Hypertension affects more than one billion people worldwide. Here we identify 113 novel loci, reporting a total of 2,103 independent genetic signals (P < 5 × 10-8) from the largest single-stage blood pressure (BP) genome-wide association study to date (n = 1,028,980 European individuals). These associations explain more than 60% of single nucleotide polymorphism-based BP heritability. Comparing top versus bottom deciles of polygenic risk scores (PRSs) reveals clinically meaningful differences in BP (16.9 mmHg systolic BP, 95% CI, 15.5-18.2 mmHg, P = 2.22 × 10-126) and more than a sevenfold higher odds of hypertension risk (odds ratio, 7.33; 95% CI, 5.54-9.70; P = 4.13 × 10-44) in an independent dataset. Adding PRS into hypertension-prediction models increased the area under the receiver operating characteristic curve (AUROC) from 0.791 (95% CI, 0.781-0.801) to 0.826 (95% CI, 0.817-0.836, ∆AUROC, 0.035, P = 1.98 × 10-34). We compare the 2,103 loci results in non-European ancestries and show significant PRS associations in a large African-American sample. Secondary analyses implicate 500 genes previously unreported for BP. Our study highlights the role of increasingly large genomic studies for precision health research.

Background: Smoking is associated with carotid intima-media thickness (C-IMT). However, knowledge about how genetics may influence this association is limited. We aimed to perform nonhypothesis driven gene-smoking interaction analyses to identify potential genetic variants, among those included in immune and metabolic platforms, that may modify the effect of smoking on carotid intima-media thickness. Methods: We used baseline data from 1551 men and 1700 women, aged 55 to 79, included in a European multi-center study. Carotid intima-media thickness maximum, the maximum of values measured at different locations of the carotid tree, was dichotomized with cut point values ≥75, respectively. Genetic data were retrieved through use of the Illumina Cardio-Metabo- and Immuno- Chips. Gene-smoking interactions were evaluated through calculations of Synergy index (S). After adjustments for multiple testing, P values of <2.4×10 −7 for S were considered significant. The models were adjusted for age, sex, education, physical activity, type of diet, and population stratification. Results: Our screening of 207 586 SNPs available for analysis, resulted in the identification of 47 significant gene-smoking synergistic interactions in relation to carotid intima-media thickness maximum. Among the significant SNPs, 28 were in protein coding genes, 2 in noncoding RNA and the remaining 17 in intergenic regions. Conclusions: Through nonhypothesis-driven analyses of gene-smoking interactions, several significant results were observed. These may stimulate further research on the role of specific genes in the process that determines the effect of smoking habits on the development of carotid atherosclerosis.

Insulin secretion is critical for glucose homeostasis, and increased levels of the precursor proinsulin relative to insulin indicate pancreatic islet beta-cell stress and insufficient insulin secretory capacity in the setting of insulin resistance. We conducted meta-analyses of genome-wide association results for fasting proinsulin from 16 European-ancestry studies in 45,861 individuals. We found 36 independent signals at 30 loci (p value < 5 × 10-8), which validated 12 previously reported loci for proinsulin and ten additional loci previously identified for another glycemic trait. Half of the alleles associated with higher proinsulin showed higher rather than lower effects on glucose levels, corresponding to different mechanisms. Proinsulin loci included genes that affect prohormone convertases, beta-cell dysfunction, vesicle trafficking, beta-cell transcriptional regulation, and lysosomes/autophagy processes. We colocalized 11 proinsulin signals with islet expression quantitative trait locus (eQTL) data, suggesting candidate genes, including ARSG, WIPI1, SLC7A14, and SIX3. The NKX6-3/ANK1 proinsulin signal colocalized with a T2D signal and an adipose ANK1 eQTL signal but not the islet NKX6-3 eQTL. Signals were enriched for islet enhancers, and we showed a plausible islet regulatory mechanism for the lead signal in the MADD locus. These results show how detailed genetic studies of an intermediate phenotype can elucidate mechanisms that may predispose one to disease.

