Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “QTLs”

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.

At least 37 records · Page 2Linked to original sources
PubMed2026

Integrated Genome-Wide Association Studies and Selection Signature Analysis Reveal the Genetic Architecture of the Wattle Trait in Goats.

Wattles are finger-like appendages on the ventral neck of goats, serving as a distinctive morphological marker for breed identification that serves potential implications for production performance. However, their genetic basis remains incompletely characterized. Here, we integrated a genome-wide association study (GWAS) and selection signature analysis to identify candidate genes and genomic regions associated with the wattle trait in goats. Using a linear mixed model, GWAS on 463 goats (23 wattled and 440 non-wattled) identified 385 quantitative trait loci (QTLs) at a 5% false discovery rate, yielding 346 candidate genes. The most significant association signal was observed on chromosome 10 (72.61-73.48 Mb), where the lead SNP (rs636481767) is located within a region containing GJD2, GREM1, and FMN1, showing strong linkage disequilibrium (r2 > 0.6) with surrounding loci. Subsequent selection signature analysis (23 wattled and 23 non-wattled) identified 83 genomic regions harboring 119 candidate genes. The strongest signals were detected at MFSD14B on chromosome 8 (FST = 0.154, log2π-ratio = 2.611) and PDLIM7 on chromosome 7 (FST = 0.144, log2π-ratio = 0.806). KEGG pathway enrichment analysis revealed that GWAS-associated genes were involved in glycosylation and immune responses, whereas selection-signature genes were enriched in DNA repair and the Hippo, Notch, and Wnt pathways. Furthermore, cross-species PheWAS revealed that human FMN1 is associated with dermatological, skeletal, and metabolic phenotypes, while porcine FMN1 is associated with backfat thickness and loin muscle depth. Overall, this study provides molecular markers of potential value for goat breeding and pinpoints key candidate genes for future functional validation of wattle development.

Animals↗
PubMed2026

Cross-tissue Mendelian randomization prioritizes RAB27B as a brain-derived candidate protein for postpartum depression.

OBJECTIVE: Postpartum depression (PPD) is one of the most common and debilitating complications of childbirth, yet the candidate proteins linking genetic risk to disease remain poorly defined. Building on recent genome-wide association studies (GWAS), we sought to integrate cross-tissue proteogenomic data to identify candidate proteins for PPD and explore therapeutic opportunities. METHODS: We conducted two-sample Mendelian randomization (MR) using genome-wide significant cis-protein QTLs from brain (n = 608 proteins), cerebrospinal fluid (CSF; n = 214), and plasma (n = 612). PPD summary statistics were obtained from FinnGen R8 (13,657 cases, 236,178 controls) and replicated in an independent GWAS. Phenome-wide association (PheWAS) was used to assess pleiotropy. Potential therapeutic targets were evaluated through DSigDB drug repurposing, molecular docking, and molecular dynamics simulations. RESULTS: Among all proteins tested, RAB27B was the only brain-derived protein surpassing Bonferroni correction (OR = 1.60; 95% CI: 1.30-1.96; P = 6.6 × 10⁻⁶), whereas no significant proteins were identified in CSF or plasma. This association was replicated in an independent GWAS (OR = 1.27; 95% CI: 1.02-1.58; P = 0.037). PheWAS identified no pleiotropic associations at genome-wide significance. In silico drug repurposing identified pregnenolone as a candidate ligand with computationally predicted stable binding to RAB27B, providing a starting point for future experimental validation. CONCLUSION: This study provides the first cross-tissue proteogenomic evidence that RAB27B is a brain-derived, reproducible candidate protein genetically associated with PPD. By extending GWAS signals to functional protein-level mechanisms and therapeutic inference, our findings nominate RAB27B and pregnenolone as promising directions for postpartum psychiatric research.

Humans↗
PubMed2026

Mitochondrial Function-Related Genes in Sleep Disorders: A Multi-Omics Mendelian Randomization Study.

