Single-Cell Transcriptomics Reveals a Conserved Metaplasia Program in Pancreatic Injury.

BACKGROUND & AIMS: Acinar to ductal metaplasia (ADM) occurs in the pancreas in response to tissue injury and is a potential precursor for adenocarcinoma. The goal of these studies was to define the populations arising from ADM, the associated transcriptional changes, and markers of disease progression. METHODS: Acinar cells were lineage-traced with enhanced yellow fluorescent protein (EYFP) to follow their fate post-injury. Transcripts of more than 13,000 EYFP+ cells were determined using single-cell RNA sequencing (scRNA-seq). Developmental trajectories were generated. Data were compared with gastric metaplasia, KrasG12D-induced neoplasia, and human pancreatitis. Results were confirmed by immunostaining and electron microscopy. KrasG12D was expressed in injury-induced ADM using several inducible Cre drivers. Surgical specimens of chronic pancreatitis from 15 patients were evaluated by immunostaining. RESULTS: scRNA-seq of ADM revealed emergence of a mucin/ductal population resembling gastric pyloric metaplasia. Lineage trajectories suggest that some pyloric metaplasia cells can generate tuft and enteroendocrine cells (EECs). Comparison with KrasG12D-induced ADM identifies populations associated with disease progression. Activation of KrasG12D expression in HNF1B+ or POU2F3+ ADM populations leads to neoplastic transformation and formation of MUC5AC+ gastric-pit-like cells. Human pancreatitis samples also harbor pyloric metaplasia with a similar transcriptional phenotype. CONCLUSIONS: Under conditions of chronic injury, acinar cells undergo a pyloric-type metaplasia to mucinous progenitor-like populations, which seed disparate tuft cell and EEC lineages. ADM-derived EEC subtypes are diverse. KrasG12D expression is sufficient to drive neoplasia when targeted to injury-induced ADM populations and offers an alternative origin for tumorigenesis. This program is conserved in human pancreatitis, providing insight into early events in pancreas diseases.

Tuft Cells Inhibit Pancreatic Tumorigenesis in Mice by Producing Prostaglandin D2.

BACKGROUND & AIMS: Development of pancreatic ductal adenocarcinoma (PDA) involves acinar to ductal metaplasia and genesis of tuft cells. It has been a challenge to study these rare cells because of the lack of animal models. We investigated the role of tuft cells in pancreatic tumorigenesis. METHODS: We performed studies with LSL-KrasG12D/+;Ptf1aCre/+ mice (KC; develop pancreatic tumors), KC mice crossed with mice with pancreatic disruption of Pou2f3 (KPouC mice; do not develop tuft cells), or mice with pancreatic disruption of the hematopoietic prostaglandin D synthase gene (Hpgds, KHC mice) and wild-type mice. Mice were allowed to age or were given caerulein to induce pancreatitis; pancreata were collected and analyzed by histology, immunohistochemistry, RNA sequencing, ultrastructural microscopy, and metabolic profiling. We performed laser-capture dissection and RNA-sequencing analysis of pancreatic tissues from 26 patients with pancreatic intraepithelial neoplasia (PanIN), 19 patients with intraductal papillary mucinous neoplasms (IPMNs), and 197 patients with PDA. RESULTS: Pancreata from KC mice had increased formation of tuft cells and higher levels of prostaglandin D2 than wild-type mice. Pancreas-specific deletion of POU2F3 in KC mice (KPouC mice) resulted in a loss of tuft cells and accelerated tumorigenesis. KPouC mice had increased fibrosis and activation of immune cells after administration of caerulein. Pancreata from KPouC and KHC mice had significantly lower levels of prostaglandin D2, compared with KC mice, and significantly increased numbers of PanINs and PDAs. KPouC and KHC mice had increased pancreatic injury after administration of caerulein, significantly less normal tissue, more extracellular matrix deposition, and higher PanIN grade than KC mice. Human PanIN and intraductal papillary mucinous neoplasm had gene expression signatures associated with tuft cells and increased expression of Hpgds messenger RNA compared with PDA. CONCLUSIONS: In mice with KRAS-induced pancreatic tumorigenesis, loss of tuft cells accelerates tumorigenesis and increases the severity of caerulein-induced pancreatic injury, via decreased production of prostaglandin D2. These data are consistent with the hypothesis that tuft cells are a metaplasia-induced tumor attenuating cell type.

