Epigenetic and transcriptional determinants of the human breast

Three-dimensional genome architecture coordinates key regulators of lineage specification in mammary epithelial cells

Although lineage-specific genes have been identified in the mammary gland, little is known about the contribution of the 3D genome organization to gene regulation in the epithelium. Here, we describe the chromatin landscape of the three major epithelial subsets through integration of long- and short-range chromatin interactions, accessibility, histone modifications, and gene expression. While basal genes display exquisite lineage specificity via distal enhancers, luminal-specific genes show widespread promoter priming in basal cells. Cell specificity in luminal progenitors is largely mediated through extensive chromatin interactions with super-enhancers in gene-body regions in addition to interactions with polycomb silencer elements. Moreover, lineage-specific transcription factors appear to be controlled through cell-specific chromatin interactivity. Finally, chromatin accessibility rather than interactivity emerged as a defining feature of the activation of quiescent basal stem cells. This work provides a comprehensive resource for understanding the role of higher-order chromatin interactions in cell-fate specification and differentiation in the adult mouse mammary gland.

The Tumor Suppressor BRCA1/2, Cancer Susceptibility and Genome Instability in Gynecological and Mammary Cancers

BRCA1 and BRCA2 germline alterations highly predispose women to breast and ovarian cancers. They are mostly found within the TNBC (Triple-Negative Breast Cancer) and the HGSOC (High-Grade Serous Ovarian Carcinoma) subsets, known by an aggressive phenotype, the lack of therapeutic targets and poor prognosis. Importantly, there is an increased risk for cervical cancer in BRCA1 and BRCA2 mutation carriers that raises questions about the link between the HPV-driven genome instability and BRCA1 and BRCA2 germline mutations. Clinical, preclinical, and in vitro studies explained the increased risk for breast and ovarian cancers by genome instability resulting from the lack or loss of many functions related to BRCA1 or BRCA2 proteins such as DNA damage repair, stalled forks and R-loops resolution, transcription regulation, cell cycle control, and oxidative stress. In this review, we decipher the relationship between BRCA1/2 alterations and genomic instability leading to gynecomammary cancers through results from patients, mice, and cell lines. Understanding the early events of BRCA1/2-driven genomic instability in gynecomammary cancers would help to find new biomarkers for early diagnosis, improve the sensitivity of emerging therapies such as PARP inhibitors, and reveal new potential therapeutic targets.

SaBrcada: Survival Intervals Prediction for Breast Cancer Patients by Dimension Raising and Age Stratification

(1) Background: Breast cancer is the second leading cause of cancer death among women. The accurate prediction of survival intervals will help physicians make informed decisions about treatment strategies or the use of palliative care. (2) Methods: Gene expression is predictive and correlates to patient prognosis. To establish a reliable prediction tool, we collected a total of 1187 RNA-seq data points from breast cancer patients (median age 58 years) in Fragments Per Kilobase Million (FPKM) format from the TCGA database. Among them, we selected 144 patients with date of death information to establish the SaBrcada-AD dataset. We first normalized the SaBrcada-AD dataset to TPM to build the survival prediction model SaBrcada. After normalization and dimension raising, we used the differential gene expression data to test eight different deep learning architectures. Considering the effect of age on prognosis, we also performed a stratified random sampling test on all ages between the lower and upper quartiles of patient age, 48 and 69 years; (3) Results: Stratifying by age 61, the performance of SaBrcada built by GoogLeNet was improved to a highest accuracy of 0.798. We also built a free website tool to provide five predicted survival periods: within six months, six months to one year, one to three years, three to five years, or over five years, for clinician reference. (4) Conclusions: We built the prediction model, SaBrcada, and the website tool of the same name for breast cancer survival analysis. Through these models and tools, clinicians will be provided with survival interval information as a basis for formulating precision medicine.

Cell type and stage specific transcriptional, chromatin and cell-cell communication landscapes in the mammary gland

The mammary gland (MG) is composed of three main epithelial lineages, the basal cells (BC), the estrogen receptor (ER) positive luminal cells (ER+ LC), and the ER negative LC (ER- LC). Defining the cell identity of each lineage and how it is modulated throughout the different stages of life is important to understand how these cells function and communicate throughout life. Here, we used transgenic mice specifically labelling ER+ LC combined to cell surface markers to isolate with high purity the 3 distinct cell lineages of the mammary gland and defined their expression profiles and chromatin landscapes by performing bulk RNAseq and ATACseq of these isolated populations in puberty, adulthood and mid-pregnancy. Our analysis identified conserved genes, ligands and transcription factor (TF) associated with a specific lineage throughout life as well as genes, ligands and TFs specific for a particular stage of the MG. In summary, our study identified genes and TF network associated with the identity, function and cell-cell communication of the different epithelial lineages of the MG at different stages of life.

RUNX1 colludes with NOTCH1 to reprogram chromatin in T cell acute lymphoblastic leukemia

Runt-related transcription factor 1 (RUNX1) is oncogenic in diverse types of leukemia and epithelial cancers where its expression is associated with poor prognosis. Current models suggest that RUNX1 cooperates with other oncogenic factors (e.g., NOTCH1, TAL1) to drive the expression of proto-oncogenes in T cell acute lymphoblastic leukemia (T-ALL) but the molecular mechanisms controlled by RUNX1 and its cooperation with other factors remain unclear. Integrative chromatin and transcriptional analysis following inhibition of RUNX1 and NOTCH1 revealed a surprisingly widespread role of RUNX1 in the establishment of global H3K27ac levels and that RUNX1 is required by NOTCH1 for cooperative transcription activation of key NOTCH1 target genes including MYC, DTX1, HES4, IL7R, and NOTCH3. Super-enhancers were preferentially sensitive to RUNX1 knockdown and RUNX1-dependent super-enhancers were disrupted following the treatment of a pan-BET inhibitor, I-BET151.

Hedgehog signaling activates a mammalian heterochronic gene regulatory network controlling differentiation timing across lineages

Many developmental signaling pathways have been implicated in lineage-specific differentiation; however, mechanisms that explicitly control differentiation timing remain poorly defined in mammals. We report that murine Hedgehog signaling is a heterochronic pathway that determines the timing of progenitor differentiation. Hedgehog activity was necessary to prevent premature differentiation of second heart field (SHF) cardiac progenitors in mouse embryos, and the Hedgehog transcription factor GLI1 was sufficient to delay differentiation of cardiac progenitors in vitro. GLI1 directly activated a de novo progenitor-specific network in vitro, akin to that of SHF progenitors in vivo, which prevented the onset of the cardiac differentiation program. A Hedgehog signaling-dependent active-to-repressive GLI transition functioned as a differentiation timer, restricting the progenitor network to the SHF. GLI1 expression was associated with progenitor status across germ layers, and it delayed the differentiation of neural progenitors in vitro, suggesting a broad role for Hedgehog signaling as a heterochronic pathway.

