SEARCH · Search PubMed

Results for “comparative genomics”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,653 records · Page 2Linked to original sources

PubMed2026

Genome-wide identification and expression analysis of the CREB/ATF family and its potential role in melanogenesis in the Manila clam (Ruditapes philippinarum).

Ruditapes philippinarum is an economically important bivalve species in China, and shell color is a trait of ecological and commercial significance. Melanin is a key determinant of shell color, and members of the CREB/ATF family have been reported to participate in melanogenesis in other organisms. In this study, members of the CREB/ATF family were systematically identified at the whole-genome level based on genomic and transcriptomic datasets, followed by analyses of their phylogenetic relationships, gene structures, and expression patterns. A total of six CREB/ATF family members were identified and classified into five subfamilies. Expression profiling and RT-qPCR validation revealed that most CREB/ATF genes were highly expressed in the mantle and displayed clear differences among shell-color phenotypes. Except for RpATF4 and RpCREBZF, most members exhibited relatively high expression levels in dark-colored shell strains, particularly in black and zebra-striped clams. Moreover, most genes showed low expression during early embryonic and larval stages but increased expression at the single-siphon spat and juvenile stages. These results suggest that the CREB/ATF family may be involved in melanin-associated shell-color regulation in R. philippinarum, providing important candidate genes and a theoretical basis for further elucidating the molecular mechanisms of shell-color formation in mollusks.

Animals↗

PubMed2026

Longitudinal whole-genome analysis of bluetongue virus identifies conserved serotype-specific genomes and distinct genomic constellations within a Colorado sheep flock (2021-2023).

Bluetongue virus (BTV) is a segmented double-stranded RNA virus of ruminants transmitted by Culicoides spp. biting midges. Although the genome consists of ten segments, classification into serotypes is primarily based on genome segment 2. However, reassortment among genomic segments is a major driver of BTV evolution and diversity. This study used longitudinal whole-genome sequencing to characterize BTV genomes collected from 2021 to 2023 within a single sheep flock in Colorado, where multiple serotypes co-circulate. Whole-genome sequences were generated from fourteen blood samples representing four serotypes: BTV-6, -11, -13, and -17. Longitudinal sampling identified multiple BTV serotypes within individual sheep across consecutive years. Tanglegram analysis comparing segment phylogenies to the segment 2 tree demonstrated incongruent topologies across all genomic segments, suggestive of reassortment or the circulation of distinct genomic constellations. Nucleotide-level comparisons revealed high sequence homology among same-serotype samples from the same year, while the greatest genetic divergence was observed among BTV-17 genomes collected in different years. Additionally, all BTV-13 genomes contained a previously undescribed nonsynonymous substitution in segment 10 predicted to extend the encoded protein by three amino acids. Together, these findings demonstrate that highly conserved BTV genomes and distinct genomic constellations can be detected at the flock level across multiple years. This longitudinal whole-genome approach reveals the genetic complexity of endemic BTV populations, including novel variants and genomic patterns consistent with reassortment that are lost with conventional serotyped-based approaches, highlighting the need to integrate whole-genome characterization into endemic BTV monitoring programs.

Animals↗

PubMed2026

Conserved host-exclusive oligonucleotide motifs enriched in pathogenic genes of human oncogenic viruses.

Comparative viral genomics can reveal sequence-level constraints influencing virus-host interactions. Relative minimal absent words (rMAWs) are short oligonucleotide motifs present in viral genomes but completely absent from the host, potentially reflecting selective pressures related to host adaptation and immune evasion. Using the EAGLE algorithm and the GRCh38 human reference genome, we systematically screened for prevalent rMAWs (prMAWs) across six major human oncogenic viruses: Epstein-Barr virus (EBV), hepatitis B virus (HBV), hepatitis C virus (HCV), human papillomavirus (HPV), human T-cell leukemia virus type 1 (HTLV-1), and human herpesvirus 8/Kaposi's sarcoma-associated herpesvirus (HHV-8/KSHV). highly conserved 11- and 12-bp prMAWs were identified in EBV, HBV, HTLV-1, and HHV-8/KSHV, with sequence prevalences ranging from 91.5% to 97.9%. Conversely, no short prMAWs were detected in HCV or HPV, likely reflecting differences in genome architecture, mutation rates, and long-term host adaptation to the human host. Importantly, the identified host-exclusive motifs exhibited non-random genomic distribution and were preferentially embedded within viral genes central to replication, persistence, immune modulation, and oncogenesis, including EBNA-1 (EBV), HBx (HBV), Tax-associated regions (HTLV-1), and lytic replication genes of HHV-8/KSHV. Notably, all detected prMAWs were enriched in GC nucleotides and exhibited marked CpG over-representation, suggesting sequence constraints associated with epigenetic regulation and viral persistence. Collectively, these highly conserved, host-exclusive signatures offer promising, candidates for sequence-directed approaches in the diagnosis, monitoring, and investigation of virus-associated cancers.

