Phenome-wide Mendelian randomization mapping the influence of the plasma proteome on complex diseases

Stroke genetics informs drug discovery and risk prediction across ancestries

Previous genome-wide association studies (GWASs) of stroke - the second leading cause of death worldwide - were conducted predominantly in populations of European ancestry1,2. Here, in cross-ancestry GWAS meta-analyses of 110,182 patients who have had a stroke (five ancestries, 33% non-European) and 1,503,898 control individuals, we identify association signals for stroke and its subtypes at 89 (61 new) independent loci: 60 in primary inverse-variance-weighted analyses and 29 in secondary meta-regression and multitrait analyses. On the basis of internal cross-ancestry validation and an independent follow-up in 89,084 additional cases of stroke (30% non-European) and 1,013,843 control individuals, 87% of the primary stroke risk loci and 60% of the secondary stroke risk loci were replicated (P < 0.05). Effect sizes were highly correlated across ancestries. Cross-ancestry fine-mapping, in silico mutagenesis analysis3, and transcriptome-wide and proteome-wide association analyses revealed putative causal genes (such as SH3PXD2A and FURIN) and variants (such as at GRK5 and NOS3). Using a three-pronged approach4, we provide genetic evidence for putative drug effects, highlighting F11, KLKB1, PROC, GP1BA, LAMC2 and VCAM1 as possible targets, with drugs already under investigation for stroke for F11 and PROC. A polygenic score integrating cross-ancestry and ancestry-specific stroke GWASs with vascular-risk factor GWASs (integrative polygenic scores) strongly predicted ischaemic stroke in populations of European, East Asian and African ancestry5. Stroke genetic risk scores were predictive of ischaemic stroke independent of clinical risk factors in 52,600 clinical-trial participants with cardiometabolic disease. Our results provide insights to inform biology, reveal potential drug targets and derive genetic risk prediction tools across ancestries.

Synthesis of Surface-Functionalized Molybdenum Disulfide Nanomaterials for Efficient Adsorption and Deep Profiling of the Human Plasma Proteome by Data-Independent Acquisition

Blood is one of the most important clinical samples for protein biomarker discovery, as it provides rich physiological and pathological information and is easy to obtain with low invasiveness. However, the discovery of protein biomarkers in the blood by mass spectrometry (MS)-based proteomic strategies has been shown to be highly challenging due to the particularly large concentration range of proteins and the strong interference by the high-abundant proteins in the blood. Therefore, developing sensitive methods for low-abundant biomarker protein identification is a key issue that has received great attention. Here, we report the synthesis and characterization of surface-functionalized magnetic molybdenum disulfide (MoS₂) for the large-scale adsorption of low-abundant plasma proteins and deep profiling by MS. MoS₂ nanomaterials resulted in the coverage of more than 3400 proteins (including a single-peptide hit) in a single LC-MS analysis without peptide prefractionation using pooled plasma samples, which were five times more than those obtained by the direct analysis of the plasma proteome. A detection limit in the low ng L^-1 range was obtained, which is rare compared with previous reports.

Systematic Mendelian randomization using the human plasma proteome to discover potential therapeutic targets for stroke

Stroke is the second leading cause of death with substantial unmet therapeutic needs. To identify potential stroke therapeutic targets, we estimate the causal effects of 308 plasma proteins on stroke outcomes in a two-sample Mendelian randomization framework and assess mediation effects by stroke risk factors. We find associations between genetically predicted plasma levels of six proteins and stroke (P ≤ 1.62 × 10-4). The genetic associations with stroke colocalize (Posterior Probability >0.7) with the genetic associations of four proteins (TFPI, TMPRSS5, CD6, CD40). Mendelian randomization supports atrial fibrillation, body mass index, smoking, blood pressure, white matter hyperintensities and type 2 diabetes as stroke risk factors (P ≤ 0.0071). Body mass index, white matter hyperintensity and atrial fibrillation appear to mediate the TFPI, IL6RA, TMPRSS5 associations with stroke. Furthermore, thirty-six proteins are associated with one or more of these risk factors using Mendelian randomization. Our results highlight causal pathways and potential therapeutic targets for stroke.

Global Biobank Meta-analysis Initiative: Powering genetic discovery across human disease

Biobanks facilitate genome-wide association studies (GWASs), which have mapped genomic loci across a range of human diseases and traits. However, most biobanks are primarily composed of individuals of European ancestry. We introduce the Global Biobank Meta-analysis Initiative (GBMI)-a collaborative network of 23 biobanks from 4 continents representing more than 2.2 million consented individuals with genetic data linked to electronic health records. GBMI meta-analyzes summary statistics from GWASs generated using harmonized genotypes and phenotypes from member biobanks for 14 exemplar diseases and endpoints. This strategy validates that GWASs conducted in diverse biobanks can be integrated despite heterogeneity in case definitions, recruitment strategies, and baseline characteristics. This collaborative effort improves GWAS power for diseases, benefits understudied diseases, and improves risk prediction while also enabling the nomination of disease genes and drug candidates by incorporating gene and protein expression data and providing insight into the underlying biology of human diseases and traits.

A practical guideline of genomics-driven drug discovery in the era of global biobank meta-analysis

Genomics-driven drug discovery is indispensable for accelerating the development of novel therapeutic targets. However, the drug discovery framework based on evidence from genome-wide association studies (GWASs) has not been established, especially for cross-population GWAS meta-analysis. Here, we introduce a practical guideline for genomics-driven drug discovery for cross-population meta-analysis, as lessons from the Global Biobank Meta-analysis Initiative (GBMI). Our drug discovery framework encompassed three methodologies and was applied to the 13 common diseases targeted by GBMI (N _mean = 1,329,242). Individual methodologies complementarily prioritized drugs and drug targets, which were systematically validated by referring previously known drug-disease relationships. Integration of the three methodologies provided a comprehensive catalog of candidate drugs for repositioning, nominating promising drug candidates targeting the genes involved in the coagulation process for venous thromboembolism and the interleukin-4 and interleukin-13 signaling pathway for gout. Our study highlighted key factors for successful genomics-driven drug discovery using cross-population meta-analyses.

Phosphoproteomics of three exercise modalities identifies canonical signaling and C18ORF25 as an AMPK substrate regulating skeletal muscle function

Exercise induces signaling networks to improve muscle function and confer health benefits. To identify divergent and common signaling networks during and after different exercise modalities, we performed a phosphoproteomic analysis of human skeletal muscle from a cross-over intervention of endurance, sprint, and resistance exercise. This identified 5,486 phosphosites regulated during or after at least one type of exercise modality and only 420 core phosphosites common to all exercise. One of these core phosphosites was S67 on the uncharacterized protein C18ORF25, which we validated as an AMPK substrate. Mice lacking C18ORF25 have reduced skeletal muscle fiber size, exercise capacity, and muscle contractile function, and this was associated with reduced phosphorylation of contractile and Ca2+ handling proteins. Expression of C18ORF25 S66/67D phospho-mimetic reversed the decreased muscle force production. This work defines the divergent and canonical exercise phosphoproteome across different modalities and identifies C18ORF25 as a regulator of exercise signaling and muscle function.

Evaluating the efficacy and mechanism of metformin targets on reducing Alzheimer's disease risk in the general population: a Mendelian randomisation study

AIMS/HYPOTHESIS

Metformin use has been associated with reduced incidence of dementia in diabetic individuals in observational studies. However, the causality between the two in the general population is unclear. This study uses Mendelian randomisation (MR) to investigate the causal effect of metformin targets on Alzheimer's disease and potential causal mechanisms in the brain linking the two.

METHODS

Genetic proxies for the effects of metformin drug targets were identified as variants in the gene for the corresponding target that associated with HbA1c level (N=344,182) and expression level of the corresponding gene (N≤31,684). The cognitive outcomes were derived from genome-wide association studies comprising 527,138 middle-aged Europeans, including 71,880 with Alzheimer's disease or Alzheimer's disease-by-proxy. MR estimates representing lifelong metformin use on Alzheimer's disease and cognitive function in the general population were generated. Effect of expression level of 22 metformin-related genes in brain cortex (N=6601 donors) on Alzheimer's disease was further estimated.

