A human breast atlas integrating single-cell proteomics and transcriptomics.

The breast is a dynamic organ whose response to physiological and pathophysiological conditions alters its disease susceptibility, yet the specific effects of these clinical variables on cell state remain poorly annotated. We present a unified, high-resolution breast atlas by integrating single-cell RNA-seq, mass cytometry, and cyclic immunofluorescence, encompassing a myriad of states. We define cell subtypes within the alveolar, hormone-sensing, and basal epithelial lineages, delineating associations of several subtypes with cancer risk factors, including age, parity, and BRCA2 germline mutation. Of particular interest is a subset of alveolar cells termed basal-luminal (BL) cells, which exhibit poor transcriptional lineage fidelity, accumulate with age, and carry a gene signature associated with basal-like breast cancer. We further utilize a medium-depletion approach to identify molecular factors regulating cell-subtype proportion in organoids. Together, these data are a rich resource to elucidate diverse mammary cell states.

Aging-Associated Alterations in Mammary Epithelia and Stroma Revealed by Single-Cell RNA Sequencing.

Aging is closely associated with increased susceptibility to breast cancer, yet there have been limited systematic studies of aging-induced alterations in the mammary gland. Here, we leverage high-throughput single-cell RNA sequencing to generate a detailed transcriptomic atlas of young and aged murine mammary tissues. By analyzing epithelial, stromal, and immune cells, we identify age-dependent alterations in cell proportions and gene expression, providing evidence that suggests alveolar maturation and physiological decline. The analysis also uncovers potential pro-tumorigenic mechanisms coupled to the age-associated loss of tumor suppressor function and change in microenvironment. In addition, we identify a rare, age-dependent luminal population co-expressing hormone-sensing and secretory-alveolar lineage markers, as well as two macrophage populations expressing distinct gene signatures, underscoring the complex heterogeneity of the mammary epithelia and stroma. Collectively, this rich single-cell atlas reveals the effects of aging on mammary physiology and can serve as a useful resource for understanding aging-associated cancer risk.

Organoid cultures from normal and cancer-prone human breast tissues preserve complex epithelial lineages.

Recently, organoid technology has been used to generate a large repository of breast cancer organoids. Here we present an extensive evaluation of the ability of organoid culture technology to preserve complex stem/progenitor and differentiated cell types via long-term propagation of normal human mammary tissues. Basal/stem and luminal progenitor cells can differentiate in culture to generate mature basal and luminal cell types, including ER+ cells that have been challenging to maintain in culture. Cells associated with increased cancer risk can also be propagated. Single-cell analyses of matched organoid cultures and native tissues by mass cytometry for 38 markers provide a higher resolution representation of the multiple mammary epithelial cell types in the organoids, and demonstrate that protein expression patterns of the tissue of origin can be preserved in culture. These studies indicate that organoid cultures provide a valuable platform for studies of mammary differentiation, transformation, and breast cancer risk.

Optofluidic real-time cell sorter for longitudinal CTC studies in mouse models of cancer.

Circulating tumor cells (CTCs) play a fundamental role in cancer progression. However, in mice, limited blood volume and the rarity of CTCs in the bloodstream preclude longitudinal, in-depth studies of these cells using existing liquid biopsy techniques. Here, we present an optofluidic system that continuously collects fluorescently labeled CTCs from a genetically engineered mouse model (GEMM) for several hours per day over multiple days or weeks. The system is based on a microfluidic cell sorting chip connected serially to an unanesthetized mouse via an implanted arteriovenous shunt. Pneumatically controlled microfluidic valves capture CTCs as they flow through the device, and CTC-depleted blood is returned back to the mouse via the shunt. To demonstrate the utility of our system, we profile CTCs isolated longitudinally from animals over 4 days of treatment with the BET inhibitor JQ1 using single-cell RNA sequencing (scRNA-Seq) and show that our approach eliminates potential biases driven by intermouse heterogeneity that can occur when CTCs are collected across different mice. The CTC isolation and sorting technology presented here provides a research tool to help reveal details of how CTCs evolve over time, allowing studies to credential changes in CTCs as biomarkers of drug response and facilitating future studies to understand the role of CTCs in metastasis.

Quantitative proteomics identify Tenascin-C as a promoter of lung cancer progression and contributor to a signature prognostic of patient survival.

