Gene expression in murine mammary epithelial stem cell-like cells shows similarities to human breast cancer gene expression

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.

Single-Cell Transcription Mapping of Murine and Human Mammary Organoids Responses to Female Hormones

During female adolescence and pregnancy, rising levels of hormones result in a cyclic source of signals that control the development of mammary tissue. While such alterations are well understood from a whole-gland perspective, the alterations that such hormones bring to organoid cultures derived from mammary glands have yet to be fully mapped. This is of special importance given that organoids are considered suitable systems to understand cross species breast development. Here we utilized single-cell transcriptional profiling to delineate responses of murine and human normal breast organoid systems to female hormones across evolutionary distinct species. Collectively, our study represents a molecular atlas of epithelial dynamics in response to estrogen and pregnancy hormones.

PAPAS promotes differentiation of mammary epithelial cells and suppresses breast carcinogenesis

Extensive remodeling of the female mammary epithelium during development and pregnancy has been linked to cancer susceptibility. The faithful response of mammary epithelial cells (MECs) to hormone signaling is key to avoiding breast cancer development. Here, we show that lactogenic differentiation of murine MECs requires silencing of genes encoding ribosomal RNA (rRNA) by the antisense transcript PAPAS. Accordingly, knockdown of PAPAS derepresses rRNA genes, attenuates the response to lactogenic hormones, and induces malignant transformation. Restoring PAPAS levels in breast cancer cells reduces tumorigenicity and lung invasion and activates many interferon-regulated genes previously linked to metastasis suppression. Mechanistically, PAPAS transcription depends on R-loop formation at the 3' end of rRNA genes, which is repressed by RNase H1 and replication protein A (RPA) overexpression in breast cancer cells. Depletion of PAPAS and upregulation of RNase H1 and RPA in human breast cancer underpin the clinical relevance of our findings.

Inhibiting SETD7 methyl-transferase activity impairs differentiation, lipid metabolism and lactogenesis in mammary epithelial cells

SETD7 (SET7/9, KMT7) is a lysine methyltransferase that targets master regulators of cell proliferation and differentiation. Here, the impact of inhibiting SETD7 catalytic activity on mammary epithelial cell differentiation was studied by focusing on genes associated with epithelial differentiation, lactogenesis, and lipid metabolism in HC11 and EpH4 cell lines. Setd7 mRNA and protein levels were induced upon lactogenic differentiation in both cell lines. Inhibition of SETD7 activity by the compound (R)-PFI-2 increased cell proliferation and downregulated E-cadherin, beta-catenin, lactoferrin, insulin-like growth factor binding protein 5, and beta-casein levels. In addition, inhibition of SETD7 activity affected the lipid profile and altered the mRNA expression of the phospholipid biosynthesis-related genes choline phosphotransferase 1, and ethanolamine-phosphate cytidylyltransferase. Altogether, the results suggest that inhibiting SETD7 catalytic activity impairs mammary epithelial and lactogenic differentiation.

Single-cell transcription mapping of murine and human mammary organoids responses to female hormones

During female adolescence and pregnancy, rising levels of hormones result in a cyclic source of signals that control the development of mammary tissue. While such alterations are well understood from a whole-gland perspective, the alterations that such hormones bring to organoid cultures derived from mammary glands have yet to be fully mapped. This is of special importance given that organoids are considered suitable systems to understand cross species breast development. Here we utilized single-cell transcriptional profiling to delineate responses of murine and human normal breast organoid systems to female hormones across evolutionary distinct species. Collectively, our study represents a molecular atlas of epithelial dynamics in response to estrogen and pregnancy hormones.

Global Down-regulation of Gene Expression Induced by Mouse Mammary Tumor Virus (MMTV) in Normal Mammary Epithelial Cells

Mouse mammary tumor virus (MMTV) is a betaretrovirus that causes breast cancer in mice. The mouse mammary epithelial cells are the most permissive cells for MMTV, expressing the highest levels of virus upon infection and being the ones later transformed by the virus due to repeated rounds of infection/superinfection and integration, leading eventually to mammary tumors. The aim of this study was to identify genes and molecular pathways dysregulated by MMTV expression in mammary epithelial cells. Towards this end, mRNAseq was performed on normal mouse mammary epithelial cells stably expressing MMTV, and expression of host genes was analyzed compared with cells in its absence. The identified differentially expressed genes (DEGs) were grouped on the basis of gene ontology and relevant molecular pathways. Bioinformatics analysis identified 12 hub genes, of which 4 were up-regulated (Angp2, Ccl2, Icam, and Myc) and 8 were down-regulated (Acta2, Cd34, Col1a1, Col1a2, Cxcl12, Eln, Igf1, and Itgam) upon MMTV expression. Further screening of these DEGs showed their involvement in many diseases, especially in breast cancer progression when compared with available data. Gene Set Enrichment Analysis (GSEA) identified 31 molecular pathways dysregulated upon MMTV expression, amongst which the PI3-AKT-mTOR was observed to be the central pathway down-regulated by MMTV. Many of the DEGs and 6 of the 12 hub genes identified in this study showed expression profile similar to that observed in the PyMT mouse model of breast cancer, especially during tumor progression. Interestingly, a global down-regulation of gene expression was observed, where nearly 74% of the DEGs in HC11 cells were repressed by MMTV expression, an observation similar to what was observed in the PyMT mouse model during tumor progression, from hyperplasia to adenoma to early and late carcinomas. Comparison of our results with the Wnt1 mouse model revealed further insights into how MMTV expression could lead to activation of the Wnt1 pathway independent of insertional mutagenesis. Thus, the key pathways, DEGs, and hub genes identified in this study can provide important clues to elucidate the molecular mechanisms involved in MMTV replication, escape from cellular anti-viral response, and potential to cause cell transformation. These data also validate the use of the MMTV-infected HC11 cells as an important model to study early transcriptional changes that could lead to mammary cell transformation.

Overexpression of Estrogen Receptor α in Mammary Glands of Aging Mice Is Associated with a Proliferative Risk Signature and Generation of Estrogen Receptor α-Positive Mammary Adenocarcinomas

