Luminal epithelial cells integrate variable responses to aging into stereotypical changes that underlie breast cancer susceptibility

  1. Rosalyn W Sayaman  Is a corresponding author
  2. Masaru Miyano
  3. Eric G Carlson
  4. Parijat Senapati
  5. Arrianna Zirbes
  6. Sundus F Shalabi
  7. Michael E Todhunter
  8. Victoria E Seewaldt
  9. Susan L Neuhausen
  10. Martha R Stampfer
  11. Dustin E Schones
  12. Mark LaBarge  Is a corresponding author
  1. University of California, San Francisco, United States
  2. City Of Hope National Medical Center, United States
  3. Lawrence Berkeley National Laboratory, United States

Abstract

Effects from aging in single cells are heterogenous, whereas at the organ- and tissue-levels aging phenotypes tend to appear as stereotypical changes. The mammary epithelium is a bilayer of two major phenotypically and functionally distinct cell lineages: luminal epithelial and myoepithelial cells. Mammary luminal epithelia exhibit substantial stereotypical changes with age that merit attention because these cells are the putative cells-of-origin for breast cancers. We hypothesize that effects from aging that impinge upon maintenance of lineage fidelity increase susceptibility to cancer initiation. We generated and analyzed transcriptomes from primary luminal epithelial and myoepithelial cells from younger <30 (y)ears old and older >55y women. In addition to age-dependent directional changes in gene expression, we observed increased transcriptional variance with age that contributed to genome-wide loss of lineage fidelity. Age-dependent variant responses were common to both lineages, whereas directional changes were almost exclusively detected in luminal epithelia and involved altered regulation of chromatin and genome organizers such as SATB1. Epithelial expression of gap junction protein GJB6 increased with age, and modulation of GJB6 expression in heterochronous co-cultures revealed that it provided a communication conduit from myoepithelial cells that drove directional change in luminal cells. Age-dependent luminal transcriptomes comprised a prominent signal that could be detected in bulk tissue during aging and transition into cancers. A machine learning classifier based on luminal-specific aging distinguished normal from cancer tissue and was highly predictive of breast cancer subtype. We speculate that luminal epithelia are the ultimate site of integration of the variant responses to aging in their surrounding tissue, and that their emergent phenotype both endows cells with the ability to become cancer-cells-of-origin and represents a biosensor that presages cancer susceptibility.

Data availability

The datasets generated and analyzed during the current study include RNA-sequencing count data publicly available as part of GSE182338 (Miyano et al., 2021; Sayaman et al., 2021; Sayaman et al., 2022; Shalabi et al., 2021; Todhunter et al., 2021). The gene expression data that support the findings of this study are available from GSE102088 (Song et al., 2017); GSE81540 (Brueffer et al., 2020; Brueffer et al., 2018; Dahlgren et al., 2021); The Cancer Genome Atlas (TCGA) Research Network: https://www.cancer.gov/tcga; and The Genotype-Tissue Expression (GTEx) Project: https://gtexportal.org/. Single-cell RNAseq data sets used for validation are available from GSE161529 (Pal et al., 2021); GSE174588 (Nee et al., 2023); and GSE198732 (Murrow et al., 2022). Analysis was conducted using standard R/Bioconductor packages and statistical tests implemented in R. All package versions, model design, and parameters are described in detail in Supplementary Methods. Summary statistics at defined significance levels are provided in Supplementary Tables; full summary statistics are provided via FigShare (https://figshare.com/s/2a7ceffccfe3f35f3ce8).

The following data sets were generated
The following previously published data sets were used

Article and author information

Author details

  1. Rosalyn W Sayaman

    Department of Laboratory Medicine, University of California, San Francisco, San Francisco, United States
    For correspondence
    rwsayaman@gmail.com
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-1343-0619
  2. Masaru Miyano

    Department of Population Sciences, City Of Hope National Medical Center, Duarte, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-1490-4743
  3. Eric G Carlson

    Department of Population Sciences, City Of Hope National Medical Center, Duarte, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-2784-1729
  4. Parijat Senapati

    Department of Diabetes Complications and Metabolism, City Of Hope National Medical Center, Duarte, United States
    Competing interests
    The authors declare that no competing interests exist.
  5. Arrianna Zirbes

    Department of Population Sciences, City Of Hope National Medical Center, Duarte, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-3849-2616
  6. Sundus F Shalabi

    Department of Population Sciences, City Of Hope National Medical Center, Duarte, United States
    Competing interests
    The authors declare that no competing interests exist.
  7. Michael E Todhunter

    Department of Population Sciences, City Of Hope National Medical Center, Duarte, United States
    Competing interests
    The authors declare that no competing interests exist.
  8. Victoria E Seewaldt

    Department of Population Sciences, City Of Hope National Medical Center, Duarte, United States
    Competing interests
    The authors declare that no competing interests exist.
  9. Susan L Neuhausen

    Department of Population Sciences, City Of Hope National Medical Center, Duarte, United States
    Competing interests
    The authors declare that no competing interests exist.
  10. Martha R Stampfer

    Biological Sciences and Engineering, Lawrence Berkeley National Laboratory, Berkeley, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-3801-5086
  11. Dustin E Schones

    Department of Diabetes Complications and Metabolism, City Of Hope National Medical Center, Duarte, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-7692-8583
  12. Mark LaBarge

    Department of Population Sciences, City Of Hope National Medical Center, Duarte, United States
    For correspondence
    mlabarge@coh.org
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-2405-4719

