Genomic and drug target evaluation of 90 cardiovascular proteins in 30,931 individuals

Targeted proteomics identifies potential biomarkers of dysglycaemia, beta cell function and insulin sensitivity in Black African men and women

AIMS/HYPOTHESIS

Using a targeted proteomics approach, we aimed to identify and validate circulating proteins associated with impaired glucose metabolism (IGM) and type 2 diabetes in a Black South African cohort. In addition, we assessed sex-specific associations between the validated proteins and pathophysiological pathways of type 2 diabetes.

METHODS

This cross-sectional study included Black South African men (n=380) and women (n=375) who were part of the Middle-Aged Soweto Cohort (MASC). Dual-energy x-ray absorptiometry was used to determine fat mass and visceral adipose tissue, and fasting venous blood samples were collected for analysis of glucose, insulin and C-peptide and for targeted proteomics, measuring a total of 184 pre-selected protein biomarkers. An OGTT was performed on participants without diabetes, and peripheral insulin sensitivity (Matsuda index), HOMA-IR, basal insulin clearance, insulin secretion (C-peptide index) and beta cell function (disposition index) were estimated. Participants were classified as having normal glucose tolerance (NGT; n=546), IGM (n=116) or type 2 diabetes (n=93). Proteins associated with dysglycaemia (IGM or type 2 diabetes) in the MASC were validated in the Swedish EpiHealth cohort (NGT, n=1706; impaired fasting glucose, n=550; type 2 diabetes, n=210).

RESULTS

We identified 73 proteins associated with dysglycaemia in the MASC, of which 34 were validated in the EpiHealth cohort. Among these validated proteins, 11 were associated with various measures of insulin dynamics, with the largest number of proteins being associated with HOMA-IR. In sex-specific analyses, IGF-binding protein 2 (IGFBP2) was associated with lower HOMA-IR in women (coefficient -0.35; 95% CI -0.44, -0.25) and men (coefficient -0.09; 95% CI -0.15, -0.03). Metalloproteinase inhibitor 4 (TIMP4) was associated with higher insulin secretion (coefficient 0.05; 95% CI 0.001, 0.11; p for interaction=0.025) and beta cell function (coefficient 0.06; 95% CI 0.02, 0.09; p for interaction=0.013) in women only. In contrast, a stronger positive association between IGFBP2 and insulin sensitivity determined using an OGTT (coefficient 0.38; 95% CI 0.27, 0.49) was observed in men (p for interaction=0.004). A posteriori analysis showed that the associations between TIMP4 and insulin dynamics were not mediated by adiposity. In contrast, most of the associations between IGFBP2 and insulin dynamics, except for insulin secretion, were mediated by either fat mass index or visceral adipose tissue in men and women. Fat mass index was the strongest mediator between IGFBP2 and insulin sensitivity (total effect mediated 40.7%; 95% CI 37.0, 43.6) and IGFBP2 and HOMA-IR (total effect mediated 39.1%; 95% CI 31.1, 43.5) in men.

CONCLUSIONS/INTERPRETATION

We validated 34 proteins that were associated with type 2 diabetes, of which 11 were associated with measures of type 2 diabetes pathophysiology such as peripheral insulin sensitivity and beta cell function. This study highlights biomarkers that are similar between cohorts of different ancestry, with different lifestyles and sociodemographic profiles. The African-specific biomarkers identified require validation in African cohorts to identify risk markers and increase our understanding of the pathophysiology of type 2 diabetes in African populations.

Fibroblast Growth Factor-23 and Risk of Cardiovascular Diseases: A Mendelian Randomization Study

BACKGROUND

Fibroblast growth factor-23 (FGF-23) is associated with a range of cardiovascular and noncardiovascular diseases in conventional epidemiological studies, but substantial residual confounding may exist. Mendelian randomization approaches can help control for such confounding.

METHODS

SCALLOP Consortium data of 19,195 participants were used to generate an FGF-23 genetic score. Data from 337,448 UK Biobank participants were used to estimate associations between higher genetically predicted FGF-23 concentration and the odds of any atherosclerotic cardiovascular disease (n=26,266 events), nonatherosclerotic cardiovascular disease (n=12,652), and noncardiovascular diseases previously linked to FGF-23. Measurements of carotid intima-media thickness and left ventricular mass were available in a subset. Associations with cardiovascular outcomes were also tested in three large case-control consortia: CARDIOGRAMplusC4D (coronary artery disease, n=181,249 cases), MEGASTROKE (stroke, n=34,217), and HERMES (heart failure, n=47,309).

RESULTS

We identified 34 independent variants for circulating FGF-23, which formed a validated genetic score. There were no associations between genetically predicted FGF-23 and any of the cardiovascular or noncardiovascular outcomes. In UK Biobank, the odds ratio (OR) for any atherosclerotic cardiovascular disease per 1-SD higher genetically predicted logFGF-23 was 1.03 (95% confidence interval [95% CI], 0.98 to 1.08), and for any nonatherosclerotic cardiovascular disease, it was 1.01 (95% CI, 0.94 to 1.09). The ORs in the case-control consortia were 1.00 (95% CI, 0.97 to 1.03) for coronary artery disease, 1.01 (95% CI, 0.95 to 1.07) for stroke, and 1.00 (95% CI, 0.95 to 1.05) for heart failure. In those with imaging, logFGF-23 was not associated with carotid or cardiac abnormalities.

CONCLUSIONS

Genetically predicted FGF-23 levels are not associated with atherosclerotic and nonatherosclerotic cardiovascular diseases, suggesting no important causal link.

PODCAST

This article contains a podcast at https://dts.podtrac.com/redirect.mp3/www.asn-online.org/media/podcast/CJASN/2023_01_10_CJN05080422.mp3.

Ten interleukins and risk of prostate cancer

Background

Interleukins (ILs) have been reported to be related to prostate cancer. The aims of this study were to estimate the levels for several key interleukins in prostate cancer and the causal effects between them.

Methods

We conducted a bi-directional two-sample Mendelian randomization (MR) study to assess the causal associations between ILs and prostate cancer. Genetic instruments and summary-level data for 10 ILs were obtained from three genome-wide association meta-analyses. Prostate cancer related data were obtained from the PRACTICAL (79,148 cases and 61,106 controls), UK Biobank (7,691 cases and 169,762 controls) and FinnGen consortium (10,414 cases and 124,994 controls), respectively.

Results

The odds ratio of prostate cancer was 0.92 (95% confidence interval (CI), 0.89, 0.96; P=1.58×10^-05) and 1.12 (95% CI, 1.07, 1.17; P=6.61×10^-07) for one standard deviation increase in genetically predicted IL-1ra and IL-6 levels, respectively. Genetically predicted levels of IL-1ß, IL-2a, IL-6ra, IL-8, IL-16, IL-17, IL-18, and IL-27 were not associated with the risk of prostate cancer. Reverse MR analysis did not find the associations between genetic liability to prostate cancer and higher levels of IL-1ra (β, -0.005; 95% CI, -0.010, 0.001; P=0.111) and IL-6 (β, 0.002; 95% CI, -0.011, 0.014; P=0.755).

Conclusion

This MR study suggests that long-term IL-6 may increase the risk of prostate cancer and IL-1ra may reduce it.

Causal effects for genetic variants of osteoprotegerin on the risk of acute myocardial infarction and coronary heart disease: A two-sample Mendelian randomization study

Although since the 1980s, the mortality of coronary heart disease(CHD) has obviously decreased due to the rise of coronary intervention, the mortality and disability of CHD were still high in some countries. Etiological studies of acute myocardial infarction(AMI) and CHD were extremely important. In this study, we used two-sample Mendelian randomization(TSMR) method to collect GWAS statistics of osteoprotegerin (OPG), AMI and CHD to reveal the causal relationship between OPG and these two diseases. In total, we identified 7 genetic variants associated with AMI and 7 genetic variants associated with CHD that were not found to be in linkage disequilibrium(LD; r ² < 0.001). Evidence of a positive effect of an OPG genetic susceptibility on AMI was discovered(IVW OR = 0.877; 95% CI = 0.787-0.977; p = 0.017; 7 SNPs) and CHD (IVW OR = 0.892; 95% CI = 0.803-0.991; p = 0.033; 7 SNPs). After removing the influence of rs1385492, we found that there was a correlation between OPG and AMI/CHD (AMI: weighted median OR = 0.818;95% CI = 0.724-0.950; p = 0.001; 6SNPs;CHD: weighted median OR = 0.842; 95% CI = 0.755-0.938; p = 1.893 × 10^-3; 6SNPs). The findings of our study indicated that OPG had a tight genetic causation association with MI or CHD. This genetic causal relationship presented us with fresh ideas for the etiology of AMI and CHD, which is an area of research that will continue in the future.

Inflammatory bowel diseases, interleukin-6 and interleukin-6 receptor subunit alpha in causal association with cerebral cortical structure: a Mendelian randomization analysis

Background

Neurological involvement and psychiatric manifestations have been documented in clinical cases of inflammatory bowel disease (IBD); however, the presence of a causal relationship remains elusive. The objective of this study is to investigate the modifications occurring in the cerebral cortex as a result of IBD.

Methods

A compendium of data extracted from a genome-wide association study (GWAS) involving a maximum of 133,380 European subjects. A series of Mendelian random analyses were applied to exclude heterogeneity and pleiotropy, ensuring the stability of the results.

Results

