Pharoah, Paul D.P.

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 ARID1A (BAF250a) is a component of the SWI/SNF chromatin modifying complex, plays an important tumour suppressor role, and is considered prognostic in several malignancies. However, in ovarian carcinomas there are contradictory reports on its relationship to outcome, immune response, and correlation with clinicopathological features. We assembled a series of 1623 endometriosis‐associated ovarian carcinomas, including 1078 endometrioid (ENOC) and 545 clear cell (CCOC) ovarian carcinomas, through combining resources of the Ovarian Tumor Tissue Analysis (OTTA) Consortium, the Canadian Ovarian Unified Experimental Resource (COEUR), local, and collaborative networks. Validated immunohistochemical surrogate assays for ARID1A mutations were applied to all samples. We investigated associations between ARID1A loss/mutation, clinical features, outcome, CD8 + tumour‐infiltrating lymphocytes (CD8 + TILs), and DNA mismatch repair deficiency (MMRd). ARID1A loss was observed in 42% of CCOCs and 25% of ENOCs. We found no associations between ARID1A loss and outcomes, stage, age, or CD8 + TIL status in CCOC. Similarly, we found no association with outcome or stage in endometrioid cases. In ENOC, ARID1A loss was more prevalent in younger patients ( p = 0.012) and was associated with MMRd ( p < 0.001) and the presence of CD8 + TILs ( p = 0.008). Consistent with MMRd being causative of ARID1A mutations, in a subset of ENOCs we also observed an association with ARID1A loss‐of‐function mutation as a result of small indels ( p = 0.035, versus single nucleotide variants). In ENOC, the association with ARID1A loss, CD8 + TILs, and age appears confounded by MMRd status. Although this observation does not explicitly rule out a role for ARID1A influence on CD8 + TIL infiltration in ENOC, given current knowledge regarding MMRd, it seems more likely that effects are dominated by the hypermutation phenotype. This large dataset with consistently applied biomarker assessment now provides a benchmark for the prevalence of ARID1A loss‐of‐function mutations in endometriosis‐associated ovarian cancers and brings clarity to the prognostic significance. © 2021 The Pathological Society of Great Britain and Ireland.

Salicylic acid (SA) has observationally been shown to decrease colorectal cancer (CRC) risk. Aspirin (acetylsalicylic acid, that rapidly deacetylates to SA) is an effective primary and secondary chemopreventive agent. Through a Mendelian randomization (MR) approach, we aimed to address whether levels of SA affected CRC risk, stratifying by aspirin use. A two-sample MR analysis was performed using GWAS summary statistics of SA (INTERVAL and EPIC-Norfolk, N = 14,149) and CRC (CCFR, CORECT, GECCO and UK Biobank, 55,168 cases and 65,160 controls). The DACHS study (4410 cases and 3441 controls) was used for replication and stratification of aspirin-use. SNPs proxying SA were selected via three methods: (1) functional SNPs that influence the activity of aspirin-metabolising enzymes; (2) pathway SNPs present in enzymes’ coding regions; and (3) genome-wide significant SNPs. We found no association between functional SNPs and SA levels. The pathway and genome-wide SNPs showed no association between SA and CRC risk (OR: 1.03, 95% CI: 0.84–1.27 and OR: 1.08, 95% CI: 0.86–1.34, respectively). Results remained unchanged upon aspirin use stratification. We found little evidence to suggest that an SD increase in genetically predicted SA protects against CRC risk in the general population and upon stratification by aspirin use.

Abstract Background: Many loci have been found to be associated with risk of epithelial ovarian cancer (EOC). However, although there is considerable variation in progression-free survival (PFS), no loci have been found to be associated with outcome at genome-wide levels of significance. Methods: We carried out a genome-wide association study (GWAS) of PFS in 2,352 women with EOC who had undergone cytoreductive surgery and standard carboplatin/paclitaxel chemotherapy. Results: We found seven SNPs at 12q24.33 associated with PFS (P &lt; 5 × 10–8), the top SNP being rs10794418 (HR = 1.24; 95% CI, 1.15–1.34; P = 1.47 × 10–8). High expression of a nearby gene, ULK1, is associated with shorter PFS in EOC, and with poor prognosis in other cancers. SNP rs10794418 is also associated with expression of ULK1 in ovarian tumors, with the allele associated with shorter PFS being associated with higher expression, and chromatin interactions were detected between the ULK1 promoter and associated SNPs in serous and endometrioid EOC cell lines. ULK1 knockout ovarian cancer cell lines showed significantly increased sensitivity to carboplatin in vitro. Conclusions: The locus at 12q24.33 represents one of the first genome-wide significant loci for survival for any cancer. ULK1 is a plausible candidate for the target of this association. Impact: This finding provides insight into genetic markers associated with EOC outcome and potential treatment options. See related commentary by Peres and Monteiro, p. 1604

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.

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.

