Proteomic Predictors of Incident Diabetes: Results From the Atherosclerosis Risk in Communities (ARIC) Study

Proteome-wide association study using cis and trans variants and applied to blood cell and lipid-related traits in the Women's Health Initiative study

In most Proteome-Wide Association Studies (PWAS), variants near the protein-coding gene (±1 Mb), also known as cis single nucleotide polymorphisms (SNPs), are used to predict protein levels, which are then tested for association with phenotypes. However, proteins can be regulated through variants outside of the cis region. An intermediate GWAS step to identify protein quantitative trait loci (pQTL) allows for the inclusion of trans SNPs outside the cis region in protein-level prediction models. Here, we assess the prediction of 540 proteins in 1002 individuals from the Women's Health Initiative (WHI), split equally into a GWAS set, an elastic net training set, and a testing set. We compared the testing r2 between measured and predicted protein levels using this proposed approach, to the testing r2 using only cis SNPs. The two methods usually resulted in similar testing r2, but some proteins showed a significant increase in testing r2 with our method. For example, for cartilage acidic protein 1, the testing r2 increased from 0.101 to 0.351. We also demonstrate reproducible findings for predicted protein association with lipid and blood cell traits in WHI participants without proteomics data and in UK Biobank utilizing our PWAS weights.

Proteomic Analyses in Diverse Populations Improved Risk Prediction and Identified New Drug Targets for Type 2 Diabetes

OBJECTIVE

Integrated analyses of plasma proteomics and genetic data in prospective studies can help assess the causal relevance of proteins, improve risk prediction, and discover novel protein drug targets for type 2 diabetes (T2D).

RESEARCH DESIGN AND METHODS

We measured plasma levels of 2,923 proteins using Olink Explore among ∼2,000 randomly selected participants from China Kadoorie Biobank (CKB) without prior diabetes at baseline. Cox regression assessed associations of individual protein with incident T2D (n = 92 cases). Proteomic-based risk models were developed with discrimination, calibration, reclassification assessed using area under the curve (AUC), calibration plots, and net reclassification index (NRI), respectively. Two-sample Mendelian randomization (MR) analyses using cis-protein quantitative trait loci identified in a genome-wide association study of CKB and UK Biobank for specific proteins were conducted to assess their causal relevance for T2D, along with colocalization analyses to examine shared causal variants between proteins and T2D.

RESULTS

Overall, 33 proteins were significantly associated (false discovery rate <0.05) with risk of incident T2D, including IGFBP1, GHR, and amylase. The addition of these 33 proteins to a conventional risk prediction model improved AUC from 0.77 (0.73-0.82) to 0.88 (0.85-0.91) and NRI by 38%, with predicted risks well calibrated with observed risks. MR analyses provided support for the causal relevance for T2D of ENTR1, LPL, and PON3, with replication of ENTR1 and LPL in Europeans using different genetic instruments. Moreover, colocalization analyses showed strong evidence (pH4 > 0.6) of shared genetic variants of LPL and PON3 with T2D.

CONCLUSIONS

Proteomic analyses in Chinese adults identified novel associations of multiple proteins with T2D with strong genetic evidence supporting their causal relevance and potential as novel drug targets for prevention and treatment of T2D.

Plasma Neuronal Growth Regulator 1 May Link Physical Activity to Reduced Risk of Type 2 Diabetes: A Proteome-Wide Study of ARIC Participants

Habitual physical activity (PA) impacts the plasma proteome and reduces the risk of developing type 2 diabetes (T2D). Using a large-scale proteome-wide approach in Atherosclerosis Risk in Communities study participants, we aimed to identify plasma proteins associated with PA and determine which of these may be causally related to lower T2D risk. PA was associated with 92 plasma proteins in discovery (P < 1.01 × 10-5), and 40 remained significant in replication (P < 5.43 × 10-4). Eighteen of these proteins were independently associated with incident T2D (P < 1.25 × 10-3), including neuronal growth regulator 1 (NeGR1; hazard ratio per SD 0.85; P = 7.5 × 10-11). Two-sample Mendelian randomization (MR) inverse variance weighted analysis indicated that higher NeGR1 reduces T2D risk (odds ratio [OR] per SD 0.92; P = 0.03) and was consistent with MR-Egger, weighted median, and weighted mode sensitivity analyses. A stronger association was observed for the single cis-acting NeGR1 genetic variant (OR per SD 0.80; P = 6.3 × 10-5). Coupled with previous evidence that low circulating NeGR1 levels promote adiposity, its association with PA and potential causal role in T2D shown here suggest that NeGR1 may link PA exposure with metabolic outcomes. Further research is warranted to confirm our findings and examine the interplay of PA, NeGR1, adiposity, and metabolic health.

ARTICLE HIGHLIGHTS

Proteomic and Metabolomic Signatures of Diet Quality in Young Adults

The assessment of "omics" signatures may contribute to personalized medicine and precision nutrition. However, the existing literature is still limited in the homogeneity of participants' characteristics and in limited assessments of integrated omics layers. Our objective was to use post-prandial metabolomics and fasting proteomics to identify biological pathways and functions associated with diet quality in a population of primarily Hispanic young adults. We conducted protein and metabolite-wide association studies and functional pathway analyses to assess the relationships between a priori diet indices, Healthy Eating Index-2015 (HEI) and Dietary Approaches to Stop Hypertension (DASH) diets, and proteins (n = 346) and untargeted metabolites (n = 23,173), using data from the MetaAIR study (n = 154, 61% Hispanic). Analyses were performed for each diet quality index separately, adjusting for demographics and BMI. Five proteins (ACY1, ADH4, AGXT, GSTA1, F7) and six metabolites (undecylenic acid, betaine, hyodeoxycholic acid, stearidonic acid, iprovalicarb, pyracarbolid) were associated with both diets (p < 0.05), though none were significant after adjustment for multiple comparisons. Overlapping proteins are involved in lipid and amino acid metabolism and in hemostasis, while overlapping metabolites include amino acid derivatives, bile acids, fatty acids, and pesticides. Enriched biological pathways were involved in macronutrient metabolism, immune function, and oxidative stress. These findings in young Hispanic adults contribute to efforts to develop precision nutrition and medicine for diverse populations.

Plasma proteins and onset of type 2 diabetes and diabetic complications: Proteome-wide Mendelian randomization and colocalization analyses

We conduct proteome-wide Mendelian randomization and colocalization analyses to decipher the associations of blood proteins with the risk of type 2 diabetes and diabetic complications. Genetic data on plasma proteome are obtained from 54,306 UK Biobank participants and 35,559 Icelanders. Summary-level data on type 2 diabetes are obtained from the DIAGRAM (DIAbetes Genetics Replication And Meta-analysis consortium) consortium (74,124 cases) and FinnGen study (33,043 cases). Data on 10 diabetic complications are obtained from FinnGen and corresponding studies. Among 1,886 proteins, genetically predicted levels of 47 plasma proteins are associated with type 2 diabetes. Eleven of these proteins have strong support of colocalization. Seventeen proteins are associated with at least one diabetic complication, although a few have colocalization support. HLA-DRA, AGER, HSPA1A, and HSPA1B are associated with most microvascular complications. This study reveals causal proteins for the onset of type 2 diabetes and diabetic complications, which enhances the understanding of molecular etiology and development of therapeutics.