Neither IBDs nor inflammatory cytokines (IL-6/IL-6Rα) were found to have a significant causality with surface area (SA) and thickness (TH) at the global level. At the regional functional brain level, Crohn's disease (CD) significantly decreased the TH of pars orbitalis (β=-0.003mm, Se=0.001mm, p_ivw =4.85×10^-4). IL-6 was observed to reduce the SA of middle temporal (β=-28.575mm², Se=6.482mm², p_ivw=1.04×10^-5) and increase the TH of fusiform (β=0.008mm, Se=0.002mm, p_ivw=8.86×10^-5) and pars opercularis (β=0.009mm, Se=0.002mm, p_ivw=2.34×10^-4). Furthermore, a causal relationship between IL-6Rα and an increase in the SA of superior frontal (β=21.132mm², Se=5.806mm², p_ivw=2.73×10^-4) and the TH of supramarginal (β=0.003mm, Se=0.0002mm, p_ivw=7.86×10^-37). All results passed sensitivity analysis and no heterogeneity and pleiotropy were detected.

Conclusion

The correlation between IBD and changes in cerebral cortical structures implies the existence of a gut-brain axis at the organismal level. It is recommended that clinical patients with IBD prioritize long-term management of inflammation, as changes at the organismal level can lead to functional pathologies. Magnetic resonance imaging (MRI) may be considered as an additional screening option for IBD.

Variants Tagging LGALS-3 Haplotype Block in Association with First Myocardial Infarction and Plasma Galectin-3 Six Months after the Acute Event

Galectin-3 is encoded by LGALS-3, located in a unique haplotype block in Caucasians. According to the Tagger server, rs4040064, rs11628437, and rs7159490 cover 82% (r² > 0.8) of the genetic variance of this HapBlock. Our aims were to examine the association of their haplotypes with first myocardial infarction (MI), changes in left ventricular echocardiographic parameters over time, and impact on plasma galectin-3 and LGALS-3 mRNA in peripheral blood mononuclear cells, both 6 months post-MI. The study group consisted of 546 MI patients and 323 controls. Gene expression was assessed in 92 patients and plasma galectin-3 in 189 patients. Rs4040064, rs11628437, rs7159490, and LGALS-3 mRNA expression were detected using TaqMan^® technology. Plasma galectin-3 concentrations were determined by the ELISA method. We found that the TGC haplotype could have a protective effect against MI (adjusted OR 0.19 [0.05-0.72], p = 0.015) and that the GAC haplotype had significantly higher galectin-3 concentrations (48.3 [37.3-59.4] ng/mL vs. 18.9 [14.5-23.4] ng/mL, p < 0.0001), both in males and compared to the referent haplotype GGC. Higher plasma Gal-3 was also associated with higher NYHA class and systolic dysfunction. Our results suggest that variants tagging LGALS-3 HapBlock could reflect plasma Gal-3 levels 6 months post-MI and may have a potential protective effect against MI in men. Further replication, validation, and functional studies are needed.

Bi-directional Mendelian randomization analysis provides evidence for the causal involvement of dysregulation of CXCL9, CCL11 and CASP8 in the pathogenesis of ulcerative colitis

BACKGROUND AND AIMS

Systemic inflammation is well-recognized to be associated with ulcerative colitis (UC), but whether these effects are causal or consequential remains unclear. We aimed to define potential causal relationship of cytokine dysregulation with different tiers of evidence.

METHODS

We firstly synthesized serum proteomic profiling data from two multi-centered observational studies, in which a panel of systemic inflammatory proteins was analyzed to examine their associations with UC risk. To further dissect observed associations, we then performed a bidirectional two-sample Mendelian randomization (TSMR) analysis from both forward and reverse directions using five genome-wide association study (GWAS) summary level data for serum proteomic profiles and the largest GWAS of 28,738 European-ancestry individuals for UC risk.

RESULTS

Pooled analysis of serum proteomic data identified 14 proteins to be associated with the risk of UC. Forward MR analysis using only cis-acting protein quantitative trait loci (cis-pQTLs) or trans-pQTLs further validated causal associations of two chemokines and the increased risk of UC: C-X-C motif chemokine ligand 9 (CXCL9) (OR, 1.45, 95% CI, 1.08-1.95, P=.012) and C-C motif chemokine ligand 11 (CCL11) (OR, 1.14, 95%CI: 1.09-1.18, P=3.89×10  -10). Using both cis- and trans-acting pQTLs, an association of caspase-8 (CASP8) (OR, 1.04, 95% CI, 1.03-1.05, P= 7.63×10  -19) was additionally identified. Reverse MR did not find any influence of genetic predisposition to UC on any of these three inflammation proteins.

CONCLUSIONS

Pre-existing elevated levels of CXCL9, CCL11 and CASP8 may play a role in the pathogenesis of UC.

Increased soluble urokinase plasminogen activator levels modulate monocyte function to promote atherosclerosis

People with kidney disease are disproportionately affected by atherosclerosis for unclear reasons. Soluble urokinase plasminogen activator receptor (suPAR) is an immune-derived mediator of kidney disease, levels of which are strongly associated with cardiovascular outcomes. We assessed suPAR's pathogenic involvement in atherosclerosis using epidemiologic, genetic, and experimental approaches. We found serum suPAR levels to be predictive of coronary artery calcification and cardiovascular events in 5,406 participants without known coronary disease. In a genome-wide association meta-analysis including over 25,000 individuals, we identified a missense variant in the plasminogen activator, urokinase receptor (PLAUR) gene (rs4760), confirmed experimentally to lead to higher suPAR levels. Mendelian randomization analysis in the UK Biobank using rs4760 indicated a causal association between genetically predicted suPAR levels and atherosclerotic phenotypes. In an experimental model of atherosclerosis, proprotein convertase subtilisin/kexin-9 (Pcsk9) transfection in mice overexpressing suPAR (suPARTg) led to substantially increased atherosclerotic plaques with necrotic cores and macrophage infiltration compared with those in WT mice, despite similar cholesterol levels. Prior to induction of atherosclerosis, aortas of suPARTg mice excreted higher levels of CCL2 and had higher monocyte counts compared with WT aortas. Aortic and circulating suPARTg monocytes exhibited a proinflammatory profile and enhanced chemotaxis. These findings characterize suPAR as a pathogenic factor for atherosclerosis acting at least partially through modulation of monocyte function.

Genetic influences of the effect of circulating inflammatory cytokines on osteoarthritis in humans

OBJECTIVE

The causal relationship between inflammatory cytokines and Osteoarthritis (OA) has not been well investigated. This study investigated the causal role of inflammatory cytokines in the risk of OA and total joint arthroplasty using the Mendelian randomization (MR) method.

METHOD

Single nucleotide polymorphisms (SNPs) robustly associated with inflammatory cytokines were used as instrumental variables. The inverse-variance weighted (IVW) method with false discovery rate (FDR) adjusted P-value (q-value) for multiple comparisons were used as the main MR method to estimate causal effects based on the summary-level data for OA (knee and hip OA, respectively) and total joint arthroplasty (TJA). Sensitivity analyses validated the robustness of the results and ensured the absence of heterogeneity and horizontal pleiotropy.

RESULTS

After FDR adjustment, macrophage colony-stimulating factor (MCSF) and vascular endothelial growth factor (VEGF) were identified as causally associated with knee OA (MCSF, odds ratio [OR]: 1.16, 95% confidence interval [CI]: 1.09-1.23, q = 5.05 × 10^-5; VEGF, OR: 1.09, 95% CI: 1.04-1.15, q = 0.011). We also observed that genetically predicted MCSF and VEGF were positively associated with the risk of TJA, and MCP3 was negatively associated with for the risk of TJA, although the effects seem fairly modest. Sensitivity analysis further excluded the influence of heterogeneity and horizontal pleiotropy.

CONCLUSIONS

Inflammatory cytokines, namely MCSF and VEGF, were causally associated with knee OA, which could enhance our understanding of inflammation in OA pathology.

Mendelian randomization and genetic colocalization infer the effects of the multi-tissue proteome on 211 complex disease-related phenotypes

BACKGROUND

Human proteins are widely used as drug targets. Integration of large-scale protein-level genome-wide association studies (GWAS) and disease-related GWAS has thus connected genetic variation to disease mechanisms via protein. Previous proteome-by-phenome-wide Mendelian randomization (MR) studies have been mainly focused on plasma proteomes. Previous MR studies using the brain proteome only reported protein effects on a set of pre-selected tissue-specific diseases. No studies, however, have used high-throughput proteomics from multiple tissues to perform MR on hundreds of phenotypes.

METHODS

Here, we performed MR and colocalization analysis using multi-tissue (cerebrospinal fluid (CSF), plasma, and brain from pre- and post-meta-analysis of several disease-focus cohorts including Alzheimer disease (AD)) protein quantitative trait loci (pQTLs) as instrumental variables to infer protein effects on 211 phenotypes, covering seven broad categories: biological traits, blood traits, cancer types, neurological diseases, other diseases, personality traits, and other risk factors. We first implemented these analyses with cis pQTLs, as cis pQTLs are known for being less prone to horizontal pleiotropy. Next, we included both cis and trans conditionally independent pQTLs that passed the genome-wide significance threshold keeping only variants associated with fewer than five proteins to minimize pleiotropic effects. We compared the tissue-specific protein effects on phenotypes across different categories. Finally, we integrated the MR-prioritized proteins with the druggable genome to identify new potential targets.

RESULTS

In the MR and colocalization analysis including study-wide significant cis pQTLs as instrumental variables, we identified 33 CSF, 13 plasma, and five brain proteins to be putative causal for 37, 18, and eight phenotypes, respectively. After expanding the instrumental variables by including genome-wide significant cis and trans pQTLs, we identified a total of 58 CSF, 32 plasma, and nine brain proteins associated with 58, 44, and 16 phenotypes, respectively. For those protein-phenotype associations that were found in more than one tissue, the directions of the associations for 13 (87%) pairs were consistent across tissues. As we were unable to use methods correcting for horizontal pleiotropy given most of the proteins were only associated with one valid instrumental variable after clumping, we found that the observations of protein-phenotype associations were consistent with a causal role or horizontal pleiotropy. Between 66.7 and 86.3% of the disease-causing proteins overlapped with the druggable genome. Finally, between one and three proteins, depending on the tissue, were connected with at least one drug compound for one phenotype from both DrugBank and ChEMBL databases.

