Couch, Fergus J.

Abstract Many somatic mutations have been detected in pancreatic ductal adenocarcinoma (PDAC), leading to the identification of some key drivers of disease progression, but the involvement of large genomic rearrangements has often been overlooked. In this study, we performed mate pair sequencing (MPseq) on genomic DNA from 24 PDAC tumors, including 15 laser-captured microdissected PDAC and 9 patient-derived xenografts, to identify genome-wide rearrangements. Large genomic rearrangements with intragenic breakpoints altering key regulatory genes involved in PDAC progression were detected in all tumors. SMAD4, ZNF521, and FHIT were among the most frequently hit genes. Conversely, commonly reported genes with copy number gains, including MYC and GATA6, were frequently observed in the absence of direct intragenic breakpoints, suggesting a requirement for sustaining oncogenic function during PDAC progression. Integration of data from MPseq, exome sequencing, and transcriptome analysis of primary PDAC cases identified limited overlap in genes affected by both rearrangements and point mutations. However, significant overlap was observed in major PDAC-associated signaling pathways, with all PDAC exhibiting reduced SMAD4 expression, reduced SMAD-dependent TGFβ signaling, and increased WNT and Hedgehog signaling. The frequent loss of SMAD4 and FHIT due to genomic rearrangements strongly implicates these genes as key drivers of PDAC, thus highlighting the strengths of an integrated genomic and transcriptomic approach for identifying mechanisms underlying disease initiation and progression. Cancer Res; 76(3); 749–61. ©2015 AACR.

Abstract The prevalence of germline pathogenic mutations in a comprehensive panel of cancer predisposition genes is not well-defined for patients with pancreatic ductal adenocarcinoma (PDAC). To estimate the frequency of mutations in a panel of 22 cancer predisposition genes, 96 patients unselected for a family history of cancer who were recruited to the Mayo Clinic Pancreatic Cancer patient registry over a 12-month period were screened by next-generation sequencing. Fourteen pathogenic mutations in 13 patients (13.5%) were identified in eight genes: four in ATM, two in BRCA2, CHEK2, and MSH6, and one in BARD1, BRCA1, FANCM, and NBN. These included nine mutations (9.4%) in established pancreatic cancer genes. Three mutations were found in patients with a first-degree relative with PDAC, and 10 mutations were found in patients with first- or second-degree relatives with breast, pancreas, colorectal, ovarian, or endometrial cancers. These results suggest that a substantial proportion of patients with PDAC carry germline mutations in predisposition genes associated with other cancers and that a better understanding of pancreatic cancer risk will depend on evaluation of families with broad constellations of tumors. These findings highlight the need for recommendations governing germline gene-panel testing of patients with pancreatic cancer. Cancer Epidemiol Biomarkers Prev; 25(1); 207–11. ©2015 AACR.

Abstract Background: A recent association study identified a common variant (rs9790517) at 4q24 to be associated with breast cancer risk. Independent association signals and potential functional variants in this locus have not been explored. Methods: We conducted a fine-mapping analysis in 55,540 breast cancer cases and 51,168 controls from the Breast Cancer Association Consortium. Results: Conditional analyses identified two independent association signals among women of European ancestry, represented by rs9790517 [conditional P = 2.51 × 10−4; OR, 1.04; 95% confidence interval (CI), 1.02–1.07] and rs77928427 (P = 1.86 × 10−4; OR, 1.04; 95% CI, 1.02–1.07). Functional annotation using data from the Encyclopedia of DNA Elements (ENCODE) project revealed two putative functional variants, rs62331150 and rs73838678 in linkage disequilibrium (LD) with rs9790517 (r2 ≥ 0.90) residing in the active promoter or enhancer, respectively, of the nearest gene, TET2. Both variants are located in DNase I hypersensitivity and transcription factor–binding sites. Using data from both The Cancer Genome Atlas (TCGA) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC), we showed that rs62331150 was associated with level of expression of TET2 in breast normal and tumor tissue. Conclusion: Our study identified two independent association signals at 4q24 in relation to breast cancer risk and suggested that observed association in this locus may be mediated through the regulation of TET2. Impact: Fine-mapping study with large sample size warranted for identification of independent loci for breast cancer risk. Cancer Epidemiol Biomarkers Prev; 24(11); 1680–91. ©2015 AACR.

