Orr, Nick

Abstract Chromosome 5p15.33 harbors several independent association signals which demonstrate antagonistic pleiotropy across cancer types, with causal mechanisms largely unresolved. To identify functional variants and enhancer elements at this locus, we performed statistical fine-mapping followed by massively parallel reporter assays (MPRA) and proliferation based CRISPRi screens. This approach identified eight multi-cancer functional variants (MCFVs) across three GWAS signals. Targeting rs421629 (part of the CLPTM1L signal marked by rs465498) with CRISPRi revealed opposing effects on TERT expression in pancreatic versus lung cancer cells, consistent with the antagonistic pleiotropy observed for this signal. Furthermore, CRISPRi nominated an intronic CLPTM1L variable number tandem repeat (VNTR) as a potent enhancer. Long-read sequencing established VNTR polymorphisms as potential causal variants for the rs465498 signal. We showed that Hippo-pathway transcription factors mediate VNTR enhancer activity in lung and pancreatic cancer cells. Together, these findings indicate that cancer susceptibility at 5p15.33 may be mediated by both SNPs and VNTRs and provide an integrated framework for resolving complex pleiotropic loci.

Genome-wide association studies (GWAS) have identified more than 200 susceptibility loci for breast cancer, but these variants explain less than a fifth of the disease risk. Although gene–environment interactions have been proposed to account for some of the remaining heritability, few studies have empirically assessed this. We obtained genotype and risk factor data from 46,060 cases and 47,929 controls of European ancestry from population-based studies within the Breast Cancer Association Consortium (BCAC). We built gene expression prediction models for 4,864 genes with a significant (P < 0.01) heritable component using the transcriptome and genotype data from the Genotype-Tissue Expression (GTEx) project. We leveraged predicted gene expression information to investigate the interactions between gene-centric genetic variation and 14 established risk factors in association with breast cancer risk, using a mixed-effects score test. After adjusting for number of tests using Bonferroni correction, no interaction remained statistically significant. The strongest interaction observed was between the predicted expression of the C13orf45 gene and age at first full-term pregnancy (PGXE = 4.44 × 10−6). In this transcriptome-informed genome-wide gene–environment interaction study of breast cancer, we found no strong support for the role of gene expression in modifying the associations between established risk factors and breast cancer risk. Our study suggests a limited role of gene–environment interactions in breast cancer risk.

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 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 We evaluated the joint associations between a new 313-variant PRS (PRS313) and questionnaire-based breast cancer risk factors for women of European ancestry, using 72 284 cases and 80 354 controls from the Breast Cancer Association Consortium. Interactions were evaluated using standard logistic regression and a newly developed case-only method for breast cancer risk overall and by estrogen receptor status. After accounting for multiple testing, we did not find evidence that per-standard deviation PRS313 odds ratio differed across strata defined by individual risk factors. Goodness-of-fit tests did not reject the assumption of a multiplicative model between PRS313 and each risk factor. Variation in projected absolute lifetime risk of breast cancer associated with classical risk factors was greater for women with higher genetic risk (PRS313 and family history) and, on average, 17.5% higher in the highest vs lowest deciles of genetic risk. These findings have implications for risk prevention for women at increased risk of breast cancer.