Mitochondrial dysfunction is linked to sleep disorders in previous report, but the potential roles of specific genes remain unclear. This study aimed to dissect different subtype-specific genetic associations and their underlying mechanisms. A multi-omics Summary-data-based Mendelian Randomization (SMR) approach was performed to identify potential causal links between mitochondrial function-related genes and sleep disorders. We integrated GWAS data from FinnGen database (the discovery set), independent GWAS datasets (covering different sleep-disorder subtypes and used for validation), and cis-QTLs (including mQTLs, eQTLs, and pQTLs) to perform systematic exploration. Specially, we performed targeted validation of tissue-specific effects, leveraging gene expression data from disease-relevant brain regions within the GTEx database. Our SMR analysis identified mitochondrial function-related genes potentially modulating sleep disorders across biological layers, initially identifying 102 genes at the methylation level, 48 at the gene expression level, and 6 at the protein abundance level. Integrative analysis subsequently prioritized DCXR and ACADVL and revealed their distinct, subtype-specific associations. DCXR exhibited a protective role in sleep apnea while ACADVL showed a paradoxical risk conferring role in daytime sleepiness. In addition, the analysis identified an epigenetic regulatory mechanism for DCXR in which its expression and protein levels are modulated by DNA methylation. Finally, validation in brain-hypothalamus tissue confirmed DCXR as a significant potential protective factor (OR = 0.929, 95% CI: 0.887-0.973, P_HEIDI = 0.999, FDR = 0.2449). Our findings implicate key mitochondrial genes, particularly DCXR and ACADVL, in the pathophysiology of specific sleep disorder subtypes, highlighting potential avenues for precision medicine. Clinical trial number: Not applicable.

Humans↗
PubMed2026

Multi-omics identification of therapeutic targets of compound sappan decoction in hepatocellular carcinoma.

BACKGROUND: Compound sappan decoction (CSD) is a multi-herbal traditional Chinese medicine formulation with clinical relevance in hepatocellular carcinoma (HCC). However, its therapeutic mechanisms remain unclear. METHODS: Bioactive compounds of CSD were identified and standardized using pharmacological and chemical databases. Potential targets were predicted via multiple target inference platforms. HCC-related genes were curated from comprehensive disease databases. Summary-data-based Mendelian randomization (SMR) was conducted to infer causal relationships between compound targets and HCC risk using large-scale quantitative trait loci (QTL) datasets and HCC genome-wide association study data. Colocalization analysis, protein-protein interaction (PPI) network construction, and GO/KEGG enrichment were performed on SMR-identified targets. Molecular docking evaluated binding affinities of representative compounds to prioritized targets. RESULTS: A total of 784 overlapping genes between predicted CSD targets and HCC-related genes were subjected to SMR analysis. Among these, 22 targets were significantly associated with HCC risk based on transcriptomic or proteomic QTLs and showed colocalization evidence. Notably, four targets (ADRB2, APOE, SYK, and PGF) were supported by both replication in an independent cohort and strong colocalization. These 22 targets were enriched in apoptosis, PI3K-Akt signaling, redox metabolism, and detoxification pathways. PPI analysis revealed central hubs including MMP9, BCL2, CASP1, and MCL1. Molecular docking demonstrated strong binding of APOE to quercetin, PGF to luteolin-7-olate, and SYK to kaempferol. CONCLUSIONS: CSD may exert therapeutic effects on HCC through modulation of genetically validated targets involved in tumor progression, inflammation, and metabolic reprogramming, supporting its potential clinical utility as an adjunctive treatment strategy. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at https://doi.org/10.1007/s12672-026-04740-8.

Caesalpinia↗
PubMed2026

Genomic diversity, inbreeding, and selection signatures in duroc, landrace, and yorkshire pigs from a long-term closed breeding system.