Humans↗

PubMed2026

Genome-wide identification and expression profiling of CSP and OBP genes in Stictocephala bisonia reveals candidate genes potentially associated with insecticide response.

Stictocephala bisonia is an important invasive agricultural pest. Due to the frequent application of insecticides in its habitat, this species is under intense selection pressure. Chemosensory proteins (CSPs) and odorant-binding proteins (OBPs) are known to play key roles in insecticide resistance, but their specific functions in S. bisonia remain unclear. In this study, we identified a total of 22 SbisCSPs and 16 SbisOBPs based on the S. bisonia genome. To screen for candidate genes potentially linked to insecticide resistance, we adopted a multi-criteria screening strategy that integrated phylogenetic analysis, molecular docking with three insecticides, and tissue-specific expression profiling. Phylogenetic analysis identified several SbisCSPs and SbisOBPs clustering with genes known to be involved in insecticide resistance, serving as an initial evolutionary filter. Molecular docking results indicated that λ-Cyhalothrin exhibited the strong predicted binding affinity with most of SbisCSPs and SbisOBPs. Subsequent qPCR validation of seven prioritized candidates revealed distinct expression patterns: SbisCSP22 was highly expressed in adults and demonstrated strong binding affinity to all three insecticides tested, suggesting a potential role in mediating multi-insecticide response. Conversely, SbisCSP17 was significantly upregulated in larvae, clustered with genes known to mediate imidacloprid resistance, and exhibited strong binding affinity to imidacloprid. Given its larval-specific expression and the soil-dwelling behavior of larvae, we hypothesize that SbisCSP17 is a key candidate gene for larvae coping with soil-treated insecticides.

Animals↗

PubMed2026

Genomic characterization and pathogenicity of ruminant Listeria monocytogenes isolates in a murine oral infection model.

Listeria monocytogenes is a major foodborne pathogen; its ruminant isolates display zoonotic characteristics, causing similar clinical signs in humans, including abortion and encephalitis. However, data on whole genome sequencing and pathogenicity of ruminant L. monocytogenes isolates remain sparse. This study aimed to analyze the genotypic characteristics of L. monocytogenes isolates from ruminants with listeriosis. Furthermore, we assessed the in vivo pathogenicity of four ruminant L. monocytogenes isolates, characterized via whole-genome sequencing-based genetic clustering, in orogastrically inoculated mice. The isolate LM18 (serotype 1/2b, ST224, SL6178) had the lowest lethal dose compared to the other three isolates including previous hypervirulence type (serotype 4b, ST1, SL1) and caused secondary bacteremia in lungs, with sustained bacterial loads in the spleen and liver. Genomic (listeria pathogenicity island -1 and -3) and virulence gene (actA and llsX) mutation analyses associated with virulence suggested from well-recognized studies could not elucidate the virulence of the isolates. SSI-1, which only exists in the isolate LM18 (serotype 1/2b, ST224, SL6178), may help L. monocytogenes survive in the gastrointestinal environment, thereby affecting its virulence. Further research should investigate the role of SSI-1 in the pathogenicity of L. monocytogenes. Moreover, additional studies utilizing larger datasets of ruminant isolates are required to validate our genotypic characterization and to obtain a comprehensive picture of further genotypic differences crucial for L. monocytogenes pathogenicity.

Animals↗

PubMed2026

Integrative genomic and transcriptomic analyses identify key regulators of skin pigmentation in Larimichthys crocea.