CONCLUSIONS

Integrating multi-tissue pQTLs with MR and the druggable genome may open doors to pinpoint novel interventions for complex traits with no effective treatments, such as ovarian and lung cancers.

Quantifying the role of transcript levels in mediating DNA methylation effects on complex traits and diseases

High-dimensional omics datasets provide valuable resources to determine the causal role of molecular traits in mediating the path from genotype to phenotype. Making use of molecular quantitative trait loci (QTL) and genome-wide association study (GWAS) summary statistics, we propose a multivariable Mendelian randomization (MVMR) framework to quantify the proportion of the impact of the DNA methylome (DNAm) on complex traits that is propagated through the assayed transcriptome. Evaluating 50 complex traits, we find that on average at least 28.3% (95% CI: [26.9%-29.8%]) of DNAm-to-trait effects are mediated through (typically multiple) transcripts in the cis-region. Several regulatory mechanisms are hypothesized, including methylation of the promoter probe cg10385390 (chr1:8'022'505) increasing the risk for inflammatory bowel disease by reducing PARK7 expression. The proposed integrative framework can be extended to other omics layers to identify causal molecular chains, providing a powerful tool to map and interpret GWAS signals.

Causal relationship between gut microbes and cardiovascular protein expression

Evidence supports associations between gut microbiota and cardiovascular protein levels in plasma. However, it is unclear whether these associations reflect a causal relationship. To reveal the causal relationship between gut microbiota and cardiovascular protein levels in plasma, we estimated their causal effects using two-sample Mendelian randomization (MR) analysis. Sensitivity analysis was also performed to assess the robustness of our results. Genome-wide association study (GWAS) of microbiomes in the MiBioGen study included 211 bacterial taxa (18,473 individuals), and GWAS of 90 cardiovascular proteins included 30,931 individuals. There were 196 bacterial taxa from five levels available for analysis. The following 14 causal relationships were identified: phylum Euryarchaeota and carbohydrate antigen 125 (β = 0.289), order Bacillales and CSF-1 (β = -0.211), genus Dorea and HSP-27 (β = 0.465), phylum Actinobacteria and IL-8 (β = 0.274), order Enterobacteriales and KIM-1 (β = -0.499), class Actinobacteria, genus Bifidobacterium, phylum Actinobacteria and LEP (β = -0.219, β = -0.201, and β = -0.221), genus Methanobrevibacter and NT-proBNP (β = 0.371), family Peptostreptococcaceae and SRC (β = 0.191), order Verrucomicrobiales, phylum Verrucomicrobia and TNF-R2 (β = 0.251 and β = 0.233), family Veillonellaceae and t-PA (β = 0.271), and class Erysipelotrichia and VEGF-D (β = 0.390). Sensitivity analysis showed no evidence of pleiotropy or heterogeneity. The results of the reverse MR analysis showed no reverse causality for any of the 13 gut microbes and 11 cardiovascular proteins. Mendelian randomization estimates provide strong evidence for a causal effect of gut microbiota-mediated alterations on cardiovascular protein expression.

Genome-wide association and Mendelian randomization study of fibroblast growth factor 21 reveals causal associations with hyperlipidemia and possibly NASH

BACKGROUND

Fibroblast growth factor 21 (FGF21) is a hepatokine that produces metabolic benefits, such as improvements of lipid profile. We performed a genome-wide association study (GWAS) to identify genetic variants associated with circulating FGF21 and investigated the causal effects of FGF21 on pertinent outcomes using Mendelian randomization (MR).

METHODS

We conducted a GWAS testing ∼7.8 million DNA sequence variants with circulating FGF21 in a discovery cohort of 6259 Swedish adults with replication in 4483 Swedish women. We then performed two-sample MR analyses of genetically predicted circulating FGF21 in relation to alcohol and nutrient intake, cardiovascular and metabolic biomarkers and diseases, and liver function biomarkers using publicly available GWAS summary statistics data.

RESULTS

Our GWAS identified multiple single-nucleotide polymorphisms with genome-wide significant associations (P < 5 × 10^-8) with circulating FGF21 on chromosomes 2 and 19 in or near the GCKR and FGF21 genes, respectively. The strongest signal at the FGF21 locus (rs2548957, β = 0.181, P < 2.18 × 10^-42) displayed in two-sample MR analyses robust associations with lower alcohol intake, lower circulating low-density lipoprotein cholesterol, apolipoprotein B, C-reactive protein, gamma-glutamyl transferase, and galectin-3 concentrations, and higher circulating insulin-like growth factor-I and alkaline phosphatase concentrations after correcting for multiple testing (P < 0.0018) whereas associations with fat mass, type 2 diabetes, and cardiovascular disease were largely null.

CONCLUSIONS

We identified robust associations of certain genetic variants in or near the GCKR and FGF21 genes with circulating FGF21 concentrations. Furthermore, our results support a strong causal effect of FGF21 on improved lipid profile, reduced alcohol consumption and C-reactive protein concentrations, and liver function biomarkers including fibrosis. We found largely null or weak positive associations with fat mass, diabetes, and cardiovascular disease as well as higher insulin-like growth factor-I concentrations, which could indicate a compensatory increase to regulate the above FGF21 resistant states in humans.

ST2 and CSF-1 as potential druggable targets of inflammatory bowel diseases: Results from two-sample Mendelian randomization study

Novel druggable targets are warranted for inflammatory bowel disease (IBD) treatment. We aimed to identify novel circulating proteins with causal associations with the risk of IBDs and provide potential therapeutic targets for IBD treatment. We performed a two-sample Mendelian randomization (MR) study to explore the associations of 55 circulating biomarkers on the risk of IBD, Crohn's disease (CD), and ulcerative colitis (UC) by leveraging the summary statistics from large genomewide association studies and protein quantitative trait loci studies. The individual estimate was pooled together by meta-analyses to estimate the causal effects of each outcome. In univariable MR, we identified several circulating proteins showed potential correlation with IBD, UC, and CD. Of note, we observed that a genetically proxied increased level of suppression of tumorigenicity 2 (ST2) was associated with an elevated risk of IBD (odds ratios [ORs] 1.133, 95% confidence interval [CI] 1.091-1.176, p < 0.0001), CD (ORs 1.188, 95% CI 1.103-1.281, p < 0.0001), and UC cohorts (ORs 1.087, 95% CI 1.050-1.125, p < 0.0001). Additionally, we observed a consistent positive correlation between the level of CSF-1 and the increased risk of IBD in individual MR, with statistically significant causal associations in the meta-analyses with ORs equal to 1.217 (IBD, 95% CI 1.115-1.328, p < 0.0001), 1.223 (CD, 95% CI 1.082-1.382, p = 0.0013), and 1.179 (UC, 95% CI 1.055-1.317, p = 0.0037). This study provided evidence for potential casual associations between circulating ST2 and CSF-1 levels, and increased risks of IBD, UC, and CD, implicating potential treatment targets for IBD and subtypes.

Proteome-wide Mendelian randomization in global biobank meta-analysis reveals multi-ancestry drug targets for common diseases

Proteome-wide Mendelian randomization (MR) shows value in prioritizing drug targets in Europeans but with limited evidence in other ancestries. Here, we present a multi-ancestry proteome-wide MR analysis based on cross-population data from the Global Biobank Meta-analysis Initiative (GBMI). We estimated the putative causal effects of 1,545 proteins on eight diseases in African (32,658) and European (1,219,993) ancestries and identified 45 and 7 protein-disease pairs with MR and genetic colocalization evidence in the two ancestries, respectively. A multi-ancestry MR comparison identified two protein-disease pairs with MR evidence in both ancestries and seven pairs with specific effects in the two ancestries separately. Integrating these MR signals with clinical trial evidence, we prioritized 16 pairs for investigation in future drug trials. Our results highlight the value of proteome-wide MR in informing the generalizability of drug targets for disease prevention across ancestries and illustrate the value of meta-analysis of biobanks in drug development.

Causal associations between CD40/CD40L and aortic diseases: A mendelian randomization study

Background: CD40 and CD40L have been reported as associated with aortic dissection (AD) and aortic aneurysm (AA), but the causality of the associations has not been established yet. Methods: We conducted a two-sample Mendelian randomization (MR) study to assess the causal inference between CD40/CD40L and aortic diseases including AD and AA. The instrumental variables (IVs) for CD40 and CD40L were selected from a high-quality protein quantitative trait loci dataset released by a genomic study involving 30,931 individuals of European ancestry. The genome-wide association studies summary statistics for AD and AA were from the FinnGen Release 7, with 288638 controls for all outcomes of interests, 680 cases for AD and 6,092 cases for AA, also from European ancestry. For AA subtypes, there were 5,881 cases of thoracic AA (TAA) and 2,434 cases of abdominal AA (AAA) respectively. Inverse-variance weighted and Wald ratio were applied for calculating causal estimates. Horizontal pleiotropy and heterogeneity were assessed using MR-Egger regression analysis and Cochran Q test, respectively. Leave-one-out analyses were further performed. Results: Three single-nucleotide polymorphisms (SNPs) for CD40 and one SNP for CD40L were selected as IVs. We found genetic proxied CD40 levels inversely associated with the risk of AD (odds ratio [OR]: 0.777, 95% confidence interval [CI]: 0.618-0.978, p = 0.031) and AA (OR: 0.905, 95% CI: 0.837-0.978, p = 0.012), consistent across TAA (both p < 0.050). There were trends of increased risks of AD and AA in the presence of CD40L while not reaching statistical significance. No significant horizontal pleiotropy or heterogeneity was observed. Conclusion: Our MR study provides evidence supporting the causal association between CD40 and the reduced risks of both AD and AA.

