Multi-ancestry polygenic risk scores for venous thromboembolism.

Venous thromboembolism (VTE) is a significant contributor to morbidity and mortality, with large disparities in incidence rates between Black and White Americans. Polygenic risk scores (PRSs) limited to variants discovered in genome-wide association studies in European-ancestry samples can identify European-ancestry individuals at high risk of VTE. However, there is limited evidence on whether high-dimensional PRS constructed using more sophisticated methods and more diverse training data can enhance the predictive ability and their utility across diverse populations. We developed PRSs for VTE using summary statistics from the International Network against Venous Thrombosis (INVENT) consortium genome-wide association studies meta-analyses of European- (71 771 cases and 1 059 740 controls) and African-ancestry samples (7482 cases and 129 975 controls). We used LDpred2 and PRS-CSx to construct ancestry-specific and multi-ancestry PRSs and evaluated their performance in an independent European- (6781 cases and 103 016 controls) and African-ancestry sample (1385 cases and 12 569 controls). Multi-ancestry PRSs with weights tuned in European-ancestry samples slightly outperformed ancestry-specific PRSs in European-ancestry test samples (e.g. the area under the receiver operating curve [AUC] was 0.609 for PRS-CSx_combinedEUR and 0.608 for PRS-CSxEUR [P = 0.00029]). Multi-ancestry PRSs with weights tuned in African-ancestry samples also outperformed ancestry-specific PRSs in African-ancestry test samples (PRS-CSxAFR: AUC = 0.58, PRS-CSx_combined AFR: AUC = 0.59), although this difference was not statistically significant (P = 0.34). The highest fifth percentile of the best-performing PRS was associated with 1.9-fold and 1.68-fold increased risk for VTE among European- and African-ancestry subjects, respectively, relative to those in the middle stratum. These findings suggest that the multi-ancestry PRS might be used to improve performance across diverse populations to identify individuals at highest risk for VTE.

Development, characterization, and replication of proteomic aging clocks: Analysis of 2 population-based cohorts.

BACKGROUND: Biological age may be estimated by proteomic aging clocks (PACs). Previous published PACs were constructed either in smaller studies or mainly in white individuals, and they used proteomic measures from only one-time point. In this study, we created de novo PACs and compared their performance to published PACs at 2 different time points in the Atherosclerosis Risk in Communities (ARIC) study of white and black participants (around 75% white and 25% black). MEDTHODS AND FINDINGS: A total of 4,712 plasma proteins were measured using SomaScan in blood samples collected in 1990 to 1992 from 11,761 midlife participants (aged 46 to 70 years) and in 2011 to 2013 from 5,183 late-life participants (aged 66 to 90 years). The de novo ARIC PACs were constructed by training them against chronological age using elastic net regression in two-thirds of healthy participants in midlife and late life and validated in the remaining one-third of healthy participants at the corresponding time point. We also computed 3 published PACs. We estimated age acceleration for each PAC as residuals after regressing each PAC on chronological age. We also calculated the change in age acceleration from midlife to late life. We examined the associations of age acceleration and change in age acceleration with mortality through 2019 from all-cause, cardiovascular disease (CVD), cancer, and lower respiratory disease (LRD) using Cox proportional hazards regression in participants (irrespective of health) after excluding the training set. The model was adjusted for chronological age, smoking, body mass index (BMI), and other confounders. We externally validated the midlife PAC using the Multi-Ethnic Study of Atherosclerosis (MESA) Exam 1 data. The ARIC PACs had a slightly stronger correlation with chronological age than published PACs in healthy participants at each time point. Associations with mortality were similar for the ARIC PACs and published PACs. For late-life and midlife age acceleration for the ARIC PACs, respectively, hazard ratios (HRs) per 1 standard deviation were 1.65 and 1.38 (both p < 0.001) for all-cause mortality, 1.37 and 1.20 (both p < 0.001) for CVD mortality, 1.21 (p = 0.028) and 1.04 (p = 0.280) for cancer mortality, and 1.68 and 1.36 (both p < 0.001) for LRD mortality. For the change in age acceleration, HRs for all-cause, CVD, and LRD mortality were comparable to the HRs for late-life age acceleration. The association between the change in age acceleration and cancer mortality was not significant. The external validation of the midlife PAC in MESA showed significant associations with mortality, as observed for midlife participants in ARIC. The main limitation is that our PACs were constructed in midlife and late-life participants. It is unknown whether these PACs could be applied to young individuals. CONCLUSIONS: In this longitudinal study, we found that the ARIC PACs and published PACs were similarly associated with an increased risk of mortality. These findings suggested that PACs show promise as biomarkers of biological age. PACs may be serve as tools to predict mortality and evaluate the effect of anti-aging lifestyle and therapeutic interventions.

