Mendelian randomization integrating GWAS and eQTL data reveals genetic determinants of complex and clinical traits

Disentangling mechanisms behind the pleiotropic effects of proximal 16p11.2 BP4-5 CNVs

Whereas 16p11.2 BP4-5 copy-number variants (CNVs) represent one of the most pleiotropic etiologies of genomic syndromes in both clinical and population cohorts, the mechanisms leading to such pleiotropy remain understudied. Identifying 73 deletion and 89 duplication carrier individuals among unrelated White British UK Biobank participants, we performed a phenome-wide association study (PheWAS) between the region's copy number and 117 complex traits and diseases, mimicking four dosage models. Forty-six phenotypes (39%) were affected by 16p11.2 BP4-5 CNVs, with the deletion-only, mirror, U-shape, and duplication-only models being the best fit for 30, 10, 4, and 2 phenotypes, respectively, aligning with the stronger deleteriousness of the deletion. Upon individually adjusting CNV effects for either body mass index (BMI), height, or educational attainment (EA), we found that sixteen testable deletion-driven associations-primarily with cardiovascular and metabolic traits-were BMI dependent, with EA playing a more subtle role and no association depending on height. Bidirectional Mendelian randomization supported that 13 out of these 16 associations were secondary consequences of the CNV's impact on BMI. For the 23 traits that remained significantly associated upon individual adjustment for mediators, matched-control analyses found that 10 phenotypes, including musculoskeletal traits, liver enzymes, fluid intelligence, platelet count, and pneumonia and acute kidney injury risk, remained associated under strict Bonferroni correction, with 10 additional nominally significant associations. These results paint a complex picture of 16p11.2 BP4-5's pleiotropic pattern that involves direct effects on multiple physiological systems and indirect co-morbidities consequential to the CNV's impact on BMI and EA, acting through trait-specific dosage mechanisms.

Prediction of causal genes at GWAS loci with pleiotropic gene regulatory effects using sets of correlated instrumental variables