RESULTS

Genetically proxied metformin use, equivalent to a 6.75 mmol/mol (1.09%) reduction on HbA1c, was associated with 4% lower odds of Alzheimer's disease (OR 0.96 [95% CI 0.95, 0.98], p=1.06×10-4) in non-diabetic individuals. One metformin target, mitochondrial complex 1 (MCI), showed a robust effect on Alzheimer's disease (OR 0.88, p=4.73×10-4) that was independent of AMP-activated protein kinase. MR of expression in brain cortex tissue showed that decreased MCI-related gene (NDUFA2) expression was associated with lower Alzheimer's disease risk (OR 0.95, p=4.64×10-4) and favourable cognitive function.

CONCLUSIONS/INTERPRETATION

Metformin use may cause reduced Alzheimer's disease risk in the general population. Mitochondrial function and the NDUFA2 gene are plausible mechanisms of action in dementia protection.

Clinlabomics: leveraging clinical laboratory data by data mining strategies

The recent global focus on big data in medicine has been associated with the rise of artificial intelligence (AI) in diagnosis and decision-making following recent advances in computer technology. Up to now, AI has been applied to various aspects of medicine, including disease diagnosis, surveillance, treatment, predicting future risk, targeted interventions and understanding of the disease. There have been plenty of successful examples in medicine of using big data, such as radiology and pathology, ophthalmology cardiology and surgery. Combining medicine and AI has become a powerful tool to change health care, and even to change the nature of disease screening in clinical diagnosis. As all we know, clinical laboratories produce large amounts of testing data every day and the clinical laboratory data combined with AI may establish a new diagnosis and treatment has attracted wide attention. At present, a new concept of radiomics has been created for imaging data combined with AI, but a new definition of clinical laboratory data combined with AI has lacked so that many studies in this field cannot be accurately classified. Therefore, we propose a new concept of clinical laboratory omics (Clinlabomics) by combining clinical laboratory medicine and AI. Clinlabomics can use high-throughput methods to extract large amounts of feature data from blood, body fluids, secretions, excreta, and cast clinical laboratory test data. Then using the data statistics, machine learning, and other methods to read more undiscovered information. In this review, we have summarized the application of clinical laboratory data combined with AI in medical fields. Undeniable, the application of Clinlabomics is a method that can assist many fields of medicine but still requires further validation in a multi-center environment and laboratory.

Genetics-informed precision treatment formulation in schizophrenia and bipolar disorder

Genetically informed drug development and repurposing is an attractive prospect for improving patient outcomes in psychiatry; however, the effectiveness of these endeavors is confounded by heterogeneity. We propose an approach that links interventions implicated by disorder-associated genetic risk, at the population level, to a framework that can target these compounds to individuals. Specifically, results from genome-wide association studies are integrated with expression data to prioritize individual "directional anchor" genes for which the predicted risk-increasing direction of expression could be counteracted by an existing drug. While these compounds represent plausible therapeutic candidates, they are not likely to be equally efficacious for all individuals. To account for this heterogeneity, we constructed polygenic scores restricted to variants annotated to the network of genes that interact with each directional anchor gene. These metrics, which we call a pharmagenic enrichment score (PES), identify individuals with a higher burden of genetic risk, localized in biological processes related to the candidate drug target, to inform precision drug repurposing. We used this approach to investigate schizophrenia and bipolar disorder and reveal several compounds targeting specific directional anchor genes that could be plausibly repurposed. These genetic risk scores, mapped to the networks associated with target genes, revealed biological insights that cannot be observed in undifferentiated genome-wide polygenic risk score (PRS). For example, an enrichment of these partitioned scores in schizophrenia cases with otherwise low PRS. In summary, genetic risk could be used more specifically to direct drug repurposing candidates that target particular genes implicated in psychiatric and other complex disorders.

Proteomic profiling platforms head to head: Leveraging genetics and clinical traits to compare aptamer- and antibody-based methods

High-throughput proteomic profiling using antibody or aptamer-based affinity reagents is used increasingly in human studies. However, direct analyses to address the relative strengths and weaknesses of these platforms are lacking. We assessed findings from the SomaScan1.3K (N = 1301 reagents), the SomaScan5K platform (N = 4979 reagents), and the Olink Explore (N = 1472 reagents) profiling techniques in 568 adults from the Jackson Heart Study and 219 participants in the HERITAGE Family Study across four performance domains: precision, accuracy, analytic breadth, and phenotypic associations leveraging detailed clinical phenotyping and genetic data. Across these studies, we show evidence supporting more reliable protein target specificity and a higher number of phenotypic associations for the Olink platform, while the Soma platforms benefit from greater measurement precision and analytic breadth across the proteome.

Phenotypic Causal Inference Using Genome-Wide Association Study Data: Mendelian Randomization and Beyond

statistics for genome-wide association studies (GWAS) are increasingly available for downstream analyses. Meanwhile, the popularity of causal inference methods has grown as we look to gather robust evidence for novel medical and public health interventions. This has led to the development of methods that use GWAS summary statistics for causal inference. Here, we describe these methods in order of their escalating complexity, from genetic associations to extensions of Mendelian randomization that consider thousands of phenotypes simultaneously. We also cover the assumptions and limitations of these approaches before considering the challenges faced by researchers performing causal inference using GWAS data. GWAS summary statistics constitute an important data source for causal inference research that offers a counterpoint to nongenetic methods when triangulating evidence. Continued efforts to address the challenges in using GWAS data for causal inference will allow the full impact of these approaches to be realized.

Proteomic discovery in diabetic kidney disease-to what end?

Emerging approaches that enable high-throughput analysis of the plasma proteome have been increasingly deployed in nephrology research. Kobayashi et al. provide a valuable addition to this literature, describing an untargeted proteomic analysis of diabetic kidney disease progression to end-stage kidney disease. This commentary places the study's findings in the context of the broader literature and outlines potential avenues toward biological insight and clinical utility.

Shared mechanisms across the major psychiatric and neurodegenerative diseases

Several common psychiatric and neurodegenerative diseases share epidemiologic risk; however, whether they share pathophysiology is unclear and is the focus of our investigation. Using 25 GWAS results and LD score regression, we find eight significant genetic correlations between psychiatric and neurodegenerative diseases. We integrate the GWAS results with human brain transcriptomes (n = 888) and proteomes (n = 722) to identify cis- and trans- transcripts and proteins that are consistent with a pleiotropic or causal role in each disease, referred to as causal proteins for brevity. Within each disease group, we find many distinct and shared causal proteins. Remarkably, 30% (13 of 42) of the neurodegenerative disease causal proteins are shared with psychiatric disorders. Furthermore, we find 2.6-fold more protein-protein interactions among the psychiatric and neurodegenerative causal proteins than expected by chance. Together, our findings suggest these psychiatric and neurodegenerative diseases have shared genetic and molecular pathophysiology, which has important ramifications for early treatment and therapeutic development.

Multi-ancestry Mendelian randomization of omics traits revealing drug targets of COVID-19 severity

BACKGROUND

Recent omic studies prioritised several drug targets associated with coronavirus disease 2019 (COVID-19) severity. However, little evidence was provided to systematically estimate the effect of drug targets on COVID-19 severity in multiple ancestries.

METHODS

In this study, we applied Mendelian randomization (MR) and colocalization approaches to understand the putative causal effects of 16,059 transcripts and 1608 proteins on COVID-19 severity in European and effects of 610 proteins on COVID-19 severity in African ancestry. We further integrated genetics, clinical and literature evidence to prioritise drug targets. Additional sensitivity analyses including multi-trait colocalization and phenome-wide MR were conducted to test for MR assumptions.

FINDINGS