Polycomb contraction differentially regulates terminal human hematopoietic differentiation programs

Background

Lifelong production of the many types of mature blood cells from less differentiated progenitors is a hierarchically ordered process that spans multiple cell divisions. The nature and timing of the molecular events required to integrate the environmental signals, transcription factor activity, epigenetic modifications, and changes in gene expression involved are thus complex and still poorly understood. To address this gap, we generated comprehensive reference epigenomes of 8 phenotypically defined subsets of normal human cord blood.

Results

We describe a striking contraction of H3K27me3 density in differentiated myelo-erythroid cells that resembles a punctate pattern previously ascribed to pluripotent embryonic stem cells. Phenotypically distinct progenitor cell types display a nearly identical repressive H3K27me3 signature characterized by large organized chromatin K27-modification domains that are retained by mature lymphoid cells but lost in terminally differentiated monocytes and erythroblasts. We demonstrate that inhibition of polycomb group members predicted to control large organized chromatin K27-modification domains influences lymphoid and myeloid fate decisions of primary neonatal hematopoietic progenitors in vitro. We further show that a majority of active enhancers appear in early progenitors, a subset of which are DNA hypermethylated and become hypomethylated and induced during terminal differentiation.

Conclusion

Primitive human hematopoietic cells display a unique repressive H3K27me3 signature that is retained by mature lymphoid cells but is lost in monocytes and erythroblasts. Intervention data implicate that control of this chromatin state change is a requisite part of the process whereby normal human hematopoietic progenitor cells make lymphoid and myeloid fate decisions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01315-1.

The dynamic alteration of transcriptional regulation by crucial TFs during tumorigenesis of gastric cancer

Background

The mechanisms of Gastric cancer (GC) initiation and progression are complicated, at least partly owing to the dynamic changes of gene regulation during carcinogenesis. Thus, investigations on the changes in regulatory networks can improve the understanding of cancer development and provide novel insights into the molecular mechanisms of cancer.

Methods

Differential co-expression analysis (DCEA), differential gene regulation network (GRN) modeling and differential regulation analysis (DRA) were integrated to detect differential transcriptional regulation events between gastric normal mucosa and cancer samples based on GSE54129 dataset. Cytological experiments and IHC staining assays were used to validate the dynamic changes of CREB1 regulated targets in different stages.

Results

A total of 1955 differentially regulated genes (DRGs) were identified and prioritized in a quantitative way. Among the top 1% DRGs, 14 out of 19 genes have been reported to be GC relevant. The four transcription factors (TFs) among the top 1% DRGs, including CREB1, BPTF, GATA6 and CEBPA, were regarded as crucial TFs relevant to GC progression. The differentially regulated links (DRLs) around the four crucial TFs were then prioritized to generate testable hypotheses on the differential regulation mechanisms of gastric carcinogenesis. To validate the dynamic alterations of gene regulation patterns of crucial TFs during GC progression, we took CREB1 as an example to screen its differentially regulated targets by using cytological and IHC staining assays. Eventually, TCEAL2 and MBNL1 were proved to be differentially regulated by CREB1 during tumorigenesis of gastric cancer.

Conclusions

By combining differential networking information and molecular cell experiments verification, testable hypotheses on the regulation mechanisms of GC around the core TFs and their top ranked DRLs were generated. Since TCEAL2 and MBNL1 have been reported to be potential therapeutic targets in SCLC and breast cancer respectively, their translation values in GC are worthy of further investigation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s10020-022-00468-7.

Functional annotation of breast cancer risk loci: current progress and future directions

Genome-wide association studies coupled with large-scale replication and fine-scale mapping studies have identified more than 150 genomic regions that are associated with breast cancer risk. Here, we review efforts to translate these findings into a greater understanding of disease mechanism. Our review comes in the context of a recently published fine-scale mapping analysis of these regions, which reported 352 independent signals and a total of 13,367 credible causal variants. The vast majority of credible causal variants map to noncoding DNA, implicating regulation of gene expression as the mechanism by which functional variants influence risk. Accordingly, we review methods for defining candidate-regulatory sequences, methods for identifying putative target genes and methods for linking candidate-regulatory sequences to putative target genes. We provide a summary of available data resources and identify gaps in these resources. We conclude that while much work has been done, there is still much to do. There are, however, grounds for optimism; combining statistical data from fine-scale mapping with functional data that are more representative of the normal "at risk" breast, generated using new technologies, should lead to a greater understanding of the mechanisms that influence an individual woman's risk of breast cancer.

Cell-type-specific epigenomic variations associated with BRCA1 mutation in pre-cancer human breast tissues

BRCA1 germline mutation carriers are predisposed to breast cancers. Epigenomic regulations have been known to strongly interact with genetic variations and potentially mediate biochemical cascades involved in tumorigenesis. Due to the cell-type specificity of epigenomic features, profiling of individual cell types is critical for understanding the molecular events in various cellular compartments within complex breast tissue. Here, we produced cell-type-specific profiles of genome-wide histone modifications including H3K27ac and H3K4me3 in basal, luminal progenitor, mature luminal and stromal cells extracted from a small pilot cohort of pre-cancer BRCA1 mutation carriers (BRCA1^mut/+) and non-carriers (BRCA1^+/+), using a low-input ChIP-seq technology that we developed. We discovered that basal and stromal cells present the most extensive epigenomic differences between mutation carriers (BRCA1^mut/+) and non-carriers (BRCA1^+/+), while luminal progenitor and mature luminal cells are relatively unchanged with the mutation. Furthermore, the epigenomic changes in basal cells due to BRCA1 mutation appear to facilitate their transformation into luminal progenitor cells. Taken together, epigenomic regulation plays an important role in the case of BRCA1 mutation for shaping the molecular landscape that facilitates tumorigenesis.

Towards resolution of the intron retention paradox in breast cancer

BACKGROUND

After many years of neglect in the field of alternative splicing, the importance of intron retention (IR) in cancer has come into focus following landmark discoveries of aberrant IR patterns in cancer. Many solid and liquid tumours are associated with drastic increases in IR, and such patterns have been pursued as both biomarkers and therapeutic targets. Paradoxically, breast cancer (BrCa) is the only tumour type in which IR is reduced compared to adjacent normal breast tissue.

