Jackson, Rebecca D.

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.

Introduction: Estrogens have been associated with risk of stroke in both men and women. Estradiol and estrone are synthesized from testosterone via the aromatase enzyme in both men and postmenopausal women. Genetic variation affecting function of the aromatase gene (CYP19A1) has been linked to elevated estradiol levels. We tested the hypothesis that these variants, as well as a modifying gene, are associated with ischemic stroke in three cohorts: the Vitamin Intervention for Stroke Prevention trial (VISP), Women’s Health Initiative (WHI), and Genetic Susceptibility to Early-Onset Stroke study (GEOS). Methods: Genomic data for individuals of European descent with ischemic stroke was compared against an equal number of controls: VISP (1723 cases - 65% men), WHI (436 cases - all postmenopausal women), GEOS (557 cases - all men age <55). Genotyping was funded through the Genomics and Randomized Trials Network (GARNET) and Gene Environment Association Studies (GENEVA) consortia. All cohorts used the Illumina Omni1-Quad chip for genotyping with imputation as needed. Eight candidate single nucleotide polymorphisms (SNPs) in CYP19A1 on chromosome 15 and an additional regulatory SNP in a modifier locus on chromosome 8 were chosen for analysis due to coding status or previous association with estradiol levels. Analyses were adjusted for age, hormone therapy use (women), and the top three principal components. Association results were determined as odds ratios for individual SNPs (effect allele = minor), with significance level of p<0.05. Results: GWAS results are tabulated. There were no significant differences in allele frequency for four SNPs in CYP19A1 among stroke patients and controls. The additional four SNPs were monomorphic in all three populations. The modifier gene SNP (rs1864729) was less common among young men with prior stroke in GEOS, but not men in VISP or in women. Conclusions: Aromatase genetic variation was not associated with risk of stroke in these three cohorts. One SNP that modifies aromatase function was less common among young, but not older men with stroke. Replication in other studies is necessary.