MR and colocalization prioritized four protein targets, FCRL3, ICAM5, ENTPD5 and OAS1 that showed effect on COVID-19 severity in European ancestry. One protein target, SERPINA1 showed a stronger effect in African ancestry but much weaker effect in European ancestry (odds ratio [OR] in Africans=0.369, 95%CI=0.203 to 0.668, P = 9.96 × 10-4; OR in Europeans=1.021, 95%CI=0.901 to 1.157, P = 0.745), which suggested that increased level of SERPINA1 will reduce COVID-19 risk in African ancestry. One protein, ICAM1 showed suggestive effect on COVID-19 severity in both ancestries (OR in Europeans=1.152, 95%CI=1.063 to 1.249, P = 5.94 × 10-4; OR in Africans=1.481, 95%CI=1.008 to 2.176; P = 0.045). The OAS1, SERPINA1 and ICAM1 effects were replicated using updated COVID-19 severity data in the two ancestries respectively, where alternative splicing events in OAS1 and ICAM1 also showed marginal effects on COVID-19 severity in Europeans. The phenome-wide MR of the prioritised targets on 622 complex traits provided information on potential beneficial effects on other diseases and suggested little evidence of adverse effects on major complications.

INTERPRETATION

Our study identified six proteins as showing putative causal effects on COVID-19 severity. OAS1 and SERPINA1 were targets of existing drugs in trials as potential COVID-19 treatments. ICAM1, ICAM5 and FCRL3 are related to the immune system. Across the six targets, OAS1 has no reliable instrument in African ancestry; SERPINA1, FCRL3, ICAM5 and ENTPD5 showed a different level of putative causal evidence in European and African ancestries, which highlights the importance of more powerful ancestry-specific GWAS and value of multi-ancestry MR in informing the effects of drug targets on COVID-19 across different populations. This study provides a first step towards clinical investigation of beneficial and adverse effects of COVID-19 drug targets.

FUNDING

No.

Mining the Plasma Proteome for Insights into the Molecular Pathology of Pulmonary Arterial Hypertension

Rationale: Pulmonary arterial hypertension (PAH) is characterized by structural remodeling of pulmonary arteries and arterioles. Underlying biological processes are likely reflected in a perturbation of circulating proteins. Objectives: To quantify and analyze the plasma proteome of patients with PAH using inherited genetic variation to inform on underlying molecular drivers. Methods: An aptamer-based assay was used to measure plasma proteins in 357 patients with idiopathic or heritable PAH, 103 healthy volunteers, and 23 relatives of patients with PAH. In discovery and replication subgroups, the plasma proteomes of PAH and healthy individuals were compared, and the relationship to transplantation-free survival in PAH was determined. To examine causal relationships to PAH, protein quantitative trait loci (pQTL) that influenced protein levels in the patient population were used as instruments for Mendelian randomization (MR) analysis. Measurements and Main Results: From 4,152 annotated plasma proteins, levels of 208 differed between patients with PAH and healthy subjects, and 49 predicted long-term survival. MR based on cis-pQTL located in proximity to the encoding gene for proteins that were prognostic and distinguished PAH from health estimated an adverse effect for higher levels of netrin-4 (odds ratio [OR], 1.55; 95% confidence interval [CI], 1.16-2.08) and a protective effect for higher levels of thrombospondin-2 (OR, 0.83; 95% CI, 0.74-0.94) on PAH. Both proteins tracked the development of PAH in previously healthy relatives and changes in thrombospondin-2 associated with pulmonary arterial pressure at disease onset. Conclusions: Integrated analysis of the plasma proteome and genome implicates two secreted matrix-binding proteins, netrin-4 and thrombospondin-2, in the pathobiology of PAH.

Genome-Wide Mapping of Plasma IgG N-Glycan Quantitative Trait Loci Identifies a Potentially Causal Association between IgG N-Glycans and Rheumatoid Arthritis

Observational studies highlight associations of IgG N-glycosylation with rheumatoid arthritis (RA); however, the causality between these conditions remains to be determined. Standard and multivariable two-sample Mendelian randomization (MR) analyses integrating a summary genome-wide association study for RA and IgG N-glycan quantitative trait loci (IgG N-glycan-QTL) data were performed to explore the potentially causal associations of IgG N-glycosylation with RA. After correcting for multiple testing (p < 2 × 10^-3), the standard MR analysis based on the inverse-variance weighted method showed a significant association of genetically instrumented IgG N-glycan (GP4) with RA (odds ratio_GP4 = 0.906, 95% confidence interval = 0.857-0.958, p = 5.246 × 10^-4). In addition, we identified seven significant associations of genetically instrumented IgG N-glycans with RA by multivariable MR analysis (p < 2 × 10^-3). Results were broadly consistent in sensitivity analyses using MR_Lasso, MR_weighted median, MR_Egger regression, and leave-one-out analysis with different instruments (all p values <0.05). There was limited evidence of pleiotropy bias (all p values > 0.05). In conclusion, our MR analysis incorporating genome-wide association studies and IgG N-glycan-QTL data revealed that IgG N-glycans were potentially causally associated with RA. Our findings shed light on the role of IgG N-glycosylation in the development of RA. Future studies are needed to validate our findings and to explore the underlying physiological mechanisms in the etiology of RA.

Identifying causal genes for depression via integration of the proteome and transcriptome from brain and blood

Genome-wide association studies (GWASs) have identified numerous risk genes for depression. Nevertheless, genes crucial for understanding the molecular mechanisms of depression and effective antidepressant drug targets are largely unknown. Addressing this, we aimed to highlight potentially causal genes by systematically integrating the brain and blood protein and expression quantitative trait loci (QTL) data with a depression GWAS dataset via a statistical framework including Mendelian randomization (MR), Bayesian colocalization, and Steiger filtering analysis. In summary, we identified three candidate genes (TMEM106B, RAB27B, and GMPPB) based on brain data and two genes (TMEM106B and NEGR1) based on blood data with consistent robust evidence at both the protein and transcriptional levels. Furthermore, the protein-protein interaction (PPI) network provided new insights into the interaction between brain and blood in depression. Collectively, four genes (TMEM106B, RAB27B, GMPPB, and NEGR1) affect depression by influencing protein and gene expression level, which could guide future researches on candidate genes investigations in animal studies as well as prioritize antidepressant drug targets.

Mendelian randomization analysis of plasma levels of CD209 and MICB proteins and the risk of varicose veins of lower extremities

Varicose veins of lower extremities (VVs) are a highly prevalent condition, the pathogenesis of which is still not fully elucidated. Mendelian randomization (MR) can provide useful preliminary information on the traits that are potentially causally related to the disease. The aim of the present study is to replicate the effects of the plasma levels of MHC class I polypeptide-related sequence B (MICB) and cluster of differentiation 209 (CD209) proteins reported in a previous hypothesis-free MR study. We conducted MR analysis using a fixed effects inverse-variance weighted meta-analysis of Wald ratios method. For MICB and CD209, we used data from a large-scale genome-wide association study (GWAS) for plasma protein levels (N = 3,301). For VVs, we used GWAS data obtained in the FinnGen project (N = 128,698), the eMERGE network (phase 3, N = 48,429), and the UK Biobank data available in the Gene ATLAS (N = 452,264). The data used in the study were obtained in individuals of European descent. The results for MICB did not pass criteria for statistical significance and replication. The results for CD209 passed all statistical significance thresholds, indicating that the genetically predicted increase in CD209 level is associated with increased risk of VVs (βMR (SE) = 0.07 (0.01), OR (95% CI) = 1.08 (1.05-1.10), P-value = 5.9 ×10-11 in the meta-analysis of three cohorts). Our findings provide further support that CD209 can potentially be involved in VVs. In future studies, independent validation of our results using data from more powerful GWASs for CD209 measured by different methods would be beneficial.

A protocol for recruiting and analyzing the disease-oriented Russian disc degeneration study (RuDDS) biobank for functional omics studies of lumbar disc degeneration