Proteomic profiling identifies novel proteins for genetic risk of severe COVID-19: the Atherosclerosis Risk in Communities Study.

BACKGROUND: Genome-wide association studies have identified six genetic variants associated with severe COVID-19, yet the mechanisms through which they may affect disease remains unclear. We investigated proteomic signatures related to COVID-19 risk variants rs657152 (ABO), rs10735079 (OAS1/OAS2/OAS3), rs2109069 (DPP9), rs74956615 (TYK2), rs2236757 (IFNAR2) and rs11385942 (SLC6A20/LZTFL1/CCR9/FYCO1/CXCR6/XCR1) as well as their corresponding downstream pathways that may promote severe COVID-19 in risk allele carriers and their potential relevancies to other infection outcomes. METHODS: A DNA aptamer-based array measured 4870 plasma proteins among 11 471 participants. Linear regression estimated associations between the COVID-19 risk variants and proteins with correction for multiple comparisons, and canonical pathway analysis was conducted. Cox regression assessed associations between proteins identified in the main analysis and risk of incident hospitalized respiratory infections (2570 events) over a 20.7-year follow-up. RESULTS: The ABO variant rs657152 was associated with 84 proteins in 7241 white participants with 24 replicated in 1671 Black participants. The TYK2 variant rs74956615 was associated with ICAM-1 and -5 in white participants with ICAM-5 replicated in Black participants. Of the 84 proteins identified in the main analysis, seven were significantly associated with incident hospitalized respiratory infections including Ephrin type-A receptor 4 (hazard ratio (HR): 0.87; P = 2.3 × 10-11) and von Willebrand factor type A (HR: 1.17; P = 1.6x10-13). CONCLUSIONS: Novel proteomics signatures and pathways for COVID-19-related risk variants TYK2 and ABO were identified. A subset of these proteins predicted greater risk of incident hospitalized pneumonia and respiratory infections. Further studies to examine these proteins in COVID-19 patients are warranted.

Whole-exome sequencing of 14 389 individuals from the ESP and CHARGE consortia identifies novel rare variation associated with hemostatic factors.

Plasma levels of fibrinogen, coagulation factors VII and VIII and von Willebrand factor (vWF) are four intermediate phenotypes that are heritable and have been associated with the risk of clinical thrombotic events. To identify rare and low-frequency variants associated with these hemostatic factors, we conducted whole-exome sequencing in 10 860 individuals of European ancestry (EA) and 3529 African Americans (AAs) from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium and the National Heart, Lung and Blood Institute's Exome Sequencing Project. Gene-based tests demonstrated significant associations with rare variation (minor allele frequency < 5%) in fibrinogen gamma chain (FGG) (with fibrinogen, P = 9.1 × 10-13), coagulation factor VII (F7) (with factor VII, P = 1.3 × 10-72; seven novel variants) and VWF (with factor VIII and vWF; P = 3.2 × 10-14; one novel variant). These eight novel rare variant associations were independent of the known common variants at these loci and tended to have much larger effect sizes. In addition, one of the rare novel variants in F7 was significantly associated with an increased risk of venous thromboembolism in AAs (Ile200Ser; rs141219108; P = 4.2 × 10-5). After restricting gene-based analyses to only loss-of-function variants, a novel significant association was detected and replicated between factor VIII levels and a stop-gain mutation exclusive to AAs (rs3211938) in CD36 molecule (CD36). This variant has previously been linked to dyslipidemia but not with the levels of a hemostatic factor. These efforts represent the largest integration of whole-exome sequence data from two national projects to identify genetic variation associated with plasma hemostatic factors.

Novel Metabolites Are Associated With Thrombosis After Cardiac Surgery in Children With Congenital Heart Disease.