Multivariate Mendelian randomization (MVMR) is a statistical technique that uses sets of genetic instruments to estimate the direct causal effects of multiple exposures on an outcome of interest. At genomic loci with pleiotropic gene regulatory effects, that is, loci where the same genetic variants are associated to multiple nearby genes, MVMR can potentially be used to predict candidate causal genes. However, consensus in the field dictates that the genetic instruments in MVMR must be independent (not in linkage disequilibrium, which is usually not possible when considering a group of candidate genes from the same locus. Here we used causal inference theory to show that MVMR with correlated instruments satisfies the instrumental set condition. This is a classical result by Brito and Pearl (2002) for structural equation models that guarantees the identifiability of individual causal effects in situations where multiple exposures collectively, but not individually, separate a set of instrumental variables from an outcome variable. Extensive simulations confirmed the validity and usefulness of these theoretical results. Importantly, the causal effect estimates remained unbiased and their variance small even when instruments are highly correlated, while bias introduced by horizontal pleiotropy or LD matrix sampling error was comparable to standard MR. We applied MVMR with correlated instrumental variable sets at genome-wide significant loci for coronary artery disease (CAD) risk using expression Quantitative Trait Loci (eQTL) data from seven vascular and metabolic tissues in the STARNET study. Our method predicts causal genes at twelve loci, each associated with multiple colocated genes in multiple tissues. We confirm causal roles for PHACTR 1 and ADAMTS 7 in arterial tissues, among others. However, the extensive degree of regulatory pleiotropy across tissues and the limited number of causal variants in each locus still require that MVMR is run on a tissue-by-tissue basis, and testing all gene-tissue pairs with cis-eQTL associations at a given locus in a single model to predict causal gene-tissue combinations remains infeasible. Our results show that within tissues, MVMR with dependent, as opposed to independent, sets of instrumental variables significantly expands the scope for predicting causal genes in disease risk loci with pleiotropic regulatory effects. However, considering risk loci with regulatory pleiotropy that also spans across tissues remains an unsolved problem.

The causal relationship between triglycerides and myocardial infarction: A two-sample Mendelian randomization

The causal relationship between triglycerides and myocardial infarction (MI) was investigated using Mendelian randomization (MR) studies. Triglycerides were the exposure factor, and MI served as the outcome variable. Inverse variance weighting was used as the main analysis method, MR-Egger, and weight median as other analysis methods for MR analysis. In addition, heterogeneity test, level multivariate analysis, and sensitivity analysis were carried out. Inverse variance weighting results showed that the increase in triglyceride level affected the incidence of MI (OR = 1.287; 95% CI = 1.185-1.398; P = 1.988 × 10-9). Consistently, the results from all 3 methods indicated a statistically significant increase in the risk of MI with higher triglyceride levels (P < .05). The results showed that patients with high triglyceride levels had a higher incidence of MI, suggesting that MI should be prevented in the high triglyceride population.

Host genetics and gut microbiota synergistically regulate feed utilization in egg-type chickens

BACKGROUND

Feed efficiency is a crucial economic trait in poultry industry. Both host genetics and gut microbiota influence feed efficiency. However, the associations between gut microbiota and host genetics, as well as their combined contributions to feed efficiency in laying hens during the late laying period, remain largely unclear.

METHODS

In total, 686 laying hens were used for whole-genome resequencing and liver transcriptome sequencing. 16S rRNA gene sequencing was conducted on gut chyme (duodenum, jejunum, ileum, and cecum) and fecal samples from 705 individuals. Bioinformatic analysis was performed by integrating the genome, transcriptome, and microbiome to screen for key genetic variations, genes, and gut microbiota associated with feed efficiency.

RESULTS

The heritability of feed conversion ratio (FCR) and residual feed intake (RFI) was determined to be 0.28 and 0.48, respectively. The ileal and fecal microbiota accounted for 15% and 10% of the FCR variance, while the jejunal, cecal, and fecal microbiota accounted for 20%, 11%, and 10% of the RFI variance. Through SMR analysis based on summary data from liver eQTL mapping and GWAS, we further identified four protein-coding genes, SUCLA2, TNFSF13B, SERTM1, and MARVELD3, that influence feed efficiency in laying hens. The SUCLA2 and TNFSF13B genes were significantly associated with SNP 1:25664581 and SNP rs312433097, respectively. SERTM1 showed significant associations with rs730958360 and 1:33542680 and is a potential causal gene associated with the abundance of Corynebacteriaceae in feces. MARVELD3 was significantly associated with the 1:135348198 and was significantly correlated with the abundance of Enterococcus in ileum. Specifically, a lower abundance of Enterococcus in ileum and a higher abundance of Corynebacteriaceae in feces were associated with better feed efficiency.

CONCLUSIONS

This study confirms that both host genetics and gut microbiota can drive variations in feed efficiency. A small portion of the gut microbiota often interacts with host genes, collectively enhancing feed efficiency. Therefore, targeting both the gut microbiota and host genetic variation by supporting more efficient taxa and selective breeding could improve feed efficiency in laying hens during the late laying period.

Dissecting the causal relationship between moderate to vigorous physical activity levels and cognitive performance: a bidirectional two-sample Mendelian randomization study

Introduction

Recent studies increasingly suggest that moderate-to-vigorous physical activity (MVPA) impacts cognitive risk. However, the bidirectional nature of this relationship warrants further exploration. To address this, we employed a Mendelian randomization (MR) approach, analyzing two distinct samples.

Methods

These analyses utilized published genome-wide association study (GWAS) summary statistics for MVPA (n = 377,234) and cognitive performance (n = 257,841). Our primary method was the inverse variance weighted (IVW) model with random effects, aiming to deduce potential causal links. Additionally, we employed supplementary methods, including MR Egger regression, Weighted median, Weighted mode, and Simple mode. For sensitivity analysis, tools like the MR Egger test, Cochran's Q, MR PRESSO, and leave-one-out (LOO) were utilized.

Results

Our findings indicate a decrease in cognitive risk with increased MVPA (Odds Ratio [OR] = 0.577, 95% Confidence Interval [CI]: 0.460-0.723, p = 1.930 × 10-6). Furthermore, enhanced cognitive levels corresponded to a reduced risk of inadequate MVPA (OR = 0.866, 95% CI: 0.839-0.895, p = 1.200 × 10-18).

Discussion

In summary, our study demonstrates that MVPA lowers cognitive risk, while poor cognitive health may impede participation in MVPA. Overall, these findings provide valuable insights for developing personalized prevention and intervention strategies in health and sports sciences.

The goldmine of GWAS summary statistics: a systematic review of methods and tools

Genome-wide association studies (GWAS) have revolutionized our understanding of the genetic architecture of complex traits and diseases. GWAS summary statistics have become essential tools for various genetic analyses, including meta-analysis, fine-mapping, and risk prediction. However, the increasing number of GWAS summary statistics and the diversity of software tools available for their analysis can make it challenging for researchers to select the most appropriate tools for their specific needs. This systematic review aims to provide a comprehensive overview of the currently available software tools and databases for GWAS summary statistics analysis. We conducted a comprehensive literature search to identify relevant software tools and databases. We categorized the tools and databases by their functionality, including data management, quality control, single-trait analysis, and multiple-trait analysis. We also compared the tools and databases based on their features, limitations, and user-friendliness. Our review identified a total of 305 functioning software tools and databases dedicated to GWAS summary statistics, each with unique strengths and limitations. We provide descriptions of the key features of each tool and database, including their input/output formats, data types, and computational requirements. We also discuss the overall usability and applicability of each tool for different research scenarios. This comprehensive review will serve as a valuable resource for researchers who are interested in using GWAS summary statistics to investigate the genetic basis of complex traits and diseases. By providing a detailed overview of the available tools and databases, we aim to facilitate informed tool selection and maximize the effectiveness of GWAS summary statistics analysis.

Extensive co-regulation of neighboring genes complicates the use of eQTLs in target gene prioritization

Identifying causal genes underlying genome-wide association studies (GWASs) is a fundamental problem in human genetics. Although colocalization with gene expression quantitative trait loci (eQTLs) is often used to prioritize GWAS target genes, systematic benchmarking has been limited due to unavailability of large ground truth datasets. Here, we re-analyzed plasma protein QTL data from 3,301 individuals of the INTERVAL cohort together with 131 eQTL Catalog datasets. Focusing on variants located within or close to the affected protein identified 793 proteins with at least one cis-pQTL where we could assume that the most likely causal gene was the gene coding for the protein. We then benchmarked the ability of cis-eQTLs to recover these causal genes by comparing three Bayesian colocalization methods (coloc.susie, coloc.abf, and CLPP) and five Mendelian randomization (MR) approaches (three varieties of inverse-variance weighted MR, MR-RAPS, and MRLocus). We found that assigning fine-mapped pQTLs to their closest protein coding genes outperformed all colocalization methods regarding both precision (71.9%) and recall (76.9%). Furthermore, the colocalization method with the highest recall (coloc.susie - 46.3%) also had the lowest precision (45.1%). Combining evidence from multiple conditionally distinct colocalizing QTLs with MR increased precision to 81%, but this was accompanied by a large reduction in recall to 7.1%. Furthermore, the choice of the MR method greatly affected performance, with the standard inverse-variance-weighted MR often producing many false positives. Our results highlight that linking GWAS variants to target genes remains challenging with eQTL evidence alone, and prioritizing novel targets requires triangulation of evidence from multiple sources.

The importance of functional genomics studies in precision rheumatology

Rheumatic diseases, those that affect the musculoskeletal system, cause significant morbidity. Among risk factors of these diseases is a significant genetic component. Recent advances in high-throughput omics techniques now allow a comprehensive profiling of patients at a genetic level through genome-wide association studies. Without functional interpretation of variants identified through these studies, clinical insight remains limited. Strategies include statistical fine-mapping that refine the list of variants in loci associated with disease, whilst colocalization techniques attempt to attribute function to variants that overlap a genetically active chromatin annotation. Functional validation using genome editing techniques can be used to further refine genetic signals and identify key pathways in cell types relevant to rheumatic disease biology. Insight gained from the combination of genetic studies and functional validation can be used to improve precision medicine in rheumatic diseases by allowing risk prediction and drug repositioning.

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.

Investigating the causal relationship between gut microbiota and gastroenteropancreatic neuroendocrine neoplasms: a bidirectional Mendelian randomization study

Background

The interaction between the intestinal flora and gastroenteropancreatic neuroendocrine neoplasms (GEP-NENs) remains poorly understood, despite the known effect of the gut microbiota on gastrointestinal adenocarcinomas. Hence, the present research aimed to determine the potential causal correlation between the intestinal flora and GEP-NENs by conducting a bidirectional Mendelian randomization (MR) analysis.

Methods

Two-sample MR analysis was conducted using the summary statistics of the gut microbiota from the MiBioGen consortium and those of GEP-NENs from the FinnGen research project. The inverse-variance weighted approach was utilized as the primary analytical method. To enhance the robustness of our findings, multiple sensitivity tests were performed, including Cochran's Q test for evaluating heterogeneity, the MR-Egger intercept test to detect horizontal pleiotropy, and the MR-PRESSO test to identify outliers and assess pleiotropy bias. Additionally, a leave-one-out analysis was performed to validate the consistency of our findings. The MR-Steiger test was also utilized to determine the causal direction in the correlation between the gut microbiota and GEP-NENs. Finally, a reverse MR analysis was performed to assess reverse causality between the intestinal flora and GEP-NENs.

Results

We identified 42 taxa of the gut microbiota that were potentially causally associated with GEP-NENs; of these taxa, 7, 8, 11, and 16 taxa were causally associated with pancreatic NENs, colorectal NENs, small intestinal NENs, and gastric NENs, respectively. After adjusting for false discovery rate (FDR) correction, we found significant causal links of Euryarchaeota with small intestinal NENs and Family XIII UCG-001 with gastric NENs. The sensitivity analyses confirmed the stability of these correlations. In the reverse MR analysis, colorectal NENs and small intestinal NENs were found to be associated with variations in 8 and 6 different taxa of the gut microbiota, respectively. After adjusting for FDR correction, no significant causal links were detected between GEP-NENs and the intestinal flora.

Conclusion

The present study reveals a potential causal association between certain taxa of the intestinal flora and GEP-NENs, thus providing new perspectives regarding the role of the intestinal flora in the development of these tumors. These insights could provide innovative approaches to screen and prevent these diseases.

A transcriptomic atlas of the human brain reveals genetically determined aspects of neuropsychiatric health

Regulation of gene expression is a vital component of neurological homeostasis. Cataloging the consequences of endogenous gene expression on the physical structure and connectivity of the brain offers a means of unifying trait-associated genetic variation with trait-associated neurological features. We perform tissue-specific transcriptome-wide association studies (TWASs) on over 3,400 neuroimaging phenotypes in the UK Biobank (N = 33,224) using our joint-tissue imputation (JTI)-TWAS method. We identify highly significant associations between predicted expression for 7,192 genes and a wide variety of measures of the brain derived from magnetic resonance imaging (MRI). Our approach generates reproducible results in internal and external replication datasets. Genetically determined expression alone is sufficient for high-fidelity reconstruction of brain structure and organization. We demonstrate complementary benefits of cross-tissue and single-tissue analyses toward an integrated neurobiology and provide evidence that gene expression outside the central nervous system provides unique insights into brain health. As an application, we provide evidence suggesting that the genetically regulated expression of schizophrenia risk genes causally affects over 73% of neurological phenotypes that are altered in individuals with schizophrenia (as identified by neuroimaging studies). Imaging features associated with neuropsychiatric traits can provide valuable insights into underlying pathophysiology. By linking neuroimaging-derived phenotypes with expression levels of specific genes, this resource represents a powerful gene prioritization schema that can improve our understanding of brain function, development, and disease. The use of multiple different cortical and subcortical atlases in the resource facilitates direct integration of these data with findings from a diverse range of clinical neuroimaging studies.

pyTWMR: transcriptome-wide Mendelian randomization in python

MOTIVATION

Mendelian randomization (MR) is a widely used approach to estimate causal effect of variation in gene expression on complex traits. Among several MR-based algorithms, transcriptome-wide summary statistics-based Mendelian Randomization approach (TWMR) enables the uses of multiple SNPs as instruments and multiple gene expression traits as exposures to facilitate causal inference in observational studies.

RESULTS

Here we present a Python-based implementation of TWMR and revTWMR. Our implementation offers GPU computational support for faster computations and robust computation mode resilient to highly correlated gene expressions and genetic variants.

AVAILABILITY AND IMPLEMENTATION

pyTWMR is available at github.com/soreshkov/pyTWMR.

Single nucleus RNA-sequencing integrated into risk variant colocalization discovers 17 cell-type-specific abdominal obesity genes for metabolic dysfunction-associated steatotic liver disease

BACKGROUND

Abdominal obesity increases the risk for non-alcoholic fatty liver disease (NAFLD), now known as metabolic dysfunction-associated steatotic liver disease (MASLD).

METHODS

To elucidate the directional cell-type level biological mechanisms underlying the association between abdominal obesity and MASLD, we integrated adipose and liver single nucleus RNA-sequencing and bulk cis-expression quantitative trait locus (eQTL) data with the UK Biobank genome-wide association study (GWAS) data using colocalization. Then we used colocalized cis-eQTL variants as instrumental variables in Mendelian randomization (MR) analyses, followed by functional validation experiments on the target genes of the cis-eQTL variants.

FINDINGS

We identified 17 colocalized abdominal obesity GWAS variants, regulating 17 adipose cell-type marker genes. Incorporating these 17 variants into MR discovers a putative tissue-of-origin, cell-type-aware causal effect of abdominal obesity on MASLD consistently with multiple MR methods without significant evidence for pleiotropy or heterogeneity. Single cell data confirm the adipocyte-enriched mean expression of the 17 genes. Our cellular experiments across human adipogenesis identify risk variant -specific epigenetic and transcriptional mechanisms. Knocking down two of the 17 genes, PPP2R5A and SH3PXD2B, shows a marked decrease in adipocyte lipidation and significantly alters adipocyte function and adipogenesis regulator genes, including DGAT2, LPL, ADIPOQ, PPARG, and SREBF1. Furthermore, the 17 genes capture a characteristic MASLD expression signature in subcutaneous adipose tissue.

INTERPRETATION

Overall, we discover a significant cell-type level effect of abdominal obesity on MASLD and trace its biological effect to adipogenesis.

FUNDING

NIH grants R01HG010505, R01DK132775, and R01HL170604; the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program (Grant No. 802825), Academy of Finland (Grants Nos. 333021), the Finnish Foundation for Cardiovascular Research the Sigrid Jusélius Foundation and the Jane and Aatos Erkko Foundation; American Association for the Study of Liver Diseases (AASLD) Advanced Transplant Hepatology award and NIH/NIDDK (P30DK41301) Pilot and Feasibility award; NIH/NIEHS F32 award (F32ES034668); Finnish Diabetes Research Foundation, Kuopio University Hospital Project grant (EVO/VTR grants 2005-2021), the Academy of Finland grant (Contract no. 138006); Academy of Finland (Grant Nos 335443, 314383, 272376 and 266286), Sigrid Jusélius Foundation, Finnish Medical Foundation, Finnish Diabetes Research Foundation, Novo Nordisk Foundation (#NNF20OC0060547, NNF17OC0027232, NNF10OC1013354) and Government Research Funds to Helsinki University Hospital; Orion Research Foundation, Maud Kuistila Foundation, Finish Medical Foundation, and University of Helsinki.

From genetic associations to genes: methods, applications, and challenges

Genome-wide association studies (GWASs) have identified numerous genetic loci associated with human traits and diseases. However, pinpointing the causal genes remains a challenge, which impedes the translation of GWAS findings into biological insights and medical applications. In this review, we provide an in-depth overview of the methods and technologies used for prioritizing genes from GWAS loci, including gene-based association tests, integrative analysis of GWAS and molecular quantitative trait loci (xQTL) data, linking GWAS variants to target genes through enhancer-gene connection maps, and network-based prioritization. We also outline strategies for generating context-dependent xQTL data and their applications in gene prioritization. We further highlight the potential of gene prioritization in drug repurposing. Lastly, we discuss future challenges and opportunities in this field.

Identifying dysregulated regions in amyotrophic lateral sclerosis through chromatin accessibility outliers

The high heritability of amyotrophic lateral sclerosis (ALS) contrasts with its low molecular diagnosis rate post-genetic testing, pointing to potential undiscovered genetic factors. To aid the exploration of these factors, we introduced EpiOut, an algorithm to identify chromatin accessibility outliers that are regions exhibiting divergent accessibility from the population baseline in a single or few samples. Annotation of accessible regions with histone chromatin immunoprecipitation sequencing and Hi-C indicates that outliers are concentrated in functional loci, especially among promoters interacting with active enhancers. Across different omics levels, outliers are robustly replicated, and chromatin accessibility outliers are reliable predictors of gene expression outliers and aberrant protein levels. When promoter accessibility does not align with gene expression, our results indicate that molecular aberrations are more likely to be linked to post-transcriptional regulation rather than transcriptional regulation. Our findings demonstrate that the outlier detection paradigm can uncover dysregulated regions in rare diseases. EpiOut is available at github.com/uci-cbcl/EpiOut.

Metabolic gene function discovery platform GeneMAP identifies SLC25A48 as necessary for mitochondrial choline import

Organisms maintain metabolic homeostasis through the combined functions of small-molecule transporters and enzymes. While many metabolic components have been well established, a substantial number remains without identified physiological substrates. To bridge this gap, we have leveraged large-scale plasma metabolome genome-wide association studies (GWAS) to develop a multiomic Gene-Metabolite Association Prediction (GeneMAP) discovery platform. GeneMAP can generate accurate predictions and even pinpoint genes that are distant from the variants implicated by GWAS. In particular, our analysis identified solute carrier family 25 member 48 (SLC25A48) as a genetic determinant of plasma choline levels. Mechanistically, SLC25A48 loss strongly impairs mitochondrial choline import and synthesis of its downstream metabolite betaine. Integrative rare variant and polygenic score analyses in UK Biobank provide strong evidence that the SLC25A48 causal effects on human disease may in part be mediated by the effects of choline. Altogether, our study provides a discovery platform for metabolic gene function and proposes SLC25A48 as a mitochondrial choline transporter.

A systems biology-based identification and in vivo functional screening of Alzheimer's disease risk genes reveal modulators of memory function

Genome-wide association studies (GWASs) have uncovered over 75 genomic loci associated with risk for late-onset Alzheimer's disease (LOAD), but identification of the underlying causal genes remains challenging. Studies of induced pluripotent stem cell (iPSC)-derived neurons from LOAD patients have demonstrated the existence of neuronal cell-intrinsic functional defects. Here, we searched for genetic contributions to neuronal dysfunction in LOAD using an integrative systems approach that incorporated multi-evidence-based gene mapping and network-analysis-based prioritization. A systematic perturbation screening of candidate risk genes in Caenorhabditis elegans (C. elegans) revealed that neuronal knockdown of the LOAD risk gene orthologs vha-10 (ATP6V1G2), cmd-1 (CALM3), amph-1 (BIN1), ephx-1 (NGEF), and pho-5 (ACP2) alters short-/intermediate-term memory function, the cognitive domain affected earliest during LOAD progression. These results highlight the impact of LOAD risk genes on evolutionarily conserved memory function, as mediated through neuronal endosomal dysfunction, and identify new targets for further mechanistic interrogation.

Unraveling the significance of AGPAT4 for the pathogenesis of endometriosis via a multi-omics approach

Endometriosis is characterized by the ectopic proliferation of endometrial cells, posing considerable diagnostic and therapeutic challenges. Our study investigates AGPAT4's involvement in endometriosis pathogenesis, aiming to unveil new therapeutic targets. Our investigation by analyzing eQTL data from GWAS for preliminary screening. Subsequently, within the GEO dataset, we utilized four machine learning algorithms to precisely identify risk-associated genes. Gene validity was confirmed through five Mendelian Randomization methods. AGPAT4 expression was measured by Single-Cell Analysis, ELISA and immunohistochemistry. We investigated AGPAT4's effect on endometrial stromal cells using RNA interference, assessing cell proliferation, invasion, and migration with CCK8, wound-healing, and transwell assays. Protein expression was analyzed by western blot, and AGPAT4 interactions were explored using AutoDock. Our investigation identified 11 genes associated with endometriosis risk, with AGPAT4 and COMT emerging as pivotal biomarkers through machine learning analysis. AGPAT4 exhibited significant upregulation in both ectopic tissues and serum samples from patients with endometriosis. Reduced expression of AGPAT4 was observed to detrimentally impact the proliferation, invasion, and migration capabilities of endometrial stromal cells, concomitant with diminished expression of key signaling molecules such as Wnt3a, β-Catenin, MMP-9, and SNAI2. Molecular docking analyses further underscored a substantive interaction between AGPAT4 and Wnt3a.Our study highlights AGPAT4's key role in endometriosis, influencing endometrial stromal cell behavior, and identifies AGPAT4 pathways as promising therapeutic targets for this condition.

Genome wide association studies on seven yield-related traits of 183 rice varieties in Bangladesh

Rice genetic diversity is regulated by multiple genes and is largely dependent on various environmental factors. Uncovering the genetic variations associated with the diversity in rice populations is the key to breed stable and high yielding rice varieties. We performed genome wide association studies (GWASs) on seven rice yielding traits (grain length, grain width, grain weight, panicle length, leaf length, leaf width, and leaf angle) based on a population of 183 rice landraces of Bangladesh. Our GWASs reveal various chromosomal regions and candidate genes that are associated with different traits in Bangladeshi rice varieties. Noteworthy was the recurrent implication of chromosome 10 in all three grain-shape-related traits (grain length, grain width, and grain weight), indicating its pivotal role in shaping rice grain morphology. Our study also underscores the involvement of transposon gene families across these three traits. For leaf related traits, chromosome 10 was found to harbor regions that are significantly associated with leaf length and leaf width. The results of these association studies support previous findings as well as provide additional insights into the genetic diversity of rice. This is the first known GWAS study on various yield-related traits in the varieties of Oryza sativa available in Bangladesh-the fourth largest rice-producing country. We believe this study will accelerate rice genetics research and breeding stable high-yielding rice in Bangladesh.

Transcriptomic imputation of genetic risk variants uncovers novel whole-blood biomarkers of Parkinson's disease

Blood-based gene expression signatures could potentially be used as biomarkers for PD. However, it is unclear whether genetically-regulated transcriptomic signatures can provide novel gene candidates for use as PD biomarkers. We leveraged on the Genotype-Tissue Expression (GTEx) database to impute whole-blood transcriptomic expression using summary statistics of three large-scale PD GWAS. A random forest classifier was used with the consensus whole-blood imputed gene signature (IGS) to discriminate between cases and controls. Outcome measures included Area under the Curve (AUC) of Receiver Operating Characteristic (ROC) Curve. We demonstrated that the IGS (n = 37 genes) is conserved across PD GWAS studies and brain tissues. IGS discriminated between cases and controls in an independent whole-blood RNA-sequencing study (1176 PD, 254 prodromal, and 860 healthy controls) with mean AUC and accuracy of 64.8% and 69.4% for PD cohort, and 78.8% and 74% for prodromal cohort. PATL2 was the top-performing imputed gene in both PD and prodromal PD cohorts, whose classifier performance varied with biological sex (higher performance for males and females in the PD and prodromal PD, respectively). Single-cell RNA-sequencing studies (scRNA-seq) of healthy humans and PD patients found PATL2 to be enriched in terminal effector CD8+ and cytotoxic CD4+ cells, whose proportions are both increased in PD patients. We demonstrated the utility of GWAS transcriptomic imputation in identifying novel whole-blood transcriptomic signatures which could be leveraged upon for PD biomarker derivation. We identified PATL2 as a potential biomarker in both clinical and prodromic PD.

Beyond guilty by association at scale: searching for causal variants on the basis of genome-wide summary statistics

Understanding the causal genetic architecture of complex phenotypes is essential for future research into disease mechanisms and potential therapies. Here, we present a novel framework for genome-wide detection of sets of variants that carry non-redundant information on the phenotypes and are therefore more likely to be causal in a biological sense. Crucially, our framework requires only summary statistics obtained from standard genome-wide marginal association testing. The described approach, implemented in open-source software, is also computationally efficient, requiring less than 15 minutes on a single CPU to perform genome-wide analysis. Through extensive genome-wide simulation studies, we show that the method can substantially outperform usual two-stage marginal association testing and fine-mapping procedures in precision and recall. In applications to a meta-analysis of ten large-scale genetic studies of Alzheimer's disease (AD), we identified 82 loci associated with AD, including 37 additional loci missed by conventional GWAS pipeline. The identified putative causal variants achieve state-of-the-art agreement with massively parallel reporter assays and CRISPR-Cas9 experiments. Additionally, we applied the method to a retrospective analysis of 67 large-scale GWAS summary statistics since 2013 for a variety of phenotypes. Results reveal the method's capacity to robustly discover additional loci for polygenic traits and pinpoint potential causal variants underpinning each locus beyond conventional GWAS pipeline, contributing to a deeper understanding of complex genetic architectures in post-GWAS analyses.

Population Neuroscience: Principles and Advances

In population neuroscience, three disciplines come together to advance our knowledge of factors that shape the human brain: neuroscience, genetics, and epidemiology (Paus, Human Brain Mapping 31:891-903, 2010). Here, I will come back to some of the background material reviewed in more detail in our previous book (Paus, Population Neuroscience, 2013), followed by a brief overview of current advances and challenges faced by this integrative approach.

Depression and hepatobiliary diseases: a bidirectional Mendelian randomization study

Background

More and more evidence suggests a close association between depression and hepatobiliary diseases, but its causal relationship is not yet clear.

Method

Using genome-wide association studies (GWAS) to summarize data, independent genetic variations associated with depression were selected as instrumental variables. Firstly, we designed a univariate Mendelian randomization (UVMR) analysis with two samples and simultaneously conducted reverse validation to evaluate the potential bidirectional causal relationship between depression and various hepatobiliary diseases. Secondly, we conducted a multivariate Mendelian randomization (MVMR) analysis on diseases closely related to depression, exploring the mediating effects of waist to hip ratio, hypertension, and daytime nap. The mediating effects were obtained through MVMR. For UVMR and MVMR, inverse variance weighted method (IVW) is considered the most important analytical method. Sensitivity analysis was conducted using Cochran'Q, MR Egger, and Leave-one-out methods.

Results

UVMR analysis showed that depression may increase the risk of non-alcoholic fatty liver disease (OR, 1.22; 95% CI, 1.03-1.46; p=0.0248) in liver diseases, while depression does not increase the risk of other liver diseases; In biliary and pancreatic related diseases, depression may increase the risk of cholelithiasis (OR, 1.26; 95% CI, 1.05-1.50; p=0.0120), chronic pancreatitis (OR, 1.61; 95% CI, 1.10-2.35; p=0.0140), and cholecystitis (OR, 1.23; 95% CI, 1.03-1.48; p=0.0250). In addition, through reverse validation, we found that non-alcoholic fatty liver disease, cholelithiasis, chronic pancreatitis, cholecystitis, or the inability to increase the risk of depression (p>0.05). The waist to hip ratio, hypertension, and daytime nap play a certain role in the process of depression leading to non-alcoholic fatty liver disease, with a mediating effect of 35.8%.

Conclusion

Depression is a susceptibility factor for non-alcoholic fatty liver disease, and the causal effect of genetic susceptibility to depression on non-alcoholic fatty liver disease is mediated by waist-hip ratio, hypertension, and daytime nap.

Screen time in the development of cardiovascular diseases: A two-sample Mendelian randomization study

BACKGROUND AND AIMS

Coronary artery disease (CAD), heart failure (HF), and ischemic heart disease (IHD) are three common cardiovascular diseases that are closely associated with metabolic activity. The global incidence and prevalence of these conditions are on the rise, primarily due to unhealthy lifestyles, aging populations, and the increasing prevalence of obesity and diabetes. Excessive screen time has emerged as a potential risk factor for various adverse health outcomes, although limited research has explored its relationship with cardiovascular disease outcomes.

METHODS AND RESULTS

A Mendelian randomization (MR) study was conducted, employing exposure-associated genetic variants as instrumental variables to explore the causal relationship between screen time use and cardiovascular disease outcomes. Single nucleotide polymorphisms (SNPs) were utilized as pooled data for the genetic variable instrument, investigating the association between screen use duration and three types of cardiovascular diseases: coronary artery disease (CAD), heart failure (HF), and ischemic heart disease (IHD). Through the MR analysis, it was revealed that the use of mobile phones and TV screens exhibited a significant causal association with the occurrence of CAD, heart failure, and IHD. However, no significant association was observed between the use of computers and these three types of cardiovascular diseases.

CONCLUSION

Our study suggests that excessive screen time use is associated with the development of cardiovascular disease. However, it should be noted that the consequences of screen time can vary depending on the reasons and purposes for its use. Implementing reasonable control over screen time, particularly for entertainment purposes, holds promise as a potential approach to mitigating cardiovascular disease.

Identification of 27 allele-specific regulatory variants in Parkinson's disease using a massively parallel reporter assay

Genome wide association studies (GWAS) have identified a number of genomic loci that are associated with Parkinson's disease (PD) risk. However, the majority of these variants lie in non-coding regions, and thus the mechanisms by which they influence disease development, and/or potential subtypes, remain largely elusive. To address this, we used a massively parallel reporter assay (MPRA) to screen the regulatory function of 5254 variants that have a known or putative connection to PD. We identified 138 loci with enhancer activity, of which 27 exhibited allele-specific regulatory activity in HEK293 cells. The identified regulatory variant(s) typically did not match the original tag variant within the PD associated locus, supporting the need for deeper exploration of these loci. The existence of allele specific transcriptional impacts within HEK293 cells, confirms that at least a subset of the PD associated regions mark functional gene regulatory elements. Future functional studies that confirm the putative targets of the empirically verified regulatory variants will be crucial for gaining a greater understanding of how gene regulatory network(s) modulate PD risk.

WITHDRAWN: The effect of ischemia on expression quantitative trait loci (eQTL) in human myocardium and insights into myocardial injury etiology

26 February, 2024. Research Square has withdrawn this preprint as it was submitted and made public without the full consent of all the authors and without the full consent of the principle investigator of the registered clinical trial. Therefore, this work should not be cited as a reference.

Sex-stratified genome-wide association and transcriptome-wide Mendelian randomization studies reveal drug targets of heart failure

The prevalence of heart failure (HF) subtypes, which are classified by left ventricular ejection fraction (LVEF), demonstrate significant sex differences. Here, we perform sex-stratified genome-wide association studies (GWASs) on LVEF and transcriptome-wide Mendelian randomization (MR) on LVEF, all-cause HF, HF with reduced ejection fraction (HFrEF), and HF with preserved ejection fraction (HFpEF). The sex-stratified GWASs of LVEF identified three sex-specific loci that were exclusively detected in the sex-stratified GWASs. Three drug target genes show sex-differential effects on HF/HFrEF via influencing LVEF, with NPR2 as the target gene for the HF drug Cenderitide under phase 2 clinical trial. Our study highlights the importance of considering sex-differential genetic effects in sex-balanced diseases such as HF and emphasizes the value of sex-stratified GWASs and MR in identifying putative genetic variants, causal genes, and candidate drug targets for HF, which is not identifiable using a sex-combined strategy.

Investigating the association between genetically proxied circulating levels of immune checkpoint proteins and cancer survival: protocol for a Mendelian randomisation analysis

INTRODUCTION

Compared with the traditional drug development pathway, investigating alternative uses for existing drugs (ie, drug repurposing) requires substantially less time, cost and resources. Immune checkpoint inhibitors are licensed for the treatment of certain breast, colorectal, head and neck, lung and melanoma cancers. These drugs target immune checkpoint proteins to reduce the suppression of T cell activation by cancer cells. As T cell suppression is a hallmark of cancer common across anatomical sites, we hypothesise that immune checkpoint inhibitors could be repurposed for the treatment of additional cancers beyond the ones already indicated.

METHODS AND ANALYSIS

We will use two-sample Mendelian randomisation to investigate the effect of genetically proxied levels of protein targets of two immune checkpoint inhibitors-programmed cell death protein 1 and programmed death ligand 1-on survival of seven cancer types (breast, colorectal, head and neck, lung, melanoma, ovarian and prostate). Summary genetic association data will be obtained from prior genome-wide association studies of circulating protein levels and cancer survival in populations of European ancestry. Various sensitivity analyses will be performed to examine the robustness of findings to potential violations of Mendelian randomisation assumptions, collider bias and the impact of alternative genetic instrument construction strategies. The impact of treatment history and tumour stage on the findings will also be investigated using summary-level and individual-level genetic data where available.

ETHICS AND DISSEMINATION

No separate ethics approval will be required for these analyses as we will be using data from previously published genome-wide association studies which individually gained ethical approval and participant consent. Results from analyses will be submitted as an open-access peer-reviewed publication and statistical code will be made freely available on the completion of the analysis.

Adjusting for genetic confounders in transcriptome-wide association studies improves discovery of risk genes of complex traits

Many methods have been developed to leverage expression quantitative trait loci (eQTL) data to nominate candidate genes from genome-wide association studies. These methods, including colocalization, transcriptome-wide association studies (TWAS) and Mendelian randomization-based methods; however, all suffer from a key problem-when assessing the role of a gene in a trait using its eQTLs, nearby variants and genetic components of other genes' expression may be correlated with these eQTLs and have direct effects on the trait, acting as potential confounders. Our extensive simulations showed that existing methods fail to account for these 'genetic confounders', resulting in severe inflation of false positives. Our new method, causal-TWAS (cTWAS), borrows ideas from statistical fine-mapping and allows us to adjust all genetic confounders. cTWAS showed calibrated false discovery rates in simulations, and its application on several common traits discovered new candidate genes. In conclusion, cTWAS provides a robust statistical framework for gene discovery.

Exploring noncoding variants in genetic diseases: from detection to functional insights

Previous studies on genetic diseases predominantly focused on protein-coding variations, overlooking the vast noncoding regions in the human genome. The development of high-throughput sequencing technologies and functional genomics tools has enabled the systematic identification of functional noncoding variants. These variants can impact gene expression, regulation, and chromatin conformation, thereby contributing to disease pathogenesis. Understanding the mechanisms that underlie the impact of noncoding variants on genetic diseases is indispensable for the development of precisely targeted therapies and the implementation of personalized medicine strategies. The intricacies of noncoding regions introduce a multitude of challenges and research opportunities. In this review, we introduce a spectrum of noncoding variants involved in genetic diseases, along with research strategies and advanced technologies for their precise identification and in-depth understanding of the complexity of the noncoding genome. We will delve into the research challenges and propose potential solutions for unraveling the genetic basis of rare and complex diseases.

Utilizing multimodal approach to identify candidate pathways and biomarkers and predicting frailty syndrome in individuals from UK Biobank

Frailty, a prevalent clinical syndrome in aging adults, is characterized by poor health outcomes, represented via a standardized frailty-phenotype (FP), and Frailty Index (FI). While the relevance of the syndrome is gaining awareness, much remains unclear about its underlying biology. Further elucidation of the genetic determinants and possible underlying mechanisms may help improve patients' outcomes allowing healthy aging.Genotype, clinical and demographic data of subjects (aged 60-73 years) from UK Biobank were utilized. FP was defined on Fried's criteria. FI was calculated using electronic-health-records. Genome-wide-association-studies (GWAS) were conducted and polygenic-risk-scores (PRS) were calculated for both FP and FI. Functional analysis provided interpretations of underlying biology. Finally, machine-learning (ML) models were trained using clinical, demographic and PRS towards identifying frail from non-frail individuals.Thirty-one loci were significantly associated with FI accounting for 12% heritability. Seventeen of those were known associations for body-mass-index, coronary diseases, cholesterol-levels, and longevity, while the rest were novel. Significant genes CDKN2B and APOE, previously implicated in aging, were reported to be enriched in lipoprotein-particle-remodeling. Linkage-disequilibrium-regression identified specific regulation in limbic-system, associated with long-term memory and cognitive-function. XGboost was established as the best performing ML model with area-under-curve as 85%, sensitivity and specificity as 0.75 and 0.8, respectively.This study provides novel insights into increased vulnerability and risk stratification of frailty syndrome via a multi-modal approach. The findings suggest frailty as a highly polygenic-trait, enriched in cholesterol-remodeling and metabolism and to be genetically associated with cognitive abilities. ML models utilizing FP and FI + PRS were established that identified frailty-syndrome patients with high accuracy.

Cerebrospinal and Brain Proteins Implicated in Neuropsychiatric and Risk Factor Traits: Evidence from Mendelian Randomization

Neuropsychiatric disorders present a global health challenge, necessitating an understanding of their molecular mechanisms for therapeutic development. Using Mendelian randomization (MR) analysis, this study explored associations between genetically predicted levels of 173 proteins in cerebrospinal fluid (CSF) and 25 in the brain with 14 neuropsychiatric disorders and risk factors. Follow-up analyses assessed consistency across plasma protein levels and gene expression in various brain regions. Proteins were instrumented using tissue-specific genetic variants, and colocalization analysis confirmed unbiased gene variants. Consistent MR and colocalization evidence revealed that lower cortical expression of low-density lipoprotein receptor-related protein 8, coupled higher abundance in the CSF and plasma, associated with lower fluid intelligence scores and decreased bipolar disorder risk. Additionally, elevated apolipoprotein-E2 and hepatocyte growth factor-like protein in the CSF and brain were related to reduced leisure screen time and lower odds of physical activity, respectively. Furthermore, elevated CSF soluble tyrosine-protein kinase receptor 1 level increased liability to attention deficit hyperactivity disorder and schizophrenia alongside lower fluid intelligence scores. This research provides genetic evidence supporting novel tissue-specific proteomic targets for neuropsychiatric disorders and their risk factors. Further exploration is necessary to understand the underlying biological mechanisms and assess their potential for therapeutic intervention.

Metabolic syndrome and risk of colorectal cancer: A Mendelian randomization study

Background

Observational studies have previously demonstrated a significant relationship among both metabolic syndrome (Mets) and colorectal cancer (CRC). Whether there is a causal link remains controversial.

Objective

To clarify whether Mets and their components have a causal effect on colorectal cancer, we have carried out a bidirectional Mendelian randomization analysis (MR).

Methods

This study started from genome-wide association data for Mets and its 5 components (hypertension, waist circumference, fasting blood glucose, serum triglycerides, and serum high-density lipoprotein cholesterol) and colorectal cancer. Mendelian randomization (MR) techniques were used in the study to examine their associations.

Results

After Benjamini-Hochberg multiple corrections, genetically predicted significant causal link exists between WC (waist circumference) and CRC. The OR was 1.35 (95 % CI: 1.08-1.69; p = 0.0096). Other Mets components (HBP, FBG, TG, HDL), on the other hand, found no evidence of a genetic link between CRC and Mets. In addition, MR results showed that CRC was not causally related to either Mets or the components. We get the same result in the validated dataset.

Conclusion

According to the bidirectional MR investigation shows a significant causal relationship among obesity and CRC in the Mets component but no causal relationship in the opposite direction.

An inverse causal relationship between serum 25-hydroxyvitamin D levels and pulmonary hypertension: A two-sample Mendelian randomization study

Observational studies have confirmed that 25-hydroxyvitamin D (25(OH)D) is associated with pulmonary hypertension (PH), but the causal association between each other is unclear. Therefore, Mendelian randomization (MR) method was performed to validate the causal association between PH and serum 25(OH)D levels. The summary data for 25(OH)D and PH were from the National Human Genome Research Institute-European Bioinformatics Institute. Catalog of human genome-wide association studies and FinnGen biobank consortium. MR analysis was utilized to explore the potential causal association between PH and 25(OH)D. To evaluate this association, inverse variance weighting was considered as the primary method. Cochran's Q test, MR-Egger intercept test, and "leave-one-out" sensitivity analyses were utilized to control the pleiotropy and heterogeneity in the study. Two-sample MR analysis revealed an inverse causal relationship between 25(OH)D and PH (odds ratio: 0.376, 95% confidence interval: 0.162-0.876, p = 2.334 × 10-2). There was no significant heterogeneity and pleiotropy. The present study confirmed the inverse causal relationship between 25(OH)D and PH. This pathway may provide another treatment pathway in PH. Further studies to elucidate this pathway is indicated.

Investigating the causal impact of gut microbiota on glioblastoma: a bidirectional Mendelian randomization study

BACKGROUND

Currently, the influence of microbiota on the occurrence, progression, and treatment of cancer is a topic of considerable research interest. Therefore, based on the theory of the gut-brain axis proved by previous studies, our objective was to uncover the causal relationship between glioblastoma and the gut microbiome using Mendelian randomization analysis.

METHODS

We conducted a bidirectional Mendelian randomization study using summary statistics of gut microbiota derived from the MiBioGen consortium, the largest database of gut microbiota. Summary statistics for glioblastoma were obtained from IEU OpenGWAS project, which included 91 cases and 218,701 controls. We assessed the presence of heterogeneity and horizontal pleiotropy in the analyzed data. We primarily employed the inverse variance weighting method to investigate the causal relationship between gut microbiota and glioblastoma after excluding cases of horizontal pleiotropy. Four other analysis methods were employed as supplementary. Excluding abnormal results based on leave-one-out sensitivity analysis. Finally, reverse Mendelian randomization analysis was performed.

RESULTS

Four genus-level taxa and one family-level taxa exhibited causal associations with glioblastoma. And these results of reverse Mendelian randomization analysis shown glioblastoma exhibited causal associations with three genus-level taxa and one family-level taxa. However, the Prevotella7(Forward, P=0.006, OR=0.34, 95%CI:0.158-0.732; Reverse, P=0.004, OR=0.972, 95%CI:0.953-0.991) shown the causal associations with glioblastoma in the bidirectional Mendelian randomization.

CONCLUSIONS

In this bidirectional Mendelian randomization study, we identified five gut microbiota species with causal associations to glioblastoma. However, additional randomized controlled trials are required to clarify the impact of gut microbiota on glioblastoma and to reveal its precise mechanisms.

Gene signatures derived from transcriptomic-causal networks stratified colorectal cancer patients for effective targeted therapy

Predictive and prognostic gene signatures derived from interconnectivity among genes can tailor clinical care to patients in cancer treatment. We identified gene interconnectivity as the transcriptomic-causal network by integrating germline genotyping and tumor RNA-seq data from 1,165 patients with metastatic colorectal cancer (CRC). The patients were enrolled in a clinical trial with randomized treatment, either cetuximab or bevacizumab in combination with chemotherapy. We linked the network to overall survival (OS) and detected novel biomarkers by controlling for confounding genes. Our data-driven approach discerned sets of genes, each set collectively stratify patients based on OS. Two signatures under the cetuximab treatment were related to wound healing and macrophages. The signature under the bevacizumab treatment was related to cytotoxicity and we replicated its effect on OS using an external cohort. We also showed that the genes influencing OS within the signatures are downregulated in CRC tumor vs. normal tissue using another external cohort. Furthermore, the corresponding proteins encoded by the genes within the signatures interact each other and are functionally related. In conclusion, this study identified a group of genes that collectively stratified patients based on OS and uncovered promising novel prognostic biomarkers for personalized treatment of CRC using transcriptomic causal networks.

Genetic insights into the association of statin and newer nonstatin drug target genes with human longevity: a Mendelian randomization analysis

BACKGROUND

It remains controversial whether the long-term use of statins or newer nonstatin drugs has a positive effect on human longevity. Therefore, this study aimed to investigate the genetic associations between different lipid-lowering therapeutic gene targets and human longevity.

METHODS

Two-sample Mendelian randomization analyses were conducted. The exposures comprised genetic variants that proxy nine drug target genes mimicking lipid-lowering effects (LDLR, HMGCR, PCKS9, NPC1L1, APOB, CETP, LPL, APOC3, and ANGPTL3). Two large-scale genome-wide association study (GWAS) summary datasets of human lifespan, including up to 500,193 European individuals, were used as outcomes. The inverse-variance weighting method was applied as the main approach. Sensitivity tests were conducted to evaluate the robustness, heterogeneity, and pleiotropy of the results. Causal effects were further validated using expression quantitative trait locus (eQTL) data.

RESULTS

Genetically proxied LDLR variants, which mimic the effects of lowering low-density lipoprotein cholesterol (LDL-C), were associated with extended lifespan. This association was replicated in the validation set and was further confirmed in the eQTL summary data of blood and liver tissues. Mediation analysis revealed that the genetic mimicry of LDLR enhancement extended lifespan by reducing the risk of major coronary heart disease, accounting for 22.8% of the mediation effect. The genetically proxied CETP and APOC3 inhibitions also showed causal effects on increased life expectancy in both outcome datasets. The lipid-lowering variants of HMGCR, PCKS9, LPL, and APOB were associated with longer lifespans but did not causally increase extreme longevity. No statistical evidence was detected to support an association between NPC1L1 and lifespan.

CONCLUSION

This study suggests that LDLR is a promising genetic target for human longevity. Lipid-related gene targets, such as PCSK9, CETP, and APOC3, might potentially regulate human lifespan, thus offering promising prospects for developing newer nonstatin therapies.

An augmented Mendelian randomization approach provides causality of brain imaging features on complex traits in a single biobank-scale dataset

Mendelian randomization (MR) is an effective approach for revealing causal risk factors that underpin complex traits and diseases. While MR has been more widely applied under two-sample settings, it is more promising to be used in one single large cohort given the rise of biobank-scale datasets that simultaneously contain genotype data, brain imaging data, and matched complex traits from the same individual. However, most existing multivariable MR methods have been developed for two-sample setting or a small number of exposures. In this study, we introduce a one-sample multivariable MR method based on partial least squares and Lasso regression (MR-PL). MR-PL is capable of considering the correlation among exposures (e.g., brain imaging features) when the number of exposures is extremely upscaled, while also correcting for winner's curse bias. We performed extensive and systematic simulations, and demonstrated the robustness and reliability of our method. Comprehensive simulations confirmed that MR-PL can generate more precise causal estimates with lower false positive rates than alternative approaches. Finally, we applied MR-PL to the datasets from UK Biobank to reveal the causal effects of 36 white matter tracts on 180 complex traits, and showed putative white matter tracts that are implicated in smoking, blood vascular function-related traits, and eating behaviors.

Proteome-wide mendelian randomization identifies causal plasma proteins in venous thromboembolism development

Genome-wide association studies (GWAS) have identified numerous risk loci for venous thromboembolism (VTE), but it is challenging to decipher the underlying mechanisms. We employed an integrative analytical pipeline to transform genetic associations to identify novel plasma proteins for VTE. Proteome-wide association studies (PWAS) were determined by functional summary-based imputation leveraging data from a genome-wide association analysis (14,429 VTE patients, 267,037 controls), blood proteomes (1348 cases), followed by Mendelian randomization, Bayesian colocalization, protein-protein interaction, and pathway enrichment analysis. Twenty genetically regulated circulating protein abundances (F2, F11, ABO, PLCG2, LRP4, PLEK, KLKB1, PROC, KNG1, THBS2, SERPINA1, RARRES2, CEL, GP6, SERPINE2, SERPINA10, OBP2B, EFEMP1, F5, and MSR1) were associated with VTE. Of these 13 proteins demonstrated Mendelian randomized correlations. Six proteins (F2, F11, PLEK, SERPINA1, RARRES2, and SERPINE2) had strong support in colocalization analysis. Utilizing multidimensional data, this study suggests PLEK, SERPINA1, and SERPINE2 as compelling proteins that may provide key hints for future research and possible diagnostic and therapeutic targets for VTE.

Leveraging molecular quantitative trait loci to comprehend complex diseases/traits from the omics perspective

Comprehending the molecular basis of quantitative genetic variation is a principal goal for complex diseases or traits. Molecular quantitative trait loci (molQTLs) have made it possible to investigate the effects of genetic variants hiding behind large-scale omics data. A deeper understanding of molQTL is urgently required in light of the multi-dimensionalization of omics data to more fully elucidate the pertinent biological mechanisms. Herein, we reviewed molQTLs with the corresponding resource from the omics perspective and further discussed the integrative strategy of GWAS-molQTL to infer their causal effects. Subsequently, we described the opportunities and challenges encountered by molQTL. The case studies showed that molQTL is essential for complex diseases and traits, whether single- or multi-omics QTLs. Overall, we highlighted the functional significance of genetic variants to employ the discovery of molQTL in complex diseases and traits.

Summary-data based Mendelian randomization identifies gene expression regulatory polymorphisms associated with bovine paratuberculosis by modulation of the nuclear factor Kappa β (NF-κß)-mediated inflammatory response

Genome-wide association studies (GWAS) have identified host genetic variants associated with paratuberculosis (PTB) susceptibility. Most of the GWAS-identified SNPs are in non-coding regions. Connecting these non-coding variants and downstream affected genes is a challenge and, up to date, only a few functional mutations or expression quantitative loci (cis-eQTLs) associated with PTB susceptibility have been identified. In the current study, the associations between imputed whole-genome sequence genotypes and whole RNA-Sequencing data from peripheral blood (PB) and ileocecal valve (ICV) samples of Spanish Holstein cows (N = 16) were analyzed with TensorQTL. This approach allowed the identification of 88 and 37 cis-eQTLs regulating the expression levels of 90 and 37 genes in PB and ICV samples, respectively (False discorey rate, FDR ≤ 0.05). Next, we applied summary-based data Mendelian randomization (SMR) to integrate the cis-eQTL dataset with GWAS data obtained from a cohort of 813 culled cattle that were classified according to the presence or absence of PTB-associated histopathological lesions in gut tissues. After multiple testing corrections (FDR ≤ 0.05), we identified two novel cis-eQTLs affecting the expression of the early growth response factor 4 (EGR4) and the bovine neuroblastoma breakpoint family member 6-like protein isoform 2 (MGC134040) that showed pleiotropic associations with the presence of multifocal and diffuse lesions in gut tissues; P = 0.002 and P = 0.017, respectively. While EGR4 acts as a brake on T-cell proliferation and cytokine production through interaction with the nuclear factor Kappa β (NF-κß), MGC134040 is a target gene of NF-κß. Our findings provide a better understanding of the genetic factors influencing PTB outcomes, confirm that the multifocal lesions are localized/confined lesions that have different underlying host genetics than the diffuse lesions, and highlight regulatory SNPs and regulated-gene targets to design future functional studies.

Molecular bases of comorbidities: present and future perspectives

Co-occurrence of diseases decreases patient quality of life, complicates treatment choices, and increases mortality. Analyses of electronic health records present a complex scenario of comorbidity relationships that vary by age, sex, and cohort under study. The study of similarities between diseases using 'omics data, such as genes altered in diseases, gene expression, proteome, and microbiome, are fundamental to uncovering the origin of, and potential treatment for, comorbidities. Recent studies have produced a first generation of genetic interpretations for as much as 46% of the comorbidities described in large cohorts. Integrating different sources of molecular information and using artificial intelligence (AI) methods are promising approaches for the study of comorbidities. They may help to improve the treatment of comorbidities, including the potential repositioning of drugs.

Association between frailty and depression: A bidirectional Mendelian randomization study

Frailty and depression were linked in observational studies, but the causality remains ambiguous. We intended to explore it using Mendelian randomization (MR). We obtained frailty genome-wide association study (GWAS) data from UK Biobank and TwinGen meta-analysis, and depression GWAS data from Psychiatric Genomics Consortium (PGC) and FinnGen (respectively recorded as PD and FD). We performed univariable and multivariable-adjusted MR with adjustments for body mass index (BMI) and physical activity (PA). Frailty was significantly associated with elevated risks of PD (OR, 1.860; 95% CI, 1.439 to 2.405; P < 0.001) and FD (OR, 1.745; 95% CI, 1.193 to 2.552; P = 0.004), and depression was meanwhile a susceptible factor for frailty (PD: β, 0.146; 95% CI, 0.086 to 0.201; P < 0.001; and FD: β, 0.112; 95% CI, 0.051 to 0.174; P < 0.001). This association was robust after adjustments for BMI or PA. Our study provides evidence of the bidirectional causal association between frailty and depression from genetic perspectives.

On the interpretation of transcriptome-wide association studies

Transcriptome-wide association studies (TWAS) aim to detect relationships between gene expression and a phenotype, and are commonly used for secondary analysis of genome-wide association study (GWAS) results. Results from TWAS analyses are often interpreted as indicating a genetic relationship between gene expression and a phenotype, but this interpretation is not consistent with the null hypothesis that is evaluated in the traditional TWAS framework. In this study we provide a mathematical outline of this TWAS framework, and elucidate what interpretations are warranted given the null hypothesis it actually tests. We then use both simulations and real data analysis to assess the implications of misinterpreting TWAS results as indicative of a genetic relationship between gene expression and the phenotype. Our simulation results show considerably inflated type 1 error rates for TWAS when interpreted this way, with 41% of significant TWAS associations detected in the real data analysis found to have insufficient statistical evidence to infer such a relationship. This demonstrates that in current implementations, TWAS cannot reliably be used to investigate genetic relationships between gene expression and a phenotype, but that local genetic correlation analysis can serve as a potential alternative.

Allelic variation in transcription factor PtoWRKY68 contributes to drought tolerance in Populus

Drought stress limits woody species productivity and influences tree distribution. However, dissecting the molecular mechanisms that underpin drought responses in forest trees can be challenging due to trait complexity. Here, using a panel of 300 Chinese white poplar (Populus tomentosa) accessions collected from different geographical climatic regions in China, we performed a genome-wide association study (GWAS) on seven drought-related traits and identified PtoWRKY68 as a candidate gene involved in the response to drought stress. A 12-bp insertion and/or deletion and three nonsynonymous variants in the PtoWRKY68 coding sequence categorized natural populations of P. tomentosa into two haplotype groups, PtoWRKY68hap1 and PtoWRKY68hap2. The allelic variation in these two PtoWRKY68 haplotypes conferred differential transcriptional regulatory activities and binding to the promoters of downstream abscisic acid (ABA) efflux and signaling genes. Overexpression of PtoWRKY68hap1 and PtoWRKY68hap2 in Arabidopsis (Arabidopsis thaliana) ameliorated the drought tolerance of two transgenic lines and increased ABA content by 42.7% and 14.3% compared to wild-type plants, respectively. Notably, PtoWRKY68hap1 (associated with drought tolerance) is ubiquitous in accessions in water-deficient environments, whereas the drought-sensitive allele PtoWRKY68hap2 is widely distributed in well-watered regions, consistent with the trends in local precipitation, suggesting that these alleles correspond to geographical adaptation in Populus. Moreover, quantitative trait loci analysis and an electrophoretic mobility shift assay showed that SHORT VEGETATIVE PHASE (PtoSVP.3) positively regulates the expression of PtoWRKY68 under drought stress. We propose a drought tolerance regulatory module in which PtoWRKY68 modulates ABA signaling and accumulation, providing insight into the genetic basis of drought tolerance in trees. Our findings will facilitate molecular breeding to improve the drought tolerance of forest trees.

Using Drosophila to identify naturally occurring genetic modifiers of amyloid beta 42- and tau-induced toxicity

Alzheimer's disease is characterized by 2 pathological proteins, amyloid beta 42 and tau. The majority of Alzheimer's disease cases in the population are sporadic and late-onset Alzheimer's disease, which exhibits high levels of heritability. While several genetic risk factors for late-onset Alzheimer's disease have been identified and replicated in independent studies, including the ApoE ε4 allele, the great majority of the heritability of late-onset Alzheimer's disease remains unexplained, likely due to the aggregate effects of a very large number of genes with small effect size, as well as to biases in sample collection and statistical approaches. Here, we present an unbiased forward genetic screen in Drosophila looking for naturally occurring modifiers of amyloid beta 42- and tau-induced ommatidial degeneration. Our results identify 14 significant SNPs, which map to 12 potential genes in 8 unique genomic regions. Our hits that are significant after genome-wide correction identify genes involved in neuronal development, signal transduction, and organismal development. Looking more broadly at suggestive hits (P < 10-5), we see significant enrichment in genes associated with neurogenesis, development, and growth as well as significant enrichment in genes whose orthologs have been identified as significantly or suggestively associated with Alzheimer's disease in human GWAS studies. These latter genes include ones whose orthologs are in close proximity to regions in the human genome that are associated with Alzheimer's disease, but where a causal gene has not been identified. Together, our results illustrate the potential for complementary and convergent evidence provided through multitrait GWAS in Drosophila to supplement and inform human studies, helping to identify the remaining heritability and novel modifiers of complex diseases.

Supervised structural learning of semiparametric regression on high-dimensional correlated covariates with applications to eQTL studies

Expression quantitative trait loci (eQTL) studies utilize regression models to explain the variance of gene expressions with genetic loci or single nucleotide polymorphisms (SNPs). However, regression models for eQTL are challenged by the presence of high dimensional non-sparse and correlated SNPs with small effects, and nonlinear relationships between responses and SNPs. Principal component analyses are commonly conducted for dimension reduction without considering responses. Because of that, this non-supervised learning method often does not work well when the focus is on discovery of the response-covariate relationship. We propose a new supervised structural dimensional reduction method for semiparametric regression models with high dimensional and correlated covariates; we extract low-dimensional latent features from a vast number of correlated SNPs while accounting for their relationships, possibly nonlinear, with gene expressions. Our model identifies important SNPs associated with gene expressions and estimates the association parameters via a likelihood-based algorithm. A GTEx data application on a cancer related gene is presented with 18 novel eQTLs detected by our method. In addition, extensive simulations show that our method outperforms the other competing methods in bias, efficiency, and computational cost.

Chromosomal deletions on 16p11.2 encompassing SH2B1 are associated with accelerated metabolic disease

New approaches are needed to treat people whose obesity and type 2 diabetes (T2D) are driven by specific mechanisms. We investigate a deletion on chromosome 16p11.2 (breakpoint 2-3 [BP2-3]) encompassing SH2B1, a mediator of leptin and insulin signaling. Phenome-wide association scans in the UK (N = 502,399) and Estonian (N = 208,360) biobanks show that deletion carriers have increased body mass index (BMI; p = 1.3 × 10-10) and increased rates of T2D. Compared with BMI-matched controls, deletion carriers have an earlier onset of T2D, with poorer glycemic control despite higher medication usage. Cystatin C, a biomarker of kidney function, is significantly elevated in deletion carriers, suggesting increased risk of renal impairment. In a Mendelian randomization study, decreased SH2B1 expression increases T2D risk (p = 8.1 × 10-6). We conclude that people with 16p11.2 BP2-3 deletions have early, complex obesity and T2D and may benefit from therapies that enhance leptin and insulin signaling.

Guidelines for performing Mendelian randomization investigations: update for summer 2023

This paper provides guidelines for performing Mendelian randomization investigations. It is aimed at practitioners seeking to undertake analyses and write up their findings, and at journal editors and reviewers seeking to assess Mendelian randomization manuscripts. The guidelines are divided into ten sections: motivation and scope, data sources, choice of genetic variants, variant harmonization, primary analysis, supplementary and sensitivity analyses (one section on robust statistical methods and one on other approaches), extensions and additional analyses, data presentation, and interpretation. These guidelines will be updated based on feedback from the community and advances in the field. Updates will be made periodically as needed, and at least every 24 months.

Transcriptomic risk scores for attention deficit/hyperactivity disorder

Attention deficit/hyperactivity disorder (ADHD) is a highly heritable neurodevelopmental disorder. We performed a transcriptome-wide association study (TWAS) using the latest genome-wide association study (GWAS) meta-analysis, in 38,691 individuals with ADHD and 186,843 controls, and 14 gene-expression reference panels across multiple brain tissues and whole blood. Based on TWAS results, we selected subsets of genes and constructed transcriptomic risk scores (TRSs) for the disorder in peripheral blood mononuclear cells of individuals with ADHD and controls. We found evidence of association between ADHD and TRSs constructed using expression profiles from multiple brain areas, with individuals with ADHD carrying a higher burden of TRSs than controls. TRSs were uncorrelated with the polygenic risk score (PRS) for ADHD and, in combination with PRS, improved significantly the proportion of variance explained over the PRS-only model. These results support the complementary predictive potential of genetic and transcriptomic profiles in blood and underscore the potential utility of gene expression for risk prediction and deeper insight in molecular mechanisms underlying ADHD.