METHODS

In this study, we have conducted a pan-cancer analysis of IR with emphasis on BrCa and its subtypes. We explored mechanisms that could cause aberrant and pathological IR and clarified why normal breast tissue has unusually high IR.

RESULTS

Strikingly, we found that aberrantly decreasing IR in BrCa can be largely attributed to normal breast tissue having the highest occurrence of IR events compared to other healthy tissues. Our analyses suggest that low numbers of IR events in breast tumours are associated with poor prognosis, particularly in the luminal B subtype. Interestingly, we found that IR frequencies negatively correlate with cell proliferation in BrCa cells, i.e. rapidly dividing tumour cells have the lowest number of IR events. Aberrant RNA-binding protein expression and changes in tissue composition are among the causes of aberrantly decreasing IR in BrCa.

CONCLUSIONS

Our results suggest that IR should be considered for therapeutic manipulation in BrCa patients with aberrantly low IR levels and that further work is needed to understand the cause and impact of high IR in other tumour types.

Therapeutic Targeting of Alternative RNA Splicing in Gastrointestinal Malignancies and Other Cancers

Recent comprehensive genomic studies including single-cell RNA sequencing and characterization have revealed multiple processes by which protein-coding and noncoding RNA processing are dysregulated in many cancers. More specifically, the abnormal regulation of mRNA and precursor mRNA (pre-mRNA) processing, which includes the removal of introns by splicing, is frequently altered in tumors, producing multiple different isoforms and diversifying protein expression. These alterations in RNA processing result in numerous cancer-specific mRNAs and pathogenically spliced events that generate altered levels of normal proteins or proteins with new functions, leading to the activation of oncogenes or the inactivation of tumor suppressor genes. Abnormally spliced pre-mRNAs are also associated with resistance to cancer treatment, and certain cancers are highly sensitive to the pharmacological inhibition of splicing. The discovery of these alterations in RNA processing has not only provided new insights into cancer pathogenesis but identified novel therapeutic vulnerabilities and therapeutic opportunities in targeting these aberrations in various ways (e.g., small molecules, splice-switching oligonucleotides (SSOs), and protein therapies) to modulate alternative RNA splicing or other RNA processing and modification mechanisms. Some of these strategies are currently progressing toward clinical development or are already in clinical trials. Additionally, tumor-specific neoantigens produced from these pathogenically spliced events and other abnormal RNA processes provide a potentially extensive source of tumor-specific therapeutic antigens (TAs) for targeted cancer immunotherapy. Moreover, a better understanding of the molecular mechanisms associated with aberrant RNA processes and the biological impact they play might provide insights into cancer initiation, progression, and metastasis. Our goal is to highlight key alternative RNA splicing and processing mechanisms and their roles in cancer pathophysiology as well as emerging therapeutic alternative splicing targets in cancer, particularly in gastrointestinal (GI) malignancies.

Eugenol modulates genomic methylation and inactivates breast cancer-associated fibroblasts through E2F1-dependent downregulation of DNMT1/DNMT3A

