Teras, Lauren R.

PURPOSE The prevalence of germline pathogenic variants (PVs) in established breast cancer predisposition genes in women in the general population over age 65 years is not well-defined. However, testing guidelines suggest that women diagnosed with breast cancer over age 65 years might have < 2.5% likelihood of a PV in a high-penetrance gene. This study aimed to establish the frequency of PVs and remaining risks of breast cancer for each gene in women over age 65 years. METHODS A total of 26,707 women over age 65 years from population-based studies (51.5% with breast cancer and 48.5% unaffected) were tested for PVs in germline predisposition gene. Frequencies of PVs and associations between PVs in each gene and breast cancer were assessed, and remaining lifetime breast cancer risks were estimated for non-Hispanic White women with PVs. RESULTS The frequency of PVs in predisposition genes was 3.18% for women with breast cancer and 1.48% for unaffected women over age 65 years. PVs in BRCA1, BRCA2, and PALB2 were found in 3.42% of women diagnosed with estrogen receptor (ER)–negative, 1.0% with ER-positive, and 3.01% with triple-negative breast cancer. Frequencies of PVs were lower among women with no first-degree relatives with breast cancer. PVs in CHEK2, PALB2, BRCA2, and BRCA1 were associated with increased risks (odds ratio = 2.9-4.0) of breast cancer. Remaining lifetime risks of breast cancer were ≥ 15% for those with PVs in BRCA1, BRCA2, and PALB2. CONCLUSION This study suggests that all women diagnosed with triple-negative breast cancer or ER-negative breast cancer should receive genetic testing and that women over age 65 years with BRCA1 and BRCA2 PVs and perhaps with PALB2 and CHEK2 PVs should be considered for magnetic resonance imaging screening.

Abstract Background: Use of menopausal hormone therapy (MHT) is associated with increased risk for breast cancer. However, the relevant mechanisms and its interaction with genetic variants are not fully understood. Methods: We conducted a genome-wide interaction analysis between MHT use and genetic variants for breast cancer risk in 27,585 cases and 34,785 controls from 26 observational studies. All women were post-menopausal and of European ancestry. Multivariable logistic regression models were used to test for multiplicative interactions between genetic variants and current MHT use. We considered interaction p-values<5x10-8 as genome-wide significant, and p-values<1x10-5 as suggestive. Linkage disequilibrium (LD)-based clumping was performed to identify independent candidate variants.Results: None of the 9.7 million genetic variants tested for interactions with MHT use reached genome-wide significance. Only 213 variants, representing 18 independent loci, had p-values<1x105. The strongest evidence was found for rs4674019 (p-value=2.27x10-7), which showed genome-wide significant interaction (p-value=3.8x10-8) with current MHT use when analysis was restricted to population-based studies only. Limiting the analyses to combined estrogen-progesterone MHT use only or to estrogen receptor (ER) positive cases did not identify any genome-wide significant evidence of interactions. Conclusions: In this large genome-wide SNP-MHT interaction study of breast cancer, we found no strong support for common genetic variants modifying the effect of MHT on breast cancer risk. These results suggest that common genetic variation has limited impact on the observed MHT–breast cancer risk association.

Abstract BackgroundCopy number variants (CNVs) are pervasive in the human genome but potential disease associations with rare CNVs have not been comprehensively assessed in large datasets. We analysed rare CNVs in genes and non-coding regions for 86,788 breast cancer cases and 76,122 controls of European ancestry with genome-wide array data.ResultsGene burden tests detected the strongest association for deletions in BRCA1 (P= 3.7E-18). Nine other genes were associated with a p-value < 0.01 including known susceptibility genes CHEK2 (P= 0.0008), ATM (P= 0.002) and BRCA2 (P= 0.008). Outside the known genes we detected associations with p-values < 0.001 for either overall or subtype-specific breast cancer at nine deletion regions and four duplication regions. Three of the deletion regions were in established common susceptibility loci.ConclusionsThis is the first genome-wide analysis of rare CNVs in a large breast cancer case-control dataset. We detected associations with exonic deletions in established breast cancer susceptibility genes. We also detected suggestive associations with non-coding CNVs in known and novel loci with large effects sizes. Larger sample sizes will be required to reach robust levels of statistical significance.