Age is a risk factor for human estrogen receptor-positive breast cancer, with highest prevalence following menopause. While transcriptome risk profiling is available for human breast cancers, it is not yet developed for prognostication for primary or secondary breast cancer development utilizing at-risk breast tissue. Both estrogen receptor α (ER) and aromatase overexpression have been linked to human breast cancer. Herein, conditional genetically engineered mouse models of estrogen receptor 1 (Esr1) and cytochrome P450 family 19 subfamily A member 1 (CYP19A1) were used to show that induction of Esr1 overexpression just before or with reproductive senescence and maintained through age 30 months resulted in significantly higher prevalence of estrogen receptor-positive adenocarcinomas than CYP19A1 overexpression. All adenocarcinomas tested showed high percentages of ER+ cells. Mammary cancer development was preceded by a persistent proliferative transcriptome risk signature initiated within 1 week of transgene induction that showed parallels to the Prosigna/Prediction Analysis of Microarray 50 human prognostic signature for early-stage human ER+ breast cancer. CYP19A1 mice also developed ER+ mammary cancers, but histology was more divided between adenocarcinoma and adenosquamous, with one ER- adenocarcinoma. Results demonstrate that, like humans, generation of ER+ adenocarcinoma in mice was facilitated by aging mice past the age of reproductive senescence. Esr1 overexpression was associated with a proliferative estrogen pathway-linked signature that preceded appearance of ER+ mammary adenocarcinomas.

Establishment and characterization of equine mammary organoids using a method translatable to other non-traditional model species

Mammary organoid (MaO) models are only available for a few traditional model organisms, limiting our ability to investigate mammary gland development and cancer across mammals. This study established equine mammary organoids (EqMaOs) from cryopreserved mammary tissue, in which mammary tissue fragments were isolated and embedded into a 3D matrix to produce EqMaOs. We evaluated viability, proliferation and budding capacity of EqMaOs at different time points during culture, showing that although the number of proliferative cells decreased over time, viability was maintained and budding increased. We further characterized EqMaOs based on expression of stem cell, myoepithelial and luminal markers, and found that EqMaOs expressed these markers throughout culture and that a bilayered structure as seen in vivo was recapitulated. We used the milk-stimulating hormone prolactin to induce milk production, which was verified by the upregulation of milk proteins, most notably β-casein. Additionally, we showed that our method is also applicable to additional non-traditional mammalian species, particularly domesticated animals such as cats, pigs and rabbits. Collectively, MaO models across species will be a useful tool for comparative developmental and cancer studies.

Aberrant RET expression impacts on normal mammary gland post-lactation transition enhancing cancer potential

RET is a receptor tyrosine kinase with oncogenic potential in the mammary epithelium. Several receptors with oncogenic activity in the breast are known to participate in specific developmental stages. We found that RET is differentially expressed during mouse mammary gland development: RET is present in lactation and its expression dramatically decreases in involution, the period during which the lactating gland returns to a quiescent state after weaning. Based on epidemiological and pre-clinical findings, involution has been described as tumor promoting. Using the Ret/MTB doxycycline-inducible mouse transgenic system, we show that sustained expression of RET in the mammary epithelium during the post-lactation transition to involution is accompanied by alterations in tissue remodeling and an enhancement of cancer potential. Following constitutive Ret expression, we observed a significant increase in neoplastic lesions in the post-involuting versus the virgin mammary gland. Furthermore, we show that abnormal RET overexpression during lactation promotes factors that prime involution, including premature activation of Stat3 signaling and, using RNA sequencing, an acute-phase inflammatory signature. Our results demonstrate that RET overexpression negatively affects the normal post-lactation transition.

Summary: We show that RET activation stimulates Stat3 signaling in mammary epithelial cell culture and in vivo during post-lactation transition, demonstrating that the RET receptor participates in the post-lactation transition priming tumorigenesis.

Tinzaparin inhibits VL30 retrotransposition induced by oxidative stress and/or VEGF in HC11 mouse progenitor mammary cells: Association between inhibition of cancer stem cell proliferation and mammosphere disaggregation

Tinzaparin is an anticoagulant and antiangiogenic drug with inhibitory properties against tumor growth. VEGF stimulates angiogenesis, while an association between reactive oxygen species (ROS) and angiogenesis is involved in tumor progression. The present study aimed to investigate the effect of tinzaparin on VL30 retrotransposition-positive mouse HC11 mammary stem-like epithelial cells, previously reported to be associated with induced mammosphere/cancer stem cell (CSC) generation and tumorigenesis. Under 24 h serum starvation, 15.2% nominal retrotransposition frequency was increased to 29%. Additionally, while treatment with 3–12 ng/ml VEGF further induced retrotransposition frequency in a dose-dependent manner (up to 40.3%), pre-incubation with tinzaparin (2 IU/ml) for 0.5–4 h reduced this frequency to 18.3% in a time-dependent manner, confirmed by analogous results in NIH3T3 fibroblasts. Treatment with 10–40 pg/ml glucose oxidase (GO) for 24 h induced HC11 cell retrotransposition in a dose-dependent manner (up to 82.5%), while a 3 h pre-incubation with tinzaparin (1 or 2 IU/ml) elicited a 13.5 or 25.5% reduction in retrotransposition, respectively. Regarding tumorigenic VL30 retrotransposition-positive HC11 cells, treatment with 2 IU/ml tinzaparin for 5 days reduced proliferation rate in a time-dependent manner (up to ~55%), and after 3 weeks, disaggregated soft agar-formed foci, as well as low-adherent mammospheres, producing single mesenchymal-like cells with a ~50% reduced retrotransposition. With respect to the VL30 retrotransposition mechanism: While 12 ng/ml VEGF increased the level of VL30 and endogenous reverse transcriptase (enRT) transcripts ~1.41- and ~1.16-fold, respectively, subsequent tinzaparin treatment reduced both endogenous/ROS- and VEGF-induced levels 1.15- and 0.40-fold (VL30) and 0.60- and 0.52-fold (enRT), respectively. To the best of our knowledge, these data demonstrate for the first time, the novel inhibition activity of tinzaparin against ROS- and VEGF-induced VL30 retrotransposition, and the proliferation and/or aggregation of mouse HC11 mammosphere/tumor-initiating CSCs, thus contributing to the inhibition of VL30 retrotransposition-induced primary tumor growth.

