Visscher, DW

Abstract Background: Despite a paucity of data on the benign breast transcriptome, reduction mammoplasty tissue is commonly used as a control for RNA biomarker discovery in breast cancer (BC). We evaluated the transcriptome of benign breast to characterize pre-menopausal and post-menopausal coding and non-coding RNA profiles. These profiles will provide insight into the metabolic landscape of benign breast and inform transcriptomic studies of BC and peri-tumoral microenvironments. Methods: With institutional approval, cryobanked fresh breast tissues from reduction mammoplasties (age 15-62y; follow-up: 4-10 y with no cancer events), with histologically-confirmed benign tissues with similar epithelial: stromal tissues were grouped into 4 age sets (< 30y, N=10; 30-39y, N=11; 40-49y, N=10; ≥50y (clinically post-menopausal), N=10). Following RNA extraction, library preparation (Illumina TruSeq Stranded mRNA kit), and sequencing (Illumina HiSeq 4000), reads were processed (MAP-RSeq v3.0.0) and aligned (STAR aligner; hg38). Differential expression (DE) analysis (edgeR 2.6.2) identified DE genes from normalized RPKM reads (absolute log2 fold change (FC) > 1 and false discovery rate (FDR) < 0.10), corrected for intra-group biases using medians (absolute FC cut-off of >1.5). Over-representation analysis [Ingenuity pathway analysis (IPA), Ingenuity® Systems] and gene set enrichment analysis [(GSEA), GeneTrail 2.0] identified significantly-enriched pathways. Results: Across the 4 age sets, 561 DE genes were identified. Compared to the post-menopausal (PM) set, the number of DE genes was highest in <30 y set (N= 372) and decreased with increasing age (N= 170, age 30-39 set and N=20, age 40-49 set), generating up-regulated (PMup) and down-regulated (PMdown) PM transcriptomic profiles. The top PMdown DE genes included RANKL, WNT4, MKI67, extracellular matrix, and lactation-related genes (lactoferrin ,MUC4, MUC16, p < 0.01). Top PMdown canonical pathways were cell cycle-related (CDK1, CCNA2, CCNB2, ESPL1, TOP2A)(p< 0.001). Top PMup genes included those involved in adipogenesis, NPY1R, NPY2R, unsaturated fatty acid synthesis and eicosanoid signaling (P< 0.001). Top PMup canonical pathways included acyl co-A hydrolysis, stearate biosynthesis, acute phase response and RXR signaling (P < 0.001). GSEA of the entire gene set (N=15,466; ranked in order of 561 DE genes) identified 16 significant pathways, functionally grouped as mitochondrial energy metabolism, proteasome, and unsaturated fatty acid biosynthesis and signaling (PPAR pathway). While coding RNAs comprised 85% of DE transcripts, lncRNA comprised ˜5%. Top DE lncRNAs and precursor miRNAs included MEG3, a tumor suppressor, LIN00092, miR22 and miR1182 (p < 0.0001). Conclusions: Pre and post-menopausal benign breast tissues have distinctive transcriptomic profiles. PM benign breast has down-regulated proliferation/cell-cycle-related pathways, and up-regulated genes in mitochondrial energy pathways, lipogenesis and inflammatory pathways. Notably, many of the top DE coding and non-coding RNAs in PM benign breast have been implicated in BC progression, highlighting their value in better understanding breast carcinogenesis and the need to characterize their functional roles in normal aging and menopause. Citation Format: Carter JM, Nair AA, Davila JI, Heinzen EP, Hoskin TL, Winham SJ, Radisky DC, Degnim AC, Visscher DW. Distinctive coding and non-coding RNA profiles of pre-menopausal and post-menopausal benign breast [abstract]. In: Proceedings of the 2018 San Antonio Breast Cancer Symposium; 2018 Dec 4-8; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2019;79(4 Suppl):Abstract nr P5-06-01.

