Single cell transcriptome atlas of mouse mammary epithelial cells across development.

BACKGROUND: Heterogeneity within the mouse mammary epithelium and potential lineage relationships have been recently explored by single-cell RNA profiling. To further understand how cellular diversity changes during mammary ontogeny, we profiled single cells from nine different developmental stages spanning late embryogenesis, early postnatal, prepuberty, adult, mid-pregnancy, late-pregnancy, and post-involution, as well as the transcriptomes of micro-dissected terminal end buds (TEBs) and subtending ducts during puberty. METHODS: The single cell transcriptomes of 132,599 mammary epithelial cells from 9 different developmental stages were determined on the 10x Genomics Chromium platform, and integrative analyses were performed to compare specific time points. RESULTS: The mammary rudiment at E18.5 closely aligned with the basal lineage, while prepubertal epithelial cells exhibited lineage segregation but to a less differentiated state than their adult counterparts. Comparison of micro-dissected TEBs versus ducts showed that luminal cells within TEBs harbored intermediate expression profiles. Ductal basal cells exhibited increased chromatin accessibility of luminal genes compared to their TEB counterparts suggesting that lineage-specific chromatin is established within the subtending ducts during puberty. An integrative analysis of five stages spanning the pregnancy cycle revealed distinct stage-specific profiles and the presence of cycling basal, mixed-lineage, and 'late' alveolar intermediates in pregnancy. Moreover, a number of intermediates were uncovered along the basal-luminal progenitor cell axis, suggesting a continuum of alveolar-restricted progenitor states. CONCLUSIONS: This extended single cell transcriptome atlas of mouse mammary epithelial cells provides the most complete coverage for mammary epithelial cells during morphogenesis to date. Together with chromatin accessibility analysis of TEB structures, it represents a valuable framework for understanding developmental decisions within the mouse mammary gland.

The genetic inheritance of the blue-eyed white phenotype in alpacas (Vicugna pacos).

White-spotting patterns in mammals can be caused by mutations in the gene KIT, whose protein is necessary for the normal migration and survival of melanocytes from the neural crest. The alpaca (Vicugna pacos) blue-eyed white (BEW) phenotype is characterized by 2 blue eyes and a solid white coat over the whole body. Breeders hypothesize that the BEW phenotype in alpacas is caused by the combination of the gene causing gray fleece and a white-spotting gene. We performed an association study using KIT flanking and intragenic markers with 40 unrelated alpacas, of which 17 were BEW. Two microsatellite alleles at KIT-related markers were significantly associated (P < 0.0001) with the BEW phenotype (bew1 and bew2). In a larger cohort of 171 related individuals, we identify an abundance of an allele (bew1) in gray animals and the occurrence of bew2 homozygotes that are solid white with pigmented eyes. Association tests accounting for population structure and familial relatedness are consistent with a proposed model where these alleles are in linkage disequilibrium with a mutation or mutations that contribute to the BEW phenotype and to individual differences in fleece color.

Ehf controls mammary alveolar lineage differentiation and is a putative suppressor of breast tumorigenesis.

The transcription factor EHF is highly expressed in the lactating mammary gland, but its role in mammary development and tumorigenesis is not fully understood. Utilizing a mouse model of Ehf deletion, herein, we demonstrate that loss of Ehf impairs mammary lobuloalveolar differentiation at late pregnancy, indicated by significantly reduced levels of milk genes and milk lipids, fewer differentiated alveolar cells, and an accumulation of alveolar progenitor cells. Further, deletion of Ehf increased proliferative capacity and attenuated prolactin-induced alveolar differentiation in mammary organoids. Ehf deletion also increased tumor incidence in the MMTV-PyMT mammary tumor model and increased the proliferative capacity of mammary tumor organoids, while low EHF expression was associated with higher tumor grade and poorer outcome in luminal A and basal human breast cancers. Collectively, these findings establish EHF as a non-redundant regulator of mammary alveolar differentiation and a putative suppressor of mammary tumorigenesis.

Identification of aberrant luminal progenitors and mTORC1 as a potential breast cancer prevention target in BRCA2 mutation carriers.

Inheritance of a BRCA2 pathogenic variant conveys a substantial life-time risk of breast cancer. Identification of the cell(s)-of-origin of BRCA2-mutant breast cancer and targetable perturbations that contribute to transformation remains an unmet need for these individuals who frequently undergo prophylactic mastectomy. Using preneoplastic specimens from age-matched, premenopausal females, here we show broad dysregulation across the luminal compartment in BRCA2mut/+ tissue, including expansion of aberrant ERBB3lo luminal progenitor and mature cells, and the presence of atypical oestrogen receptor (ER)-positive lesions. Transcriptional profiling and functional assays revealed perturbed proteostasis and translation in ERBB3lo progenitors in BRCA2mut/+ breast tissue, independent of ageing. Similar molecular perturbations marked tumours bearing BRCA2-truncating mutations. ERBB3lo progenitors could generate both ER+ and ER- cells, potentially serving as cells-of-origin for ER-positive or triple-negative cancers. Short-term treatment with an mTORC1 inhibitor substantially curtailed tumorigenesis in a preclinical model of BRCA2-deficient breast cancer, thus uncovering a potential prevention strategy for BRCA2 mutation carriers.

Chromosome-Level Alpaca Reference Genome VicPac3.1 Improves Genomic Insight Into the Biology of New World Camelids.

