Power, false discovery rate and Winner's Curse in eQTL studies

Disentangling the divergent causal pathways underlying the association between body mass index and bone mineral density: a comprehensive Mendelian randomization study

BACKGROUND

While the protective role of body mass index (BMI) in bone mass has been well-documented, the divergent associations between BMI and estimated bone mineral density (eBMD), attributed to its highly heterogeneous nature, remain insufficiently understood.

METHODS

Leveraging the hitherto largest genome-wide summary statistics, we conducted a two-sample Mendelian randomization (MR) to re-evaluate the effect of genetically predicted BMI on eBMD. Then, MR-Clust was applied to examine the potential presence of distinct causal pathways underlying the BMI-eBMD link. Utilizing tissue-partitioned MR, we estimated the distinct effects of separated tissue-specific subcomponents of BMI on eBMD, further supplemented by multivariable MR of body composition phenotypes on eBMD.

RESULTS

We reconfirmed the significant positive association between genetically predicted BMI and eBMD (βIVW = 0.13, P value = 1.28 × 10-34). Potential distinct causal pathways contributing to the observed total effect were identified by MR-Clust, with some exerting a protective effect while others leading to its deterioration. Tissue-partitioned MR suggested a marginally independent protective association between skeletal muscle-tissue instrumented BMI and eBMD (βIVW = 0.14, P value = 4.98 × 10-2) after accounting for adipose-tissue instrumented BMI, which was supported by the independent association between genetically predicted lean mass and eBMD after accounting for other body composition phenotypes.

CONCLUSIONS

Our results shed preliminary insights into the intricate relationship between obesity and bone mass, highlighting divergent causal pathways underlying the association between BMI and eBMD. Our findings emphasize the potential importance of precision obesity management over merely a general indicator as BMI in future public health strategies for osteoporosis prevention.

TEX10: A Novel Drug Target and Potential Therapeutic Direction for Sleep Apnea Syndrome

Background

Sleep apnea syndrome (SAS) is a prevalent sleep disorder strongly associated with obesity, metabolic dysregulation, and cardiovascular diseases. While its underlying pathophysiological mechanisms remain incompletely understood, genetic factors likely play a pivotal role in SAS pathogenesis. This study investigates the causal relationships between potential drug target genes and SAS using multiple statistical approaches, aiming to provide novel insights for targeted therapeutic development.

Methods

We conducted a comprehensive genetic analysis integrating multiple methodologies to investigate gene-SAS relationships. Using publicly available GWAS and eQTL databases, we performed Mendelian Randomization (MR) analysis with the inverse variance weighted (IVW) method, validated by weighted median and MR-Egger approaches. Summary-data-based MR (SMR) analysis, coupled with HEIDI testing, assessed direct gene expression-SAS associations while controlling for linkage disequilibrium (LD). Colocalization analysis evaluated the probability of shared causal variants between SNPs, gene expression, and SAS. Statistical significance was determined using Benjamini-Hochberg multiple testing correction (FDR < 0.05). Additionally, mediation analysis explored TEX10's influence on SAS through metabolic intermediates including BMI, waist circumference, and HDL cholesterol.

Results

We identified 18 candidate drug target genes significantly associated with SAS, with MAPKAPK3, TNXB, MPHOSPH8, and TEX10 showing consistent associations across multiple analyses. TEX10, in particular, exhibited significant associations with SAS risk in blood, cerebral cortex, hippocampus, and basal ganglia (PP.H4 > 0.9). Mediation analysis suggested that TEX10 might influence SAS risk indirectly through BMI, waist circumference, and HDL cholesterol levels.

Conclusion

Our study identified multiple potential therapeutic targets causally linked to SAS, with TEX10 emerging as a key candidate gene. These findings advance our understanding of SAS pathogenesis and offer promising directions for personalized diagnostics and targeted therapies.

Multiomic QTL mapping reveals phenotypic complexity of GWAS loci and prioritizes putative causal variants

Most GWAS loci are presumed to affect gene regulation; however, only ∼43% colocalize with expression quantitative trait loci (eQTLs). To address this colocalization gap, we map eQTLs, chromatin accessibility QTLs (caQTLs), and histone acetylation QTLs (haQTLs) using molecular samples from three early developmental-like tissues. Through colocalization, we annotate 10.4% (n = 540) of GWAS loci in 15 traits by QTL phenotype, temporal specificity, and complexity. We show that integration of chromatin QTLs results in a 2.3-fold higher annotation rate of GWAS loci because they capture distal GWAS loci missed by eQTLs, and that 5.4% (n = 13) of GWAS colocalizing eQTLs are early developmental specific. Finally, we utilize the iPSCORE multiomic QTLs to prioritize putative causal variants overlapping transcription factor motifs to elucidate the potential genetic underpinnings of 296 GWAS-QTL colocalizations.

reg-eQTL: Integrating transcription factor effects to unveil regulatory variants

Regulatory single-nucleotide variants (rSNVs) in noncoding regions of the genome play a crucial role in gene transcription by altering transcription factor (TF) binding, chromatin states, and other epigenetic modifications. Existing expression quantitative trait locus (eQTL) methods identify genomic loci associated with gene-expression changes, but they often fall short in pinpointing causal variants. We introduce reg-eQTL, a computational method that incorporates TF effects and interactions with genetic variants into eQTL analysis. This approach provides deeper insights into the regulatory mechanisms, bringing us one step closer to identifying potential causal variants by uncovering how TFs interact with SNVs to influence gene expression. This method defines a trio consisting of a genetic variant, a target gene, and a TF and tests its impact on gene transcription. In comprehensive simulations, reg-eQTL shows improved power of detecting rSNVs with low population frequency, weak effects, and synergetic interaction with TF as compared to traditional eQTL methods. Application of reg-eQTL to GTEx data from lung, brain, and whole-blood tissues uncovered regulatory trios that include eQTLs and increased the number of eQTLs shared across tissue types. Regulatory networks constructed on the basis of these trios reveal intricate gene regulation across tissue types.

Expression quantitative trait loci associated with performance traits, blood biochemical parameters, and cytokine profile in pigs

Identifying expression Quantitative Trait Loci (eQTL) and functional candidate variants associated with blood biochemical parameters can contribute to the understanding of genetic mechanisms underlying phenotypic variation in complex traits in pigs. We identified eQTLs through gene expression levels in muscle and liver tissues of Large White pigs. The identified eQTL were then tested for association with biochemical parameters, cytokine profiles, and performance traits of pigs. A total of 41,759 SNPs and 15,093 and 15,516 expression gene levels from muscle and liver tissues, respectively, enabled the identification of 1,199 eQTL. The eQTL identified related the SNP rs345667860 as significantly associated with interleukin-6 and interleukin-18 in liver tissue, while the rs695637860 SNP was associated with aspartate aminotransferase and interleukin-6, and rs337362164 was associated with high-density lipoprotein of the blood serum. In conclusion, the identification of three eQTL significantly associated with aspartate aminotransferase and cytokine levels in both serum and liver tissues suggests a potential role for these variants in modulating immune function and overall health in production pigs. Further research is needed to validate these findings and explore their potential for improving pig health and productivity.

The contribution of genetic determinants of blood gene expression and splicing to molecular phenotypes and health outcomes