Abstract Background: Differences in ER+ and ER- breast cancer tumor biology are well-documented, but little is known about the field of background benign breast in which those cancers arise. We evaluated the transcriptome of benign breast tissues from women with concurrent ipsilateral ER+ and ER- breast cancers (BC) to characterize coding and non-coding RNA profiles. This pilot study provides insight into the transcriptomic landscapes of benign breast tissues in patients with BC, and the microenvironment of the at-risk benign tissue. Methods: With institutional approval, cryobanked breast tissues from patients with concurrent ipsilateral ER+ BC (benign ER+BC, N=14) or ER- BC (benign ER-BC, N=10), were selected for benign tissues with similar epithelial:stromal ratios and grouped into pre (preM) and post-menopausal (PM) groups. Following RNA sequencing (Illumina TruSeq Stranded mRNA kit & Illumina HiSeq 4000), reads were processed (MAP-RSeq v3.0.0) and aligned (STAR aligner; hg38). Differential expression (DE) analysis (edgeR 2.6.2) identified DE genes from normalized RPKM reads (absolute log2 fold change (FC) > 1 and false discovery rate (FDR) < 0.10), corrected for intra-group biases. Over-representation analysis [Ingenuity pathway analysis (IPA), Ingenuity® Systems] and gene set enrichment analysis [(GSEA), GeneTrail 2.0] identified significantly-enriched pathways. Results: In the PM group, there were 144 DE transcripts between benign ER+BC and benign ER-BC, including coding RNAs (40%), antisense RNAs (35%) and lncRNAs (7%). In contrast, the preM group had no significantly DE genes between benign ER+BC and benign ER-BC. In the PM DE coding gene set, the top DE transcripts in benign ER+BC included many genes implicated in BC development# or ER+ BC progression (*) (e.g. up-regulated: KAAG1*, DNAJB7/HSP40*, TMEM151B, ZBTB32# p < 0.001 ; down-regulated: CPB1*, FOS, TPPP3#, CLEC3B#p < 0.001). Top canonical pathways altered in benign ER+BC included MAPK, PI3K, and acute phase response pathways (p<0.05). GSEA of the entire gene set (N=15,223; ranked in order of 144 DE genes) identified 72 altered pathways (P < 0.005); those with the highest normalization enrichment scores (NES) (> 0.4) functionally grouped as immune function-related (T cell function and antigen presentation). Depleted pathways with NES > 0.4, (N=6) functionally grouped into proteasome-related, fatty acid biosynthesis and mitochondrial energy metabolism. Among the non-coding DE gene set, notably, the entire DE antisense RNA gene set (N=51 transcripts) was up-regulated in benign ER+ BC compared to benign ER-BC (P< 0.001) with a subset (N=11) showing marked up-regulation (> 4 log2FC). Among the DE antisense RNAs, 70% have reported roles in carcinogenesis or BC progression (e.g. KRT7-AS, NAV2-AS2, CCDC144NL-AS1, RP11-66B24.4, HSA-MIR4454). Conclusion: The benign breast transcriptome differs between postmenopausal women with ER+ vs ER- BC, with distinctive coding and non-coding RNA signatures. In postmenopausal women with ER+ BC, benign breast expresses a unique antisense RNA set and is enriched in genes implicated in BC development or progression. These data provide insight into at-risk benign breast and facilitate identification of potential biomarkers of carcinogenesis. Citation Format: Carter JM, Nair AA, Davila JI, Heinzen EP, Hoskin TL, Winham SJ, Radisky DC, Visscher DW, Degnim AC. A unique coding and non-coding benign breast transcriptome in post-menopausal ER+ breast cancer [abstract]. In: Proceedings of the 2018 San Antonio Breast Cancer Symposium; 2018 Dec 4-8; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2019;79(4 Suppl):Abstract nr P3-08-10.

