SEARCH · PubMed Health

Results for “genome-wide analysis”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 19 recordsLinked to original sources

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

Genome-wide characterization of heat shock protein genes reveals thermal stress-responsive candidates in Litopenaeus vannamei.

Heat shock proteins (HSPs) are conserved molecular chaperones involved in protein folding, refolding, aggregation prevention, and degradation of damaged proteins. However, the genomic organization and thermal responsiveness of HSP genes in the Pacific white shrimp (Litopenaeus vannamei) remain incompletely understood. Here, we performed a genome-wide analysis of the HSP gene family and examined its phylogenetic relationships, structural features, duplication patterns, sequence variation, interaction networks, and transcriptional responses to acute heat stress. A total of 34 HSP genes were identified and classified into the HSP90, HSP70, HSP40/DNAJ, HSP60, and small HSP families. Phylogenetic, motif, gene structure, synteny, and subcellular localization analyses revealed evolutionary conservation and structural diversification among family members. Three duplicated gene pairs were identified, comprising two segmental duplications and one tandem duplication. All pairs exhibited Ka/Ks ratios below 1, consistent with purifying selection of varying strength. Sequence analysis identified 295 nonsynonymous single-nucleotide polymorphisms, of which 12 were consistently predicted to be deleterious by multiple algorithms. Protein-protein interaction analysis indicated enrichment of protein-folding and cellular stress-response functions. RT-qPCR analysis showed significant induction of HSPA4, HSP90AA1, TRAP1, BiP, and DNAJA1 after 6, 12, and 24 h of exposure to 34 °C, whereas DNAJC3 was significantly induced only at 12 h. All six genes reached their highest transcript abundance at 12 h. These findings may provide a genomic framework for HSP genes in L. vannamei and identify candidate genes and variants associated with thermal stress responses.

Animals↗

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

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

Peripheral pain threshold, glycaemic status, and LAMP3 genetic variation: A community-based analysis.

Diabetic polyneuropathy is a common complication of diabetes, yet substantial inter-individual variation in peripheral pain perception suggests underlying genetic influences. This population-based study investigated clinical, metabolic, and genetic determinants of pain threshold using intraepidermal electrical stimulation in 906 participants from the Iwaki Health Promotion Project 2017. Genome-wide association analysis identified 12 loci showing suggestive associations, among which a missense variant in LAMP3 (rs482912) was prioritized as a biologically plausible candidate. Phenotype-stratified analyses showed that individuals carrying the CT or CC genotypes had lower PINT indices than those with the TT genotype, indicating reduced pain thresholds. Notably, the CC genotype retained an association with lower pain threshold using intraepidermal electrical stimulation under conditions of metabolic stress, including impaired glucose tolerance, elevated HbA1c, and obesity, whereas this association was attenuated in the presence of hypertension. Single-cell RNA sequencing analysis of human skin revealed that LAMP3-positive mature dendritic cells, enriched in immunoregulatory molecules, exhibited transcriptional enrichment of inflammatory, antigen-presenting, and nociception-related pathways, including NF-κB, JAK-STAT, cytokine signaling, and neuroimmune sensitization cascades. Autopsy-based skin analysis further demonstrated genotype-associated differences in dermal LAMP3-positive cell infiltration and CD8-positive T-cell abundance, while CD4-positive T-cell abundance and intraepidermal nerve fiber density remained unchanged across genotypes. Taken together, these findings suggest a potential association between LAMP3 variation and individual differences in peripheral pain threshold and provide biological context supporting a role for neuroimmune interactions in early sensory modulation under metabolic stress. Given the suggestive genetic evidence and indirect mechanistic data, these observations should be interpreted as exploratory and hypothesis-generating.

Humans↗

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, structural characterization, and evolutionary analysis of growth-related gene families in African catfish (Clarias gariepinus).