Repurposing drugs to treat cardiovascular disease in the era of precision medicine

Drug repurposing is the use of a given therapeutic agent for indications other than that for which it was originally designed or intended. The concept is appealing because of potentially lower development costs and shorter timelines than are needed to produce a new drug. To date, drug repurposing for cardiovascular indications has been opportunistic and driven by knowledge of disease mechanisms or serendipitous observation rather than by systematic endeavours to match an existing drug to a new indication. Innovations in two areas of personalized medicine - computational approaches to associate drug effects with disease signatures and predictive model systems to screen drugs for disease-modifying activities - support efforts that together create an efficient pipeline to systematically repurpose drugs to treat cardiovascular disease. Furthermore, new experimental strategies that guide the medicinal chemistry re-engineering of drugs could improve repurposing efforts by tailoring a medicine to its new indication. In this Review, we summarize the historical approach to repurposing and discuss the technological advances that have created a new landscape of opportunities.

Identification of 969 protein quantitative trait loci in an African American population with kidney disease attributed to hypertension

Investigations into the causal underpinnings of disease processes can be aided by the incorporation of genetic information. Genetic studies require populations varied in both ancestry and prevalent disease in order to optimize discovery and ensure generalizability of findings to the global population. Here, we report the genetic determinants of the serum proteome in 466 African Americans with chronic kidney disease attributed to hypertension from the richly phenotyped African American Study of Kidney Disease and Hypertension (AASK) study. Using the largest aptamer-based protein profiling platform to date (6,790 proteins or protein complexes), we identified 969 genetic associations with 900 unique proteins; including 52 novel cis (local) associations and 379 novel trans (distant) associations. The genetic effects of previously published cis-protein quantitative trait loci (pQTLs) were found to be highly reproducible, and we found evidence that our novel genetic signals colocalize with gene expression and disease processes. Many trans- pQTLs were found to reflect associations mediated by the circulating cis protein, and the common trans-pQTLs are enriched for processes involving extracellular vesicles, highlighting a plausible mechanism for distal regulation of the levels of secreted proteins. Thus, our study generates a valuable resource of genetic associations linking variants to protein levels and disease in an understudied patient population to inform future studies of drug targets and physiology.

Causal risk factors for asthma in Mendelian randomization studies: A systematic review and meta-analysis

BACKGROUND

Several risk factors for asthma have been proposed; however, the causality of these associations is sometimes unclear. Mendelian randomization is a powerful epidemiological approach that can help elucidate the causality of risk factors. The aim of the present study was to identify causal risk factors for asthma through Mendelian Randomization studies.

METHODS

A systematic search of PubMed and EMBASE was conducted, to identify studies investigating risk factors for asthma or respiratory allergies through Mendelian Randomization. When two or more studies investigated the same risk factor a meta-analysis was conducted. Of 239 studies initially identified, 41 were included.

RESULTS

A causal association between adiposity and adult asthma risk was found in 10 out of 12 studies with a summary risk ratio of 1.05 per kg/m2 increase in BMI (95% CI: 1.03-1.07). Puberty timing (n = 3), alcohol (n = 2), and linoleic acid (n = 1) had causal effects on asthma risk, while vitamins/minerals (n = 6) showed no consistent effect on asthma. The effect of smoking on adult asthma conflicted between studies. Several of the significant associations of asthma with immune related proteins (n = 5) and depression (n = 2) investigated through multiple traits analyses could generally benefit from replications in independent datasets.

CONCLUSION

This systematic review and meta-analysis found evidence for causal effects of adiposity, puberty timing, linoleic acid, alcohol, immune related proteins, and depression on risk of asthma.

Role of Inflammatory Factors in Mediating the Effect of Lipids on Nonalcoholic Fatty Liver Disease: A Two-Step, Multivariable Mendelian Randomization Study

AIMS/HYPOTHESIS

20-80% of Nonalcoholic Fatty Liver Disease (NAFLD) have been observed to have dyslipidemia. Nevertheless, the probable mechanism of dyslipidemia's effect on NAFLD remains unclear. Mendelian randomization (MR) was utilized to investigate the relationship between lipids, inflammatory factors, and NAFLD; and also, to determine the proportion mediated by interleukin-17(IL-17) and interleukin-1β(IL-1β) for the effect between lipids and NAFLD.

METHODS

Summary statistics of traits were obtained from the latest and largest genome-wide association study (GWAS). The UK Biobank provided a summary of lipid statistics, which comprised up to 500,000 participants of European descent. And NAFLD GWAS summary statistics were obtained from the FinnGen Biobank which included a total sample size of 218,792 participants of European ancestry. In order to gain an overall picture of how lipids affect NAFLD, MR with two samples was carried out. Multivariable MR determined lipids direct effects on NAFLD after adjusting for inflammatory factors, namely IL-1β, interleukin-6(IL-6), interleukin-16(IL-16), IL-17, and interleukin-18(IL-18); those lipids comprise HDL cholesterol (HDL-C), LDL cholesterol (LDL-C), triglycerides (TGs), apolipoprotein A1 (ApoA1), and apolipoprotein B (ApoB). For the purpose of determining the MR impact, an inverse variance weighted (IVW) meta-analysis of each Wald Ratio was carried out, while other methods were also performed for sensitivity analysis.

RESULTS