Abstract Background: The LAR subtype is a genomically distinct subset of TNBC. Using a large cohort of non-metastatic TNBC patients (pts) with long term follow-up, we sought to further characterize the clinicopathologic features and outcomes of LAR vs non-LAR TNBC. Methods: From a cohort of 9982 women with surgically-treated non-metastatic breast cancer, 605 met criteria for TNBC (ER/PR<1% and HER2-negative) by central pathology. RNA extracted from 304 FFPE tumor specimens using the HighPure RNA extraction kit was subjected to TruSeq RNA Access library preparation and sequencing on a HiSeq2500. Adequate RNA was available for 283 pts. Tumors were classified as LAR or non-LAR using a shrunken centroid model, CABAL (Clustering Among BAsal and Luminal androgen receptor). In addition to previously described analyses [Leon-Ferre et al, Breast Cancer Res Treat 2017], immunohistochemical (IHC) androgen receptor (AR) staining was performed and the impact of various parameters on invasive disease-free survival (IDFS) and overall survival (OS) was assessed using Cox proportional hazards models. Results: 58 (20%) tumors were classified as LAR and 225 (80%) as non-LAR. Compared to non-LAR, LAR pts were older (mean age 65 vs 54) and more often postmenopausal (79%vs53%), both p=0.01. Apocrine histology was more common among LAR tumors (21%vs0%), which were also lower grade (grade3: 69%vs95%) and had lower Ki-67 (Ki-67>15%: 64%vs82%), all p<0.01. Additionally, LAR tumors had lower median stromal tumor infiltrating lymphocytes (TILs, 20%vs25%) and were less frequently lymphocyte-predominant [≥50% stromal or intratumoral TILs (19%vs32%)], although neither reached statistical significance. AR IHC was available for 223 of 283 tumors. Median AR IHC score in LAR was 65% (range 0-100%) vs 0% (range 0-90%) in non-LAR. T/N stage, surgery type, and receipt of adjuvant chemotherapy (AdjCT) or radiotherapy were similar between LAR and non-LAR. LAR pts had shorter IDFS and OS compared to non-LAR (5.6 vs 11.8 yrs and 10.8 vs 20.8 yrs, respectively), although this did not reach statistical significance. Test of proportional hazard assumption was not significant for IDFS or OS (p = 0.30 and 0.09). IDFS estimates were numerically higher in LAR vs non-LAR (80.2%vs70.5%,p = 0.92) at 3yrs post-diagnosis; whereas the opposite was true (40.9%vs55.6%,p = 0.07) after 10yrs. OS estimates at 3 and 5yrs were similar between LAR and non-LAR, but at 10yrs OS was inferior in LAR (40.9%vs66.4%,p = 0.24). In a univariate analysis including both LAR and non-LAR, older age, higher N stage, lower TILs and absence of AdjCT were associated with poorer IDFS and OS. In a multivariate analysis, higher N stage and absence of AdjCT remained associated with both poorer IDFS and OS; while lower stromal TILs were associated with poorer IDFS (p=0.01), and with a trend towards poorer OS (p=0.07). Conclusions: LAR TNBCs occurred in older women, were lower grade, and had lower TIL density than nonLAR tumors. While significant differences in IDFS or OS were not demonstrated, LAR pts exhibited a numerically lower risk of a disease event at 3yrs, but higher risk by 10yrs compared to nonLAR pts. In the entire cohort, higher N stage, absence of AdjCT and lower TILs were independently associated with poorer outcomes. Citation Format: Leon-Ferre RA, Polley M-Y, Liu H, Kalari KR, Boughey JC, Liu MC, Cafourek V, Negron V, Ingle JN, Thompson KJ, Tang X, Barman P, Carlson E, Visscher DW, Carter JC, Couch FJ, Goetz MP. Characteristics, outcomes and prognostic factors of luminal androgen receptor (LAR) triple-negative breast cancer (TNBC) [abstract]. In: Proceedings of the 2018 San Antonio Breast Cancer Symposium; 2018 Dec 4-8; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2019;79(4 Suppl):Abstract nr P3-08-01.

