Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Viral 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.

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

217 records · Page 3Linked to original sources
PubMed2026

A multi-model genome-wide association study identifies genetic variants underlying resistance to Largemouth Bass Ranavirus (LMBV) in Micropterus salmoides.

Largemouth bass (Micropterus salmoides) is an economically important freshwater aquaculture species, yet recurrent outbreaks of Largemouth Bass Ranavirus (LMBV) continue to impair production and cause substantial losses. The genetic basis of host variation in LMBV resistance remains insufficiently characterized. Here, we applied a multi-model genome-wide association study (GWAS) to identify loci associated with resistance following a controlled challenge with the LMBV-23PY strain. Whole-genome resequencing was performed for 146 phenotyped fish, including 72 susceptible and 74 resistant individuals. After stringent quality control, 877,262 high-quality variants were retained and tested using six GWAS models. Across binary survival status and survival time phenotypes, 32 shared suggestive variants were consistently detected across models, representing suggestive loci for LMBV-23PY resistance. Genes within ±50 kb of these loci were annotated, and functional enrichment highlighted immune- and redox-related biological processes. Three prioritized candidates-GSTT3L (glutathione S-transferase theta-3-like), CGRP2 (calcitonin gene-related peptide 2), and NPPC (natriuretic peptide C)-were associated with pathways involved in oxidative stress responses and immune regulation. Collectively, these results provide insight into the genetic architecture of LMBV-23PY resistance in largemouth bass and identify suggestive variants and associated candidate genes for downstream validation, functional interrogation, and the development of marker-assisted and genome-enabled breeding strategies.

Animals↗
PubMed2026

Comparative genomic and proteomic analysis reveals orthogroup structured evolution of tick protease inhibitors.

Protease inhibitors (PIs) play central roles in regulating endogenous proteolysis and host-parasite interactions in ticks. However, the evolutionary architecture underlying their diversification across tick lineages remains insufficiently resolved. Here, we performed a genome-wide comparative analysis of predicted proteomes from 14 tick species to systematically characterize PI repertoires. In total, 4931 putative PIs were identified and grouped into 20 families using the MEROPS classification system. Further, PI families such as Antistasin, WAP-type, and Pacifastin, which have not previously been systematically reported in tick genomes, were classified. Orthogroup inference demonstrated that PI expansion is structured at the level of evolutionary lineages rather than uniformly across families. By stratifying orthogroups according to duplication burden and taxonomic conservation, we identified a broadly conserved single-copy core under strong purifying selection. Motif level analysis of serpin reactive center loops further revealed conservation of inhibitory specificity within single copy orthogroups and diversification of key functional residues in duplication-associated lineages. Integration of secretion prediction and tissue-resolved proteomics from Hyalomma anatolicum and Rhipicephalus microplus demonstrated that evolutionary stratification is reflected at the protein level. Together, these findings provide an orthogroup-resolved evolutionary framework linking duplication dynamics, molecular evolution, and tissue-level protein deployment. This integrative approach offers a systematic basis for prioritizing conserved and diversified PI lineages for future functional and anti-tick intervention studies.

Animals↗
PubMed2026

Routine methods misidentify Serratia spp.: Limitations of MALDI-TOF MS revealed by whole-genome sequencing.

Accurate species-level identification within the genus Serratia remains challenging due to extensive phenotypic overlap and high genomic relatedness among closely related and recently described taxa. This study presents an evaluation of routine and genome-based identification approaches applied to clinical Serratia isolates, integrating phenotypic assays, MALDI-TOF MS (Bruker Daltonics), 16S rRNA gene sequencing, and Whole-Genome Sequencing (WGS). A total of 103 isolates collected from a teaching hospital were analyzed. WGS was performed on a subset of isolates. Conventional biochemical methods classified all isolates as Serratia marcescens, whereas MALDI-TOF MS identified 60.1% as S. marcescens, 11.6% as S. ureilytica, and 28.1% just at the genus level. Peak analysis from MALDI-TOF MS revealed specific peaks associated with S. marcescens and S. ureilytica, but limited discriminatory power. WGS of six isolates initially identified as S. ureilytica by MALDI-TOF MS revealed reclassification as Serratia sarumanii (n = 5) and Serratia montpellierensis (n = 1), supported by Average Nucleotide Identity (ANI), Average Amino Acid Identity (AAI), and Digital DNA-DNA Hybridization (dDDH) thresholds. In contrast, 16S rRNA analysis showed limited species-level resolution. Phylogenomic and SNP-based analyses confirmed these classifications with strong support. Overall, this study underscores the critical role of high-resolution genomic approaches for precise species identification and highlights the need for continuous expansion and curation of MALDI-TOF MS reference databases to support reliable clinical diagnostics and epidemiological surveillance of emerging Serratia species.