Active cancer-associated fibroblasts (CAFs) are major components of the tumor microenvironment, which promote carcinogenesis and modulate response to therapy. Therefore, targeting these cells or reducing their paracrine pro-carcinogenic effects could be of great therapeutic value. To this end, we sought to investigate the effect of eugenol, a natural phenolic molecule, on active breast CAFs. We have shown that decitabine (5-Aza-2'-deoxycytidine, DAC) and eugenol inhibit the expression of the DNA methyltransferase genes DNMT1 and DNMT3A at both the protein and mRNA levels in breast CAF cells. While the effect of eugenol was persistent, DAC had only a transient inhibitory effect on the mRNA level of both DNMT genes. Furthermore, eugenol and DAC suppressed the invasive/migratory and proliferative potential of CAF cells as well as their paracrine pro-carcinogenic effects both in vitro and in humanized orthotopic tumor xenografts. Interestingly, these inhibitory effects of decitabine and eugenol were mediated through E2F1 downregulation. Indeed, ectopic expression of E2F1 upregulated both genes and attenuated the effects of eugenol. Additionally, we provide clear evidence that eugenol, like DAC, strongly modulates the methylation pattern in active CAF cells, through methylating several oncogenes and demethylating various important tumor suppressor genes, which affected their mRNA expression levels. Importantly, the E2F1 promoter was also hypermethylated and the gene downregulated in response to eugenol. Together, these findings show that the active features of breast CAF cells can be normalized through eugenol-dependent targeting of DNMT1/DNMT3A and the consequent modulation in gene methylation.

Epigenomic programming in early fetal brain development

Aim: To provide a comprehensive understanding of gene regulatory networks in the developing human brain and a foundation for interpreting pathogenic deregulation. Materials & methods: We generated reference epigenomes and transcriptomes of dissected brain regions and primary neural progenitor cells (NPCs) derived from cortical and ganglionic eminence tissues of four normal human fetuses. Results: Integration of these data across developmental stages revealed a directional increase in active regulatory states, transcription factor activities and gene transcription with developmental stage. Consistent with differences in their biology, NPCs derived from cortical and ganglionic eminence regions contained common, region specific, and gestational week specific regulatory states. Conclusion: We provide a high-resolution regulatory network for NPCs from different brain regions as a comprehensive reference for future studies.

The changing paradigm of intron retention: regulation, ramifications and recipes

Intron retention (IR) is a form of alternative splicing that has long been neglected in mammalian systems although it has been studied for decades in non-mammalian species such as plants, fungi, insects and viruses. It was generally assumed that mis-splicing, leading to the retention of introns, would have no physiological consequence other than reducing gene expression by nonsense-mediated decay. Relatively recent landmark discoveries have highlighted the pivotal role that IR serves in normal and disease-related human biology. Significant technical hurdles have been overcome, thereby enabling the robust detection and quantification of IR. Still, relatively little is known about the cis- and trans-acting modulators controlling this phenomenon. The fate of an intron to be, or not to be, retained in the mature transcript is the direct result of the influence exerted by numerous intrinsic and extrinsic factors at multiple levels of regulation. These factors have altered current biological paradigms and provided unexpected insights into the transcriptional landscape. In this review, we discuss the regulators of IR and methods to identify them. Our focus is primarily on mammals, however, we broaden the scope to non-mammalian organisms in which IR has been shown to be biologically relevant.

Epigenetics in Inflammatory Breast Cancer: Biological Features and Therapeutic Perspectives

Evidence has emerged implicating epigenetic alterations in inflammatory breast cancer (IBC) origin and progression. IBC is a rare and rapidly progressing disease, considered the most aggressive type of breast cancer (BC). At clinical presentation, IBC is characterized by diffuse erythema, skin ridging, dermal lymphatic invasion, and peau d'orange aspect. The widespread distribution of the tumor as emboli throughout the breast and intra- and intertumor heterogeneity is associated with its poor prognosis. In this review, we highlighted studies documenting the essential roles of epigenetic mechanisms in remodeling chromatin and modulating gene expression during mammary gland differentiation and the development of IBC. Compiling evidence has emerged implicating epigenetic changes as a common denominator linking the main risk factors (socioeconomic status, environmental exposure to endocrine disruptors, racial disparities, and obesity) with IBC development. DNA methylation changes and their impact on the diagnosis, prognosis, and treatment of IBC are also described. Recent studies are focusing on the use of histone deacetylase inhibitors as promising epigenetic drugs for treating IBC. All efforts must be undertaken to unravel the epigenetic marks that drive this disease and how this knowledge could impact strategies to reduce the risk of IBC development and progression.

LDAPred: A Method Based on Information Flow Propagation and a Convolutional Neural Network for the Prediction of Disease-Associated lncRNAs

Long non-coding RNAs (lncRNAs) play a crucial role in the pathogenesis and development of complex diseases. Predicting potential lncRNA-disease associations can improve our understanding of the molecular mechanisms of human diseases and help identify biomarkers for disease diagnosis, treatment, and prevention. Previous research methods have mostly integrated the similarity and association information of lncRNAs and diseases, without considering the topological structure information among these nodes, which is important for predicting lncRNA-disease associations. We propose a method based on information flow propagation and convolutional neural networks, called LDAPred, to predict disease-related lncRNAs. LDAPred not only integrates the similarities, associations, and interactions among lncRNAs, diseases, and miRNAs, but also exploits the topological structures formed by them. In this study, we construct a dual convolutional neural network-based framework that comprises the left and right sides. The embedding layer on the left side is established by utilizing lncRNA, miRNA, and disease-related biological premises. On the right side of the frame, multiple types of similarity, association, and interaction relationships among lncRNAs, diseases, and miRNAs are calculated based on information flow propagation on the bi-layer networks, such as the lncRNA-disease network. They contain the network topological structure and they are learned by the right side of the framework. The experimental results based on five-fold cross-validation indicate that LDAPred performs better than several state-of-the-art methods. Case studies on breast cancer, colon cancer, and osteosarcoma further demonstrate LDAPred's ability to discover potential lncRNA-disease associations.

Mammary stem cells and progenitors: targeting the roots of breast cancer for prevention

Breast cancer prevention is daunting, yet not an unsurmountable goal. Mammary stem and progenitors have been proposed as the cells-of-origin in breast cancer. Here, we present the concept of limiting these breast cancer precursors as a risk reduction approach in high-risk women. A wealth of information now exists for phenotypic and functional characterization of mammary stem and progenitor cells in mouse and human. Recent work has also revealed the hormonal regulation of stem/progenitor dynamics as well as intrinsic lineage distinctions between mammary epithelial populations. Leveraging these insights, molecular marker-guided chemoprevention is an achievable reality.

Transcriptional regulation of normal human mammary cell heterogeneity and its perturbation in breast cancer

The mammary gland in adult women consists of biologically distinct cell types that differ in their surface phenotypes. Isolation and molecular characterization of these subpopulations of mammary cells have provided extensive insights into their different transcriptional programs and regulation. This information is now serving as a baseline for interpreting the heterogeneous features of human breast cancers. Examination of breast cancer mutational profiles further indicates that most have undergone a complex evolutionary process even before being detected. The consequent intra-tumoral as well as inter-tumoral heterogeneity of these cancers thus poses major challenges to deriving information from early and hence likely pervasive changes in potential therapeutic interest. Recently described reproducible and efficient methods for generating human breast cancers de novo in immunodeficient mice transplanted with genetically altered primary cells now offer a promising alternative to investigate initial stages of human breast cancer development. In this review, we summarize current knowledge about key transcriptional regulatory processes operative in these partially characterized subpopulations of normal human mammary cells and effects of disrupting these processes in experimentally produced human breast cancers.

Overlooked roles of DNA damage and maternal age in generating human germline mutations

The textbook view that most germline mutations in mammals arise from replication errors is indirectly supported by the fact that there are both more mutations and more cell divisions in the male than in the female germline. When analyzing large de novo mutation datasets in humans, we find multiple lines of evidence that call that view into question. Notably, despite the drastic increase in the ratio of male to female germ cell divisions after the onset of spermatogenesis, even young fathers contribute three times more mutations than young mothers, and this ratio barely increases with parental age. This surprising finding points to a substantial contribution of damage-induced mutations. Indeed, C-to-G transversions and CpG transitions, which together constitute over one-fourth of all base substitution mutations, show genomic distributions and sex-specific age dependencies indicative of double-strand break repair and methylation-associated damage, respectively. Moreover, we find evidence that maternal age at conception influences the mutation rate both because of the accumulation of damage in oocytes and potentially through an influence on the number of postzygotic mutations in the embryo. These findings reveal underappreciated roles of DNA damage and maternal age in the genesis of human germline mutations.

Early Parity Epigenetic Footprint of FOXA1 Gene Body in Normal Breast Tissue of Iranian Women

BACKGROUND

Young age at first full-term pregnancy (FFTP) is an important factor in breast cancer risk reduction. It is postulated that this protective effect is the result of stable molecular signatures imprinted by physiological process of pregnancy, but the molecular mechanism of this protective role is unclear. The aim of the current study was to identify the effect of early FFTP on methylation status of FOXA1 gene body. FOXA1 is an essential transcription factor for mammary gland development and estrogen responsiveness of breast tissue.

METHODS

Fresh frozen normal breast tissues (n = 51) were collected from Iranian women who underwent cosmetic mammoplasty (27 nulliparous women and 24 parous women who have experienced first pregnancy before the age of 25). DNA was extracted and then methylated DNA immunoprecipitation (MeDIP) real-time PCR was used to assess FOXA1 gene body methylation.

RESULTS

Our results revealed that FOXA1 methylation level is significantly higher in early parous compared with nulliparous group (p = 0.041).

CONCLUSION

Our study provides new hint about the association between early FFTP and epigenetic modifications within gene body of FOXA1 in normal breast tissue. More investigation is required for clarifying molecular mechanisms underlying this association in order to develop breast cancer prevention strategies.

MicroRNA-196a is regulated by ER and is a prognostic biomarker in ER+ breast cancer

Background

MicroRNAs are potent post-transcriptional regulators involved in all hallmarks of cancer. Mir-196a is transcribed from two loci and has been implicated in a wide range of developmental and pathogenic processes, with targets including Hox, Fox, Cdk inhibitors and annexins. Genetic variants and altered expression of MIR196A are associated with risk and progression of multiple cancers including breast cancer, however little is known about the regulation of the genes encoding this miRNA, nor the impact of variants therein.

Methods

Genomic data and chromatin interaction analysis were used to discover functional promoter and enhancer elements for MIR196A. Expression data were used to associate MIR196A with mechanisms of resistance, breast cancer subtypes and prognosis.

Results

Here we demonstrate that MIR196A displays complex and dynamic expression patterns, in part controlled by long-range transcriptional regulation between promoter and enhancer elements bound by ERα. Expression of this miRNA is significantly increased in drug-resistant models of hormone-receptor positive disease. The expression of MIR196A also proves to be a robust prognostic factor for patients with advanced and post-menopausal ER+ disease.

Conclusion

This work sheds light on the normal and abnormal regulation of MIR196A and provides a novel stratification method for therapeutically resistant breast cancer.

DNA methylation/hydroxymethylation regulate gene expression and alternative splicing during terminal granulopoiesis

Aim: To determine whether epigenetic modifications of DNA regulate gene expression and alternative splicing during terminal granulopoiesis. Materials & methods: Using whole genome bisulfite sequencing, reduced representation hydroxymethylation profiling and mRNA sequencing, we compare changes in DNA methylation, DNA hydroxymethylation, gene expression and alternative splicing in mouse promyelocytes and granulocytes. Results & conclusion: We show reduced DNA methylation at the promoters and enhancers of key granulopoiesis genes, indicating a regulatory role in the activation of lineage-specific genes during differentiation. Notably, increased DNA hydroxymethylation in exons is associated with preferential inclusion of specific exons in granulocytes. Overall, DNA methylation and hydroxymethylation changes at particular genomic loci may play specific roles in gene regulation or alternative splicing during terminal granulopoiesis. Data deposition: Whole genome bisulfite sequencing of mouse promyelocytes and granulocytes: Gene Expression Omnibus (GSE85517); mRNA sequencing of mouse promyelocytes and granulocytes: Gene Expression Omnibus (GSE48307); reduced representation 5-hydroxymethylation profiling of mouse promyelocytes and granulocytes: Bioproject (PRJNA495696).

BRCA1-Dependent Transcriptional Regulation: Implication in Tissue-Specific Tumor Suppression

Germ-line mutations in breast cancer susceptibility gene 1 (BRCA1) predominantly predispose women to breast and ovarian cancers. BRCA1 is best known for its functions in maintenance of genomic integrity including repairing DNA double-strand breaks through homologous recombination and suppressing DNA replication stress. However, whether these universally important BRCA1 functions in maintenance of genomic stability are sufficient to account for its tissue-specific tumor-suppressing function remains unclear. Accumulating evidence indicates that there are previously underappreciated roles of BRCA1 in transcriptional regulation and chromatin remodeling. In this review, we discuss the functional significance of interactions between BRCA1 and various transcription factors, its role in epigenetic regulation and chromatin dynamics, and BRCA1-dependent crosstalk between the machineries of transcription and genome integrity. Furthermore, we propose a model of how transcriptional regulation could contribute to tissue-dependent tumor-suppressing function of BRCA1.

Silencing of retrotransposon-derived imprinted gene RTL1 is the main cause for postimplantational failures in mammalian cloning

Substantial rates of fetal loss plague all in vitro procedures involving embryo manipulations, including human-assisted reproduction, and are especially problematic for mammalian cloning where over 90% of reconstructed nuclear transfer embryos are typically lost during pregnancy. However, the epigenetic mechanism of these pregnancy failures has not been well described. Here we performed methylome and transcriptome analyses of pig induced pluripotent stem cells and associated cloned embryos, and revealed that aberrant silencing of imprinted genes, in particular the retrotransposon-derived RTL1 gene, is the principal epigenetic cause of pregnancy failure. Remarkably, restoration of RTL1 expression in pig induced pluripotent stem cells rescued fetal loss. Furthermore, in other mammals, including humans, low RTL1 levels appear to be the main epigenetic cause of pregnancy failure.

Differentiation-state plasticity is a targetable resistance mechanism in basal-like breast cancer

Intratumoral heterogeneity in cancers arises from genomic instability and epigenomic plasticity and is associated with resistance to cytotoxic and targeted therapies. We show here that cell-state heterogeneity, defined by differentiation-state marker expression, is high in triple-negative and basal-like breast cancer subtypes, and that drug tolerant persister (DTP) cell populations with altered marker expression emerge during treatment with a wide range of pathway-targeted therapeutic compounds. We show that MEK and PI3K/mTOR inhibitor-driven DTP states arise through distinct cell-state transitions rather than by Darwinian selection of preexisting subpopulations, and that these transitions involve dynamic remodeling of open chromatin architecture. Increased activity of many chromatin modifier enzymes, including BRD4, is observed in DTP cells. Co-treatment with the PI3K/mTOR inhibitor BEZ235 and the BET inhibitor JQ1 prevents changes to the open chromatin architecture, inhibits the acquisition of a DTP state, and results in robust cell death in vitro and xenograft regression in vivo.

Resistance to therapy can be driven by intratumoral heterogeneity. Here, the authors show that drug tolerant persistent cell populations emerge during treatment, and these emergent populations arise through epigenetically mediated cell state transitions rather than sub population selection.

Identification of protein kinase inhibitors to reprogram breast cancer cells

Direct reversion of cancers into normal-like tissues is an ideal strategy for cancer treatment. Recent reports have showed that defined transcription factors can induce reprogramming of cancer cells into pluripotent stem cells, supporting this notion. Here, we have developed a reprogramming method that uses a conceptually unique strategy for breast cancer cell treatment. We have screened a kinase inhibitor library and found that Rho-associated protein kinase (ROCK) and mammalian target of rapamycin (mTOR) kinase inhibitors can substitute for all transcription factors to be sufficient to reprogram breast cancer cells into progenitor cells. Furthermore, ROCK–mTOR inhibitors could reprogram breast cancer cells to another terminal lineage-adipogenic cells. Genome-wide transcriptional analysis shows that the induced fat-like cells have a profile different from breast cancer cells and similar to that of normal adipocytes. In vitro and in vivo tumorigenesis assays have shown that induced fat-like cells lose proliferation and tumorigenicity. Moreover, reprogramming treatment with ROCK–mTOR inhibitors prevents breast cancer local recurrence in mice. Currently, ROCK–mTOR inhibitors are already used as antitumor drugs in patients, thus, this reprogramming strategy has significant potential to move rapidly toward clinical trials for breast cancer treatment.

MiR-29b-1-5p is altered in BRCA1 mutant tumours and is a biomarker in basal-like breast cancer

Depletion of BRCA1 protein in mouse mammary glands results in defects in lactational development and increased susceptibility to mammary cancer. Extensive work has focussed on the role of BRCA1 in the normal breast and in the development of breast cancer, the cell of origin for BRCA1 tumours and the protein-coding genes altered in BRCA1 deficient cells. However, the role of non-coding RNAs in BRCA1-deficient cells is poorly understood. To evaluate miRNA expression in BRCA1 deficient mammary cells, RNA sequencing was performed on the mammary glands of Brca1 knockout mice. We identified 140 differentially expressed miRNAs, 9 of which were also differentially expressed in human BRCA1 breast tumours or familial non-BRCA1 patients and during normal gland development. We show that BRCA1 binds to putative cis-elements in promoter regions of the miRNAs with the potential to regulate their expression, and that four miRNAs (miR-29b-1-5p, miR-664, miR-16-2 and miR-744) significantly stratified the overall survival of basal-like tumours. Importantly the prognostic value of miR-29b-1-5p was higher in significance than several commonly used clinical biomarkers. These results emphasise the role of Brca1 in modulating expression of miRNAs and highlights the potential for BRCA1 regulated miRNAs to be informative biomarkers associated with BRCA1 loss and survival in breast cancer.

Epigenomics of mammary gland development

Differentiation of stem cells into highly specialised cells requires gene expression changes brought about by remodelling of the chromatin architecture. During this lineage-commitment process, the majority of DNA needs to be packaged into inactive heterochromatin, allowing only a subset of regulatory elements to remain open and functionally required genes to be expressed. Epigenetic mechanisms such as DNA methylation, post-translational modifications to histone tails, and nucleosome positioning all potentially contribute to the changes in higher order chromatin structure during differentiation. The mammary gland is a particularly useful model to study these complex epigenetic processes since the majority of its development is postnatal, the gland is easily accessible, and development occurs in a highly reproducible manner. Inappropriate epigenetic remodelling can also drive tumourigenesis; thus, insights into epigenetic remodelling during mammary gland development advance our understanding of breast cancer aetiology. We review the current literature surrounding DNA methylation and histone modifications in the developing mammary gland and its implications for breast cancer.

High-Resolution Single-Cell DNA Methylation Measurements Reveal Epigenetically Distinct Hematopoietic Stem Cell Subpopulations

Increasing evidence of functional and transcriptional heterogeneity in phenotypically similar cells examined individually has prompted interest in obtaining parallel methylome data. We describe the development and application of such a protocol to index-sorted murine and human hematopoietic cells that are highly enriched in their content of functionally defined stem cells. Utilizing an optimized single-cell bisulfite sequencing protocol, we obtained quantitative DNA methylation measurements of up to 5.7 million CpGs in single hematopoietic cells. In parallel, we developed an analytical strategy (PDclust) to define single-cell DNA methylation states through pairwise comparisons of single-CpG methylation measurements. PDclust revealed that a single-cell epigenetic state can be described by a small (<1%) stochastically sampled fraction of CpGs and that these states are reflective of cell identity and state. Using relationships revealed by PDclust, we derive near complete methylomes for epigenetically distinct subpopulations of hematopoietic cells enriched for functional stem cell content.

Mammary molecular portraits reveal lineage-specific features and progenitor cell vulnerabilities

Casey et al. integrate epigenomic, transcriptomic, and proteomic profiling of primary basal and luminal mammary cells to identify master epigenetic regulators of the mammary epithelium and uncover stem and progenitor cell vulnerabilities. They develop a pipeline to identify drugs that abrogate progenitor cell activity in normal and high-risk breast cancer patient samples in vitro and in vivo.

The mammary epithelium depends on specific lineages and their stem and progenitor function to accommodate hormone-triggered physiological demands in the adult female. Perturbations of these lineages underpin breast cancer risk, yet our understanding of normal mammary cell composition is incomplete. Here, we build a multimodal resource for the adult gland through comprehensive profiling of primary cell epigenomes, transcriptomes, and proteomes. We define systems-level relationships between chromatin–DNA–RNA–protein states, identify lineage-specific DNA methylation of transcription factor binding sites, and pinpoint proteins underlying progesterone responsiveness. Comparative proteomics of estrogen and progesterone receptor–positive and –negative cell populations, extensive target validation, and drug testing lead to discovery of stem and progenitor cell vulnerabilities. Top epigenetic drugs exert cytostatic effects; prevent adult mammary cell expansion, clonogenicity, and mammopoiesis; and deplete stem cell frequency. Select drugs also abrogate human breast progenitor cell activity in normal and high-risk patient samples. This integrative computational and functional study provides fundamental insight into mammary lineage and stem cell biology.

MiR-519d suppresses breast cancer tumorigenesis and metastasis via targeting MMP3

Breast cancer (BC) is the most common cause of death in women throughout the world. Although microRNAs (miRNAs) have been identified as novel regulators in carcinogenesis, there are still abundant hidden treasure needed to be excavated. In the present study, we found that miR-519d expression was remarkably decreased in both human BC tissues and MCF-7 cells. CCK8 and 5-Ethynyl-2'-deoxyuridine (EdU) assays were used to evaluate cell proliferation. Wound-healing and transwell assays were performed for detection of cell migration and invasion. The results demonstrated miR-519d overexpression dramatically suppressed MCF-7 cells proliferation, migration and invasion. While downregulation of miR-519d by miR-519d inhibitor substantially increased MCF-7 cell carcinogenesis. Further analysis identified Matrix Metalloproteinase-3 (MMP3) as a direct target of miR-519d. QRT-PCR and western blot results indicated the correlative expression of miR-519d and MMP3 in BC tissues and MCF-7 cells. In summary, our data uncovered the novel molecular interaction between miR-519d and MMP3, indicating a therapeutic strategy of miR-519d for BC.

Analysis of Normal Human Mammary Epigenomes Reveals Cell-Specific Active Enhancer States and Associated Transcription Factor Networks

The normal adult human mammary gland is a continuous bilayered epithelial system. Bipotent and myoepithelial progenitors are prominent and unique components of the outer (basal) layer. The inner (luminal) layer includes both luminal-restricted progenitors and a phenotypically separable fraction that lacks progenitor activity. We now report an epigenomic comparison of these three subsets with one another, with their associated stromal cells, and with three immortalized, non-tumorigenic human mammary cell lines. Each genome-wide analysis contains profiles for six histone marks, methylated DNA, and RNA transcripts. Analysis of these datasets shows that each cell type has unique features, primarily within genomic regulatory regions, and that the cell lines group together. Analyses of the promoter and enhancer profiles place the luminal progenitors in between the basal cells and the non-progenitor luminal subset. Integrative analysis reveals networks of subset-specific transcription factors.

Vitamin C induces specific demethylation of H3K9me2 in mouse embryonic stem cells via Kdm3a/b

Background

Histone methylation patterns regulate gene expression and are highly dynamic during development. The erasure of histone methylation is carried out by histone demethylase enzymes. We had previously shown that vitamin C enhances the activity of Tet enzymes in embryonic stem (ES) cells, leading to DNA demethylation and activation of germline genes.

Results

We report here that vitamin C induces a remarkably specific demethylation of histone H3 lysine 9 dimethylation (H3K9me2) in naïve ES cells. Vitamin C treatment reduces global levels of H3K9me2, but not other histone methylation marks analyzed, as measured by western blot, immunofluorescence and mass spectrometry. Vitamin C leads to widespread loss of H3K9me2 at large chromosomal domains as well as gene promoters and repeat elements. Vitamin C-induced loss of H3K9me2 occurs rapidly within 24 h and is reversible. Importantly, we found that the histone demethylases Kdm3a and Kdm3b are required for vitamin C-induced demethylation of H3K9me2. Moreover, we show that vitamin C-induced Kdm3a/b-mediated H3K9me2 demethylation and Tet-mediated DNA demethylation are independent processes at specific loci. Lastly, we document Kdm3a/b are partially required for the upregulation of germline genes by vitamin C.

Conclusions

These results reveal a specific role for vitamin C in histone demethylation in ES cells and document that DNA methylation and H3K9me2 cooperate to silence germline genes in pluripotent cells.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-017-0143-3) contains supplementary material, which is available to authorized users.

