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.

Whole genome sequencing based analysis of inflammation biomarkers in the Trans-Omics for Precision Medicine (TOPMed) consortium.

Inflammation biomarkers can provide valuable insight into the role of inflammatory processes in many diseases and conditions. Sequencing based analyses of such biomarkers can also serve as an exemplar of the genetic architecture of quantitative traits. To evaluate the biological insight, which can be provided by a multi-ancestry, whole-genome based association study, we performed a comprehensive analysis of 21 inflammation biomarkers from up to 38 465 individuals with whole-genome sequencing from the Trans-Omics for Precision Medicine (TOPMed) program (with varying sample size by trait, where the minimum sample size was n = 737 for MMP-1). We identified 22 distinct single-variant associations across 6 traits-E-selectin, intercellular adhesion molecule 1, interleukin-6, lipoprotein-associated phospholipase A2 activity and mass, and P-selectin-that remained significant after conditioning on previously identified associations for these inflammatory biomarkers. We further expanded upon known biomarker associations by pairing the single-variant analysis with a rare variant set-based analysis that further identified 19 significant rare variant set-based associations with 5 traits. These signals were distinct from both significant single variant association signals within TOPMed and genetic signals observed in prior studies, demonstrating the complementary value of performing both single and rare variant analyses when analyzing quantitative traits. We also confirm several previously reported signals from semi-quantitative proteomics platforms. Many of these signals demonstrate the extensive allelic heterogeneity and ancestry-differentiated variant-trait associations common for inflammation biomarkers, a characteristic we hypothesize will be increasingly observed with well-powered, large-scale analyses of complex traits.

A genetic association study of circulating coagulation factor VIII and von Willebrand factor levels.

ABSTRACT: Coagulation factor VIII (FVIII) and its carrier protein von Willebrand factor (VWF) are critical to coagulation and platelet aggregation. We leveraged whole-genome sequence data from the Trans-Omics for Precision Medicine (TOPMed) program along with TOPMed-based imputation of genotypes in additional samples to identify genetic associations with circulating FVIII and VWF levels in a single-variant meta-analysis, including up to 45 289 participants. Gene-based aggregate tests were implemented in TOPMed. We identified 3 candidate causal genes and tested their functional effect on FVIII release from human liver endothelial cells (HLECs) and VWF release from human umbilical vein endothelial cells. Mendelian randomization was also performed to provide evidence for causal associations of FVIII and VWF with thrombotic outcomes. We identified associations (P < 5 × 10-9) at 7 new loci for FVIII (ST3GAL4, CLEC4M, B3GNT2, ASGR1, F12, KNG1, and TREM1/NCR2) and 1 for VWF (B3GNT2). VWF, ABO, and STAB2 were associated with FVIII and VWF in gene-based analyses. Multiphenotype analysis of FVIII and VWF identified another 3 new loci, including PDIA3. Silencing of B3GNT2 and the previously reported CD36 gene decreased release of FVIII by HLECs, whereas silencing of B3GNT2, CD36, and PDIA3 decreased release of VWF by HVECs. Mendelian randomization supports causal association of higher FVIII and VWF with increased risk of thrombotic outcomes. Seven new loci were identified for FVIII and 1 for VWF, with evidence supporting causal associations of FVIII and VWF with thrombotic outcomes. B3GNT2, CD36, and PDIA3 modulate the release of FVIII and/or VWF in vitro.

Predicted leukocyte telomere length and risk of myeloid neoplasms.