Spectrometry, Mass, Matrix-Assisted Laser Desorpti↗
PubMed2026

Comparative analysis of genomic variations among different Cdo1 paralogs for salinity-adaptation in oysters.

Under rapid climate change and anthropogenic activities, oysters, a global aquaculture species, are subjected to exacerbated culturing environments, especially for those living in in-shore estuarine species, such as Suminoe oysters Crassostrea ariakensis. This study aims to investigate the molecular mechanisms of salinity adaptation of C. ariakensis. We performed an expression genome-wide association study (eGWAS) to compare genetic regulation among 5 paralogous copies of a key salinity-related gene, cysteine dioxygenase 1 (Cdo1). A total of 40 significant eSNPs with 82 adjacent eGenes were identified in 2 copies (Cdo1_26639 and Cdo1_1666). We identified only trans-eSNPs for Cdo1_26639 and more cis-eSNPs for Cdo1_1666, and different eGenes for these 2 Cdo1 copies, which indicated that the expressional regulation of these paralogs may undergo distinct pathways. We identified 3 eGenes that exhibited identical expression patterns with Cdo1_26639 and Cdo1_1666, including 6-Pgdh, Trapp and tandem copy of Cdo1_27337. The expression correlation between Cdo1 copies and eGenes was enhanced under salinity stresses, suggesting the crucial role of eGenes in regulating Cdo1's expression in response to salinity changes. Our results provide comprehensive identification and comparison of eSNPs across different paralogous copies of one gene, along with insights into the molecular mechanisms underlying salinity tolerance, and genetic markers for breeding salinity-resistant oysters.

Animals↗
PubMed2026

Systematic modular engineering of genome-integrated Escherichia coli MG1655 for high-level 2'-fucosyllactose production.

2'-Fucosyllactose (2'-FL), the most abundant human milk oligosaccharide (HMO), has attracted considerable interest for its prebiotic and immunomodulatory functions, with broad applications in infant nutrition. In this study, we report the development of a high-yield, genome-integrated 2'-FL-producing strain based on Escherichia coli MG1655 through systematic modular optimization. Starting from a single-copy BKHT strain (MGC06), we first optimized the copy number of the α-1,2-fucosyltransferase (α-1,2-FT) gene BKHT. Subsequently, the GDP-L-fucose supply was enhanced through coordinated genomic integration of the gene clusters cpsG-cpsB and gmd-fcl, while the multidrug efflux transporter gene mdfA was integrated to improve product export and strain robustness. BKHT copy number was then re-evaluated in the optimized background, with four copies yielding the highest production. The final engineered strain, harboring all genetic modifications stably integrated into the chromosome, produced 17.18 g/L 2'-FL in shake-flask culture. In fed-batch fermentation using a 5-L bioreactor, this strain achieved a titer of 154.12 g/L after 60 h, with a productivity of 2.57 g/L/h. Notably, throughout the entire fermentation process, no antibiotics or inducers were supplemented, underscoring the genetic stability and regulatory compliance of this plasmid-free system. To our knowledge, this represents the highest 2'-FL titer reported to date, positioning our engineered strain as a promising candidate for commercial 2'-FL production.

Escherichia coli↗
PubMed2026

Genome-wide identification of the HSP70 superfamily in tropical sea cucumber Stichopus monotuberculatus and their expression analysis under low-salinity stress.