Remodeling of Zn2+ homeostasis upon differentiation of mammary epithelial cells

Zinc is the second most abundant transition metal in humans and an essential nutrient required for growth and development of newborns. During lactation, mammary epithelial cells differentiate into a secretory phenotype, uptake zinc from blood circulation, and export it into mother's milk. At the cellular level, many zinc-dependent cellular processes, such as transcription, metabolism of nutrients, and proliferation are involved in the differentiation of mammary epithelial cells. Using mouse mammary epithelial cells as a model system, we investigated the remodeling of zinc homeostasis during differentiation induced by treatment with the lactogenic hormones cortisol and prolactin. RNA-Seq at different stages of differentiation revealed changes in global gene expression, including genes encoding zinc-dependent proteins and regulators of zinc homeostasis. Increases in mRNA levels of three zinc homeostasis genes, Slc39a14 (ZIP14) and metallothioneins (MTs) I and II were induced by cortisol but not by prolactin. The cortisol-induced increase was partially mediated by the nuclear glucocorticoid receptor signaling pathway. An increase in the cytosolic labile Zn2+ pool was also detected in lactating mammary cells, consistent with upregulation of MTs. We found that the zinc transporter ZIP14 was important for the expression of a major milk protein, whey acid protein (WAP), as knockdown of ZIP14 dramatically decreased WAP mRNA levels. In summary, our study demonstrated remodeling of zinc homeostasis upon differentiation of mammary epithelial cells resulting in changes in cytosolic Zn2+ and differential expression of zinc homeostasis genes, and these changes are important for establishing the lactation phenotype.

External and internal EGFR-activating signals drive mammary epithelial cells proliferation and viability

During puberty, the mammary gland undergoes an intense growth, dependent on the interplay between the Epidermal Growth Factor Receptor (EGFR) in the stroma and different mammary epithelial receptors. We hypothesize that EGFR expressed in the mammary epithelium also has a role in puberty and the epithelial cells can self-sustain by EGFR-mediated autocrine signaling. We adopted mammary cell lines from different species, as in vitro model for the epithelium, and we observed that EGFR-signaling positively affects their survival and proliferation. Once deprived of external growth factors, mammary cells still showed strong Erk 1/2 phosphorylation, abolished upon EGFR inhibition, coupled with a further reduction in survival and proliferation. Based on gene expression analysis, three EGFR-ligands (AREG, EREG and HBEGF) are likely to mediate this autocrine signaling. In conclusion, internal EGFR-activating signals sustain mammary epithelial cell proliferation and survival in vitro.

VL30 retrotransposition is associated with induced EMT, CSC generation and tumorigenesis in HC11 mouse mammary stem‑like epithelial cells

Retrotransposons copy their sequences via an RNA intermediate, followed by reverse transcription into cDNA and random insertion, into a new genomic locus. New retrotransposon copies may lead to cell transformation and/or tumorigenesis through insertional mutagenesis. Methylation is a major defense mechanism against retrotransposon RNA expression and retrotransposition in differentiated cells, whereas stem cells are relatively hypo-methylated. Epithelial-to-mesenchymal transition (EMT), which transforms normal epithelial cells into mesenchymal-like cells, also contributes to tumor progression and tumor metastasis. Cancer stem cells (CSCs), a fraction of undifferentiated tumor-initiating cancer cells, are reciprocally related to EMT. In the present study, the outcome of long terminal repeat (LTR)-Viral-Like 30 (VL30) retrotransposition was examined in mouse mammary stem-like/progenitor HC11 epithelial cells. The transfection of HC11 cells with a VL30 retrotransposon, engineered with an EGFP-based retrotransposition cassette, elicited a higher retrotransposition frequency in comparison to differentiated J3B1A and C127 mouse mammary cells. Fluorescence microscopy and PCR analysis confirmed the specificity of retrotransposition events. The differentiated retrotransposition-positive cells retained their epithelial morphology, while the respective HC11 cells acquired mesenchymal features associated with the loss of E-cadherin, the induction of N-cadherin, and fibronectin and vimentin protein expression, as well as an increased transforming growth factor (TGF)-β1, Slug, Snail-1 and Twist mRNA expression. In addition, they were characterized by cell proliferation in low serum, and the acquisition of CSC-like properties indicated by mammosphere formation under anchorage-independent conditions. Mammospheres exhibited an increased Nanog and Oct4 mRNA expression and a CD44⁺/CD24^−/low antigenic phenotype, as well as self-renewal and differentiation capacity, forming mammary acini-like structures. DNA sequencing analysis of retrotransposition-positive HC11 cells revealed retrotransposed VL30 copies integrated at the vicinity of EMT-, cancer type- and breast cancer-related genes. The inoculation of these cells into Balb/c mice produced cytokeratin-positive tumors containing pancytokeratin-positive cells, indicative of cell invasion features. On the whole, the findings of the present study demonstrate, for the first time, to the best of our knowledge, that stem-like epithelial HC11 cells are amenable to VL30 retrotransposition associated with the induction of EMT and CSC generation, leading to tumorigenesis.

A Primary Cell and Organoid Platform for Evaluating Pharmacological Responses in Mammary Epithelial Cells

An essential process in predicting the in vivo pharmacological activity of a candidate molecule involves the evaluation of target responses using established model systems. While these models largely comprise immortalized cells, which are often serially passaged as monolayers on uniformly stiff substrates and are modified to overexpress one or more components of the pathway-of-interest, the importance of cell identity, heterogeneity, and three-dimensional (3D) context to target response is gaining increasing attention. Here, we assess intracellular calcium responses in mouse mammary epithelial cells in three distinct model systems: 3D primary organoids, 2D primary epithelial cells, and 2D immortalized cells. Specifically, we assess intracellular calcium responses to a number of extracellular signals implicated in the regulation of basal (or myoepithelial) cell function. These findings provide further insights into cell type and context-specific pharmacological responses in mammary epithelial cells and highlight the opportunities and challenges in the adoption of architecturally complex and heterogeneous in vitro assays in pharmacological research.

Replication-dependent histone isoforms: a new source of complexity in chromatin structure and function