The somatotropic axis encompassing growth hormone (GH), insulin-like growth factor (IGF), myostatin (MSTN), and prolactin (PRL) signalling cascades is the master regulator of somatic growth, metabolism, and development in vertebrates. African catfish (Clarias gariepinus), a commercially pivotal aquaculture species, now possesses a chromosome-level reference genome (CGAR_prim_01v2); however, a systematic, genome-wide characterization spanning all five interconnected growth-related gene families has not previously been undertaken in this species. Here, we identified and characterized 15 growth-related genes spanning gh1, ghra, ghrb, Igf1, Igf2a, Igf2b, igf1ra, Igf1rb, Igf2r, Mstna, Mstnb, prl, prlra, prlrb, and smtlb distributed across 13 chromosomes. Complete one-to-one orthology with zebrafish confirmed strong dosage-balance conservation across >120 million years of teleost divergence. Physicochemical analysis resolved a clear biochemical dichotomy between compact, basic secreted ligands (19.88-45.81 kDa; pI up to 10.02) and large, acidic, heavily glycosylated membrane receptors (56.82-270.80 kDa; pI 4.85-5.97). Phylogenetic analysis confirmed 3R whole-genome duplication origins for all paralog pairs, while synteny analysis revealed a disruption of the ancestral gh1-prl chromosomal block in C. gariepinus, a finding that warrants further comparative and functional investigation. This genomic atlas provides the sequence and structural information including exon-intron boundaries, domain architecture, and chromosomal coordinates needed as a prerequisite for future marker-assisted selection and CRISPR-based myostatin-editing efforts in African catfish aquaculture, though translation into applied breeding outcomes will require subsequent functional and expression studies.

Animals↗

PubMed2026

Identification and characterization of G protein-coupled receptors in the nocturnal halictid bee Megalopta genalis.

G protein-coupled receptors (GPCRs) are one of the largest families of membrane proteins in insects, regulating vision, neural signal transduction, and various physiological behaviors. Megalopta genalis exhibits a unique facultatively eusocial lifestyle and possesses adaptations for nocturnal activity; however, its GPCR family has not yet been systematically characterized. In this study, we performed genome-wide identification, phylogenetic analysis, and expression profiling of GPCRs in M. genalis by integrating genomic annotation and transcriptomic analysis. The results showed that a total of 99 GPCRs were identified in the genome of M. genalis, which were classified into four major families. Here, we show that M. genalis has undergone lineage-specific GPCR repertoire remodeling, marked by the expansion of novel orphan receptors and the systematic loss of multiple receptor subtypes, such as the neuropeptide receptors MIP-R and NPFR. Moreover, opsins have formed a diverse array of combinations and non-GPCR odorant receptors have undergone significant expansion via tandem duplication. Together, these features may represent part of the molecular repertoire associated with the adaptation of M. genalis to a nocturnal lifestyle. Furthermore, transcriptomic analysis revealed distinct spatiotemporal expression divergence within each of the Mth/Mthl and Fz GPCR families, suggesting functional specialization across development and adult tissues. This study provides the first systematic identification and initial functional characterization of GPCRs in M. genalis, revealing an evolutionary pattern characterized by the coexistence of contraction and expansion within the GPCR family. These findings lay a foundation for further studies aimed at elucidating the roles of these GPCRs in regulating M. genalis physiology and behavior.

Animals↗

PubMed2026

Genome-wide scans reveal candidate genes associated with wing morph differentiation in Tetrix japonica.

Wing dimorphism is an important dispersal-related trait in insects, but its genomic basis remains poorly understood in pygmy grasshoppers. Here, we integrated genome-wide single-nucleotide polymorphism (SNP) analyses, population structure inference, selection scans, and functional annotation to investigate genomic differentiation between long- and short-winged Tetrix japonica. Principal component analysis (PCA), ADMIXTURE, and phylogenetic analyses revealed weak genome-wide separation between morphs, indicating differentiation on a largely shared genetic background. Genome-wide scans based on the fixation index (FST), nucleotide diversity ratios, and Tajima's D, using 50-kb non-overlapping windows and empirical top-5% outlier thresholds, identified multiple candidate regions across seven chromosomes. The broader long- and short-winged candidate sets spanned 9.35 Mb and 9.37 Mb and directly overlapped 82 and 77 genes, respectively. Candidate genes were associated with signaling/hormone regulation, membrane transport, metabolism, cytoskeletal organization, extracellular matrix structure, and development. Short-winged candidate genes were significantly enriched for ABC-type transporter activity and ATP hydrolysis activity. Because all individuals originated from a single laboratory-maintained population with weak genome-wide structure, these regions should be regarded as candidate loci from a screening-stage analysis that require validation in independent populations and by functional assays, rather than as confirmed targets of selection.

Animals↗

PubMed2026

