Cunningham, Julie M.

Abstract Introduction: Current models to predict breast cancer risk do not differentiate risk for estrogen receptor (ER) positive and negative breast cancer (BC), despite growing evidence that these tumors are biologically very different. We hypothesized that women with ER+ BC cancers have different clinical risk factors and histologic findings on prior benign breast biopsies than those with ER- BC. Methods: After IRB approval, we examined associations of age at benign biopsy and histologic features of the benign biopsy with ER status of incident BCs within the Mayo Benign Breast Disease Cohort. Benign biopsy slides were reviewed for extent of lobular involution and degree of epithelial proliferation by a single breast pathologist blinded to BC events. Invasive BCs occurring within 15 years after benign biopsy were classified as ER+ if ER staining was >1%. BC case-only associations with ER status were evaluated using multivariate logistic regression. Full-cohort hazard ratios (HR) and 95% confidence intervals (CI) for risk of ER-specific subtypes were estimated using Cox proportion hazards regression. Results: Among 13,410 women undergoing a benign breast biopsy from 1967-2001, 656 invasive BCs (459 ER+, 106 ER, 106 unknown) occurred within 15 years. Women who developed ER+ and ER- BCs were similar in age at the time of their prior benign breast biopsy (p=0.34). Although benign biopsies in women who later developed ER+ BC were more likely to show complete involution (23% vs 15% for ER- BC), this was not statistically significant (p=0.06). However, the degree of epithelial proliferation was significantly associated with ER status of later BCs (p=0.001), with ER+ BCs more likely than ER- BCs to have had a prior biopsy with atypical hyperplasia (16% vs 8%), and ER+ BCs less likely than ER- BCs to have had a prior biopsy with proliferative disease without atypia (33% vs 52%); this association remained after multivariate adjustment (p=0.003). We further pursued the association of epithelial proliferation with differential risk of ER+ and ER- BC in our overall cohort of 13,410 women (Table 1). Compared to women with non-proliferative disease, women with proliferative disease +/- atypia had ∼2-fold hazard ratios for ER- BC, whereas hazard ratios for ER+ BC were higher in women with atypical hyperplasia (∼4-fold) compared to proliferative disease. Hazard Ratio of ER+ and ER- BC within 15 years of BBD based upon degree of epithelial proliferation at benign biopsyCancer TypeGroupNN EventsHazard Ratio (95% CI)p-valueInvasive ER- ≤15 yearsNP8478421.00 (ref)<0.0001 PD4229542.63 (1.76, 3.93) AH70382.60 (1.22, 5.54) Invasive ER+ ≤15 yearsNP84782321.00 (ref)<0.0001 PD42291511.33 (1.09, 1.64) AH703744.41 (3.39, 5.73) NP=nonproliferative disease, PD=proliferative disease without atypia; AH=atypical hyperplasia Conclusion: ER+ and ER- breast cancers appear to have different features on prior benign breast biopsy, with atypical hyperplasia showing increased risk for both types of breast cancer, but a greater risk for ER+ tumors. Citation Format: Amy C Degnim, Derek C Radisky, Robert A Vierkant, Ryan D Frank, Marlene H Frost, Vernon S Pankratz, Celine M Vachon, Tanya L Hoskin, Julie M Cunningham, Chen Wang, Jean-Pierre Kocher, Teresa M Allers, Joanne L Johnson, Tina J Hieken, Karthik Ghosh, Lynn C Hartmann, Daniel W Visscher. Histologic features of benign breast biopsy tissue and association with ER positive and ER negative breast cancer in the Mayo BBD cohort study [abstract]. In: Proceedings of the Thirty-Seventh Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2014 Dec 9-13; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2015;75(9 Suppl):Abstract nr P6-10-06.