BACKGROUND: Thrombosis is a major complication after cardiac surgery in children with congenital heart disease. The mechanisms underlying thrombosis development remain poorly understood. We aimed to identify novel circulating metabolites before cardiac surgery that are associated with thrombosis after surgery in children with congenital heart disease. METHODS: In this prospective cohort study, all blood samples were drawn right before surgical incision and after the induction of anesthesia, and plasma was separated immediately under 4 °C. Untargeted metabolomic data were measured by Metabolon in plasma from children (age range, 0 days-18 years) with congenital heart disease undergoing cardiac surgery. The primary outcome was thrombosis within 30 days of surgery or before discharge. Associations of individual metabolites with thrombosis were assessed with logistic regression with false discovery rate correction for multiple comparison and adjustment for clinical characteristics; elastic net regression was used to select a prediction model. RESULTS: Out of 1115 metabolites measured in samples from 203 children, 776 met the quality control criteria. In total, 25 children (12.3%) developed thrombosis. Among the 776 metabolites, 175 were significantly associated with thrombosis (false discovery rate Q<0.05). The top 3 metabolites showing the strongest associations with thrombosis were eicosapentaenoate, stearidonate, and andro steroid monosulfate C19H28O6S (false discovery rate, 0.01 for all). Pathway analysis showed that the pathways of nicotinate and nicotinamide metabolism and glycerophospholipid metabolism were enriched (false discovery rate, 0.003 for both) and had significant impact on the development of thrombosis. In elastic net regression analysis, the area under the receiver operating-characteristic curve of a prediction model for thrombosis was 0.969 in the training sample (70% of the total sample) and 0.833 in the testing sample (the remaining 30%). CONCLUSIONS: We have identified promising novel metabolites and metabolic pathways associated with thrombosis. Future studies are warranted to confirm these findings and examine the mechanistic pathways to thrombosis.

Proteomics Analysis of Genetic Liability of Abdominal Aortic Aneurysm Identifies Plasma Neogenin and Kit Ligand: The ARIC Study.

BACKGROUND: Genome-wide association studies have reported 23 gene loci related to abdominal aortic aneurysm (AAA)-a potentially lethal condition characterized by a weakened dilated vessel wall. This study aimed to identify proteomic signatures and pathways related to these risk loci to better characterize AAA genetic susceptibility. METHODS: Plasma concentrations of 4870 proteins were determined using a DNA aptamer-based array. Linear regression analysis estimated the associations between the 23 risk alleles and plasma protein levels with adjustments for potential confounders in a race-stratified analysis of 1671 Black and 7241 White participants. Significant proteins were then evaluated for their prediction of clinical AAA (454 AAA events in 11 064 individuals), and those significantly associated with AAA were further interrogated using Mendelian randomization analysis. RESULTS: Risk variants proximal to PSRC1-CELSR2-SORT1, PCIF1-ZNF335-MMP9, RP11-136O12.2/TRIB1, ZNF259/APOA5, IL6R, PCSK9, LPA, and APOE were associated with 118 plasma proteins in Whites and 59 were replicated in Black participants. Novel associations with clinical AAA incidence were observed for kit ligand (HR, 0.59 [95% CI, 0.42-0.82] for top versus first quintiles) and neogenin (HR, 0.64 [95% CI, 0.46-0.88]) over a median 21.2-year follow-up; neogenin was also associated with ultrasound-detected asymptomatic AAA (N=4295; 57 asymptomatic AAA cases). Mendelian randomization inverse variance weighted estimates suggested that AAA risk is promoted by lower levels of kit ligand (OR per SD=0.67; P=1.4×10-5) and neogenin (OR per SD=0.50; P=0.03). CONCLUSIONS: Low levels of neogenin and kit ligand may be novel risk factors for AAA development in potentially causal pathways. These findings provide insights and potential targets to reduce AAA susceptibility.

Antithrombin, Protein C, and Protein S: Genome and Transcriptome-Wide Association Studies Identify 7 Novel Loci Regulating Plasma Levels.

BACKGROUND: Antithrombin, PC (protein C), and PS (protein S) are circulating natural anticoagulant proteins that regulate hemostasis and of which partial deficiencies are causes of venous thromboembolism. Previous genetic association studies involving antithrombin, PC, and PS were limited by modest sample sizes or by being restricted to candidate genes. In the setting of the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium, we meta-analyzed across ancestries the results from 10 genome-wide association studies of plasma levels of antithrombin, PC, PS free, and PS total. METHODS: Study participants were of European and African ancestries, and genotype data were imputed to TOPMed, a dense multiancestry reference panel. Each of the 10 studies conducted a genome-wide association studies for each phenotype and summary results were meta-analyzed, stratified by ancestry. Analysis of antithrombin included 25 243 European ancestry and 2688 African ancestry participants, PC analysis included 16 597 European ancestry and 2688 African ancestry participants, PSF and PST analysis included 4113 and 6409 European ancestry participants. We also conducted transcriptome-wide association analyses and multiphenotype analysis to discover additional associations. Novel genome-wide association studies and transcriptome-wide association analyses findings were validated by in vitro functional experiments. Mendelian randomization was performed to assess the causal relationship between these proteins and cardiovascular outcomes. RESULTS: Genome-wide association studies meta-analyses identified 4 newly associated loci: 3 with antithrombin levels (GCKR, BAZ1B, and HP-TXNL4B) and 1 with PS levels (ORM1-ORM2). transcriptome-wide association analyses identified 3 newly associated genes: 1 with antithrombin level (FCGRT), 1 with PC (GOLM2), and 1 with PS (MYL7). In addition, we replicated 7 independent loci reported in previous studies. Functional experiments provided evidence for the involvement of GCKR, SNX17, and HP genes in antithrombin regulation. CONCLUSIONS: The use of larger sample sizes, diverse populations, and a denser imputation reference panel allowed the detection of 7 novel genomic loci associated with plasma antithrombin, PC, and PS levels.