PURPOSE This study assessed the joint association of pathogenic variants (PVs) in breast cancer (BC) predisposition genes and polygenic risk scores (PRS) with BC in the general population. METHODS A total of 26,798 non-Hispanic white BC cases and 26,127 controls from predominately population-based studies in the Cancer Risk Estimates Related to Susceptibility consortium were evaluated for PVs in BRCA1, BRCA2, ATM, CHEK2, PALB2, BARD1, BRIP1, CDH1, and NF1. PRS based on 105 common variants were created using effect estimates from BC genome-wide association studies; the performance of an overall BC PRS and estrogen receptor–specific PRS were evaluated. The odds of BC based on the PVs and PRS were estimated using penalized logistic regression. The results were combined with age-specific incidence rates to estimate 5-year and lifetime absolute risks of BC across percentiles of PRS by PV status and first-degree family history of BC. RESULTS The estimated lifetime risks of BC among general-population noncarriers, based on 10th and 90th percentiles of PRS, were 9.1%-23.9% and 6.7%-18.2% for women with or without first-degree relatives with BC, respectively. Taking PRS into account, more than 95% of BRCA1, BRCA2, and PALB2 carriers had > 20% lifetime risks of BC, whereas, respectively, 52.5% and 69.7% of ATM and CHEK2 carriers without first-degree relatives with BC, and 78.8% and 89.9% of those with a first-degree relative with BC had > 20% risk. CONCLUSION PRS facilitates personalization of BC risk among carriers of PVs in predisposition genes. Incorporating PRS into BC risk estimation may help identify > 30% of CHEK2 and nearly half of ATM carriers below the 20% lifetime risk threshold, suggesting the addition of PRS may prevent overscreening and enable more personalized risk management approaches.

Abstract Background: Use of menopausal hormone therapy (MHT) is associated with increased risk for breast cancer. However, the relevant mechanisms and its interaction with genetic variants are not fully understood. Methods: We conducted a genome-wide interaction analysis between MHT use and genetic variants for breast cancer risk in 27,585 cases and 34,785 controls from 26 observational studies. All women were post-menopausal and of European ancestry. Multivariable logistic regression models were used to test for multiplicative interactions between genetic variants and current MHT use. We considered interaction p-values<5x10-8 as genome-wide significant, and p-values<1x10-5 as suggestive. Linkage disequilibrium (LD)-based clumping was performed to identify independent candidate variants.Results: None of the 9.7 million genetic variants tested for interactions with MHT use reached genome-wide significance. Only 213 variants, representing 18 independent loci, had p-values<1x105. The strongest evidence was found for rs4674019 (p-value=2.27x10-7), which showed genome-wide significant interaction (p-value=3.8x10-8) with current MHT use when analysis was restricted to population-based studies only. Limiting the analyses to combined estrogen-progesterone MHT use only or to estrogen receptor (ER) positive cases did not identify any genome-wide significant evidence of interactions. Conclusions: In this large genome-wide SNP-MHT interaction study of breast cancer, we found no strong support for common genetic variants modifying the effect of MHT on breast cancer risk. These results suggest that common genetic variation has limited impact on the observed MHT–breast cancer risk association.

In this study we aim to examine gene–environment interactions (GxEs) between genes involved with estrogen metabolism and environmental factors related to estrogen exposure. GxE analyses were conducted with 1970 Korean breast cancer cases and 2052 controls in the case-control study, the Seoul Breast Cancer Study (SEBCS). A total of 11,555 SNPs from the 137 candidate genes were included in the GxE analyses with eight established environmental factors. A replication test was conducted by using an independent population from the Breast Cancer Association Consortium (BCAC), with 62,485 Europeans and 9047 Asians. The GxE tests were performed by using two-step methods in GxEScan software. Two interactions were found in the SEBCS. The first interaction was shown between rs13035764 of NCOA1 and age at menarche in the GE|2df model (p-2df = 1.2 × 10−3). The age at menarche before 14 years old was associated with the high risk of breast cancer, and the risk was higher when subjects had homozygous minor allele G. The second GxE was shown between rs851998 near ESR1 and height in the GE|2df model (p-2df = 1.1 × 10−4). Height taller than 160 cm was associated with a high risk of breast cancer, and the risk increased when the minor allele was added. The findings were not replicated in the BCAC. These results would suggest specificity in Koreans for breast cancer risk.

Abstract Background: Lipid traits have been inconsistently linked to risk of non-Hodgkin lymphoma (NHL). We examined the association of genetically predicted lipid traits with risk of diffuse large B-cell lymphoma (DLBCL), chronic lymphocytic leukemia (CLL), follicular lymphoma (FL), and marginal zone lymphoma (MZL) using Mendelian randomization (MR) analysis. Methods: Genome-wide association study data from the InterLymph Consortium were available for 2,661 DLBCLs, 2,179 CLLs, 2,142 FLs, 824 MZLs, and 6,221 controls. SNPs associated (P &lt; 5 × 10−8) with high-density lipoprotein (HDL, n = 164), low-density lipoprotein (LDL, n = 137), total cholesterol (TC, n = 161), and triglycerides (TG, n = 123) were used as instrumental variables (IV), explaining 14.6%, 27.7%, 16.8%, and 12.8% of phenotypic variation, respectively. Associations between each lipid trait and NHL subtype were calculated using the MR inverse variance–weighted method, estimating odds ratios (OR) per standard deviation and 95% confidence intervals (CI). Results: HDL was positively associated with DLBCL (OR = 1.14; 95% CI, 1.00–1.30) and MZL (OR = 1.09; 95% CI, 1.01–1.18), while TG was inversely associated with MZL risk (OR = 0.90; 95% CI, 0.83–0.99), all at nominal significance (P &lt; 0.05). A positive trend was observed for HDL with FL risk (OR = 1.08; 95% CI, 0.99–1.19; P = 0.087). No associations were noteworthy after adjusting for multiple testing. Conclusions: We did not find evidence of a clear or strong association of these lipid traits with the most common NHL subtypes. While these IVs have been previously linked to other cancers, our findings do not support any causal associations with these NHL subtypes. Impact: Our results suggest that prior reported inverse associations of lipid traits are not likely to be causal and could represent reverse causality or confounding.