Replication-dependent histones are expressed in a cell cycle regulated manner and supply the histones necessary to support DNA replication. In mammals, the replication-dependent histones are encoded by a family of genes that are located in several clusters. In humans, these include 16 genes for histone H2A, 22 genes for histone H2B, 14 genes for histone H3, 14 genes for histone H4 and 6 genes for histone H1. While the proteins encoded by these genes are highly similar, they are not identical. For many years, these genes were thought to encode functionally equivalent histone proteins. However, several lines of evidence have emerged that suggest that the replication-dependent histone genes can have specific functions and may constitute a novel layer of chromatin regulation. This Survey and Summary reviews the literature on replication-dependent histone isoforms and discusses potential mechanisms by which the small variations in primary sequence between the isoforms can alter chromatin function. In addition, we summarize the wealth of data implicating altered regulation of histone isoform expression in cancer.

A miR-206 regulated gene landscape enhances mammary epithelial differentiation

miR‐206 is known to suppress breast cancer. However, while it is expressed in mammary stem cells, its function in such nontumor cells is not well understood. Here, we explore the role of miR‐206 in undifferentiated, stem‐like mammary cells using the murine mammary differentiation model HC11, genome‐wide gene expression analysis, and functional assays. We describe the miR‐206‐regulated gene landscape and propose a network whereby miR‐206 suppresses tumor development. We functionally demonstrate that miR‐206 in nontumor stem‐like cells induces a G1–S cell cycle arrest, and reduces colony formation and epithelial‐to‐mesenchymal transition markers. Finally, we show that addition of miR‐206 accelerates the mammary differentiation process along with related accumulation of lipids. We conclude that miR‐206 impacts a network of signaling pathways, and acts as a regulator of proliferation, stemness, and mammary cell differentiation in nontumor stem‐like mammary cells. Our study provides a broad insight into the breast cancer suppressive functions of miR‐206.

Comprehensive profiling of transcriptional networks specific for lactogenic differentiation of HC11 mammary epithelial stem-like cells

The development of mammary gland as a lactogenic tissue is a highly coordinated multistep process. The epithelial cells of lactiferous tubules undergo profound changes during the developmental window of puberty, pregnancy, and lactation. Several hormones including estrogen, progesterone, glucocorticoids and prolactin act in concert, and orchestrate the development of mammary gland. Understanding the gene regulatory networks that coordinate proliferation and differentiation of HC11 Mammary Epithelial stem-like Cells (MEC) under the influence of lactogenic hormones is critical for elucidating the mechanism of lactogenesis in detail. In this study, we analyzed transcriptome profiles of undifferentiated MEC (normal) and compared them with Murine Embryonic Stem Cells (ESC) using next-generation mRNA sequencing. Further, we analyzed the transcriptome output during lactogenic differentiation of MEC following treatment with glucocorticoids (primed state) and both glucocorticoids and prolactin together (prolactin state). We established stage-specific gene regulatory networks in ESC and MEC (normal, priming and prolactin states). We validated the top up-and downregulated genes in each stage of differentiation of MEC by RT-PCR and found that they are comparable with that of RNA-seq data. HC11 MEC display decreased expression of Pou5f1 and Sox2, which is crucial for the differentiation of MEC, which otherwise ensure pluripotency to ESC. Cited4 is induced during priming and is involved in milk secretion. MEC upon exposure to both glucocorticoids and prolactin undergo terminal differentiation, which is associated with the expression of several genes, including Xbp1 and Cbp that are required for cell growth and differentiation. Our study also identified differential expression of transcription factors and epigenetic regulators in each stage of lactogenic differentiation. We also analyzed the transcriptome data for the pathways that are selectively activated during lactogenic differentiation. Further, we found that selective expression of chromatin modulators (Dnmt3l, Chd9) in response to glucocorticoids suggests a highly coordinated stage-specific lactogenic differentiation of MEC.

Effects of miRNAs on functions of breast cancer stem cells and treatment of breast cancer

Breast cancer is one of the most common malignancies for women, which accounts for 30% of all female malignancies. The formation of breast cancer stem cells (BCSCs) is attributed to the acquisition of stemness of tumor cells. With self-renewal potential, these stem cells are insensitive to either radiotherapy or chemotherapy but are significant in regulating tumor behaviors and drug resistance. MicroRNA (miRNA) is a kind of noncoding small RNA for negatively regulating gene expressions. Research findings suggest that many miRNAs specifically regulate the expression of target genes and signal pathways of BCSCs. They play an important role in self-renewal, growth, and metastasis of breast cancer cells as potential targets for treating breast cancer. These signal pathways include phosphatase and tensin homolog deleted on chromosome 10-phosphatidylinositol 3-kinase/Akt, Wnt/β-catenin, Notch, and so on. This paper reviews the progress of research about miRNAs in self-renewal, metastasis, epithelial-mesenchymal transition and metastasis, mediation of resistance to chemotherapies, and treatment of breast cancer.

R-spondin3 Is Associated with Basal-Progenitor Behavior in Normal and Tumor Mammary Cells

R-spondin3 (RSPO3) is a member of a family of secreted proteins that enhance Wnt signaling pathways in diverse processes, including cancer. However, the role of RSPO3 in mammary gland and breast cancer development remains unclear. In this study, we show that RSPO3 is expressed in the basal stem cell-enriched compartment of normal mouse mammary glands but is absent from committed mature luminal cells in which exogenous RSPO3 impairs lactogenic differentiation. RSPO3 knockdown in basal-like mouse mammary tumor cells reduced canonical Wnt signaling, epithelial-to-mesenchymal transition-like features, migration capacity, and tumor formation in vivo Conversely, RSPO3 overexpression, which was associated with some LGR and RUNX factors, highly correlated with the basal-like subtype among patients with breast cancer. Thus, we identified RSPO3 as a novel key modulator of breast cancer development and a potential target for treatment of basal-like breast cancers.Significance: These findings identify RSPO3 as a potential therapetuic target in basal-like breast cancers.Graphical Abstract: http://cancerres.aacrjournals.org/content/canres/78/16/4497/F1.large.jpg Cancer Res; 78(16); 4497-511. ©2018 AACR.