Heat shock proteins (HSPs) are a group of evolutionarily conserved molecular chaperones that serve as indispensable core regulators in preserving cellular homeostasis and orchestrating organismal stress responses. The tropical sea cucumber Stichopus monotuberculatus, a high-value aquaculture species, is sensitive to fluctuations in environmental salinity-a challenge that has emerged as a critical bottleneck limiting its large-scale commercial cultivation. However, no systematic investigation has been conducted to characterize the HSP70 superfamily in S. monotuberculatus and elucidate its functional roles in salinity adaptation. In the present study, we performed a comprehensive genome-wide scan and identified 19 HSP70 superfamily genes in the S. monotuberculatus genome, with the HSP70IV subfamily showing remarkable gene expansion, containing 8 distinct copies. Phylogenetic analysis, conserved motif identification, and gene structure characterization demonstrated high evolutionary conservation within each HSP subfamily. These genes were unevenly distributed across the chromosomes of S. monotuberculatus, and prediction of cis-acting elements revealed that their upstream regulatory regions were enriched with numerous functional elements associated with stress response and immune regulation. Salinity stress experiments revealed that under severe low-salinity conditions (18‰), the expression levels of SmHSPA14L and multiple HSP70IV subfamily members were significantly elevated, while SmHYOU1D was significantly downregulated; in contrast, only subtle changes were detected in the expression of most HSP70 genes under moderate low-salinity stress (24‰). These findings strongly suggest that HSP70 genes, particularly the expanded HSP70IV subfamily, may act as key modulators in the low-salinity stress response. This work provides valuable insight into the molecular mechanisms underlying salinity adaptation in tropical sea cucumbers.

Animals↗
PubMed2026

Integrating genomic distance analyses in the description of a new family, genus, and species of sponge-associated antipatharians (black corals).

Antipatharians (black corals) are among the least studied coral groups, with much of their diversity still undescribed. Here, we present an integrative morphological, phylogenomic and genomic distance study of deep-sea antipatharians sampled in high seas areas of the North Pacific Ocean and from New Zealand's Exclusive Economic Zone. These corals grow on hexactinellid sponges - a unique characteristic in the order Antipatharia. Using a dataset of ultra-conserved elements and exons, combined with morphological analyses, we reconstruct phylogenomic relationships and formally describe a new family (Eidikopathidae fam. nov.), a new genus (Eidikopathesgen. nov.), and two new species (E. korallispongiasp. nov., E. zealandkoralliasp. nov.). Morphologically, the new family is distinguished by a corallum consisting of a network of loose branches that fuse with the sponge skeletal framework. Phylogenomic analyses recovered consistent topologies with strong nodal support, corroborating the distinct evolutionary placement of this sponge-associated lineage. Pairwise genomic distances estimated using the Tamura-Nei model were concordant with patristic genomic distances, identifying Pteridopathidae as the genetically closest family to Eidikopathidae fam. nov., followed by Myriopathidae and Stylopathidae, which were recovered as sister families in the phylogeny. This pattern shows that genomic distance complements, rather than simply mirrors, tree topology by quantifying accumulated sequence divergence among lineages. Together, these results provide the first genomic distance framework for Antipatharia, offering a baseline for future systematic, evolutionary, and biodiversity studies on this fundamental shallow, mesophotic and deep-sea coral group.

Animals↗
PubMed2026

Regional genomic analysis of lineage distribution and transferable multidrug resistance among chicken-associated Salmonella Kentucky isolates in China.