Maintenance of telomere length has long been established to play a role in the biology of cancer and several studies suggest that it may be especially important in myeloid malignancies. To overcome potential bias in confounding and reverse causation of observational studies, we use both a polygenic risk score (PRS) and inverse-variance weighted (IVW) Mendelian randomization (MR) analyses to estimate the relationship between genetically predicted leukocyte telomere length (LTL) and acute myeloid leukemia (AML) risk in 498 cases and 2099 controls and myelodysplastic syndrome (MDS) risk in 610 cases and 1759 controls. Genetic instruments derived from four recent studies explaining 1.23-4.57% of telomere variability were considered. We used multivariable logistic regression to estimate odds ratios (OR, 95% confidence intervals [CI]) as the measure of association between individual single-nucleotide polymorphisms and myeloid malignancies. We observed a significant association between a PRS of longer predicted LTL and AML using three genetic instruments (OR = 4.03 per ~1200 base pair [bp] increase in LTL, 95% CI: 1.65, 9.85 using Codd et al. [Codd, V., Nelson, C.P., Albrecht, E., Mangino, M., Deelen, J., Buxton, J.L., Hottenga, J.J., Fischer, K., Esko, T., Surakka, I. et al. (2013) Identification of seven loci affecting mean telomere length and their association with disease. Nat. Genet., 45, 422-427 427e421-422.], OR = 3.48 per one-standard deviation increase in LTL, 95% CI: 1.74, 6.97 using Li et al. [Li, C., Stoma, S., Lotta, L.A., Warner, S., Albrecht, E., Allione, A., Arp, P.P., Broer, L., Buxton, J.L., Alves, A.D.S.C. et al. (2020) Genome-wide association analysis in humans links nucleotide metabolism to leukocyte telomere length. Am. J. Hum. Genet., 106, 389-404.] and OR = 2.59 per 1000 bp increase in LTL, 95% CI: 1.03, 6.52 using Taub et al. [Taub, M.A., Conomos, M.P., Keener, R., Iyer, K.R., Weinstock, J.S., Yanek, L.R., Lane, J., Miller-Fleming, T.W., Brody, J.A., Raffield, L.M. et al. (2022) Genetic determinants of telomere length from 109,122 ancestrally diverse whole-genome sequences in TOPMed. Cell Genom., 2.] genetic instruments). MR analyses further indicated an association between LTL and AML risk (PIVW ≤ 0.049) but not MDS (all PIVW ≥ 0.076). Findings suggest variation in genes relevant to telomere function and maintenance may be important in the etiology of AML but not MDS.

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.

Whole genome sequence analysis of platelet traits in the NHLBI Trans-Omics for Precision Medicine (TOPMed) initiative.

Platelets play a key role in thrombosis and hemostasis. Platelet count (PLT) and mean platelet volume (MPV) are highly heritable quantitative traits, with hundreds of genetic signals previously identified, mostly in European ancestry populations. We here utilize whole genome sequencing (WGS) from NHLBI's Trans-Omics for Precision Medicine initiative (TOPMed) in a large multi-ethnic sample to further explore common and rare variation contributing to PLT (n = 61 200) and MPV (n = 23 485). We identified and replicated secondary signals at MPL (rs532784633) and PECAM1 (rs73345162), both more common in African ancestry populations. We also observed rare variation in Mendelian platelet-related disorder genes influencing variation in platelet traits in TOPMed cohorts (not enriched for blood disorders). For example, association of GP9 with lower PLT and higher MPV was partly driven by a pathogenic Bernard-Soulier syndrome variant (rs5030764, p.Asn61Ser), and the signals at TUBB1 and CD36 were partly driven by loss of function variants not annotated as pathogenic in ClinVar (rs199948010 and rs571975065). However, residual signal remained for these gene-based signals after adjusting for lead variants, suggesting that additional variants in Mendelian genes with impacts in general population cohorts remain to be identified. Gene-based signals were also identified at several genome-wide association study identified loci for genes not annotated for Mendelian platelet disorders (PTPRH, TET2, CHEK2), with somatic variation driving the result at TET2. These results highlight the value of WGS in populations of diverse genetic ancestry to identify novel regulatory and coding signals, even for well-studied traits like platelet traits.

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.

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.

Cross-Ancestry Investigation of Venous Thromboembolism Genomic Predictors.

