A proteome-wide association study identifies putative causal proteins for breast cancer risk.

BACKGROUND: Genome-wide association studies (GWAS) have identified more than 200 breast cancer risk-associated genetic loci, yet the causal genes and biological mechanisms for most loci remain elusive. Proteins, as final gene products, are pivotal in cellular function. In this study, we conducted a proteome-wide association study (PWAS) to identify proteins in breast tissue related to breast cancer risk. METHODS: We profiled the proteome in fresh frozen breast tissue samples from 120 cancer-free European-ancestry women from the Susan G. Komen Tissue Bank (KTB). Protein expression levels were log2-transformed then normalized via quantile and inverse-rank transformations. GWAS data were also generated for these 120 samples. These data were used to build statistical models to predict protein expression levels via cis-genetic variants using the elastic net method. The prediction models were then applied to the GWAS summary statistics data of 133,384 breast cancer cases and 113,789 controls to assess the associations of genetically predicted protein expression levels with breast cancer risk overall and its subtypes using the S-PrediXcan method. RESULTS: A total of 6388 proteins were detected in the normal breast tissue samples from 120 women with a high detection false discovery rate (FDR) p value < 0.01. Among the 5820 proteins detected in more than 80% of participants, prediction models were successfully built for 2060 proteins with R > 0.1 and P < 0.05. Among these 2060 proteins, five proteins were significantly associated with overall breast cancer risk at an FDR p value < 0.1. Among these five proteins, the corresponding genes for proteins COPG1, DCTN3, and DDX6 were located at least 1 Megabase away from the GWAS-identified breast cancer risk variants. COPG1 was associated with an increased risk of breast cancer with a p value of 8.54 × 10-4. Both DCTN3 and DDX6 were associated with a decreased risk of breast cancer with p values of 1.01 × 10-3 and 3.25 × 10-4, respectively. The corresponding genes for the remaining two proteins, LSP1 and DNAJA3, were located in previously GWAS-identified breast cancer risk loci. After adjusting for GWAS-identified risk variants, the association for DNAJA3 was still significant (p value of 9.15 × 10-5 and adjusted p value of 1.94 × 10-4). However, the significance for LSP1 became weaker with a p value of 0.62. Stratification analyses by breast cancer subtypes identified three proteins, SMARCC1, LSP1, and NCKAP1L, associated with luminal A, luminal B, and ER-positive breast cancer. NCKAP1L was located at least 1Mb away from the GWAS-identified breast cancer risk variants. After adjusting for GWAS-identified breast cancer risk variants, the association for protein LSP1 was still significant (adjusted p value of 6.43 × 10-3 for luminal B subtype). CONCLUSION: We conducted the first breast-tissue-based PWAS and identified seven proteins associated with breast cancer, including five proteins not previously implicated. These findings help improve our understanding of the underlying genetic mechanism of breast cancer development.

Host, reproductive, and lifestyle factors in relation to quantitative histologic metrics of the normal breast.