Duroc (DD), Landrace (LL), and Yorkshire (YY) are among the most widely used commercial pig breeds, having undergone intense long-term selection within closed breeding systems. This study presents a comprehensive genomic analysis of genetic diversity, inbreeding patterns, and selection signatures in DD, LL, and YY populations that have been subject to close breeding for over 15 years. Genomic and pedigree data were available for 1,088 animals (DD = 348, LL = 276, YY = 464), genotyped using the GenoBaits® Porcine 100 K SNP panel. Principal component analysis and genetic diversity metrics revealed distinct population structures among the three breeds. Pairwise genetic differentiation supported this pattern, with DD showing the greatest divergence from LL (0.34 ± 0.24) and YY (0.33 ± 0.24), while LL and YY were more closely related (FST = 0.22 ± 0.19). Linkage disequilibrium (LD) analysis further confirmed these differences, as DD exhibited the highest average r² (0.34), followed by LL (0.28) and YY (0.25). Within-breed genetic diversity metrics, including observed heterozygosity (HO: 0.37 in DD, 0.39 in LL, 0.38 in YY), expected heterozygosity (HE: 0.36 in DD, 0.37 in LL, 0.38 in YY), and minor allele frequency (MAF: 0.27 in DD, 0.28 in LL, 0.29 in YY), indicated greater genetic variability in LL and YY compared to DD. Runs of homozygosity (ROH) analyses revealed different patterns of autozygosity, with DD exhibiting more long ROH indicative of recent inbreeding, while YY harbored a higher number of short ROH, suggestive of more ancient demographic events. ROH-based inbreeding coefficients (FROH) consistently exceeded pedigree-based estimates (FPED) across all breeds, highlighting the presence of recent or unrecorded inbreeding that pedigree data may not fully capture. According to Generation Proxy Selection Mapping (GPSM), 17, 1, and 12 significant SNPs were detected in DD, LL, and YY, respectively. Functional annotation of ROH islands and GPSM-significant loci revealed both breed-specific and overlapping QTLs related to traits such as growth, reproduction, and carcass. In general, the findings of this study contribute to a deeper understanding of the genomic consequences of long-term closed breeding and provide reference information to support consideration of breeding strategies that balance continued selection for productivity with the maintenance of genetic diversity in modern commercial pig populations.

Animals↗
PubMed2024

Crosstalk between epitranscriptomic and epigenomic modifications and its implication in human diseases.

Crosstalk between N6-methyladenosine (m6A) and epigenomes is crucial for gene regulation, but its regulatory directionality and disease significance remain unclear. Here, we utilize quantitative trait loci (QTLs) as genetic instruments to delineate directional maps of crosstalk between m6A and two epigenomic traits, DNA methylation (DNAme) and H3K27ac. We identify 47 m6A-to-H3K27ac and 4,733 m6A-to-DNAme and, in the reverse direction, 106 H3K27ac-to-m6A and 61,775 DNAme-to-m6A regulatory loci, with differential genomic location preference observed for different regulatory directions. Integrating these maps with complex diseases, we prioritize 20 genome-wide association study (GWAS) loci for neuroticism, depression, and narcolepsy in brain; 1,767 variants for asthma and expiratory flow traits in lung; and 249 for coronary artery disease, blood pressure, and pulse rate in muscle. This study establishes disease regulatory paths, such as rs3768410-DNAme-m6A-asthma and rs56104944-m6A-DNAme-hypertension, uncovering locus-specific crosstalk between m6A and epigenomic layers and offering insights into regulatory circuits underlying human diseases.

Humans↗
PubMed2025

Identification of novel type 1 and type 2 diabetes genes by co-localization of human islet eQTL and GWAS variants with colocRedRibbon.

Over 1,000 genetic variants have been associated with diabetes by genome-wide association studies (GWASs), but for most, their functional impact is unknown; only 7% alter gene expression in pancreatic islets in expression quantitative trait locus (eQTL) studies. To fill this gap, we developed a co-localization pipeline, colocRedRibbon, that prefilters eQTLs by the direction of effect on gene expression and shortlists overlapping eQTL and GWAS variants prior to co-localization. Applying colocRedRibbon to recent diabetes and glycemic trait GWASs, we identified 292 co-localizing gene regions, including 24 co-localizations for type 1 diabetes and 268 for type 2 diabetes and glycemic traits, representing a 4-fold increase. A low-frequency type 2 diabetes protective variant increases islet MYO5C expression, and a type 1 diabetes protective variant increases FUT2 expression. These novel co-localizations advance the understanding of diabetes genetics and its impact on human islet biology. colocRedRibbon has broad applicability to co-localize GWASs and various QTLs.

Humans↗
PubMed2026

Integrating GWAS and Transcriptome Analysis Identifies Candidate Genes for Kernel Starch Quality Traits in Maize.