The yellow body coloration of large yellow croaker (Larimichthys crocea) constitutes a crucial economic trait, yet its underlying genetic regulatory mechanisms remain poorly understood. This study systematically elucidated the molecular basis of body color variation by integrating genome resequencing and skin transcriptome analyses, combined with the contextual analysis of key pigmentation-related genes and phenotypic histological validation. 200 phenotyped individuals (including yellow-selected lines, F1 progeny, and normal control groups, all derived from a well-characterized aquaculture stock) identified 39 significantly associated SNPs (-log₁₀(P) ≥ 6), mapping to multiple candidate genes. These genes were significantly enriched in pathways related to pigment deposition (GO:0033059), melanosome organization (GO:0032438), melanogenesis, and tyrosine metabolism. Cross-developmental stage transcriptome analysis revealed 2395 differentially expressed genes (DEGs). Multi-omics integration identified eight overlapping candidate genes, including tyrp1, slc45a2, oca2, and dgat2, among which tyrp1 was prioritized for in-depth validation based on its core regulatory role in eumelanin synthesis, significant SNP association signal, and consistent downregulation in transcriptomic data. Experimental validation demonstrated that the g.895C > T mutation in exon 2 of tyrp1b was strongly significantly associated with the yellow phenotype: the frequency of mutant genotypes (TT/CT) reached 92.86%in the yellow-selected group, whereas the control group exclusively exhibited the wild-type genotype (CC). qPCR confirmed significantly downregulated tyrp1b expression in the skin of yellow individuals, consistent with the transcriptome trend. Histological and stereomicroscopic observations of skin tissues further validated the physiological basis of the yellow phenotype, revealing a significant reduction in melanophore number and abnormal melanosome morphology in yellow-phenotype individuals, accompanied by increased xanthophore density. These results suggest that tyrp1b mutation is strongly associated with the yellow phenotype. However, the presence of a wild-type CC individual in the yellow group indicates that this mutation is not strictly required for yellow coloration, suggesting that other genetic or environmental factors may also contribute to the phenotype, Additionally, downregulation of the carotenoid metabolism gene bco2 coupled with upregulation of xdh, together with the functional changes of slc45a2 and oca2, may synergistically promote xanthophore pigment deposition, contributing to the yellow phenotype. As melanin synthesis in large yellow croaker relies on the conserved tyrosinase pathway and transporter proteins, mutations in associated genes (tyrp1b, slc45a2, oca2) represent a primary underlying cause for the loss of melanin-based coloration and transition to a yellow phenotype in L. crocea. These findings provide key molecular targets and a theoretical foundation for molecular breeding of body color in this species, and also enrich the understanding of xanthism regulatory mechanisms in teleosts.

Animals↗

PubMed2026

Genomic identification and functional characterization of the nuclear receptor gene family in relation to sex determination and gonad development in the Pacific oyster (Crassostrea gigas).

Nuclear receptors (NRs) are a large superfamily of transcription factors that control a wide range of physiological processes by modulating the expression of downstream target genes. Numerous studies have confirmed that NR family members play critical and conserved roles in sex determination and gonadal development across metazoans. However, in mollusks, systematic characterization of NRs and their potential functions in gonadal regulation remain largely unexplored. In this study, 46 NR gene family members in the Pacific oyster (Crassostrea gigas) were identified and assigned to eight subfamilies. All NR family members contain at least one of the two core domains (DNA-binding domain, DBD; ligand-binding domain, LBD), and conserved exon-intron structures were observed within the same subgroup, indicating their evolutionary conservation. Furthermore, expression profiling revealed high expression of CgNR2F, CgNR5A1-1, and CgNR0B1 in undifferentiated gonads, suggesting their potential involvement in sex determination. CgNR1A and CgNR2E5 were specifically expressed in female gonads and exhibited female-biased expression patterns, indicating a putative role in ovarian development. Moreover, CgNR3A and CgNR3B showed high expression levels during the undifferentiated stage and early male development stage, implying their possible participation in male gonadal development and gametogenesis. These results expand the understanding of the NR gene family in C. gigas and help elucidate the potential functions of NR genes in sex determination and gonadal development.

Animals↗

PubMed2026

Cost-Effectiveness and the Economics of Genomic Testing and Molecularly Matched Therapies.

