Vierkant, Robert A.

Abstract BackgroundFormalin-fixed, paraffin-embedded (FFPE) tissues have many advantages for identification of risk biomarkers, including wide availability and potential for extended follow-up endpoints. However, RNA derived from archival FFPE samples has limited quality. Here we identified parameters that determine which FFPE samples have the potential for successful RNA extraction, library preparation, and generation of usable RNAseq data.MethodsWe optimized library preparation protocols designed for use with FFPE samples using seven FFPE and Fresh Frozen replicate pairs, and tested optimized protocols using a study set of 130 FFPE biopsies from women with benign breast disease. Metrics from RNA extraction and preparation procedures were collected and compared with bioinformatics sequencing summary statistics. Finally, a decision tree model was built to learn the relationship between pre-sequencing lab metrics and qc pass/fail status as determined by bioinformatics metrics.ResultsSamples that failed bioinformatics qc tended to have low median sample-wise correlation within the cohort (Spearman correlation < 0.75), low number of reads mapped to gene regions (< 25 million), or low number of detectable genes (11,400 # of detected genes with CPM > 2). The median RNA concentration and pre-capture library Qubit values for qc failed samples were 18.7 ng/ul and 1.68 ng/ul respectively, which were significantly lower than those of qc pass samples (40.8 ng/ul and 5.58 ng/ul). We built a decision tree model based on input RNA concentration, input library qubit values, and achieved an F score of 0.803 in predicting QC status (pass/fail) of FFPE samples. ConclusionsWe provide a bioinformatics quality control recommendation for FFPE samples from breast tissue by evaluating bioinformatic and sample metrics. Our results suggest a minimum concentration of 25 ng/ul FFPE-extracted RNA for library preparation and 1.8 ng/ul pre-capture library output to achieve adequate RNA-seq data for downstream bioinformatics analysis.

Abstract Background: Approximately one million U.S women are diagnosed with benign breast disease (BBD) annually, which increases breast cancer (BC) risk between 1.5- and four-fold, depending on the pathologic characteristics of the lesion and other factors. BBD lesions are associated with inflammation and fibrosis, suggesting that inflammatory pathways may play a role in progression of BBD to BC. Aspirin and other non-steroidal anti-inflammatory drugs (NSAIDs) are well tolerated, widely used drugs that block inflammation by inhibiting cyclooxygenase enzymes, which lowers prostaglandin synthesis, and by exerting numerous additional anti-inflammatory effects. NSAIDs are inversely associated with several types of cancer, including BC in some, but not all studies, and have been linked to lower BC risk among BBD patients in prior publications. We examined the association between NSAID use and BC risk, overall and by self-reported history of BBD, using data from the American Cancer Society Cancer Prevention Study II (CPS-II). Methods: The Cancer Prevention Study-II (CPS-II) Nutrition Cohort is a study of approximately 189,000 U.S. men and women followed prospectively for cancer incidence and mortality. Study participants filled out a detailed questionnaire at baseline in 1992/1993 and follow-up surveys in 1997 and every two years until 2017. On each follow-up survey, participants reported type, quantity, and duration of NSAID use as well as diagnoses of BBD. Incident BCs were either self-reported and validated via medical records or cancer registry linkage, or identified via National Death Index linkage. Women in the analytic cohort were followed from date of 1997 survey until diagnosis of BC, reported use of tamoxifen, report of cancer other than BC, date of death, loss to follow-up or end of follow-up. Hazard ratios (HRs) and 95% confidence intervals (CIs) assessing associations of BBD and NSAID use with BC risk were estimated using Cox proportional hazards regression analysis. Statistical interactions between BBD and NSAID were evaluated by fitting and testing the corresponding cross-product terms, along with the main effects of each. Both NSAID use and BBD were modeled as time-dependent variables, allowing exposure status to change over the course of the multiple surveys returned. Results: Of the 60,517 women included in the study, 3,899 (6.4%) developed BC. Mean follow-up time was 12.7 years. A total of 23,661 women (39%) reported a diagnosis of BBD at some point during the study period, and 54,656 (90%) reported some NSAID use. Women self-reporting BBD were 46% more likely to develop BC than those without self-reported BBD (HR 1.46, 95% CI 1.36-1.57). Overall, BC risk was not increased among NSAID users versus non-users (HR 0.95, 95% CI 0.88-1.02). In contrast, BC risk was reduced among current NSAID users with a history of BBD (HR 0.87, 95% CI 0.78-0.97), but not among NSAID users who did not report a history of BBD (HR 1.02, 95% CI 0.92-1.13, p-interaction=0.04). Risk did not differ across number of pills per month (RR=0.84 for &lt;15 pills/month, 0.89 for 15-29, 0.88 for 30-44, 0.78 for 45-59, 0.88 for 60+). Conclusion: In our study, NSAID use had a statistically-significant protective effect on BC risk among women with BBD, suggesting its potential use as a chemopreventive agent in this high-risk group of women, but not among those who did not report a history of BBD. Citation Format: Robert A Vierkant, Matthew Masters, Lauren R Teras, Mark E Sherman. Breast cancer (BC) risk among patients with benign breast disease (BBD) by NSAID use [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr P2-11-01.

Abstract Background:Breast terminal duct lobular units (TDLUs), the source of most breast cancer (BC) precursors, are shaped by age-related involution, a gradual process, and postpartum involution (PPI), a dramatic inflammatory process that restores baseline microanatomy after weaning. Dysregulated PPI is implicated in the pathogenesis of postpartum BCs. We propose that assessment of TDLUs in the postpartum period may have value in risk estimation, but characteristics of these tissues in relation to epidemiological factors are incompletely described.Methods: Using validated Artificial Intelligence and morphometric methods, we analyzed digitized images of tissue sections of normal breast tissues stained with hematoxylin and eosin from donors <45 years from the Komen Tissue Bank (180 parous and 545 nulliparous). Metrics assessed by AI, included: TDLU count; adipose tissue fraction; mean acini count/TDLU; mean dilated acini; mean average acini area; mean “capillary” area; mean epithelial area; mean ratio of epithelial area versus intralobular stroma; mean mononuclear cell count (surrogate of immune cells); mean fat area proximate to TDLUs and TDLU area. We compared epidemiologic characteristics collected via questionnaire by parity status and race, using a Wilcoxon rank sum test or Fisher’s exact test. Histologic features were compared between nulliparous and parous women (overall and by time between last birth and donation [ recent birth: ≤5 years versus remote birth: >5 years]) using multivariable regression models.Results:Normal breast tissues of parous women contained significantly higher TDLU counts and acini counts, more frequent dilated acini, higher mononuclear cell counts in TDLUs and smaller acini area per TDLU than nulliparas (all multivariable analyses p<0.001). Differences in TDLU counts and average acini size persisted for >5 years postpartum whereas increases in immune cells were most marked <5 years of a birth. Relationships were suggestively modified by several other factors, including demographic and reproductive characteristics, ethanol consumption and breastfeeding duration. Conclusions:Our study identified sustained expansion of TDLU numbers and reduced average acini area among parous versus nulliparous women and notable increases in immune responses within five years following childbirth. Further, we show that quantitative characteristics of normal breast samples vary with demographic features and BC risk factors.

NT-proBNP, soluble ST2 (sST2), and galectin-3 are biomarkers of cardiac dysfunction that have been proposed as identifiers of patients experiencing asymptomatic cardiac dysfunction after anthracycline-based chemotherapy. This study aimed to compare the proportion of breast cancer (BC) survivors with elevated serum levels of these three putative biomarkers by prior receipt of anthracycline (yes vs. no). Five-hundred-eighty survivors of BC who had received anthracycline-based chemotherapy were matched by age and time between diagnosis and serum storage to 580 who had not. Cardiac biomarker levels were analyzed using immunoassays. Analyses were carried out using linear and logistic regression models. Anthracycline recipients had higher values of NT-proBNP than non-recipients (mean 116.0 ng/L vs. 97.0 ng/L, respectively; p < 0.001). Values for ST2 and galectin-3 did not significantly differ by receipt of anthracycline. After further adjustment for age at breast cancer diagnosis, ethnicity, and receipt of trastuzumab, associations between receipt of anthracycline and higher NT-proBNP persisted (p < 0.001), showing that NT-proBNP may be a biomarker of cardiovascular toxicity after receipt of anthracycline-based chemotherapy. Further research to assess the clinical utility of NT-proBNP testing after receipt of anthracycline is recommended. sST2 and galectin-3 do not appear to differentiate between anthracycline recipients and non-recipients amongst breast cancer survivors.

Abstract Benign breast disease (BBD) is an established risk factor for development of breast cancer (BC). However, little is known about the prevalence of somatic genetic alterations at the time of BBD diagnosis and relationship to subsequent BC. To characterize DNA variants in BBD and investigate association with BC risk, we isolated DNA from formalin-fixed paraffin-embedded (FFPE) BBD biopsies and assessed mutations using a panel of 93 breast cancer predisposition genes (Qiagen GeneRead). The study set consisted of 120 BBD patients from three age-matched groups defined based on 16 years of BC follow-up: BBD-controls (cancer-free after 16 years), BBD-ER+, or BBD-ER- breast cancer (within 16 years). Bioinformatics quality checks included comparisons of paired fresh-frozen and FFPE replicates to reduce false variant calls induced by FFPE artifacts. Gene-level variant burden differences among study groups were assessed using 12 statistical approaches, adjusted for patient age, year of biopsy, and histology. To investigate possible immune implications of increased variant burden, extent of leukocytes in the BBD biopsies was assessed using CD45 immunohistochemistry. To assess variant profile changes during the progression from BBD to BC, we also performed targeted sequencing of 12 BBD-matched subsequent tumor samples (six ER+ and six ER-) developed within 5 years of BBD diagnosis. Gene-level association results showed consistently skewed distributions suggesting a higher variant burden in BBD controls: among the top 10 genes associated with case-status (p&lt;0.05), all had BC-risk odds ratio estimates less than 1.0. Immune infiltrates, assessed by CD45 H-scores, were positively associated with overall variant burden (r=0.42, p=0.0031), and noticeably higher in cancer-free controls (r=0.50, p=0.005). Comparison of mutation variants in 12 BBD cases that developed cancer within 5 years of diagnosis and their subsequent tumors revealed that more than half of the non-SNP variants were shared between the BBD and the subsequent tumor (54.6%), whereas less than one quarter were unique to the tumor (20.6%) or unique to the BBD (24.9%). Through BBD-BC paired variant analysis, three subjects were found with newly emergent clones uniquely appearing in the tumors and four subjects had loss of clones in the tumor that were present in the BBD. These results show that somatic variant burden is lower in patients with BBD that progress to BC, mutation profiles of BC progression are heterogeneous, and BBD tissues with reduced variant burden also show reduced immune infiltrate. This implicates a potential protective role for immune activation in BBD. Future directions will define the characteristics of the potentially protective immune response and its role in variant heterogeneity in a larger sample set, integrated with paired RNA sequencing to better understand the pathways involved and the relationship to immune infiltration. Citation Format: Stacey J. Winham, Julie M. Cunningham, Yuanhang Liu, Aditya V. Bhagwate, Ethan P. Heinzen, Samantha J. McDonough, Melody L. Stallings-Mann, Robert A. Vierkant, Jodi M. Carter, Mark E. Sherman, Derek C. Radisky, Amy C. Degnim, Chen Wang. Somatic variant burdenin benign breast disease and association with risk of subsequent breast cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2274.

Abstract Background: Benign breast disease (BBD) is a risk factor for breast cancer (BC); however, little is known about the genetic alterations present at the time of BBD diagnosis, and how these relate to risk of incident BC. Methods: A subset of a long-term BBD cohort was selected to examine DNA variation across three BBD groups (42 future estrogen receptor-positive (ER+) BC, 36 future estrogen receptor-negative (ER-) BC and 42 controls cancer-free for at least 16 years post-BBD). DNA extracted from archival formalin fixed, paraffin-embedded (FFPE) tissue blocks was tested and analyzed for presence of DNA alterations using a targeted panel of 93 BC-associated genes. To address artifacts frequently observed in FFPE tissues (e.g., C>T changes), we applied three filtering strategies based on alternative allele frequencies and nucleotide substitution context. Gene-level associations were performed using two burden tests, SKAT-O and logistic regression, and adjusted for clinical and technical covariates including epithelial percentage, histologic impression, patient age and year of BBD biopsy. Results: After filtering, the variant frequency of SNPs in our samples was highly consistent with population allele frequencies reported in 1KG/ExAC (0.986, p <1e-16). Ten genes were found to be associated with later cancer status by four of 12 association methods, with nominal p-values <0.05: MED12, MSH2, BRIP1, PMS1, GATA3, MUC16, FAM175A, EXT2, MLH1 and TGFB1. Additionally, analysis revealed 11 and 19 gene-level associations with p-values <0.05 and OR<1 for all cancer cases and ER- cases, respectively. Comparing between the three case groups, BBD ER+ cases were closer to controls in mutation profile, while BBD ER- cases were distinct. Notably, the variant burden was significantly higher in controls than in either ER+ or ER- cases. CD45 expression was associated with mutational burden (p=0.003).Conclusions: Somatic mutations were more frequent in benign breast tissue from women who did not develop cancer, opening questions of clonal diversity or immune-mediated restraint on future cancer development. CD45 expression was positively associated with mutational burden, most strongly in controls. Further studies in both normal and premalignant tissues are needed to better understand the role of somatic gene mutations and their contribution to future cancer development.

Abstract Background: Women with benign breast disease (BBD) are at increased risk of breast cancer (BC). Classic studies based on film-based mammographic screening and pathology diagnosis of surgical biopsies conducted in the 1980s established a hierarchy of increasing BC risk: non-proliferative (NP) BBD, proliferative BBD without atypia (PDWA) and atypical hyperplasia (AH). Given changes in epidemiological BC risk factors and introduction of percutaneous core needle biopsy (CNB) in mid-1990s, and later, digital mammography, we hypothesized that the patient characteristics and relative frequency of BBD diagnoses have changed over time. Accordingly, we performed a longitudinal analysis of the frequency of patient characteristics and BBD diagnoses in the Mayo BBD cohort. Methods: Utilizing the Mayo Clinic Surgical and Pathology Indices, women ages 18 to 85 who had a BBD biopsy between 1/1/67 and 12/31/13 were identified. Breast pathologists reviewed biopsies masked to diagnoses of incident BC diagnosed in follow-up. Demographic characteristics and BC events were obtained by query of institutional data sources and participant surveys. Trends were evaluated for the following eras: 1: pre-mammogram (1967-1981), 2: pre-CNB (1982-1992), 3: CNB Transition (1993-2001), and 4: CNB (2002-2013). Demographics were formally compared across eras using chi-square tests for categorical variables and analyses of variance (ANOVAs) for continuous variables. Results: From 1967-2013, the cohort includes 19,582 unique women with BBD. The frequency of CNB increased from eras 1-4: 0.04%, 0.6%, 51.3 %, and 88.9%, respectively. Mean age at BBD diagnosis was younger in era 1 (48.0 years) vs eras 2-4 (53.2, 52.0, and 51.8, respectively, p&lt;0.001). The percentage of biopsies diagnosed as PDWA increased from era 1-4 (25.7%, 34.3%, 35.2%, 46.2%, p&lt;0.001), as did the percentage with AH (2.4%, 5.1%, 8.6%, 12.3%, p&lt;0.001). Over eras 1-4, the percentages of women with a strong family history of BC increased (9.9%, 12.7%, 17.1%, and 29.0%, p&lt;0.001) as did mean BMI (24.8, 26.4, 27.4, and 28.6, p&lt;0.001). With a median follow-up of 10.9 years, 1,719 breast cancers have developed, with increasing proportion of noninvasive (DCIS-only) disease across eras 1-4: 15.0%, 21.1%, 21.2%, and 33.2%, p&lt; 0.001. Conclusions: Analysis of this large, single institution BBD cohort for the 46 year period 1967-2013 demonstrates that BC risk factors among BBD patients has changed over time, with subjects demonstrating increasing age, BMI, and family history, and that the percentages of BBD classified as PDWA and AH have increased. Impact on BC risk will be further investigated. Citation Format: Amy Degnim, Karthik Ghosh, Jodi M Carter, Robert A Vierkant, Matthew R Jensen, Stacey J Winham, Tanya L Hoskin, Marlene Frost, Teresa M Allers, Denise L Gehling, Mindy J Kern, Laura M Pacheco-Spann, Celine M Vachon, Derek C Radisky, Daniel W Visscher, Mark E Sherman. Benign breast disease: Temporal trends from 1967 to 2013 [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr PS7-11.

Abstract Background: Women with benign breast disease (BBD) experience an increased risk of breast cancer (BC). Histologic classification of BBD, as non-proliferative disease (NP), proliferative disease without atypia (PDWA) or atypical hyperplasia (AH), stratifies groups of patients into progressively higher categories of BC risk. However, this classification does not comprehensively assess the proliferative state of the epithelium throughout the biopsy. In addition, while AH is considered the most high-risk class of BBD, it is not always a highly proliferative lesion; atypical ductal lesions may reflect focal cytologic and architectural changes. We evaluated the association of an alternative classification of BBD severity and BC risk based on subjective grading of: 1) the maximal degree of epithelial proliferation and 2) multifocality of epithelial proliferation. Methods: Pathologists reviewed biopsies from participants aged 18 to 85 years in the Mayo BBD cohort (2002-2013), masked to BC outcomes, ascertained via questionnaires, tumor registry data and medical record review. Biopsies were classified as NP, PDWA or AH and semi-quantitatively scored for: 1) maximal degree of epithelial proliferation within a focus (DP) (0-3; none to severe) and 2) multifocality of proliferation (MP) (0-3; none to multiple foci). DP and MP scores were also summed to give a DP+MP score (0-6). Associations of DP and MP with BC risk were examined using Cox proportional hazards regression analyses, adjusting for age at BBD biopsy. Women were followed from date of initial biopsy to date of BC, death or last follow-up. Results: Of the 1529 assessable biopsies, 544 (35.6%) were classified as NP, 708 (46.3%) as PDWA and 277 (18.1%) as AH. Both DP and MP scores had significant positive correlation with increasing BBD severity (DP: r=0.51, p&lt; 0.001; MP: r=0.52, p&lt; 0.001). Mean (SD) DP scores were 0.6 (0.6) for NP, 1.6 (0.9) for PDWA, and 1.8 (0.7) for AH (ANOVA p&lt;0.001). Mean (SD) for MP scores were 0.6 (0.6) for NP, 1.4 (0.8) for PDWA, and 1.8 (0.8) for AH (ANOVA p&lt;0.001). Mean (SD) for DP+MP scores were 1.2 (1.2) for NP, 2.9 (1.5) for PDWA, and 3.6 (1.2) for AH (ANOVA p&lt;0.001). With median follow-up of 8.8 years for controls and 5.3 years for cases, 10.6% of the women in the cohort developed BC. Compared to those with DP scores of 0, women with DP scores of 3 had significantly increased BC risk (HR 1.42, 95% CI: 1.16, 1.74, p=0.003). MP was associated with a non-significant increase in BC risk for scores of 3 versus 0 (HR: 1.20, 95% CI: 0.97,1.49, p=0.11). DP+MP scores of 6 conferred the highest BC risk (HR (score 6 vs. 0): 1.62, 95% CI 1.18,2.21, p=0.02). Results did not substantively differ after adjusting for BBD severity as NP, PDWA or AH. Conclusions: In this preliminary analysis within the Mayo BBD cohort, both proliferative degree (DP) and multifocality (MP) scores were correlated with histologic severity of BBD. DP and DP+MP scores were each associated with increased BC risk. We conclude that improved characterization of epithelial proliferation in BBD biopsies may enable refined prediction of individual BC risk. Citation Format: Jodi M Carter, Matthew R Jensen, Robert A Vierkant, Stacey J Winham, Tanya L Hoskin, Marlene Frost, Karthik Ghosh, Derek C Radisky, Amy C Degnim, Mark E Sherman. Epithelial proliferation score as an independent breast cancer risk predictor in benign breast disease [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr PD10-10.

AbstractTumor‐infiltrating lymphocytes in colorectal cancer (CRC) predict better survival. However, associations between T‐lymphocyte count in histologically normal tissues from patients with CRC and survival remain uncertain. We examined associations of CD3+ T‐cells in colorectal tumor and histologically normal tissues with CRC‐specific and all‐cause mortality in the prospective Iowa Women's Health Study. Tissue microarrays were constructed using paraffin‐embedded colorectal tissue samples from 464 women with tumor tissues and 314 women with histologically normal tissues (55–69 years at baseline) diagnosed with incident CRC from 1986 to 2002 and followed through 2014 (median follow‐up 20.5 years). Three tumor and two histologically normal tissue cores for each patient were immunostained using CD3+ antibody and quantified, and the counts were averaged across the cores in each tissue. Cox proportional hazards regression estimated hazard ratios (HR) and 95% confidence interval (CI) for CRC‐specific and all‐cause mortality. After adjustment for age at diagnosis, body mass index, smoking status, tumor grade, and stage, HRs (95% CI) for the highest versus lowest tertile of tumor CD3+ score were 0.59 (0.38–0.89) for CRC‐specific mortality and 0.82 (0.63–1.05) for all‐cause mortality; for histologically normal CD3+ score, the corresponding HRs (95% CI) were 0.47 (0.19–1.17) and 0.50 (0.27–0.90), respectively. The CD3+ score combining the tumor and histologically normal scores was inversely associated with CRC‐specific and all‐cause mortality. Although the association between tumor CD3+ score and all‐cause mortality was not significant, both higher CD3+ T‐lymphocyte counts in tumor and histologically normal scores tended to be associated with lower CRC‐specific and all‐cause mortality.

Abstract Background: Accumulating evidence suggests a relationship between endometrial cancer and ovarian cancer. Independent genome-wide association studies (GWAS) for endometrial cancer and ovarian cancer have identified 16 and 27 risk regions, respectively, four of which overlap between the two cancers. We aimed to identify joint endometrial and ovarian cancer risk loci by performing a meta-analysis of GWAS summary statistics from these two cancers. Methods: Using LDScore regression, we explored the genetic correlation between endometrial cancer and ovarian cancer. To identify loci associated with the risk of both cancers, we implemented a pipeline of statistical genetic analyses (i.e., inverse-variance meta-analysis, colocalization, and M-values) and performed analyses stratified by subtype. Candidate target genes were then prioritized using functional genomic data. Results: Genetic correlation analysis revealed significant genetic correlation between the two cancers (rG = 0.43, P = 2.66 × 10−5). We found seven loci associated with risk for both cancers (PBonferroni &lt; 2.4 × 10−9). In addition, four novel subgenome-wide regions at 7p22.2, 7q22.1, 9p12, and 11q13.3 were identified (P &lt; 5 × 10−7). Promoter-associated HiChIP chromatin loops from immortalized endometrium and ovarian cell lines and expression quantitative trait loci data highlighted candidate target genes for further investigation. Conclusions: Using cross-cancer GWAS meta-analysis, we have identified several joint endometrial and ovarian cancer risk loci and candidate target genes for future functional analysis. Impact: Our research highlights the shared genetic relationship between endometrial cancer and ovarian cancer. Further studies in larger sample sets are required to confirm our findings.