Maize (Zea mays L.) starch quality is a complex trait with significant implications for grain processing and industrial applications. However, the genetic basis underlying starch quality, particularly for gelatinization and thermodynamic properties, remains poorly understood. In this study, we evaluated 12 starch quality traits, including seven gelatinization characteristics, four thermodynamic traits, and kernel starch content (KSC) in a diverse panel of 335 maize inbred lines. Considerable phenotypic variation was observed for all traits. A total of 228 quantitative trait loci (QTLs) were significantly associated with 12 starch quality traits through genome-wide association studies (GWAS). By integrating a dynamic transcriptome analysis of two maize inbred lines with contrasting starch quality, we identified 60 candidate genes. One gene, waxy1, encoding a starch synthase, was found to be associated with enthalpy of gelatinization (ΔHgel) and pasting temperature (Ptemp). Six variants in waxy1 contributed to natural variation in ΔHgel and Ptemp, and a cost-effective InDel and two PARMS-based molecular markers were developed and validated in 144 maize inbred lines, enabling efficient marker-assisted selection. Our findings provide key genes and molecular markers for high-quality maize breeding with improved starch properties.

Zea mays↗
PubMed2025

Integrative omics of the genetic basis for wheat WUE and drought resilience reveal the function of TaMYB7-A1.

Improving wheat drought resilience and water use efficiency (WUE) is critical for sustaining productivity under increasing water scarcity. Here, we integrate genome-wide association study (GWAS), expression quantitative trait locus (eQTL) mapping, population-transcriptome analysis, and summary-data-based mendelian randomization (SMR), followed by functional validation using indexed EMS mutants and transgenic lines, to systematically identify key WUE regulators. GWAS across water conditions in 228 accessions identifies 73 quantitative trait loci (QTLs) for WUE-traits. Transcriptome profiling of 110 diverse accessions reveals 28 drought-responsive modules. eQTL mapping uncovers 146,966 regulatory variants, including condition-specific hotspots associated with key drought-related pathways. Integrative analysis underscores 85 high-confidence candidate genes, notably TaMYB7-A1. Overexpression of TaMYB7-A1 enhances photosynthesis, WUE, root development, and grain yield under drought condition by activating TaPIP2;2-B1 (water transport), TaRD20-D1 (stomatal regulation), and TaABCB4-B1 (root growth), reflecting reduced water loss and improved physiological resilience. Our study presents a comprehensive regulatory map and robust targets for wheat drought adaptation and resilient cultivar breeding.

Triticum↗
PubMed2026

Proinsulin regulators identified with CRISPR screen and in vivo mouse QTL mapping.

Altered proinsulin levels in β-cells and bloodstream are hallmarks of diabetes and other diseases, but our knowledge about the proinsulin regulators remains limited. Here we perform a genome-wide CRISPR screen to identify 84 proinsulin regulators that alter intracellular proinsulin/insulin ratio in a mouse β-cell line. The proinsulin regulators are distinct from the insulin regulators from a previous orthogonal CRISPR screen. Functional annotation of the proinsulin regulators highlights Golgi as the primary organelle for proinsulin storage and regulation. Trafficking towards the Golgi increases the intra-cellular proinsulin/insulin ratio, while trafficking away from the Golgi, including exocytosis and Golgi-to-ER retrograde transport, decreases the intracellular proinsulin levels. We also map mouse quantitative trait loci (QTLs) associated with plasma proinsulin levels and use the CRISPR screen results to pinpoint the causal genes within the QTL loci. Interestingly, protein disulfide isomerase Pdia6 is the strongest hit from both CRISPR screen and the in vivo QTL mapping. Knocking down Pdia6 significantly reduce proinsulin accumulation in Golgi and secretory granules. Intriguingly, Pdia6-depletion in both human and mouse β-cells does not affect the folding status of proinsulin but causes significantly impaired proinsulin production through a UPR-independent mechanism. Taken together, our genetic profiles provide mechanistic insights into the regulation of proinsulin/insulin homeostasis.

Animals↗
PubMed2026

Analysis of deep-resequencing data of 984 soybean accessions reveals structural variations underlying agronomic traits.