RNA sequencing of murine mammary epithelial stem-like cells (HC11) undergoing lactogenic differentiation and its comparison with embryonic stem cells

OBJECTIVES

Understanding of transcriptional networks specifying HC11 murine mammary epithelial stem cell-like cells (MEC) in comparison with embryonic stem cells (ESCs) and their rewiring, under the influence of glucocorticoids (GC) and prolactin (PRL) hormones, is critical for elucidating the mechanism of lactogenesis. In this data note, we provide RNA sequencing data from murine MECs and ESCs, MECs treated with steroid hormone alone and in combination with PRL. This data could help in understanding temporal dynamics of mRNA transcription that impact the process of lactogenesis associated with mammary gland development. Further integration of these data sets with existing datasets of cells derived from various stages of mammary gland development and different types of breast tumors, should pave the way for effective prognosis and to develop therapies for breast cancer.

DATA DESCRIPTION

We have generated RNA-sequencing data representing steady-state levels of mRNAs from murine ESCs, normal MECs (N), MECs primed (P) with hydrocortisone (HC) alone and in combination with PRL hormone by using Illumina sequencing platform. We have generated ~ 58 million reads for ESCs with an average length of ~ 100 nt and an average 115 million good quality mapped reads with an average length of ~ 150 nt for different stages of MECs differentiation.

Basolateral pressure challenges mammary epithelial cell monolayer integrity, in vitro

Mammary gland epithelium is physiologically exposed to variations of hydrostatic pressure due to accumulation of milk and removal by suckling and mechanical milking. Integrity of the mammary gland epithelium primarily relies on the tight junction. To analyze pressure-induced effects on the tight junction, we established a modified Ussing chamber and tested the hypothesis if hydrostatic pressure on the basal side of the epithelium is able to affect barrier properties in a mammary epithelial cell model, in vitro. Therefore, a conventional Ussing chamber was modified by an additional tube system to apply hydrostatic pressure. Monolayers of the mammary epithelial cell line HC11 were mounted in the modified Ussing chambers and incubated with increasing basal hydrostatic pressure. Transepithelial resistance and short circuit current were recorded and compared to controls. Hydrostatic pressure was stably applied and incubation steps of 30 min were technically feasible, leading to a decrease of transepithelial resistance and an increase of short circuit current in all monolayers. In a series of experiments simulating the physiological exposure time by short intervals of 5 min, these electrophysiological findings were also observed, and monolayer integrity was not significantly perturbed as analyzed by fluorescence immunohistochemistry selectively staining tight junction proteins. Moreover, electrophysiology demonstrated reversibility of effects. In conclusion, the modified Ussing chamber is an adequate method to analyze the effects of hydrostatic pressure on epithelial cell monolayers, in vitro. Both, the reduction of transepithelial resistance and the increase of short circuit current may be interpreted as protective reactions.

Distance in cancer gene expression from stem cells predicts patient survival

The degree of histologic cellular differentiation of a cancer has been associated with prognosis but is subjectively assessed. We hypothesized that information about tumor differentiation of individual cancers could be derived objectively from cancer gene expression data, and would allow creation of a cancer phylogenetic framework that would correlate with clinical, histologic and molecular characteristics of the cancers, as well as predict prognosis. Here we utilized mRNA expression data from 4,413 patient samples with 7 diverse cancer histologies to explore the utility of ordering samples by their distance in gene expression from that of stem cells. A differentiation baseline was obtained by including expression data of human embryonic stem cells (hESC) and human mesenchymal stem cells (hMSC) for solid tumors, and of hESC and CD34+ cells for liquid tumors. We found that the correlation distance (the degree of similarity) between the gene expression profile of a tumor sample and that of stem cells orients cancers in a clinically coherent fashion. For all histologies analyzed (including carcinomas, sarcomas, and hematologic malignancies), patients with cancers with gene expression patterns most similar to that of stem cells had poorer overall survival. We also found that the genes in all undifferentiated cancers of diverse histologies that were most differentially expressed were associated with up-regulation of specific oncogenes and down-regulation of specific tumor suppressor genes. Thus, a stem cell-oriented phylogeny of cancers allows for the derivation of a novel cancer gene expression signature found in all undifferentiated forms of diverse cancer histologies, that is competitive in predicting overall survival in cancer patients compared to previously published prediction models, and is coherent in that gene expression was associated with up-regulation of specific oncogenes and down-regulation of specific tumor suppressor genes associated with regulation of the multicellular state.

Genome-wide effects of MELK-inhibitor in triple-negative breast cancer cells indicate context-dependent response with p53 as a key determinant

Triple-negative breast cancer (TNBC) is an aggressive, highly recurrent breast cancer subtype, affecting approximately one-fifth of all breast cancer patients. Subpopulations of treatment-resistant cancer stem cells within the tumors are considered to contribute to disease recurrence. A potential druggable target for such cells is the maternal embryonic leucine-zipper kinase (MELK). MELK expression is upregulated in mammary stem cells and in undifferentiated cancers, where it correlates with poor prognosis and potentially mediates treatment resistance. Several MELK inhibitors have been developed, of which one, OTSSP167, is currently in clinical trials. In order to better understand how MELK and its inhibition influence TNBC, we verified its anti-proliferative and apoptotic effects in claudin-low TNBC cell lines MDA-MB-231 and SUM-159 using MTS assays and/or trypan blue viability assays together with analysis of PARP cleavage. Then, using microarrays, we explored which genes were affected by OTSSP167. We demonstrate that different sets of genes are regulated in MDA-MB-231 and SUM-159, but in both cell lines genes involved in cell cycle, mitosis and protein metabolism and folding were regulated. We identified p53 (TP53) as a potential upstream regulator of the regulated genes. Using western blot we found that OTSSP167 downregulates mutant p53 in all tested TNBC cell lines (MDA-MB-231, SUM-159, and BT-549), but upregulates wild-type p53 in the luminal A subtype MCF-7 cell line. We propose that OTSSP167 might have context-dependent or off-target effects, but that one consistent mechanism of action could involve the destabilization of mutant p53.

Runx Genes in Breast Cancer and the Mammary Lineage