Conventional measurements of fasting and postprandial blood glucose levels investigated in genome-wide association studies (GWAS) cannot capture the effects of DNA variability on ‘around the clock’ glucoregulatory processes. Here we show that GWAS meta-analysis of glucose measurements under nonstandardized conditions (random glucose (RG)) in 476,326 individuals of diverse ancestries and without diabetes enables locus discovery and innovative pathophysiological observations. We discovered 120 RG loci represented by 150 distinct signals, including 13 with sex-dimorphic effects, two cross-ancestry and seven rare frequency signals. Of these, 44 loci are new for glycemic traits. Regulatory, glycosylation and metagenomic annotations highlight ileum and colon tissues, indicating an underappreciated role of the gastrointestinal tract in controlling blood glucose. Functional follow-up and molecular dynamics simulations of lower frequency coding variants in glucagon-like peptide-1 receptor ( GLP1R ), a type 2 diabetes treatment target, reveal that optimal selection of GLP-1R agonist therapy will benefit from tailored genetic stratification. We also provide evidence from Mendelian randomization that lung function is modulated by blood glucose and that pulmonary dysfunction is a diabetes complication. Our investigation yields new insights into the biology of glucose regulation, diabetes complications and pathways for treatment stratification.

Objective Neonatal lupus erythematosus (NLE) may develop after transplacental transfer of maternal autoantibodies with cardiac manifestations (congenital heart block, CHB) including atrioventricular block, atrial and ventricular arrhythmias, and cardiomyopathies. The association with anti-Ro/SSA antibodies is well established, but a recurrence rate of only 12%–16% despite persisting maternal autoantibodies suggests that additional factors are required for CHB development. Here, we identify fetal genetic variants conferring risk of CHB and elucidate their effects on cardiac function. Methods A genome-wide association study was performed in families with at least one case of CHB. Gene expression was analysed by microarrays, RNA sequencing and PCR and protein expression by western blot, immunohistochemistry, immunofluorescence and flow cytometry. Calcium regulation and connectivity were analysed in primary cardiomyocytes and cells induced from pleuripotent stem cells. Fetal heart performance was analysed by Doppler/echocardiography. Results We identified DNAJC6 as a novel fetal susceptibility gene, with decreased cardiac expression of DNAJC6 associated with the disease risk genotype. We further demonstrate that fetal cardiomyocytes deficient in auxilin, the protein encoded by DNAJC6, have abnormal connectivity and Ca2+ homoeostasis in culture, as well as decreased cell surface expression of the Cav1.3 calcium channel. Doppler echocardiography of auxilin-deficient fetal mice revealed cardiac NLE abnormalities in utero, including abnormal heart rhythm with atrial and ventricular ectopias, as well as a prolonged atrioventricular time intervals. Conclusions Our study identifies auxilin as the first genetic susceptibility factor in NLE modulating cardiac function, opening new avenues for the development of screening and therapeutic strategies in CHB.

Objective: To identify causal pathophysiological mechanisms for atherosclerosis and incident cardiovascular events using protein measurements. Approach and Results: Carotid artery atherosclerosis was assessed by ultrasound, and 86 cardiovascular-related proteins were measured using the Olink CVD-I panel in 7 Swedish prospective studies (11 754 individuals). The proteins were analyzed in relation to intima-media thickness in the common carotid artery (IMT-CCA), plaque occurrence, and incident cardiovascular events (composite end point of myocardial infarction or ischemic stroke) using a discovery/replication approach in different studies. After adjustments for traditional cardiovascular risk factors, 11 proteins remained significantly associated with IMT-CCA in the replication stage, whereas 9 proteins were replicated for plaque occurrence and 17 proteins for incident cardiovascular events. NT-proBNP (N-terminal pro-B-type natriuretic peptide) and MMP (matrix metalloproteinase)-12 were associated with both IMT-CCA and incident events, but the overlap was considerably larger between plaque occurrence and incident events, including MMP-12, TIM-1 (T-cell immunoglobulin and mucin domain 1), GDF (growth/differentiation factor)-15, IL (interleukin)-6, U-PAR (urokinase plasminogen activator surface receptor), LOX-1 (lectin-like oxidized LDL [low-density lipoprotein] receptor 1), and TRAIL-R2 (TNF [tumor necrosis factor]-related apoptosis-inducing ligand receptor 2). Only MMP-12 was associated with IMT-CCA, plaque, and incident events with a positive and concordant direction of effect. However, a 2-sample Mendelian randomization analysis suggested that increased MMP-12 may be protective against ischemic stroke ( P =5.5×10 −7 ), which is in the opposite direction of the observational analyses. Conclusions: The present meta-analysis discovered several proteins related to carotid atherosclerosis that partly differed in their association with IMT-CCA, plaque, and incident atherosclerotic disease. Mendelian randomization analysis for the top finding, MMP-12, suggests that the increased levels of MMP-12 could be a consequence of atherosclerotic burden rather than the opposite chain of events.