Abstract Panel-based germline genetic testing of established breast and ovarian cancer predisposition genes is widely used in high-risk clinics to inform on risks of breast and ovarian cancer. While most protein truncating mutations identified in predisposition genes are associated with moderate to high risks of these cancers, variants of uncertain significance (VUS), in the form of missense and intronic variants, have an undefined relevance to cancer. Many hundreds of VUS in BRCA1 and BRCA2 identified by clinical genetic testing over the last 20 years still have uncertain clinical relevance. In addition, VUS in other known predisposition genes are being identified by current panel-testing in up to 10% of individuals without BRCA1 or BRCA2 mutations. These VUS cause problems in cancer risk estimation and clinical management of patients undergoing clinical genetic testing. In an effort to evaluate the clinical significance of VUS in cancer predisposition genes we have developed functional assays that characterize the impact of the variants on the functions of the proteins encoded by the predisposition genes. In particular, a homology directed repair reporter assay providing a quantitative measure of homologous recombination repair of DNA double strand breaks has been developed for characterization of VUS in homologous recombination associated genes. When applied to 13 known pathogenic and 30 known neutral variants from the DNA binding domain of BRCA2, this assay exhibits 100% specificity and 100% sensitivity. Based on these results a model estimating the probability of pathogenicity of each BRCA2 VUS was developed. The HDR assay and probability model has now been applied to 230 BRCA2 VUS identified through genetic testing. Here we identify 80 BRCA2 variants that are deficient in HDR activity. We further combined these results with prior probabilities of pathogenicity derived from the sequence based Align GVGD in silico prediction model to estimate the posterior probability of pathogenicity for each VUS. Our results suggest that in the absence of family based segregation information, results from the HDR functional assay and in silico prediction models can be used to predict the pathogenicity of BRCA2 VUS. The assay has also been used to characterize VUS in the PALB2 and RAD51C predisposition genes, leading to identification of a subset of missense variants that are deficient in HDR repair activity and may be associated with increased risk of breast and ovarian cancer. These results will prove useful for improved risk assessment of carriers of these BRCA2, PALB2 and RAD51C variants and their family members and for validation of other approaches to mutation classification. Citation Format: Fergus J. Couch, Hermela Shimelis, Andreas Schroeder, Chunling Hu, Emily Hallberg, Gary Lipton, Edwin Iversen, Noralane M. Lindor. Characterization of variants of uncertain significance (VUS) in breast and ovarian cancer predisposition genes. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4685. doi:10.1158/1538-7445.AM2015-4685