A full understanding of RUNX gene function in different epithelial lineages has been thwarted by the lethal phenotypes observed when constitutively knocking out these mammalian genes. However temporal expression of the Runx genes throughout the different phases of mammary gland development is indicative of a functional role in this tissue. A few studies have emerged describing how these genes impact on the fate of mammary epithelial cells by regulating lineage differentiation and stem/progenitor cell potential, with implications for the transformed state. The importance of the RUNX/CBFβ core factor binding complex in breast cancer has very recently been highlighted with both RUNX1 and CBFβ appearing in a comprehensive gene list of predicted breast cancer driver mutations. Nonetheless, the evidence to date shows that the RUNX genes can have dualistic outputs with respect to promoting or constraining breast cancer phenotypes, and that this may be aligned to individual subtypes of the clinical disease. We take this opportunity to review the current literature on RUNX and CBFβ in the normal and neoplastic mammary lineage while appreciating that this is likely to be the tip of the iceberg in our knowledge.

Nanomedicine-Mediated Therapies to Target Breast Cancer Stem Cells

Accumulating evidences have suggested the existence of breast cancer stem cells (BCSCs), which possess the potential of both self-renewal and differentiation. The origin of BCSCs might have relationship to the development of normal mammary stem cells. BCSCs are believed to play a key role in the initiation, recurrence and chemo-/radiotherapy resistances of breast cancer. Therefore, elimination of BCSCs is crucial for breast cancer therapy. However, conventional chemo and radiation therapies cannot eradicate BCSCs effectively. Fortunately, nanotechnology holds great potential for specific and efficient anti-BCSCs treatment. “Smart” nanocarriers can distinguish BCSCs from the other breast cancer cells and selectively deliver therapeutic agents to the BCSCs. Emerging findings suggest that BCSCs in breast cancer could be successfully inhibited and even eradicated by functionalized nanomedicines. In this review, we focus on origin of BCSCs, strategies used to target BCSCs, and summarize the nanotechnology-based delivery systems that have been applied for eliminating BCSCs in breast cancer.

Runx2 contributes to the regenerative potential of the mammary epithelium

Although best known for its role in bone development and associated structures the transcription factor RUNX2 is expressed in a wide range of lineages, including those of the mammary gland. Previous studies have indicated that Runx2 can regulate aspects of mammary cell function and influence the properties of cancer cells. In this study we investigate the role of Runx2 in the mammary stem/progenitor population and its relationship with WNT signalling. Results show that RUNX2 protein is differentially expressed throughout embryonic and adult development of the murine mammary gland with high levels of expression in mammary stem-cell enriched cultures. Importantly, functional analysis reveals a role for Runx2 in mammary stem/progenitor cell function in in vitro and in vivo regenerative assays. Furthermore, RUNX2 appears to be associated with WNT signalling in the mammary epithelium and is specifically upregulated in mouse models of WNT-driven breast cancer. Overall our studies reveal a novel function for Runx2 in regulating mammary epithelial cell regenerative potential, possibly acting as a downstream target of WNT signalling.

miR-200a inhibits migration of triple-negative breast cancer cells through direct repression of the EPHA2 oncogene

Triple-negative breast cancer (TNBC) is characterized by aggressiveness and affects 10-20% of breast cancer patients. Since TNBC lacks expression of ERα, PR and HER2, existing targeted treatments are not effective and the survival is poor. In this study, we demonstrate that the tumor suppressor microRNA miR-200a directly regulates the oncogene EPH receptor A2 (EPHA2) and modulates TNBC migration. We show that EPHA2 expression is correlated with poor survival specifically in basal-like breast cancer and that its expression is repressed by miR-200a through direct interaction with the 3'UTR of EPHA2. This regulation subsequently affects the downstream activation of AMP-activated protein kinase (AMPK) and results in decreased cell migration of TNBC. We establish that miR-200a directs cell migration in a dual manner; in addition to regulating the well-characterized E-cadherin pathway it also regulates a EPHA2 pathway. The miR-200a-EPHA2 axis is a novel mechanism highlighting the possibility of utilizing miR-200a delivery to target TNBC metastases.

miR-206 inhibits cell migration through direct targeting of the actin-binding protein coronin 1C in triple-negative breast cancer

Patients with triple-negative breast cancer (TNBC) have an overall poor prognosis, which is primarily due to a high metastatic capacity of these tumors. Novel therapeutic approaches to target the signaling pathways that promote metastasis are desirable, in order to improve the outcome for these patients. A loss of function of a microRNA, miR-206, is related to increased metastasis potential in breast cancers but the mechanism is not known. In this study, we show that miR-206 was decreased in TNBC clinical tumor samples and cell lines whereas one of its predicted targets, actin-binding protein CORO1C, was increased. Expression of miR-206 significantly reduced proliferation and migration while repressing CORO1C mRNA and protein levels. We demonstrate that miR-206 interacts with the 3'-untranslated region (3'-UTR) of CORO1C and regulates this gene post-transcriptionally. This post-transcriptional regulation was dependent on two miR-206-binding sites within the 3'-UTR of CORO1C and was relieved by mutations of corresponding sites. Further, silencing of CORO1C reduced tumor cell migration and affected the actin skeleton and cell morphology, similar to miR-206 expression, but did not reduce proliferation. In accordance with this, overexpression of CORO1C rescued the inhibitory effect of miR-206 on cell migration. Our findings suggest that miR-206 represses tumor cell migration through direct targeting of CORO1C in TNBC cells which modulates the actin filaments. This pathway is a novel mechanism that offers a mechanistic basis through which the metastatic potential of TNBC tumors could be targeted.

Fatty acid and phospholipid biosynthetic pathways are regulated throughout mammary epithelial cell differentiation and correlate to breast cancer survival