Background The association of social isolation or lack of social network ties in older adults is unknown. This knowledge gap is important since the risk of heart failure (HF) and social isolation increase with age. The study examines whether social isolation is associated with incident HF in older women, and examines depressive symptoms as a potential mediator and age and race and ethnicity as effect modifiers. Methods and Results This study included 44 174 postmenopausal women of diverse race and ethnicity from the WHI (Women's Health Initiative) study who underwent annual assessment for HF adjudication from baseline enrollment (1993–1998) through 2018. We conducted a mediation analysis to examine depressive symptoms as a potential mediator and further examined effect modification by age and race and ethnicity. Incident HF requiring hospitalization was the main outcome. Social isolation was a composite variable based on marital/partner status, religious ties, and community ties. Depressive symptoms were assessed using CES‐D (Center for Epidemiology Studies‐Depression). Over a median follow‐up of 15.0 years, we analyzed data from 36 457 women, and 2364 (6.5%) incident HF cases occurred; 2510 (6.9%) participants were socially isolated. In multivariable analyses adjusted for sociodemographic, behavioral, clinical, and general health/functioning; socially isolated women had a higher risk of incident HF than nonisolated women (HR, 1.23; 95% CI, 1.08–1.41). Adding depressive symptoms in the model did not change this association (HR, 1.22; 95% CI, 1.07–1.40). Neither race and ethnicity nor age moderated the association between social isolation and incident HF. Conclusions Socially isolated older women are at increased risk for developing HF, independent of traditional HF risk factors. Registration URL: http://www.clinicaltrials.gov ; Unique identifier: NCT00000611.

Abstract We previously identified five single nucleotide polymorphisms (SNPs) at four susceptibility loci for diffuse large B-cell lymphoma (DLBCL) in individuals of European ancestry through a large genome-wide association study (GWAS). To further elucidate genetic susceptibility to DLBCL, we sought to validate two loci at 3q13.33 and 3p24.1 that were suggestive in the original GWAS with additional genotyping. In the meta-analysis (5662 cases and 9237 controls) of the four original GWAS discovery scans and three replication studies, the 3q13.33 locus (rs9831894; minor allele frequency [MAF] = 0.40) was associated with DLBCL risk [odds ratio (OR) = 0.83, P = 3.62 × 10−13]. rs9831894 is in linkage disequilibrium (LD) with additional variants that are part of a super-enhancer that physically interacts with promoters of CD86 and ILDR1. In the meta-analysis (5510 cases and 12 817 controls) of the four GWAS discovery scans and four replication studies, the 3p24.1 locus (rs6773363; MAF = 0.45) was also associated with DLBCL risk (OR = 1.20, P = 2.31 × 10−12). This SNP is 29 426-bp upstream of the nearest gene EOMES and in LD with additional SNPs that are part of a highly lineage-specific and tumor-acquired super-enhancer that shows long-range interaction with AZI2 promoter. These loci provide additional evidence for the role of immune function in the etiology of DLBCL, the most common lymphoma subtype.

Abstract Background: We previously identified 5 SNPs at 4 susceptibility loci for diffuse large B-cell lymphoma (DLBCL) in individuals of European ancestry through a large genome-wide association study (Cerhan et al, Nat Gen 2014;46:1233-8); however, much of the heritability remains unexplained. To further elucidate genetic susceptibility to DLBCL, we sought to validate 2 loci at 3q13.33 and 3p24.1 that were suggestive in the original report. Methods: We selected two SNPs rs9831894 (3q13.33) and rs6773363 (3p24.1) for de novo replication from studies at the Mayo Clinic, MD Anderson, and Memorial Sloan Kettering. Logistic regression was used to estimate odds ratios (ORs), using a log-additive model, and adjusting for age, gender and significant eigenvectors for each genotyping center separately. For both SNPs, we then conducted a meta-analysis of all replication studies (n=3) and our original GWAS studies (n=4), encompassing in total 5662 cases and 9237 controls for rs9831894, and 5510 cases and 12,817 controls for rs6773363. Meta-analysis was conducted using the fixed-effects inverse variance method based on the β estimates and standard errors from each study. Results: In a meta-analysis of the 4 original GWAS scans (3856 cases and 7666 controls), rs9831894 (MAF=0.40) was associated with DLBCL risk (OR=0.84, P=4.52 x 10-9) and this replicated in a meta-analysis of the 3 studies with de novo genotyping (OR=0.80, P=4.17 x 10-5); the overall meta-analysis showed a strong association with DLBCL risk (OR=0.83, P=3.62 x 10-13). This locus maps near a plausible candidate gene, CD86, a protein coding gene and a member of the immunoglobulin superfamily that encodes a type I membrane protein. Binding of CD86 with cytotoxic T-lymphocyte-associated protein 4 negatively regulates T-cell activation and diminishes the immune response, while binding of CD86 with CD28 antigen is a costimulatory signal for activation of the T-cell. In a meta-analysis of the same 4 GWAS scans, rs6773363 (MAF=0.45) was tentatively associated with DLBCL risk (OR=1.17, P= 3.68 x 10-7) and this replicated in a meta-analysis of the studies with de novo genotyping (OR=1.27, P=3.78 x 10-7); the overall meta-analysis showed a strong association with DLBCL risk (OR=1.20, P=2.31 x 10-12). This locus also maps near a plausible candidate gene, eomesodermin (EOMES), a transcription factor crucial for embryonic development of the central nervous system and also putatively involved in T-cell differentiation in viral infection defense. Conclusion: In this follow-up analysis of our initial GWAS, we have identified two additional loci associated with risk of DLBCL, the most common lymphoma subtype. These loci provide additional evidence for the role of immune function in the etiology of DLBCL. Citation Format: Geffen Kleinstern, Michelle Hildebrand, Vijai Joseph, Sonja I. Berndt, Hervé Ghesquières, James McKay, Sophia S. Wang, Alexandra Nieters, David Cox, Alain Monnereau, Angela R. Brooks-Wilson, Qing Lan, Mads Melbye, Rebecca D. Jackson, Lauren R. Teras, Mark P. Purdue, Claire M. Vajdic, Demetrius Albanes, Anne Zeleniuch-Jacquotte, Simon Crouch, Yawei Zhang, Susan L. Slager, Xifeng Wu, Karin E. Smedby, Gilles Salles, Christine F. Skibola, Nathaniel Rothman, Stephen J. Chanock, James R. Cerhan. Inherited variants at 3q13.33 and 3p24.1 influences risk of diffuse large B-cell lymphoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 221.

Abstract Background Previous genome-wide association studies (GWAS) have identified 42 loci (P < 5 × 10−8) associated with risk of colorectal cancer (CRC). Expanded consortium efforts facilitating the discovery of additional susceptibility loci may capture unexplained familial risk. Methods We conducted a GWAS in European descent CRC cases and control subjects using a discovery–replication design, followed by examination of novel findings in a multiethnic sample (cumulative n = 163 315). In the discovery stage (36 948 case subjects/30 864 control subjects), we identified genetic variants with a minor allele frequency of 1% or greater associated with risk of CRC using logistic regression followed by a fixed-effects inverse variance weighted meta-analysis. All novel independent variants reaching genome-wide statistical significance (two-sided P < 5 × 10−8) were tested for replication in separate European ancestry samples (12 952 case subjects/48 383 control subjects). Next, we examined the generalizability of discovered variants in East Asians, African Americans, and Hispanics (12 085 case subjects/22 083 control subjects). Finally, we examined the contributions of novel risk variants to familial relative risk and examined the prediction capabilities of a polygenic risk score. All statistical tests were two-sided. Results The discovery GWAS identified 11 variants associated with CRC at P < 5 × 10−8, of which nine (at 4q22.2/5p15.33/5p13.1/6p21.31/6p12.1/10q11.23/12q24.21/16q24.1/20q13.13) independently replicated at a P value of less than .05. Multiethnic follow-up supported the generalizability of discovery findings. These results demonstrated a 14.7% increase in familial relative risk explained by common risk alleles from 10.3% (95% confidence interval [CI] = 7.9% to 13.7%; known variants) to 11.9% (95% CI = 9.2% to 15.5%; known and novel variants). A polygenic risk score identified 4.3% of the population at an odds ratio for developing CRC of at least 2.0. Conclusions This study provides insight into the architecture of common genetic variation contributing to CRC etiology and improves risk prediction for individualized screening.

Objective Determinants of the increased risk of diffuse large B-cell lymphoma (DLBCL) in SLE are unclear. Using data from a recent lymphoma genome-wide association study (GWAS), we assessed whether certain lupus-related single nucleotide polymorphisms (SNPs) were also associated with DLBCL. Methods GWAS data on European Caucasians from the International Lymphoma Epidemiology Consortium (InterLymph) provided a total of 3857 DLBCL cases and 7666 general-population controls. Data were pooled in a random-effects meta-analysis. Results Among the 28 SLE-related SNPs investigated, the two most convincingly associated with risk of DLBCL included the CD40 SLE risk allele rs4810485 on chromosome 20q13 (OR per risk allele=1.09, 95% CI 1.02 to 1.16, p=0.0134), and the HLA SLE risk allele rs1270942 on chromosome 6p21.33 (OR per risk allele=1.17, 95% CI 1.01 to 1.36, p=0.0362). Of additional possible interest were rs2205960 and rs12537284. The rs2205960 SNP, related to a cytokine of the tumour necrosis factor superfamily TNFSF4, was associated with an OR per risk allele of 1.07, 95% CI 1.00 to 1.16, p=0.0549. The OR for the rs12537284 (chromosome 7q32, IRF5 gene) risk allele was 1.08, 95% CI 0.99 to 1.18, p=0.0765. Conclusions These data suggest several plausible genetic links between DLBCL and SLE.

Genome-wide association studies have identified >50 common variants associated with kidney function, but these variants do not fully explain the variation in eGFR. We performed a two-stage meta-analysis of associations between genotypes from the Illumina exome array and eGFR on the basis of serum creatinine (eGFRcrea) among participants of European ancestry from the CKDGen Consortium ( n Stage1 : 111,666; n Stage2 : 48,343). In single-variant analyses, we identified single nucleotide polymorphisms at seven new loci associated with eGFRcrea ( PPM1J , EDEM3, ACP1, SPEG, EYA4, CYP1A1 , and ATXN2L ; P Stage1 <3.7×10 −7 ), of which most were common and annotated as nonsynonymous variants. Gene-based analysis identified associations of functional rare variants in three genes with eGFRcrea, including a novel association with the SOS Ras/Rho guanine nucleotide exchange factor 2 gene, SOS2 ( P =5.4×10 −8 by sequence kernel association test). Experimental follow-up in zebrafish embryos revealed changes in glomerular gene expression and renal tubule morphology in the embryonic kidney of acp1- and sos2 -knockdowns. These developmental abnormalities associated with altered blood clearance rate and heightened prevalence of edema. This study expands the number of loci associated with kidney function and identifies novel genes with potential roles in kidney formation.

Abstract Background: Diffuse large B-cell lymphoma (DLBCL) is the most common lymphoma subtype and is clinically aggressive. Relatives of DLBCL patients are at elevated risk for DLBCL. In small genome-wide association studies (GWAS) in which all NHL subtypes were combined, no conclusive loci for DLBCL have been previously identified in individuals of European background, whereas a recent study conducted in East Asia identified a locus at 3q27. Methods: To identify genetic susceptibility loci for DLBCL, we conducted a meta-analysis of three new GWAS and one prior scan, totaling 3,857 DLBCL cases and 7,666 controls of European ancestry. Imputation was used to combine all the data using the 1000 Genomes Project release v3 and IMPUTE2. Based on these findings, we selected 6 promising SNPs for genotyping in an additional 202 cases and 3,431 controls, which were then included in a second meta-analysis. SNPTEST was used to estimate odds ratios (ORs), using the additive model and adjusting for age, gender and significant eigenvectors. Results: In our multi-stage analysis, five independent SNPs in four loci achieved genome-wide significance marked by rs116446171 at 6p25.3 (EXOC2; OR=2.33, P=1.60x10-20), rs13255292 and rs4733601 at 8q24.21 (PVT1; OR=1.23, P=4.86x10-12 and OR=1.18, P=1.63x10-8, respectively), rs2523607 at 6p21.33 (HLA-B; OR=1.34, 3.35x10-9) and rs2681416 at 3q13.33 (CD86; OR=1.20, P=6.04x10-9). The two 8q24.21 SNPs displayed minimal linkage disequilibrium (r2=0.03) and conditional analysis supported their independence. We did not observe a notable signal for a locus previously reported for DLBCL on 3q27 in East Asia, rs6773854 (OR=1.06, P=0.81). We estimated that common SNPs, including but not limited to the loci discovered in this study, explain approximately 16% of the variance for DLBCL risk overall. Discussion: The susceptibility locus at 6p25.3 maps near a plausible DLBCL candidate gene, EXOC2 (exocyst complex component 2), which interacts with Ral proteins, and disruption of the Ral-exocyst regulatory node impacts cellular and developmental processes associated with malignant transformation and progression. The two 8q24.21 variants are approximately 1Mb telomeric to the 8q24 region linked with multiple cancers including CLL. Both variants are positioned in close proximity to PVT1, which is a non-coding RNA implicated in the MYC activation, and MYC is known to be deregulated in some DLBCLs. The locus at 3q13.33 maps to CD86 (B7-2), a B-7 family member that along with CD80 (B7-1) is essential in cognate interactions between T- and B-lymphocytes; both play key roles in immunodeficiency, autoimmune diseases, and in anti-tumor immunity. Finally, we had multiple signals in the HLA region, the strongest of which was at 6p21.33 near HLA-B. HLA-B encodes the HLA class I heavy chain paralogue, which heterodimerizes with a light chain (β2 microglobulin) and plays a central role in presenting intracellularly processed self or foreign antigens to CD8+ cytotoxic T lymphocytes. Class I molecules have been extensively linked to the risk of a variety of immune-mediated diseases and cancers including DLBCL. Conclusions: These data provide substantial new evidence for genetic susceptibility to this B-cell malignancy, and point towards pathways involved in immune recognition and immune function in the pathogenesis of DLBCL. Citation Format: James R. Cerhan, Sonja I. Berndt, Joseph Vijai, Hervé Ghesquières, James McKay, Sophia S. Wang, Zhaoming Wang, Meredith Yeager, Alexandra Nieters, David Cox, Alain Monnereau, Christopher Flowers, Anneclaire J. De Roos, Angela R. Brooks-Wilson, Qing Lan, Gianluca Severi, Mads Melbye, Rebecca D. Jackson, Lauren R. Teras, Mark Purdue, Claire Vajdic, Demetrius Albanes, Kimberly A. Bertrand, Anne Zeleniuch-Jacquotte, Simon Crouch, Yawei Zhang, Paolo Vineis, Susan L. Slager, Karin E. Smedby, Gilles Salles, Christine F. Skibola, Nathaniel Rothman, Stephen J. Chanock, on behalf of the NHL GWAS Project. Genome-wide association study identifies multiple susceptibility loci for diffuse large B-cell lymphoma. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr LB-272. doi:10.1158/1538-7445.AM2014-LB-272