Interleukin 6 (IL-6) is a multifunctional cytokine with both pro- and anti-inflammatory properties with a heritability estimate of up to 61%. The circulating levels of IL-6 in blood have been associated with an increased risk of complex disease pathogenesis. We conducted a two-staged, discovery and replication meta genome-wide association study (GWAS) of circulating serum IL-6 levels comprising up to 67 428 (ndiscovery = 52 654 and nreplication = 14 774) individuals of European ancestry. The inverse variance fixed effects based discovery meta-analysis, followed by replication led to the identification of two independent loci, IL1F10/IL1RN rs6734238 on chromosome (Chr) 2q14, (Pcombined = 1.8 × 10-11), HLA-DRB1/DRB5 rs660895 on Chr6p21 (Pcombined = 1.5 × 10-10) in the combined meta-analyses of all samples. We also replicated the IL6R rs4537545 locus on Chr1q21 (Pcombined = 1.2 × 10-122). Our study identifies novel loci for circulating IL-6 levels uncovering new immunological and inflammatory pathways that may influence IL-6 pathobiology.

Glycemic traits are used to diagnose and monitor type 2 diabetes and cardiometabolic health. To date, most genetic studies of glycemic traits have focused on individuals of European ancestry. Here we aggregated genome-wide association studies comprising up to 281,416 individuals without diabetes (30% non-European ancestry) for whom fasting glucose, 2-h glucose after an oral glucose challenge, glycated hemoglobin and fasting insulin data were available. Trans-ancestry and single-ancestry meta-analyses identified 242 loci (99 novel; P < 5 × 10−8), 80% of which had no significant evidence of between-ancestry heterogeneity. Analyses restricted to individuals of European ancestry with equivalent sample size would have led to 24 fewer new loci. Compared with single-ancestry analyses, equivalent-sized trans-ancestry fine-mapping reduced the number of estimated variants in 99% credible sets by a median of 37.5%. Genomic-feature, gene-expression and gene-set analyses revealed distinct biological signatures for each trait, highlighting different underlying biological pathways. Our results increase our understanding of diabetes pathophysiology by using trans-ancestry studies for improved power and resolution. A trans-ancestry meta-analysis of GWAS of glycemic traits in up to 281,416 individuals identifies 99 novel loci, of which one quarter was found due to the multi-ancestry approach, which also improves fine-mapping of credible variant sets.

Homeostatic control of blood glucose requires different physiological responses in the fasting and post-prandial states. We reasoned that glucose measurements under non-standardised conditions (random glucose; RG) may capture diverse glucoregulatory processes more effectively than previous genome-wide association studies (GWAS) of fasting glycaemia or after standardised glucose loads. Through GWAS meta-analysis of RG in 493,036 individuals without diabetes of diverse ethnicities we identified 128 associated loci represented by 162 distinct signals, including 14 with sex-dimorphic effects, 9 discovered through trans-ethnic analysis, and 70 novel signals for glycaemic traits. Novel RG loci were particularly enriched in expression in the ileum and colon, indicating a prominent role for the gastrointestinal tract in the control of blood glucose. Functional studies and molecular dynamics simulations of coding variants of GLP1R , a well-established type 2 diabetes treatment target, provided a genetic framework for optimal selection of GLP-1R agonist therapy. We also provided new evidence from Mendelian randomisation that lung function is modulated by blood glucose and that pulmonary dysfunction is a diabetes complication. Thus, our approach based on RG GWAS provided wide-ranging insights into the biology of glucose regulation, diabetes complications and the potential for treatment stratification.

A Correction to this paper has been published: https://doi.org/10.1038/s41467-021-21276-3