This work combined gene and protein expression, gas chromatography-flame ionization detector, and hydrophilic interaction liquid chromatography-tandem mass spectrometry to compare lipid metabolism changes in undifferentiated/proliferating vs. functionally differentiated mammary epithelial cells (MECs) and to study their correlation to breast cancer survival. Sixty-eight genes involved in lipid metabolism were changed in MEC differentiation. Differentiated cells showed induction of Elovl6 (2-fold), Scd1 (4-fold), and Fads2 (2-fold), which correlated with increased levels of C16:1 n-7 and C18:1 n-9 (1.5-fold), C20:3 n-6 (2.5-fold), and C20:4 n-6 (6-fold) fatty acids (FAs) and more phospholipids (PLs) containing these species. Further, increased expression (2- to 3-fold) of genes in phosphatidylethanolamine (PE) de novo biosynthesis resulted in a 20% PE increase. Proliferating/undifferentiated cells showed higher C16:0 (1.7-fold) and C18:2 n-6 (4.2-fold) levels and more PLs containing C16:0 FAs [PC(16:0/16:1), PG(16:0/18:2), PG(16:0/18:1), and SM(16:0/18:0)]. Kaplan-Meier analysis of data from 3455 patients with breast cancer disclosed a positive correlation for 59% of genes expressed in differentiated MECs with better survival. PE biosynthesis and FA oxidation correlated with better prognosis in patients with breast cancer, including the basal-like subtype. Therefore, genes involved in mammary gland FA and PL metabolism and their resulting molecular species reflect the cellular proliferative ability and differentiation state and deserve further studies as potential markers of breast cancer progression

Mammary morphogenesis and regeneration require the inhibition of EMT at terminal end buds by Ovol2 transcriptional repressor

Epithelial cells possess remarkable plasticity, having the ability to become mesenchymal cells through alterations in adhesion and motility (epithelial-to-mesenchymal transition [EMT]). However, how epithelial plasticity is kept in check in epithelial cells during tissue development and regeneration remains to be fully understood. Here we show that restricting the EMT of mammary epithelial cells by transcription factor Ovol2 is required for proper morphogenesis and regeneration. Deletion of Ovol2 blocks mammary ductal morphogenesis, depletes stem and progenitor cell reservoirs, and leads epithelial cells to undergo EMT in vivo to become nonepithelial cell types. Ovol2 directly represses myriad EMT inducers, and its absence switches response to TGF-β from growth arrest to EMT. Furthermore, forced expression of the repressor isoform of Ovol2 is able to reprogram metastatic breast cancer cells from a mesenchymal to an epithelial state. Our findings underscore the critical importance of exquisitely regulating epithelial plasticity in development and cancer.

SAP domain-dependent Mkl1 signaling stimulates proliferation and cell migration by induction of a distinct gene set indicative of poor prognosis in breast cancer patients

Background

The main cause of death of breast cancer patients is not the primary tumor itself but the metastatic disease. Identifying breast cancer-specific signatures for metastasis and learning more about the nature of the genes involved in the metastatic process would 1) improve our understanding of the mechanisms of cancer progression and 2) reveal new therapeutic targets. Previous studies showed that the transcriptional regulator megakaryoblastic leukemia-1 (Mkl1) induces tenascin-C expression in normal and transformed mammary epithelial cells. Tenascin-C is known to be expressed in metastatic niches, is highly induced in cancer stroma and promotes breast cancer metastasis to the lung.

Methods

Using HC11 mammary epithelial cells overexpressing different Mkl1 constructs, we devised a subtractive transcript profiling screen to identify the mechanism by which Mkl1 induces a gene set co-regulated with tenascin-C. We performed computational analysis of the Mkl1 target genes and used cell biological experiments to confirm the effect of these gene products on cell behavior. To analyze whether this gene set is prognostic of accelerated cancer progression in human patients, we used the bioinformatics tool GOBO that allowed us to investigate a large breast tumor data set linked to patient data.

Results

We discovered a breast cancer-specific set of genes including tenascin-C, which is regulated by Mkl1 in a SAP domain-dependent, serum response factor-independent manner and is strongly implicated in cell proliferation, cell motility and cancer. Downregulation of this set of transcripts by overexpression of Mkl1 lacking the SAP domain inhibited cell growth and cell migration. Many of these genes are direct Mkl1 targets since their promoter-reporter constructs were induced by Mkl1 in a SAP domain-dependent manner. Transcripts, most strongly reduced in the absence of the SAP domain were mechanoresponsive. Finally, expression of this gene set is associated with high-proliferative poor-outcome classes in human breast cancer and a strongly reduced survival rate for patients independent of tumor grade.

Conclusions

This study highlights a crucial role for the transcriptional regulator Mkl1 and its SAP domain during breast cancer progression. We identified a novel gene set that correlates with bad prognosis and thus may help in deciding the rigor of therapy.

MicroRNAs and triple negative breast cancer

Triple Negative Breast Cancer (TNBC) is a very aggressive tumor subtype, which still lacks specific markers for an effective targeted therapy. Despite the common feature of negativity for the three most relevant receptors (ER, PgR and HER2), TNBC is a very heterogeneous disease where different subgroups can be recognized, and both gene and microRNA profiling studies have recently been carried out to dissect the different molecular entities. Moreover, several microRNAs playing a crucial role in triple negative breast cancer biology have been identified, providing the experimental basis for a possible therapeutic application. Indeed, the causal involvement of microRNAs in breast cancer and the possible use of these small noncoding RNA molecules as biomarkers has been extensively studied with promising results. Their application as therapeutic tools might represent an innovative approach, especially for a tumor subgroup still lacking an efficient and specific therapy such as TNBC. In this review, we summarize our knowledge on the most important microRNAs described in TNBC.

Cancer research by means of tissue engineering--is there a rationale?

Tissue engineering (TE) has evoked new hopes for the cure of organ failure and tissue loss by creating functional substitutes in the laboratory. Besides various innovations in the context of Regenerative Medicine (RM), TE also provided new technology platforms to study mechanisms of angiogenesis and tumour cell growth as well as potentially tumour cell spreading in cancer research. Recent advances in stem cell technology--including embryonic and adult stem cells and induced pluripotent stem cells--clearly show the need to better understand all relevant mechanisms to control cell growth when such techniques will be administered to patients. Such TE-Cancer research models allow us to investigate the interactions that occur when replicating physiological and pathological conditions during the initial phases of replication, morphogenesis, differentiation and growth under variable given conditions. Tissue microenvironment has been extensively studied in many areas of TE and it plays a crucial role in cell signalling and regulation of normal and malignant cell functions. This article is intended to give an overview on some of the most recent developments and possible applications of TE and RM methods with regard to the improvement of cancer research with TE platforms. The synthesis of TE with innovative methods of molecular biology and stem-cell technology may help investigate and potentially modulate principal phenomena of tumour growth and spreading, as well as tumour-related angiogenesis. In the future, these models have the potential to investigate the optimal materials, culture conditions and material structure to propagate tumour growth.