Reactivation of multipotency by oncogenic PIK3CA induces breast tumour heterogeneity.

Breast cancer is the most frequent cancer in women and consists of heterogeneous types of tumours that are classified into different histological and molecular subtypes. PIK3CA and P53 (also known as TP53) are the two most frequently mutated genes and are associated with different types of human breast cancers. The cellular origin and the mechanisms leading to PIK3CA-induced tumour heterogeneity remain unknown. Here we used a genetic approach in mice to define the cellular origin of Pik3ca-derived tumours and the impact of mutations in this gene on tumour heterogeneity. Surprisingly, oncogenic Pik3ca(H1047R) mutant expression at physiological levels in basal cells using keratin (K)5-CreER(T2) mice induced the formation of luminal oestrogen receptor (ER)-positive/progesterone receptor (PR)-positive tumours, while its expression in luminal cells using K8-CReER(T2) mice gave rise to luminal ER(+)PR(+) tumours or basal-like ER(-)PR(-) tumours. Concomitant deletion of p53 and expression of Pik3ca(H1047R) accelerated tumour development and induced more aggressive mammary tumours. Interestingly, expression of Pik3ca(H1047R) in unipotent basal cells gave rise to luminal-like cells, while its expression in unipotent luminal cells gave rise to basal-like cells before progressing into invasive tumours. Transcriptional profiling of cells that underwent cell fate transition upon Pik3ca(H1047R) expression in unipotent progenitors demonstrated a profound oncogene-induced reprogramming of these newly formed cells and identified gene signatures characteristic of the different cell fate switches that occur upon Pik3ca(H1047R) expression in basal and luminal cells, which correlated with the cell of origin, tumour type and different clinical outcomes. Altogether our study identifies the cellular origin of Pik3ca-induced tumours and reveals that oncogenic Pik3ca(H1047R) activates a multipotent genetic program in normally lineage-restricted populations at the early stage of tumour initiation, setting the stage for future intratumoural heterogeneity. These results have important implications for our understanding of the mechanisms controlling tumour heterogeneity and the development of new strategies to block PIK3CA breast cancer initiation.

The influence of tamoxifen on normal mouse mammary gland homeostasis.

INTRODUCTION: Lineage tracing using inducible genetic labeling has emerged to be a powerful method for interrogating the developmental fate of cells in intact tissues. A common induction mechanism is the use of tamoxifen-dependent Cre recombinase (CreER and CreERT2), but the effects of tamoxifen at doses normally used in lineage-tracing studies on normal adult mammary gland homeostasis are not known. METHODS: We used flow cytometry and immunostaining of intact glands to determine whether varying doses of tamoxifen skew the distribution and the apoptosis and proliferation status of different types of mammary epithelial cells in vivo. We also examined how tamoxifen influences the number of progenitor and mammary repopulating units (MRUs). RESULTS: Our results indicate that ≥5 mg/25 g body weight of tamoxifen induces a transient increase in cell proliferation and in the number of basal cells in the adult mammary epithelium up to 7 days after tamoxifen administration. However, in the medium term (3 weeks), all doses of tamoxifen≥1 mg/25 g body weight result in a decrease in the number of basal and EpCAM+CD49b- luminal cells and a decrease in progenitor cell function. Tamoxifen at doses≥5 mg/25 g body weight induced a transient increase in caspase-3-mediated apoptotic cell death within the mammary epithelium. However, mammary epithelial cell numbers in all subpopulations were restored to their original levels by 8 weeks. No long-lasting effects of tamoxifen on MRU numbers or on pubertal ductal development were observed. CONCLUSION: Tamoxifen can skew the distribution of mammary cell types in a dose-dependent manner, and thus caution must be taken when interpreting lineage-tracing studies using high doses of tamoxifen, particularly when short-duration analyses of a quantitative nature are being performed.

Cancer Stem Cells: Concepts, Challenges, and Opportunities for Cancer Therapy.