Abstract Introduction: Peripheral blood leukocyte (PBL) telomere length is associated with colon cancer (CC) risk, but the relationship of PBL telomere length to survival in CC patients receiving chemotherapy is unknown. Methods: We measured relative telomere length (RTL) by qPCR of PBL DNA from 2,748 cases (53% male, 85% white, median age 58 yrs) in the North American Phase III Intergroup Trial N0147 that compared disease-free survival (DFS) among completely resected Stage III CC treated with mFOLFOX6 +/- cetuximab. Tumor markers were assessed, and 12% were BRAF mutant; 35% Kras mutant, and 12% deficient mismatch repair (dMMR). Linear and non-linear prognostic associations of PBL RTL to the primary endpoint of DFS were assessed by Cox models. Optimal categories of RTL were developed using cumulative martingale residual plots, designated from shortest to longest RTL (Table 1). Results: Never smokers had longer PBL RTL (p = 0.0006) than past or current smokers. Those with diabetes mellitus (DM) had longer PBL RTL than those without DM (p = 0.0165). PBL RTL was not associated with any tumor markers (MMR status; p = 0.7174; BRAF mutant vs wild-type (wt); p =0.4039; Kras mutant vs wt; p= 0.1149). PBL RTL was significantly associated with DFS in a non-linear form (p = 0.0389). Cases in the longer RTL category had more favorable DFS than those with shorter or longest RTL, (Table), which remained after adjusting for age, sex , T/N stage, grade, tumor markers, and treatment (p = 0.0230). We observed no significant interaction (p = 0.1490) between RTL categories and treatment. Conclusion: We observed a non-linear association between RTL and DFS in stage III CC, suggesting that an optimal range telomere length can influence tumor behavior in CC patients. Future studies of telomere maintenance mechanisms may help to explain the association of RTL with improved DFS. RTL is associated with smoking and DM in CC cases but further studies are needed to determine if either tobacco use or DM affect survival from CC by modulating PBL telomere length. Citation Format: Nivedita Basu, Ruth A. Johnson, Qian Shi, Garth D. Nelson, Julie M. Cunningham, Richard M. Goldberg, Frank A. Sinicrope, Daniel J. Sargent, Steven R. Alberts, Gloria M. Petersen, Stephen N. Thibodeau, Halcyon G. Skinner, Kristin Litzelman, Songwon Seo, Russell J. Vanderboom, Ronald E. Gangnon, Corinne D. Engelman, David N. Rider, Lisa A. Boardman. Telomeres and survival in stage III colon cancers: Findings from NCCTG N0147 (Alliance) study. [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 5102. doi:10.1158/1538-7445.AM2014-5102

Abstract Mammographic density measures adjusted for age and body mass index (BMI) are heritable predictors of breast cancer risk, but few mammographic density-associated genetic variants have been identified. Using data for 10,727 women from two international consortia, we estimated associations between 77 common breast cancer susceptibility variants and absolute dense area, percent dense area and absolute nondense area adjusted for study, age, and BMI using mixed linear modeling. We found strong support for established associations between rs10995190 (in the region of ZNF365), rs2046210 (ESR1), and rs3817198 (LSP1) and adjusted absolute and percent dense areas (all P < 10−5). Of 41 recently discovered breast cancer susceptibility variants, associations were found between rs1432679 (EBF1), rs17817449 (MIR1972-2: FTO), rs12710696 (2p24.1), and rs3757318 (ESR1) and adjusted absolute and percent dense areas, respectively. There were associations between rs6001930 (MKL1) and both adjusted absolute dense and nondense areas, and between rs17356907 (NTN4) and adjusted absolute nondense area. Trends in all but two associations were consistent with those for breast cancer risk. Results suggested that 18% of breast cancer susceptibility variants were associated with at least one mammographic density measure. Genetic variants at multiple loci were associated with both breast cancer risk and the mammographic density measures. Further understanding of the underlying mechanisms at these loci could help identify etiologic pathways implicated in how mammographic density predicts breast cancer risk. Cancer Res; 75(12); 2457–67. ©2015 AACR.

Key PointsThe NF-κB subunits p52 and RelB increase IRF4 promoter activity and expression in PTCL cells. A positive feedback loop involving CD30, NF-κB, and IRF4 drives PTCL cell proliferation and can be blocked by NF-κB inhibitors.

Abstract Purpose: To analyze the effect of germline mutations in BRCA1 and BRCA2 on mortality in patients with ovarian cancer up to 10 years after diagnosis. Experimental Design: We used unpublished survival time data for 2,242 patients from two case–control studies and extended survival time data for 4,314 patients from previously reported studies. All participants had been screened for deleterious germline mutations in BRCA1 and BRCA2. Survival time was analyzed for the combined data using Cox proportional hazard models with BRCA1 and BRCA2 as time-varying covariates. Competing risks were analyzed using Fine and Gray model. Results: The combined 10-year overall survival rate was 30% [95% confidence interval (CI), 28%–31%] for non-carriers, 25% (95% CI, 22%–28%) for BRCA1 carriers, and 35% (95% CI, 30%–41%) for BRCA2 carriers. The HR for BRCA1 was 0.53 at time zero and increased over time becoming greater than one at 4.8 years. For BRCA2, the HR was 0.42 at time zero and increased over time (predicted to become greater than 1 at 10.5 years). The results were similar when restricted to 3,202 patients with high-grade serous tumors and to ovarian cancer–specific mortality. Conclusions: BRCA1/2 mutations are associated with better short-term survival, but this advantage decreases over time and in BRCA1 carriers is eventually reversed. This may have important implications for therapy of both primary and relapsed disease and for analysis of long-term survival in clinical trials of new agents, particularly those that are effective in BRCA1/2 mutation carriers. Clin Cancer Res; 21(3); 652–7. ©2014 AACR.

ObjectiveMutations in Charcot–Marie–Tooth disease (CMT) genes are the cause of rare familial forms of polyneuropathy. Whether allelic variability in CMT genes is also associated with common forms of polyneuropathy—considered “acquired” in medical parlance—is unknown. Chemotherapy‐induced peripheral neuropathy (CIPN) occurs commonly in cancer patients and is individually unpredictable. We used CIPN as a clinical model to investigate the association of non‐CMT polyneuropathy with CMT genes.MethodsA total of 269 neurologically asymptomatic cancer patients were enrolled in the clinical trial Alliance N08C1 to receive the neurotoxic drug paclitaxel, while undergoing prospective assessments for polyneuropathy. Forty‐nine CMT genes were analyzed by targeted massively parallel sequencing of genomic DNA from patient blood.ResultsA total of 119 (of 269) patients were identified from the 2 ends of the polyneuropathy phenotype distribution: patients that were most and least susceptible to paclitaxel polyneuropathy. The CMT gene PRX was found to be deleteriously mutated in patients who were susceptible to CIPN but not in controls (p = 8 × 10−3). Genetic variation in another CMT gene, ARHGEF10, was highly significantly associated with CIPN (p = 5 × 10−4). Three nonsynonymous recurrent single nucleotide variants contributed to the ARHGEF10 signal: rs9657362, rs2294039, and rs17683288. Of these, rs9657362 had the strongest effect (odds ratio = 4.8, p = 4 × 10−4).InterpretationThe results reveal an association of CMT gene allelic variability with susceptibility to CIPN. The findings raise the possibility that other acquired polyneuropathies may also be codetermined by genetic etiological factors, of which some may be related to genes already known to cause the phenotypically related Mendelian disorders of CMT. Ann Neurol 2014;76:727–737

Abstract Tumor infiltrating regulatory T cells (Tregs) appear to play a role in survival in a number of solid tumors, including ovarian cancer. We assessed genetic variation via tag SNPs (N=749) in 25 Treg associated genes (CD28, CTLA4, FOXP3, IDO1, IL10, IL10RA, IL15, 1L17RA, IL23A, IL23R, IL2RA, IL6, IL6R, IL8, LGALS1, LGALS9, MAP3K8, STAT5A, STAT5B, TGFB1, TGFB2, TGFB3, TGFBR1, TGRBR2, and TGFBR3) in relation to ovarian cancer survival. SNPs were selected based on r2 ≥ 0.8, MAF ≥ 0.05, 5 kb upstream and downstream from the gene. We had SNP and outcome data for ovarian cancer from 10,084 combined ovarian cancer cases, including 5,248 high grade serous, 1,452 endometrioid, 661 mucinous, and 795 clear cell cases of European descent across multiple studies from the Ovarian Cancer Association Consortium (OCAC). Cox proportional hazards regression modeling was used to estimate per-allele hazard ratios (HRs) and 95% confidence intervals (CIs) for associations with overall survival, adjusting for age at diagnosis, population substructure PCs, study site, histology (for the all case analysis), tumor stage, tumor grade and oral contraceptive use. The strongest associations were found for endometrioid histologic subtype and IL2RA SNPs (maximum r2=0.67) rs11256497 [HR=1.42, 95% CI: 1.22-1.64; p=5.7 x 10−6], rs791587 [HR=1.36, 95% CI: 1.17-1.57; p=6.2 x 10−5], rs2476491 [HR=1.40, 95% CI: 1.19-1.64; p=5.6 x 10−5], and rs10795763 [HR=1.35, 95% CI: 1.17-1.57; p=7.9 x 10−5], as well as clear cell ovarian cancer and CTLA4 SNP rs231775 [HR=0.67, 95% CI: 0.54-0.82; p=9.3 x 10−5]. No other associations were observed at p<1 x 10−4. Overall, this study provided evidence that variation in genes related to Tregs, particularly IL2RA and CTLA4, influence ovarian cancer survival. Citation Format: Bridget Charbonneau, Keith Knutson, Robert A. Vierkant, Julie M. Cunningham, Zachary C. Fogarty, Kimberly R. Kalli, Matthew J. Maurer, Ann L. Oberg, Brooke L. Fridley, Claudia C. Preston, Paul D.P. Pharoah, Susan J. Ramus, Catherine M. Phelan, Kunle Odunsi, Kirsten Moysich, Ellen L. Goode. Polymorphisms in regulatory T cell related genes and ovarian cancer survival. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 4844. doi:10.1158/1538-7445.AM2013-4844

Abstract Background. Transmembrane transport (GO:0055085) is the process whereby a solute/hormone/ion/iron is transported from one side of a membrane to the other. Disruption of these processes leads to defects in homeostasis which contributes to cancer risk. We hypothesized that germline single nucleotide polymorphisms (SNPs) in the cellular transport genes are associated with EOC risk and histopathology subtype. Methods. We genotyped 305 SNPs from 120 cellular transport related genes in 14,525 cases and 23,448 controls of European ancestry, 387 African, and 2,388 Asian ancestries from 43 studies in the Ovarian Cancer Association Consortium (OCAC). A custom Illumina iSelect designed for the Collaborative Oncological Gene-environment Study (COGS) was used. For women of European ancestry, both invasive cancers combined and the four main histological subtypes (serous (n=8,369), mucinous (n=943), endometroid (n=2,067) and clear cell (n=1,024) carcinoma) were analyzed, while for women of African and Asian ancestry only the serous subtype was analyzed adjusting for study site and first 5 principal components. Odds ratios and 95% confidence intervals were calculated using unconditional logistic regression. We fitted a log-additive genetic model. Results. In women of European ancestry, the strongest evidence of an association for all invasive EOCs was the iron transport gene, HEPH rs17216603 (OR=0.85, 95%CI=0.78-0.93, P=2.8x10−4), which was also the most significant association for serous histological subtype (OR=0.81, 95%CI=0.73-0.91, P=2.0x10−4). For endometrioid EOC, the strongest association was observed in the glucuronidation gene, UGT1A5 rs11563251 (OR=0.82, 95%CI=0.73-0.92, P=6.6x10−4). For clear cell and mucinous EOC, the thyroid hormone transport protein, SERPINA7 rs1804495 showed the strongest associations (OR=1.06, 95%CI=1.06-1.38, P=0.0041) and (OR=0.85, 95%CI=0.73-0.98, P=0.0407), respectively. For women of African ancestry, the strongest association for serous EOC was the glutathione conjugation gene, MGST rs6488840 (OR=0.55, 95%CI=0.37-0.82, P=0.0035). This SNP was not significant in women of European or Asian ancestry. In women of Asian ancestry, the HEPH rs17216603 (OR=1.45, 95%CI=1.15-1.83, P=0.0019) demonstrated the strongest association, although the directionality of risk was reversed compared to European women. Conclusions. Our results suggest that variation in transmembrane transport genes may influence EOC risk and histopathologic subtype in women of European, African and Asian ancestries, highlighting the differences in the disease etiology across populations. Additional epidemiological and functional studies are warranted to elucidate translational and clinical utility of these associations. Citation Format: Ganna Chornokur, Jonathan Tyrer, Hui-Yi Lin, Gang Han, Xiaotao Qu, Chen Zhihua, Ya-Yu Tsai, Ellen L. Goode, Julie M. Cunningham, Edwin Iversen, Susan J. Ramus, Andrew Berchuk, Joellen M. Schildkraut, Alvaro Monteiro, Simon A. Gayther, Steven A. Narod, Paul Pharoah, Thomas A. Sellers, Catherine M. Phelan. Variation in transmembrane transport genes influence epithelial ovarian cancer risk and histopathologic subtype. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 3644. doi:10.1158/1538-7445.AM2013-3644

Abstract Introduction: Symptom control is an important objective in the care of patients with cancer and a critical goal for long-term survivors of antineoplastic treatments. Commonly used chemotherapies, including paclitaxel, are associated with a toxic effect on the peripheral nervous system in a substantial fraction of patients, termed chemotherapy induced peripheral neuropathy (CIPN), that is often debilitating. CIPN cannot be predicted from clinical parameters in neurologically normal subjects, suggesting a genetic basis as the root cause. Patients with hereditary neuropathy (HN) due to one of ~50 recognized genes are prone to developing severe CIPN even if the hereditary predisposition is clinically silent. Assessment of CIPN has usually been an afterthought of cancer treatment studies using crude symptom measurement instruments. N08C1 and N08CA are recent clinical trials_the latter analyzed only in early 2013_focusing on extensive phenotyping by patient reported outcomes (PRO) to support genetic investigations of paclitaxel-related CIPN. Massively parallel (“nextgen”) DNA sequencing is the genetic technology of choice and represents an as-yet unrealized potential to study cancer symptom control. Methods: Exome sequencing_analysis of 20,794 genes at nucleotide resolution_was performed on blood of CIPN cases (n=73) and controls (n=46) from N08C1. The primary endpoint was the association of HN genes with CIPN. Bioinformatic analysis determined SNV and other genetic variations across 6.2x107 nucleotides. The association of genes with CIPN was statistically examined by the over-dispersion test C-alpha. The role of singleton SNV was assessed by a burden test. N08CA served as validation. Results: 386,644 SNV were detected with an estimated sensitivity of 90% and specificity of 99.9%; the orthogonal validation rate (PCR of >1000 genotypes) was 99.7%. The association of the minor allele (C>T) at rs7349683 in the gene EPHA5, previously reported by Baldwin et al. (2012) for CALGB40101, was confirmed in N08C1. Two hereditary neuropathy genes were identified as new CIPN candidate genes: Firstly, ARHGEF10 was highly significant with p=5x10-4 based on single gene testing and p<0.02 after Bonferroni correction for multiple testing of 41 hereditary neuropathy genes with non-synonymous SNV. Three SNV were predictive, rs17683288, rs2294039 and rs9657362, the latter with an odds ratio (OR)=3.62. Validation in an independent patient group (n=75) from another paclitaxel CIPN trial, N08CA, confirmed the impact of ARHGEF10 on CIPN (p<0.025; rs9657362 OR=2.51). Secondly, peraxin (PRX) was disrupted by non-synonymous singleton SNV in n=8 cases and n=0 controls (OR=6.4; unadjusted p=0.04). 7 of the 8 mutations were also scored as deleterious by both Polyphen 2 and SIFT. Because PRX SNV were non-recurrent validation was performed in N08C1 by Sanger sequencing. Several neuro-genes across the exome were noted as additional CIPN gene candidates. Conclusions: ARHGEF10 and PRX associations with CIPN suggest that CIPN and HN might share genetic roots in a subset of patients. Construct validity of N08C1 is supported by confirmation of a published result from CALGB40101 relating EPHA5 to CIPN. Germline exome sequencing in cancer patients is highly informative for understanding complex pharmacogenomic traits including outcomes in symptom control research. Citation Format: Andreas S. Beutler, Amit Kulkarni, Rahul Kanwar, Rui Qin, Julie M. Cunningham, Terry M. Therneau, Charles L. Loprinzi. Sequencing symptom control: Results from the Alliance N08C1 and N08CA genetics of chemotherapy neuropathy trials. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr LB-196. doi:10.1158/1538-7445.AM2013-LB-196