Lumbar intervertebral disc degeneration (DD) disease is one of the main risk factors for low back pain and a leading cause of population absenteeism and disability worldwide. Despite a variety of biological studies, lumbar DD is not yet fully understood, partially because there are only few studies that use systematic and integrative approaches. This urges the need for studies that integrate different omics (including genomics and transcriptomics) measured on samples within a single cohort. This protocol describes a disease-oriented Russian disc degeneration study (RuDDS) biobank recruitment and analyses aimed to facilitate further omics studies of lumbar DD integrating genomic, transcriptomic and glycomic data. A total of 1,100 participants aged over 18 with available lumbar MRI scans, medical histories and biological material (whole blood, plasma and intervertebral disc tissue samples from surgically treated patients) will be enrolled during the three-year period from two Russian clinical centers. Whole blood, plasma and disc tissue specimens will be used for genotyping with genome-wide SNP-arrays, glycome profiling and RNA sequencing, respectively. Omics data will be further used for a genome-wide association study of lumbar DD with in silico functional annotation, analysis of plasma glycome and lumbar DD disease interactions and transcriptomic data analysis including an investigation of differential expression patterns associated with lumbar DD disease. Statistical tests applied in each of the analyses will meet the standard criteria specific to the attributed study field. In a long term, the results of the study will expand fundamental knowledge about lumbar DD development and contribute to the elaboration of novel personalized approaches for disease prediction and therapy. Additionally to the lumbar disc degeneration study, a RuDDS cohort could be used for other genetic studies, as it will have unique omics data.

Trial registration numberNCT04600544.

Combining evidence from Mendelian randomization and colocalization: Review and comparison of approaches

Mendelian randomization and colocalization are two statistical approaches that can be applied to summarized data from genome-wide association studies (GWASs) to understand relationships between traits and diseases. However, despite similarities in scope, they are different in their objectives, implementation, and interpretation, in part because they were developed to serve different scientific communities. Mendelian randomization assesses whether genetic predictors of an exposure are associated with the outcome and interprets an association as evidence that the exposure has a causal effect on the outcome, whereas colocalization assesses whether two traits are affected by the same or distinct causal variants. When considering genetic variants in a single genetic region, both approaches can be performed. While a positive colocalization finding typically implies a non-zero Mendelian randomization estimate, the reverse is not generally true: there are several scenarios which would lead to a non-zero Mendelian randomization estimate but lack evidence for colocalization. These include the existence of distinct but correlated causal variants for the exposure and outcome, which would violate the Mendelian randomization assumptions, and a lack of strong associations with the outcome. As colocalization was developed in the GWAS tradition, typically evidence for colocalization is concluded only when there is strong evidence for associations with both traits. In contrast, a non-zero estimate from Mendelian randomization can be obtained despite only nominally significant genetic associations with the outcome at the locus. In this review, we discuss how the two approaches can provide complementary information on potential therapeutic targets.

Plasma proteome analyses in individuals of European and African ancestry identify cis-pQTLs and models for proteome-wide association studies

Improved understanding of genetic regulation of the proteome can facilitate identification of the causal mechanisms for complex traits. We analyzed data on 4,657 plasma proteins from 7,213 European American (EA) and 1,871 African American (AA) individuals from the Atherosclerosis Risk in Communities study, and further replicated findings on 467 AA individuals from the African American Study of Kidney Disease and Hypertension study. Here, we identified 2,004 proteins in EA and 1,618 in AA, with most overlapping, which showed associations with common variants in cis-regions. Availability of AA samples led to smaller credible sets and notable number of population-specific cis-protein quantitative trait loci. Elastic Net produced powerful models for protein prediction in both populations. An application of proteome-wide association studies to serum urate and gout implicated several proteins, including IL1RN, revealing the promise of the drug anakinra to treat acute gout flares. Our study demonstrates the value of large and diverse ancestry study to investigate the genetic mechanisms of molecular phenotypes and their relationship with complex traits.

Identifying factors contributing to increased susceptibility to COVID-19 risk: a systematic review of Mendelian randomization studies

Background

To summarize modifiable factors for coronavirus disease 2019 (COVID-19) suggested by Mendelian randomization studies.

Methods

In this systematic review, we searched PubMed, EMBASE and MEDLINE, from inception to 15 November 2021, for Mendelian randomization studies in English. We selected studies that assessed associations of genetically predicted exposures with COVID-19-related outcomes (severity, hospitalization and susceptibility). Risk of bias of the included studies was evaluated based on the consideration of the three main assumptions for instrumental variable analyses.

Results

We identified 700 studies through systematic search, of which 50 Mendelian randomization studies were included. Included studies have explored a wide range of socio-demographic factors, lifestyle attributes, anthropometrics and biomarkers, predisposition to diseases and druggable targets in COVID-19 risk. Mendelian randomization studies suggested that increases in smoking, obesity and inflammatory factors were associated with higher risk of COVID-19. Predisposition to ischaemic stroke, combined bipolar disorder and schizophrenia, attention-deficit and hyperactivity disorder, chronic kidney disease and idiopathic pulmonary fibrosis was potentially associated with higher COVID-19 risk. Druggable targets, such as higher protein expression of histo-blood group ABO system transferase (ABO), interleukin (IL)-6 and lower protein expression of 2′-5′ oligoadenylate synthetase 1 (OAS1) were associated with higher risk of COVID-19. There was no strong genetic evidence supporting the role of vitamin D, glycaemic traits and predisposition to cardiometabolic diseases in COVID-19 risk.

Conclusion

This review summarizes modifiable factors for intervention (e.g. smoking, obesity and inflammatory factors) and proteomic signatures (e.g. OAS1 and IL-6) that could help identify drugs for treating COVID-19.

Mendelian Randomization Analysis Reveals Causal Effects of Plasma Proteome on Body Composition Traits

CONTEXT

Observational studies have demonstrated associations between plasma proteins and obesity, but evidence of causal relationship remains to be studied.

OBJECTIVE

We aimed to evaluate the causal relationship between plasma proteins and body composition.

METHODS

We conducted a 2-sample Mendelian randomization (MR) analysis based on the genome-wide association study (GWAS) summary statistics of 23 body composition traits and 2656 plasma proteins. We then performed hierarchical cluster analysis to evaluate the structure and pattern of the identified causal associations, and we performed gene ontology enrichment analysis to explore the functional relevance of the identified proteins.

RESULTS

We identified 430 putatively causal effects of 96 plasma proteins on 22 body composition traits (except obesity status) with strong MR evidence (P < 2.53 × 10 - 6, at a Bonferroni-corrected threshold). The top 3 causal associations are follistatin (FST) on trunk fat-free mass (Beta = -0.63, SE = 0.04, P = 2.00 × 10-63), insulin-like growth factor-binding protein 1 (IGFBP1) on trunk fat-free mass (Beta = -0.54, SE = 0.03, P = 1.79 × 10-57) and r-spondin-3 (RSPO3) on WHR (waist circumference/hip circumference) (Beta = 0.01, SE = 4.47 × 10-4, P = 5.45 × 10-60), respectively. Further clustering analysis and pathway analysis demonstrated that the pattern of causal effect to fat mass and fat-free mass may be different.

CONCLUSION

Our findings may provide evidence for causal relationships from plasma proteins to various body composition traits and provide basis for further targeted functional studies.

Genetics of osteopontin in patients with chronic kidney disease: The German Chronic Kidney Disease study