Cancer stem cells (CSCs) are a subpopulation of cancer cells with self-renewal capacity, that fuel tumor growth and contribute to the heterogeneous nature of tumors. First identified in hematological malignancies, CSC populations have to date been proposed in solid tumors in various organs. In vitro and in vivo assays, mouse genetic models, and more recently single-cell sequencing technologies and other '-omics' methodologies have not only facilitated the identification of novel CSC populations but also revealed and clarified novel properties of CSCs. Increasingly, both cell-autonomous and CSC niche factors are recognized as important contributors of CSC properties. The deepened understanding of CSC properties and characteristics would enable and facilitate the rational design of CSC-specific therapeutics that would, ideally, have high selectivity for cancer cells, eliminate tumor bulk, and prevent tumor recurrence. Addressing these issues would form some of the key challenges of the CSC research field in the coming years.

Single-Cell Chromatin Analysis of Mammary Gland Development Reveals Cell-State Transcriptional Regulators and Lineage Relationships.

Technological improvements enable single-cell epigenetic analyses of organ development. We reasoned that high-resolution single-cell chromatin accessibility mapping would provide needed insight into the epigenetic reprogramming and transcriptional regulators involved in normal mammary gland development. Here, we provide a single-cell resource of chromatin accessibility for murine mammary development from the peak of fetal mammary stem cell (fMaSC) functional activity in late embryogenesis to the differentiation of adult basal and luminal cells. We find that the chromatin landscape within individual cells predicts both gene accessibility and transcription factor activity. The ability of single-cell chromatin profiling to separate E18 fetal mammary cells into clusters exhibiting basal-like and luminal-like chromatin features is noteworthy. Such distinctions were not evident in analyses of droplet-based single-cell transcriptomic data. We present a web application as a scientific resource for facilitating future analyses of the gene regulatory networks involved in mammary development.

Single-Cell Transcriptomes Distinguish Stem Cell State Changes and Lineage Specification Programs in Early Mammary Gland Development.

The mammary gland consists of cells with gene expression patterns reflecting their cellular origins, function, and spatiotemporal context. However, knowledge of developmental kinetics and mechanisms of lineage specification is lacking. We address this significant knowledge gap by generating a single-cell transcriptome atlas encompassing embryonic, postnatal, and adult mouse mammary development. From these data, we map the chronology of transcriptionally and epigenetically distinct cell states and distinguish fetal mammary stem cells (fMaSCs) from their precursors and progeny. fMaSCs show balanced co-expression of factors associated with discrete adult lineages and a metabolic gene signature that subsides during maturation but reemerges in some human breast cancers and metastases. These data provide a useful resource for illuminating mammary cell heterogeneity, the kinetics of differentiation, and developmental correlates of tumorigenesis.

Stem and progenitor cell division kinetics during postnatal mouse mammary gland development.

The cycling properties of mammary stem and progenitor cells is not well understood. To determine the division properties of these cells, we administered synthetic nucleosides for varying periods of time to mice at different stages of postnatal development and monitored the rate of uptake of these nucleosides in the different mammary cell compartments. Here we show that most cell division in the adult virgin gland is restricted to the oestrogen receptor-expressing luminal cell lineage. Our data also demonstrate that the oestrogen receptor-expressing, milk and basal cell subpopulations have telomere lengths and cell division kinetics that are not compatible with these cells being hierarchically organized; instead, our data indicate that in the adult homeostatic gland, each cell type is largely maintained by its own restricted progenitors. We also observe that transplantable stem cells are largely quiescent during oestrus, but are cycling during dioestrus when progesterone levels are high.

Mammary stem cells have myoepithelial cell properties.

Contractile myoepithelial cells dominate the basal layer of the mammary epithelium and are considered to be differentiated cells. However, we observe that up to 54% of single basal cells can form colonies when seeded into adherent culture in the presence of agents that disrupt actin-myosin interactions, and on average, 65% of the single-cell-derived basal colonies can repopulate a mammary gland when transplanted in vivo. This indicates that a high proportion of basal myoepithelial cells can give rise to a mammary repopulating unit (MRU). We demonstrate that myoepithelial cells, flow-sorted using two independent myoepithelial-specific reporter strategies, have MRU capacity. Using an inducible lineage-tracing approach we follow the progeny of myoepithelial cells that express α-smooth muscle actin and show that they function as long-lived lineage-restricted stem cells in the virgin state and during pregnancy.