The development of high-quality chromosomally assigned reference genomes constitutes a key feature for understanding genome architecture of a species and is critical for the discovery of the genetic blueprints of traits of biological significance. South American camelids serve people in extreme environments and are important fiber and companion animals worldwide. Despite this, the alpaca reference genome lags far behind those available for other domestic species. Here we produced a chromosome-level improved reference assembly for the alpaca genome using the DNA of the same female Huacaya alpaca as in previous assemblies. We generated 190X Illumina short-read, 8X Pacific Biosciences long-read and 60X Dovetail Chicago® chromatin interaction scaffolding data for the assembly, used testis and skin RNAseq data for annotation, and cytogenetic map data for chromosomal assignments. The new assembly VicPac3.1 contains 90% of the alpaca genome in just 103 scaffolds and 76% of all scaffolds are mapped to the 36 pairs of the alpaca autosomes and the X chromosome. Preliminary annotation of the assembly predicted 22,462 coding genes and 29,337 isoforms. Comparative analysis of selected regions of the alpaca genome, such as the major histocompatibility complex (MHC), the region involved in the Minute Chromosome Syndrome (MCS) and candidate genes for high-altitude adaptations, reveal unique features of the alpaca genome. The alpaca reference genome VicPac3.1 presents a significant improvement in completeness, contiguity and accuracy over VicPac2 and is an important tool for the advancement of genomics research in all New World camelids.

Intraclonal Plasticity in Mammary Tumors Revealed through Large-Scale Single-Cell Resolution 3D Imaging.

Breast tumors are inherently heterogeneous, but the evolving cellular organization through neoplastic progression is poorly understood. Here we report a rapid, large-scale single-cell resolution 3D imaging protocol based on a one-step clearing agent that allows visualization of normal tissue architecture and entire tumors at cellular resolution. Imaging of multicolor lineage-tracing models of breast cancer targeted to either basal or luminal progenitor cells revealed profound clonal restriction during progression. Expression profiling of clones arising in Pten/Trp53-deficient tumors identified distinct molecular signatures. Strikingly, most clones harbored cells that had undergone an epithelial-to-mesenchymal transition, indicating widespread, inherent plasticity. Hence, an integrative pipeline that combines lineage tracing, 3D imaging, and clonal RNA sequencing technologies offers a comprehensive path for studying mechanisms underlying heterogeneity in whole tumors.

Foxp1 Is Indispensable for Ductal Morphogenesis and Controls the Exit of Mammary Stem Cells from Quiescence.

Long-lived quiescent mammary stem cells (MaSCs) are presumed to coordinate the dramatic expansion of ductal epithelium that occurs through the different phases of postnatal development, but little is known about the molecular regulators that underpin their activation. We show that ablation of the transcription factor Foxp1 in the mammary gland profoundly impairs ductal morphogenesis, resulting in a rudimentary tree throughout life. Foxp1-deficient glands were highly enriched for quiescent Tspan8hi MaSCs, which failed to become activated even in competitive transplantation assays, thus highlighting a cell-intrinsic defect. Foxp1 deletion also resulted in aberrant expression of basal genes in luminal cells, inferring a role in cell-fate decisions. Notably, Foxp1 was uncovered as a direct repressor of Tspan8 in basal cells, and deletion of Tspan8 rescued the defects in ductal morphogenesis elicited by Foxp1 loss. Thus, a single transcriptional regulator Foxp1 can control the exit of MaSCs from dormancy to orchestrate differentiation and development.

Identification of quiescent and spatially restricted mammary stem cells that are hormone responsive.

Despite accumulating evidence for a mammary differentiation hierarchy, the basal compartment comprising stem cells remains poorly characterized. Through gene expression profiling of Lgr5+ basal epithelial cells, we identify a new marker, Tetraspanin8 (Tspan8). Fractionation based on Tspan8 and Lgr5 expression uncovered three distinct mammary stem cell (MaSC) subsets in the adult mammary gland. These exist in a largely quiescent state but differ in their reconstituting ability, spatial localization, and their molecular and epigenetic signatures. Interestingly, the deeply quiescent MaSC subset (Lgr5+Tspan8hi) resides within the proximal region throughout life, and has a transcriptome strikingly similar to that of claudin-low tumours. Lgr5+Tspan8hi cells appear to originate from the embryonic mammary primordia before switching to a quiescent state postnatally but can be activated by ovarian hormones. Our findings reveal an unexpected degree of complexity within the adult MaSC compartment and identify a dormant subset poised for activation in response to physiological stimuli.

Construction of developmental lineage relationships in the mouse mammary gland by single-cell RNA profiling.

The mammary epithelium comprises two primary cellular lineages, but the degree of heterogeneity within these compartments and their lineage relationships during development remain an open question. Here we report single-cell RNA profiling of mouse mammary epithelial cells spanning four developmental stages in the post-natal gland. Notably, the epithelium undergoes a large-scale shift in gene expression from a relatively homogeneous basal-like program in pre-puberty to distinct lineage-restricted programs in puberty. Interrogation of single-cell transcriptomes reveals different levels of diversity within the luminal and basal compartments, and identifies an early progenitor subset marked by CD55. Moreover, we uncover a luminal transit population and a rare mixed-lineage cluster amongst basal cells in the adult mammary gland. Together these findings point to a developmental hierarchy in which a basal-like gene expression program prevails in the early post-natal gland prior to the specification of distinct lineage signatures, and the presence of cellular intermediates that may serve as transit or lineage-primed cells.