The biological mechanisms through which most nonprotein-coding genetic variants affect disease risk are unknown. To investigate gene-regulatory mechanisms, we mapped blood gene expression and splicing quantitative trait loci (QTLs) through bulk RNA sequencing in 4,732 participants and integrated protein, metabolite and lipid data from the same individuals. We identified cis-QTLs for the expression of 17,233 genes and 29,514 splicing events (in 6,853 genes). Colocalization analyses revealed 3,430 proteomic and metabolomic traits with a shared association signal with either gene expression or splicing. We quantified the relative contribution of the genetic effects at loci with shared etiology, observing 222 molecular phenotypes significantly mediated by gene expression or splicing. We uncovered gene-regulatory mechanisms at disease loci with therapeutic implications, such as WARS1 in hypertension, IL7R in dermatitis and IFNAR2 in COVID-19. Our study provides an open-access resource on the shared genetic etiology across transcriptional phenotypes, molecular traits and health outcomes in humans ( https://IntervalRNA.org.uk ).

Secure and federated quantitative trait loci mapping with privateQTL

Understanding the relationship between genotypes and phenotypes is crucial for advancing personalized medicine. Expression quantitative trait loci (eQTL) mapping plays a significant role by correlating genetic variants to gene expression levels. Despite the progress made by large-scale projects, eQTL mapping still faces challenges in statistical power and privacy concerns. Multi-site studies can increase sample sizes but are hindered by privacy issues. We present privateQTL, a novel framework leveraging secure multi-party computation for secure and federated eQTL mapping. When tested in a real-world scenario with data from different studies, privateQTL outperformed meta-analysis by accurately correcting for covariates and batch effect and retaining higher accuracy and precision for both eGene-eVariant mapping and effect size estimation. In addition, privateQTL is modular and scalable, making it adaptable for other molecular phenotypes and large-scale studies. Our results indicate that privateQTL is a practical solution for privacy-preserving collaborative eQTL mapping.

Stimulating Wnt signaling reveals context-dependent genetic effects on gene regulation in primary human neural progenitors

Gene regulatory effects have been difficult to detect at many non-coding loci associated with brain-related traits, likely because some genetic variants have distinct functions in specific contexts. To explore context-dependent gene regulation, we measured chromatin accessibility and gene expression after activation of the canonical Wnt pathway in primary human neural progenitors (n = 82 donors). We found that TCF/LEF motifs and brain-structure-associated and neuropsychiatric-disorder-associated variants were enriched within Wnt-responsive regulatory elements. Genetically influenced regulatory elements were enriched in genomic regions under positive selection along the human lineage. Wnt pathway stimulation increased detection of genetically influenced regulatory elements/genes by 66%/53% and enabled identification of 397 regulatory elements primed to regulate gene expression. Stimulation also increased identification of shared genetic effects on molecular and complex brain traits by up to 70%, suggesting that genetic variant function during neurodevelopmental patterning can lead to differences in adult brain and behavioral traits.

BMI Interacts with the Genome to Regulate Gene Expression Globally, with Emphasis in the Brain and Gut

Genome-wide association studies identify common genomic variants associated with disease across a population. Individual environmental effects are often not included, despite evidence that environment mediates genomic regulation of higher order biology. Body mass index (BMI) is associated with complex disorders across clinical specialties, yet has not been modeled as a genomic environment. Here, we tested for expression quantitative trait (eQTL) loci that contextually regulate gene expression across the BMI spectrum using an interaction approach. We parsed the impact of cell type, enhancer interactions, and created novel BMI-dynamic gene expression predictor models. We found that BMI main effects associated with endocrine gene expression, while interactive variant-by-BMI effects impacted gene expression in the brain and gut. Cortical BMI-dynamic loci were experimentally dysregulated by inflammatory cytokines in an in vitro system. Using BMI-dynamic models, we identify novel genes in nitric oxide signaling pathways in the nucleus accumbens significantly associated with depression and smoking. While neither genetics nor BMI are sufficient as standalone measures to capture the complexity of downstream cellular consequences, including environment powers disease gene discovery.

Private information leakage from single-cell count matrices

The increase in publicly available human single-cell datasets, encompassing millions of cells from many donors, has significantly enhanced our understanding of complex biological processes. However, the accessibility of these datasets raises significant privacy concerns. Due to the inherent noise in single-cell measurements and the scarcity of population-scale single-cell datasets, recent private information quantification studies have focused on bulk gene expression data sharing. To address this gap, we demonstrate that individuals in single-cell gene expression datasets are vulnerable to linking attacks, where attackers can infer their sensitive phenotypic information using publicly available tissue or cell-type-specific expression quantitative trait loci (eQTLs) information. We further develop a method for genotype prediction and genotype-phenotype linking that remains effective without relying on eQTL information. We show that variants from one study can be exploited to uncover private information about individuals in another study.

A brief guide to analyzing expression quantitative trait loci

Molecular quantitative trait locus (molQTL) mapping has emerged as an important approach for elucidating the functional consequences of genetic variants and unraveling the causal mechanisms underlying diseases or complex traits. However, the variety of analysis tools and sophisticated methodologies available for molQTL studies can be overwhelming for researchers with limited computational expertise. Here, we provide a brief guideline with a curated list of methods and software tools for analyzing expression quantitative trait loci, the most widely studied type of molQTL.

Transcriptome-wide Mendelian randomization during CD4+ T cell activation reveals immune-related drug targets for cardiometabolic diseases

Immunity has shown potentials in informing drug development for cardiometabolic diseases, such as type 2 diabetes (T2D) and coronary artery disease (CAD). Here, we performed a transcriptome-wide Mendelian randomization (MR) study to estimate the putative causal effects of 11,021 gene expression profiles during CD4+ T cells activation on the development of T2D and CAD. Robust MR and colocalization evidence was observed for 162 genes altering T2D risk and 80 genes altering CAD risk, with 12% and 16% respectively demonstrating CD4+ T cell specificity. We observed temporal causal patterns during T cell activation in 69 gene-T2D pairs and 34 gene-CAD pairs. These genes were eight times more likely to show robust genetic evidence. We further identified 25 genes that were targets for drugs under clinical investigation, including LIPA and GCK. This study provides evidence to support immune-to-metabolic disease connections, and prioritises immune-mediated drug targets for cardiometabolic diseases.

Quantitative trait locus mapping in placenta: A comparative study of chorionic villus and birth placenta

The placenta, a pivotal player in the prenatal environment, holds crucial insights into early developmental pathways and future health outcomes. In this study, we explored genetic molecular regulation in chorionic villus samples (CVS) from the first trimester and placenta tissue at birth. We assessed quantitative trait locus (QTL) mapping on DNA methylation and gene expression data in a Finnish cohort of 574 individuals. We found more QTLs in birth placenta than in first-trimester placenta. Nevertheless, a substantial amount of associations overlapped in their effects and showed consistent direction in both tissues, with increasing molecular genetic effects from early pregnancy to birth placenta. The identified QTLs in birth placenta were most enriched in genes with placenta-specific expression. Conducting a phenome-wide-association study (PheWAS) on the associated SNPs, we observed numerous overlaps with genome-wide association study (GWAS) hits (spanning 57 distinct traits and 23 SNPs), with notable enrichments for immunological, skeletal, and respiratory traits. The QTL-SNP rs1737028 (chr6:29737993) presented with the highest number of GWAS hits. This SNP was related to HLA-G expression via DNA methylation and was associated with various immune, respiratory, and psychiatric traits. Our findings implicate increasing genetic molecular regulation during the course of pregnancy and support the involvement of placenta gene regulation, particularly in immunological traits. This study presents a framework for understanding placenta-specific gene regulation during pregnancy and its connection to health-related traits.

Smoking-informed methylation and expression QTLs in human brain and colocalization with smoking-associated genetic loci

Smoking is a leading cause of preventable morbidity and mortality. Smoking is heritable, and genome-wide association studies (GWASs) of smoking behaviors have identified hundreds of significant loci. Most GWAS-identified variants are noncoding with unknown neurobiological effects. We used genome-wide genotype, DNA methylation, and RNA sequencing data in postmortem human nucleus accumbens (NAc) to identify cis-methylation/expression quantitative trait loci (meQTLs/eQTLs), investigate variant-by-cigarette smoking interactions across the genome, and overlay QTL evidence at smoking GWAS-identified loci to evaluate their regulatory potential. Active smokers (N = 52) and nonsmokers (N = 171) were defined based on cotinine biomarker levels and next-of-kin reporting. We simultaneously tested variant and variant-by-smoking interaction effects on methylation and expression, separately, adjusting for biological and technical covariates and correcting for multiple testing using a two-stage procedure. We found >2 million significant meQTL variants (padj < 0.05) corresponding to 41,695 unique CpGs. Results were largely driven by main effects, and five meQTLs, mapping to NUDT12, FAM53B, RNF39, and ADRA1B, showed a significant interaction with smoking. We found 57,683 significant eQTL variants for 958 unique eGenes (padj < 0.05) and no smoking interactions. Colocalization analyses identified loci with smoking-associated GWAS variants that overlapped meQTLs/eQTLs, suggesting that these heritable factors may influence smoking behaviors through functional effects on methylation/expression. One locus containing MUSTN1 and ITIH4 colocalized across all data types (GWAS, meQTL, and eQTL). In this first genome-wide meQTL map in the human NAc, the enriched overlap with smoking GWAS-identified genetic loci provides evidence that gene regulation in the brain helps explain the neurobiology of smoking behaviors.

Genetics of cell-type-specific post-transcriptional gene regulation during human neurogenesis

The function of some genetic variants associated with brain-relevant traits has been explained through colocalization with expression quantitative trait loci (eQTL) conducted in bulk postmortem adult brain tissue. However, many brain-trait associated loci have unknown cellular or molecular function. These genetic variants may exert context-specific function on different molecular phenotypes including post-transcriptional changes. Here, we identified genetic regulation of RNA editing and alternative polyadenylation (APA) within a cell-type-specific population of human neural progenitors and neurons. More RNA editing and isoforms utilizing longer polyadenylation sequences were observed in neurons, likely due to higher expression of genes encoding the proteins mediating these post-transcriptional events. We also detected hundreds of cell-type-specific editing quantitative trait loci (edQTLs) and alternative polyadenylation QTLs (apaQTLs). We found colocalizations of a neuron edQTL in CCDC88A with educational attainment and a progenitor apaQTL in EP300 with schizophrenia, suggesting that genetically mediated post-transcriptional regulation during brain development leads to differences in brain function.

A Deep Dive into Statistical Modeling of RNA Splicing QTLs Reveals New Variants that Explain Neurodegenerative Disease

Genome-wide association studies (GWAS) have identified thousands of putative disease causing variants with unknown regulatory effects. Efforts to connect these variants with splicing quantitative trait loci (sQTLs) have provided functional insights, yet sQTLs reported by existing methods cannot explain many GWAS signals. We show current sQTL modeling approaches can be improved by considering alternative splicing representation, model calibration, and covariate integration. We then introduce MAJIQTL, a new pipeline for sQTL discovery. MAJIQTL includes two new statistical methods: a weighted multiple testing approach for sGene discovery and a model for sQTL effect size inference to improve variant prioritization. By applying MAJIQTL to GTEx, we find significantly more sGenes harboring sQTLs with functional significance. Notably, our analysis implicates the novel variant rs582283 in Alzheimer's disease. Using antisense oligonucleotides, we validate this variant's effect by blocking the implicated YBX3 binding site, leading to exon skipping in the gene MS4A3.

Mapping Epigenetic Gene Variant Dynamics: Comparative Analysis of Frequency, Functional Impact and Trait Associations in African and European Populations

Epigenetic modifications influence gene expression levels, impact organismal traits, and play a role in the development of diseases. Therefore, variants in genes involved in epigenetic processes are likely to be important in disease susceptibility, and the frequency of variants may vary between populations with African and European ancestries. Here, we analyse an integrated dataset to define the frequencies, associated traits, and functional impact of epigenetic gene variants among individuals of African and European ancestry represented in the UK Biobank. We find that the frequencies of 88.4% of epigenetic gene variants significantly differ between these groups. Furthermore, we find that the variants are associated with many traits and diseases, and some of these associations may be population-specific owing to allele frequency differences. Additionally, we observe that variants associated with traits are significantly enriched for quantitative trait loci that affect DNA methylation, chromatin accessibility, and gene expression. We find that methylation quantitative trait loci account for 71.2% of the variants influencing gene expression. Moreover, variants linked to biomarker traits exhibit high correlation. We therefore conclude that epigenetic gene variants associated with traits tend to differ in their allele frequencies among African and European populations and are enriched for QTLs.

Sources of gene expression variation in a globally diverse human cohort

Genetic variation that influences gene expression and splicing is a key source of phenotypic diversity1-5. Although invaluable, studies investigating these links in humans have been strongly biased towards participants of European ancestries, which constrains generalizability and hinders evolutionary research. Here to address these limitations, we developed MAGE, an open-access RNA sequencing dataset of lymphoblastoid cell lines from 731 individuals from the 1000 Genomes Project6, spread across 5 continental groups and 26 populations. Most variation in gene expression (92%) and splicing (95%) was distributed within versus between populations, which mirrored the variation in DNA sequence. We mapped associations between genetic variants and expression and splicing of nearby genes (cis-expression quantitative trait loci (eQTLs) and cis-splicing QTLs (sQTLs), respectively). We identified more than 15,000 putatively causal eQTLs and more than 16,000 putatively causal sQTLs that are enriched for relevant epigenomic signatures. These include 1,310 eQTLs and 1,657 sQTLs that are largely private to underrepresented populations. Our data further indicate that the magnitude and direction of causal eQTL effects are highly consistent across populations. Moreover, the apparent 'population-specific' effects observed in previous studies were largely driven by low resolution or additional independent eQTLs of the same genes that were not detected. Together, our study expands our understanding of human gene expression diversity and provides an inclusive resource for studying the evolution and function of human genomes.

eQTLs identify regulatory networks and drivers of variation in the individual response to sepsis

Sepsis is a clinical syndrome of life-threatening organ dysfunction caused by a dysregulated response to infection, for which disease heterogeneity is a major obstacle to developing targeted treatments. We have previously identified gene-expression-based patient subgroups (sepsis response signatures [SRS]) informative for outcome and underlying pathophysiology. Here, we aimed to investigate the role of genetic variation in determining the host transcriptomic response and to delineate regulatory networks underlying SRS. Using genotyping and RNA-sequencing data on 638 adult sepsis patients, we report 16,049 independent expression (eQTLs) and 32 co-expression module (modQTLs) quantitative trait loci in this disease context. We identified significant interactions between SRS and genotype for 1,578 SNP-gene pairs and combined transcription factor (TF) binding site information (SNP2TFBS) and predicted regulon activity (DoRothEA) to identify candidate upstream regulators. Overall, these approaches identified putative mechanistic links between host genetic variation, cell subtypes, and the individual transcriptomic response to infection.

Genome-wide association and expression quantitative trait loci in cattle reveals common genes regulating mammalian fertility

Most genetic variants associated with fertility in mammals fall in non-coding regions of the genome and it is unclear how these variants affect fertility. Here we use genome-wide association summary statistics for Heifer puberty (pubertal or not at 600 days) from 27,707 Bos indicus, Bos taurus and crossbred cattle; multi-trait GWAS signals from 2119 indicine cattle for four fertility traits, including days to calving, age at first calving, pregnancy status, and foetus age in weeks (assessed by rectal palpation of the foetus); and expression quantitative trait locus for whole blood from 489 indicine cattle, to identify 87 putatively functional genes affecting cattle fertility. Our analysis reveals a significant overlap between the set of cattle and previously reported human fertility-related genes, impling the existence of a shared pool of genes that regulate fertility in mammals. These findings are crucial for developing approaches to improve fertility in cattle and potentially other mammals.

Genetic influence on vascular smooth muscle cell apoptosis

Vascular smooth muscle cell (VSMC) proliferation, migration, and apoptosis play important roles in many physiological processes and pathological conditions. To identify genetic influences on VSMC behavior, we measured these traits and undertook genome-wide association studies in primary umbilical artery-derived VSMCs from >2000 individuals. Although there were no genome-wide significant associations for VSMC proliferation or migration, genetic variants at two genomic loci (7p15.3 and 7q32.3) showed highly significant associations with VSMC apoptosis (P = 1.95 × 10-13 and P = 7.47 × 10-9, respectively). The lead variant at the 7p51.3 locus was associated with increased expression of the GSDME and PALS2 genes in VSMCs. Knockdown of GSDME or PALS2 in VSMCs attenuated apoptotic cell death. A protein co-immunoprecipitation assay indicated that GSDME complexed with PALS2. PALS2 knockdown attenuated activated caspase-3 and GSDME fragmentation, whilst GSDME knockdown also reduced activated caspase-3. These findings provide new insights into the genetic regulation of VSMC apoptosis, with potential utility for therapeutic development.

Control of false discoveries in grouped hypothesis testing for eQTL data

BACKGROUND

Expression quantitative trait locus (eQTL) analysis aims to detect the genetic variants that influence the expression of one or more genes. Gene-level eQTL testing forms a natural grouped-hypothesis testing strategy with clear biological importance. Methods to control family-wise error rate or false discovery rate for group testing have been proposed earlier, but may not be powerful or easily apply to eQTL data, for which certain structured alternatives may be defensible and may enable the researcher to avoid overly conservative approaches.

RESULTS

In an empirical Bayesian setting, we propose a new method to control the false discovery rate (FDR) for grouped hypotheses. Here, each gene forms a group, with SNPs annotated to the gene corresponding to individual hypotheses. The heterogeneity of effect sizes in different groups is considered by the introduction of a random effects component. Our method, entitled Random Effects model and testing procedure for Group-level FDR control (REG-FDR), assumes a model for alternative hypotheses for the eQTL data and controls the FDR by adaptive thresholding. As a convenient alternate approach, we also propose Z-REG-FDR, an approximate version of REG-FDR, that uses only Z-statistics of association between genotype and expression for each gene-SNP pair. The performance of Z-REG-FDR is evaluated using both simulated and real data. Simulations demonstrate that Z-REG-FDR performs similarly to REG-FDR, but with much improved computational speed.

CONCLUSION

Our results demonstrate that the Z-REG-FDR method performs favorably compared to other methods in terms of statistical power and control of FDR. It can be of great practical use for grouped hypothesis testing for eQTL analysis or similar problems in statistical genomics due to its fast computation and ability to be fit using only summary data.

Multi-omic QTL mapping in early developmental tissues reveals phenotypic and temporal complexity of regulatory variants underlying GWAS loci

Most GWAS loci are presumed to affect gene regulation, however, only ∼43% colocalize with expression quantitative trait loci (eQTLs). To address this colocalization gap, we identify eQTLs, chromatin accessibility QTLs (caQTLs), and histone acetylation QTLs (haQTLs) using molecular samples from three early developmental (EDev) tissues. Through colocalization, we annotate 586 GWAS loci for 17 traits by QTL complexity, QTL phenotype, and QTL temporal specificity. We show that GWAS loci are highly enriched for colocalization with complex QTL modules that affect multiple elements (genes and/or peaks). We also demonstrate that caQTLs and haQTLs capture regulatory variations not associated with eQTLs and explain ∼49% of the functionally annotated GWAS loci. Additionally, we show that EDev-unique QTLs are strongly depleted for colocalizing with GWAS loci. By conducting one of the largest multi-omic QTL studies to date, we demonstrate that many GWAS loci exhibit phenotypic complexity and therefore, are missed by traditional eQTL analyses.

TBK1, a prioritized drug repurposing target for amyotrophic lateral sclerosis: evidence from druggable genome Mendelian randomization and pharmacological verification in vitro

BACKGROUND

There is a lack of effective therapeutic strategies for amyotrophic lateral sclerosis (ALS); therefore, drug repurposing might provide a rapid approach to meet the urgent need for treatment.

METHODS

To identify therapeutic targets associated with ALS, we conducted Mendelian randomization (MR) analysis and colocalization analysis using cis-eQTL of druggable gene and ALS GWAS data collections to determine annotated druggable gene targets that exhibited significant associations with ALS. By subsequent repurposing drug discovery coupled with inclusion criteria selection, we identified several drug candidates corresponding to their druggable gene targets that have been genetically validated. The pharmacological assays were then conducted to further assess the efficacy of genetics-supported repurposed drugs for potential ALS therapy in various cellular models.

RESULTS

Through MR analysis, we identified potential ALS druggable genes in the blood, including TBK1 [OR 1.30, 95%CI (1.19, 1.42)], TNFSF12 [OR 1.36, 95%CI (1.19, 1.56)], GPX3 [OR 1.28, 95%CI (1.15, 1.43)], TNFSF13 [OR 0.45, 95%CI (0.32, 0.64)], and CD68 [OR 0.38, 95%CI (0.24, 0.58)]. Additionally, we identified potential ALS druggable genes in the brain, including RESP18 [OR 1.11, 95%CI (1.07, 1.16)], GPX3 [OR 0.57, 95%CI (0.48, 0.68)], GDF9 [OR 0.77, 95%CI (0.67, 0.88)], and PTPRN [OR 0.17, 95%CI (0.08, 0.34)]. Among them, TBK1, TNFSF12, RESP18, and GPX3 were confirmed in further colocalization analysis. We identified five drugs with repurposing opportunities targeting TBK1, TNFSF12, and GPX3, namely fostamatinib (R788), amlexanox (AMX), BIIB-023, RG-7212, and glutathione as potential repurposing drugs. R788 and AMX were prioritized due to their genetic supports, safety profiles, and cost-effectiveness evaluation. Further pharmacological analysis revealed that R788 and AMX mitigated neuroinflammation in ALS cell models characterized by overly active cGAS/STING signaling that was induced by MSA-2 or ALS-related toxic proteins (TDP-43 and SOD1), through the inhibition of TBK1 phosphorylation.

CONCLUSIONS

Our MR analyses provided genetic evidence supporting TBK1, TNFSF12, RESP18, and GPX3 as druggable genes for ALS treatment. Among the drug candidates targeting the above genes with repurposing opportunities, FDA-approved drug-R788 and AMX served as effective TBK1 inhibitors. The subsequent pharmacological studies validated the potential of R788 and AMX for treating specific ALS subtypes through the inhibition of TBK1 phosphorylation.

The genetic basis of autoimmunity seen through the lens of T cell functional traits

Autoimmune disease heritability is enriched in T cell-specific regulatory regions of the genome. Modern-day T cell datasets now enable association studies between single nucleotide polymorphisms (SNPs) and a myriad of molecular phenotypes, including chromatin accessibility, gene expression, transcriptional programs, T cell antigen receptor (TCR) amino acid usage, and cell state abundances. Such studies have identified hundreds of quantitative trait loci (QTLs) in T cells that colocalize with genetic risk for autoimmune disease. The key challenge facing immunologists today lies in synthesizing these results toward a unified understanding of the autoimmune T cell: which genes, cell states, and antigens drive tissue destruction?

Transgressive segregation, hopeful monsters, and phenotypic selection drove rapid genetic gains and breakthroughs in predictive breeding for quantitative resistance to Macrophomina in strawberry

Two decades have passed since the strawberry (Fragaria x ananassa) disease caused by Macrophomina phaseolina, a necrotrophic soilborne fungal pathogen, began surfacing in California, Florida, and elsewhere. This disease has since become one of the most common causes of plant death and yield losses in strawberry. The Macrophomina problem emerged and expanded in the wake of the global phase-out of soil fumigation with methyl bromide and appears to have been aggravated by an increase in climate change-associated abiotic stresses. Here we show that sources of resistance to this pathogen are rare in gene banks and that the favorable alleles they carry are phenotypically unobvious. The latter were exposed by transgressive segregation and selection in populations phenotyped for resistance to Macrophomina under heat and drought stress. The genetic gains were immediate and dramatic. The frequency of highly resistant individuals increased from 1% in selection cycle 0 to 74% in selection cycle 2. Using GWAS and survival analysis, we found that phenotypic selection had increased the frequencies of favorable alleles among 10 loci associated with resistance and that favorable alleles had to be accumulated among four or more of these loci for an individual to acquire resistance. An unexpectedly straightforward solution to the Macrophomina disease resistance breeding problem emerged from our studies, which showed that highly resistant cultivars can be developed by genomic selection per se or marker-assisted stacking of favorable alleles among a comparatively small number of large-effect loci.

Molecular quantitative trait loci in reproductive tissues impact male fertility in cattle

Breeding bulls are well suited to investigate inherited variation in male fertility because they are genotyped and their reproductive success is monitored through semen analyses and thousands of artificial inseminations. However, functional data from relevant tissues are lacking in cattle, which prevents fine-mapping fertility-associated genomic regions. Here, we characterize gene expression and splicing variation in testis, epididymis, and vas deferens transcriptomes of 118 mature bulls and conduct association tests between 414,667 molecular phenotypes and 21,501,032 genome-wide variants to identify 41,156 regulatory loci. We show broad consensus in tissue-specific and tissue-enriched gene expression between the three bovine tissues and their human and murine counterparts. Expression- and splicing-mediating variants are more than three times as frequent in testis than epididymis and vas deferens, highlighting the transcriptional complexity of testis. Finally, we identify genes (WDR19, SPATA16, KCTD19, ZDHHC1) and molecular phenotypes that are associated with quantitative variation in male fertility through transcriptome-wide association and colocalization analyses.

Identification of eQTLs using different sets of single nucleotide polymorphisms associated with carcass and body composition traits in pigs

BACKGROUND

Mapping expression quantitative trait loci (eQTLs) in skeletal muscle tissue in pigs is crucial for understanding the relationship between genetic variation and phenotypic expression of carcass traits in meat animals. Therefore, the primary objective of this study was to evaluate the impact of different sets of single nucleotide polymorphisms (SNP), including scenarios removing SNPs pruned for linkage disequilibrium (LD) and SNPs derived from SNP chip arrays and RNA-seq data from liver, brain, and skeletal muscle tissues, on the identification of eQTLs in the Longissimus lumborum tissue, associated with carcass and body composition traits in Large White pigs. The SNPs identified from muscle mRNA were combined with SNPs identified in the brain and liver tissue transcriptomes, as well as SNPs from the GGP Porcine 50 K SNP chip array. Cis- and trans-eQTLs were identified based on the skeletal muscle gene expression level, followed by functional genomic analyses and statistical associations with carcass and body composition traits in Large White pigs.

RESULTS

The number of cis- and trans-eQTLs identified across different sets of SNPs (scenarios) ranged from 261 to 2,539 and from 29 to 13,721, respectively. Furthermore, 6,180 genes were modulated by eQTLs in at least one of the scenarios evaluated. The eQTLs identified were not significantly associated with carcass and body composition traits but were significantly enriched for many traits in the "Meat and Carcass" type QTL. The scenarios with the highest number of cis- (n = 304) and trans- (n = 5,993) modulated genes were the unpruned and LD-pruned SNP set scenarios identified from the muscle transcriptome. These genes include 84 transcription factor coding genes.

CONCLUSIONS

After LD pruning, the set of SNPs identified based on the transcriptome of the skeletal muscle tissue of pigs resulted in the highest number of genes modulated by eQTLs. Most eQTLs are of the trans type and are associated with genes influencing complex traits in pigs, such as transcription factors and enhancers. Furthermore, the incorporation of SNPs from other genomic regions to the set of SNPs identified in the porcine skeletal muscle transcriptome contributed to the identification of eQTLs that had not been identified based on the porcine skeletal muscle transcriptome alone.

Sources of gene expression variation in a globally diverse human cohort

Genetic variation influencing gene expression and splicing is a key source of phenotypic diversity. Though invaluable, studies investigating these links in humans have been strongly biased toward participants of European ancestries, diminishing generalizability and hindering evolutionary research. To address these limitations, we developed MAGE, an open-access RNA-seq data set of lymphoblastoid cell lines from 731 individuals from the 1000 Genomes Project spread across 5 continental groups and 26 populations. Most variation in gene expression (92%) and splicing (95%) was distributed within versus between populations, mirroring variation in DNA sequence. We mapped associations between genetic variants and expression and splicing of nearby genes (cis-eQTLs and cis-sQTLs, respective), identifying >15,000 putatively causal eQTLs and >16,000 putatively causal sQTLs that are enriched for relevant epigenomic signatures. These include 1310 eQTLs and 1657 sQTLs that are largely private to previously underrepresented populations. Our data further indicate that the magnitude and direction of causal eQTL effects are highly consistent across populations and that apparent "population-specific" effects observed in previous studies were largely driven by low resolution or additional independent eQTLs of the same genes that were not detected. Together, our study expands understanding of gene expression diversity across human populations and provides an inclusive resource for studying the evolution and function of human genomes.

eQTL mapping in fetal-like pancreatic progenitor cells reveals early developmental insights into diabetes risk

The impact of genetic regulatory variation active in early pancreatic development on adult pancreatic disease and traits is not well understood. Here, we generate a panel of 107 fetal-like iPSC-derived pancreatic progenitor cells (iPSC-PPCs) from whole genome-sequenced individuals and identify 4065 genes and 4016 isoforms whose expression and/or alternative splicing are affected by regulatory variation. We integrate eQTLs identified in adult islets and whole pancreas samples, which reveal 1805 eQTL associations that are unique to the fetal-like iPSC-PPCs and 1043 eQTLs that exhibit regulatory plasticity across the fetal-like and adult pancreas tissues. Colocalization with GWAS risk loci for pancreatic diseases and traits show that some putative causal regulatory variants are active only in the fetal-like iPSC-PPCs and likely influence disease by modulating expression of disease-associated genes in early development, while others with regulatory plasticity likely exert their effects in both the fetal and adult pancreas by modulating expression of different disease genes in the two developmental stages.

An atlas of genome-wide gene expression and metabolite associations and possible mediation effects towards body mass index

Investigating the cross talk of different omics layers is crucial to understand molecular pathomechanisms of metabolic diseases like obesity. Here, we present a large-scale association meta-analysis of genome-wide whole blood and peripheral blood mononuclear cell (PBMC) gene expressions profiled with Illumina HT12v4 microarrays and metabolite measurements from dried blood spots (DBS) characterized by targeted liquid chromatography tandem mass spectrometry (LC-MS/MS) in three large German cohort studies with up to 7706 samples. We found 37,295 associations comprising 72 amino acids (AA) and acylcarnitine (AC) metabolites (including ratios) and 8579 transcripts. We applied this catalogue of associations to investigate the impact of associating transcript-metabolite pairs on body mass index (BMI) as an example metabolic trait. This is achieved by conducting a comprehensive mediation analysis considering metabolites as mediators of gene expression effects and vice versa. We discovered large mediation networks comprising 27,023 potential mediation effects within 20,507 transcript-metabolite pairs. Resulting networks of highly connected (hub) transcripts and metabolites were leveraged to gain mechanistic insights into metabolic signaling pathways. In conclusion, here, we present the largest available multi-omics integration of genome-wide transcriptome data and metabolite data of amino acid and fatty acid metabolism and further leverage these findings to characterize potential mediation effects towards BMI proposing candidate mechanisms of obesity and related metabolic diseases. KEY MESSAGES: Thousands of associations of 72 amino acid and acylcarnitine metabolites and 8579 genes expand the knowledge of metabolome-transcriptome associations. A mediation analysis of effects on body mass index revealed large mediation networks of thousands of obesity-related gene-metabolite pairs. Highly connected, potentially mediating hub genes and metabolites enabled insight into obesity and related metabolic disease pathomechanisms.

Smoking-informed methylation and expression QTLs in human brain and colocalization with smoking-associated genetic loci

Smoking is a leading cause of preventable morbidity and mortality. Smoking is heritable, and genome-wide association studies (GWAS) of smoking behaviors have identified hundreds of significant loci. Most GWAS-identified variants are noncoding with unknown neurobiological effects. We used genome-wide genotype, DNA methylation, and RNA sequencing data in postmortem human nucleus accumbens (NAc) to identify cis-methylation/expression quantitative trait loci (meQTLs/eQTLs), investigate variant-by-cigarette smoking interactions across the genome, and overlay QTL evidence at smoking GWAS-identified loci to evaluate their regulatory potential. Active smokers (N=52) and nonsmokers (N=171) were defined based on cotinine biomarker levels and next-of-kin reporting. We simultaneously tested variant and variant-by-smoking interaction effects on methylation and expression, separately, adjusting for biological and technical covariates and using a two-stage multiple testing approach with eigenMT and Bonferroni corrections. We found >2 million significant meQTL variants (padj<0.05) corresponding to 41,695 unique CpGs. Results were largely driven by main effects; five meQTLs, mapping to NUDT12, FAM53B, RNF39, and ADRA1B, showed a significant interaction with smoking. We found 57,683 significant eQTLs for 958 unique eGenes (padj<0.05) and no smoking interactions. Colocalization analyses identified loci with smoking-associated GWAS variants that overlapped meQTLs/eQTLs, suggesting that these heritable factors may influence smoking behaviors through functional effects on methylation/expression. One locus containing MUSTIN1 and ITIH4 colocalized across all data types (GWAS + meQTL + eQTL). In this first genome-wide meQTL map in the human NAc, the enriched overlap with smoking GWAS-identified genetic loci provides evidence that gene regulation in the brain helps explain the neurobiology of smoking behaviors.

Master regulator activity QTL protocol to implicate regulatory pathways potentially mediating GWAS signals using eQTL data

Here, we present a protocol to identify transcriptional regulators potentially mediating downstream biological effects of germline variants associated with complex traits of interest, which enables functional hypothesis generation independent of colocalizing expression quantitative trait loci (eQTLs). We describe steps for tissue-/cell-type-specific co-expression network modeling, expression regulator activity inference, and identification of representative phenotypic master regulators. Finally, we detail activity QTL and eQTL analyses. This protocol requires genotype, expression, and relevant covariables and phenotype data from existing eQTL datasets. For complete details on the use and execution of this protocol, please refer to Hoskins et al.1.

Genetics of cell-type-specific post-transcriptional gene regulation during human neurogenesis

The function of some genetic variants associated with brain-relevant traits has been explained through colocalization with expression quantitative trait loci (eQTL) conducted in bulk post-mortem adult brain tissue. However, many brain-trait associated loci have unknown cellular or molecular function. These genetic variants may exert context-specific function on different molecular phenotypes including post-transcriptional changes. Here, we identified genetic regulation of RNA-editing and alternative polyadenylation (APA), within a cell-type-specific population of human neural progenitors and neurons. More RNA-editing and isoforms utilizing longer polyadenylation sequences were observed in neurons, likely due to higher expression of genes encoding the proteins mediating these post-transcriptional events. We also detected hundreds of cell-type-specific editing quantitative trait loci (edQTLs) and alternative polyadenylation QTLs (apaQTLs). We found colocalizations of a neuron edQTL in CCDC88A with educational attainment and a progenitor apaQTL in EP300 with schizophrenia, suggesting genetically mediated post-transcriptional regulation during brain development lead to differences in brain function.

ClipperQTL: ultrafast and powerful eGene identification method

A central task in expression quantitative trait locus (eQTL) analysis is to identify cis-eGenes (henceforth "eGenes"), i.e., genes whose expression levels are regulated by at least one local genetic variant. Among the existing eGene identification methods, FastQTL is considered the gold standard but is computationally expensive as it requires thousands of permutations for each gene. Alternative methods such as eigenMT and TreeQTL have lower power than FastQTL. In this work, we propose ClipperQTL, which reduces the number of permutations needed from thousands to 20 for data sets with large sample sizes (> 450) by using the contrastive strategy developed in Clipper; for data sets with smaller sample sizes, it uses the same permutation-based approach as FastQTL. We show that ClipperQTL performs as well as FastQTL and runs about 500 times faster if the contrastive strategy is used and 50 times faster if the conventional permutation-based approach is used. The R package ClipperQTL is available at https://github.com/heatherjzhou/ClipperQTL.

Robust identification of regulatory variants (eQTLs) using a differential expression framework developed for RNA-sequencing

BACKGROUND

A gap currently exists between genetic variants and the underlying cell and tissue biology of a trait, and expression quantitative trait loci (eQTL) studies provide important information to help close that gap. However, two concerns that arise with eQTL analyses using RNA-sequencing data are normalization of data across samples and the data not following a normal distribution. Multiple pipelines have been suggested to address this. For instance, the most recent analysis of the human and farm Genotype-Tissue Expression (GTEx) project proposes using trimmed means of M-values (TMM) to normalize the data followed by an inverse normal transformation.

RESULTS

In this study, we reasoned that eQTL analysis could be carried out using the same framework used for differential gene expression (DGE), which uses a negative binomial model, a statistical test feasible for count data. Using the GTEx framework, we identified 35 significant eQTLs (P < 5 × 10^-8) following the ANOVA model and 39 significant eQTLs (P < 5 × 10^-8) following the additive model. Using a differential gene expression framework, we identified 930 and six significant eQTLs (P < 5 × 10^-8) following an analytical framework equivalent to the ANOVA and additive model, respectively. When we compared the two approaches, there was no overlap of significant eQTLs between the two frameworks. Because we defined specific contrasts, we identified trans eQTLs that more closely resembled what we expect from genetic variants showing complete dominance between alleles. Yet, these were not identified by the GTEx framework.

CONCLUSIONS

Our results show that transforming RNA-sequencing data to fit a normal distribution prior to eQTL analysis is not required when the DGE framework is employed. Our proposed approach detected biologically relevant variants that otherwise would not have been identified due to data transformation to fit a normal distribution.

Multi-Omics Studies in Historically Excluded Populations: The Road to Equity

Over the past few decades, genomewide association studies (GWASs) have identified the specific genetics variants contributing to many complex diseases by testing millions of genetic variations across the human genome against a variety of phenotypes. However, GWASs are limited in their ability to uncover mechanistic insight given that most significant associations are found in non-coding region of the genome. Furthermore, the lack of diversity in studies has stymied the advance of precision medicine for many historically excluded populations. In this review, we summarize most popular multi-omics approaches (genomics, transcriptomics, proteomics, and metabolomics) related to precision medicine and highlight if diverse populations have been included and how their findings have advance biological understanding of disease and drug response. New methods that incorporate local ancestry have been to improve the power of GWASs for admixed populations (such as African Americans and Latinx). Because most signals from GWAS are in the non-coding region, other machine learning and omics approaches have been developed to identify the potential causative single-nucleotide polymorphisms and genes that explain these phenotypes. These include polygenic risk scores, expression quantitative trait locus mapping, and transcriptome-wide association studies. Analogous protein methods, such as proteins quantitative trait locus mapping, proteome-wide association studies, and metabolomic approaches provide insight into the consequences of genetic variation on protein abundance. Whereas, integrated multi-omics studies have improved our understanding of the mechanisms for genetic association, we still lack the datasets and cohorts for historically excluded populations to provide equity in precision medicine and pharmacogenomics.

Fine mapping spatiotemporal mechanisms of genetic variants underlying cardiac traits and disease

The causal variants and genes underlying thousands of cardiac GWAS signals have yet to be identified. Here, we leverage spatiotemporal information on 966 RNA-seq cardiac samples and perform an expression quantitative trait locus (eQTL) analysis detecting eQTLs considering both eGenes and eIsoforms. We identify 2,578 eQTLs associated with a specific developmental stage-, tissue- and/or cell type. Colocalization between eQTL and GWAS signals of five cardiac traits identified variants with high posterior probabilities for being causal in 210 GWAS loci. Pulse pressure GWAS loci are enriched for colocalization with fetal- and smooth muscle- eQTLs; pulse rate with adult- and cardiac muscle- eQTLs; and atrial fibrillation with cardiac muscle- eQTLs. Fine mapping identifies 79 credible sets with five or fewer SNPs, of which 15 were associated with spatiotemporal eQTLs. Our study shows that many cardiac GWAS variants impact traits and disease in a developmental stage-, tissue- and/or cell type-specific fashion.

Epigenetic and Genetic Population Structure is Coupled in a Marine Invertebrate

Delineating the relative influence of genotype and the environment on DNA methylation is critical for characterizing the spectrum of organism fitness as driven by adaptation and phenotypic plasticity. In this study, we integrated genomic and DNA methylation data for two distinct Olympia oyster (Ostrea lurida) populations while controlling for within-generation environmental influences. In addition to providing the first characterization of genome-wide DNA methylation patterns in the oyster genus Ostrea, we identified 3,963 differentially methylated loci between populations. Our results show a clear coupling between genetic and epigenetic patterns of variation, with 27% of variation in interindividual methylation differences explained by genotype. Underlying this association are both direct genetic changes in CpGs (CpG-SNPs) and genetic variation with indirect influence on methylation (mQTLs). When comparing measures of genetic and epigenetic population divergence at specific genomic regions this relationship surprisingly breaks down, which has implications for the methods commonly used to study epigenetic and genetic coupling in marine invertebrates.

Disentangling the aetiological pathways between body mass index and site-specific cancer risk using tissue-partitioned Mendelian randomisation

BACKGROUND

Body mass index (BMI) is known to influence the risk of various site-specific cancers, however, dissecting which subcomponents of this heterogenous risk factor are predominantly responsible for driving disease effects has proven difficult to establish. We have leveraged tissue-specific gene expression to separate the effects of distinct phenotypes underlying BMI on the risk of seven site-specific cancers.

METHODS

SNP-exposure estimates were weighted in a multivariable Mendelian randomisation analysis by their evidence for colocalization with subcutaneous adipose- and brain-tissue-derived gene expression using a recently developed methodology.

RESULTS

Our results provide evidence that brain-tissue-derived BMI variants are predominantly responsible for driving the genetically predicted effect of BMI on lung cancer (OR: 1.17; 95% CI: 1.01-1.36; P = 0.03). Similar findings were identified when analysing cigarettes per day as an outcome (Beta = 0.44; 95% CI: 0.26-0.61; P = 1.62 × 10-6), highlighting a possible shared aetiology or mediator effect between brain-tissue BMI, smoking and lung cancer. Our results additionally suggest that adipose-tissue-derived BMI variants may predominantly drive the effect of BMI and increased risk for endometrial cancer (OR: 1.71; 95% CI: 1.07-2.74; P = 0.02), highlighting a putatively important role in the aetiology of endometrial cancer.

CONCLUSIONS

The study provides valuable insight into the divergent underlying pathways between BMI and the risk of site-specific cancers.

Analysis of merged transcriptomic and genomic datasets to identify genes and pathways underlying residual feed intake in growing pigs

Improvement of feed efficiency (FE) in pigs is an important milestone in order to reduce the economic and environmental impact of pig production. The goal of finding biomarkers for FE has persisted for decades. However, due to the complexity of the FE trait, these goals have still not been met. Here, we search for quantitative trait loci (QTL), candidate genes, and biological pathways associated with FE using both genotype and RNA-seq data. We obtained genotype and colon epithelium RNA-seq data for 375 and 96 pigs, respectively. In total, a genome-wide association study (GWAS) and differential expression (DE) analysis led to detection of three QTL on SSC9 and 17 DE-genes associated with FE. Possible intersection points between genes located in QTL and DE-genes were found on levels of transcription factor-target interaction. Moreover, cis-eQTL analysis revealed associations between genotype and expression levels of three DE-genes and three genes located in the GWAS QTLs, which may establish the connection between genotype and phenotype through DE. Finally, single nucleotide polymorphism calling using RNA-seq data for genes located in GWAS QTLs revealed 53 polymorphisms of which eleven were missense variants.

A comparison of the genes and genesets identified by GWAS and EWAS of fifteen complex traits

Identifying genomic regions pertinent to complex traits is a common goal of genome-wide and epigenome-wide association studies (GWAS and EWAS). GWAS identify causal genetic variants, directly or via linkage disequilibrium, and EWAS identify variation in DNA methylation associated with a trait. While GWAS in principle will only detect variants due to causal genes, EWAS can also identify genes via confounding, or reverse causation. We systematically compare GWAS (N > 50,000) and EWAS (N > 4500) results of 15 complex traits. We evaluate if the genes or gene ontology terms flagged by GWAS and EWAS overlap, and find substantial overlap for diastolic blood pressure, (gene overlap P = 5.2 × 10^-6; term overlap P = 0.001). We superimpose our empirical findings against simulated models of varying genetic and epigenetic architectures and observe that in most cases GWAS and EWAS are likely capturing distinct genesets. Our results indicate that GWAS and EWAS are capturing different aspects of the biology of complex traits.

Current challenges in understanding the role of enhancers in disease

Enhancers play a central role in the spatiotemporal control of gene expression and tend to work in a cell-type-specific manner. In addition, they are suggested to be major contributors to phenotypic variation, evolution and disease. There is growing evidence that enhancer dysfunction due to genetic, structural or epigenetic mechanisms contributes to a broad range of human diseases referred to as enhanceropathies. Such mechanisms often underlie the susceptibility to common diseases, but can also play a direct causal role in cancer or Mendelian diseases. Despite the recent gain of insights into enhancer biology and function, we still have a limited ability to predict how enhancer dysfunction impacts gene expression. Here we discuss the major challenges that need to be overcome when studying the role of enhancers in disease etiology and highlight opportunities and directions for future studies, aiming to disentangle the molecular basis of enhanceropathies.

The Awesome Power of Human Genetics of Infectious Disease

Since the identification of sickle cell trait as a heritable form of resistance to malaria, candidate gene studies, linkage analysis paired with sequencing, and genome-wide association (GWA) studies have revealed many examples of genetic resistance and susceptibility to infectious diseases. GWA studies enabled the identification of many common variants associated with small shifts in susceptibility to infectious diseases. This is exemplified by multiple loci associated with leprosy, malaria, HIV, tuberculosis, and coronavirus disease 2019 (COVID-19), which illuminate genetic architecture and implicate pathways underlying pathophysiology. Despite these successes, most of the heritability of infectious diseases remains to be explained. As the field advances, current limitations may be overcome by applying methodological innovations such as cellular GWA studies and phenome-wide association (PheWA) studies as well as by improving methodological rigor with more precise case definitions, deeper phenotyping, increased cohort diversity, and functional validation of candidate loci in the laboratory or human challenge studies.

Mapping responsive genomic elements to heat stress in a maize diversity panel

BACKGROUND

Many plant species exhibit genetic variation for coping with environmental stress. However, there are still limited approaches to effectively uncover the genomic region that regulates distinct responsive patterns of the gene across multiple varieties within the same species under abiotic stress.

RESULTS

By analyzing the transcriptomes of more than 100 maize inbreds, we reveal many cis- and trans-acting eQTLs that influence the expression response to heat stress. The cis-acting eQTLs in response to heat stress are identified in genes with differential responses to heat stress between genotypes as well as genes that are only expressed under heat stress. The cis-acting variants for heat stress-responsive expression likely result from distinct promoter activities, and the differential heat responses of the alleles are confirmed for selected genes using transient expression assays. Global footprinting of transcription factor binding is performed in control and heat stress conditions to document regions with heat-enriched transcription factor binding occupancies.

CONCLUSIONS

Footprints enriched near proximal regions of characterized heat-responsive genes in a large association panel can be utilized for prioritizing functional genomic regions that regulate genotype-specific responses under heat stress.

Ecological Speciation Promoted by Divergent Regulation of Functional Genes Within African Cichlid Fishes

Rapid ecological speciation along depth gradients has taken place repeatedly in freshwater fishes, yet molecular mechanisms facilitating such diversification are typically unclear. In Lake Masoko, an African crater lake, the cichlid Astatotilapia calliptera has diverged into shallow-littoral and deep-benthic ecomorphs with strikingly different jaw structures within the last 1,000 years. Using genome-wide transcriptome data, we explore two major regulatory transcriptional mechanisms, expression and splicing-QTL variants, and examine their contributions to differential gene expression underpinning functional phenotypes. We identified 7,550 genes with significant differential expression between ecomorphs, of which 5.4% were regulated by cis-regulatory expression QTLs, and 9.2% were regulated by cis-regulatory splicing QTLs. We also found strong signals of divergent selection on differentially expressed genes associated with craniofacial development. These results suggest that large-scale transcriptome modification plays an important role during early-stage speciation. We conclude that regulatory variants are important targets of selection driving ecologically relevant divergence in gene expression during adaptive diversification.

Evaluating molecular fingerprint-based models of drug side effects against a statistical control

There are many machine learning models that use molecular fingerprints of drugs to predict side effects. Characterizing their skill is necessary for understanding their usefulness in pharmaceutical development. Here, we analyze a statistical control of side effect prediction skill, develop a pipeline for benchmarking models, and evaluate how well existing models predict side effects identified in pharmaceutical documentation. We demonstrate that molecular fingerprints are useful for ranking drugs by their likelihood to cause a given side effect. However, the predictions for one or more drugs overall benefit only marginally from molecular fingerprints when ranking the likelihoods of many possible side effects, and display at most modest overall skill at identifying the side effects that do and do not occur.

Genetic susceptibility loci for Chlamydia trachomatis endometrial infection influence expression of genes involved in T cell function, tryptophan metabolism and epithelial integrity

Objectives

Identify genetic loci of enhanced susceptibility to Chlamydial trachomatis (Ct) upper genital tract infection in women.

Methods

We performed an integrated analysis of DNA genotypes and blood-derived mRNA profiles from 200 Ct-exposed women to identify expression quantitative trait loci (eQTL) and determine their association with endometrial chlamydial infection using a mediation test. We further evaluated the effect of a lead eQTL on the expression of CD151 by immune cells from women with genotypes associated with low and high whole blood expression of CD151, respectively.

Results

We identified cis-eQTLs modulating mRNA expression of 81 genes (eGenes) associated with altered risk of ascending infection. In women with endometrial infection, eGenes involved in proinflammatory signaling were upregulated. Downregulated eGenes included genes involved in T cell functions pivotal for chlamydial control. eGenes encoding molecules linked to metabolism of tryptophan, an essential chlamydial nutrient, and formation of epithelial tight junctions were also downregulated in women with endometrial infection. A lead eSNP rs10902226 was identified regulating CD151, a tetrospanin molecule important for immune cell adhesion and migration and T cell proliferation. Further in vitro experiments showed that women with a CC genotype at rs10902226 had reduced rates of endometrial infection with increased CD151 expression in whole blood and T cells when compared to women with a GG genotype.

Conclusions

We discovered genetic variants associated with altered risk for Ct ascension. A lead eSNP for CD151 is a candidate genetic marker for enhanced CD4 T cell function and reduced susceptibility.

Effects of Coronary Artery Disease-Associated Variants on Vascular Smooth Muscle Cells

BACKGROUND

Genome-wide association studies have identified many genetic loci that are robustly associated with coronary artery disease (CAD). However, the underlying biological mechanisms are still unknown for most of these loci, hindering the progress to medical translation. Evidence suggests that the genetic influence on CAD susceptibility may act partly through vascular smooth muscle cells (VSMCs).

METHODS

We undertook genotyping, RNA sequencing, and cell behavior assays on a large bank of VSMCs (n>1499). Expression quantitative trait locus and splicing quantitative trait locus analyses were performed to identify genes with an expression that was influenced by CAD-associated variants. To identify candidate causal genes for CAD, we ascertained colocalizations of VSMC expression quantitative trait locus signals with CAD association signals by performing causal variants identification in associated regions analysis and the summary data-based mendelian randomization test. Druggability analysis was then performed on the candidate causal genes. CAD risk variants were tested for associations with VSMC proliferation, migration, and apoptosis. Collective effects of multiple CAD-associated variants on VSMC behavior were estimated by polygenic scores.

RESULTS

Approximately 60% of the known CAD-associated variants showed statistically significant expression quantitative trait locus or splicing quantitative trait locus effects in VSMCs. Colocalization analyses identified 84 genes with expression quantitative trait locus signals that significantly colocalized with CAD association signals, identifying them as candidate causal genes. Druggability analysis indicated that 38 of the candidate causal genes were druggable, and 13 had evidence of drug-gene interactions. Of the CAD-associated variants tested, 139 showed suggestive associations with VSMC proliferation, migration, or apoptosis. A polygenic score model explained up to 5.94% of variation in several VSMC behavior parameters, consistent with polygenic influences on VSMC behavior.

CONCLUSIONS

This comprehensive analysis shows that a large percentage of CAD loci can modulate gene expression in VSMCs and influence VSMC behavior. Several candidate causal genes identified are likely to be druggable and thus represent potential therapeutic targets.

Deleterious mutations and the rare allele burden on rice gene expression

Deleterious genetic variation is maintained in populations at low frequencies. Under a model of stabilizing selection, rare (and presumably deleterious) genetic variants are associated with increase or decrease in gene expression from some intermediate optimum. We investigate this phenomenon in a population of largely Oryza sativa ssp. indica rice landraces under normal unstressed wet and stressful drought field conditions. We include single nucleotide polymorphisms, insertion/deletion mutations, and structural variants in our analysis and find a stronger association between rare variants and gene expression outliers under the stress condition. We also show an association of the strength of this rare variant effect with linkage, gene expression levels, network connectivity, local recombination rate, and fitness consequence scores, consistent with the stabilizing selection model of gene expression.