Elucidating mechanisms of genetic cross-disease associations at the PROCR vascular disease locus

Figures

• 

  Fig. 1. Schematic overview of the study design to elucidate molecular underpinnings of cross-disease associations. Credits: The immune response, Big Picture (https://www.stem.org.uk/rx34vg).

• 

  Fig. 2. Association of PROCR-219Gly with a range of health outcomes and circulating cardiovascular biomarkers. a Phenome-wide association scan of PROCR-p.S219G (rs867186) across 1402 broad electronic health record-derived ICD-codes from the UK Biobank. Unadjusted P values were obtained from the PheWeb portal. b Forest plot showing the associations of the minor (G) allele of rs867186 genotype with different cardiovascular conditions. Association statistics for VTE outcomes were obtained from the INVENT consortium (a) or UK Biobank (b). Data are presented as odds ratios with 95% confidence intervals (horizontal lines). Box sizes are proportional to inverse-variance weights. For each phenotypic subgroup, associations are ordered by their effect size. P values were obtained from the published GWAS. Associations that passed correction for multiple testing in this analysis (P = 0.05/13 traits = 3.85 × 10−3) are highlighted in red. The number of cases and controls for each association is shown in the forest plot. Supplementary Data 1 provides the association statistics for all traits, as well as data sources and references. c Forest plot showing the associations of rs867186-G with clinical…

• 

  Fig. 3. Statistical colocalization of cardiovascular outcomes and traits at the PROCR locus. a Regional association plots at the PROCR gene locus, showing the genetic association with coagulation factor VII, protein C, activated protein C, DVT and CAD. Unadjusted P values were obtained from the published GWAS. Details about the statistical analysis and source of the data are given in the Methods section. Color key indicates r2 with the respective lead variants in the GWAS. b Plot showing the colocalization posterior probabilities explained by each of the genetic variants at the chr20q11.22 locus tested in the colocalization analysis.

• 

  Fig. 4. Effect of rs867186 genotype on plasma biomarkers and EPCR expression on HUVECs. a Boxplots showing the distribution of plasma biomarker levels as a function of rs867186 genotype in up to 52 individuals. We measured plasma levels of sEPCR (n = 52 individuals across the three genotypic groups), APC (n = 52) and TAT complex (n = 51) using immunoassays, and PC levels (n = 51) using a chromogenic assay. All measurements were done with three technical replicates. The boxplots show the interquartile range in the box with the median as a horizontal line. Whiskers extend to ±1.5 times the interquartile range. Dashed lines indicate the fitted linear regression model for biomarker~genotype. P values for the additive regression model are indicated. b Boxplots showing the distribution of EPCR levels on HUVECs homozygous for the rs867186-G-allele or A-allele (n = 3 cell lines per genotypic group). Data show mean fluorescence intensity values of EPCR on untreated HUVECs, normalized to mean fluorescence intensity values of homozygotes of the rs867186-A-allele. The boxplot shows the interquartile range in the box with the median as a horizontal line. Whiskers extend to ±1.5 times the…

• 

  Fig. 5. Effect of APC on cell adhesion molecule expression and leukocyte–endothelial cell adhesion. a Barplots showing gene expression levels of ICAM1 and PROCR in human umbilical vein endothelial cells (HUVECs) and human coronary artery endothelial cells (HCAECs) relative to the control condition (i.e., 0 nM APC; indicated with a dashed line). Cells were co-incubated with 1 ng/ml TNF-α and varying concentrations of APC (0, 0.1, 1, 10, 100 nM) for 24 h. Data are shown for n = 5 (ICAM1) and n = 4 (PROCR) biological replicates in HUVECs and n = 4 (ICAM1) and n = 5 (PROCR) biological replicates in HCAECs. Error bars show standard deviations of the means. The blue and green lines indicate the fitted linear regression model for gene expression level~log(APC concentration). P values for the F-test of the linear regression model are shown. For each biological replicate, three technical replicates were averaged. b Barplots showing mean cell adhesion events using static adhesion assays with PMA-stimulated monocytic cells (U937) and TNF-α-activated HUVECs or HCAECs. Cells were co-incubated with 1 ng/ml TNF-α and 100 nM APC for 24 h (Methods). Data are shown for n = 5 and n = 4 biological…

• 

  Fig. 6. Proposed molecular mechanism underlying the PROCR-p.S219G variant. Credits: Icons were made by Pixel perfect from https://www.flaticon.com/.