Attenuation of RNA polymerase II pausing mitigates BRCA1-associated R-loop accumulation and tumorigenesis

Most BRCA1-associated breast tumours are basal-like yet originate from luminal progenitors. BRCA1 is best known for its functions in double-strand break repair and resolution of DNA replication stress. However, it is unclear whether loss of these ubiquitously important functions fully explains the cell lineage-specific tumorigenesis. In vitro studies implicate BRCA1 in elimination of R-loops, DNA-RNA hybrid structures involved in transcription and genetic instability. Here we show that R-loops accumulate preferentially in breast luminal epithelial cells, not in basal epithelial or stromal cells, of BRCA1 mutation carriers. Furthermore, R-loops are enriched at the 5' end of those genes with promoter-proximal RNA polymerase II (Pol II) pausing. Genetic ablation of Cobra1, which encodes a Pol II-pausing and BRCA1-binding protein, ameliorates R-loop accumulation and reduces tumorigenesis in Brca1-knockout mouse mammary epithelium. Our studies show that Pol II pausing is an important contributor to BRCA1-associated R-loop accumulation and breast cancer development.

Intron retention is regulated by altered MeCP2-mediated splicing factor recruitment

While intron retention (IR) is considered a widely conserved and distinct mechanism of gene expression control, its regulation is poorly understood. Here we show that DNA methylation directly regulates IR. We also find reduced occupancy of MeCP2 near the splice junctions of retained introns, mirroring the reduced DNA methylation at these sites. Accordingly, MeCP2 depletion in tissues and cells enhances IR. By analysing the MeCP2 interactome using mass spectrometry and RNA co-precipitation, we demonstrate that decreased MeCP2 binding near splice junctions facilitates IR via reduced recruitment of splicing factors, including Tra2b, and increased RNA polymerase II stalling. These results suggest an association between IR and a slower rate of transcription elongation, which reflects inefficient splicing factor recruitment. In summary, our results reinforce the interdependency between alternative splicing involving IR and epigenetic controls of gene expression.