Salmonella enterica serovar Kentucky is an important multidrug-resistant foodborne pathogen in the poultry meat supply chain. Although recent broader genomic studies have elucidated the population structure and epidemiological significance of major lineages in China (e.g., ST198 and ST314), the regional dynamics within local poultry supply chains remain insufficiently characterized. In this study, 31 chicken meat-derived isolates from Shanghai and 39 publicly available genomes from China were analyzed using antimicrobial susceptibility testing, whole-genome sequencing, phylogenetic analysis, conjugation experiments, and complete sequencing of representative plasmids. This enabled a systematic characterization of the molecular epidemiological features of the population and the mechanisms underlying resistance dissemination. Population genomic analysis revealed a lineage composition markedly different from the global epidemiological pattern: ST314 was the predominant sequence type among the Shanghai chicken-derived isolates (74.2%), whereas the internationally recognized high-risk clone ST198 accounted for only 25.8% of the local isolates. However, risk stratification analysis indicated that although ST198 was detected less frequently, it carried a significantly greater burden of acquired resistance genes and therefore represented a higher-risk resistant lineage. Functional and structural validation further elucidated the molecular basis of resistance dissemination within this high-risk lineage. Conjugation experiments confirmed the co-transfer of a multidrug resistance module carrying blaTEM-1 and blaCTX-M-267 to the recipient strain Escherichia coli J53. Complete plasmid analysis revealed that these two β-lactam resistance genes were co-localized on a 242-kb transferable plasmid flanked by Tn1331, Tn3, and multiple transposase-associated elements, thereby providing a structural basis for their horizontal transfer. This study provides important molecular epidemiological evidence for lineage-specific surveillance and risk-stratified control of resistant Salmonella in the poultry meat supply chain and further underscores the need for continuous monitoring of mobile genetic elements within a One Health framework.

Animals↗
PubMed2026

Machine learning-ready genomic biomarkers: ATF3 polymorphisms predict postoperative analgesic demand through AI-compatible phenotyping.

PURPOSE: To determine whether ATF3 polymorphisms can serve as genetic biomarkers for machine learning-based precision analgesia by establishing a genotype-phenotype association suitable for predictive modeling of postoperative opioid requirements. METHODS: In a prospective cohort of 167 adults undergoing abdominal surgery, ATF3 SNPs rs3122721 and rs3125293 were genotyped. A structured dataset architecture was developed to represent genetic profiles as input features for supervised learning models, enabling translational analysis of genotype‑dependent opioid consumption over 72 h. RESULTS: Patients with homozygous genotypes of the ATF3 SNPs had significantly higher opioid requirements than non‑carriers, despite reporting similar subjective pain scores. This consistent genotype‑dependent pattern provided a clinically relevant phenotype suitable for integration into predictive algorithms. CONCLUSION: ATF3 genotyping offers a promising biomarker for computationally informed precision analgesia. By linking genomic variability to clinically meaningful outcomes within a structured clinical and genomic framework, this approach supports the future development of risk-stratified clinical decision-support systems to optimize postoperative pain management.Trial registration ChiCTR1900021991, registered 30 April 2019. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at https://doi.org/10.1007/s13755-026-00480-9.

ATF3↗
PubMed2026

Comparative analyses of olfactory receptor repertoires in Schizothorax fish based on the chromosome-level genomes: Implications for regulatory roles of dietary differentiation and ploidy variation.

The olfactory receptor (OR) genes constitute the molecular basis of fish olfaction, mediating survival behaviors and environmental adaptation while coevolving with habitat-driven evolution. Schizothorax, a cyprinid genus endemic to the Qinghai-Tibetan Plateau, exhibits remarkable dietary divergence and ploidy variation in response to plateau environmental changes, which presumably facilitates the adaptive evolution of OR genes. However, the evolutionary patterns of OR genes associated with trophic divergence and ploidy variation in this genus remain unclear. In this study, three species were selected: the herbivorous diploid S. macropogon, the carnivorous diploid S. lantsangensis, and the herbivorous tetraploid S. curvilabiatus. S. macropogon possessed 142 OR genes (92.25% functional), primarily located on chromosomes 14 and 24, with the fewest sequence clusters. Such compact gene repertoire and highly overlapping chromosomal clusters indicated specialization for a herbivorous olfactory niche. S. lantsangensis contained 127 OR genes (93.70% functional), concentrated on chromosomes 4 and 5, with fewer sequence clusters and a scattered distribution, reflecting evolution of OR genes under carnivorous feeding habits. The herbivorous tetraploid S. curvilabiatus exhibited striking features: 316 OR genes (94.30% functional), the most subfamilies, unique ε and κ OR subfamilies, and species-specific motifs. These characteristics revealed that ploidy, rather than herbivory, dominated OR gene evolution. In conclusion, dietary differentiation and ploidy variation together drove olfactory adaptive evolution in Schizothorax, providing new insights into vertebrate OR gene ecological adaptation.

Animals↗
PubMed2026