Osteopontin (OPN), encoded by SPP1, is a phosphorylated glycoprotein predominantly synthesized in kidney tissue. Increased OPN mRNA and protein expression correlates with proteinuria, reduced creatinine clearance, and kidney fibrosis in animal models of kidney disease. But its genetic underpinnings are incompletely understood. We therefore conducted a genome-wide association study (GWAS) of OPN in a European chronic kidney disease (CKD) population. Using data from participants of the German Chronic Kidney Disease (GCKD) study (N = 4,897), a GWAS (minor allele frequency [MAF]≥1%) and aggregated variant testing (AVT, MAF<1%) of ELISA-quantified serum OPN, adjusted for age, sex, estimated glomerular filtration rate (eGFR), and urinary albumin-to-creatinine ratio (UACR) was conducted. In the project, GCKD participants had a mean age of 60 years (SD 12), median eGFR of 46 mL/min/1.73m2 (p25: 37, p75: 57) and median UACR of 50 mg/g (p25: 9, p75: 383). GWAS revealed 3 loci (p<5.0E-08), two of which replicated in the population-based Young Finns Study (YFS) cohort (p<1.67E-03): rs10011284, upstream of SPP1 encoding the OPN protein and related to OPN production, and rs4253311, mapping into KLKB1 encoding prekallikrein (PK), which is processed to kallikrein (KAL) implicated through the kinin-kallikrein system (KKS) in blood pressure control, inflammation, blood coagulation, cancer, and cardiovascular disease. The SPP1 gene was also identified by AVT (p = 2.5E-8), comprising 7 splice-site and missense variants. Among others, downstream analyses revealed colocalization of the OPN association signal at SPP1 with expression in pancreas tissue, and at KLKB1 with various plasma proteins in trans, and with phenotypes (bone disorder, deep venous thrombosis) in human tissue. In summary, this GWAS of OPN levels revealed two replicated associations. The KLKB1 locus connects the function of OPN with PK, suggestive of possible further post-translation processing of OPN. Further studies are needed to elucidate the complex role of OPN within human (patho)physiology.

A proteomic surrogate for cardiovascular outcomes that is sensitive to multiple mechanisms of change in risk

A reliable, individualized, and dynamic surrogate of cardiovascular risk, synoptic for key biologic mechanisms, could shorten the path for drug development, enhance drug cost-effectiveness and improve patient outcomes. We used highly multiplexed proteomics to address these objectives, measuring about 5000 proteins in each of 32,130 archived plasma samples from 22,849 participants in nine clinical studies. We used machine learning to derive a 27-protein model predicting 4-year likelihood of myocardial infarction, stroke, heart failure, or death. The 27 proteins encompassed 10 biologic systems, and 12 were associated with relevant causal genetic traits. We independently validated results in 11,609 participants. Compared to a clinical model, the ratio of observed events in quintile 5 to quintile 1 was 6.7 for proteins versus 2.9 for the clinical model, AUCs (95% CI) were 0.73 (0.72 to 0.74) versus 0.64 (0.62 to 0.65), c-statistics were 0.71 (0.69 to 0.72) versus 0.62 (0.60 to 0.63), and the net reclassification index was +0.43. Adding the clinical model to the proteins only improved discrimination metrics by 0.01 to 0.02. Event rates in four predefined protein risk categories were 5.6, 11.2, 20.0, and 43.4% within 4 years; median time to event was 1.71 years. Protein predictions were directionally concordant with changed outcomes. Adverse risks were predicted for aging, approaching an event, anthracycline chemotherapy, diabetes, smoking, rheumatoid arthritis, cancer history, cardiovascular disease, high systolic blood pressure, and lipids. Reduced risks were predicted for weight loss and exenatide. The 27-protein model has potential as a "universal" surrogate end point for cardiovascular risk.

Genetic regulation of post-translational modification of two distinct proteins

Post-translational modifications diversify protein functions and dynamically coordinate their signalling networks, influencing most aspects of cell physiology. Nevertheless, their genetic regulation or influence on complex traits is not fully understood. Here, we compare the genetic regulation of the same PTM of two proteins - glycosylation of transferrin and immunoglobulin G (IgG). By performing genome-wide association analysis of transferrin glycosylation, we identify 10 significantly associated loci, 9 of which were not reported previously. Comparing these with IgG glycosylation-associated genes, we note protein-specific associations with genes encoding glycosylation enzymes (transferrin - MGAT5, ST3GAL4, B3GAT1; IgG - MGAT3, ST6GAL1), as well as shared associations (FUT6, FUT8). Colocalisation analyses of the latter suggest that different causal variants in the FUT genes regulate fucosylation of the two proteins. Glycosylation of these proteins is thus genetically regulated by both shared and protein-specific mechanisms.

Enhancer promoter interactome and Mendelian randomization identify network of druggable vascular genes in coronary artery disease

Coronary artery disease (CAD) is a multifactorial disorder, which is partly heritable. Herein, we implemented a mapping of CAD-associated candidate genes by using genome-wide enhancer-promoter conformation (H3K27ac-HiChIP) and expression quantitative trait loci (eQTL). Enhancer-promoter anchor loops from human coronary artery smooth muscle cells (HCASMC) explained 22% of the heritability for CAD. 3D enhancer-promoter genome mapping of CAD-genes in HCASMC was enriched in vascular eQTL genes. By using colocalization and Mendelian randomization analyses, we identified 58 causal candidate vascular genes including some druggable targets (MAP3K11, CAMK1D, PDGFD, IPO9 and CETP). A network analysis of causal candidate genes was enriched in TGF beta and MAPK pathways. The pharmacologic inhibition of causal candidate gene MAP3K11 in vascular SMC reduced the expression of athero-relevant genes and lowered cell migration, a cardinal process in CAD. Genes connected to enhancers are enriched in vascular eQTL and druggable genes causally associated with CAD.

Systematic review of Mendelian randomization studies on risk of cancer

BACKGROUND

We aimed to map and describe the current state of Mendelian randomization (MR) literature on cancer risk and to identify associations supported by robust evidence.

METHODS

We searched PubMed and Scopus up to 06/10/2020 for MR studies investigating the association of any genetically predicted risk factor with cancer risk. We categorized the reported associations based on a priori designed levels of evidence supporting a causal association into four categories, namely robust, probable, suggestive, and insufficient, based on the significance and concordance of the main MR analysis results and at least one of the MR-Egger, weighed median, MRPRESSO, and multivariable MR analyses. Associations not presenting any of the aforementioned sensitivity analyses were not graded.

RESULTS

We included 190 publications reporting on 4667 MR analyses. Most analyses (3200; 68.6%) were not accompanied by any of the assessed sensitivity analyses. Of the 1467 evaluable analyses, 87 (5.9%) were supported by robust, 275 (18.7%) by probable, and 89 (6.1%) by suggestive evidence. The most prominent robust associations were observed for anthropometric indices with risk of breast, kidney, and endometrial cancers; circulating telomere length with risk of kidney, lung, osteosarcoma, skin, thyroid, and hematological cancers; sex steroid hormones and risk of breast and endometrial cancer; and lipids with risk of breast, endometrial, and ovarian cancer.

CONCLUSIONS

Despite the large amount of research on genetically predicted risk factors for cancer risk, limited associations are supported by robust evidence for causality. Most associations did not present a MR sensitivity analysis and were thus non-evaluable. Future research should focus on more thorough assessment of sensitivity MR analyses and on more transparent reporting.

Coding and regulatory variants are associated with serum protein levels and disease

Circulating proteins can be used to diagnose and predict disease-related outcomes. A deep serum proteome survey recently revealed close associations between serum protein networks and common disease. In the current study, 54,469 low-frequency and common exome-array variants were compared to 4782 protein measurements in the serum of 5343 individuals from the AGES Reykjavik cohort. This analysis identifies a large number of serum proteins with genetic signatures overlapping those of many diseases. More specifically, using a study-wide significance threshold, we find that 2021 independent exome array variants are associated with serum levels of 1942 proteins. These variants reside in genetic loci shared by hundreds of complex disease traits, highlighting serum proteins' emerging role as biomarkers and potential causative agents of a wide range of diseases.

Little genomic support for Cyclophilin A-matrix metalloproteinase-9 pathway as a therapeutic target for cognitive impairment in APOE4 carriers