BACKGROUND: Venous thromboembolism (VTE) is a life-threatening vascular event with environmental and genetic determinants. Recent VTE genome-wide association studies (GWAS) meta-analyses involved nearly 30 000 VTE cases and identified up to 40 genetic loci associated with VTE risk, including loci not previously suspected to play a role in hemostasis. The aim of our research was to expand discovery of new genetic loci associated with VTE by using cross-ancestry genomic resources. METHODS: We present new cross-ancestry meta-analyzed GWAS results involving up to 81 669 VTE cases from 30 studies, with replication of novel loci in independent populations and loci characterization through in silico genomic interrogations. RESULTS: In our genetic discovery effort that included 55 330 participants with VTE (47 822 European, 6320 African, and 1188 Hispanic ancestry), we identified 48 novel associations, of which 34 were replicated after correction for multiple testing. In our combined discovery-replication analysis (81 669 VTE participants) and ancestry-stratified meta-analyses (European, African, and Hispanic), we identified another 44 novel associations, which are new candidate VTE-associated loci requiring replication. In total, across all GWAS meta-analyses, we identified 135 independent genomic loci significantly associated with VTE risk. A genetic risk score of the significantly associated loci in Europeans identified a 6-fold increase in risk for those in the top 1% of scores compared with those with average scores. We also identified 31 novel transcript associations in transcriptome-wide association studies and 8 novel candidate genes with protein quantitative-trait locus Mendelian randomization analyses. In silico interrogations of hemostasis and hematology traits and a large phenome-wide association analysis of the 135 GWAS loci provided insights to biological pathways contributing to VTE, with some loci contributing to VTE through well-characterized coagulation pathways and others providing new data on the role of hematology traits, particularly platelet function. Many of the replicated loci are outside of known or currently hypothesized pathways to thrombosis. CONCLUSIONS: Our cross-ancestry GWAS meta-analyses identified new loci associated with VTE. These findings highlight new pathways to thrombosis and provide novel molecules that may be useful in the development of improved antithrombosis treatments.

Allelic Heterogeneity at the CRP Locus Identified by Whole-Genome Sequencing in Multi-ancestry Cohorts.

Whole-genome sequencing (WGS) can improve assessment of low-frequency and rare variants, particularly in non-European populations that have been underrepresented in existing genomic studies. The genetic determinants of C-reactive protein (CRP), a biomarker of chronic inflammation, have been extensively studied, with existing genome-wide association studies (GWASs) conducted in >200,000 individuals of European ancestry. In order to discover novel loci associated with CRP levels, we examined a multi-ancestry population (n = 23,279) with WGS (∼38× coverage) from the Trans-Omics for Precision Medicine (TOPMed) program. We found evidence for eight distinct associations at the CRP locus, including two variants that have not been identified previously (rs11265259 and rs181704186), both of which are non-coding and more common in individuals of African ancestry (∼10% and ∼1% minor allele frequency, respectively, and rare or monomorphic in 1000 Genomes populations of East Asian, South Asian, and European ancestry). We show that the minor (G) allele of rs181704186 is associated with lower CRP levels and decreased transcriptional activity and protein binding in vitro, providing a plausible molecular mechanism for this African ancestry-specific signal. The individuals homozygous for rs181704186-G have a mean CRP level of 0.23 mg/L, in contrast to individuals heterozygous for rs181704186 with mean CRP of 2.97 mg/L and major allele homozygotes with mean CRP of 4.11 mg/L. This study demonstrates the utility of WGS in multi-ethnic populations to drive discovery of complex trait associations of large effect and to identify functional alleles in noncoding regulatory regions.

Predicted leukocyte telomere length and risk of germ cell tumours.