Abstract Background: While TNBC remains the most aggressive type of breast cancer (BC), substantial heterogeneity in biology and outcomes exists among TNBC subtypes. Historically, risk stratification of TNBC has been based on anatomic factors such as tumor size, nodal involvement and presence of distant metastases. However, these features alone fail to accurately predict outcomes. Tumor immune infiltration (sTILs) and distribution of immune cell subsets in the perip heral blood (NLR) have emerged as variables reported to be associated with outcomes in TNBC. We sought to evaluate whether NLR and sTILs provided independent prognostic information in TNBC. Methods: From a cohort of 9,982 women who underwent BC surgery at Mayo Clinic, Rochester, MN between Jan 1985 and Dec 2012, we identified 605 centrally-confirmed TNBC tumors. Patients (pts) with prior BC, bilateral BC, non-invasive disease, stage IV, neoadjuvant therapy, endocrine therapy, or adenoid cystic histology were excluded. For eligible tumors, clinical and pathologic variables were evaluated, including peripheral blood NLR and central assessment of sTILs per the 2014 International TILs Working Group recommendations. We calculated the Pearson correlation coefficient (PCC) between NLR and sTILs and constructed Cox Proportional Hazards Models to evaluate their association with invasive-disease free (IDFS) and overall survival (OS). NLR and sTILs were both analyzed as continuous variables. Results: Most pts had T1-2 (95%) and N0-1 disease (86%). Median OS follow-up was 10.6yrs. Median IDFS was 12yrs (95%CI 10.2-16.7) and median OS was 18.8yrs (95%CI 15.6-20.8). NLR and sTILs were available in 408 and 599 pts, respectively. The median NLR and sTIL content were 2.29 (0.14-10.50) and 20% (0-90%), respectively. NLR and sTILs were poorly correlated (PCC 0.0237). On univariate analysis (UVA), a higher NLR was associated with worse IDFS (HR 1.13; 95%CI 1.02-1.26, p=0.02) and OS (HR 1.17; 95%CI 1.04-1.31, p=0.01). Each 1% increment in sTILs was associated with improved IDFS (HR 0.99; 95%CI 0.98-0.99, p<0.001) and OS (HR 0.99, 95%CI 0.98-1.00, p<0.001). Among pts with high sTILs (≥20%), a higher NLR remained significantly associated with worse IDFS (HR 1.21; 95%CI 1.05-1.38, p=0.007) and OS (HR 1.25; 95%CI 1.09-1.44, p=0.001). In contrast, among pts with low sTILs (<20%), NLR was not associated with IDFS (HR 1.07; 95%CI 0.89-1.28, p=0.49) or OS (HR 1.07; 95%CI 0.88-1.30, p=0.49). The interaction test between NLR and sTILs did not reach statistical significance. A multivariate analysis (MVA; including age, menopausal status, histologic subtype, grade, tumor size, nodal stage, Ki-67, NLR, sTILs, adjuvant chemotherapy, type of surgery and adjuvant radiation) showed that sTILs remained independently associated with IDFS (HR 0.99, 95%CI 0.97-1.0, p=0.019) and OS (HR 0.99, 95% CI 0.97-1.0, p=0.044), whereas NLR did not. Conclusions: A lower NLR and a higher sTIL content were each associated with improved IDFS and OS among pts with nonmetastatic TNBC on UVA. However, when evaluated on a MVA, only sTILs remained independently associated with IDFS and OS. Our data suggest that the effect of sTILs on outcomes may not be modified by the NLR. Citation Format: Leon-Ferre RA, Polley M-Y, Liu H, Gilbert J, Cafourek V, Hillman D, Negron V, Boughey JC, Liu MC, Ingle JN, Kalari K, Couch FJ, Visscher DW, Goetz MP. Prognostic value of the neutrophil-to-lymphocyte ratio (NLR) and its relation to stromal tumor infiltrating lymphocytes (sTILs) in triple negative breast cancer (TNBC) [abstract]. In: Proceedings of the 2017 San Antonio Breast Cancer Symposium; 2017 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2018;78(4 Suppl):Abstract nr P3-05-06.

Abstract Purpose: Triple negative breast cancer (TNBC) is an aggressive breast cancer subtype with substantial risks of disease recurrence. While cytotoxic chemotherapy is commonly administered and reduces recurrence, disease outcomes vary considerably and few prognostic tools are available for risk stratification for TNBC patients. We constructed and validated clinical calculators for invasive-disease free survival (IDFS) and overall survival (OS) for TNBC and compared their performance against AJCC-based models which include only tumor size and nodal status. Methods: From a surgical cohort of 9,982 patients who underwent breast cancer surgery at Mayo Clinic between January 1985 and December 2012, 605 centrally reviewed TNBC patients were identified and used to construct Cox models for IDFS and OS. Patients treated with neoadjuvant chemotherapy were excluded. Variables considered included age, menopausal status, tumor size, nodal status, Nottingham grade, type of breast surgery (mastectomy vs. lumpectomy), adjuvant radiation therapy, adjuvant chemotherapy, Ki67, stromal tumor infiltrating lymphocytes (sTILs), and neutrophil-to-lymphocyte ratio (NLR). Missing values were imputed using single imputation with all variables (including outcomes) included in the imputation model. Backward step-down procedure was used for model selections. The final models were internally validated for calibration and discrimination using bootstrapping methods and compared with AJCC-based models. Results: For both IDFS and OS, higher sTIL's, less extensive nodal involvement, use of adjuvant chemotherapy, and lower NLR were significant predictors of improved clinical outcomes. Premenopausal status and younger age were additionally predictive of improved IDFS and OS, respectively. Models for IDFS and OS have good calibration and are associated with bias-corrected C-indices of 0.68 and 0.71, respectively, as compared with C-indices of 0.59 and 0.62 for AJCC-based models. Conclusions: Our data indicate that a clinical calculator that includes sTIL's, NLR, menopausal status, age, nodal involvement as well as chemotherapy use can provide significantly greater prediction of clinical risk than tumor size and nodal status alone. These tools may be used to identify TNBC patients at elevated risk of disease relapse and to aid physician's communication with patients regarding their long-term disease outlook and planning treatment strategies. External validation is required to further evaluate broader applicability of this tool, which was developed utilizing a single-institutional experience. Citation Format: Polley M-YC, Leon-Ferre RA, Liu H, Gilbert J, Cafourek V, Hillman DW, Negron V, Boughey JC, Liu MC, Ingle JN, Kalari K, Couch F, Visscher DW, Goetz MP. Mayo clinic TNBC outcome calculator: A clinical calculator to predict disease relapse and survival in women with triple-negative breast cancer [abstract]. In: Proceedings of the 2017 San Antonio Breast Cancer Symposium; 2017 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2018;78(4 Suppl):Abstract nr P1-06-07.