Cost-effectiveness analysis of precision oncology can help guide value-driven care. Next-generation sequencing is increasingly cost-efficient over single gene testing because diagnostic algorithms require multiple individual gene tests to determine biomarker status. Matched targeted therapy is often not cost-effective due to the high cost associated with drug treatment. However, genomic profiling can promote cost-effective care by identifying patients who are unlikely to benefit from therapy. Additional applications of genomic profiling such as universal testing for hereditary cancer syndromes and germline testing in patients with cancer may represent cost-effective approaches compared with traditional history-based diagnostic methods.

Humans↗

PubMed2026

Divergent trajectories of genome architecture and chromosome evolution in ferns and angiosperms.

Ferns and angiosperms represent the two largest vascular plant lineages but exhibit striking genomic and ecological contrasts. We investigated whether differences in genome size, chromosome architecture, GC content, and stomatal traits reveal divergent evolutionary trajectories between these lineages. We assembled the most comprehensive dataset to date, integrating genome size, chromosome number and size, GC content, and stomatal traits for over 1100 fern species and compared it with an extensive angiosperm dataset. Ferns exhibited markedly lower variability and c. 16-fold slower rates of chromosome size evolution than angiosperms. A persistent positive relationship between genome size and chromosome number in ferns suggests limited cytological post-polyploid diploidization. While ferns generally possess larger stomata, this difference disappears after accounting for genome size, indicating that nucleotypic constraints, rather than lineage-specific physiology, dictate stomatal dimensions. Both groups share a unimodal GC-genome size relationship peaking at c. 14 Gbp. Larger fern chromosomes imply lower genome-wide recombination rates, potentially limiting genetic reshuffling and adaptive potential. Our results highlight fundamentally divergent evolutionary trajectories, likely shaped by meiotic symmetry in ferns and meiotic asymmetry, possibly centromere drive, and post-polyploid diploidization in angiosperms, defining the functional and genomic landscapes of these lineages across deep evolutionary timescales.

Genome, Plant↗

PubMed2026

Plant species identification by genome skimming across the vascular plant tree of life.

Accurate species identification is essential for biodiversity conservation and sustainable use, yet standard plant DNA barcoding often fails to achieve species-level resolution. We present a large-scale empirical evaluation of genome skimming as a tool to improve plant species discrimination. Using standardised data from 1969 individuals representing 475 species from 32 genera across major lineages of the vascular plant tree of life, we compare conventional plastid + internal transcribed spacer (ITS) barcodes with genome skimming approaches. Standard barcoding using rbcL, matK, trnH-psbA and ITS resolved about half of species (49.3%), with six genera showing < 25% species discrimination. By contrast, genome skimming enabled the recovery of complete plastid genomes, yielding 57.6% species discrimination. It also generated sufficient nuclear genomic data for additional resolution from k-mer analysis, achieving 66.8% species discrimination - an average gain of 17.5% over standard barcodes - while eliminating cases of extreme failure (< 25% resolution). The recovery of complete plastomes and ribosomal DNAs from genome skims also ensures backward compatibility with existing barcode datasets. Our results demonstrate that genome skimming provides data that substantially improves species-level resolution across diverse plant lineages and offers a scalable, high-throughput approach for building comprehensive reference resources to support global biodiversity initiatives.

DNA Barcoding, Taxonomic↗

PubMed2026

Genome-wide characterisation of the myosin light chain gene family in Chinese perch (Siniperca chuatsi) and its expression patterns in muscle fibre types and injury response.

The Class II myosin light chain (myl) genes in Chinese perch (Siniperca chuatsi) have not yet been systematically characterised, and relationships with muscle fibre specification, development, and injury-associated remodelling remain unclear. In this study, fast and slow muscle fibres were initially distinguished using myofibrillar ATPase histochemistry. Subsequently, genome-wide mining identified 16 Class II myl genes, comprising eight essential and eight regulatory light-chain subunits. Their conserved-domain features, chromosomal distribution, phylogenetic relationships and expression profiles were analysed. Transcriptomic profiling showed that summed myl transcript abundance was higher in fast muscle than in slow muscle, accounting for 67.2% of the pooled myl transcript pool across the two muscle types (paired t-test, raw P = 0.036). mylpfa, myl1 and mylz3 were the major fast-muscle-associated genes, whereas myl10, myl2b and myl13 were preferentially expressed in slow muscle at the transcript level. These patterns support these genes as candidate fibre-type-associated expression markers. Developmental profiling identified stage-associated myl expression patterns, including a possible expression shift between mylpfb and mylpfa. In the descriptive injury-repair time course (d0-d7), FPKM profiles indicated that fast-muscle-associated genes (mylpfa, mylz3 and myl1) were lower at d1 and recovered by d3, whereas several slow-muscle-associated genes showed biphasic transcript-level increases. The slow-muscle-associated RLC gene mylpfb showed a delayed expression peak at d7. Notably, the embryonic isoform myl6l showed a modest increase from approximately 2 FPKM at d0 to 4-5 FPKM after injury, suggesting a possible injury-associated expression pattern that requires further validation. Together, these findings provide a genome-wide description of the Chinese perch myl gene family and identify candidate fibre-type-associated genes and descriptive injury-associated isoform expression patterns.