BACKGROUND: Genetically predicted leukocyte telomere length (LTL) has been evaluated in several studies of childhood and adult cancer. We test whether genetically predicted longer LTL is associated with germ cell tumours (GCT) in children and adults. METHODS: Paediatric GCT samples were obtained from a Children's Oncology Group study and state biobank programs in California and Michigan (N = 1413 cases, 1220 biological parents and 1022 unrelated controls). Replication analysis included 396 adult testicular GCTs (TGCT) and 1589 matched controls from the UK Biobank. Mendelian randomisation was used to look at the association between genetically predicted LTL and GCTs and TERT variants were evaluated within GCT subgroups. RESULTS: We identified significant associations between TERT variants reported in previous adult TGCT GWAS in paediatric GCT: TERT/rs2736100-C (OR = 0.82; P = 0.0003), TERT/rs2853677-G (OR = 0.80; P = 0.001), and TERT/rs7705526-A (OR = 0.81; P = 0.003). We also extended these findings to females and tumours outside the testes. In contrast, we did not observe strong evidence for an association between genetically predicted LTL by other variants and GCT risk in children or adults. CONCLUSION: While TERT is a known susceptibility locus for GCT, our results suggest that LTL predicted by other variants is not strongly associated with risk in either children or adults.

Genetic drivers of Cushing's disease: Frequency and associated phenotypes.

PURPOSE: Cushing's disease (CD) is often explained by a single somatic sequence change. Germline defects, however, often go unrecognized. We aimed to determine the frequency and associated phenotypes of genetic drivers of CD in a large cohort. METHODS: We studied 245 unrelated patients with CD (139 female, 56.7%), including 230 (93.9%) pediatric and 15 (6.1%) adult patients. Germline exome sequencing was performed in 184 patients; tumor exome sequencing was also done in 27 of them. A total of 43 germline samples and 92 tumor samples underwent Sanger sequencing of specific genes. Rare variants of uncertain significance, likely pathogenic (LP), or pathogenic variants in CD-associated genes, were identified. RESULTS: Germline variants (13 variants of uncertain significance, 8 LP, and 11 pathogenic) were found in 8 of 19 patients (42.1%) with positive family history and in 23 of 226 sporadic patients (10.2%). Somatic variants (1 LP and 7 pathogenic) were found in 20 of 119 tested individuals (16.8%); one of them had a coexistent germline defect. Altogether, variants of interest were identified at the germline level in 12.2% of patients, at the somatic level in 7.8%, and coexisting germline and somatic variants in 0.4%, accounting for one-fifth of the cohort. CONCLUSION: We report an estimate of the contribution of multiple germline and somatic genetic defects underlying CD in a single cohort.

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.

Progressive reduction in circulating levels of carotenoids and other micronutrients in patients with chronic pancreatitis.