Retinal microstructural alterations as early phenotypes of depression in radiogenomics analysis.

BACKGROUND: With the increasing prevalence of depression, there is an urgent clinical need for early screening in depression. The retina offers a promising window for early screening in depression due to its rapid, non-invasive, objective, eye-brain correlated characteristics, but previous research has yielded conflicting alterations in retinal microstructure in depression. METHODS: We screened retinal optical coherence tomography and brain magnetic resonance imaging data in the UK Biobank to enroll 23,225 participants for retinal study of depression occurrence, and 1475 participants for the eye-brain association study. We also used genetic data (ID: ebi-a-GCST90014267 and ukb-d-20,448) from the Integrative Epidemiology Unit Open Genome-Wide Association Study for Mendelian randomization analysis. We used Cox regression to assess the association between retinal microstructure and depression risk, Mendelian randomization to infer causality, and mediation analysis to explore retina-brain pathway association. RESULTS: The Cox regression analysis showed that retinal ganglion cell-inner plexiform layer (GCIPL) thickness remained a significant predictor of depression. The Mendelian randomization analysis indicated a positive statistical association between GCIPL thickness and depression. Moreover, there was a significant positive correlation (all p < 0.001) between the volume of specific depression-related brain regions and the GCIPL thickness. Adjusting for age, sex, and head size, the mediation analysis provided preliminary evidence for a potential anatomical pathway linking retinal GCIPL thickness to depression-related brain regions through primary visual cortex and secondary visual cortex volumes. CONCLUSION: Thickened retinal GCIPL is a potential early phenotype of depression and has a potential association pathway with depression-related brain regions using a radiogenomics approach.

Humans↗

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

Integrative multi-omics reveals a fibroblast-centered, ZFHX3-prioritized regulatory framework linking sick sinus syndrome and atrial fibrillation.

OBJECTIVE: To define shared genetic and multi-scale mechanisms underlying comorbidity between sick sinus syndrome (SSS) and atrial fibrillation (AF). METHODS: We integrated genome-wide association study (GWAS) summary statistics for SSS and AF with Genotype-Tissue Expression (GTEx) expression and splicing quantitative trait loci (eQTL/sQTL), atrial single-cell and spatial transcriptomics, and epigenomics. We identified trait-relevant tissues and pathways, prioritized shared cell types, quantified genome-wide and local genetic sharing, detected joint loci by cross-trait meta-analysis, and linked loci to regulatory programs via colocalization and cell-prioritized co-expression networks. RESULTS: Both traits showed strongest enrichment in cardiac tissue, especially Heart Atrial Appendage. Fibroblasts from the left atrial appendage were consistently prioritized as the key shared cell population. SSS and AF displayed significant positive genome-wide genetic correlation, with multiple locally shared regions, including six major loci. Cross-trait meta-analysis identified eight joint-phenotype SNPs implicating four susceptibility genes. ZFHX3 was the leading tissue-cell-gene candidate, acting as a hub in fibroblast co-expression modules and colocalizing with cardiac regulatory signals. CONCLUSION: Shared liability for SSS and AF is highly tissue- and cell-specific, converging on regulatory networks in atrial appendage fibroblasts, with ZFHX3 serving as a central mechanistic and biomarker node.

Humans↗

PubMed2026

Genetic evidence for a causal relationship between melatonin metabolism and depression.

To investigate the causal relevance of melatonin metabolism, which provides the biological basis for circulating melatonin levels, to specific depression symptom subtypes, we performed a targeted systematic review of melatonin metabolism pathways in the human brain and liver. Using two-sample Mendelian randomization (MR), we assessed the causal effects of metabolism pathways and/or individual genes on major depressive disorder (MDD) and nine symptom subtypes derived from Patient Health Questionnaire-9 (PHQ-9). Instrumental variables (IVs) were expression quantitative trait loci (eQTL) for eight individual genes, one synthesis route, and three degradation routes. Results were assessed using Bayesian colocalization and phenome-wide association analyses. At the pathway-level, the genetically proxied synthesis-route signal was associated with PHQ-9 Assessment 5 (PHQ9A5, OR: 0.89, 95% CI: 0.85-0.93), but sensitivity analyses suggested this association was primarily driven by TPH1 and may reflect serotonin-related biology. In contrast, higher brain melatonin degradation raised the risk of both PHQ9A1 (OR: 1.03, 95% CI: 1.02-1.04) and PHQ9A7 (OR: 1.03, 95% CI: 1.02-1.03). Within degradation, up-regulation of the kynurenine sub-pathway increased the odds of PHQ9A3 (OR: 1.05, 95% CI: 1.02-1.07), PHQ9A4 (OR = 1.04, 95% CI: 1.02-1.06) and PHQ9A7 (OR: 1.05, 95% CI: 1.02-1.07). Gene-level analyses were largely concordant, except for SULT1A1, whose higher expression was genetically protective for PHQ9A3 but risk-increased for PHQ9A1 and PHQ9A4. Overall, these results demonstrate that melatonin metabolism exerts symptom-specific and pathway-specific causal effects on depression. A stratified view of melatonin's role may help optimize the application of exogenous melatonin supplementation.