BACKGROUND: Emerging data indicate that variations in quantitative epithelial and stromal tissue composition and their relative abundance in benign breast biopsies independently impact risk of future invasive breast cancer. To gain further insights into breast cancer etiopathogenesis, we investigated associations between epidemiological factors and quantitative tissue composition metrics of the normal breast. METHODS: The study participants were 4108 healthy women ages 18-75 years who voluntarily donated breast tissue to the US-based Susan G. Komen Tissue Bank (KTB; 2008-2019). Using high-accuracy machine learning algorithms, we quantified the percentage of epithelial, stromal, adipose, and fibroglandular tissue, as well as the proportion of fibroglandular tissue that is epithelium relative to stroma (i.e., epithelium-to-stroma proportion, ESP) on digitized hematoxylin and eosin (H&E)-stained normal breast biopsy specimens. Data on epidemiological factors were obtained from participants using a detailed questionnaire administered at the time of tissue donation. Associations between epidemiological factors and square root transformed tissue metrics were investigated using multivariable linear regression models. RESULTS: With increasing age, the amount of stromal, epithelial, and fibroglandular tissue declined and adipose tissue increased, while that of ESP demonstrated a bimodal pattern. Several epidemiological factors were associated with individual tissue composition metrics, impacting ESP as a result. Compared with premenopausal women, postmenopausal women had lower ESP [ß (95% Confidence Interval (CI)) = -0.28 (- 0.43, - 0.13); P < 0.001] with ESP peaks at 30-40 years and 60-70 years among pre- and postmenopausal women, respectively. Pregnancy [ß (95%CI) vs nulligravid = 0.19 (0.08, 0.30); P < 0.001] and increasing number of live births (P-trend < 0.001) were positively associated with ESP, while breastfeeding was inversely associated with ESP [ß (95%CI) vs no breastfeeding = -0.15 (- 0.29, - 0.01); P = 0.036]. A positive family history of breast cancer (FHBC) [ß (95%CI) vs no FHBC = 0.14 (0.02-0.26); P = 0.02], being overweight or obese [ß (95%CI) vs normal weight = 0.18 (0.06-0.30); P = 0.004 and 0.32 (0.21-0.44); P < 0.001, respectively], and Black race [ß (95%CI) vs White = 0.12 (- 0.005, 0.25); P = 0.06] were positively associated with ESP. CONCLUSION: Our findings revealed that cumulative exposure to etiological factors over the lifespan impacts normal breast tissue composition metrics, individually or jointly, to alter their dynamic equilibrium, with potential implications for breast cancer susceptibility and tumor etiologic heterogeneity.

Towards defining morphologic parameters of normal parous and nulliparous breast tissues by artificial intelligence.

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.

Serum hormone levels and normal breast histology among premenopausal women.

PURPOSE: Breast terminal duct lobular units (TDLUs) are the main source of breast cancer (BC) precursors. Higher serum concentrations of hormones and growth factors have been linked to increased TDLU numbers and to elevated BC risk, with variable effects by menopausal status. We assessed associations of circulating factors with breast histology among premenopausal women using artificial intelligence (AI) and preliminarily tested whether parity modifies associations. METHODS: Pathology AI analysis was performed on 316 digital images of H&E-stained sections of normal breast tissues from Komen Tissue Bank donors ages ≤ 45 years to assess 11 quantitative metrics. Associations of circulating factors with AI metrics were assessed using regression analyses, with inclusion of interaction terms to assess effect modification. RESULTS: Higher prolactin levels were related to larger TDLU area (p < 0.001) and increased presence of adipose tissue proximate to TDLUs (p < 0.001), with less significant positive associations for acini counts (p = 0.012), dilated acini (p = 0.043), capillary area (p = 0.014), epithelial area (p = 0.007), and mononuclear cell counts (p = 0.017). Testosterone levels were associated with increased TDLU counts (p < 0.001), irrespective of parity, but associations differed by adipose tissue content. AI data for TDLU counts generally agreed with prior visual assessments. CONCLUSION: Among premenopausal women, serum hormone levels linked to BC risk were also associated with quantitative features of normal breast tissue. These relationships were suggestively modified by parity status and tissue composition. We conclude that the microanatomic features of normal breast tissue may represent a marker of BC risk.

A risk-associated Active transcriptome phenotype expressed by histologically normal human breast tissue and linked to a pro-tumorigenic adipocyte population.