The relationship between serum manganese concentration with all-cause and cause-specific mortality: a retrospective and population-based cross-sectional study.

BACKGROUND: The study aimed to explore the association between manganese concentration and all-cause, cardiovascular disease (CVD)-related, and cancer-related mortality in the general population of the United States. METHODS: We integrated the data from the National Health and Nutrition Examination Survey from 2011 to 2018. A total of 9,207 subjects were selected based on the inclusion and exclusion criteria. The relationship between manganese concentration and all-cause, CVD-related, and cancer-related mortality was analyzed by constructing a Cox proportional hazard regression model and a restricted cubic spline (RCS) plot. Additionally, subgroup analyses stratified by age, sex, race/ethnicity, hypertension, diabetes mellitus (DM), chronic heart disease, chronic heart failure, angina pectoris, heart attack, stroke, and BMI were further performed. RESULTS: In the full adjusted model, compared with the lowest quartile, the adjusted hazard ratios with 95% confidence intervals (CIs) for all-cause, CVD-related, and cancer-related mortality across manganese quartiles were (1.11 (0.87,1.41), 0.96 (0.74, 1.23), and 1.23 (0.96, 1.59); P-value for trend =0.041), (0.86 (0.54, 1.37), 0.87 (0.55, 1.40), and 1.07 (0.67, 1.72); P-value for trend =0.906), and (1.45 (0.92, 2.29), 1.14 (0.70, 1.88), and 1.26 (0.75, 2.11); P-value for trend =0.526), respectively. The RCS curve shown a U-shaped association between manganese concentration and all-cause mortality and CVD-related mortality (P-value for nonlinear <0.05). However, there was an increase and then a decrease in the link between manganese concentration and cancer-related mortality (P-value for nonlinear <0.05). Manganese exposure was positively correlated with sex (correlation coefficient, r =0.19, P-value <0.001) and negatively correlated with age (correlation coefficient, r =-0.11, P-value <0.001) and serum creatinine (correlation coefficient, r =-0.12, P-value <0.001), respectively. CONCLUSIONS: Our findings suggest that elevated serum manganese concentrations are associated with all-cause and CVD-related mortality in the U.S. population and that maintenance of serum manganese between 8.67-9.23 µg/L may promote public health.

Proteomic analysis of diabetes genetic risk scores identifies complement C2 and neuropilin-2 as predictors of type 2 diabetes: the Atherosclerosis Risk in Communities (ARIC) Study.

AIMS/HYPOTHESIS: Genetic predisposition to type 2 diabetes is well-established, and genetic risk scores (GRS) have been developed that capture heritable liabilities for type 2 diabetes phenotypes. However, the proteins through which these genetic variants influence risk have not been thoroughly investigated. This study aimed to identify proteins and pathways through which type 2 diabetes risk variants may influence pathophysiology. METHODS: Using a proteomics data-driven approach in a discovery sample of 7241 White participants in the Atherosclerosis Risk in Communities Study (ARIC) cohort and a replication sample of 1674 Black ARIC participants, we interrogated plasma levels of 4870 proteins and four GRS of specific type 2 diabetes phenotypes related to beta cell function, insulin resistance, lipodystrophy, BMI/blood lipid abnormalities and a composite score of all variants combined. RESULTS: Twenty-two plasma proteins were identified in White participants after Bonferroni correction. Of the 22 protein-GRS associations that were statistically significant, 10 were replicated in Black participants and all but one were directionally consistent. In a secondary analysis, 18 of the 22 proteins were found to be associated with prevalent type 2 diabetes and ten proteins were associated with incident type 2 diabetes. Two-sample Mendelian randomisation indicated that complement C2 may be causally related to greater type 2 diabetes risk (inverse variance weighted estimate: OR 1.65 per SD; p=7.0 × 10-3), while neuropilin-2 was inversely associated (OR 0.44 per SD; p=8.0 × 10-3). CONCLUSIONS/INTERPRETATION: Identified proteins may represent viable intervention or pharmacological targets to prevent, reverse or slow type 2 diabetes progression, and further research is needed to pursue these targets.

Comparison of Proteomic Measurements Across Platforms in the Atherosclerosis Risk in Communities (ARIC) Study.