Melatonin↗

PubMed2026

Prioritizing Parkinson's disease risk-associated mitochondrial candidate genes via multi-omics integrative analysis.

BACKGROUND: Mitochondrial dysfunction has been implicated in Parkinson's disease (PD), but the genetically regulated mitochondrial genes associated with PD risk remain incompletely defined. METHODS: We conducted a summary-data-based genetic epidemiology study integrating summary-based Mendelian randomization (SMR), Heterogeneity in dependent instruments (HEIDI) filtering, and Bayesian colocalization to prioritize mitochondrial-related molecular features associated with PD risk. Mitochondrial-related genes were defined using MitoCarta3.0. Genetically predicted gene expression and plasma protein abundance were evaluated using expression quantitative trait loci (eQTL) data from eQTLGen and GTEx v8, and protein quantitative trait loci (pQTL) data was assessed using International Parkinson's Disease Genomics Consortium (IPDGC) as the discovery genome-wide association study (GWAS) and FinnGen as the replication dataset. Prespecified QTL analyses were interpreted using FDR correction, HEIDI filtering, and colocalization support. DNA methylation QTL analysis, mitochondrial phenotype MR, and single-nucleus RNA-seq analysis were performed as complementary analyses. RESULTS: In the primary eQTL analysis, higher genetically predicted TTC19 expression was associated with lower PD risk (OR = 0.80, 95% CI: 0.74-0.87, PPH4 = 0.80), whereas higher MALSU1 expression was associated with increased PD risk (OR = 2.21, 95% CI: 1.59-3.06, PPH4 = 0.96). Both associations survived FDR correction, passed HEIDI filtering, and showed colocalization support. GTEx whole-blood data supported the direction of the TTC19 association. No mitochondrial protein reached significance after FDR correction and colocalization filtering in the primary pQTL analysis. Complementary methylation analysis highlighted cg06270993 as an exploratory regulatory signal for MALSU1. CONCLUSIONS: This MR-colocalization study prioritizes TTC19 and MALSU1 as genetically supported mitochondrial-related candidate genes associated with PD risk. Further validation is required to define their functional roles in PD pathogenesis.

Humans↗

PubMed2026

Analyzing salinity tolerance in grass carp (Ctenopharyngodon idella): Insights from genome-wide association study and genomic selection.

Grass carp (Ctenopharyngodon idella) is one of the most widely cultured freshwater fish species globally. However, the expansion of its farming scale faces severe limitation owing to freshwater scarcity; therefore, the development of strains with greater salinity tolerance is key for expanding production using brackish water resources. To investigate the genetic basis of salinity tolerance in grass carp, a genome-wide association study (GWAS) was conducted using 200 individuals representing extreme phenotypes, namely salinity-tolerant and salinity-sensitive groups. In total, 17 single nucleotide polymorphisms (SNPs) related to salinity tolerance were detected, which were distributed across 11 chromosomes. Through gene annotation, 38 candidate genes were obtained from these loci. Enrichment analysis revealed these candidate genes are primarily implicated in key biological processes, including osmotic regulation, energy metabolism, and stress responses. Analyses of different SNP densities revealed that the 5 K SNP density panel can balance prediction accuracy and computational efficiency. The BayesA model achieved the highest prediction accuracy under the GWAS_Evenly selection strategy, with substantial reductions in mean absolute error and mean square error. This study reveals the genetic mechanisms of salinity tolerance in grass carp, which might be optimized through genomic selection, and provides insights for selectively breeding new varieties with greater salinity tolerance.

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

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