Understanding why individuals with severe mental illness (Schizophrenia, Bipolar Disorder and Major Depressive Disorder) have increased risk of cardiometabolic disease (including obesity, type 2 diabetes and cardiovascular disease), and identifying those at highest risk of cardiometabolic disease are important priority areas for researchers. For individuals with European ancestry we explored whether genetic variation could identify sub-groups with different metabolic profiles. Loci associated with schizophrenia, bipolar disorder and major depressive disorder from previous genome-wide association studies and loci that were also implicated in cardiometabolic processes and diseases were selected. In the IMPROVE study (a high cardiovascular risk sample) and UK Biobank (general population sample) multidimensional scaling was applied to genetic variants implicated in both psychiatric and cardiometabolic disorders. Visual inspection of the resulting plots used to identify distinct clusters. Differences between these clusters were assessed using chi-squared and Kruskall-Wallis tests. In IMPROVE, genetic loci associated with both schizophrenia and cardiometabolic disease (but not bipolar disorder or major depressive disorder) identified three groups of individuals with distinct metabolic profiles. This grouping was replicated within UK Biobank, with somewhat less distinction between metabolic profiles. This work focused on individuals of European ancestry and is unlikely to apply to more genetically diverse populations. Overall, this study provides proof of concept that common biology underlying mental and physical illness may help to stratify subsets of individuals with different cardiometabolic profiles.

The relationship between intake of saturated fats and subclinical atherosclerosis, as well as the possible influence of genetic variants, is poorly understood and investigated. We aimed to investigate this relationship, with a hypothesis that it would be positive, and to explore whether genetics may modulate it, using data from a European cohort including 3,407 participants aged 54–79 at high risk of cardiovascular disease. Subclinical atherosclerosis was assessed by carotid intima-media thickness (C-IMT), measured at baseline and after 30 months. Logistic regression (OR; 95% CI) was employed to assess the association between high intake of food rich in saturated fat (vs. low) and: (1) the mean and the maximum values of C-IMT in the whole carotid artery (C-IMT mean , C-IMT max ), in the bifurcation (Bif-), the common (CC-) and internal (ICA-) carotid arteries at baseline (binary, cut-point ≥ 75th), and (2) C-IMT progression (binary, cut-point > zero). For the genetic-diet interaction analyses, we considered 100,350 genetic variants. We defined interaction as departure from additivity of effects. After age- and sex-adjustment, high intake of saturated fat was associated with increased C-IMT mean (OR:1.27;1.06–1.47), CC-IMT mean (OR:1.22;1.04–1.44) and ICA-IMT mean (OR:1.26;1.07–1.48). However, in multivariate analysis results were no longer significant. No clear associations were observed between high intake of saturated fat and risk of atherosclerotic progression. There was no evidence of interactions between high intake of saturated fat and any of the genetic variants considered, after multiple testing corrections. High intake of saturated fats was not independently associated with subclinical atherosclerosis. Moreover, we did not identify any significant genetic-dietary fat interactions in relation to risk of subclinical atherosclerosis.

Background Use of targeted exome-arrays with common, rare variants and functionally enriched variation has led to discovery of new genes contributing to population variation in risk factors. Plasminogen activator-inhibitor 1 (PAI-1), tissue plasminogen activator (tPA), and the plasma product D-dimer are important components of the fibrinolytic system. There have been few large-scale genome-wide or exome-wide studies of PAI-1, tPA, and D-dimer. Objectives We sought to discover new genetic loci contributing to variation in these traits using an exome-array approach. Methods Cohort-level analyses and fixed effects meta-analyses of PAI-1 (n = 15 603), tPA (n = 6876,) and D-dimer (n = 19 306) from 12 cohorts of European ancestry with diverse study design were conducted, including single-variant analyses and gene-based burden testing. Results Five variants located in NME7, FGL1, and the fibrinogen locus, all associated with D-dimer levels, achieved genome-wide significance (P < 5 x 10(-8)). Replication was sought for these 5 variants, as well as 45 well-imputed variants with P < 1 x 10(-4) in the discovery using an independent cohort. Replication was observed for three out of the five significant associations, including a novel and uncommon (0.013 allele frequency) coding variant p.Trp256Leu in FGL1 (fibrinogen-like-1) with increased plasma D-dimer levels. Additionally, a candidate-gene approach revealed a suggestive association for a coding variant (rs143202684-C) in SERPINB2, and suggestive associations with consistent effect in the replication analysis include an intronic variant (rs11057830-A) in SCARB1 associated with increased D-dimer levels. Conclusion This work provides new evidence for a role of FGL1 in hemostasis.