Therapeutic targets for halting the progression of Alzheimer's disease pathology are lacking. Recent evidence suggests that APOE4, but not APOE3, activates the Cyclophilin-A matrix metalloproteinase-9 (CypA-MMP9) pathway, leading to an accelerated breakdown of the blood-brain barrier (BBB) and thereby causing neuronal and synaptic dysfunction. Furthermore, blockade of the CypA-MMP9 pathway in APOE4 knock-in mice restores BBB integrity and subsequently normalizes neuronal and synaptic function. Thus, CypA has been suggested as a potential target for treating APOE4 mediated neurovascular injury and the resulting neuronal dysfunction and degeneration. The odds of drug targets passing through clinical trials are greatly increased if they are supported by genomic evidence. We found little evidence to suggest that CypA or MMP9 affects the risk of Alzheimer's disease or cognitive impairment using two-sample Mendelian randomization and polygenic risk score analysis in humans. This casts doubt on whether they are likely to represent effective drug targets for cognitive impairment in human APOE4 carriers.

Circulating inflammatory cytokines and risk of five cancers: a Mendelian randomization analysis

BACKGROUND

Epidemiological and experimental evidence has linked chronic inflammation to cancer aetiology. It is unclear whether associations for specific inflammatory biomarkers are causal or due to bias. In order to examine whether altered genetically predicted concentration of circulating cytokines are associated with cancer development, we performed a two-sample Mendelian randomisation (MR) analysis.

METHODS

Up to 31,112 individuals of European descent were included in genome-wide association study (GWAS) meta-analyses of 47 circulating cytokines. Single nucleotide polymorphisms (SNPs) robustly associated with the cytokines, located in or close to their coding gene (cis), were used as instrumental variables. Inverse-variance weighted MR was used as the primary analysis, and the MR assumptions were evaluated in sensitivity and colocalization analyses and a false discovery rate (FDR) correction for multiple comparisons was applied. Corresponding germline GWAS summary data for five cancer outcomes (breast, endometrial, lung, ovarian, and prostate), and their subtypes were selected from the largest cancer-specific GWASs available (cases ranging from 12,906 for endometrial to 133,384 for breast cancer).

RESULTS

There was evidence of inverse associations of macrophage migration inhibitory factor with breast cancer (OR per SD = 0.88, 95% CI 0.83 to 0.94), interleukin-1 receptor antagonist with endometrial cancer (0.86, 0.80 to 0.93), interleukin-18 with lung cancer (0.87, 0.81 to 0.93), and beta-chemokine-RANTES with ovarian cancer (0.70, 0.57 to 0.85) and positive associations of monokine induced by gamma interferon with endometrial cancer (3.73, 1.86 to 7.47) and cutaneous T-cell attracting chemokine with lung cancer (1.51, 1.22 to 1.87). These associations were similar in sensitivity analyses and supported in colocalization analyses.

CONCLUSIONS

Our study adds to current knowledge on the role of specific inflammatory biomarker pathways in cancer aetiology. Further validation is needed to assess the potential of these cytokines as pharmacological or lifestyle targets for cancer prevention.

Review of Mendelian Randomization Studies on Endometrial Cancer

Endometrial cancer (EC) is a common gynecological cancer. In some parts of the world, the incidence and mortality of EC are on the rise. Understanding the risk factors of EC is necessary to prevent the occurrence of this disease. Observational studies have revealed the association between certain modifiable environmental risk factors and EC risk. However, due to unmeasured confounding, measurement errors, and reverse causality, observational studies sometimes have limited ability to judge robust causal inferences. In recent years, Mendelian randomization (MR) analysis has received extensive attention, providing valuable insights for cancer-related research, and is expected to identify potential therapeutic interventions. In MR analysis, genetic variation (alleles are randomly assigned during meiosis and are usually independent of environmental or lifestyle factors) is used instead of modifiable exposure to study the relationship between risk factors and disease. Therefore, MR analysis can make causal inference about exposure and disease risk. This review briefly describes the key principles and assumptions of MR analysis; summarizes published MR studies on EC; focuses on the correlation between different risk factors and EC risks; and discusses the application of MR methods in EC research. The results of MR studies on EC showed that type 2 diabetes, uterine fibroids, higher body mass index, higher plasminogen activator inhibitor-1 (PAI-1), higher fasting insulin, early insulin secretion, longer telomere length, higher testosterone and higher plasma cortisol levels are associated with increased risk of EC. In contrast, later age of menarche, higher circulatory tumor necrosis factor, higher low-density lipoprotein cholesterol, and higher sex hormone-binding globulin levels are associated with reduced risk of EC. In general, despite some limitations, MR analysis still provides an effective way to explore the causal relationship between different risk factors and EC.

Mendelian randomization of genetically independent aging phenotypes identifies LPA and VCAM1 as biological targets for human aging

Length and quality of life are important to us all, yet identification of promising drug targets for human aging using genetics has had limited success. In the present study, we combine six European-ancestry genome-wide association studies of human aging traits-healthspan, father and mother lifespan, exceptional longevity, frailty index and self-rated health-in a principal component framework that maximizes their shared genetic architecture. The first principal component (aging-GIP1) captures both length of life and indices of mental and physical wellbeing. We identify 27 genomic regions associated with aging-GIP1, and provide additional, independent evidence for an effect on human aging for loci near HTT and MAML3 using a study of Finnish and Japanese survival. Using proteome-wide, two-sample, Mendelian randomization and colocalization, we provide robust evidence for a detrimental effect of blood levels of apolipoprotein(a) and vascular cell adhesion molecule 1 on aging-GIP1. Together, our results demonstrate that combining multiple aging traits using genetic principal components enhances the power to detect biological targets for human aging.

Clinically Relevant Circulating Protein Biomarkers for Type 1 Diabetes: Evidence From a Two-Sample Mendelian Randomization Study

OBJECTIVE

To identify circulating proteins influencing type 1 diabetes susceptibility using Mendelian randomization (MR).

RESEARCH DESIGN AND METHODS

We used a large-scale two-sample MR study, using cis genetic determinants (protein quantitative trait loci [pQTL]) of up to 1,611 circulating proteins from five large genome-wide association studies, to screen for causal associations of these proteins with type 1 diabetes risk in 9,684 case subjects with type 1 diabetes and 15,743 control subjects. Further, pleiotropy-robust MR methods were used in sensitivity analyses using both cis and trans-pQTL.

RESULTS

We found that a genetically predicted SD increase in signal regulatory protein gamma (SIRPG) level was associated with increased risk of type 1 diabetes risk (MR odds ratio [OR] 1.66 [95% 1.36-2.03]; P = 7.1 × 10-7). The risk of type 1 diabetes increased almost twofold per genetically predicted standard deviation (SD) increase in interleukin-27 Epstein-Barr virus-induced 3 (IL27-EBI3) protein levels (MR OR 1.97 [95% CI 1.48-2.62]; P = 3.7 × 10-6). However, an SD increase in chymotrypsinogen B1 (CTRB1) was associated with decreased risk of type 1 diabetes (MR OR 0.84 [95% CI 0.77-0.90]; P = 6.1 × 10-6). Sensitivity analyses using MR methods testing for pleiotropy while including trans-pQTL showed similar results. While the MR-Egger suggested no pleotropic effect (P value MR-Egger intercept = 0.31), there was evidence of pleiotropy in MR-PRESSO (P value global test = 0.006).

CONCLUSIONS

We identified three novel circulating protein biomarkers associated with type 1 diabetes risk using an MR approach. These biomarkers are promising targets for development of drugs and/or of screening tools for early prediction of type 1 diabetes.

Target Discovery for Drug Development Using Mendelian Randomization

Making drug development more efficient by identifying promising drug targets can contribute to resource savings. Identifying promising drug targets using human genetic approaches can remove barriers related to translation. In addition, genetic information can be used to identify potentially causal relationships between a drug target and disease. Mendelian randomization (MR) is a class of approaches used to identify causal associations between pairs of genetically predicted traits using data from human genetic studies. MR can be used to prioritize candidate drug targets by predicting disease outcomes and adverse events that could result from the manipulation of a drug target. The theory behind MR is reviewed, including a discussion of MR assumptions, different MR analytical methods, tests for violations of assumptions, and MR methods that can be robust to some violations of MR assumptions. A protocol to perform two-sample MR (2SMR) with summary genome-wide association study (GWAS) results is described. An example of 2SMR examining the causal relationship between low-density lipoprotein (LDL) and coronary artery disease (CAD) is provided as an illustration of the protocol.