Abstract Genome-wide association studies (GWAS) have identified several susceptibility regions that confer risk of two or more cancers. A region at 19p13 is associated with both breast and ovarian cancer risk in the general population as well as in BRCA1/BRCA2 mutation carriers. Using the iCOGs custom genotyping array, we analyzed 438 SNPs in this region in 7,435 ER-negative breast cancer cases and 42,599 controls from BCAC, and 15,347 ovarian cancer cases and 30,845 controls from OCAC and 15,252 BRCA1 mutation carriers. We identified two significant peaks of association. Peak 1, consisting of 13 highly correlated SNPs (r2>0.95), spans a 19.4kb region harboring the BABAM1, ANKLE1 and ABHD8 genes. Peak 1 SNPs are associated with serous ovarian cancer (top hit rs4808075, P = 9.2×10-20), ER-negative breast cancer (rs61494113, P = 1.1×10-13) and breast cancer in BRCA1 mutation carriers (rs61494113, P = 7.7×10-16). Peak 2 spans over 100kb and contains 7 SNPs associated with breast cancer in BRCA1 mutation carriers only (rs4464206, P = 8.9×10-8). Multiple assays were used to establish the putative functional roles of risk SNPs and candidate susceptibility genes in this region. Expression quantitative trait locus (eQTL) analyses identified significant genotype-gene expression correlations between genotype at the Peak 1 SNP rs480816 and ABHD8 expression in ovarian tumors (P = 3.0×10-5) and normal breast tissues (P = 2.8×10-3). The Peak 2 SNP, rs7246243, was significantly associated with BABAM1 expression (P = 5×10-3) in normal breast tissue only. Overexpression of BABAM1 in normal ovarian and breast cells did not modulate neoplastic phenotypes, whereas overexpression of ABHD8 induced a significant reduction in migration in both cell types in vitro. Differentially expressed genes in ABHD8 models, identified by RNAseq, were associated with cellular movement (P<1.82×10-44). Chromosome conformation capture (3C) identified an interaction between the ABHD8 promoter and a region containing 4 risk SNPs (rs4808075, rs10419397, rs56069439 and rs4808076) in breast and ovarian normal and cancer cells but only rs56069439 coincides with epigenetic enhancer marks in breast and ovarian cells. CRISPR/Cas9 genome editing was used to delete a 57bp region containing rs56069439 and expression of BABAM1, ANKLE1 and ABHD8 were measured. Taken together, these data suggest a common functional mechanism underlying risk of breast and ovarian cancer at this locus, mediated by rs56069439 variants differentially regulating expression of ABHD8, a gene involved in cell motility. Citation Format: Kate Lawrenson, Siddhartha Kar, Karoline Kuckenbaeker, Stacey Edwards, Qiyuan Li, Jonathan Tyrer, Jonathan Beesley, Kyriaki Michailidou, Susan Ramus, Alison Dunning, Fergus Couch, Alvaro Monteiro, Jacques Simard, Matthew Freedman, Douglas Easton, Georgia Chenevix-Trench, Paul Pharoah, Antonis Antoniou, Simon Gayther, OCAC, BCAC and CIMBA. Common functional mechanisms underlying pleiotropy at the 19p13.1 breast and ovarian cancer cusceptibility locus. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 2783. doi:10.1158/1538-7445.AM2015-2783

Abstract Distinct subtypes of triple negative (TN) breast cancer have currently been identified by tumor expression profiling and are thought to be influenced by the tumor microenvironment. However, little is known about the relationship between distinguishing histopathologic features of TN tumors, which reflect aspects of both tumor behavior and tumor microenvironment, and molecular TN subtypes. We analyzed whole genome expression profiles of 593 TN tumors from the Triple Negative Breast Cancer Consortium (TNBCC) using the Illumina Whole Genome DASL array. Histopathologic review of mitotic index, grade, trabecular pattern, growth pattern, involution, presence of ductal carcinoma in situ, and degree of lymphoid infiltrate, necrosis, fibrosis, and lobulitis was completed for 310 of these tumors. Clusters were identified both agnostically and using genes associated with these histopathologic features, and subtypes were evaluated for enriched pathways and associations with histopathologic features. Agnostic clustering identified four stable clusters which represented a basal-like (BL) signature, a luminal signature with overexpression of androgen receptor and targets, an immune signature, and a stromal (STR) signature. Overall, the TN tumors tended to have little or no necrosis (91%), high grade (88%), a high mitotic index (80%), fibrosis (77%), an epicentric growth pattern (62%), incomplete lobular involution (55%), trabecular structures (37%), lymphocytic infiltration (25%), adjacent ductal carcinoma in situ (24%), and variable amounts of lobulitis. Using a combined total of 2,776 unique probes individually associated with the ten histopathologic features measured, we identified six distinct TN subtypes, with an additional BL and STR signature compared to the four agnostic subtypes. Our findings suggest that histopathologic features that reflect heterogeneity in tumor cell structure and behavior as well as the tumor microenvironment, particularly the presence of lymphocytic infiltrate, are critical for TN subtype identification. Further work is needed to validate these signatures and to better understand the immunologic pathways associated with each subtype. Citation Format: Kristen Purrington, Daniel Visscher, Drakoulis Yannoukakos, Jane Carpenter, Heli Nevanlinna, Arto Mannerma, Xianshu Wang, Graham Giles, Wei Zheng, Angela Cox, Hiltrud Brauch, Ute Hamann, Diana Eccles, Celine M. Vachon, Fergus J. Couch. Genes associated with histopathologic features of triple negative breast tumors predict molecular subtypes. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 1924. doi:10.1158/1538-7445.AM2015-1924

