SEARCH · Search PubMed

Results for “Genome-Wide Association Study”

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.

643 recordsLinked to original sources

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

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

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

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

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

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

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

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

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 analysis of the plant-specific PLATZ gene family in Taraxacum kok-saghyz and its roles in response to drought and salt tolerance.

Abiotic stress severely limits plant growth and productivity. Taraxacum kok-saghyz Rodin (TKS), known for its environmental resilience, represents a valuable resource for identifying stress-tolerant genes to improve stress-adaptive crops. Plant AT-rich protein and zinc-binding protein (PLATZ) transcription factors serve as core regulators of plant growth, developmental processes, and adaptive responses to various stress conditions; however, they remain uncharacterized in TKS. Here, we identified 10 TksPLATZ genes through a whole-genome analysis. Phylogenetically, these genes were grouped into five distinct evolutionary branches. Promoter sequence analysis revealed multiple types of cis-acting regulatory elements that are connected with hormonal signal responses and environmental stress adaptation. Integrated analysis of transcriptome datasets and RT-qPCR validation demonstrated that TksPLATZ genes display tissue-specific expression profiles and show distinct responsive patterns to drought and salt stress treatments. Among them, TksPLATZ1, TksPLATZ2 and TksPLATZ7 were markedly induced under both stressors and were selected for further functional study. We demonstrated that TksPLATZ1, TksPLATZ2 and TksPLATZ7 localize to the cell nucleus and act as transcriptional activators and repressors, respectively. Phenotypic data from overexpression experiments in plants confirm that heterologous expression of TksPLATZ1, TksPLATZ2, and TksPLATZ7 enhances the tolerance of Arabidopsis to salt and osmotic stress. These findings provide valuable genetic resources for improving plant tolerance to environmental stresses.

Salt Tolerance↗

PubMed2026

Genome-wide insights into the evolutionary and demographic history of the red alga Mazzaella laminarioides: Evidence for speciation with ancient migration along the southeast Pacific coast.

The mechanisms driving lineage divergence in red algae remain unexplored, despite the group's remarkable diversity and ancient evolutionary history. The red alga Mazzaella laminarioides, a Chilean intertidal species complex composed of three parapatric cryptic lineages (North, Center, South), offers a valuable system to evaluate these processes, as its life history combines severe dispersal limitation with a haploid-diploid cycle that may influence the emergence of reproductive barriers. We reconstructed its evolutionary history using whole-genome sequencing and nuclear genome assembly of representative individuals from each lineage. Phylogenomic analyses based on 1,507 single-copy orthologs recovered three deeply divergent lineages with limited nuclear discordance consistent with incomplete lineage sorting. For both splits, demographic modelling was most consistent with an Ancient Migration scenario, although support over strict isolation was moderate, suggesting that divergence may have begun with low asymmetric ancestral gene flow followed by subsequent loss of connectivity, demographic bottlenecks, and later population expansion. Coding sequence analyses revealed lineage-specific dN/dS heterogeneity; only one South-lineage locus passed FDR correction (metaxin-1, mitochondrial protein import), with two further South-lineage candidates in chlorophyll and heme biosynthesis falling below the FDR threshold. Together, these signals suggest that divergent selective pressures on energy acquisition may have contributed to divergence at the southern end of the distribution. These results add to the small but growing body of whole-genome data for red algae and, alongside recent macroalgal studies, suggest that ancestral connectivity could be a recurrent feature of lineage divergence even in marine organisms with extremely restricted dispersal.

Rhodophyta↗

PubMed2026

Cross-tissue multi-omics integration highlights BPHL and mitochondrial targets in Alzheimer's disease.

BACKGROUND: Mitochondrial dysfunction is a hallmark of Alzheimer's disease (AD), yet specific molecular targets remain to be fully characterized. METHODS: A summary-data-based Mendelian randomization (SMR) framework integrated AD genome-wide association study (GWAS) statistics (39,918 cases) with blood DNA methylation quantitative trait loci (mQTL), gene expression (eQTL), and protein (pQTL) data for 1136 mitochondria-related genes. Associations were assessed using Bayesian colocalization and HEIDI testing. Tissue relevance was evaluated in four brain regions (hippocampus, amygdala, cortex, frontal cortex) using GTEx and external transcriptomic datasets. RESULTS: Screening identified eight candidates supported across blood mQTL and eQTL layers. Stepwise central nervous system (CNS) evaluation singled out biphenyl hydrolase-like (BPHL) as the consistent candidate. Higher genetically predicted BPHL expression was associated with reduced AD risk across the hippocampus (OR=0.920, 95% CI 0.873-0.970), amygdala (OR=0.925, 95%CI 0.880-0.973), cortex (OR=0.943, 95% CI 0.908-0.978), and frontal cortex (OR=0.938, 95%CI 0.901-0.976). These findings aligned with protein-protein interactions connecting BPHL to respiratory complexes and lower BPHL expression in independent AD brains. Functional enrichment converged on oxidative phosphorylation pathways. CONCLUSIONS: By integrating multi-omics data with tissue-specific validation, this study nominates BPHL as a consistent protective candidate in the brain. These findings provide genetic support for mitochondrial molecular perturbations in AD, offering insights for future validation.