BACKGROUND: Previous studies have identified and validated a risk-associated Active transcriptome phenotype commonly expressed in the cancer-adjacent and histologically normal epithelium, stroma, and adipose containing peritumor microenvironment of clinically established invasive breast cancers, conferring a 2.5- to 3-fold later risk of dying from recurrent breast cancer. Expression of this Active transcriptome phenotype has not yet been evaluated in normal breast tissue samples unassociated with any benign or malignant lesions; however, it has been associated with increased peritumor adipocyte composition. METHODS: Detailed histologic and transcriptomic (RNAseq) analyses were performed on normal breast biopsy samples from 151 healthy, parous, non-obese (mean BMI = 29.60 ± 7.92) women, ages 27-66 who donated core breast biopsy samples to the Komen Tissue Bank, and whose average breast cancer risk estimate (Gail score) at the time of biopsy (1.27 ± 1.34) would not qualify them for endocrine prevention therapy. RESULTS: Full genome RNA sequencing (RNAseq) identified 52% (78/151) of these normal breast samples as expressing the Active breast phenotype. While Active signature genes were found to be most variably expressed in mammary adipocytes, donors with the Active phenotype had no difference in BMI but significantly higher Gail scores (1.46 vs. 1.18; p = 0.007). Active breast samples possessed 1.6-fold more (~ 80%) adipocyte nuclei, larger cross-sectional adipocyte areas (p < 0.01), and 0.5-fold fewer stromal and epithelial cell nuclei (p < 1e-6). Infrequent low-level expression of cancer gene hotspot mutations was detected but not enriched in the Active breast samples. Active samples were enriched in gene sets associated with adipogenesis and fat metabolism (FDR q ≤ 10%), higher signature scores for cAMP-dependent lipolysis known to drive breast cancer progression, white adipose tissue browning (Wilcoxon p < 0.01), and genes associated with adipocyte activation (leptin, adiponectin) and remodeling (CAV1, BNIP3), adipokine growth factors (IGF-1, FGF2), and pro-inflammatory fat signaling (IKBKG, CCL13). CONCLUSIONS: The risk-associated Active transcriptome phenotype first identified in cancer-adjacent breast tissues also occurs commonly in healthy women without breast disease who do not qualify for breast cancer chemoprevention, and independently of breast expressed cancer-associated mutations. The risk-associated Active phenotype appears driven by a pro-tumorigenic adipocyte microenvironment that can predate breast cancer development.

DNA methylation age is elevated in breast tissue of healthy women.

BACKGROUND: Limited evidence suggests that female breast tissue ages faster than other parts of the body according to an epigenetic biomarker of aging known as the "epigenetic clock." However, it is unknown whether breast tissue samples from healthy women show a similar accelerated aging effect relative to other tissues, and what could drive this acceleration. The goal of this study is to validate our initial finding of advanced DNA methylation (DNAm) age in breast tissue, by directly comparing it to that of peripheral blood tissue from the same individuals, and to do a preliminary assessment of hormonal factors that could explain the difference. METHODS: We utilized n = 80 breast and 80 matching blood tissue samples collected from 40 healthy female participants of the Susan G. Komen Tissue Bank at the Indiana University Simon Cancer Center who donated these samples at two time points spaced at least a year apart. DNA methylation levels (Illumina 450K platform) were used to estimate the DNAm age. RESULTS: DNAm age was highly correlated with chronological age in both peripheral blood (r = 0.94, p < 0.0001) and breast tissues (r = 0.86, p < 0.0001). A measure of epigenetic age acceleration (age-adjusted DNAm Age) was substantially increased in breast relative to peripheral blood tissue (p = 1.6 × 10-11). The difference between DNAm age of breast and blood decreased with advancing chronologic age (r = -0.53, p = 4.4 × 10-4). CONCLUSIONS: Our data clearly demonstrate that female breast tissue has a higher epigenetic age than blood collected from the same subject. We also observe that the degree of elevation in breast diminishes with advancing age. Future larger studies will be needed to examine associations between epigenetic age acceleration and cumulative hormone exposure.

Characterizing the heterogeneity of triple-negative breast cancers using microdissected normal ductal epithelium and RNA-sequencing.