Introduction. Acquired severe aplastic anemia (SAA) is a life-threatening disorder characterized by severe progressive pancytopenia and hypocellular bone marrow. The etiology of acquired SAA is not understood but believed to be related to abnormal immune responses to environmental exposures. We conducted a genome-wide association study (GWAS) to identify common germline variants associated with SAA. Methods. We identified 895 patients with SAA who underwent related or unrelated hematopoietic cell transplant (HCT) with clinical data and pre-HCT blood samples available in the Center for International Blood and Marrow Transplant Research (CIBMTR) database and biorepository. Pre-HCT DNA was extracted from blood of patients with SAA and genome-wide genotyping was conducted using the Illumina OmniExpress array. We excluded 93 inherited bone marrow failure cases. The SAA cases were grouped cases into discovery and validation sets based on time of batch sample receipt. Analyses were limited to patients of European ancestry based on principal component analyses to minimize the potential effect of population stratification. Controls were genomically matched and selected from previously scanned cancer-free subjects at the Cancer Genomics Research Laboratory, Division of Cancer Epidemiology and Genetics, NCI. The final analysis included 534 acquired SAA cases (359 in the discovery set and 175 in the validation set), and 2,455 controls (1,396 in the discovery set, and 1,059 in the validation set). Results. Patients with SAA in this study received HCT between 1989-2015 at a median age of 21 years, 56% were male, and the median time between SAA diagnosis and HCT was 11 months. Strong genome-wide association signals were identified across the human leukocyte antigen (HLA) genes encoded at the major histocompatibility complex (MHC) on chromosome 6p21 (Figure 1). The top SNP was located in the P4 binding pocket of the HLA class II gene HLA-DPB1(rs1042151A>G, p.Met105Val, pooled-odds ratio [OR]=1.75, 95% confidence interval [CI]=1.50-2.03, p=1.94x10-13). The expression of HLA-DP in CD19+ cells from 175 healthy donors was significantly different by rs1042151 A>G genotype (p=2.04x10-6) (Figure 2). A second SNP near HLA-B, rs28367832G>A, also reached genome-wide significance (pooled-OR=1.49, 95% CI=1.22-1.78, p=7.27x10-9). Copy-number variant analysis and next generation sequencing also identified somatic, clonal copy-neutral loss-of-heterozygosity affecting class I HLA genes in 8.6% of the SAA cases and none of the controls. Conclusion. This SAA GWAS identified strong association signals between common germline genetic variants in HLA class I and II genes and SAA. The main SNP is associated with changes in HLA-DP expression suggesting a key role for this locus in SAA etiology. This study adds further evidence to the connection between SAA and immune dysregulation. Disclosures Cerhan: Janssen: Membership on an entity's Board of Directors or advisory committees; NanoString: Research Funding; Celgene: Research Funding. Lee:Incyte: Research Funding; Syndax: Research Funding; Amgen: Research Funding; Novartis: Research Funding; Takeda: Research Funding; Kadmon: Consultancy, Research Funding; Pfizer: Consultancy, Research Funding; AstraZeneca: Research Funding.

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.

AbstractEpidemiologic studies show an increased risk of non‐Hodgkin lymphoma (NHL) in patients with autoimmune disease (AD), due to a combination of shared environmental factors and/or genetic factors, or a causative cascade: chronic inflammation/antigen‐stimulation in one disease leads to another. Here we assess shared genetic risk in genome‐wide‐association‐studies (GWAS).Secondary analysis of GWAS of NHL subtypes (chronic lymphocytic leukemia, diffuse large B‐cell lymphoma, follicular lymphoma, and marginal zone lymphoma) and ADs (rheumatoid arthritis, systemic lupus erythematosus, and multiple sclerosis). Shared genetic risk was assessed by (a) description of regional genetic of overlap, (b) polygenic risk score (PRS), (c)"diseasome", (d)meta‐analysis.Descriptive analysis revealed few shared genetic factors between each AD and each NHL subtype. The PRS of ADs were not increased in NHL patients (nor vice versa). In the diseasome, NHLs shared more genetic etiology with ADs than solid cancers (p = .0041). A meta‐analysis (combing AD with NHL) implicated genes of apoptosis and telomere length.This GWAS‐based analysis four NHL subtypes and three ADs revealed few weakly‐associated shared loci, explaining little total risk. This suggests common genetic variation, as assessed by GWAS in these sample sizes, may not be the primary explanation for the link between these ADs and NHLs.