Intron retention is a conserved mechanism that controls gene expression but its regulation is poorly understood. Here, the authors provide evidence that DNA methylation regulates intron retention and find reduced MeCP2 occupancy and splicing factor recruitment near affected splice junctions.

Intron retention as a component of regulated gene expression programs

Intron retention has long been an exemplar of regulated splicing with case studies of individual events serving as models that provided key mechanistic insights into the process of splicing control. In organisms such as plants and budding yeast, intron retention is well understood as a major mechanism of gene expression regulation. In contrast, in mammalian systems, the extent and functional significance of intron retention have, until recently, remained greatly underappreciated. Technical challenges to the global detection and quantitation of transcripts with retained introns have often led to intron retention being overlooked or dismissed as “noise”. Now, however, with the wealth of information available from high-throughput deep sequencing, combined with focused computational and statistical analyses, we are able to distinguish clear intron retention patterns in various physiological and pathological contexts. Several recent studies have demonstrated intron retention as a central component of gene expression programs during normal development as well as in response to stress and disease. Furthermore, these studies revealed various ways in which intron retention regulates protein isoform production, RNA stability and translation efficiency, and rapid induction of expression via post-transcriptional splicing of retained introns. In this review, we highlight critical findings from these transcriptomic studies and discuss commonalties in the patterns prevalent in intron retention networks at the functional and regulatory levels.