Genomic structural variants (SVs) are major sources of genetic variation and have profound impacts on phenotypic traits. However, their functional effects remain largely unexplored in soybean. Here, we resequence 940 soybean accessions. Together with 44 publicly available datasets, we identify 602,281 SVs. Using a graph-based genome, we detect an additional 58,760 presence/absence variations (PAVs) that broadly affect gene expression. Population genomic analyses reveal that SVs serve as a core driving force for soybean domestication and improvement. Integrating SVs with QTLs for oil and protein content, and performing GWAS on 27 traits, we identify key functional SVs. These include transposable element insertions altering seed coat color, multiple insertions within a cytochrome P450 gene modifying flower and hypocotyl color, and a GmMATE1 deletion enhancing seed size. Together, our study establishes a comprehensive SV map of soybean, offering a valuable resource for dissecting the genetic basis of complex traits to accelerate molecular breeding.

Glycine max↗
PubMed2026

EST-SSR based genetic polymorphism among Lablab (Lablab purpureus L. Sweet) accessions contrasting for drought stress at seedling stage.

Lablab is a multipurpose and the most drought-tolerant (DT) crop compared with its relatives. Despite its potential, Lablab is still an underutilized crop with a lack of improved varieties in many countries. The DT (D349, D147, HA4, D363, D352, D359, D348, D311, D55 and D250) and drought-susceptible (DS) (D271, D66, D106, D6, D26, D255, D28, D186, D95, and D258) accessions were earlier identified according to their morphological and biochemical responses to moisture stress at the seedling stage. These accessions were used to establish genetic polymorphism among the accessions contrasting for drought stress based on the Expressed Sequence Tag-Simple Sequence Repeats (EST-SSR) markers. The CTAB protocol was employed for the genomic DNA extraction. After DNA quality and quantity verification, the PCR was conducted using 16 EST-SSR primer pairs specific to the Lablab. The products were separated through the horizontal polyacrylamide gel electrophoresis (hPAGE). Discriminating ability of the markers and primers' efficiency were evaluated based on various genetic parameters. Principal Coordinate Analysis (PCoA) was performed to estimate the distance matrix among the population and among the accessions. While cluster analysis was processed to trace the genetic relationship among the accessions, dendrogram was constructed to decipher their genetic relationship. Analysis of Molecular Variance (AMOVA) was finally computed to quantify the diversity level and genetic relationship among the population, and among the accessions. A low polymorphism (GD = 0.19) was observed between the DT and DS accessions, likely due to limited discriminatory power of the EST-SSR markers. However, the PCoA, cluster analysis and AMOVA identified DT (D147, HA4, and D349) and DS (D106, D95, and D271) accessions as strongly contrasting populations under drought stress, with D147, HA4, D349, D363, D359, D352, and D348 further recommended as DT accessions. Given the low polymorphism observed, further validation using more informative molecular markers and advanced genomic approaches is recommended to improve the identification of drought-tolerance genes and related QTLs to support Lablab breeding programs.

Expressed Sequence Tags↗
PubMed2026

Genome-wide association identifies and validates genomic region controlling grain yield and agronomic traits in extra-early orange maize inbred lines under drought.

In order to meet the expected maize yield by 2050, breeders must work to improve breeding program efficiency by intensifying the implementation of new and improved technologies such as marker-assisted selection (MAS). Dissecting the genomic regions associated with drought tolerance is the first step forward in MAS program deployment for maize improvement under drought stress. Genome-wide association studies (GWAS) were used to investigate and identify quantitative trait loci (QTLs) associated with six traits under drought stress. One hundred and eighty-seven extra-early orange maize inbred lines were evaluated under managed drought stress at Ikenne, in Nigeria, during the 2022 and 2023 dry seasons. The materials were also genotyped using 9355 DArTseq SNP markers and analyzed using the enriched compressed mixed linear model (ECMLM). Enriched compressed mixed linear model was used for association-trait analysis. The ECMLM-based GWAS identified 45 candidate genomic loci associated with the six traits, including five for grain yield, with R2 ranging from 8.79 to 25.3%. Independent validation using the multi-locus 3VmrMLM approach confirmed seven high-confidence genomic loci consistently detected by both methods across grain yield, anthesis-silking interval, ear aspect, and ears per plant, providing additional statistical support for these genomic regions. Candidate gene annotation identified biologically relevant genes underlying the validated loci, including Zm00001eb238250 (protein-serine/threonine phosphatase), Zm00001eb040940 (trehalose-phosphatase), Zm00001eb117820 (homeobox protein knotted-1-like 4), Zm00001eb145560 (zinc ion-binding protein), and Zm00001eb294180 (WRKY DNA-binding domain protein), suggesting their potential roles in drought adaptation and grain productivity. These findings improve our understanding of the genetic architecture of drought tolerance in extra-early orange maize and provide valuable genomic resources for accelerating drought-resilient maize breeding.