Abstract Background Women with benign breast disease (BBD) have an increased risk for breast cancer (BC). Almost 30% of all BCs develop in women with prior BBD. Information about the features of the expected BCs after BBD would enable individualized surveillance and prevention strategies for these women. We sought to characterize BCs developing in a large cohort of women with BBD. Methods Our cohort includes 13,485 women who underwent breast biopsy for mammographic or palpable concern between 1967-2001. Biopsy slides were reviewed and classified as nonproliferative (NP), proliferative disease without atypia (PDWA), or atypical hyperplasia (AH). BCs were identified by follow-up questionnaires, medical records, and Tumor Registry data. BC tissues were obtained and reviewed. Results With a median follow-up of 15.8 years, 1273 women developed BC. Most BCs were invasive (81%), among which 61% were ductal, 13% mixed ductal/lobular, and 14% lobular. Two thirds were intermediate or high-grade, and 29% were node positive. Cancer characteristics were similar across the three histologic categories of BBD, with similar frequency of DCIS and invasive disease, tumor size, time to invasive BC, histologic type of BC, nodal positivity and HER2 positivity. Women with AH had a higher frequency of ER+ BC (91%) vs women with PD (80%) or NP (85%), p=0.02. Conclusion A substantial proportion of all BCs develop in women with prior BBD. We show that the majority of BCs after BBD are invasive tumors of ductal type, with a substantial proportion node positive. 84% of cancers were ER-positive. Chemoprevention should be strongly encouraged in higher-risk women with BBD. Citation Format: Visscher DW, Frost MH, Hartmanan LC, Frank RD, Vierkant RA, McCullough AE, Winham SJ, Vachon C, Ghosh K, Brandt KR, Farrell AM, Tarabishy Y, Hieken TJ, Haddad TC, Kraft RA, Radisky DC, Degnim AC. Clinicopathologic features of breast cancers that develop in women with previous benign breast disease. [abstract]. In: Proceedings of the Thirty-Eighth Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2015 Dec 8-12; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2016;76(4 Suppl):Abstract nr P6-10-19.