Characterizing Genetic Susceptibility to Breast Cancer in Women of African Ancestry

Background: Genome-wide association studies have identified approximately 100 common genetic variants associated with breast cancer risk, the majority of which were discovered in women of European ancestry. Because of different patterns of linkage disequilibrium, many of these genetic markers may not represent signals in populations of African ancestry.Methods: We tested 74 breast cancer risk variants and conducted fine-mapping of these susceptibility regions in 6,522 breast cancer cases and 7,643 controls of African ancestry from three genetic consortia (AABC, AMBER, and ROOT).Results: Fifty-four of the 74 variants (73%) were found to have ORs that were directionally consistent with those previously reported, of which 12 were nominally statistically significant (P < 0.05). Through fine-mapping, in six regions (3p24, 12p11, 14q13, 16q12/FTO, 16q23, 19p13), we observed seven markers that better represent the underlying risk variant for overall breast cancer or breast cancer subtypes, whereas in another two regions (11q13, 16q12/TOX3), we identified suggestive evidence of signals that are independent of the reported index variant. Overlapping chromatin features and regulatory elements suggest that many of the risk alleles lie in regions with biological functionality.Conclusions: Through fine-mapping of known susceptibility regions, we have revealed alleles that better characterize breast cancer risk in women of African ancestry.Impact: The risk alleles identified represent genetic markers for modeling and stratifying breast cancer risk in women of African ancestry. Cancer Epidemiol Biomarkers Prev; 26(7); 1016-26. ©2017 AACR.