Affinity-based proteomics assays identify particular profiles of the circulating histidine-rich glycoprotein to increase with age and elevated levels to indicate a higher risk of all-cause mortality. Despite recognizing aging as a common risk factor of many human diseases, little is known about its molecular traits. To identify age-associated proteins circulating in human blood, we screened 156 individuals aged 50–92 using exploratory and multiplexed affinity proteomics assays. Profiling eight additional study sets (N = 3,987), performing antibody validation, and conducting a meta-analysis revealed a consistent age association (P = 6.61 × 10−6) for circulating histidine-rich glycoprotein (HRG). Sequence variants of HRG influenced how the protein was recognized in the immunoassays. Indeed, only the HRG profiles affected by rs9898 were associated with age and predicted the risk of mortality (HR = 1.25 per SD; 95% CI = 1.12–1.39; P = 6.45 × 10−5) during a follow-up period of 8.5 yr after blood sampling (IQR = 7.7–9.3 yr). Our affinity proteomics analysis found associations between the particular molecular traits of circulating HRG with age and all-cause mortality. The distinct profiles of this multipurpose protein could serve as an accessible and informative indicator of the physiological processes related to biological aging.

Background/Aim The association between alcohol consumption and subclinical atherosclerosis is still unclear. Using data from a European multicentre study, we assess subclinical atherosclerosis and its 30-month progression by carotid intima-media thickness (C-IMT) measurements, and correlate this information with self-reported data on alcohol consumption. Methods Between 2002–2004, 1772 men and 1931 women aged 54–79 years with at least three risk factors for cardiovascular disease (CVD) were recruited in Italy, France, Netherlands, Sweden, and Finland. Self-reported alcohol consumption, assessed at baseline, was categorized as follows: none (0 g/d), very-low (0 − 5 g/d), low (> 5 to ≤ 10 g/d), moderate (> 10 to ≤ 20 g/d for women, > 10 to ≤ 30 g/d for men) and high (> 20 g/d for women, > 30 g/d for men). C-IMT was measured in millimeters at baseline and after 30 months. Measurements consisted of the mean and maximum values of the common carotids (CC), internal carotid artery (ICA), and bifurcations (Bif) and whole carotid tree. We used quantile regression to describe the associations between C-IMT measures and alcohol consumption categories, adjusting for sex, age, physical activity, education, smoking, diet, and latitude. Results Adjusted differences between median C-IMT values in different levels of alcohol consumption (vs. very-low) showed that moderate alcohol consumption was associated with lower C-IMT max [− 0.17(95%CI − 0.32; − 0.02)], and Bif-IMT mean [− 0.07(95%CI − 0.13; − 0.01)] at baseline and decreasing C-IMT mean [− 0.006 (95%CI − 0.011; − 0.000)], Bif-IMT mean [− 0.016(95%CI − 0.027; − 0.005)], ICA-IMT mean [− 0.009(95% − 0.016; − 0.002)] and ICA-IMT max [− 0.016(95%: − 0.032; − 0.000)] after 30 months. There was no evidence of departure from linearity in the association between alcohol consumption and C-IMT. Conclusion In this European population at high risk of CVD, findings show an inverse relation between moderate alcohol consumption and carotid subclinical atherosclerosis and its 30-month progression, independently of several potential confounders.

Circulating proteins are vital in human health and disease and are frequently used as biomarkers for clinical decision-making or as targets for pharmacological intervention. Here, we map and replicate protein quantitative trait loci (pQTL) for 90 cardiovascular proteins in over 30,000 individuals, resulting in 451 pQTLs for 85 proteins. For each protein, we further perform pathway mapping to obtain trans-pQTL gene and regulatory designations. We substantiate these regulatory findings with orthogonal evidence for trans-pQTLs using mouse knockdown experiments (ABCA1 and TRIB1) and clinical trial results (chemokine receptors CCR2 and CCR5), with consistent regulation. Finally, we evaluate known drug targets, and suggest new target candidates or repositioning opportunities using Mendelian randomization. This identifies 11 proteins with causal evidence of involvement in human disease that have not previously been targeted, including EGF, IL-16, PAPPA, SPON1, F3, ADM, CASP-8, CHI3L1, CXCL16, GDF15 and MMP-12. Taken together, these findings demonstrate the utility of large-scale mapping of the genetics of the proteome and provide a resource for future precision studies of circulating proteins in human health.