Animals↗

PubMed2026

A conserved distal-tail helical extension defines a tailspike attachment architecture in Gram-negative siphophages.

Rapid growth of bacteriophage genome collections has outpaced functional annotation of tail-tip proteins, limiting comparative analysis of host-recognition structures. Starting from a shared distal-tail gene organization in the Salmonella phages 9NA and Jersey, I developed a morphogenetic bioinformatic framework integrating gene synteny, sequence comparison, profile hidden Markov model (HMM) screening, structural evidence, structure-aware searching, and AlphaFold modeling. Comparison with the experimentally characterized lambda and Sf11 tail assemblies identified a predominantly alpha-helical C-terminal extension of the distal-tail (DT) protein associated with tailspike attachment, termed the distal-tail helical extension (DT-helix). Screening 541,986 proteins from 5167 complete NCBI RefSeq tailed-phage genomes, followed by evidence-based evaluation of sequence, genomic context, and structural architecture, identified 165 curated DT-helical-extension-associated phages. Their DT proteins segregated into six sequence groups. In the four principal multi-member groups, cognate tailspikes showed group-specific conservation in proximal N-terminal regions but substantially greater downstream diversity, consistent with sequence constraint at the DT-tailspike attachment boundary. A complementary ProstT5/Foldseek search supported the established groups but revealed no convincing additional highly divergent family. Together with the experimentally characterized Sf11 attachment interface, these findings define a recurrent morphogenetic architecture linking conserved distal-tail scaffolds to more variable receptor-binding proteins across siphophages infecting Gram-negative bacteria. Although universal exchangeability is not established, the identified scaffold-receptor-binding boundaries provide a framework for molecular characterization and rational phage engineering. Accession-level information for the 165 curated phages is available through PhageTailDB.

Viral Tail Proteins↗

PubMed2026

Integrative modeling of the genome structure and dynamics in fission yeast.

Genome organization in the nucleus is highly structured and dynamic. Recent advances in genomic technology have enabled the measurement of genome-wide architecture and locus-specific motion, yielding contact maps and live-cell trajectories. However, these outcomes are derived from different modalities and are not directly comparable, with their quantitative integration being a key challenge. Here we establish a genome-wide live-cell imaging platform in fission yeast Schizosaccharomyces pombe, tracking 131 chromosomal loci, along with the spindle pole body (SPB) and nucleolus, to construct a quantitative map of locus dynamics. By integrating these dynamics with contact data through polymer modeling of Hi-C data, we build a physics-based "digital twin" of the S. pombe genome consistent with the spatiotemporal dynamics of interphase chromatin. We validate it against genome-wide mobility patterns and known architectural features, including centromere and telomere clustering. The model also identifies distinct dynamical regimes: centromere- and telomere-proximal loci relax within [Formula: see text]150 s, whereas the remaining loci relax within [Formula: see text]70 s. We measure semiperiodic dynamics of SPB motion, including a characteristic peak near 225 s and [Formula: see text] fluctuations. We use the model with SPB-directed forcing to show how these low-frequency components propagate through the genome to drive genome-wide chromatin displacements. Together, this predictive physics-based modeling framework integrates genome structure and dynamics to reveal how nuclear mechanical driving forces shape chromosome motion, linking mechanically driven chromatin responses to genome maintenance and regulation.

Schizosaccharomyces↗

PubMed2026

Ramu stunt virus genome reveals previously unreported segments and nucleocapsid domain duplication in Mechlorovirus.

