School of Mathematics and Statistics - Research Publications
Permanent URI for this collection
Search Results
1 - 10 of 12
-
ItemR code and downstream analysis objects for the scRNA-seq atlas of normal and tumorigenic human breast tissue(NATURE PORTFOLIO, 2022-03-23)Breast cancer is a common and highly heterogeneous disease. Understanding cellular diversity in the mammary gland and its surrounding micro-environment across different states can provide insight into cancer development in the human breast. Recently, we published a large-scale single-cell RNA expression atlas of the human breast spanning normal, preneoplastic and tumorigenic states. Single-cell expression profiles of nearly 430,000 cells were obtained from 69 distinct surgical tissue specimens from 55 patients. This article extends the study by providing quality filtering thresholds, downstream processed R data objects, complete cell annotation and R code to reproduce all the analyses. Data quality assessment measures are presented and details are provided for all the bioinformatic analyses that produced results described in the study.
-
ItemIn vivo genome-editing screen identifies tumor suppressor genes that cooperate with Trp53 loss during mammary tumorigenesis(WILEY, 2022-01-26)Breast cancer is a heterogeneous disease that comprises multiple histological and molecular subtypes. To gain insight into mutations that drive breast tumorigenesis, we describe a pipeline for the identification and validation of tumor suppressor genes. Based on an in vivo genome-wide CRISPR/Cas9 screen in Trp53+/- heterozygous mice, we identified tumor suppressor genes that included the scaffold protein Axin1, the protein kinase A regulatory subunit gene Prkar1a, as well as the proof-of-concept genes Pten, Nf1, and Trp53 itself. Ex vivo editing of primary mammary epithelial organoids was performed to further interrogate the roles of Axin1 and Prkar1a. Increased proliferation and profound changes in mammary organoid morphology were observed for Axin1/Trp53 and Prkar1a/Trp53 double mutants compared to Pten/Trp53 double mutants. Furthermore, direct in vivo genome editing via intraductal injection of lentiviruses engineered to express dual short-guide RNAs revealed that mutagenesis of Trp53 and either Prkar1a, Axin1, or Pten markedly accelerated tumor development compared to Trp53-only mutants. This proof-of-principle study highlights the application of in vivo CRISPR/Cas9 editing for uncovering cooperativity between defects in tumor suppressor genes that elicit mammary tumorigenesis.
-
ItemSingle cell transcriptome atlas of mouse mammary epithelial cells across development(BMC, 2021-06-29)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.
-
ItemCanonical PRC2 function is essential for mammary gland development and affects chromatin compaction in mammary organoids(PUBLIC LIBRARY SCIENCE, 2018-08-01)Distinct transcriptional states are maintained through organization of chromatin, resulting from the sum of numerous repressive and active histone modifications, into tightly packaged heterochromatin versus more accessible euchromatin. Polycomb repressive complex 2 (PRC2) is the main mammalian complex responsible for histone 3 lysine 27 trimethylation (H3K27me3) and is integral to chromatin organization. Using in vitro and in vivo studies, we show that deletion of Suz12, a core component of all PRC2 complexes, results in loss of H3K27me3 and H3K27 dimethylation (H3K27me2), completely blocks normal mammary gland development, and profoundly curtails progenitor activity in 3D organoid cultures. Through the application of mammary organoids to bypass the severe phenotype associated with Suz12 loss in vivo, we have explored gene expression and chromatin structure in wild-type and Suz12-deleted basal-derived organoids. Analysis of organoids led to the identification of lineage-specific changes in gene expression and chromatin structure, inferring cell type-specific PRC2-mediated gene silencing of the chromatin state. These expression changes were accompanied by cell cycle arrest but not lineage infidelity. Together, these data indicate that canonical PRC2 function is essential for development of the mammary gland through the repression of alternate transcription programs and maintenance of chromatin states.
-
ItemConstruction of developmental lineage relationships in the mouse mammary gland by single-cell RNA profiling(NATURE PUBLISHING GROUP, 2017-11-20)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.
-
ItemTranscriptome analyses of mouse and human mammary cell subpopulations reveal multiple conserved genes and pathways(BMC, 2010-01-01)INTRODUCTION: Molecular characterization of the normal epithelial cell types that reside in the mammary gland is an important step toward understanding pathways that regulate self-renewal, lineage commitment, and differentiation along the hierarchy. Here we determined the gene expression signatures of four distinct subpopulations isolated from the mouse mammary gland. The epithelial cell signatures were used to interrogate mouse models of mammary tumorigenesis and to compare with their normal human counterpart subsets to identify conserved genes and networks. METHODS: RNA was prepared from freshly sorted mouse mammary cell subpopulations (mammary stem cell (MaSC)-enriched, committed luminal progenitor, mature luminal and stromal cell) and used for gene expression profiling analysis on the Illumina platform. Gene signatures were derived and compared with those previously reported for the analogous normal human mammary cell subpopulations. The mouse and human epithelial subset signatures were then subjected to Ingenuity Pathway Analysis (IPA) to identify conserved pathways. RESULTS: The four mouse mammary cell subpopulations exhibited distinct gene signatures. Comparison of these signatures with the molecular profiles of different mouse models of mammary tumorigenesis revealed that tumors arising in MMTV-Wnt-1 and p53-/- mice were enriched for MaSC-subset genes, whereas the gene profiles of MMTV-Neu and MMTV-PyMT tumors were most concordant with the luminal progenitor cell signature. Comparison of the mouse mammary epithelial cell signatures with their human counterparts revealed substantial conservation of genes, whereas IPA highlighted a number of conserved pathways in the three epithelial subsets. CONCLUSIONS: The conservation of genes and pathways across species further validates the use of the mouse as a model to study mammary gland development and highlights pathways that are likely to govern cell-fate decisions and differentiation. It is noteworthy that many of the conserved genes in the MaSC population have been considered as epithelial-mesenchymal transition (EMT) signature genes. Therefore, the expression of these genes in tumor cells may reflect basal epithelial cell characteristics and not necessarily cells that have undergone an EMT. Comparative analyses of normal mouse epithelial subsets with murine tumor models have implicated distinct cell types in contributing to tumorigenesis in the different models.
-
ItemDeaf-I regulates epithelial cell proliferation and side-branching in the mammary gland(BMC, 2008-10-01)BACKGROUND: The transcription factor DEAF-1 has been identified as a high affinity binding partner of the LIM-only protein LMO4 that plays important roles in mammary gland development and breast cancer. Here we investigated the influence of DEAF-1 on human and mouse mammary epithelial cells both in vitro and in vivo and identified a potential target gene. RESULTS: Overexpression of DEAF-1 in human breast epithelial MCF10A cells enhanced cell proliferation in the mammary acini that develop in 3D cultures. To investigate the effects of Deaf-1 on mammary gland development and oncogenesis, we generated MMTV-Deaf-1 transgenic mice. Increased ductal side-branching was observed in young virgin mammary glands, accompanied by augmented cell proliferation. In addition, the ratio of the progesterone receptor isoforms PRA and PRB, previously implicated in regulating ductal side-branching, was altered. Affymetrix gene profiling studies revealed Rac3 as a potential target gene and quantitative RT-PCR analysis confirmed that Rac3 was upregulated by Deaf-1 in immortalized mouse mammary epithelial cells. Furthermore, MMTV-Deaf-1 transgenic mammary glands were found to have elevated levels of Rac3 mRNA, suggesting that it is a bona fide target. CONCLUSION: We have demonstrated that overexpression of Deaf-1 enhances the proliferation of human breast epithelial cells in vitro and mouse epithelial cells in vivo. Transgenic mammary glands overexpressing Deaf-1 exhibited a modest side-branching phenotype, accompanied by an increase in the number of BrdU-positive cells and a decrease in the proportion of PRA-expressing cells. Although proliferation was enhanced in Deaf-1 transgenic mice, overexpression of this gene was not sufficient to induce the formation of mammary tumors. In addition, our studies identified Rac3, encoding a small Rho-like GTPase, as a potential target of Deaf-1 in mouse mammary epithelial cells.
-
ItemIntegration of microRNA signatures of distinct mammary epithelial cell types with their gene expression and epigenetic portraits(BMC, 2015-06-18)INTRODUCTION: MicroRNAs (miRNAs) have been implicated in governing lineage specification and differentiation in multiple organs; however, little is known about their specific roles in mammopoiesis. We have determined the global miRNA expression profiles of functionally distinct epithelial subpopulations in mouse and human mammary tissue, and compared these to their cognate transcriptomes and epigenomes. Finally, the human miRNA signatures were used to interrogate the different subtypes of breast cancer, with a view to determining miRNA networks deregulated during oncogenesis. METHODS: RNA from sorted mouse and human mammary cell subpopulations was subjected to miRNA expression analysis using the TaqMan MicroRNA Array. Differentially expressed (DE) miRNAs were correlated with gene expression and histone methylation profiles. Analysis of miRNA signatures of the intrinsic subtypes of breast cancer in The Cancer Genome Atlas (TCGA) database versus those of normal human epithelial subpopulations was performed. RESULTS: Unique miRNA signatures characterized each subset (mammary stem cell (MaSC)/basal, luminal progenitor, mature luminal, stromal), with a high degree of conservation across species. Comparison of miRNA and transcriptome profiles for the epithelial subtypes revealed an inverse relationship and pinpointed key developmental genes. Interestingly, expression of the primate-specific miRNA cluster (19q13.4) was found to be restricted to the MaSC/basal subset. Comparative analysis of miRNA signatures with H3 lysine modification maps of the different epithelial subsets revealed a tight correlation between active or repressive marks for the top DE miRNAs, including derepression of miRNAs in Ezh2-deficient cellular subsets. Interrogation of TCGA-identified miRNA profiles with the miRNA signatures of different human subsets revealed specific relationships. CONCLUSIONS: The derivation of global miRNA expression profiles for the different mammary subpopulations provides a comprehensive resource for understanding the interplay between miRNA networks and target gene expression. These data have highlighted lineage-specific miRNAs and potential miRNA-mRNA networks, some of which are disrupted in neoplasia. Furthermore, our findings suggest that key developmental miRNAs are regulated by global changes in histone modification, thus linking the mammary epigenome with genome-wide changes in the expression of genes and miRNAs. Comparative miRNA signature analyses between normal breast epithelial cells and breast tumors confirmed an important linkage between luminal progenitor cells and basal-like tumors.
-
ItemA pooled shRNA screen for regulators of primary mammary stem and progenitor cells identifies roles for Asap1 and Prox1(BMC, 2015-04-03)BACKGROUND: The molecular regulators that orchestrate stem cell renewal, proliferation and differentiation along the mammary epithelial hierarchy remain poorly understood. Here we have performed a large-scale pooled RNAi screen in primary mouse mammary stem cell (MaSC)-enriched basal cells using 1295 shRNAs against genes principally involved in transcriptional regulation. METHODS: MaSC-enriched basal cells transduced with lentivirus pools carrying shRNAs were maintained as non-adherent mammospheres, a system known to support stem and progenitor cells. Integrated shRNAs that altered culture kinetics were identified by next generation sequencing as relative frequency changes over time. RNA-seq-based expression profiling coupled with in vitro progenitor and in vivo transplantation assays was used to confirm a role for candidate genes in mammary stem and/or progenitor cells. RESULTS: Utilizing a mammosphere-based assay, the screen identified several candidate regulators. Although some genes had been previously implicated in mammary gland development, the vast majority of genes uncovered have no known function within the mammary gland. RNA-seq analysis of freshly purified primary mammary epithelial populations and short-term cultured mammospheres was used to confirm the expression of candidate regulators. Two genes, Asap1 and Prox1, respectively implicated in breast cancer metastasis and progenitor cell function in other systems, were selected for further analysis as their roles in the normal mammary gland were unknown. Both Prox1 and Asap1 were shown to act as negative regulators of progenitor activity in vitro, and Asap1 knock-down led to a marked increase in repopulating activity in vivo, implying a role in stem cell activity. CONCLUSIONS: This study has revealed a number of novel genes that influence the activity or survival of mammary stem and/or progenitor cells. Amongst these, we demonstrate that Prox1 and Asap1 behave as negative regulators of mammary stem/progenitor function. Both of these genes have also been implicated in oncogenesis. Our findings provide proof of principle for the use of short-term cultured primary MaSC/basal cells in functional RNAi screens.
-
ItemDifferential methylation analysis of reduced representation bisulfite sequencing experiments using edgeR.(F1000 Research Ltd, 2017)Cytosine methylation is an important DNA epigenetic modification. In vertebrates, methylation occurs at CpG sites, which are dinucleotides where a cytosine is immediately followed by a guanine in the DNA sequence from 5' to 3'. When located in the promoter region of a gene, DNA methylation is often associated with transcriptional silencing of the gene. Aberrant DNA methylation is associated with the development of various diseases such as cancer. Bisulfite sequencing (BS-seq) is the current "gold-standard" technology for high-resolution profiling of DNA methylation. Reduced representation bisulfite sequencing (RRBS) is an efficient form of BS-seq that targets CpG-rich DNA regions in order to save sequencing costs. A typical bioinformatics aim is to identify CpGs that are differentially methylated (DM) between experimental conditions. This workflow demonstrates that differential methylation analysis of RRBS data can be conducted using software and methodology originally developed for RNA-seq data. The RNA-seq pipeline is adapted to methylation by adding extra columns to the design matrix to account for read coverage at each CpG, after which the RRBS and RNA-seq pipelines are almost identical. This approach is statistically natural and gives analysts access to a rich collection of analysis tools including generalized linear models, gene set testing and pathway analysis. The article presents a complete start to finish case study analysis of RRBS profiles of different cell populations from the mouse mammary gland using the Bioconductor package edgeR. We show that lineage-committed cells are typically hyper-methylated compared to progenitor cells and this is true on all the autosomes but not the sex chromosomes. We demonstrate a strong negative correlation between methylation of promoter regions and gene expression as measured by RNA-seq for the same cell types, showing that methylation is a regulatory mechanism involved in epithelial linear commitment.