Abstract Background Women with atypical hyperplasia (AH) are at an approximately four-fold increased risk of subsequent breast cancer (BC). Mammographic breast density (MBD) is a well-established risk factor for BC, but its contribution to BC risk in women with AH remains an open question. We previously reported no association between MBD [measured by Wolfe's parenchymal pattern (PP)] and BC risk in a cohort of 147 women with AH. Here, we present results in an expanded cohort of 459 women diagnosed with AH between 1985 and 2001. Methods The Mayo Clinic Benign Breast Disease Cohort includes 13,485 women who had benign core and/or excisional biopsy 1967-2001. Biopsy tissues were reviewed by our study pathologist to determine presence of AH. MBD was available from clinical records starting in 1985, coded as PP (the standard for 1985-1996) or BI-RADS (1997-2001) density criteria. The original four-level PP (N1-fatty, P1-ductal prominence <25% of breast, P2-ductal prominence >25%, DY-dysplasia) and BI-RADS (fatty, scattered densities, heterogeneously dense, extremely dense) measures were re-categorized as low, moderate or high MBD by combining the middle two categories for each. BC events and clinical information were obtained by questionnaires, medical records and the Mayo Clinic Tumor Registry. Women were followed from benign biopsy to date of BC, death or last contact. Standardized incidence ratios (SIRs) were generated overall and within subgroups defined by density measure (PP vs. BI-RADS), number of atypical foci, and BMI by dividing the observed number of BCs by population-based expected values. Cox regression was used to estimate MBD hazard ratios after adjustment for demographic and clinical variables. Results Of the 551 women diagnosed with AH between 1985 and 2001, 459 (83%) had MBD data within 1 year prior to biopsy. Of these, 68 (15%) had low, 221 (48%) had moderate, and 170 (37%) had high MBD, respectively. Over a median follow-up of 11.7 years, 80 BCs were observed. SIRs for breast cancer did not differ significantly across density categories, overall or within any subgroups examined (see Table). Cox regression adjusting for age, BMI and density measure (PP vs. BI-RADS) also failed to identify an association with MBD (p=0.55). Low MBDModerate MBDHigh MBD N / BCsSIR (95% CI)N / BCsSIR (95% CI)N / BCsSIR (95% CI)P-valueOverall68/123.5 (1.8,6.1)221/393.6 (2.5,4.9)170/293.4 (2.3,4.8)0.97MBD Measure PP59/113.6 (1.8,6.5)85/153.0 (1.7,5.0)130/243.3 (2.2,5.0)0.87BI-RADS9/12.7 (0.1,14.7)136/244.0 (2.6,6.0)40/53.4 (1.1,7.9)0.90No. Atypical Foci 147/62.3 (0.9,5.1)123/182.8 (1.7,4.4)96/163.5 (2.0,5.8)0.63214/47.5 (2.0,19.1)58/134.8 (2.6,8.3)41/62.5 (0.9,5.5)0.213+7/26.7 (0.8,24.0)40/84.4 (1.9,8.6)33/74.0 (1.0,8.2)0.83BMI <2525/21.5 (0.2,5.4)75/184.9 (2.9,7.7)101/173.4 (2.0,5.5)0.1625-2919/54.9 (1.6,11.3)68/92.8 (1.3,5.2)36/42.1 (0.6,5.4)0.4530+23/54.8 (1.6,11.2)76/123.0 (1.6,5.3)32/74.0 (1.6,8.2)0.67SIRs compare observed numer of BCs to expected using Iowa SEER data. Analyses account for the effects of age and calendar period. P-value is test of heterogeneity in SIRs across columns. Conclusions We found no evidence of an association between MBD and subsequent BC in women with AH. Citation Format: Vierkant RA, Degnim AC, Hartmann LC, Frank RD, Radisky DC, Visscher DW, Frost MH, Winham SJ, Ghosh K, Vachon CM. No evidence of association between mammographic breast density and risk of breast cancer in women with atypical hyperplasia. [abstract]. In: Proceedings of the Thirty-Eighth Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2015 Dec 8-12; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2016;76(4 Suppl):Abstract nr P6-09-05.

Abstract Introduction: CD4+ and CD8+ T cells are important effector cells in the adaptive immune response against cancers. In premalignant proliferative breast tissues, the role of the immune system remains relatively undefined. In this study, we sought to investigate and quantify the presence of CD4 and CD8 expressing immune cells in benign breast tissues and their association with breast cancer risk. Method: Archived breast tissue samples from women with benign breast disease (BBD) or no clinical breast disease [Komen Tissue Bank (KTB)] were obtained for this age-matched case-control study. BBD samples included BBD cases (subsequently developed breast cancer) and BBD controls). Tissue sections were characterized for non-proliferative or proliferative disease and degree of lobular involution. CD4+ and CD8+ cell densities (cells/mm2) were obtained by digital image quantitation for up to 10 individual lobules per sample and a per-sample estimate of cell density was calculated as the median cell density across lobules within a sample. Statistical analysis was performed using Wilcoxon signed-rank and Kruskal-Wallis tests. Results: Among 267 women (median age 52, range 35-75) comprising 89 age-matched triplets, 2417 lobules were evaluated [783 KTB, 804 BBD case, 830 BBD control]. 1372 (57%) lobules were normal and 1045 (43%) were fibrocystic. CD4+ cells were present in 1903 (79%) lobules while CD8+ cells were present in 2214 (92%) lobules. CD4 and CD8 densities showed a moderate positive correlation (r = 0.39) with each other, but overall, lobules showed a lower density of CD4+ than CD8+ cells (mean difference -54 cells/mm2, 95% CI: -68 to -39 cells/mm2). At the per sample level (see Table), BBD controls showed significantly higher levels of both CD4+ (p=0.009) and CD8+ cells (p=0.0001) compared to KTB samples. In contrast, the BBD cases only showed elevated CD8 (p=0.007) and not CD4 (p=0.21) cell infiltrates compared to KTB samples. Comparing BBD cases to BBD controls, BBD cases showed significantly lower CD8 cell density as compared to BBD controls (p=0.002) and cases had lower CD4 cell density but the difference was not significant (p=0.07). Despite these differences in the cell densities, the ratio of CD4:CD8+ cells did not differ significantly across groups (p=0.83). The CD4:CD8+ ratio also did not differ significantly by histologic impression (p=0.74), degree of involution (p=0.27), or age (p=0.32). KTBBBD ControlsBBD CasesMedian (IQR)(n=89)(n=89)(n=89)CD4+ cells/mm240 (17-95)60 (25-124)49 (21-92)CD8+ cells/mm285 (49-131)157 (84-225)134 (80-168)Ratio (CD4:CD8)*0.45 (0.14-0.96)0.41 (0.15-0.81)0.40 (0.19-0.70)*Ratio is missing for n=7 samples (3 KTB, 2 BBD control, 2 BBD case) with undefined ratio due to zero CD8+ cells for all lobules in the sample. Conclusion: Our findings suggest that BBD is associated with increased infiltration of T cells into the breast tissue. Furthermore, these immune responses appear more robust in BBD patients at lower risk of developing breast cancer. Citation Format: Pena Jimenez A, Hoskin TL, Arshad MA, Visscher DW, Brahmbhatt RD, Miller EE, Stallings Mann ML, Carter JM, Murphy LM, Winham SJ, Radisky DC, Denison LA, Knutson KA, Degnim AC. CD4 and CD8 T cell densities are increased in benign breast disease compared to normal breast tissue. [abstract]. In: Proceedings of the Thirty-Eighth Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2015 Dec 8-12; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2016;76(4 Suppl):Abstract nr P4-04-09.

