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.

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.

Single-cell RNA sequencing of developing mouse mammary epithelia reveals the timing of lineage specification. Giraddi et al. find that fetal mammary stem cells co-express factors that define distinct lineages in their progeny and bear functionally relevant metabolic program signatures that change with differentiation and are resurrected in human breast cancers and metastases.

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.

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.