Triple-negative breast cancers (TNBCs) are a heterogeneous set of tumors defined by an absence of actionable therapeutic targets (ER, PR, and HER-2). Microdissected normal ductal epithelium from healthy volunteers represents a novel comparator to reveal insights into TNBC heterogeneity and to inform drug development. Using RNA-sequencing data from our institution and The Cancer Genome Atlas (TCGA) we compared the transcriptomes of 94 TNBCs, 20 microdissected normal breast tissues from healthy volunteers from the Susan G. Komen for the Cure Tissue Bank, and 10 histologically normal tissues adjacent to tumor. Pathway analysis comparing TNBCs to optimized normal controls of microdissected normal epithelium versus classic controls composed of adjacent normal tissue revealed distinct molecular signatures. Differential gene expression of TNBC compared with normal comparators demonstrated important findings for TNBC-specific clinical trials testing targeted agents; lack of over-expression for negative studies and over-expression in studies with drug activity. Next, by comparing each individual TNBC to the set of microdissected normals, we demonstrate that TNBC heterogeneity is attributable to transcriptional chaos, is associated with non-silent DNA mutational load, and explains transcriptional heterogeneity in addition to known molecular subtypes. Finally, chaos analysis identified 146 core genes dysregulated in >90 % of TNBCs revealing an over-expressed central network. In conclusion, use of microdissected normal ductal epithelium from healthy volunteers enables an optimized approach for studying TNBC and uncovers biological heterogeneity mediated by transcriptional chaos.

Outdoor air pollution and histologic composition of normal breast tissue.

BACKGROUND: Biologic pathways underlying the association between outdoor air pollution and breast cancer risk are poorly understood. Breast tissue composition may reflect cumulative exposure to breast cancer risk factors and has been associated with breast cancer risk among patients with benign breast disease. Herein, we evaluated whether fine particulate matter (PM2.5) was associated with the histologic composition of normal breast tissue. METHODS: Machine-learning algorithms were applied to digitized hematoxylin and eosin-stained biopsies of normal breast tissue to quantify the epithelium, stroma, adipose and total tissue area from 3,977 individuals aged 18-75 years from a primarily Midwestern United States population who donated breast tissue samples to the Susan G. Komen Tissue Bank (2009-2019). Annual levels of PM2.5 were assigned to each woman's residential address based on year of tissue donation. We applied predictive k-means to assign participants to clusters with similar PM2.5 chemical composition and used linear regression to examine the cross-sectional associations between a 5-µg/m3 increase in PM2.5 and square root-transformed proportions of epithelium, stroma, adipose, and epithelium-to-stroma proportion [ESP], overall and by PM2.5 cluster. RESULTS: Higher residential PM2.5 was associated with lower proportion of breast stromal tissue [ß = -0.93, 95% confidence interval: (-1.52, -0.33)], but was not related to the proportion of epithelium [ß = -0.11 (-0.34, 0.11)]. Although PM2.5 was not associated with ESP overall [ß = 0.24 (-0.16, 0.64)], the association significantly differed by PM2.5 chemical composition (p-interaction = 0.04), with a positive association evident only among an urban, Midwestern cluster with higher concentrations of nitrate (NO3-) and ammonium (NH4+) [ß = 0.49 (0.03, 0.95)]. CONCLUSIONS: Our findings are consistent with a possible role of PM2.5 in breast cancer etiology and suggest that changes in breast tissue composition may be a potential pathway by which outdoor air pollution impacts breast cancer risk. This study further underscores the importance of considering heterogeneity in PM2.5 composition and its impact on breast carcinogenesis.

The impact of reproductive factors on DNA methylation-based telomere length in healthy breast tissue.

Estrogen promotes breast tissue proliferation and telomerase activation. We investigated the effects of reproductive history on cell cycling and telomere length using a DNA methylation-based estimate of telomere length (DNAmTL) in breast and blood from healthy women donors. We demonstrate that DNAmTL is shorter in breast than in blood, and that nulliparous women have longer age-adjusted DNAmTL in both breast and blood, potentially explaining their higher risk of breast cancer.

The Effects of Lifetime Estrogen Exposure on Breast Epigenetic Age.