The extracellular microenvironment is an integral component of normal and diseased tissues that is poorly understood owing to its complexity. To investigate the contribution of the microenvironment to lung fibrosis and adenocarcinoma progression, two pathologies characterized by excessive stromal expansion, we used mouse models to characterize the extracellular matrix (ECM) composition of normal lung, fibrotic lung, lung tumors, and metastases. Using quantitative proteomics, we identified and assayed the abundance of 113 ECM proteins, which revealed robust ECM protein signatures unique to fibrosis, primary tumors, or metastases. These analyses indicated significantly increased abundance of several S100 proteins, including Fibronectin and Tenascin-C (Tnc), in primary lung tumors and associated lymph node metastases compared with normal tissue. We further showed that Tnc expression is repressed by the transcription factor Nkx2-1, a well-established suppressor of metastatic progression. We found that increasing the levels of Tnc, via CRISPR-mediated transcriptional activation of the endogenous gene, enhanced the metastatic dissemination of lung adenocarcinoma cells. Interrogation of human cancer gene expression data revealed that high TNC expression correlates with worse prognosis for lung adenocarcinoma, and that a three-gene expression signature comprising TNC, S100A10, and S100A11 is a robust predictor of patient survival independent of age, sex, smoking history, and mutational load. Our findings suggest that the poorly understood ECM composition of the fibrotic and tumor microenvironment is an underexplored source of diagnostic markers and potential therapeutic targets for cancer patients.

Clonal dynamics following p53 loss of heterozygosity in Kras-driven cancers.

Although it has become increasingly clear that cancers display extensive cellular heterogeneity, the spatial growth dynamics of genetically distinct clones within developing solid tumours remain poorly understood. Here we leverage mosaic analysis with double markers (MADM) to trace subclonal populations retaining or lacking p53 within oncogenic Kras-initiated lung and pancreatic tumours. In both models, p53 constrains progression to advanced adenocarcinomas. Comparison of lineage-related p53 knockout and wild-type clones reveals a minor role of p53 in suppressing cell expansion in lung adenomas. In contrast, p53 loss promotes both the initiation and expansion of low-grade pancreatic intraepithelial neoplasia (PanINs), likely through differential expression of the p53 regulator p19ARF. Strikingly, lineage-related cells are often dispersed in lung adenomas and PanINs, contrasting with more contiguous growth of advanced subclones. Together, these results support cancer type-specific suppressive roles of p53 in early tumour progression and offer insights into clonal growth patterns during tumour development.

Foxa2 and Cdx2 cooperate with Nkx2-1 to inhibit lung adenocarcinoma metastasis.

Despite the fact that the majority of lung cancer deaths are due to metastasis, the molecular mechanisms driving metastatic progression are poorly understood. Here, we present evidence that loss of Foxa2 and Cdx2 synergizes with loss of Nkx2-1 to fully activate the metastatic program. These three lineage-specific transcription factors are consistently down-regulated in metastatic cells compared with nonmetastatic cells. Knockdown of these three factors acts synergistically and is sufficient to promote the metastatic potential of nonmetastatic cells to that of naturally arising metastatic cells in vivo. Furthermore, silencing of these three transcription factors is sufficient to account for a significant fraction of the gene expression differences between the nonmetastatic and metastatic states in lung adenocarcinoma, including up-regulated expression of the invadopodia component Tks5long, the embryonal proto-oncogene Hmga2, and the epithelial-to-mesenchymal mediator Snail. Finally, analyses of tumors from a genetically engineered mouse model and patients show that low expression of Nkx2-1, Foxa2, and Cdx2 strongly correlates with more advanced tumors and worse survival. Our findings reveal that a large part of the complex transcriptional network in metastasis can be controlled by a small number of regulatory nodes that function redundantly, and loss of multiple nodes is required to fully activate the metastatic program.

Differential Tks5 isoform expression contributes to metastatic invasion of lung adenocarcinoma.

Metastasis accounts for the vast majority of cancer-related deaths, yet the molecular mechanisms that drive metastatic spread remain poorly understood. Here we report that Tks5, which has been linked to the formation of proteolytic cellular protrusions known as invadopodia, undergoes an isoform switch during metastatic progression in a genetically engineered mouse model of lung adenocarcinoma. Nonmetastatic primary tumor-derived cells predominantly expressed a short isoform, Tks5short, while metastatic primary tumor- and metastasis-derived cells acquired increased expression of the full-length isoform Tks5long. This elevation of Tks5long to Tks5short ratio correlated with a commensurate increase in invadopodia activity in metastatic cells compared with nonmetastatic cells. Further characterization of these isoforms by knockdown and overexpression experiments demonstrated that Tks5long promoted invadopodia in vitro and increased metastasis in transplant models and an autochthonous model of lung adenocarcinoma. Conversely, Tks5short decreased invadopodia stability and proteolysis, acting as a natural dominant-negative inhibitor to Tks5long. Importantly, high Tks5long and low Tks5short expressions in human lung adenocarcinomas correlated with metastatic disease and predicted worse survival of early stage patients. These data indicate that tipping the Tks5 isoform balance to a high Tks5long to Tks5short ratio promotes invadopodia-mediated invasion and metastasis.