Ramu stunt virus (RmSV), a member of the genus Mechlorovirus within the family Phenuiviridae, was previously described as a six-segmented RNA virus infecting sugarcane. In this study, we re-examined type material and additional isolates using high-throughput sequencing and RT-PCR validation, revealing that RmSV possesses a nine-segmented genome, making it the largest reported in the Phenuiviridae. This expanded architecture includes duplicated RNA segments (RNA 2a and RNA 2b) encoding nucleocapsid-like proteins and two novel segments (RNA 7 and RNA 8). Comparative analysis showed that RNA 2a and 2b share about 84% amino acid identity, while RNA 5 encodes a third nucleocapsid homolog, indicating unprecedented domain redundancy. Structural modeling confirmed that all three nucleocapsid proteins maintain a conserved fold despite low sequence identity, with electrostatic mapping suggesting differential RNA-binding potential. Additionally, RNA 6 encodes a hypothetical protein structurally similar to the rice stripe virus disease-specific S-protein, implicating a role in symptom development. Transcript abundance analysis revealed RNA 6 as the most highly expressed segment across isolates. These findings revise the genomic composition of RmSV, highlight mechanisms of genome plasticity and adaptive evolution in plant-infecting bunyaviruses, and underscore practical implications for diagnostic assay design, resistance breeding, and biosecurity surveillance.

Genome, Viral↗

PubMed2026

A chromosome-level, haplotype-resolved genome assembly for the barn owl, Tyto alba.

Recent advances in long-read sequencing have enabled near telomere-to-telomere (T2T) assemblies across diverse taxa. However, avian genomes remain challenging due to numerous microchromosomes, small, typically < 20Mb, DNA molecules that are gene-, GC-, and repeat-rich. As a consequence, microchromosomes are often missing from genome assemblies. Here, we present a chromosome-level, haplotype-resolved genome assembly for the Western barn owl (Tyto alba). Using a trio-binning strategy with Illumina parental reads combined with PacBio HiFi and Oxford Nanopore Technologies data, we generated two phased contig sets. These were scaffolded into 40 linkage groups using a linkage map. Comparative analyses identified unplaced HiFi scaffolds corresponding to microchromosomes, which we integrated into six additional microchromosomes using long reads information. The two assemblies present 46 chromosomes, matching the karyotype of the species. They exhibit strong synteny between parental haplotypes, except for a ∼38 Mb complex region on chromosome 7 containing nested inversions. This high-quality reference provides a haplotype-resolved and chromosome-level genome for Strigiformes, enabling fine-scale studies of structural variation and avian genome evolution.

Tyto alba↗

PubMed2026

Exploratory identification and cellular functional characterization of ppiabl as a candidate gene associated with growth traits in Paralichthys olivaceus.

The Japanese flounder (Paralichthys olivaceus) is an important mariculture species. However, the genetic mechanisms underlying its growth traits remain poorly understood. To explore the genetic basis of growth variation, whole-genome resequencing was performed in a cultured cohort of 60 individuals, followed by exploratory genome-wide association analysis and candidate-gene prioritization. The results revealed heritability estimates of 0.40 for body weight and 0.24 for body length, with substantial overlap in associated loci between the two traits. Exploratory association and variant-annotation analyses prioritized ppiabl, which carries a nonconservative missense variant, as a candidate gene for further investigation. Tissue expression analysis showed that ppiabl was highly expressed in muscle tissue. This gene encodes a protein belonging to the conserved peptidyl-prolyl cis-trans isomerase family. In Japanese flounder primary muscle cells, ppiabl knockdown was associated with altered expression of growth-related genes and an increased G1-phase fraction, whereas overexpression produced changes in the opposite direction. In line with this, fast-growing individuals were found to have significantly larger muscle fiber areas than slow-growing ones. These findings suggest that ppiabl may be involved in muscle-related cellular processes associated with growth variation in Japanese flounder, although its contribution to whole-animal growth requires further validation. Overall, this exploratory study prioritizes ppiabl as a candidate gene potentially associated with growth-related cellular processes in Japanese flounder, although validation in larger independent populations and in vivo models is required.

Animals↗

PubMed2026

Genome-wide identification and expression profiling of HSD3B and SDR42E1 genes in the Pacific oyster (Crassostrea gigas): potential associations with gonadal development.