Putative Candidate Drug Targets for Sarcopenia-Related Traits Identified Through Mendelian Randomization Analysis of the Blood Proteome

Purpose: The increasing prevalence of sarcopenia remains an ongoing challenge to health care systems worldwide. The lack of treatments encouraged the discovery of human proteomes to find potential therapeutic targets. As one of the major components of the human proteome, plasma proteins are functionally connected with various organs of the body to regulate biological processes and mediate overall homeostasis, which makes it crucial in various complex processes such as aging and chronic diseases. By performing a systematic causal analysis of the plasma proteome, we attempt to reveal the etiological mechanism and discover drug targets for sarcopenia. Methods: By using data from four genome-wide association studies for blood proteins and the UK Biobank data for sarcopenia-related traits, we applied two-sample Mendelian randomization (MR) analysis to evaluate 310 plasma proteins as possible causal mediators of sarcopenia-related traits: appendicular lean mass (ALM) and handgrip strength (right and left). Then we performed a two-sample bidirectional Mendelian randomization analysis for the identified putatively causal proteins to assess potential reverse causality that the trait values may influence protein levels. Finally, we performed phenome-wide MR analysis of the identified putatively causal proteins for 784 diseases to test the possible side effects of these proteins on other diseases. Results: Five plasma proteins were identified as putatively causal mediators of sarcopenia-related traits. Specifically, leukocyte immunoglobulin-like receptor subfamily B member 2 (LILRB2), asporin (ASPN), and contactin-2 (CNTN2) had potential causal effects on appendicular lean mass, and ecto-ADP-ribosyltransferase 4 (ART4) and superoxide dismutase 2 (SOD2) had putative causal effects on the handgrip strength, respectively. None of the five putatively causal proteins had a reverse causality relationship with sarcopenia-related traits, and no side effects on other diseases were identified. Conclusion: We identified five plasma proteins that may serve as putatively potential novel drug targets for sarcopenia. Our study attested to the value of two-sample MR analysis in identifying and prioritizing putatively potential therapeutic targets for complex diseases.

Finding genetically-supported drug targets for Parkinson's disease using Mendelian randomization of the druggable genome

Parkinson's disease is a neurodegenerative movement disorder that currently has no disease-modifying treatment, partly owing to inefficiencies in drug target identification and validation. We use Mendelian randomization to investigate over 3,000 genes that encode druggable proteins and predict their efficacy as drug targets for Parkinson's disease. We use expression and protein quantitative trait loci to mimic exposure to medications, and we examine the causal effect on Parkinson's disease risk (in two large cohorts), age at onset and progression. We propose 23 drug-targeting mechanisms for Parkinson's disease, including four possible drug repurposing opportunities and two drugs which may increase Parkinson's disease risk. Of these, we put forward six drug targets with the strongest Mendelian randomization evidence. There is remarkably little overlap between our drug targets to reduce Parkinson's disease risk versus progression, suggesting different molecular mechanisms. Drugs with genetic support are considerably more likely to succeed in clinical trials, and we provide compelling genetic evidence and an analysis pipeline to prioritise Parkinson's disease drug development.

Large-Scale Plasma Protein Profiling of Incident Myocardial Infarction, Ischemic Stroke, and Heart Failure

Background

We recently reported a link between plasma levels of 2 of 84 cardiovascular disease (CVD)–related proteins and the 3 major CVDs, myocardial infarction, ischemic stroke, and heart failure. The present study investigated whether measurement of almost 10 times the number of proteins could lead to discovery of additional risk markers for CVD.

Methods and Results

We measured 742 proteins using the proximity extension assay in 826 male participants of ULSAM (Uppsala Longitudinal Study of Adult Men) who were free from CVD at the age of 70 years. Cox proportional hazards models were adjusted for age only, as well as all traditional risk factors. During a 12.5‐year median follow‐up (maximal, 22.0 years), 283 incident CVDs occurred. Forty‐one proteins were significantly (false discovery rate <0.05) related to the combined end point of incident CVD, with N‐terminal pro–brain natriuretic peptide as the top finding, while 53 proteins were related to incident myocardial infarction. A total of 13 and 16 proteins were significantly related to incident ischemic stroke and heart failure, respectively. Growth differentiation factor 15, 4‐disulfide core domain protein 2, and kidney injury molecule were related to all of the 3 major CVD outcomes. A lasso selection of 11 proteins improved discrimination of incident CVD by 5.0% (P=0.0038).

Conclusions

Large‐scale proteomics seem useful for the discovery of new risk markers for CVD and to improve risk prediction in an elderly population of men. Further studies are needed to replicate the findings in independent samples of both men and women of different ages.

Large-scale integration of the plasma proteome with genetics and disease

The plasma proteome can help bridge the gap between the genome and diseases. Here we describe genome-wide association studies (GWASs) of plasma protein levels measured with 4,907 aptamers in 35,559 Icelanders. We found 18,084 associations between sequence variants and levels of proteins in plasma (protein quantitative trait loci; pQTL), of which 19% were with rare variants (minor allele frequency (MAF) < 1%). We tested plasma protein levels for association with 373 diseases and other traits and identified 257,490 associations. We integrated pQTL and genetic associations with diseases and other traits and found that 12% of 45,334 lead associations in the GWAS Catalog are with variants in high linkage disequilibrium with pQTL. We identified 938 genes encoding potential drug targets with variants that influence levels of possible biomarkers. Combining proteomics, genomics and transcriptomics, we provide a valuable resource that can be used to improve understanding of disease pathogenesis and to assist with drug discovery and development.

Crosstalk between Statins and Cancer Prevention and Therapy: An Update

The importance of statins in cancer has been discussed in many studies. They are known for their anticancer properties against solid tumors of the liver or lung, as well as diffuse cancers, such as multiple myeloma or leukemia. Currently, the most commonly used statins are simvastatin, rosuvastatin and atorvastatin. The anti-tumor activity of statins is largely related to their ability to induce apoptosis by targeting cancer cells with high selectivity. Statins are also involved in the regulation of the histone acetylation level, the disturbance of which can lead to abnormal activity of genes involved in the regulation of proliferation, differentiation and apoptosis. As a result, tumor growth and its invasion may be promoted, which is associated with a poor prognosis. High levels of histone deacetylases are observed in many cancers; therefore, one of the therapeutic strategies is to use their inhibitors. Combining statins with histone deacetylase inhibitors can induce a synergistic anticancer effect.

Synergistic insights into human health from aptamer- and antibody-based proteomic profiling

Affinity-based proteomics has enabled scalable quantification of thousands of protein targets in blood enhancing biomarker discovery, understanding of disease mechanisms, and genetic evaluation of drug targets in humans through protein quantitative trait loci (pQTLs). Here, we integrate two partly complementary techniques-the aptamer-based SomaScan^® v4 assay and the antibody-based Olink assays-to systematically assess phenotypic consequences of hundreds of pQTLs discovered for 871 protein targets across both platforms. We create a genetically anchored cross-platform proteome-phenome network comprising 547 protein-phenotype connections, 36.3% of which were only seen with one of the two platforms suggesting that both techniques capture distinct aspects of protein biology. We further highlight discordance of genetically predicted effect directions between assays, such as for PILRA and Alzheimer's disease. Our results showcase the synergistic nature of these technologies to better understand and identify disease mechanisms and provide a benchmark for future cross-platform discoveries.

Mapping the proteo-genomic convergence of human diseases

Characterization of the genetic regulation of proteins is essential for understanding disease etiology and developing therapies. We identified 10,674 genetic associations for 3892 plasma proteins to create a cis-anchored gene-protein-disease map of 1859 connections that highlights strong cross-disease biological convergence. This proteo-genomic map provides a framework to connect etiologically related diseases, to provide biological context for new or emerging disorders, and to integrate different biological domains to establish mechanisms for known gene-disease links. Our results identify proteo-genomic connections within and between diseases and establish the value of cis-protein variants for annotation of likely causal disease genes at loci identified in genome-wide association studies, thereby addressing a major barrier to experimental validation and clinical translation of genetic discoveries.