BACKGROUND: The plasma proteome can be quantified using different types of highly multiplexed technologies, including aptamer-based and proximity-extension immunoassay methods. There has been limited characterization of how these protein measurements correlate across platforms and with absolute measures from targeted immunoassays. METHODS: We assessed the comparability of (a) highly multiplexed aptamer-based (SomaScan v4; Somalogic) and proximity-extension immunoassay (OLINK Proseek® v5003; Olink) methods in 427 Atherosclerosis Risk in Communities (ARIC) Study participants (Visit 5, 2011-2013), and (b) 18 of the SomaScan protein measurements against targeted immunoassays in 110 participants (55 cardiovascular disease cases, 55 controls). We calculated Spearman correlations (r) between the different measurements and compared associations with case-control status. RESULTS: There were 417 protein comparisons (366 unique proteins) between the SomaScan and Olink platforms. The average correlation was r = 0.46 (range: -0.21 to 0.97; 79 [19%] with r ≥ 0.8). For the comparison of SomaScan and targeted immunoassays, 6 of 18 assays (growth differentiation factor 15 [GDF15], interleukin-1 receptor-like 1 [ST2], interstitial collagenase [MMP1], adiponectin, leptin, and resistin) had good correlations (r ≥ 0.8), 2 had modest correlations (0.5 ≤ r < 0.8; osteopontin and interleukin-6 [IL6]), and 10 were poorly correlated (r < 0.5; metalloproteinase inhibitor 1 [TIMP1], stromelysin-1 [MMP3], matrilysin [MMP7], C-C motif chemokine 2 [MCP1], interleukin-10 [IL10], vascular cell adhesion protein 1 [VCAM1], intercellular adhesion molecule 1 [ICAM1], interleukin-18 [IL18], tumor necrosis factor [TNFα], and visfatin) overall. Correlations for SomaScan and targeted immunoassays were similar according to case status. CONCLUSIONS: There is variation in the quantitative measurements for many proteins across aptamer-based and proximity-extension immunoassays (approximately 1/2 showing good or modest correlation and approximately 1/2 poor correlation) and also for correlations of these highly multiplexed technologies with targeted immunoassays. Design and interpretation of protein quantification studies should be informed by the variation across measurement techniques for each protein.

Sequencing Analysis at 8p23 Identifies Multiple Rare Variants in DLC1 Associated with Sleep-Related Oxyhemoglobin Saturation Level.

Average arterial oxyhemoglobin saturation during sleep (AvSpO2S) is a clinically relevant measure of physiological stress associated with sleep-disordered breathing, and this measure predicts incident cardiovascular disease and mortality. Using high-depth whole-genome sequencing data from the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) project and focusing on genes with linkage evidence on chromosome 8p23,1,2 we observed that six coding and 51 noncoding variants in a gene that encodes the GTPase-activating protein (DLC1) are significantly associated with AvSpO2S and replicated in independent subjects. The combined DLC1 association evidence of discovery and replication cohorts reaches genome-wide significance in European Americans (p = 7.9 × 10-7). A risk score for these variants, built on an independent dataset, explains 0.97% of the AvSpO2S variation and contributes to the linkage evidence. The 51 noncoding variants are enriched in regulatory features in a human lung fibroblast cell line and contribute to DLC1 expression variation. Mendelian randomization analysis using these variants indicates a significant causal effect of DLC1 expression in fibroblasts on AvSpO2S. Multiple sources of information, including genetic variants, gene expression, and methylation, consistently suggest that DLC1 is a gene associated with AvSpO2S.

A Mendelian randomization of γ' and total fibrinogen levels in relation to venous thromboembolism and ischemic stroke.