Abstract Background: Tumor-infiltrating lymphocytes (TIL) have prognostic and potentially predictive significance in the (neo)adjuvant treatment of high-risk breast cancer. However, quantitative TIL measurement is not routinely performed. It is unclear why some tumors attract large quantities of TIL while others do not. We sought to confirm the association between TIL and pathologic complete response rate (pCR) and to further use next generation sequencing (NGS) to identify genes and gene pathways associated with the presence/absence of TIL. Methods: We studied 140 women with high risk stage I-III breast cancer, enrolled in the Breast Cancer Genome Guided Therapy Study (BEAUTY), obtaining serial biopsies for DNA/RNA sequencing and MRI imaging to assess response to neoadjuvant chemotherapy (NAC) with taxane (+/- trastuzumab+/-pertuzumab for HER2+ disease) followed by AC or (F)EC. Diagnostic pre-NAC core needle biopsies and surgical resection specimens post-NAC were available from 110 patients. Stromal TIL were semi-quantitated on a scale of 1-4 (with 1: ≤10/hpf, 2: subtle infiltrate >10/hpf, 3: moderate infiltrate readily visible at low power magnification, 4: dense infiltrate with innumerable lymphocytes). For this analysis, low TIL was defined as scores of 1-2 vs. high defined as 3-4. Using pre-NAC biopsies, RNAseq was performed using the Illumina HiSeq2000 and the Mayo Analysis Pipeline for RNAseq (MAP-Rseq) for quality control, sequence alignment, and gene counts. The quantity of TIL was associated with transcripts across the transcriptome after conditional quantile normalization. Differentially expressed genes were obtained using EdgeR analysis, using a false discovery rate of 0.05, and pathways were evaluated using GAGE methods. Results: The pCR and residual cancer burden (RCB)-0/I rates by stromal TIL status within each molecular subtype are presented in the table. A diverse spectrum of 1344 genes with differential expression between tumors with high vs. low stromal TIL was identified. The genes with >2.0-fold change (FC) and p<1e-09 included S100A7 (4.49 FC), LCN2 (2.48 FC), and ART3 (2.82 FC) (genes known to be involved in immune regulation), as well as TDRD1 (2.71 FC) (a gene related to ERG [ETS-related gene] expression). In addition, the "regulation of actin cytoskeleton" pathway was upregulated in tumors with high TIL, while the "Hedgehog signaling" and "Wnt signaling" pathways were downregulated. Molecular SubtypeStromal TILspCR rate n (%)RCB-0/I rateLuminal AHigh------Luminal ALow0/9 (0%)0/9 (0%)Luminal BHigh1/9 (11.1%)1/8 (12.5%)Luminal BLow3/24 (12.5%)6/23 (26.1%)ER+/HER2+High3/9 (33.3%)4/9 (44.4%)ER+/HER2+Low1/6 (16.7%)1/6 (16.7%)ER-/HER2+High8/9 (88.9%)7/7 (100%)ER-/HER2+Low4/8 (50.0%)6/8 (75.0%)Triple NegativeHigh10/19 (52.6%)13/19 (68.4%)Triple NegativeLow7/14 (50.0%)9/13 (69.2%) Conclusions: We identified genes and gene pathways associated with high TIL expression in breast tumors prior to NAC that provide insight into the interactions between TIL and tumors. TIL can be easily semi-quantitated on H&E and along with these novel biomarkers, may contribute to the personalization of breast cancer therapy. Citation Format: Moyer AM, Boughey JC, Kalari KR, Suman VJ, McLaughlin SA, Moreno-Aspitia A, Northfelt DW, Gray RJ, Sinnwell JP, Carlson EE, Dockter TJ, Jones KN, Felten SJ, Conners AL, Wieben ED, Ingle JN, Wang L, Weinshilboum RM, Visscher DW, Goetz MP. Differential mRNA expression patterns in breast tumors with high vs. low quantity of stromal tumor–Infiltrating lymphocytes. [abstract]. In: Proceedings of the Thirty-Eighth Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2015 Dec 8-12; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2016;76(4 Suppl):Abstract nr P4-04-05.