We discovered a positive association between genetically predicted TGs levels and a 45.5% elevated risk of NAFLD, while genetically predicted IL-1β [(IVW: OR 1.315 (1.060-1.630), p = 0.012) and IL-17 [(IVW: OR 1.468 (1.035-2.082), p = 0.032] were positively associated with 31.5% and 46.8% increased risk of NAFLD, respectively. Moreover, TG was positively associated with 10.5% increased risk of IL-1β and 17.3% increased risk of IL-17. The proportion mediated by IL-17 and IL-1β respectively and both was 2.6%, 3.1%, 14.1%.

CONCLUSION

Genetically predicted TGs, IL-1β, and IL-17 were positively associated with increased risk of NAFLD, with evidence that IL-1β and IL-17 mediated TGs effect upon NAFLD risk. It indicated that early diet management, weight management, lipid-lowering and anti-inflammatory treatment should be carried out for patients with hyperlipidemia to prevent the NAFLD.

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.

Constructing an atlas of associations between polygenic scores from across the human phenome and circulating metabolic biomarkers

BACKGROUND

Polygenic scores (PGS) are becoming an increasingly popular approach to predict complex disease risk, although they also hold the potential to develop insight into the molecular profiles of patients with an elevated genetic predisposition to disease.

METHODS

We sought to construct an atlas of associations between 125 different PGS derived using results from genome-wide association studies and 249 circulating metabolites in up to 83,004 participants from the UK Biobank.

RESULTS

As an exemplar to demonstrate the value of this atlas, we conducted a hypothesis-free evaluation of all associations with glycoprotein acetyls (GlycA), an inflammatory biomarker. Using bidirectional Mendelian randomization, we find that the associations highlighted likely reflect the effect of risk factors, such as adiposity or liability towards smoking, on systemic inflammation as opposed to the converse direction. Moreover, we repeated all analyses in our atlas within age strata to investigate potential sources of collider bias, such as medication usage. This was exemplified by comparing associations between lipoprotein lipid profiles and the coronary artery disease PGS in the youngest and oldest age strata, which had differing proportions of individuals undergoing statin therapy. Lastly, we generated all PGS-metabolite associations stratified by sex and separately after excluding 13 established lipid-associated loci to further evaluate the robustness of findings.

CONCLUSIONS

We envisage that the atlas of results constructed in our study will motivate future hypothesis generation and help prioritize and deprioritize circulating metabolic traits for in-depth investigations. All results can be visualized and downloaded at http://mrcieu.mrsoftware.org/metabolites_PRS_atlas.

FUNDING

This work is supported by funding from the Wellcome Trust, the British Heart Foundation, and the Medical Research Council Integrative Epidemiology Unit.

Constructing an atlas of associations between polygenic scores from across the human phenome and circulating metabolic biomarkers

Background

Polygenic scores (PGS) are becoming an increasingly popular approach to predict complex disease risk, although they also hold the potential to develop insight into the molecular profiles of patients with an elevated genetic predisposition to disease.

Methods

We sought to construct an atlas of associations between 125 different PGS derived using results from genome-wide association studies and 249 circulating metabolites in up to 83,004 participants from the UK Biobank.

Results

As an exemplar to demonstrate the value of this atlas, we conducted a hypothesis-free evaluation of all associations with glycoprotein acetyls (GlycA), an inflammatory biomarker. Using bidirectional Mendelian randomization, we find that the associations highlighted likely reflect the effect of risk factors, such as adiposity or liability towards smoking, on systemic inflammation as opposed to the converse direction. Moreover, we repeated all analyses in our atlas within age strata to investigate potential sources of collider bias, such as medication usage. This was exemplified by comparing associations between lipoprotein lipid profiles and the coronary artery disease PGS in the youngest and oldest age strata, which had differing proportions of individuals undergoing statin therapy. Lastly, we generated all PGS-metabolite associations stratified by sex and separately after excluding 13 established lipid-associated loci to further evaluate the robustness of findings.

Conclusions

We envisage that the atlas of results constructed in our study will motivate future hypothesis generation and help prioritize and deprioritize circulating metabolic traits for in-depth investigations. All results can be visualized and downloaded at http://mrcieu.mrsoftware.org/metabolites_PRS_atlas.

Funding

This work is supported by funding from the Wellcome Trust, the British Heart Foundation, and the Medical Research Council Integrative Epidemiology Unit.

Genetic Liability to Rheumatoid Arthritis in Relation to Coronary Artery Disease and Stroke Risk

OBJECTIVE

To assess the causality of the associations of rheumatoid arthritis (RA) with coronary artery disease (CAD) and stroke using the Mendelian randomization approach.

METHODS

Independent single-nucleotide polymorphisms strongly associated with RA (n = 70) were selected as instrumental variables from a genome-wide association meta-analysis including 14,361 RA patients and 43,923 controls of European ancestry. Summary-level data for CAD, all stroke, any ischemic stroke and its subtypes, intracerebral hemorrhage (ICH), and subarachnoid hemorrhage were obtained from meta-analyses of genetic studies, international genetic consortia, the UK Biobank, and the FinnGen consortium. We obtained summary-level data for common cardiovascular risk factors and related inflammatory biomarkers to assess possible mechanisms.

RESULTS

Genetic liability to RA was associated with an increased risk of CAD and ICH. For a 1-unit increase in log odds of RA, the combined odds ratios were 1.02 (95% confidence interval [1.01, 1.03]; P = 0.003) for CAD and 1.05 (95% confidence interval [1.02, 1.08]; P = 0.001) for ICH. Genetic liability to RA was associated with increased levels of tumor necrosis factor and C-reactive protein (CRP). The association with CAD was attenuated after adjustment for genetically predicted CRP levels. There were no associations of genetic liability to RA with the other studied outcomes.

CONCLUSION

This study found that genetic liability to RA was associated with an increased risk of CAD and ICH and that the association with CAD might be mediated by CRP. The heightened cardiovascular risk should be actively monitored and managed in RA patients, and this may include dampening systemic inflammation.

Immune system-wide Mendelian randomization and triangulation analyses support autoimmunity as a modifiable component in dementia-causing diseases

Immune system and blood-brain barrier dysfunction are implicated in the development of Alzheimer's and other dementia-causing diseases, but their causal role remains unknown. We performed Mendelian randomization for 1,827 immune system- and blood-brain barrier-related biomarkers and identified 127 potential causal risk factors for dementia-causing diseases. Pathway analyses linked these biomarkers to amyloid-β, tau and α-synuclein pathways and to autoimmunity-related processes. A phenome-wide analysis using Mendelian randomization-based polygenic risk score in the FinnGen study (n = 339,233) for the biomarkers indicated shared genetic background for dementias and autoimmune diseases. This association was further supported by human leukocyte antigen analyses. In inverse-probability-weighted analyses that simulate randomized controlled drug trials in observational data, anti-inflammatory methotrexate treatment reduced the incidence of Alzheimer's disease in high-risk individuals (hazard ratio compared with no treatment, 0.64, 95% confidence interval 0.49-0.88, P = 0.005). These converging results from different lines of human research suggest that autoimmunity is a modifiable component in dementia-causing diseases.

Coupled myovascular expansion directs cardiac growth and regeneration

Heart regeneration requires multiple cell types to enable cardiomyocyte (CM) proliferation. How these cells interact to create growth niches is unclear. Here, we profile proliferation kinetics of cardiac endothelial cells (CECs) and CMs in the neonatal mouse heart and find that they are spatiotemporally coupled. We show that coupled myovascular expansion during cardiac growth or regeneration is dependent upon VEGF-VEGFR2 signaling, as genetic deletion of Vegfr2 from CECs or inhibition of VEGFA abrogates both CEC and CM proliferation. Repair of cryoinjury displays poor spatial coupling of CEC and CM proliferation. Boosting CEC density after cryoinjury with virus encoding Vegfa enhances regeneration. Using Mendelian randomization, we demonstrate that circulating VEGFA levels are positively linked with human myocardial mass, suggesting that Vegfa can stimulate human cardiac growth. Our work demonstrates the importance of coupled CEC and CM expansion and reveals a myovascular niche that may be therapeutically targeted for heart regeneration.

Serum IL-1ra Is Associated with but Has No Genetic Link to Type 1 Diabetes

Interleukin-1 antagonism is a proposed biomarker and potential therapy for the delay and/or treatment of type 1 diabetes (T1D). We evaluated the role of circulating interleukin-1 receptor antagonist (IL-1ra) in a prospectively monitored cohort of T1D patients. In order to determine a mechanistic association between IL-1ra and T1D, we performed co-localization analyses between serum IL-1ra protein quantitative trait loci and T1D genome-wide analysis studies. Adjusting for human leukocyte antigen (HLA) genotypes, first degree relative status, gender, and age, serum levels of IL-1ra were lower in subjects who progressed to T1D compared to the controls (p = 0.023). Our results suggest that females have higher levels of IL-1ra compared to males (p = 0.005). The 2q14.1 region associated with serum IL-1ra levels is not associated with a risk of developing T1D. Our data suggest that IL-1 antagonism by IL-1ra is not an effective therapy in T1D, but IL-1ra may be a biomarker for progression to T1D.

Circulating vascular endothelial growth factor and cancer risk: A bidirectional mendelian randomization

In observational studies, circulating vascular endothelial growth factor (VEGF) has been reported to be associated with certain types of cancer. The purpose of this study was to verify whether there is a causal relationship between circulating VEGF and different types of cancer and the direction of the causal relationship. Summary statistical data were obtained from the corresponding genome-wide association studies (GWASs) to investigate the causal relationship between circulating VEGF and the risk of several cancers, including breast cancer, ovarian cancer, lung cancer, colorectal cancer, anus and anal canal cancer, prostate cancer, esophageal cancer, kidney cancer, bladder cancer, thyroid cancer, malignant neoplasm of the brain and malignant neoplasm of the liver and intrahepatic bile ducts. A two-sample bidirectional Mendelian randomization (MR) analysis and sensitivity tests were used to evaluate the validity of causality. A causal relationship was detected between circulating VEGF and colorectal cancer (OR 1.21, 95% CI 1.11–1.32, p < 0.000) and colon adenocarcinoma (OR 1.245, 95% CI 1.10–1.412, p < 0.000). Suggestive evidence of association was detected in VEGF on malignant neoplasms of the rectum (OR 1.16, 95% CI 1.00–1.34, p = 0.049). No causal relationship was found between circulating VEGF and other types of cancer, nor was there a reverse causal relationship from tumors to VEGF (p > 0.05). Circulating VEGF has a causal relationship with specific types of cancer. Our findings highlight and confirm the importance of circulating VEGF in the prevention and treatment of colorectal cancer.

Leukotriene A4 Hydrolase and Hepatocyte Growth Factor Are Risk Factors of Sudden Cardiac Death Due to First-Ever Myocardial Infarction

Patients at a high risk for sudden cardiac death (SCD) without previous history of cardiovascular disease remain a challenge to identify. Atherosclerosis and prothrombotic states involve inflammation and non-cardiac tissue damage that may play active roles in SCD development. Therefore, we hypothesized that circulating proteins implicated in inflammation and tissue damage are linked to the future risk of SCD. We conducted a prospective nested case–control study of SCD cases with verified myocardial infarction (N = 224) and matched controls without myocardial infarction (N = 224), aged 60 ± 10 years time and median time to event was 8 years. Protein concentrations (N = 122) were measured using a proximity extension immunoassay. The analyses revealed 14 proteins significantly associated with an increased risk of SCD, from which two remained significant after adjusting for smoking status, systolic blood pressure, BMI, cholesterol, and glucose levels. We identified leukotriene A4 hydrolase (LTA4H, odds ratio 1.80, corrected confidence interval (CIcorr) 1.02–3.17) and hepatocyte growth factor (HGF; odds ratio 1.81, CIcorr 1.06–3.11) as independent risk markers of SCD. Elevated LTA4H may reflect increased systemic and pulmonary neutrophilic inflammatory processes that can contribute to atherosclerotic plaque instability. Increased HGF levels are linked to obesity-related metabolic disturbances that are more prevalent in SCD cases than the controls.

Considerations of biomarker application for cancer continuum in the era of precision medicine

Purpose of the review

The goal of this review is to highlight emerging biomarker research by the key phases of the cancer continuum and outline the methodological considerations for biomarker application.

Recent findings

While biomarkers have an established role in targeted therapy and to some extent, disease monitoring, their role in early detection and survivorship remains to be elucidated. With the advent of omics technology, the discovery of biomarkers has been accelerated exponentially, therefore careful consideration to ensure an unbiased study design and robust validity is crucial.

Summary

The rigor of biomarker research holds the key to the success of precision health care. The potential clinical utility and the feasibility of implementation should be central to future biomarker research study design.

Circulating adipokine levels and preeclampsia: A bidirectional Mendelian randomization study

Background: Several observational studies have demonstrated that significantly rising circulating adipokine levels are pervasive in preeclampsia or eclampsia disorder (or preeclampsia toxemia (PET)). However, it remains unclear whether this relationship is causal. In this study, we sought to elucidate the causal effects of circulating adipokine levels on PET.

Methods: Summary-level data and independent genetic variants strongly associated with common adipokine molecule (adiponectin, leptin, resistin, sOB-R, and PAI-1) levels were drawn from public genome-wide association study (GWASs). Additionally, the corresponding effects between instrumental variables and PET outcomes were acquired from the FinnGen consortium, including 4,743 cases and 136,325 controls of European ancestry. Subsequently, an inverse-variance weighted (IVW) approach was applied for the principal two-sample Mendelian randomization (MR) and multivariable MR (MVMR) analyses. Various complementary sensitivity analyses were then carried out to determine the robustness of our models.

Results: The results of the IVW method did not reveal any causal relationship shared across genetically predisposed adipokine levels and PET risk (for adiponectin, OR = 0.86, 95% CI: 0.65–1.13, p = 0.274). Additionally, no significant associations were identified after taking into account five circulating adipokines in MVMR research. Complementary sensitivity analysis also supported no significant associations between them. In the reverse MR analysis, genetically predicted PET risk showed a suggestive association with elevating PAI-1 levels by the IVW method (Beta = 0.120, 95% CI: 0.014, 0.227, p = 0.026). Furthermore, there were no strong correlations between genetic liability to PET and other adipokine levels (p > 0.05).

Conclusion: Our MR study did not provide robust evidence supporting the causal role of common circulating adipokine levels in PET, whereas genetically predicted PET may instrumentally affect PAI-1 levels. These findings suggest that PAI-1 may be a useful biomarker for monitoring the diagnosis or therapy of PET rather than a therapeutic target for PET.

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.

Integrated analyses of growth differentiation factor-15 concentration and cardiometabolic diseases in humans

Growth differentiation factor-15 (GDF15) is a stress response cytokine that is elevated in several cardiometabolic diseases and has attracted interest as a potential therapeutic target. To further explore the association of GDF15 with human disease, we conducted a broad study into the phenotypic and genetic correlates of GDF15 concentration in up to 14,099 individuals. Assessment of 772 traits across 6610 participants in FINRISK identified associations of GDF15 concentration with a range of phenotypes including all-cause mortality, cardiometabolic disease, respiratory diseases and psychiatric disorders, as well as inflammatory markers. A meta-analysis of genome-wide association studies (GWAS) of GDF15 concentration across three different assay platforms (n=14,099) confirmed significant heterogeneity due to a common missense variant (rs1058587; p.H202D) in GDF15, potentially due to epitope-binding artefacts. After conditioning on rs1058587, statistical fine mapping identified four independent putative causal signals at the locus. Mendelian randomisation (MR) analysis found evidence of a causal relationship between GDF15 concentration and high-density lipoprotein (HDL) but not body mass index (BMI). Using reverse MR, we identified a potential causal association of BMI on GDF15 (IVW pFDR = 0.0040). Taken together, our data derived from human population cohorts do not support a role for moderately elevated GDF15 concentrations as a causal factor in human cardiometabolic disease but support its role as a biomarker of metabolic stress.

Genetic determinants of metabolic biomarkers and their associations with cardiometabolic traits in Hispanic/Latino adolescents

BACKGROUND

Metabolic regulation plays a significant role in energy homeostasis, and adolescence is a crucial life stage for the development of cardiometabolic disease (CMD). This study aims to investigate the genetic determinants of metabolic biomarkers-adiponectin, leptin, ghrelin, and orexin-and their associations with CMD risk factors.

METHODS

We characterized the genetic determinants of the biomarkers among Hispanic/Latino adolescents of the Santiago Longitudinal Study (SLS) and identified the cumulative effects of genetic variants on adiponectin and leptin using biomarker polygenic risk scores (PRS). We further investigated the direct and indirect effect of the biomarker PRS on downstream body fat percent (BF%) and glycemic traits using structural equation modeling.

RESULTS

We identified putatively novel genetic variants associated with the metabolic biomarkers. A substantial amount of biomarker variance was explained by SLS-specific PRS, and the prediction was improved by including the putatively novel loci. Fasting blood insulin and insulin resistance were associated with PRS for adiponectin, leptin, and ghrelin, and BF% was associated with PRS for adiponectin and leptin. We found evidence of substantial mediation of these associations by the biomarker levels.

CONCLUSIONS

The genetic underpinnings of metabolic biomarkers can affect the early development of CMD, partly mediated by the biomarkers.

IMPACT

This study characterized the genetic underpinnings of four metabolic hormones and investigated their potential influence on adiposity and insulin biology among Hispanic/Latino adolescents. Fasting blood insulin and insulin resistance were associated with polygenic risk score (PRS) for adiponectin, leptin, and ghrelin, with evidence of some degree of mediation by the biomarker levels. Body fat percent (BF%) was also associated with PRS for adiponectin and leptin. This provides important insight on biological mechanisms underlying early metabolic dysfunction and reveals candidates for prevention efforts. Our findings also highlight the importance of ancestrally diverse populations to facilitate valid studies of the genetic architecture of metabolic biomarker levels.

Causal analysis between altered levels of interleukins and obstructive sleep apnea

Background

Inflammation proteins including interleukins (ILs) have been reported to be related to obstructive sleep apnea (OSA). The aims of this study were to estimate the levels for several key interleukins in OSA and the causal effects between them.

Method

Weighted mean difference (WMD) was used to compare the expression differences of interleukins between OSA and control, and the changed levels during OSA treatments in the meta-analysis section. A two-sample Mendelian randomization (MR) was used to estimate the causal directions and effect sizes between OSA risks and interleukins. The inverse-variance weighting (IVW) was used as the primary method followed by several other MR methods including MR Egger, Weighted median, and MR-Robust Adjusted Profile Score as sensitivity analysis.

Results

Nine different interleukins—IL-1β, IL-2, IL-4, IL-6, IL-8, IL-12, IL-17, IL-18, and IL-23—were elevated in OSA compared with control to varying degrees, ranging from 0.82 to 100.14 pg/ml, and one interleukin, IL-10, was decreased by 0.77 pg/ml. Increased IL-1β, IL-6, and IL-8 rather than IL-10 can be reduced in OSA by effective treatments. Further, the MR analysis of the IVW method showed that there was no significant evidence to support the causal relationships between OSA and the nine interleukins—IL-1β, IL-2, IL-4, IL-5, IL-6, IL-8, IL-10, IL-17, and IL-18. Among them, the causal effect of OSA on IL-5 was almost significant [estimate: 0.267 (−0.030, 0.564), p = 0.078]. These results were consistent in the sensitivity analysis.

Conclusions

Although IL-1β, IL-2, IL-4, IL-6, IL-8, IL-12, IL-17, IL-18, and IL-23 were increasing and IL-10 was reducing in OSA, no significant causal relationships were observed between them by MR analysis. Further research is needed to test the causality of OSA risk on elevated IL-5 level.

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.

Mendelian randomization study for the roles of IL-18 and IL-1 receptor antagonist in the development of inflammatory bowel disease

BACKGROUND AND AIMS

IL-1 and IL-18 play important roles in intestine barrier integrity maintenance and inflammatory response. However, their net effects on the risk of IBD are still inconclusive. Here, we used Mendelian randomization (MR) approaches to investigate the causal associations of IL-18 and IL-1Ra (receptor antagonist) on the risks of IBD and subtypes.

METHODS

For IL-18, both three-sample and two-sample MR approaches were used for the causal inferences. In three-sample MR, three single nucleotide polymorphisms (SNPs) and the effect values were extracted from two quantitative trait loci (pQTL) datasets with non-overlapping populations. In two-sample MR, we extracted genetic instruments information from the same larger pQTL dataset. For IL-1Ra, we applied the two-sample MR method with summary-statistics from the larger pQTL dataset. Summary-level results of three large IBD/CD/UC genome-wide association studies in European ancestry were employed. Inverse-variance weighted method, various sensitivity analyses and meta-analysis were performed to give causal estimates, detect heterogeneity and correct for outliers.

RESULTS

We observed consistent positive causal effects of IL-18 on all three major outcomes using three-sample MR, with meta-analyses odds ratios (ORs) equal to 1.240 (IBD), 1.199 (CD) and 1.274 (UC) respectively. The two-sample MR demonstrated similar results. Moreover, genetically predicted IL-1Ra is inversely associated with the risk of IBD/UC/CD with ORs equal to 0.915 (IBD), 0.902 (CD) and 0.899 (UC) respectively in meta-analyses.

CONCLUSIONS

This study suggested genetically predicted IL-18 and IL-1Ra level are causally associated with an increased and decreased risk of IBD and subtypes.

Circulating Hepatocyte Growth Factor Reflects Activation of Vascular Repair in Response to Stress

•

HGF is released by stressed human vascular cells and promotes vascular cell repair responses in autocrine and/or paracrine ways.

•

Subjects with a low capacity to express HGF in response to systemic stress have an increased cardiovascular risk.

•

Human atherosclerotic plaques with a low content of HGF have a more unstable phenotype.

•

The present study shows that subjects with a low ability to express HGF in response to metabolic stress have an increased risk to suffer cardiovascular events.

Hepatocyte growth factor (HGF) is released by stressed human vascular cells and promotes vascular cell repair responses in both autocrine and paracrine ways. Subjects with a low capacity to express HGF in response to systemic stress have an increased cardiovascular risk. Human atherosclerotic plaques with a low content of HGF have a more unstable phenotype. The present study shows that subjects with a low ability to express HGF in response to metabolic stress have an increased risk to suffer myocardial infarction and stroke.

Plaque characteristics and biomarkers predicting regression and progression of carotid atherosclerosis

The factors that influence the atherosclerotic disease process in high-risk individuals remain poorly understood. Here, we used a combination of vascular imaging, risk factor assessment, and biomarkers to identify factors associated with 3-year change in carotid disease severity in a cohort of high-risk subjects treated with preventive therapy (n = 865). The results show that changes in intima-media thickness (IMT) are most pronounced in the carotid bulb. Progression of bulb IMT demonstrates independent associations with baseline bulb IMT, the plaque gray scale median (GSM), and the plasma level of platelet-derived growth factor (PDGF) (standardized β-coefficients and 95% confidence interval [CI] −0.14 [−0.06 to −0.02] p = 0.001, 0.15 [0.02–0.07] p = 0.001, and 0.20 [0.03–0.07] p < 0.001, respectively). Plasma PDGF correlates with the plaque GSM (0.23 [0.15–0.29] p < 0.001). These observations provide insight into the atherosclerotic process in high-risk subjects by showing that progression primarily occurs in fibrotic plaques and is associated with increased levels of PDGF.

•

High age, male gender, and smoking increases risk of carotid disease progression

•

Plaques that progress are more echogenic, indicating an increased degree of fibrosis

•

Progression is associated with high plasma levels of pro-fibrotic growth factors

•

Regression is most common in large, less fibrotic plaques

Gene-based whole genome sequencing meta-analysis of 250 circulating proteins in three isolated European populations

OBJECTIVE

Deep sequencing offers unparalleled access to rare variants in human populations. Understanding their role in disease is a priority, yet prohibitive sequencing costs mean that many cohorts lack the sample size to discover these effects on their own. Meta-analysis of individual variant scores allows the combination of rare variants across cohorts and study of their aggregated effect at the gene level, boosting discovery power. However, the methods involved have largely not been field-tested. In this study, we aim to perform the first meta-analysis of gene-based rare variant aggregation optimal tests, applied to the human cardiometabolic proteome.

METHODS

Here, we carry out this analysis across MANOLIS, Pomak and ORCADES, three isolated European cohorts with whole-genome sequencing (total N = 4,422). We examine the genetic architecture of 250 proteomic traits of cardiometabolic relevance. We use a containerised pipeline to harmonise variant lists across cohorts and define four sets of qualifying variants. For every gene, we interrogate protein-damaging variants, exonic variants, exonic and regulatory variants, and regulatory only variants, using the CADD and Eigen scores to weigh variants according to their predicted functional consequence. We perform single-cohort rare variant analysis and meta-analyse variant scores using the SMMAT package.

RESULTS

We describe 5 rare variant pQTLs (RV-pQTL) which pass our stringent significance threshold (7.45 × 10-11) and quality control procedure. These were split between four cis signals for MARCO, TEK, MMP2 and MPO, and one trans association for GDF2 in the SERPINA11 gene. We show that the cis-MPO association, which was not detectable using the single-point data alone, is driven by 5 missense and frameshift variants. These include rs140636390 and rs119468010, which are specific to MANOLIS and ORCADES, respectively. We show how this kind of signal could improve the predictive accuracy of genetic factors in common complex disease such as stroke and cardiovascular disease.

CONCLUSIONS

Our proof-of-concept study demonstrates the power of gene-based meta-analyses for discovering disease-relevant associations complementing common-variant signals by incorporating population-specific rare variation.

Increased proteolytic cleavage of osteoglycin is associated with a stable plaque phenotype and lower risk of cardiovascular events

BACKGROUND AND AIMS

Extracellular matrix (ECM) remodeling is one of the key components in the formation of vulnerable atherosclerotic plaques and cardiovascular events. We recently showed that the full-length ECM-proteoglycan osteoglycin was associated with plaque vulnerability and future cardiovascular events. In the present study, we aimed to investigate the association of cleaved osteoglycin with plaque phenotype.

METHODS

Two-hundred human carotid plaques were analyzed by immunohistochemistry. Cleaved osteoglycin and active caspase-3 were assessed by ELISA. ECM components (collagen, elastin and glycosaminoglycans) were assessed by colorimetric assays in plaque tissue homogenates. Matrix metalloproteinases (MMPs) were assessed using Milliplex. MMP-cleavage of osteoglycin and its effect on apoptosis were studied in vitro. Cardiovascular events were recorded during follow-up using national registries.

RESULTS

Plaque levels of cleaved osteoglycin were significantly higher in asymptomatic plaques and correlated to α-actin plaque area, collagen, elastin and inversely to lipids, active. caspase-3 and a histological vulnerability index. Cleaved osteoglycin correlated to several MMPs, especially MMP-12, which was also shown to cleave osteoglycin in vitro. In vitro cleavage of osteoglycin was also associated with less smooth muscle cell apoptosis. Patients with high plaque levels of cleaved osteoglycin had a significantly lower risk to suffer from future cardiovascular events.

CONCLUSIONS

The current study shows that cleaved osteoglycin is associated with a stable plaque phenotype and lower risk for future cardiovascular events. Potentially due to reduced cell apoptosis and ability to retain LDL. These results indicate that targeting the cleavage of osteoglycin may be a potential therapeutic strategy to stabilize plaques.

Proteomic Profiles of Body Mass Index and Waist-to-Hip Ratio and Their Role in Incidence of Diabetes

CONTEXT

It is unclear to what extent the plasma proteome of abdominal fat distribution differs from that of body mass index, and whether the differences have clinical implications.

OBJECTIVE

To evaluate the difference between the plasma proteomic profiles of body mass index (BMI) and waist-to-hip ratio (WHR), and then examine the identified BMI- or WHR-specific proteins in relation to incidence of diabetes.

METHODS

Data were obtained from the Malmö Diet and Cancer-Cardiovascular Cohort study in the general community. Participants (n = 4203) with no previous diabetes (aged 57.2 ± 6.0 years, 37.8% men) were included. Plasma proteins (n = 136) were measured by the Proseek proximity extension method. BMI- and WHR-specific proteins were identified at baseline using a 2-step iterative resampling approach to optimize internal replicability followed by β coefficient comparisons. The identified proteins were considered internally replicated and were then studied in relation to incident diabetes by Cox proportional hazards regression analysis. The main outcome measure was incident diabetes over a mean follow-up of 20.3 ± 5.9 years.

RESULTS

After excluding 21 overlapping proteins and proteins that did not show significantly different associations with BMI vs WHR, 10 internally replicated proteins were found to be specific to BMI, and 22 were found to be specific to WHR (false discovery rate-adjusted P < .05). Of the WHR-specific proteins, 18 remained associated with diabetes risk after multivariate adjustments, whereas none of the BMI-specific proteins showed associations with diabetes risk.

CONCLUSION

Abdominal fat distribution was associated with some unique characteristics of the plasma proteome that potentially could be related to its additional risk of diabetes beyond general obesity.

Interleukin-6 Receptor Blockade can Increase the Risk of Nonalcoholic Fatty Liver Disease: Indications From Mendelian Randomization

Background: Interleukin-6 receptor (IL-6R) blockade has been approved for inflammation-associated diseases and whether it is effective in treating non-alcoholic fatty liver disease (NAFLD) is still unknown.

Methods: A target-based Mendelian randomization was performed to appraise whether inhibiting the IL-6 signaling pathway via IL-6R blockade can reduce the risk of NAFLD. The previously established genetic proxy SNP rs2228145 was mainly used to appraise the therapeutic effects and the genetic-predicted circulating IL-6 level was treated as the exposure with ∼30,000 samples. The genetic association between SNP rs2228145 (A > C) and NAFLD was obtained from non-FinnGen GWAS (1,483 cases and 17,781controls) and FinnGen GWAS (894 cases and 217,898 controls). The causal effects were estimated using a Wald ratio method and were combined using a fixed-effects meta-analysis. Furthermore, the SNP rs12048091 was employed as another proxy in the sensitivity analysis.

Results: The positive control analysis suggested the SNP rs2228145 can mimic the effects of IL-6R blockade where inhibiting IL-6 signaling can reduce the risk of rheumatoid arthritis [OR = 0.68 (0.58, 0.80)] and coronary heart disease [OR = 0.75 (0.68, 0.84)]. This Mendelian randomization analysis suggested that IL-6R blockade can adversely increase the risk of NAFLD in the non-FinnGen GWAS [OR = 1.99 (1.27, 3.13)] while not significant in the FinnGen consortium. The fixed-effects meta-analysis indicated inhibiting the IL-6 signaling pathway can reduce the risk of NAFLD [OR = 1.80 (1.26, 2.57)]. When including SNP rs12048091 as the genetic instrument, the meta-analysis using two genetic variants also indicated a similar effect on NAFLD [OR = 1.83 (1.32, 2.53)]. There was no heterogeneity in the whole analysis.

Conclusion: Our Mendelian randomization suggested inhibiting the IL-6 signaling pathway via IL-6R blockade might increase the risk of NAFLD, suggesting IL-6R should play a protective role in NAFLD.

Identification of novel genes influencing eosinophil-specific protein levels in asthma families

BACKGROUND

Eosinophils play a key role in the asthma allergic response by releasing cytotoxic molecules such as eosinophil cationic protein (ECP) and eosinophil-derived neurotoxin (EDN) that generate epithelium damages.

OBJECTIVE

To identify genetic variants influencing ECP and EDN levels in asthma-ascertained families.

METHODS

We performed univariate and bivariate genome-wide association analyses of ECP and EDN levels in 1,018 subjects from EGEA study with follow-up in 153 subjects from SLSJ study and combined the results of these two studies through meta-analysis. We then conducted Bayesian statistical fine-mapping together with quantitative trait locus and functional annotation analyses to identify the most likely functional genetic variants and candidate genes.

RESULTS

We identified five genome-wide significant loci (P<5x10^-8) including seven distinct signals associated with ECP and/or EDN levels. The genes targeted by our fine-mapping and functional search include RNASE2 and RNASE3 (14q11) which encode EDN and ECP respectively and four other genes which regulate ECP/EDN levels. These four genes were the following: JAK1 (1p31) a transcription factor with a key role in the immune response and a potential therapeutic target for eosinophilic asthma, ARHGAP25 (2p13) involved in leukocyte recruitment to inflammatory sites, NDUFA4 (7p21) encoding a component of the mitochondrial respiratory chain and involved in cellular response to stress and CTSL (9q22) involved in immune response, extra-cellular remodeling and allergic inflammation.

CONCLUSION

This study demonstrates that the analysis of specific phenotypes produced by eosinophils allows identifying genes with a major role in allergic response and inflammation and offering potential therapeutic targets for asthma.

A Mendelian randomization study investigating the causal role of inflammation on Parkinson’s disease

There is increasing evidence for inflammation as a determinant in the pathogenesis of Parkinson’s disease, but its role in parkinsonian neurodegeneration remains elusive. It is not clear whether inflammatory cascades are causes or consequences of dopamine neuron death. In the present study, we aim to perform an in-depth statistical investigation of the causal relationship between inflammation and Parkinson’s disease using a two-sample Mendelian randomization design.

Genetic instruments were selected using summary-level data from the largest genome-wide association studies to date (sample size ranging from 13 955 to 204 402 individuals) conducted on a European population for the following inflammation biomarkers: C-reactive protein, interleukin-6, interleukin 1 receptor antagonist and tumour necrosis factor α. Genetic association data on Parkinson’s disease (56 306 cases and 1 417 791 controls) and age at onset of Parkinson’s disease (28 568 cases) were obtained from the International Parkinson’s Disease Genomics Consortium. On primary analysis, causal associations were estimated on sets of strong (P-value < 5 × 10⁻⁸; F-statistic > 10) and independent (linkage disequilibrium r² < 0.001) genetic instruments using the inverse-variance weighted method. In sensitivity analysis, we estimated causal effects using robust Mendelian randomization methods and after removing pleiotropic genetic variants. Reverse causation was also explored. We repeated the analysis on different data sources for inflammatory biomarkers to check the consistency of the findings.

In all the three data sources selected for interleukin-6, we found statistical evidence for an earlier age at onset of Parkinson’s disease associated with increased interleukin-6 concentration [years difference per 1 log-unit increase = −2.364, 95% confidence interval (CI) = −4.789–0.060; years difference per 1 log-unit increase = −2.011, 95% CI = −3.706 to −0.317; years difference per 1 log-unit increase = −1.569, 95% CI = −2.891 to −0.247]. We did not observe any statistical evidence for causal effects of C-reactive protein, interleukin 1 receptor antagonist and tumour necrosis factor α on both Parkinson’s disease and its age at onset. Results after excluding possible pleiotropic genetic variants were consistent with findings from primary analyses. When investigating reverse causation, we did not find evidence for a causal effect of Parkinson’s disease or age at onset on any biomarkers of inflammation.

We found evidence for a causal association between the onset of Parkinson’s disease and interleukin-6. The findings of this study suggest that the pro-inflammatory activity of the interleukin-6 cytokine could be a determinant of prodromal Parkinson’s disease.

The heart-brain axis: A proteomics study of meditation on the cardiovascular system of Tibetan Monks

Background

There have been mixed reports on the beneficial effects of meditation in cardiovascular disease (CVD), which is widely considered the leading cause of death worldwide.

Methods

To clarify the role of meditation in modulating the heart-brain axis, we implemented an extreme phenotype strategy, i.e., Tibetan monks (BMI > 30) who practised 19.20 ± 7.82 years of meditation on average and their strictly matched non-meditative Tibetan controls. Hypothesis-free advanced proteomics strategies (Data Independent Acquisition and Targeted Parallel Reaction Monitoring) were jointly applied to systematically investigate and target the plasma proteome underlying meditation. Total cholesterol, low-density lipoprotein cholesterol  (LDL-C), apolipoprotein B (Apo B) and lipoprotein (a) [Lp(a)] as the potential cardiovascular risk factors were evaluated. Heart rate variability (HRV) was assessed by electrocardiogram.

Findings

Obesity, hypertension, and reduced HRV is offset by long-term meditation. Notably, meditative monks have blood pressure and HRV comparable to their matched Tibetan controls. Meditative monks have a protective plasma proteome, related to decreased atherosclerosis, enhanced glycolysis, and oxygen release, that confers resilience to the development of CVD. In addition, clinical risk factors in plasma were significantly decreased in monks compared with controls, including total cholesterol, LDL-C, Apo B, and Lp(a).

Interpretation

To our knowledge, this work is the first well-controlled proteomics investigation of long-term meditation, which opens up a window for individuals characterized by a sedentary lifestyle to improve their cardiovascular health with an accessible method practised for more than two millennia.

Funding

See the Acknowledgements section.

Connecting Genomics and Proteomics to Identify Protein Biomarkers for Adult and Youth-Onset Type 2 Diabetes: A Two-Sample Mendelian Randomization Study

Type 2 diabetes shows an increasing prevalence in both adults and children. Identification of biomarkers for both youth and adult-onset type 2 diabetes is crucial for development of screening tools or drug targets. In this study, using two-sample Mendelian randomization (MR), we identified 22 circulating proteins causally linked to adult type 2 diabetes and 11 proteins with suggestive evidence for association with youth-onset type 2 diabetes. Among these, colocalization analysis further supported a role in type 2 diabetes for C-type mannose receptor 2 (MR odds ratio [OR] 0.85 [95% CI 0.79-0.92] per genetically predicted SD increase in protein level), MANS domain containing 4 (MR OR 0.90 [95% CI 0.88-0.92]), sodium/potassium-transporting ATPase subunit β2 (MR OR 1.10 [95% CI 1.06-1.15]), endoplasmic reticulum oxidoreductase 1β (MR OR 1.09 [95% CI 1.05-1.14]), spermatogenesis-associated protein 20 (MR OR 1.12 [95% CI 1.06-1.18]), haptoglobin (MR OR 0.96 [95% CI 0.94-0.98]), and α1-3-N-acetylgalactosaminyltransferase and α1-3-galactosyltransferase (MR OR 1.04 [95% CI 1.03-1.05]). Our findings support a causal role in type 2 diabetes for a set of circulating proteins, which represent promising type 2 diabetes drug targets.

Causal Associations Between Circulating Adipokines and Cardiovascular Disease: A Mendelian Randomization Study

CONTEXT

Observational studies have suggested associations between adipokines and cardiovascular disease (CVD), but the roles of certain adipokines remain controversial, and these associations have not yet been ascertained causally.

OBJECTIVE

To investigate whether circulating adipokines causally affect the risk of CVD using 2-sample Mendelian randomization (MR).

METHODS

Independent genetic variants strongly associated with adiponectin, resistin, chemerin, and retinol binding protein 4 (RBP4) were selected from public genome-wide association studies. Summary-level statistics for CVD, including coronary artery disease (CAD), myocardial infarction, atrial fibrillation (AF), heart failure (HF), and stroke and its subtypes were collected. The inverse-variance weighted and Wald ratio methods were used for the MR estimates. The MR pleiotropy residual sum and outlier, weighted median, MR-Egger, leave-one-out analysis, MR Steiger, and colocalization analyses were used in the sensitivity analysis.

RESULTS

Genetically predicted resistin levels were positively associated with AF risk (odds ratio [OR] 1.09; 95% confidence interval [CI], 1.04-1.13; P = 4.1 × 10-5), which was attenuated to null after adjusting for blood pressure. We observed suggestive associations between higher genetically predicted chemerin levels and an increased risk of CAD (OR 1.27; 95% CI, 1.01-1.60; P = 0.040), higher genetically predicted RBP4 levels and an increased risk of HF (OR 1.14; 95% CI, 1.02-1.27; P = 0.024). There was no causal association between genetically predicted adiponectin levels and CVD risk.

CONCLUSIONS

Our findings reveal the causal association between resistin and AF, probably acting through blood pressure, and suggest potential causal associations between chemerin and CAD, RBP4, and HF.

Genetic analyses identify pleiotropy and causality for blood proteins and highlight Wnt/β-catenin signalling in migraine

Migraine is a common complex disorder with a significant polygenic SNP heritability (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${h}_{{SNP}}^{2}$$\end{document}hSNP2). Here we utilise genome-wide association study (GWAS) summary statistics to study pleiotropy between blood proteins and migraine under the polygenic model. We estimate \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${h}_{{SNP}}^{2}$$\end{document}hSNP2 for 4625 blood protein GWASs and identify 325 unique proteins with a significant \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${h}_{{SNP}}^{2}$$\end{document}hSNP2 for use in subsequent genetic analyses. Pleiotropy analyses link 58 blood proteins to migraine risk at genome-wide, gene and/or single-nucleotide polymorphism levels—suggesting shared genetic influences or causal relationships. Notably, the identified proteins are largely distinct from migraine GWAS loci. We show that higher levels of DKK1 and PDGFB, and lower levels of FARS2, GSTA4 and CHIC2 proteins have a significant causal effect on migraine. The risk-increasing effect of DKK1 is particularly interesting—indicating a role for downregulation of β-catenin-dependent Wnt signalling in migraine risk, suggesting Wnt activators that restore Wnt/β-catenin signalling in brain could represent therapeutic tools against migraine.

Understanding of the causes and treatment of migraine is incomplete. Here, the authors detect pleiotropic genetic effects and causal relationships between migraine and 58 proteins that are largely distinct from migraine-associated loci identified by GWAS.

An optimal variant to gene distance window derived from an empirical definition of cis and trans protein QTLs

Background

A genome-wide association study (GWAS) correlates variation in the genotype with variation in the phenotype across a cohort, but the causal gene mediating that impact is often unclear. When the phenotype is protein abundance, a reasonable hypothesis is that the gene encoding that protein is the causal gene. However, as variants impacting protein levels can occur thousands or even millions of base pairs from the gene encoding the protein, it is unclear at what distance this simple hypothesis breaks down.

Results

By making the simple assumption that cis-pQTLs should be distance dependent while trans-pQTLs are distance independent, we arrive at a simple and empirical distance cutoff separating cis- and trans-pQTLs. Analyzing a recent large-scale pQTL study (Pietzner in Science 374:eabj1541, 2021) we arrive at an estimated distance cutoff of 944 kilobasepairs (95% confidence interval: 767–1,161) separating the cis and trans regimes.

Conclusions

We demonstrate that this simple model can be applied to other molecular GWAS traits. Since much of biology is built on molecular traits like protein, transcript and metabolite abundance, we posit that the mathematical models for cis and trans distance distributions derived here will also apply to more complex phenotypes and traits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04706-x.