Abstract Background: Prospective validation of breast cancer risk models integrating classical risk factors and genetic variants is required for risk-stratified prevention and screening strategies. The objective of this study was to validate a breast cancer risk model integrating classical risk factors and a 313-variant polygenic risk score (PRS) in multiple prospective cohort studies, and to project five-year risk of breast cancer in six different countries. Methods: The study population included 7,529 cases and 230,103 controls from 14 prospective cohort studies in Australia, Germany, the Netherlands, Sweden, UK, and USA. We used the Individualized Coherent Absolute Risk Estimator (iCARE) tool for risk model building, validation, and risk projection. Expected five-year risk of invasive or in situ breast cancer was compared to observed risk, overall and within deciles of expected risk using goodness of fit statistics. We evaluated calibration of the relative risk through meta-analysis across cohorts, and of the absolute risk within each cohort. Model discrimination was evaluated using the area under the curve (AUC), and percentages of women crossing risk thresholds. Projections of five-year risk distributions were estimated for women of European ancestry aged 50-70 years in the general populations of these six countries. Results: Analysis showed overall good calibration of the integrated iCARE-based model relative risk for both women younger than 50 years (χ2=14.9, P=0.09) and aged 50 years or older (χ2=14.6, P=0.10), with a small overestimation of risk for women in the highest decile of expected risk (RR = 3.5 expected vs 2.3 (95% CI 1.6 to 3.2) observed for women &lt;50 years; and 2.8 expected vs 2.3 (95% CI 2.0 to 2.7) observed for women 50+ years). The age-adjusted AUCs for the integrated model were 63.1 (95% CI 60.9 to 65.3) and 62.9 (95% CI 61.8 to 64.0), for the two age groups respectively. The calibration of absolute risk showed substantial variation across cohorts, particularly for the older group, but had no systematic bias. Model based projections in the general populations showed that compared to the population average, women in the 1st and 99th percentiles of the integrated risk score had relative risks 0.19 and 3.56 respectively. The proportion of women of European ancestry aged 50-70 years with a five-year risk greater than 3% (threshold for consideration of risk-lowering drugs by U.S. Preventive Services Task Force) ranged from 7.1% in Germany to 18.2% in the US, which corresponds to ~5.5 million women in the US. Conclusions: Five-year risk predictions from a model with classical risk factors and PRS are well calibrated and provide substantial risk stratification across multiple cohorts in six different countries. Further studies are needed to evaluate the clinical utility of the validated model for risk stratified screening and prevention of breast cancer. Citation Format: Parichoy Pal Choudhury, Amber Wilcox, Chi Gao, Brian Carter, Anika Husing, Mark Brook, Mikael Eriksson, Kara Martin, Chris Scott, Min Shi, Thomas Ahearn, Michael Jones, Nick Orr, Minouk Schoemaker, Kamila Czene, Jenny Chang-Claude, Jacques Simard, Doug Easton, Marjanka K. Schmidt, Dale Sandler, Clarice R. Weinberg, Celine Vachon, Roger Milne, Per Hall, Anthony Swerdlow, Rudolph Kaaks, Myrto Barrdahl, Mia Gaudet, Antonis Antoniou, Peter Kraft, Montserrat Garcia-Closas, Nilanjan Chatterjee. Validation of breast cancer risk model incorporating classical risk factors and polygenic risk scores in 14 prospective cohort studies in 6 countries [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 962.

AbstractFemale Hodgkin lymphoma (HL) patients treated with chest radiotherapy (RT) have a very high risk of breast cancer. The contribution of genetic factors to this risk is unclear. We therefore examined 211 155 germline single-nucleotide polymorphisms (SNPs) for gene-radiation interaction on breast cancer risk in a case-only analysis including 327 breast cancer patients after chest RT for HL and 4671 first primary breast cancer patients. Nine SNPs showed statistically significant interaction with RT on breast cancer risk (false discovery rate, &lt;20%), of which 1 SNP in the PVT1 oncogene attained the Bonferroni threshold for statistical significance. A polygenic risk score (PRS) composed of these SNPs (RT-interaction-PRS) and a previously published breast cancer PRS (BC-PRS) derived in the general population were evaluated in a case-control analysis comprising the 327 chest-irradiated HL patients with breast cancer and 491 chest-irradiated HL patients without breast cancer. Patients in the highest tertile of the RT-interaction-PRS had a 1.6-fold higher breast cancer risk than those in the lowest tertile. Remarkably, we observed a fourfold increased RT-induced breast cancer risk in the highest compared with the lowest decile of the BC-PRS. On a continuous scale, breast cancer risk increased 1.4-fold per standard deviation of the BC-PRS, similar to the effect size found in the general population. This study demonstrates that genetic factors influence breast cancer risk after chest RT for HL. Given the high absolute breast cancer risk in radiation-exposed women, these results can have important implications for the management of current HL survivors and future patients.

Background The rarity of mutations in PALB2, CHEK2 and ATM make it difficult to estimate precisely associated cancer risks. Population-based family studies have provided evidence that at least some of these mutations are associated with breast cancer risk as high as those associated with rare BRCA2 mutations. We aimed to estimate the relative risks associated with specific rare variants in PALB2, CHEK2 and ATM via a multicentre case-control study. Methods We genotyped 10 rare mutations using the custom iCOGS array: PALB2 c.1592delT, c.2816T>G and c.3113G>A, CHEK2 c.349A>G, c.538C>T, c.715G>A, c.1036C>T, c.1312G>T, and c.1343T>G and ATM c.7271T>G. We assessed associations with breast cancer risk (42 671 cases and 42 164 controls), as well as prostate (22 301 cases and 22 320 controls) and ovarian (14 542 cases and 23 491 controls) cancer risk, for each variant. Results For European women, strong evidence of association with breast cancer risk was observed for PALB2 c.1592delT OR 3.44 (95% CI 1.39 to 8.52, p=7.1×10−5), PALB2 c.3113G>A OR 4.21 (95% CI 1.84 to 9.60, p=6.9×10−8) and ATM c.7271T>G OR 11.0 (95% CI 1.42 to 85.7, p=0.0012). We also found evidence of association with breast cancer risk for three variants in CHEK2, c.349A>G OR 2.26 (95% CI 1.29 to 3.95), c.1036C>T OR 5.06 (95% CI 1.09 to 23.5) and c.538C>T OR 1.33 (95% CI 1.05 to 1.67) (p≤0.017). Evidence for prostate cancer risk was observed for CHEK2 c.1343T>G OR 3.03 (95% CI 1.53 to 6.03, p=0.0006) for African men and CHEK2 c.1312G>T OR 2.21 (95% CI 1.06 to 4.63, p=0.030) for European men. No evidence of association with ovarian cancer was found for any of these variants. Conclusions This report adds to accumulating evidence that at least some variants in these genes are associated with an increased risk of breast cancer that is clinically important.

Abstract Purpose: Expression of the tumor suppressor gene E-cadherin is diminished in lobular breast cancers and has been implicated in epithelial mesenchymal transition. We assessed risk factor associations for breast cancer stratified by low vs. high E-cadherin protein expression in a pooled analysis within the Breast Cancer Association Consortium (BCAC) studies. Methods: E-cadherin tumor tissue staining was performed centrally at the NCI on formalin-fixed paraffin-embedded tissue microarray (TMA) sections representing 6,010 breast cancer patients from 12 US and European BCAC studies. TMAs were digitally scanned and scored using the SlidePath Digital Image Hub (Leica Biosystems, Wetzlar, Germany). For 5,896 cancers with evaluable tumors, E-cadherin was visually scored as estimated percent of positive cells times stain intensity (0, 1+, 2+, 3+) (score range 0-300). E-cadherin low was defined as tumors with a score &lt; 100. Risk factor associations for low vs. high E-cadherin expressing tumors were evaluated by logistic regression, adjusted for age and study site, and stratified by ER status and histologic subtype. To assess the consistency in results by study, meta-analyses were performed using the random-effects model. All statistical tests were two-sided. Results: E-cadherin low cancers comprised 20% of tumors and were associated with lobular histology, well/moderately differentiated cancers, &gt; 2cm in size, and HER2-negative status (χ2, P &lt; 0.003 for all factors). E-cadherin low status was associated with family history of breast cancer (FH) (OR = 1.32, 95% CI = 1.09-1.60, P-het = 0.005) and ever use of menopausal hormones (OR = 1.26, 95% CI = 1.02-1.56, P-het = 0.03). Study specific meta-plots showed consistent effects across studies for menopausal hormone therapy (I2 = 0.0%, p = 0.64); however, E-cadherin loss by FH did show evidence of heterogeneity by study (I2 = 54%, P = 0.03). Differences in E-cadherin expression remained significant for FH, and menopausal hormone use when further adjusted for ER (P &lt; 0.05). We also found relationships with E-cadherin loss to vary by BMI, number of live births, age at first birth, and oral contraceptive use in stratified analysis by ER status and histologic subtype. Conclusion: This large pooled analysis shows that breast cancer risk factor associations may differ by E-cadherin expression independent of ER status, suggesting that it may represent a marker of etiologic heterogeneity. Citation Format: Hisani N. Horne, Hannah Oh, Mark E. Sherman, Maya Palakal, Stephen H. Hewitt, Marjanka Schmidt, Javier Benitez, Roger Milne, Hermann Brenner, Heli Nevanlinna, Arto Mannermaa, Georgia Chenevix-Trench, Fergus Couch, Peter Devilee, Diana Eccles, Maartje Hooning, Anthony J. Swerdlow, Nick Orr, Melissa A. Troester, Renata Cora, Paul D. Pharoah, Montserrat Garcia-Closas, Jonine D. Figueroa. Breast cancer risk factor associations by loss of E-cadherin tumor tissue expression: A pooled analysis of 5,896 cases in 12 studies from the Breast Cancer Association Consortium (BCAC). [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 3451.