Sex steroids are lipid-soluble signaling molecules that regulate sex differentiation, reproductive development and physiological homeostasis in animals. 3β-Hydroxysteroid dehydrogenase/Δ5-Δ4 isomerase (3β-HSD) is a key steroidogenic enzyme, whereas SDR42E1, an extended short-chain dehydrogenase/reductase, has been implicated in sterol- and steroid-related metabolism. However, the composition, evolutionary relationships and expression patterns of the HSD3B- and SDR42E1-related genes in bivalve gonadal development remain poorly characterized. In this study, five PF01073-containing genes, comprising three CgHsd3b and two CgSdr42e1 genes, were identified in the Pacific oyster Crassostrea gigas. Phylogenetic analysis separated the proteins into HSD3B-related and SDR42E1-related groups, and gene-structure and motif analyses indicated subfamily-level divergence. All five proteins retained the SDR domain but differed in exon-intron structure and motif composition. Each contained the extended-SDR TGxxGxxG motif, whereas exact classical [ST]GxxxGxG and NNAG motifs were absent. Tyr- and Lys-equivalent residues were conserved, while the HSD3B1 Ser-equivalent position contained Thr in two C. gigas proteins and Ser in one. These features support their classification as extended-SDR proteins but do not establish enzymatic activity or substrate specificity. The three CgHsd3b genes were dispersed on one chromosome, whereas CgSdr42e1-1 and CgSdr42e1-2 were adjacent on another chromosome, suggesting a possible local duplication event for the CgSdr42e1 pair. Public RNA-seq data showed distinct tissue- and gonadal-stage expression patterns, with several genes displaying gonad-biased or female-stage-associated expression. Independent RT-qPCR profiling of the representative genes CgHsd3b-3 and CgSdr42e1-1 detected stage-dependent expression, although tissue rankings differed from those in the public RNA-seq datasets. These differences may reflect the use of independent biological samples, tissue composition, normalization procedures, and platform-specific measurements. Because enzymatic assays, metabolite measurements, cellular localization, and functional perturbation were not performed, the results identify candidate genes whose expression is associated with gonadal development rather than demonstrating regulatory roles. This study provides a comparative framework for future functional investigation of sterol- and steroid-related metabolism in bivalves.

Animals↗

PubMed2026

Genomic and food-safety evaluation of Staphylococcus chromogenes in Chinese dairy milk.

Non-aureus staphylococci and mammaliicocci (NASM) cause mastitis and may contaminate milk and dairy products. Milk samples (n = 1916) from cows with subclinical or clinical mastitis (SCM and CM, respectively) were collected from 28 large-scale (> 500 lactating cows) Chinese dairy farms. Overall, 999 NASM isolates representing 19 species were identified by MALDI-TOF MS and cpn60 sequencing, with Staphylococcuschromogenes, Mammaliicoccus sciuri and Staphylococcus haemolyticus being most prevalent. Antimicrobial resistance (AMR) was determined with disc diffusion; non-susceptible to penicillin was most common (SCM, 30% and CM, 29%) whereas cefoxitin non-susceptible NASM accounted for 8-10% of isolates; among these, 12.5% carried mecA but none carried mecC. Galleria mellonella was used to assess virulence of 78 strains of S. chromogenes, a dominant species; subsequently, 32 strains, representing higher- and lower-virulence in the Galleria model, were selected for whole-genome sequencing and comparative genomics. S. chromogenes isolates from CM had higher virulence (p < 0.05) than those from SCM. The 32 genomes comprised 20 sequence types, indicating high genetic diversity. No robust genomic marker of Galleria virulence phenotype was identified in this selected WGS subset. Acquired resistance genes (n = 5) were detected, including a first report of fusC in S. chromogenes; the fusC-positive isolate had an elevated fusidic acid MIC (8 mg/L). Although S. chromogenes persisted in milk at 4 °C, pasteurization (64 °C for 30 min) reduced viable counts to below detection. This study provided new insights into the prevalence, AMR, genomic diversity, and dairy-chain relevance of milk-derived NASM, particularly S. chromogenes. However, the genomic findings were based on an intentionally selected WGS subset and should be interpreted as hypothesis-generating rather than population-representative.

Animals↗