Evolutionary architecture and lineage-specific diversification of Forkhead box transcription factors in Perna viridis.

The Forkhead box (Fox) transcription factors are evolutionarily conserved regulators of development, cell cycle, and apoptosis across metazoans. This study provides the first comprehensive genome-wide analysis of the Fox gene family in the Asian green mussel (Perna viridis). We identified 28 Fox genes distributed across 10 chromosomes. Comparative analysis reveals the absence of the FoxI, FoxQ1, FoxR and FoxS subfamily, consistent with other bivalves and indicative of lineage-specific gene loss during molluscan evolution. Notably, gene duplications in the FoxAB, FoxD, FoxH, FoxN1-4, FoxQ2 and FoxQD subfamilies may reflect functional diversification associated with environmental adaptation. Exon-intron structural variability, including intron loss in several paralogues, suggests structural diversification and potential regulatory variation. Phylogenetic reconstruction confirmed the monophyly of core Fox classes while highlighting divergent expansion patterns in lophotrochozoans. Selection analyses showed strong purifying selection across duplicated Fox paralogs, supporting functional conservation after lineage-specific expansion. Gene Ontology enrichment linked Fox genes to stress response, apoptosis, and transcriptional regulation. By integrating phylogenetic, structural, and transcriptomic analyses, this study provides a genomic framework for understanding Fox gene organisation, evolution, and tissue-associated expression patterns in Perna viridis and establishes a comparative resource for future functional studies in bivalves.

Animals↗
PubMed2026

Genomic epidemiology of clinically critical antibiotic resistance in Salmonella enterica causing bloodstream infections across six Chinese provinces, 1994-2023.

Clinically critical antibiotic-resistant Salmonella enterica (S. enterica) causing bloodstream infections remains a public health challenge. Here, we aim to reveal the emergence and trends of clinically important antibiotic resistance in S. enterica causing bloodstream infections using 833 isolates from six Chinese provincial-level administrative areas during 1994-2023. We identified 48 serovars and 64 sequence types (STs). Overall, 8.52% of 833 isolates were resistant or had decreased susceptibility to ciprofloxacin, 4.32% and 6.84% reported resistance or decreased susceptibility to third- and fourth-generation cephalosporins (3GCs and 4GCs), 1.80% reported resistance to fosfomycin, and 2.16% reported resistance to azithromycin. Across these six regions, azithromycin and fosfomycin resistance is increasing, as is decreased susceptibility or resistance to ciprofloxacin, 3GCs, and 4GCs, especially among younger children and elderly people. Clinically prioritized antibiotic resistance also varies by region, serovar, and age group. S. Paratyphi A genotype 2.3.3 strains are mainly divided into 2 lineages distributed in Guangxi and Shanghai. Within the scope of this passive surveillance dataset, S. Typhi genotype 4.3.1.2.1 was identified as the earliest documented case among the collected isolates. Our retrospective and longitudinal genomic epidemiology study provides critical data for the formulation of treatment guidelines and policies for bloodstream infections and for the monitoring and control of antimicrobial resistance.

Humans↗
PubMed2026

Genome-wide characterization of the TGF-β superfamily identifies bmp15, gdf9, and gsdf as sex-biased candidate regulators of gonadal differentiation in the synchronous hermaphrodite Plectropomus leopardus.

The transforming growth factor-β (TGF-β) superfamily plays conserved roles in vertebrate reproduction and gonadal sex differentiation. However, its genomic repertoire and sex-biased expression patterns remain unclear in the leopard coral grouper (Plectropomus leopardus), a species with synchronous hermaphroditism. Here, we performed a genome-wide identification of the TGF-β superfamily, identifying 42 genes from the chromosome-level genome. Phylogenetic and synteny analyses indicated that segmental duplication under purifying selection contributed to family expansion. Expression profiling across multiple tissues and four gonadal developmental stages (undifferentiated, 120 dph; early differentiated, 15 months; mature testis, 3 years; mature ovary, 3 years) identified eight gonad-enriched genes, among which bmp15 and gdf9 exhibited pronounced female-biased expression, with transcripts localized exclusively to the oocyte cytoplasm, particularly in stage II-III oocytes. In contrast, gsdf showed male-biased expression and was localized in spermatogenic cells of the testis. These reciprocal expression patterns indicate that bmp15/gdf9 and gsdf are candidate factors associated with gonadal sex differentiation. Our study provides the first comprehensive characterization of the TGF-β superfamily in P. leopardus and highlights bmp15, gdf9, and gsdf as candidate sex-differentiation factors in this hermaphroditic species.