Pairwise effects between lipid GWAS genes modulate lipid plasma levels and cellular uptake

Complex traits are characterized by multiple genes and variants acting simultaneously on a phenotype. However, studying the contribution of individual pairs of genes to complex traits has been challenging since human genetics necessitates very large population sizes, while findings from model systems do not always translate to humans. Here, we combine genetics with combinatorial RNAi (coRNAi) to systematically test for pairwise additive effects (AEs) and genetic interactions (GIs) between 30 lipid genome-wide association studies (GWAS) genes. Gene-based burden tests from 240,970 exomes show that in carriers with truncating mutations in both, APOB and either PCSK9 or LPL ("human double knock-outs") plasma lipid levels change additively. Genetics and coRNAi identify overlapping AEs for 12 additional gene pairs. Overlapping GIs are observed for TOMM40/APOE with SORT1 and NCAN. Our study identifies distinct gene pairs that modulate plasma and cellular lipid levels primarily via AEs and nominates putative drug target pairs for improved lipid-lowering combination therapies.

Integrative analysis of the plasma proteome and polygenic risk of cardiometabolic diseases

Cardiometabolic diseases are frequently polygenic in architecture, comprising a large number of risk alleles with small effects spread across the genome1-3. Polygenic scores (PGS) aggregate these into a metric representing an individual's genetic predisposition to disease. PGS have shown promise for early risk prediction4-7 and there is an open question as to whether PGS can also be used to understand disease biology8. Here, we demonstrate that cardiometabolic disease PGS can be used to elucidate the proteins underlying disease pathogenesis. In 3,087 healthy individuals, we found that PGS for coronary artery disease, type 2 diabetes, chronic kidney disease and ischaemic stroke are associated with the levels of 49 plasma proteins. Associations were polygenic in architecture, largely independent of cis and trans protein quantitative trait loci and present for proteins without quantitative trait loci. Over a follow-up of 7.7 years, 28 of these proteins associated with future myocardial infarction or type 2 diabetes events, 16 of which were mediators between polygenic risk and incident disease. Twelve of these were druggable targets with therapeutic potential. Our results demonstrate the potential for PGS to uncover causal disease biology and targets with therapeutic potential, including those that may be missed by approaches utilizing information at a single locus.

Disentangling sex differences in the shared genetic architecture of posttraumatic stress disorder, traumatic experiences, and social support with body size and composition

There is a well-known association of traumatic experiences and posttraumatic stress disorder (PTSD) with body size and composition, including consistent differences between sexes. However, the biology underlying these associations is unclear. To understand the genetic underpinnings of this complex relationship, we investigated genome-wide datasets informative of African and European ancestries from the Psychiatric Genomic Consortium, the UK Biobank, the GIANT Consortium, and the Million Veteran Program. We used genome-wide association statistics to estimate sex-specific genetic correlations (r g ) of traumatic experiences, social support, and PTSD with multiple anthropometric traits. After multiple testing corrections (false discovery rate, FDR q < 0.05), we observed 58 significant r g relationships in females (e.g., childhood physical abuse and body mass index, BMI r g  = 0.245, p = 3.88 × 10-10) and 21 significant r g relationships in males (e.g., been involved in combat or exposed to warzone and leg fat percentage; r g  = 0.405, p = 4.42 × 10-10). We performed causal inference analyses of these genetic overlaps using Mendelian randomization and latent causal variable approaches. Multiple female-specific putative causal relationships were observed linking body composition/size with PTSD (e.g., leg fat percentage→PTSD; beta = 0.319, p = 3.13 × 10-9), traumatic experiences (e.g., childhood physical abuse→waist circumference; beta = 0.055, p = 5.07 × 10-4), and childhood neglect (e.g., "someone to take you to doctor when needed as a child"→BMI; beta = -0.594, p = 1.09 × 10-5). In males, we observed putative causal effects linking anthropometric-trait genetic liabilities to traumatic experiences (e.g., BMI→childhood physical abuse; beta = 0.028, p = 8.19 × 10-3). Some of these findings were replicated in individuals of African descent although the limited sample size available did not permit us to conduct a sex-stratified analysis in this ancestry group. In conclusion, our findings provide insights regarding sex-specific causal networks linking anthropometric traits to PTSD, traumatic experiences, and social support.

Validation of lipid-related therapeutic targets for coronary heart disease prevention using human genetics

Drug target Mendelian randomization (MR) studies use DNA sequence variants in or near a gene encoding a drug target, that alter the target's expression or function, as a tool to anticipate the effect of drug action on the same target. Here we apply MR to prioritize drug targets for their causal relevance for coronary heart disease (CHD). The targets are further prioritized using independent replication, co-localization, protein expression profiles and data from the British National Formulary and clinicaltrials.gov. Out of the 341 drug targets identified through their association with blood lipids (HDL-C, LDL-C and triglycerides), we robustly prioritize 30 targets that might elicit beneficial effects in the prevention or treatment of CHD, including NPC1L1 and PCSK9, the targets of drugs used in CHD prevention. We discuss how this approach can be generalized to other targets, disease biomarkers and endpoints to help prioritize and validate targets during the drug development process.

Advancing the use of genome-wide association studies for drug repurposing

Genome-wide association studies (GWAS) have revealed important biological insights into complex diseases, which are broadly expected to lead to the identification of new drug targets and opportunities for treatment. Drug development, however, remains hampered by the time taken and costs expended to achieve regulatory approval, leading many clinicians and researchers to consider alternative paths to more immediate clinical outcomes. In this Review, we explore approaches that leverage common variant genetics to identify opportunities for repurposing existing drugs, also known as drug repositioning. These approaches include the identification of compounds by linking individual loci to genes and pathways that can be pharmacologically modulated, transcriptome-wide association studies, gene-set association, causal inference by Mendelian randomization, and polygenic scoring.

Effects of adiposity on the human plasma proteome: observational and Mendelian randomisation estimates

BACKGROUND

Variation in adiposity is associated with cardiometabolic disease outcomes, but mechanisms leading from this exposure to disease are unclear. This study aimed to estimate effects of body mass index (BMI) on an extensive set of circulating proteins.

METHODS

We used SomaLogic proteomic data from up to 2737 healthy participants from the INTERVAL study. Associations between self-reported BMI and 3622 unique plasma proteins were explored using linear regression. These were complemented by Mendelian randomisation (MR) analyses using a genetic risk score (GRS) comprised of 654 BMI-associated polymorphisms from a recent genome-wide association study (GWAS) of adult BMI. A disease enrichment analysis was performed using DAVID Bioinformatics 6.8 for proteins which were altered by BMI.

RESULTS

Observationally, BMI was associated with 1576 proteins (P < 1.4 × 10-5), with particularly strong evidence for a positive association with leptin and fatty acid-binding protein-4 (FABP4), and a negative association with sex hormone-binding globulin (SHBG). Observational estimates were likely confounded, but the GRS for BMI did not associate with measured confounders. MR analyses provided evidence for a causal relationship between BMI and eight proteins including leptin (0.63 standard deviation (SD) per SD BMI, 95% CI 0.48-0.79, P = 1.6 × 10-15), FABP4 (0.64 SD per SD BMI, 95% CI 0.46-0.83, P = 6.7 × 10-12) and SHBG (-0.45 SD per SD BMI, 95% CI -0.65 to -0.25, P = 1.4 × 10-5). There was agreement in the magnitude of observational and MR estimates (R2 = 0.33) and evidence that proteins most strongly altered by BMI were enriched for genes involved in cardiovascular disease.

CONCLUSIONS

This study provides evidence for a broad impact of adiposity on the human proteome. Proteins strongly altered by BMI include those involved in regulating appetite, sex hormones and inflammation; such proteins are also enriched for cardiovascular disease-related genes. Altogether, results help focus attention onto new proteomic signatures of obesity-related disease.