RUNX2 in mammary gland development and breast cancer

Runx2 is best known as an essential factor in osteoblast differentiation and bone development but, like many other transcription factors involved in development, is known to operate over a much wider tissue range. Our understanding of these other aspects of Runx2 function is still at a relatively early stage and the importance of its role in cell fate decisions and lineage maintenance in non-osseous tissues is only beginning to emerge. One such tissue is the mammary gland, where Runx2 is known to be expressed and participate in the regulation of mammary specific genes. Furthermore, differential and temporal expression of this gene is observed during mammary epithelial differentiation in vivo, strongly indicative of an important functional role. Although the precise nature of that role remains elusive, preliminary evidence hints at possible involvement in the regulation of mammary stem and/or progenitor cells. As with many genes important in regulating cell fate, RUNX2 has also been linked to metastatic cancer where in some established breast cell lines, retention of expression is associated with a more invasive phenotype. More recently, expression analysis has been extended to primary breast cancers where high levels of RUNX2 align with a specific subtype of the disease. That RUNX2 expression correlates with the so called "Triple Negative" subtype is particularly interesting given the known cross talk between Runx2 and estrogen receptor signaling pathways. This review summaries our current understanding of Runx2 in mammary gland development and cancer, and postulates a role that may link both these processes.

EXP-PAC: providing comparative analysis and storage of next generation gene expression data

Microarrays and more recently RNA sequencing has led to an increase in available gene expression data. How to manage and store this data is becoming a key issue. In response we have developed EXP-PAC, a web based software package for storage, management and analysis of gene expression and sequence data. Unique to this package is SQL based querying of gene expression data sets, distributed normalization of raw gene expression data and analysis of gene expression data across experiments and species. This package has been populated with lactation data in the international milk genomic consortium web portal (http://milkgenomics.org/). Source code is also available which can be hosted on a Windows, Linux or Mac APACHE server connected to a private or public network (http://mamsap.it.deakin.edu.au/~pcc/Release/EXP_PAC.html).

MicroRNA-regulated gene networks during mammary cell differentiation are associated with breast cancer

MicroRNAs (miRNAs) play pivotal roles in stem cell biology, differentiation and oncogenesis and are of high interest as potential breast cancer therapeutics. However, their expression and function during normal mammary differentiation and in breast cancer remain to be elucidated. In order to identify which miRNAs are involved in mammary differentiation, we thoroughly investigated miRNA expression during functional differentiation of undifferentiated, stem cell-like, murine mammary cells using two different large-scale approaches followed by qPCR. Significant changes in expression of 21 miRNAs were observed in repeated rounds of mammary cell differentiation. The majority, including the miR-200 family and known tumor suppressor miRNAs, was upregulated during differentiation. Only four miRNAs, including oncomiR miR-17, were downregulated. Pathway analysis indicated complex interactions between regulated miRNA clusters and major pathways involved in differentiation, proliferation and stem cell maintenance. Comparisons with human breast cancer tumors showed the gene profile from the undifferentiated, stem-like stage clustered with that of poor-prognosis breast cancer. A common nominator in these groups was the E2F pathway, which was overrepresented among genes targeted by the differentiation-induced miRNAs. A subset of miRNAs could further discriminate between human non-cancer and breast cancer cell lines, and miR-200a/miR-200b, miR-146b and miR-148a were specifically downregulated in triple-negative breast cancer cells. We show that miR-200a/miR-200b can inhibit epithelial-mesenchymal transition (EMT)-characteristic morphological changes in undifferentiated, non-tumorigenic mammary cells. Our studies propose EphA2 as a novel and important target gene for miR-200a. In conclusion, we present evidentiary data on how miRNAs are involved in mammary cell differentiation and indicate their related roles in breast cancer.

MYBL2 (B-MYB) in cervical cancer: putative biomarker

INTRODUCTION

Cervical cancer is the second most common cancer affecting women worldwide. It is characterized by chromosomal aberrations and alteration in the expression levels of many cell cycle regulatory proteins. MYBL2 is a member of the MYB proto-oncogene family that encodes DNA binding proteins. These proteins are involved in cell proliferation and control of cellular differentiation.

MATERIALS AND METHODS

Four established cervical cancer cell lines were examined and compared with normal cervix using gene expression profiling and comparative genomic hybridization, and results were correlated to identify potential novel cervical cancer biomarkers. Results were validated using TaqMan polymerase chain reaction, and the potential role of MYBL2 as a clinical biomarker was then evaluated by immunohistochemistry on 30 tissue samples.

RESULTS

MYBL2 was found to be overexpressed in the cervical cancer cell lines by gene expression profiling, and this result was confirmed using TaqMan polymerase chain reaction. Analysis of comparative genomic hybridization data indicated that chromosome 20q13.1, which encodes the MYBL2 gene, was amplified in the human papillomavirus (HPV) type 16-positive CaSki and SiHa cell lines but not in the HPV-18-positive HeLa or HPV-negative C33A cell lines.

DISCUSSION

Although MYBL2 staining was predominantly absent in normal cervical epithelium, strong staining (score of 2 or 3) was identified in all cases of cervical intraepithelial neoplasia, cervical glandular intraepithelial neoplasia, and invasive cancer on immunohistochemistry. In addition, strong staining of a population of diffusely scattered single cells is identified. We postulate that these may represent so-called cancer stem-like cells.

Key signalling nodes in mammary gland development and cancer: Myc

Myc has been intensely studied since its discovery more than 25 years ago. Insight has been gained into Myc's function in normal physiology, where its role appears to be organ specific, and in cancer where many mechanisms contribute to aberrant Myc expression. Numerous signals and pathways converge on Myc, which in turn acts on a continuously growing number of identified targets, via transcriptional and nontranscriptional mechanisms. This review will concentrate on Myc as a signaling mediator in the mammary gland, discussing its regulation and function during normal development, as well as its activation and roles in breast cancer.