Animals↗
PubMed2026

PGR expression as a pharmacogenomic companion biomarker to GENE70-derived genomic risk in ER-positive/HER2-negative breast cancer.

BACKGROUND: The biology of the estrogen receptor-positive (ER+) and human epidermal growth factor receptor 2-negative (HER2-) breast cancers is heterogeneous even when they are categorized by their risk via genomics. Transcriptomic PGR expression reflects endocrine pathway activity and may provide complementary biological information within established GENE70-derived genomic-risk categories. Whether this molecular marker improves the biological interpretation of genomic-risk stratification beyond conventional clinicopathological assessment remains uncertain. OBJECTIVES: The aim of this study was to determine whether transcriptomic PGR expression provides complementary biological and prognostic information within reconstructed GENE70-derived genomic-risk categories and refines the characterization of endocrine-related tumour biology in ER-positive/HER2-negative breast cancer. METHODS: This study analysed publicly available transcriptomic and clinical data from three cohorts: METABRIC (discovery cohort), GSE96058/SCAN-B cohort (validation cohort) and TCGA-BRCA cohort (molecular validation cohort). The GENE70-derived genomic-risk score was reconstructed for each cohort using matched genes. Cox regression, Kaplan-Meier analysis and subgroup comparisons were used to assess relationships between PGR expression, clinicopathologic variables, molecular features and survival outcomes. RESULTS: Across the three independent cohorts, low transcriptomic PGR expression was consistently associated with higher GENE70-derived genomic risk, increased MKI67 expression, reduced ESR1 expression and enrichment of the Luminal B subtype. Survival findings differed between cohorts. In the discovery METABRIC cohort, transcriptomic PGR expression showed heterogeneous associations with survival, particularly within GENE70-derived high-risk subgroups, whereas the external GSE96058/SCAN-B validation cohort demonstrated consistent associations between low PGR expression and poorer overall survival in both the overall ER-positive/HER2-negative population and GENE70-derived high-risk subgroups. CONCLUSION: These findings suggest that transcriptomic PGR provides complementary biological and prognostic information within GENE70-derived genomic-risk categories. However, because treatment response was not evaluated in the present study, the findings should not be interpreted as evidence of predictive or pharmacogenomic utility and prospective studies incorporating treatment-response analyses are required before such applications can be established.

Humans↗
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

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

Genome-wide SNP data support species boundaries in sympatric Polylepis Ruiz & Pav. (Rosaceae) species from Bolivia and Ecuador.

Species delimitation in the South American genus Polylepis is notoriously challenging due to high morphological similarity and phenotypic plasticity, likely driven by hybridization and gene flow. Previous phylogenetic studies suggested that genetic structure aligns more strongly with geography than with taxonomy, questioning existing species concepts and hampering conservation efforts. We used double-digest RAD sequencing (ddRADseq) to generate genome-wide SNP data for 11 Polylepis species sampled across multiple localities in Bolivia and Ecuador. Population genetic analyses, phylogenetic inference, and network approaches were combined to assess whether genetic structure aligns more closely with taxonomy or geography. Morphologically defined species formed largely cohesive genetic lineages across regions, with species identity explaining substantially more genetic variation than locality. While localized admixture and reticulation were detected among closely related taxa, widespread species showed strong genetic cohesion and clear separation from congeners. Our results indicate that the sampled Polylepis species from Bolivia and Ecuador maintain distinct genetic identities despite localized signals consistent with gene flow. This genome-wide support for current taxonomy highlights Polylepis as a valuable model for studying speciation under gene flow and indicates that multiple geographic sampling will be essential in reconstructing a robust phylogeny of the genus, with important implications for conservation planning in Andean montane forests.

Bolivia↗
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↗