Objective An understanding of the etiologic heterogeneity of colorectal cancer (CRC) is critical for improving precision prevention, including individualized screening recommendations and the discovery of novel drug targets and repurposable drug candidates for chemoprevention. Known differences in molecular characteristics and environmental risk factors among tumors arising in different locations of the colorectum suggest partly distinct mechanisms of carcinogenesis. The extent to which the contribution of inherited genetic risk factors for CRC differs by anatomical subsite of the primary tumor has not been examined. Design To identify new anatomical subsite-specific risk loci, we performed genome-wide association study (GWAS) meta-analyses including data of 48 214 CRC cases and 64 159 controls of European ancestry. We characterised effect heterogeneity at CRC risk loci using multinomial modelling. Results We identified 13 loci that reached genome-wide significance (p<5×10 −8 ) and that were not reported by previous GWASs for overall CRC risk. Multiple lines of evidence support candidate genes at several of these loci. We detected substantial heterogeneity between anatomical subsites. Just over half (61) of 109 known and new risk variants showed no evidence for heterogeneity. In contrast, 22 variants showed association with distal CRC (including rectal cancer), but no evidence for association or an attenuated association with proximal CRC. For two loci, there was strong evidence for effects confined to proximal colon cancer. Conclusion Genetic architectures of proximal and distal CRC are partly distinct. Studies of risk factors and mechanisms of carcinogenesis, and precision prevention strategies should take into consideration the anatomical subsite of the tumour.

Abstract Background: Evidence for aspirin's chemopreventative properties on colorectal cancer (CRC) is substantial, but its mechanism of action is not well-understood. We combined a proteomic approach with Mendelian randomization (MR) to identify possible new aspirin targets that decrease CRC risk. Methods: Human colorectal adenoma cells (RG/C2) were treated with aspirin (24 hours) and a stable isotope labeling with amino acids in cell culture (SILAC) based proteomics approach identified altered protein expression. Protein quantitative trait loci (pQTLs) from INTERVAL (N = 3,301) and expression QTLs (eQTLs) from the eQTLGen Consortium (N = 31,684) were used as genetic proxies for protein and mRNA expression levels. Two-sample MR of mRNA/protein expression on CRC risk was performed using eQTL/pQTL data combined with CRC genetic summary data from the Colon Cancer Family Registry (CCFR), Colorectal Transdisciplinary (CORECT), Genetics and Epidemiology of Colorectal Cancer (GECCO) consortia and UK Biobank (55,168 cases and 65,160 controls). Results: Altered expression was detected for 125/5886 proteins. Of these, aspirin decreased MCM6, RRM2, and ARFIP2 expression, and MR analysis showed that a standard deviation increase in mRNA/protein expression was associated with increased CRC risk (OR: 1.08, 95% CI, 1.03–1.13; OR: 3.33, 95% CI, 2.46–4.50; and OR: 1.15, 95% CI, 1.02–1.29, respectively). Conclusions: MCM6 and RRM2 are involved in DNA repair whereby reduced expression may lead to increased DNA aberrations and ultimately cancer cell death, whereas ARFIP2 is involved in actin cytoskeletal regulation, indicating a possible role in aspirin's reduction of metastasis. Impact: Our approach has shown how laboratory experiments and population-based approaches can combine to identify aspirin-targeted proteins possibly affecting CRC risk.