We used a target-centric strategy to identify transporter proteins upregulated in pancreatic ductal adenocarcinoma (PDAC) as potential targets for a functional imaging probe to complement existing anatomical imaging approaches. We performed transcriptomic profiling (microarray and RNASeq) on histologically confirmed primary PDAC tumors and normal pancreas tissue from 33 patients, including five patients whose tumors were not visible on computed tomography. Target expression was confirmed with immunohistochemistry on tissue microarrays from 94 PDAC patients. The best imaging target identified was SLC6A14 (a neutral and basic amino acid transporter). SLC6A14 was overexpressed at the transcriptional level in all patients and expressed at the protein level in 95% of PDAC tumors. Very little is known about the role of SLC6A14 in PDAC and our results demonstrate that this target merits further investigation as a candidate transporter for functional imaging of PDAC.

Abstract Purpose: Mammographic breast density is an established and strong risk factor for breast cancer. Recently, a number of common genetic susceptibility loci have been identified as risk factors for breast cancer. We present the first report on the relative contribution of 76 validated breast cancer susceptibility loci, in the context of a polygenic risk score (PRS), to the breast density and breast cancer association. We also examine whether the PRS improves prediction of the Breast Cancer Surveillance Consortium (BCSC) 5-year risk model above and beyond breast density and clinical factors included in the model. Methods: The study population included 1643 cases and 2397 controls from three independent epidemiologic studies: the Mayo Mammography Health Study (MMHS), the Mayo Clinic Breast Cancer Study (MCBCS), and the Bavarian Breast Cancer Cases and Control Study (BBCC). Data collected on patients in each of the studies included clinical risk factor, BI-RADS breast density [a) almost entirely fat; b) scattered fibroglandular densities; c) heterogeneously dense; and d) extremely dense] and genotypes for the 76 breast cancer susceptibility loci known at the time of the study. We formed a PRS from the reported per-SNP odds ratios (OR) for 76 known breast cancer susceptibility loci, and evaluated whether BI-RADS density and the PRS were independent risk factors for breast cancer, when adjusted for age, 1/BMI and study. We also incorporated the PRS (OR) into the BCSC 5-year risk model and estimated 5-year risk for the MMHS nested case-control study of 339 invasive cases, 765 controls. Results: BI-RADS density (p<0.0001) and PRS (p<0.0001) were independent risk factors for breast cancer that together showed greater discrimination of risk (AUC=0.69) than density (AUC=0.66; ΔAUC=0.029) or PRS score alone (AUC=0.68; ΔAUC=0.013; p<0.001). Relative to those with scattered fibroglandular densities and average (2nd quartile) PRS, women with extremely dense breasts and in the highest PRS quartile, had a 2.7 fold (95%CI: 1.7-4.1) increased risk of breast cancer, while those with fatty breasts and in the lowest PRS quartile had a reduced risk (OR=0.30, 95%CI: 0.18-0.51). Incorporation of the PRS into the BCSC risk model improved discrimination (ΔAUC=0.031, p=0.001), for a net reclassification improvement of 20% (95%CI: 11%-28%), split equally among cases (9%) and controls (11%). Conclusion: BI-RADS density and the PRS are both important risk factors for breast cancer that can be included in breast cancer risk models to improve prediction of this disease. Using these models to identify high and low-risk risk groups will facilitate improved tailored screening and primary prevention interventions. Citation Format: Celine M. Vachon, V. Shane Pankratz, Christopher G. Scott, Lothar Haeberle, Elad Ziv, Matthew R. Jensen, Kathleen R. Brandt, Dana H. Whaley, Janet E. Olson, Katharina Heusinger, Carolin C. Hack, Sebastian M. Jud, Matthias W. Beckmann, Jeffrey A. Tice, Kristen S. Purrington, Thomas A. Sellers, Karla Kerlikowske, Peter A. Fasching, Fergus J. Couch. The contribution of common breast cancer susceptibility loci to the breast density and breast cancer association and the Breast Cancer Surveillance Consortium (BCSC) risk model. [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 940. doi:10.1158/1538-7445.AM2014-940