Alzheimer Disease↗

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

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

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

Stage-specific ROMO1 in rheumatoid arthritis: predictive immune insights into the MIF pathway and HLA-DR/IL2RA axis via integrated GWAS, transcriptomic, single-cell, and spatial profiling.

Emerging evidence links reactive oxygen species modulator 1 (ROMO1), a key mitochondrial ROS regulator, to rheumatoid arthritis (RA) pathogenesis. However, its exact mechanism remains elusive given the conflicting evidence about its specific function. We used a four-level integrative framework combining multi-omics data and literature‑supported mechanistic inference. At the genetic level, Mendelian randomization (MR) was performed to explore potential causal relationships between ROMO1, IL2RA, HLA-DR, MIF, and RA risk, followed by differential expression analysis and machine learning-based feature selection to identify key mROS genes. The temporal expression dynamics of ROMO1 were assessed in RA progression. At the cellular and tissue levels, we integrated single-cell RNA sequencing and spatial transcriptomics to map cell-type-specific expression and synovial localization of ROMO1-related immune cells and pathways. Finally, our multi-omics findings were contextualized with literature-supported mechanistic inference. (1) MR results were consistent with a potential protective effect of ROMO1 on RA (OR = 0.52) and its potential regulation of risk factors IL2RA (OR = 0.46) and HLA-DR (OR = 0.40). Conversely, IL2RA (OR = 1.42), HLA-DR (OR = 1.88), and MIF (OR = 1.17) were positively associated with RA risk. Additionally, ROMO1 was identified as a top candidate diagnostic predictor with stage-specific dynamics: downregulated in the early but upregulated in the late/remission stages. (2) Single-cell RNA sequencing showed ROMO1's cell-specific expression in CD14+ HLA-DR+ CD74+ monocytes and CD4+ IL2RA+ T cells. Cell communication analysis further suggested that these cells may participate in MIF pathway regulation. Spatial transcriptomics subsequently identified that ROMO1-related cells localized to synovial pathological regions, with MIF pathway changes correlated with RA progression. (3) Finally, literature-supported mechanistic inference suggests that ROMO1 may modulate mROS levels to promote anti-inflammatory M2 macrophage polarization, which could theoretically contribute to reduced systemic inflammation and the alleviation of multi-organ decline in RA. This integrated multi-omics investigation, supported by literature-based mechanistic inference, suggests ROMO1 as a stage-dependent biomarker candidate and potential immune regulator in RA.

Humans↗

PubMed2026

Introgression shapes the genomic conflict landscape of Malus, providing evidence for a reticulate backbone in a woody crop lineage.

Phylogenomic discordance is widespread across plants, but its evolutionary significance is often obscured when conflict is treated primarily as analytical noise rather than as evidence of underlying processes. In woody lineages in particular, incomplete lineage sorting, introgression, and genome duplication can interact over long timescales to produce complex genomic histories that are not adequately summarized by a strictly bifurcating tree. Here, we use Malus as a model woody genus to investigate how these processes structure conflict across a genus-scale, accession-based phylogenomic framework. Using broad taxon sampling, hundreds of nuclear loci, plastid genomes, and genome-wide SNP summaries, we reconstruct a robust nuclear backbone for sampled Malus lineages and evaluate where discordance is concentrated and which processes best explain it. Nuclear analyses resolve eight major clades, whereas conflict is non-random and localized to recurrent hotspots rather than evenly distributed across the tree. Cytonuclear discordance is similarly concentrated, especially around Clade H, represented by sampled accessions of M. tschonoskii, where localized plastid-nuclear disagreement is consistent with candidate plastid capture or organellar introgression. Multiple complementary analyses further indicate that the strongest conflict is not explained by ILS alone, but instead reflects lineage-structured introgression, while polyploid complexes represent additional localized sources of evolutionary complexity. Together, these results provide evidence for a reticulate genomic backbone in Malus and show how integrating nuclear, plastid, and genome-wide conflict analyses can help distinguish background discordance from process-specific signals in woody plant radiations. Several lineage-level reticulation hypotheses identified here should now be tested with broader population-level sampling and curated reference accessions.

Malus↗