Evidence of steady-state fibroblast subtypes in the normal human breast as cells-of-origin for perturbed-state fibroblasts in breast cancer

BACKGROUND

Human breast cancer most frequently originates within a well-defined anatomical structure referred to as the terminal duct lobular unit (TDLU). This structure is endowed with its very own lobular fibroblasts representing one out of two steady-state fibroblast subtypes-the other being interlobular fibroblasts. While cancer-associated fibroblasts (CAFs) are increasingly appreciated as covering a spectrum of perturbed states, we lack a coherent understanding of their relationship-if any-with the steady-state fibroblast subtypes. To address this, we here established two autologous CAF lines representing inflammatory CAFs (iCAFs) and myofibroblast CAFs (myCAFs) and compared them with already established interlobular- and lobular fibroblasts with respect to their origin and impact on tumor formation.

METHODS

Primary breast tumor-derived CAFs were transduced to express human telomerase reverse transcriptase (hTERT) and sorted into CD105low and CD105high populations using fluorescence-activated cell sorting (FACS). The two populations were tested for differentiation similarities to iCAF and myCAF states through transcriptome-wide RNA-Sequencing (RNA-Seq) including comparison to an available iCAF-myCAF cell state atlas. Inference of origin in interlobular and lobular fibroblasts relied on RNA-Seq profiles, immunocytochemistry and growth characteristics. Osteogenic differentiation and bone formation assays in culture and in vivo were employed to gauge for origin in bone marrow-derived mesenchymal stem cells (bMSCs). Functional characteristics were assessed with respect to contractility in culture and interaction with tumor cells in mouse xenografts. The cells' gene expression signatures were tested for association with clinical outcome of breast cancer patients using survival data from The Cancer Genome Atlas database.

RESULTS

We demonstrate that iCAFs have properties in common with interlobular fibroblasts while myCAFs and lobular fibroblasts are related. None of the CAFs qualify as bMSCs as revealed by lack of critical performance in bone formation assays. Functionally, myCAFs and lobular fibroblasts are almost equally tumor promoting as opposed to iCAFs and interlobular fibroblasts. A myCAF gene signature is found to associate with poor breast cancer-specific survival.

CONCLUSIONS

We propose that iCAFs and myCAFs originate in interlobular and lobular fibroblasts, respectively, and more importantly, that the tumor-promoting properties of lobular fibroblasts render the TDLU an epicenter for breast cancer evolution.

Serum insulin-like growth factor (IGF)-I and IGF binding protein-3 in relation to terminal duct lobular unit involution of the normal breast in Caucasian and African American women: The Susan G. Komen Tissue Bank

Lesser degrees of terminal duct lobular unit (TDLU) involution, as reflected by higher numbers of TDLUs and acini/TDLU, are associated with elevated breast cancer risk. In rodent models, the insulin-like growth factor (IGF) system regulates involution of the mammary gland. We examined associations of circulating IGF measures with TDLU involution in normal breast tissues among women without precancerous lesions. Among 715 Caucasian and 283 African American (AA) women who donated normal breast tissue samples to the Komen Tissue Bank between 2009 and 2012 (75% premenopausal), serum concentrations of IGF-I and binding protein (IGFBP)-3 were quantified using enzyme-linked immunosorbent assay. Hematoxilyn and eosin-stained tissue sections were assessed for numbers of TDLUs ("TDLU count"). Zero-inflated Poisson regression models with a robust variance estimator were used to estimate relative risks (RRs) for association of IGF measures (tertiles) with TDLU count by race and menopausal status, adjusting for potential confounders. AA (vs. Caucasian) women had higher age-adjusted mean levels of serum IGF-I (137 vs. 131 ng/mL, p = 0.07) and lower levels of IGFBP-3 (4165 vs. 4684 ng/mL, p < 0.0001). Postmenopausal IGFBP-3 was inversely associated with TDLU count among AA (RR_T3vs.T1  = 0.49, 95% CI = 0.28-0.84, p-trend = 0.04) and Caucasian (RR_T3vs.T1 =0.64, 95% CI = 0.42-0.98, p-trend = 0.04) women. In premenopausal women, higher IGF-I:IGFBP-3 ratios were associated with higher TDLU count in Caucasian (RR_T3vs.T1 =1.33, 95% CI = 1.02-1.75, p-trend = 0.04), but not in AA (RR_T3vs.T1 =0.65, 95% CI = 0.42-1.00, p-trend = 0.05), women. Our data suggest a role of the IGF system, particularly IGFBP-3, in TDLU involution of the normal breast, a breast cancer risk factor, among Caucasian and AA women.

Ages at menarche- and menopause-related genetic variants in relation to terminal duct lobular unit involution in normal breast tissue

Reduced levels of terminal duct lobular unit (TDLU) involution, as reflected by higher numbers of TDLUs and acini per TDLU, have been associated with higher breast cancer risk. Younger age at menarche and older age at menopause have been previously related to lower levels of TDLU involution. To determine a possible genetic link, we examined whether single-nucleotide polymorphisms (SNPs) previously established in genome-wide association studies (GWAS) for ages at menarche and menopause are associated with TDLU involution. We conducted a pooled analysis of 862 women from two studies. H&E tissue sections were assessed for numbers of TDLUs and acini/TDLU. Poisson regression models were used to estimate associations of 36 menarche- and 21 menopause-SNPs with TDLU counts, acini counts/TDLU, and the product of these two measures, adjusting for age and study site. Fourteen percent of evaluated SNPs (eight SNPs) were associated with TDLU counts at p < 0.05, suggesting an enrichment of associations with TDLU counts. However, only menopause-SNPs had >50 % that were either significantly or nonsignificantly associated with TDLU measures in the directions consistent with their relationships shown in GWAS. Among ten SNPs that were statistically significantly associated with at least one TDLU involution measure (p < 0.05), seven SNPs (rs466639: RXRG; rs2243803: SLC14A2; rs2292573: GAB2; rs6438424: 3q13.32; rs7606918: METAP1D; rs11668344: TMEM150B; rs1635501: EXO1) were associated in the consistent directions. Our data suggest that the loci associated with ages at menarche and menopause may influence TDLU involution, suggesting some shared genetic mechanisms. However, larger studies are needed to confirm the results.