BACKGROUND: Although micronutrients modulate immunity and inflammation, it remains elusive whether they are implicated in the development and progression of chronic pancreatitis (CP). This study aimed to investigate differences in the circulating levels of selected carotenoids and vitamins between CP and controls and trends in the levels of these micronutrients across controls, early CP, and definite CP. METHODS: Demographic and lifestyle data were extracted from medical records for 53 patients with CP (13 early and 38 definite) and obtained using a questionnaire for 52 controls. Plasma ß-carotene, lycopene, cryptoxanthin, zeaxanthin, and α-tocopherol and serum 25(OH)D, folate, IL-6, TNF-α, and MCP-1 were measured with state-of-the-art methods. RESULTS: The levels of all micronutrients (except folate) were significantly lower in CP than in controls. There was a progressive decrease in the levels of these micronutrients across controls, early CP, and definite CP (all p values for trend: ≤0.0012); e.g., plasma lycopene was 36.6, 21.5, and 14.5 µg/dL for controls, early CP, and definite CP, respectively. After adjustment for confounders, there were strong, inverse associations between the levels of all micronutrients (except folate) and CP (e.g., OR (95% CI) for ≥ median vs.

Whole genome sequence association with E-selectin levels reveals loss-of-function variant in African Americans.

E-selectin mediates the rolling of circulating leukocytes during inflammatory processes. Previous genome-wide association studies in European and Asian individuals have identified the ABO locus associated with E-selectin levels. Using Trans-Omics for Precision Medicine whole genome sequencing data in 2249 African Americans (AAs) from the Jackson Heart Study, we examined genome-wide associations with soluble E-selectin levels. In addition to replicating known signals at ABO, we identified a novel association of a common loss-of-function, missense variant in Fucosyltransferase 6 (FUT6; rs17855739,p.Glu274Lys, P = 9.02 × 10-24) with higher soluble E-selectin levels. This variant is considerably more common in populations of African ancestry compared to non-African ancestry populations. We replicated the association of FUT6 p.Glu274Lys with higher soluble E-selectin in an independent population of 748 AAs from the Women's Health Initiative and identified an additional pleiotropic association with vitamin B12 levels. Despite the broad role of both selectins and fucosyltransferases in various inflammatory, immune and cancer-related processes, we were unable to identify any additional disease associations of the FUT6 p.Glu274Lys variant in an electronic medical record-based phenome-wide association scan of over 9000 AAs.

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.

Exome sequencing of child-parent trios with bladder exstrophy: Findings in 26 children.

Bladder exstrophy (BE) is a rare, lower ventral midline defect with the bladder and part of the urethra exposed. The etiology of BE is unknown but thought to be influenced by genetic variation with more recent studies suggesting a role for rare variants. As such, we conducted paired-end exome sequencing in 26 child/mother/father trios. Three children had rare (allele frequency ≤ 0.0001 in several public databases) inherited variants in TSPAN4, one with a loss-of-function variant and two with missense variants. Two children had loss-of-function variants in TUBE1. Four children had rare missense or nonsense variants (one per child) in WNT3, CRKL, MYH9, or LZTR1, genes previously associated with BE. We detected 17 de novo missense variants in 13 children and three de novo loss-of-function variants (AKR1C2, PRRX1, PPM1D) in three children (one per child). We also detected rare compound heterozygous loss-of-function variants in PLCH2 and CLEC4M and rare inherited missense or loss-of-function variants in additional genes applying autosomal recessive (three genes) and X-linked recessive inheritance models (13 genes). Variants in two genes identified may implicate disruption in cell migration (TUBE1) and adhesion (TSPAN4) processes, mechanisms proposed for BE, and provide additional evidence for rare variants in the development of this defect.

Common α-globin variants modify hematologic and other clinical phenotypes in sickle cell trait and disease.

Co-inheritance of α-thalassemia has a significant protective effect on the severity of complications of sickle cell disease (SCD), including stroke. However, little information exists on the association and interactions for the common African ancestral α-thalassemia mutation (-α3.7 deletion) and ß-globin traits (HbS trait [SCT] and HbC trait) on important clinical phenotypes such as red blood cell parameters, anemia, and chronic kidney disease (CKD). In a community-based cohort of 2,916 African Americans from the Jackson Heart Study, we confirmed the expected associations between SCT, HbC trait, and the -α3.7 deletion with lower mean corpuscular volume/mean corpuscular hemoglobin and higher red blood cell count and red cell distribution width. In addition to the recently recognized association of SCT with lower estimated glomerular filtration rate and glycated hemoglobin (HbA1c), we observed a novel association of the -α3.7 deletion with higher HbA1c levels. Co-inheritance of each additional copy of the -α3.7 deletion significantly lowered the risk of anemia and chronic kidney disease among individuals with SCT (P-interaction = 0.031 and 0.019, respectively). Furthermore, co-inheritance of a novel α-globin regulatory variant was associated with normalization of red cell parameters in individuals with the -α3.7 deletion and significantly negated the protective effect of α-thalassemia on stroke in 1,139 patients with sickle cell anemia from the Cooperative Study of Sickle Cell Disease (CSSCD) (P-interaction = 0.0049). Functional assays determined that rs11865131, located in the major alpha-globin enhancer MCS-R2, was the most likely causal variant. These findings suggest that common α- and ß-globin variants interact to influence hematologic and clinical phenotypes in African Americans, with potential implications for risk-stratification and counseling of individuals with SCD and SCT.

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.