Abstract Background: Tumorigenesis involves both genetic and epigenetic changes. Epigenetic alterations are reversible and are promising cancer therapeutic targets. Decitabine (5-aza-2'-deoxycytidine), a DNA methyltransferase inhibitor, is FDA approved for hematological malignancies. However, the effect of decitabine in breast cancer is not completely understood. Previous reports indicated that one decitabine mechanism involves regulation of protein levels for DNMT1, the major DNA methyltransferase that methylates hemimethylated CpG di-nucleotides in DNA. However, the E3 ligase involved in this process has not been identified. Whether decitabine also regulates DNMT3A and 3B in a similar fashion remains unclear. Therefore, our goals were to 1) understand mechanisms underlying decitabine action, 2) test the antitumor activity of decitabine in breast cancer models and 3) identify biomarkers associated with response to decitabine. Methods and Results: Western blots of breast cancer cell lines showed that DNMT1, DNMT3A, and DNMT3B protein levels decreased following decitabine treatment without a reduction in mRNA levels. Bioinformatic analysis of DNA methyltransferase sequences revealed a potential TRAF6 binding motif, and the interaction with TRAF6 (TNF receptor-associated factor 6) was confirmed by IP. TRAF6 functions as an E3 ligase. To determine whether TRAF6 might be the E3 ligase responsible for the degradation of DNMTs after decitabine treatment, we knocked down TRAF6 by RNA interference or knocked out the TRAF6 gene by CRISPR/Cas9. Down regulation of TRAF6 attenuated DNMT ubiquitination and increased DNMT protein levels, suggesting that TRAF6 might mediate proteasome-dependent degradation of all three DNMTs. This was further confirmed by reconstituting the knockout cells with WT and a TRAF6-C70A mutant, followed by assessing DNMT protein levels. Global DNA methylation was also increased after TRAF6 depletion and was confirmed in TRAF6 knock out cells in which DNMT levels were unaffected by decitabine. Cell cytotoxicity and colony forming assays showed that TRAF6 knockout cells were resistant to decitabine, suggesting that a major decitabine mechanism of action is through the regulation of TRAF6 which, in turn, degrades DNMTs, leading to decreased global methylation. Finally, decitabine significantly induced TRAF6 at both mRNA and protein levels, a process that might create positive feedback leading to increased degradation of DNMT proteins upon decitabine treatment. Based on these results, we further hypothesized that levels of the three DNMTs might influence decitabine response. Using 18 breast cancer patient derived xenograft (PDX) models, we found a wide range of DNMT protein levels regardless of ER/HER2 status. DNMT levels in the PDX models were directly associated with sensitivity to decitabine treatment, confirming our hypothesis. Conclusion: Our data showed that decitabine might be an effective agent for treating breast cancer and revealed a novel mechanism underlying decitabine treatment. Baseline DNMT protein levels may serve as a biomarker for predicting decitabine drug response. Citation Format: Yu J, Qin B, Boughey JC, Moyer AM, Visscher DW, Sinnwell JP, Yin P, Thompson KJ, Docter TJ, Kalari KR, Suman VJ, Wieben ED, Felten SJ, Conners AL, Jones KN, McLaughlin SA, Copland JA III, Moreno Aspitia A, Northfelt DW, Gray RJ, Ingle JN, Lou Z, Weinshilboum R, Goetz MP, Wang L. Regulation of DNA methyltransferases via TRAF6 determines breast cancer response to decitabine. [abstract]. In: Proceedings of the Thirty-Eighth Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2015 Dec 8-12; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2016;76(4 Suppl):Abstract nr P3-07-51.