Funding

National Cancer Institute (T32CA221709)

  • Rosalyn W Sayaman

Yvonne Craig-Aldrich Fund for Cancer Research

  • Mark LaBarge

City of Hope Center for Cancer and Aging

  • Mark LaBarge

National Cancer Institute (K01CA279498)

  • Rosalyn W Sayaman

National Cancer Institute (U01CA244109)

  • Mark LaBarge

National Institute on Aging (R33AG059206)

  • Mark LaBarge

National Institute of Biomedical Imaging and Bioengineering (R01EB024989)

  • Mark LaBarge

National Cancer Institute (R01CA237602)

  • Mark LaBarge

Department of Defense (BC141351)

  • Mark LaBarge

Department of Defense (BC181737)

  • Mark LaBarge

Conrad N. Hilton Foundation

  • Mark LaBarge

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Ethics

Human subjects: The FACS-enriched luminal epithelial (LEP) and MEP myoepithelial (MEP) cells from finite lifespan, non-immortalized human mammary epithelial cells (HMECs) grown to 4th passage serve as our experimental model system. HMECs were cultured from breast tissue organoids collected from reduction mammoplasties (RM). Deidentified surgical discard tissue was obtained non-prospectively with consent for research and publication under a Lawrence Berkeley National Laboratory (Berkeley, CA) approved IRB 22997, or from City of Hope under IRB 17185 for sample distribution and collection.

Copyright

This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.

Metrics

  • 242
    views
  • 49
    downloads
  • 0
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Rosalyn W Sayaman
  2. Masaru Miyano
  3. Eric G Carlson
  4. Parijat Senapati
  5. Arrianna Zirbes
  6. Sundus F Shalabi
  7. Michael E Todhunter
  8. Victoria E Seewaldt
  9. Susan L Neuhausen
  10. Martha R Stampfer
  11. Dustin E Schones
  12. Mark LaBarge
(2024)
Luminal epithelial cells integrate variable responses to aging into stereotypical changes that underlie breast cancer susceptibility
eLife 13:e95720.
https://doi.org/10.7554/eLife.95720

Share this article

https://doi.org/10.7554/eLife.95720

Further reading

    1. Cancer Biology
    Bruno Bockorny, Lakshmi Muthuswamy ... Senthil K Muthuswamy
    Tools and Resources

    Pancreatic cancer has the worst prognosis of all common tumors. Earlier cancer diagnosis could increase survival rates and better assessment of metastatic disease could improve patient care. As such, there is an urgent need to develop biomarkers to diagnose this deadly malignancy. Analyzing circulating extracellular vesicles (cEVs) using ‘liquid biopsies’ offers an attractive approach to diagnose and monitor disease status. However, it is important to differentiate EV-associated proteins enriched in patients with pancreatic ductal adenocarcinoma (PDAC) from those with benign pancreatic diseases such as chronic pancreatitis and intraductal papillary mucinous neoplasm (IPMN). To meet this need, we combined the novel EVtrap method for highly efficient isolation of EVs from plasma and conducted proteomics analysis of samples from 124 individuals, including patients with PDAC, benign pancreatic diseases and controls. On average, 912 EV proteins were identified per 100 µL of plasma. EVs containing high levels of PDCD6IP, SERPINA12, and RUVBL2 were associated with PDAC compared to the benign diseases in both discovery and validation cohorts. EVs with PSMB4, RUVBL2, and ANKAR were associated with metastasis, and those with CRP, RALB, and CD55 correlated with poor clinical prognosis. Finally, we validated a seven EV protein PDAC signature against a background of benign pancreatic diseases that yielded an 89% prediction accuracy for the diagnosis of PDAC. To our knowledge, our study represents the largest proteomics profiling of circulating EVs ever conducted in pancreatic cancer and provides a valuable open-source atlas to the scientific community with a comprehensive catalogue of novel cEVs that may assist in the development of biomarkers and improve the outcomes of patients with PDAC.

    1. Cancer Biology
    2. Evolutionary Biology
    Lingjie Zhang, Tong Deng ... Chung-I Wu
    Research Article

    Tumorigenesis, like most complex genetic traits, is driven by the joint actions of many mutations. At the nucleotide level, such mutations are cancer-driving nucleotides (CDNs). The full sets of CDNs are necessary, and perhaps even sufficient, for the understanding and treatment of each cancer patient. Currently, only a small fraction of CDNs is known as most mutations accrued in tumors are not drivers. We now develop the theory of CDNs on the basis that cancer evolution is massively repeated in millions of individuals. Hence, any advantageous mutation should recur frequently and, conversely, any mutation that does not is either a passenger or deleterious mutation. In the TCGA cancer database (sample size n=300–1000), point mutations may recur in i out of n patients. This study explores a wide range of mutation characteristics to determine the limit of recurrences (i*) driven solely by neutral evolution. Since no neutral mutation can reach i*=3, all mutations recurring at i≥3 are CDNs. The theory shows the feasibility of identifying almost all CDNs if n increases to 100,000 for each cancer type. At present, only <10% of CDNs have been identified. When the full sets of CDNs are identified, the evolutionary mechanism of tumorigenesis in each case can be known and, importantly, gene targeted therapy will be far more effective in treatment and robust against drug resistance.