Subtype-specific CpG island shore methylation and mutation patterns in 30 breast cancer cell lines

Background

Aberrant epigenetic modifications, including DNA methylation, are key regulators of gene activity in tumorigenesis. Breast cancer is a heterogeneous disease, and large-scale analyses indicate that tumor from normal and benign tissues, as well as molecular subtypes of breast cancer, can be distinguished based on their distinct genomic, transcriptomic, and epigenomic profiles. In this study, we used affinity-based methylation sequencing data in 30 breast cancer cell lines representing functionally distinct cancer subtypes to investigate methylation and mutation patterns at the whole genome level.

Results

Our analysis revealed significant differences in CpG island (CpGI) shore methylation and mutation patterns among breast cancer subtypes. In particular, the basal-like B type, a highly aggressive form of the disease, displayed distinct CpGI shore hypomethylation patterns that were significantly associated with downstream gene regulation. We determined that mutation rates at CpG sites were highly correlated with DNA methylation status and observed distinct mutation rates among the breast cancer subtypes. These findings were validated by using targeted bisulfite sequencing of differentially expressed genes (n=85) among the cell lines.

Conclusions

Our results suggest that alterations in DNA methylation play critical roles in gene regulatory process as well as cytosine substitution rates at CpG sites in molecular subtypes of breast cancer.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-016-0356-2) contains supplementary material, which is available to authorized users.

The Ties That Bind: Mapping the Dynamic Enhancer-Promoter Interactome

Coupling chromosome conformation capture to molecular enrichment for promoter-containing DNA fragments enables the systematic mapping of interactions between individual distal regulatory sequences and their target genes. In this Minireview, we describe recent progress in the application of this technique and related complementary approaches to gain insight into the lineage- and cell-type-specific dynamics of interactions between regulators and gene promoters.

Mapping of Variable DNA Methylation Across Multiple Cell Types Defines a Dynamic Regulatory Landscape of the Human Genome

DNA methylation is an important epigenetic modification involved in many biological processes and diseases. Many studies have mapped DNA methylation changes associated with embryogenesis, cell differentiation, and cancer at a genome-wide scale. Our understanding of genome-wide DNA methylation changes in a developmental or disease-related context has been steadily growing. However, the investigation of which CpGs are variably methylated in different normal cell or tissue types is still limited. Here, we present an in-depth analysis of 54 single-CpG-resolution DNA methylomes of normal human cell types by integrating high-throughput sequencing-based methylation data. We found that the ratio of methylated to unmethylated CpGs is relatively constant regardless of cell type. However, which CpGs made up the unmethylated complement was cell-type specific. We categorized the 26,000,000 human autosomal CpGs based on their methylation levels across multiple cell types to identify variably methylated CpGs and found that 22.6% exhibited variable DNA methylation. These variably methylated CpGs formed 660,000 variably methylated regions (VMRs), encompassing 11% of the genome. By integrating a multitude of genomic data, we found that VMRs enrich for histone modifications indicative of enhancers, suggesting their role as regulatory elements marking cell type specificity. VMRs enriched for transcription factor binding sites in a tissue-dependent manner. Importantly, they enriched for GWAS variants, suggesting that VMRs could potentially be implicated in disease and complex traits. Taken together, our results highlight the link between CpG methylation variation, genetic variation, and disease risk for many human cell types.

Cell Fate Decisions During Breast Cancer Development

During the formation of breast cancer, many genes become altered as cells evolve progressively from normal to a pre-malignant to a malignant state of growth. How mutations in genes lead to specific subtypes of human breast cancer is only partially understood. Here we review how initial genetic or epigenetic alterations within mammary epithelial cells (MECs) can alter cell fate decisions and put pre-malignant cells on a path towards cancer development with specific phenotypes. Understanding the early stages of breast cancer initiation and progression and how normal developmental processes are hijacked during transformation has significant implications for improving early detection and prevention of breast cancer. In addition, insights gleaned from this understanding may also be important for developing subtype-specific treatment options.

Barcoding reveals complex clonal dynamics of de novo transformed human mammary cells

Most human breast cancers have diversified genomically and biologically by the time they become clinically evident. Early events involved in their genesis and the cellular context in which these events occur have thus been difficult to characterize. Here we present the first formal evidence of the shared and independent ability of basal cells and luminal progenitors, isolated from normal human mammary tissue and transduced with a single oncogene (KRAS(G12D)), to produce serially transplantable, polyclonal, invasive ductal carcinomas within 8 weeks of being introduced either subrenally or subcutaneously into immunodeficient mice. DNA barcoding of the initial cells revealed a dramatic change in the numbers and sizes of clones generated from them within 2 weeks, and the first appearance of many 'new' clones in tumours passaged into secondary recipients. Both primary and secondary tumours were phenotypically heterogeneous and primary tumours were categorized transcriptionally as 'normal-like'. This system challenges previous concepts that carcinogenesis in normal human epithelia is necessarily a slow process requiring the acquisition of multiple driver mutations. It also presents the first description of initial events that accompany the genesis and evolution of malignant human mammary cell populations, thereby contributing new understanding of the rapidity with which heterogeneity in their properties can develop.

Integrative analysis of 111 reference human epigenomes

The reference human genome sequence set the stage for studies of genetic variation and its association with human disease, but epigenomic studies lack a similar reference. To address this need, the NIH Roadmap Epigenomics Consortium generated the largest collection so far of human epigenomes for primary cells and tissues. Here we describe the integrative analysis of 111 reference human epigenomes generated as part of the programme, profiled for histone modification patterns, DNA accessibility, DNA methylation and RNA expression. We establish global maps of regulatory elements, define regulatory modules of coordinated activity, and their likely activators and repressors. We show that disease- and trait-associated genetic variants are enriched in tissue-specific epigenomic marks, revealing biologically relevant cell types for diverse human traits, and providing a resource for interpreting the molecular basis of human disease. Our results demonstrate the central role of epigenomic information for understanding gene regulation, cellular differentiation and human disease.