Zea mays↗
PubMed2026

Haplotype-specific expression of a terpene synthase underlies linalool variation in the grapevine cultivar Riesling.

Grapevine cultivars vary widely in monoterpenoid content, yet the genetic and regulatory mechanisms underlying this variation remain poorly characterized beyond highly aromatic Muscat types. We profiled free volatiles and monoterpenoid glycosides in a Riesling × Cabernet Sauvignon F1 mapping population, revealing extensive variation and transgressive segregation consistent with multigenic control. QTL mapping identified 70 significant loci associated with 48 volatile compounds and monoterpene glycosides, including two major QTLs explaining 33.6% and 33.4% of phenotypic variance in (3S)-linalool accumulation. Integration of haplotype-resolved transcriptomics with metabolite data, enabled by a chromosome-scale diploid Riesling genome assembly, resolved a (3S)-linalool/nerolidol synthase cluster on chromosome 10 and identified VviTPS54 as the strongest candidate underlying linalool variation. VviTPS54 exhibited haplotype-specific expression strongly correlated with (3S)-linalool accumulation across genotypes, while no QTL was detected at the 1-deoxy-D-xylulose-5-phosphate synthase 1 (VviDXS1) locus previously identified in Muscat cultivars. In addition, VviDXS1 expression was not correlated with terpene levels, indicating that regulatory variation within terpene synthase clusters, rather than methylerythritol phosphate (MEP) pathway flux, drives monoterpenoid composition in this population. These results establish regulatory variation of terpene synthases as a key mechanism underlying monoterpenoid diversity in grapevine and demonstrate that resolving such variation requires haplotype-phased genome assemblies coupled with haplotype-resolved transcriptomics to detect allele-specific expression differences at complex, heterozygous loci.

Grapevine↗
PubMed2026

Integrative multi-omics quantitative trait loci prioritize CASP7 as a candidate protective gene for cataract.

Cataracts are the leading cause of vision loss worldwide. Despite surgery being the only effective treatment, its economic burden highlights the necessity of exploring the pathogenesis of cataracts. In this study, we analyzed 4 large-scale GWAS (genome-wide association study) datasets for cataracts and performed SMR analysis along with heterogeneity in dependent instruments (HEIDI) testing to explore the effects of methylation, expression, and protein QTLs on cataracts. We further validated shared genetic variants through COLOC analysis. Additionally, we searched datasets related to cataracts from the Gene Expression Omnibus (GEO) database for differentially expressed genes (DEGs) and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) enrichment analyses. By integrating summary-based Mendelian randomization (SMR) results with bioinformatics findings, CASP7 showed a consistent protective-direction association with cataract risk (mQTL: OR [95% CI] = 0.959 [0.941-0.977], FDR-adjusted P = .039; eQTL: OR [95% CI] = 0.897 [0.860-0.937], FDR-adjusted P = .0046; pQTL: OR [95% CI] = 0.597 [0.483-0.738], FDR-adjusted P = .00083). GEO-based analyses provided transcriptomic support for CASP7 involvement in cataract-related lens biology. These findings prioritize CASP7 as a genetically supported candidate protective gene associated with cataract risk. Because this study is based on public summary-level and transcriptomic datasets, the results should be interpreted cautiously and require functional validation in human lens-relevant systems.

Quantitative Trait Loci↗
PubMed2025

Dynamic neuro-immune regulation of psychiatric risk loci in human neurons.