Fibrinogen is a key component of the coagulation cascade, and variation in its circulating levels may contribute to thrombotic diseases, such as venous thromboembolism (VTE) and ischemic stroke. Gamma prime (γ') fibrinogen is an isoform of fibrinogen that has anticoagulant properties. We applied 2-sample Mendelian randomization (MR) to estimate the causal effect of total circulating fibrinogen and its isoform, γ' fibrinogen, on risk of VTE and ischemic stroke subtypes using summary statistics from genome-wide association studies. Genetic instruments for γ' fibrinogen and total fibrinogen were selected, and the inverse-variance weighted MR approach was used to estimate causal effects in the main analysis, complemented by sensitivity analyses that are more robust to the inclusion of pleiotropic variants, including MR-Egger, weighted median MR, and weighted mode MR. The main inverse-variance weighted MR estimates based on a combination of 16 genetic instruments for γ' fibrinogen and 75 genetic instruments for total fibrinogen indicated a protective effect of higher γ' fibrinogen and higher total fibrinogen on VTE risk. There was also a protective effect of higher γ' fibrinogen levels on cardioembolic and large artery stroke risk. Effect estimates were consistent across sensitivity analyses. Our results provide evidence to support effects of genetically determined γ' fibrinogen on VTE and ischemic stroke risk. Further research is needed to explore mechanisms underlying these effects and their clinical applications.

Association of Differential Leukocyte Count With Incident Abdominal Aortic Aneurysm Over 22.5 Years: The ARIC Study.

OBJECTIVE: Leukocytes contribute to the development of abdominal aortic aneurysm (AAA). We evaluated whether associations of differential leukocyte counts with AAA persist after accounting for traditional risk factors of AAA. APPROACH AND RESULTS: Among 11 217 adults from the Atherosclerosis Risk in Communities Study, we evaluated associations of differential leukocyte counts at baseline (1987­1989) with incident AAAs over a median follow-up of 22.5 years, using Cox proportional hazards regression. Each differential leukocyte count was categorized into 5 groups­below normal, tertiles within the normal range, and above normal, with the first tertile serving as the referent. We identified 377 incident AAAs through 2011, using hospital discharge diagnoses, linked Medicare records, or death certificates. At baseline, higher neutrophil, monocyte, and eosinophil counts were associated with higher risk of AAA, independent of smoking, other differential leukocyte counts, and other traditional risk factors. The association with incident AAA was the strongest for above normal neutrophil count, with an adjusted hazard ratio (95% CI) of 2.17 (1.29­3.64). Below normal neutrophil, lymphocyte, eosinophil and basophil counts were associated with higher risk of AAA with adjusted hazard ratio (95% CI) between 1.86 (1.04­3.35) and 1.62 (1.10­2.39). CONCLUSIONS: Higher neutrophil, monocyte, and eosinophil counts in midlife are associated with higher risk of AAA, even after accounting for traditional risk factors such as smoking, obesity, and atherosclerosis. This suggests the need to identify nontraditional risk factors and treatment strategies to mitigate the residual risk of AAA conferred by midlife inflammation. Whether immunosuppression is associated with higher risk of AAA needs further investigation.

Multi-ethnic genome-wide association analyses of white blood cell and platelet traits in the Population Architecture using Genomics and Epidemiology (PAGE) study.

BACKGROUND: Circulating white blood cell and platelet traits are clinically linked to various disease outcomes and differ across individuals and ancestry groups. Genetic factors play an important role in determining these traits and many loci have been identified. However, most of these findings were identified in populations of European ancestry (EA), with African Americans (AA), Hispanics/Latinos (HL), and other races/ethnicities being severely underrepresented. RESULTS: We performed ancestry-combined and ancestry-specific genome-wide association studies (GWAS) for white blood cell and platelet traits in the ancestrally diverse Population Architecture using Genomics and Epidemiology (PAGE) Study, including 16,201 AA, 21,347 HL, and 27,236 EA participants. We identified six novel findings at suggestive significance (P < 5E-8), which need confirmation, and independent signals at six previously established regions at genome-wide significance (P < 2E-9). We confirmed multiple previously reported genome-wide significant variants in the single variant association analysis and multiple genes using PrediXcan. Evaluation of loci reported from a Euro-centric GWAS indicated attenuation of effect estimates in AA and HL compared to EA populations. CONCLUSIONS: Our results highlighted the potential to identify ancestry-specific and ancestry-agnostic variants in participants with diverse backgrounds and advocate for continued efforts in improving inclusion of racially/ethnically diverse populations in genetic association studies for complex traits.

A genome-wide association study identifies new loci for factor VII and implicates factor VII in ischemic stroke etiology.

Factor VII (FVII) is an important component of the coagulation cascade. Few genetic loci regulating FVII activity and/or levels have been discovered to date. We conducted a meta-analysis of 9 genome-wide association studies of plasma FVII levels (7 FVII activity and 2 FVII antigen) among 27 495 participants of European and African ancestry. Each study performed ancestry-specific association analyses. Inverse variance weighted meta-analysis was performed within each ancestry group and then combined for a trans-ancestry meta-analysis. Our primary analysis included the 7 studies that measured FVII activity, and a secondary analysis included all 9 studies. We provided functional genomic validation for newly identified significant loci by silencing candidate genes in a human liver cell line (HuH7) using small-interfering RNA and then measuring F7 messenger RNA and FVII protein expression. Lastly, we used meta-analysis results to perform Mendelian randomization analysis to estimate the causal effect of FVII activity on coronary artery disease, ischemic stroke (IS), and venous thromboembolism. We identified 2 novel (REEP3 and JAZF1-AS1) and 6 known loci associated with FVII activity, explaining 19.0% of the phenotypic variance. Adding FVII antigen data to the meta-analysis did not result in the discovery of further loci. Silencing REEP3 in HuH7 cells upregulated FVII, whereas silencing JAZF1 downregulated FVII. Mendelian randomization analyses suggest that FVII activity has a positive causal effect on the risk of IS. Variants at REEP3 and JAZF1 contribute to FVII activity by regulating F7 expression levels. FVII activity appears to contribute to the etiology of IS in the general population.

Genomic and transcriptomic association studies identify 16 novel susceptibility loci for venous thromboembolism.

Venous thromboembolism (VTE) is a significant contributor to morbidity and mortality. To advance our understanding of the biology contributing to VTE, we conducted a genome-wide association study (GWAS) of VTE and a transcriptome-wide association study (TWAS) based on imputed gene expression from whole blood and liver. We meta-analyzed GWAS data from 18 studies for 30 234 VTE cases and 172 122 controls and assessed the association between 12 923 718 genetic variants and VTE. We generated variant prediction scores of gene expression from whole blood and liver tissue and assessed them for association with VTE. Mendelian randomization analyses were conducted for traits genetically associated with novel VTE loci. We identified 34 independent genetic signals for VTE risk from GWAS meta-analysis, of which 14 are newly reported associations. This included 11 newly associated genetic loci (C1orf198, PLEK, OSMR-AS1, NUGGC/SCARA5, GRK5, MPHOSPH9, ARID4A, PLCG2, SMG6, EIF5A, and STX10) of which 6 replicated, and 3 new independent signals in 3 known genes. Further, TWAS identified 5 additional genetic loci with imputed gene expression levels differing between cases and controls in whole blood (SH2B3, SPSB1, RP11-747H7.3, RP4-737E23.2) and in liver (ERAP1). At some GWAS loci, we found suggestive evidence that the VTE association signal for novel and previously known regions colocalized with expression quantitative trait locus signals. Mendelian randomization analyses suggested that blood traits may contribute to the underlying risk of VTE. To conclude, we identified 16 novel susceptibility loci for VTE; for some loci, the association signals are likely mediated through gene expression of nearby genes.

A large-scale exome array analysis of venous thromboembolism.

Although recent Genome-Wide Association Studies have identified novel associations for common variants, there has been no comprehensive exome-wide search for low-frequency variants that affect the risk of venous thromboembolism (VTE). We conducted a meta-analysis of 11 studies comprising 8,332 cases and 16,087 controls of European ancestry and 382 cases and 1,476 controls of African American ancestry genotyped with the Illumina HumanExome BeadChip. We used the seqMeta package in R to conduct single variant and gene-based rare variant tests. In the single variant analysis, we limited our analysis to the 64,794 variants with at least 40 minor alleles across studies (minor allele frequency [MAF] ~0.08%). We confirmed associations with previously identified VTE loci, including ABO, F5, F11, and FGA. After adjusting for multiple testing, we observed no novel significant findings in single variant or gene-based analysis. Given our sample size, we had greater than 80% power to detect minimum odds ratios greater than 1.5 and 1.8 for a single variant with MAF of 0.01 and 0.005, respectively. Larger studies and sequence data may be needed to identify novel low-frequency and rare variants associated with VTE risk.

Genome-Wide Association Transethnic Meta-Analyses Identifies Novel Associations Regulating Coagulation Factor VIII and von Willebrand Factor Plasma Levels.

BACKGROUND: Factor VIII (FVIII) and its carrier protein von Willebrand factor (VWF) are associated with risk of arterial and venous thrombosis and with hemorrhagic disorders. We aimed to identify and functionally test novel genetic associations regulating plasma FVIII and VWF. METHODS: We meta-analyzed genome-wide association results from 46 354 individuals of European, African, East Asian, and Hispanic ancestry. All studies performed linear regression analysis using an additive genetic model and associated ≈35 million imputed variants with natural log-transformed phenotype levels. In vitro gene silencing in cultured endothelial cells was performed for candidate genes to provide additional evidence on association and function. Two-sample Mendelian randomization analyses were applied to test the causal role of FVIII and VWF plasma levels on the risk of arterial and venous thrombotic events. RESULTS: We identified 13 novel genome-wide significant ( P≤2.5×10-8) associations, 7 with FVIII levels ( FCHO2/TMEM171/TNPO1, HLA, SOX17/RP1, LINC00583/NFIB, RAB5C-KAT2A, RPL3/TAB1/SYNGR1, and ARSA) and 11 with VWF levels ( PDHB/PXK/KCTD6, SLC39A8, FCHO2/TMEM171/TNPO1, HLA, GIMAP7/GIMAP4, OR13C5/NIPSNAP, DAB2IP, C2CD4B, RAB5C-KAT2A, TAB1/SYNGR1, and ARSA), beyond 10 previously reported associations with these phenotypes. Functional validation provided further evidence of association for all loci on VWF except ARSA and DAB2IP. Mendelian randomization suggested causal effects of plasma FVIII activity levels on venous thrombosis and coronary artery disease risk and plasma VWF levels on ischemic stroke risk. CONCLUSIONS: The meta-analysis identified 13 novel genetic loci regulating FVIII and VWF plasma levels, 10 of which we validated functionally. We provide some evidence for a causal role of these proteins in thrombotic events.

Ancestry-specific associations identified in genome-wide combined-phenotype study of red blood cell traits emphasize benefits of diversity in genomics.

BACKGROUND: Quantitative red blood cell (RBC) traits are highly polygenic clinically relevant traits, with approximately 500 reported GWAS loci. The majority of RBC trait GWAS have been performed in European- or East Asian-ancestry populations, despite evidence that rare or ancestry-specific variation contributes substantially to RBC trait heritability. Recently developed combined-phenotype methods which leverage genetic trait correlation to improve statistical power have not yet been applied to these traits. Here we leveraged correlation of seven quantitative RBC traits in performing a combined-phenotype analysis in a multi-ethnic study population. RESULTS: We used the adaptive sum of powered scores (aSPU) test to assess combined-phenotype associations between ~ 21 million SNPs and seven RBC traits in a multi-ethnic population (maximum n = 67,885 participants; 24% African American, 30% Hispanic/Latino, and 43% European American; 76% female). Thirty-nine loci in our multi-ethnic population contained at least one significant association signal (p < 5E-9), with lead SNPs at nine loci significantly associated with three or more RBC traits. A majority of the lead SNPs were common (MAF > 5%) across all ancestral populations. Nineteen additional independent association signals were identified at seven known loci (HFE, KIT, HBS1L/MYB, CITED2/FILNC1, ABO, HBA1/2, and PLIN4/5). For example, the HBA1/2 locus contained 14 conditionally independent association signals, 11 of which were previously unreported and are specific to African and Amerindian ancestries. One variant in this region was common in all ancestries, but exhibited a narrower LD block in African Americans than European Americans or Hispanics/Latinos. GTEx eQTL analysis of all independent lead SNPs yielded 31 significant associations in relevant tissues, over half of which were not at the gene immediately proximal to the lead SNP. CONCLUSION: This work identified seven loci containing multiple independent association signals for RBC traits using a combined-phenotype approach, which may improve discovery in genetically correlated traits. Highly complex genetic architecture at the HBA1/2 locus was only revealed by the inclusion of African Americans and Hispanics/Latinos, underscoring the continued importance of expanding large GWAS to include ancestrally diverse populations.

Discovery, fine-mapping, and conditional analyses of genetic variants associated with C-reactive protein in multiethnic populations using the Metabochip in the Population Architecture using Genomics and Epidemiology (PAGE) study.

C-reactive protein (CRP) is a circulating biomarker indicative of systemic inflammation. We aimed to evaluate genetic associations with CRP levels among non-European-ancestry populations through discovery, fine-mapping and conditional analyses. A total of 30 503 non-European-ancestry participants from 6 studies participating in the Population Architecture using Genomics and Epidemiology study had serum high-sensitivity CRP measurements and ∼200 000 single nucleotide polymorphisms (SNPs) genotyped on the Metabochip. We evaluated the association between each SNP and log-transformed CRP levels using multivariate linear regression, with additive genetic models adjusted for age, sex, the first four principal components of genetic ancestry, and study-specific factors. Differential linkage disequilibrium patterns between race/ethnicity groups were used to fine-map regions associated with CRP levels. Conditional analyses evaluated for multiple independent signals within genetic regions. One hundred and sixty-three unique variants in 12 loci in overall or race/ethnicity-stratified Metabochip-wide scans reached a Bonferroni-corrected P-value <2.5E-7. Three loci have no (HACL1, OLFML2B) or only limited (PLA2G6) previous associations with CRP levels. Six loci had different top hits in race/ethnicity-specific versus overall analyses. Fine-mapping refined the signal in six loci, particularly in HNF1A. Conditional analyses provided evidence for secondary signals in LEPR, IL1RN and HNF1A, and for multiple independent signals in CRP and APOE. We identified novel variants and loci associated with CRP levels, generalized known CRP associations to a multiethnic study population, refined association signals at several loci and found evidence for multiple independent signals at several well-known loci. This study demonstrates the benefit of conducting inclusive genetic association studies in large multiethnic populations.