BACKGROUND AND AIMS:Genetic loci associated with CHD show different relationships with intima-media thickness in the common carotid artery (IMT-CCA) and in the bulb (IMT-bulb). We evaluated if IMT-CCA and IMT-bulb differ also with respect to circulating protein profiles and risk of incident atherosclerotic disease. METHODS:In three Swedish cohorts (MDC, IMPROVE, PIVUS, total n > 7000), IMT-CCA and IMT-bulb were assessed by ultrasound at baseline, and 86 cardiovascular-related proteins were analyzed. In the PIVUS study only, IMT-CCA and IMT-bulb were investigated in relation to incident atherosclerotic disease over 10 years of follow-up. RESULTS:In a meta-analysis of the analysis performed separately in the cohorts, three proteins, matrix metalloproteinase-12 (MMP-12), hepatocyte growth factor (HGF) and N-terminal pro-B-type natriuretic peptide (NT-proBNP), were associated with IMT-CCA when adjusted for traditional cardiovascular risk factors. Five proteins were associated with IMT-bulb (MMP-12, growth/differentiation factor 15 (GDF-15), osteoprotegerin, growth hormone and renin). Following adjustment for cardiovascular risk factors, IMT-bulb was significantly more closely related to incident stroke or myocardial infarction (total number of cases, 111) than IMT-CCA in the PIVUS study (HR 1.51 for 1 SD, 95%CI 1.21-1.87, p < 0.001 vs HR 1.17, 95%CI 0.93-1.47, p = 0.16). MMP-12 levels were related to this combined end-point (HR 1.30, 95%CI 1.08-1.56, p = 0.0061). CONCLUSIONS:Elevated levels of MMP-12 were associated with both IMT-CCA and IMT-bulb, but other proteins were significantly related to IMT in only one of these locations. The finding that IMT-bulb was more closely related to incident atherosclerotic disease than IMT-CCA emphasizes a difference between these measurements of IMT.

Background: Coronary artery disease (CAD) is accelerated in subjects with type 2 diabetes mellitus (T2D). Methods: To test whether this reflects differential genetic influences on CAD risk in subjects with T2D, we performed a systematic assessment of genetic overlap between CAD and T2D in 66 643 subjects (27 708 with CAD and 24 259 with T2D). Variants showing apparent association with CAD in stratified analyses or evidence of interaction were evaluated in a further 117 787 subjects (16 694 with CAD and 11 537 with T2D). Results: None of the previously characterized CAD loci was found to have specific effects on CAD in T2D individuals, and a genome-wide interaction analysis found no new variants for CAD that could be considered T2D specific. When we considered the overall genetic correlations between CAD and its risk factors, we found no substantial differences in these relationships by T2D background. Conclusions: This study found no evidence that the genetic architecture of CAD differs in those with T2D compared with those without T2D.

The genes regulating circulating levels of soluble gp130 (sgp130), the antagonist of the inflammatory response in atherosclerosis driven by interleukin 6, are largely unknown. Aims of the present study were to identify genetic loci associated with circulating sgp130 and to explore the potential association between variants associated with sgp130 and markers of subclinical atherosclerosis. The study is based on IMPROVE ( n = 3703), a cardiovascular multicentre study designed to investigate the determinants of carotid intima media thickness, a measure of subclinical atherosclerosis. Genomic DNA was genotyped by the CardioMetaboChip and ImmunoChip. About 360,842 SNPs were tested for association with log-transformed sgp130, using linear regression adjusted for age, gender, and population stratification using PLINK v1.07. A p value of 1 × 10 −5 was chosen as threshold for significance value. In an exploratory analysis, SNPs associated with sgp130 were tested for association with c-IMT measures. We identified two SNPs significantly associated with sgp130 levels and 24 showing suggestive association with sgp130 levels. One SNP (rs17688225) on chromosome 14 was positively associated with sgp130 serum levels ( β = 0.03 SE = 0.007, p = 4.77 × 10 −5 ) and inversely associated with c-IMT (c-IMT mean–max β = −0.001 SE = 0.005, p = 0.0342). Our data indicate that multiple loci regulate sgp130 levels and suggest a possible common pathway between sgp130 and c-IMT measures.

Circulating proteins are vital in human health and disease and are frequently used as biomarkers for clinical decision-making or as targets for pharmacological intervention. By mapping and replicating protein quantitative trait loci (pQTL) for 90 cardiovascular proteins in over 30,000 individuals, we identified 467 pQTLs for 85 proteins. The pQTLs were used in combination with other sources of information to evaluate known drug targets, and suggest new target candidates or repositioning opportunities, underpinned by a) causality assessment using Mendelian randomization, b) pathway mapping using -pQTL gene assignments, and c) protein-centric polygenic risk scores enabling matching of plausible target mechanisms to sub-groups of individuals enabling precision medicine.