The prenatal environment influences neurodevelopment and subsequent clinical trajectories for psychiatric outcomes in childhood and adolescence. Yet it remains unclear if the impact of maternal and fetal immune activation varies with distinct polygenic risk profiles. Therefore, here we catalogue genotype and environment (GxE) interactions, contrasting allele-specific regulatory activity between inflammatory contexts. We report a cue-specific neuronal massively parallel reporter assay (MPRA) of 220 loci from genome-wide association study (GWAS) linked to ten brain traits/disorders, empirically dissecting the impact of interleukin-6 (IL-6) and interferon-alpha (IFNα) on transcriptional activity. Of 1,469 active candidate regulatory risk elements (MPRA-active CRSs) across three conditions, we identify 316 with dynamic variant-specific effects (MPRA-QTLs) in human induced pluripotent stem cell (hiPSC)-derived glutamatergic neurons. Broadly, across hundreds of variants, neuronal immune-mediated regulatory activity is driven by differences in transcription factor binding and chromatin accessibility, the gene targets of which show pleiotropic enrichments for brain, metabolic, and immune disorders. Dynamic genetic regulation mediates immune effects, informing our understanding of mechanisms governing pleiotropy and variable penetrance. Understanding neurodevelopmental GxE interactions will inform mental health trajectories and resolve mechanisms mediating prenatal risk.

dynamic expression quantitative trait loci↗
PubMed2019

Unique genetic basis of the distinct antibiotic potency of high acetic acid production in the probiotic yeast Saccharomyces cerevisiae var. boulardii.

The yeast Saccharomyces boulardii has been used worldwide as a popular, commercial probiotic, but the basis of its probiotic action remains obscure. It is considered conspecific with budding yeast Saccharomyces cerevisiae, which is generally used in classical food applications. They have an almost identical genome sequence, making the genetic basis of probiotic potency in S. boulardii puzzling. We now show that S. boulardii produces at 37°C unusually high levels of acetic acid, which is strongly inhibitory to bacterial growth in agar-well diffusion assays and could be vital for its unique application as a probiotic among yeasts. Using pooled-segregant whole-genome sequence analysis with S. boulardii and S. cerevisiae parent strains, we succeeded in mapping the underlying QTLs and identified mutant alleles of SDH1 and WHI2 as the causative alleles. Both genes contain a SNP unique to S. boulardii (sdh1 F317Y and whi2 S287*) and are fully responsible for its high acetic acid production. S. boulardii strains show different levels of acetic acid production, depending on the copy number of the whi2 S287* allele. Our results offer the first molecular explanation as to why S. boulardii could exert probiotic action as opposed to S. cerevisiae They reveal for the first time the molecular-genetic basis of a probiotic action-related trait in S. boulardii and show that antibacterial potency of a probiotic microorganism can be due to strain-specific mutations within the same species. We suggest that acquisition of antibacterial activity through medium acidification offered a selective advantage to S. boulardii in its ecological niche and for its application as a probiotic.

Acetic Acid↗
PubMed2026

Genome-wide association studies of plant traits and functional analysis of leaf development-related genes in citrus.

Labor-saving and high-light-efficiency tree architecture is a key breeding objective for woody fruit trees like citrus. However, population genetics information on these traits remains limited. In this study, tree architecture, thorn, and leaf traits were evaluated in 353 F2 progeny derived from a cross between Clementine mandarin and precocious trifoliate orange-an early-flowering variety. A random subset of 300 offspring was sequenced for a genome-wide association study (GWAS), which detected 10 216 significantly associated SNPs and defined several major quantitative trait loci (QTLs) for the target traits. Subsequent bulked segregant analysis (BSA) and GWAS on individuals with extreme compound leaf phenotypes mapped the causal gene(s) to a 0.8 Mb region (22.15-22.95 Mb) on chromosome 4. Genetic analysis across multiple hybrid combinations confirmed that the compound leaf trait in trifoliate orange is dominantly inherited and follows Mendelian segregation. Transcriptome profiling of parental leaves at different developmental stages identified a KNOX gene, CiKNAT6, as a candidate. Further validation using CAPS markers and Hi-Tom sequencing demonstrated tight linkage between an InDel polymorphism in CiKNAT6 and leaf shape across diverse citrus species and the F2 population, with co-segregation observed for the compound leaf trait. Due to alternative splicing producing seven splice variants, the CiKNAT6 DNA sequence was selected for genetic transformation experiments. Functional analysis revealed that the Clementine mandarin allele of CiKNAT6 is non-functional owing to an InDel, whereas ectopic expression of the trifoliate orange allele in tobacco and lemon induced leaf curling and reduced leaf size. CRISPR-Cas9 knockout of CiKNAT6 in trifoliate orange resulted in increased leaf area. These findings provide valuable genetic resources and insights for future studies on tree architecture and leaf morphology.

Plant Leaves↗