BACKGROUND: Estrogens are thought to contribute to breast cancer risk through cell cycling and accelerated breast aging. We hypothesize that lifetime estrogen exposure drives early epigenetic breast aging observed in healthy women. In this study, we examined associations between hormonal factors and epigenetic aging measures in healthy breast tissues. METHODS: We extracted DNA from breast tissue specimens from 192 healthy female donors to the Susan G. Komen Tissue Bank at the Indiana University Simon Cancer Center. Methylation experiments were performed using the Illumina EPIC 850K array platform. Age-adjusted regression models were used to examine for associations between factors related to estrogen exposure and five DNA methylation-based estimates: Grim age, pan-tissue age, Hannum age, phenotypic age, and skin and blood clock age. RESULTS: Women were aged 19-90 years, with 95 premenopausal, and 97 nulliparous women. The age difference (Grim age - chronologic age) was higher at earlier ages close to menarche. We found significant associations between earlier age at menarche and age-adjusted accelerations according to the Grim clock, the skin and blood clock, and between higher body mass index (BMI) and age-adjusted accelerations in the Grim clock, Hannum clock, phenotypic clock, and skin and blood clock. CONCLUSIONS: Earlier age at menarche and higher BMI are associated with elevations in DNA methylation-based age estimates in healthy breast tissues, suggesting that cumulative estrogen exposure drives breast epigenetic aging. IMPACT: Epigenetic clock measures may help advance inquiry into the relationship between accelerated breast tissue aging and an elevated incidence of breast cancer in younger women.

Pregnancy Hypertension and a Commonly Inherited IGF1R Variant (rs2016347) Reduce Breast Cancer Risk by Enhancing Mammary Gland Involution.

BACKGROUND: Terminal duct lobular units (TDLUs) are the anatomic sites of breast cancer initiation, and breast tissue involution resulting in lower TDLU counts has been associated with decreased breast cancer risk. The insulin-like growth factor (IGF) pathway plays a role in breast involution, and systemic changes in this developmental pathway occur with hypertensive disorders of pregnancy (HDP), which have also been associated with lower breast cancer risk, especially in women carrying a functional variant of IGF1R SNP rs2016347. We proposed that this breast cancer protective effect might be explained by increased breast tissue involution. MATERIALS AND METHODS: We conducted a retrospective cohort study utilizing the Komen Tissue Bank, which collects breast tissue core biopsies from women without a history of breast cancer. Eighty white non-Hispanic women with a history of HDP were selected along with 120 nonexposed participants, and after genotyping for rs2016347, TDLU parameters were histologically measured blinded to participant characteristics from fixed biopsy sections. RESULTS: Stratified models by HDP status demonstrated that among HDP+ participants, those carrying two T alleles of rs2016347 had a decrease in TDLU counts of 53.2% when compared to those with no T alleles (p=0.049). Trend analysis demonstrated a multiplicative decrease in counts of 31.6% per T allele (p=0.050). Although no statistically significant interaction was seen between HDP status and T alleles, interaction terms showed increasingly negative values reaching a p value of 0.124 for HDP × 2T alleles. CONCLUSIONS: The observed statistically significant decrease in TDLU counts signifies increased breast epithelial involution in women with prior HDP who inherited the TT genotype of IGF1R SNP rs2016347. The increasing degree of breast involution with greater rs2016347 T allele copy number is consistent with the known progressive reduction in IGF1R expression in breast and other normal tissues.

The Susan G. Komen for the Cure Tissue Bank at the IU Simon Cancer Center: a unique resource for defining the "molecular histology" of the breast.

"Molecular histology" of the breast may be conceptualized as encompassing the normative ranges of histologic structure and marker expression in normal breast tissues in relation to a woman's age and life experiences. Studies of molecular histology can aid our understanding of early events in breast carcinogenesis and provide data for comparison with diseased breast tissues. Until recently, lack of epidemiologically annotated, optimally prepared normal breast tissues obtained from healthy women presented a barrier to breast cancer research. The Komen Tissue Bank at Indiana University (Indianapolis, IN) is a unique biorepository that was developed to overcome this limitation. The Bank enrolls healthy donors who provide questionnaire data, blood, and up to four breast biopsies, which are prepared as both formalin-fixed, paraffin-embedded and frozen tissues. The resource is accessible to researchers worldwide through a proposal submission, review, and approval process. As of November 2010, the Bank had collected specimens and information from 1,174 donors. In this review, we discuss the importance of studying normal breast tissues, assess the strengths and limitations of studying normal tissues obtained from different sources, and summarize the features of the Komen Tissue Bank. As research projects are completed, results will be posted on the Bank's website.