Abstract Background: It is not known whether modifiable lifestyle factors that predict survival after invasive breast cancer differ by subtype. Methods: We analyzed data for 121,435 women diagnosed with breast cancer from 67 studies in the Breast Cancer Association Consortium with 16,890 deaths (8,554 breast cancer specific) over 10 years. Cox regression was used to estimate associations between risk factors and 10-year all-cause mortality and breast cancer–specific mortality overall, by estrogen receptor (ER) status, and by intrinsic-like subtype. Results: There was no evidence of heterogeneous associations between risk factors and mortality by subtype (Padj &gt; 0.30). The strongest associations were between all-cause mortality and BMI ≥30 versus 18.5–25 kg/m2 [HR (95% confidence interval (CI), 1.19 (1.06–1.34)]; current versus never smoking [1.37 (1.27–1.47)], high versus low physical activity [0.43 (0.21–0.86)], age ≥30 years versus &lt;20 years at first pregnancy [0.79 (0.72–0.86)]; &gt;0–&lt;5 years versus ≥10 years since last full-term birth [1.31 (1.11–1.55)]; ever versus never use of oral contraceptives [0.91 (0.87–0.96)]; ever versus never use of menopausal hormone therapy, including current estrogen–progestin therapy [0.61 (0.54–0.69)]. Similar associations with breast cancer mortality were weaker; for example, 1.11 (1.02–1.21) for current versus never smoking. Conclusions: We confirm associations between modifiable lifestyle factors and 10-year all-cause mortality. There was no strong evidence that associations differed by ER status or intrinsic-like subtype. Impact: Given the large dataset and lack of evidence that associations between modifiable risk factors and 10-year mortality differed by subtype, these associations could be cautiously used in prognostication models to inform patient-centered care.

Abstract Improving risk prediction and prevention strategies through identifying new susceptibility genes for non-high-grade serous ovarian cancer (non-HGS) would represent an important advance in reducing incidence and mortality. Epithelial ovarian cancer (EOC) has five main histotypes that have distinct pathologies, molecular changes, clinical characteristics, and tissues of origin. They are classified as high-grade serous (HGS), low-grade serous (LGS), clear-cell (CCC), endometrioid (END), and mucinous (MUC). A family history of breast or ovarian cancer is the strongest single risk factor for EOC. Germline mutations in the high-penetrance susceptibility genes BRCA1 and BRCA2 confer EOC risks of 44% and 17% by age 80, respectively. Only 40% of the excess familial risk is known, suggesting there are many other susceptibility genes yet unidentified. The non-HGS histotypes are poorly studied due to limited sample size within individual studies. Thus, families of cases with non-HGS tumors have few risk prediction options. Whole-exome sequencing (WES) analysis was performed on 251 non-HGS cases (56 LGS, 55 CCC, 117 END, 23 MUC). Cases were screened negative for BRCA1/2 mutations and selected for a family history of ovarian or breast cancer, or young onset (&lt;45 years) of ovarian cancer. WES identified 1,278 genes with rare predicted truncating mutation in at least one case. Predicted truncating mutations in cases were compared to 171,584 controls from the Exome Aggregation Consortium (ExAC) and the Genome Aggregation Database (gnomAD). Gene-set enrichment analysis showed enrichment for genes involved in the DNA repair pathway (p=3.08 × 10−7). Twenty-five candidate genes (including 4 DNA repair genes) were selected for validation by targeted sequencing. Five genes important in HGS (BRIP1, FANCM, PALB2, RAD51C, RAD51D) as well as two additional candidates (XRCC2 and XRCC3) were also sequenced. Targeted sequencing was performed on 1,779 cases (669 END, 356 CCC, 327 LGS, 427 MUC) and 1,863 controls from studies in the Ovarian Cancer Association Consortium (OCAC). Library preparation and target enrichment was performed using the Fluidigm Juno System and Illumina sequencing was performed on a NovaSeq 6000 Sequencing System. A case-control analysis of predicted truncating variants in the 32 genes was performed. Using a minimum alternate allele frequency of 30%, we identified a higher frequency of mutations in non-HGS cases than controls in ERCC6 (p=0.034) and IL31RA (p=0.034). ERCC6 is involved in DNA repair and IL31RA is a cytokine receptor. Sanger sequencing validation of variants to confirm these results is ongoing. Several other genes showed suggestive higher frequencies of mutations in cases than controls. A larger case control analysis will be performed to confirm those findings. Identifying novel susceptibility genes for non-HGS may have clinical impact by reducing disease-associated mortality through improving risk prediction, identifying prevention strategies, and developing new targeted treatments. Citation Format: Marina Pavanello, Ed Dicks, Honglin Song, Amir Ariff, Adelyn Bolithon, Maria P. Intermaggio, Mark Pinese, Kirsten Moysich, Kunle O. Odunsi, Ellen Goode, David D. Bowtell, Peter Fasching, Jennifer A. Doherty, Francesmary Modugno, Susanne K. Kjær, Penelope M. Webb, Anna Wu, Anna deFazio, Ovarian Cancer Association Consortium, Paul James, Deepak Subramanian, Ian Campbell, Simon A. Gayther, Paul D.P. Pharoah, Susan J. Ramus. Germline mutations in new susceptibility genes for non-high-grade serous ovarian cancer [abstract]. In: Proceedings of the AACR Special Conference on Advances in Ovarian Cancer Research; 2019 Sep 13-16, 2019; Atlanta, GA. Philadelphia (PA): AACR; Clin Cancer Res 2020;26(13_Suppl):Abstract nr B37.