Epigenetic and transcriptional determinants of the human breast.

While significant effort has been dedicated to the characterization of epigenetic changes associated with prenatal differentiation, relatively little is known about the epigenetic changes that accompany post-natal differentiation where fully functional differentiated cell types with limited lifespans arise. Here we sought to address this gap by generating epigenomic and transcriptional profiles from primary human breast cell types isolated from disease-free human subjects. From these data we define a comprehensive human breast transcriptional network, including a set of myoepithelial- and luminal epithelial-specific intronic retention events. Intersection of epigenetic states with RNA expression from distinct breast epithelium lineages demonstrates that mCpG provides a stable record of exonic and intronic usage, whereas H3K36me3 is dynamic. We find a striking asymmetry in epigenomic reprogramming between luminal and myoepithelial cell types, with the genomes of luminal cells harbouring more than twice the number of hypomethylated enhancer elements compared with myoepithelial cells.

Intermediate DNA methylation is a conserved signature of genome regulation.

The role of intermediate methylation states in DNA is unclear. Here, to comprehensively identify regions of intermediate methylation and their quantitative relationship with gene activity, we apply integrative and comparative epigenomics to 25 human primary cell and tissue samples. We report 18,452 intermediate methylation regions located near 36% of genes and enriched at enhancers, exons and DNase I hypersensitivity sites. Intermediate methylation regions average 57% methylation, are predominantly allele-independent and are conserved across individuals and between mouse and human, suggesting a conserved function. These regions have an intermediate level of active chromatin marks and their associated genes have intermediate transcriptional activity. Exonic intermediate methylation correlates with exon inclusion at a level between that of fully methylated and unmethylated exons, highlighting gene context-dependent functions. We conclude that intermediate DNA methylation is a conserved signature of gene regulation and exon usage.

Regulatory network decoded from epigenomes of surface ectoderm-derived cell types.

Developmental history shapes the epigenome and biological function of differentiated cells. Epigenomic patterns have been broadly attributed to the three embryonic germ layers. Here we investigate how developmental origin influences epigenomes. We compare key epigenomes of cell types derived from surface ectoderm (SE), including keratinocytes and breast luminal and myoepithelial cells, against neural crest-derived melanocytes and mesoderm-derived dermal fibroblasts, to identify SE differentially methylated regions (SE-DMRs). DNA methylomes of neonatal keratinocytes share many more DMRs with adult breast luminal and myoepithelial cells than with melanocytes and fibroblasts from the same neonatal skin. This suggests that SE origin contributes to DNA methylation patterning, while shared skin tissue environment has limited effect on epidermal keratinocytes. Hypomethylated SE-DMRs are in proximity to genes with SE relevant functions. They are also enriched for enhancer- and promoter-associated histone modifications in SE-derived cells, and for binding motifs of transcription factors important in keratinocyte and mammary gland biology. Thus, epigenomic analysis of cell types with common developmental origin reveals an epigenetic signature that underlies a shared gene regulatory network.

Functional DNA methylation differences between tissues, cell types, and across individuals discovered using the M&M algorithm.

DNA methylation plays key roles in diverse biological processes such as X chromosome inactivation, transposable element repression, genomic imprinting, and tissue-specific gene expression. Sequencing-based DNA methylation profiling provides an unprecedented opportunity to map and compare complete DNA methylomes. This includes one of the most widely applied technologies for measuring DNA methylation: methylated DNA immunoprecipitation followed by sequencing (MeDIP-seq), coupled with a complementary method, methylation-sensitive restriction enzyme sequencing (MRE-seq). A computational approach that integrates data from these two different but complementary assays and predicts methylation differences between samples has been unavailable. Here, we present a novel integrative statistical framework M&M (for integration of MeDIP-seq and MRE-seq) that dynamically scales, normalizes, and combines MeDIP-seq and MRE-seq data to detect differentially methylated regions. Using sample-matched whole-genome bisulfite sequencing (WGBS) as a gold standard, we demonstrate superior accuracy and reproducibility of M&M compared to existing analytical methods for MeDIP-seq data alone. M&M leverages the complementary nature of MeDIP-seq and MRE-seq data to allow rapid comparative analysis between whole methylomes at a fraction of the cost of WGBS. Comprehensive analysis of nineteen human DNA methylomes with M&M reveals distinct DNA methylation patterns among different tissue types, cell types, and individuals, potentially underscoring divergent epigenetic regulation at different scales of phenotypic diversity. We find that differential DNA methylation at enhancer elements, with concurrent changes in histone modifications and transcription factor binding, is common at the cell, tissue, and individual levels, whereas promoter methylation is more prominent in reinforcing fundamental tissue identities.

DNA hypomethylation within specific transposable element families associates with tissue-specific enhancer landscape.

Transposable element (TE)-derived sequences comprise half of the human genome and DNA methylome and are presumed to be densely methylated and inactive. Examination of genome-wide DNA methylation status within 928 TE subfamilies in human embryonic and adult tissues identified unexpected tissue-specific and subfamily-specific hypomethylation signatures. Genes proximal to tissue-specific hypomethylated TE sequences were enriched for functions important for the relevant tissue type, and their expression correlated strongly with hypomethylation within the TEs. When hypomethylated, these TE sequences gained tissue-specific enhancer marks, including monomethylation of histone H3 at lysine 4 (H3K4me1) and occupancy by p300, and a majority exhibited enhancer activity in reporter gene assays. Many such TEs also harbored binding sites for transcription factors that are important for tissue-specific functions and showed evidence of evolutionary selection. These data suggest that sequences derived from TEs may be responsible for wiring tissue type-specific regulatory networks and may have acquired tissue-specific epigenetic regulation.

The clonal and mutational evolution spectrum of primary triple-negative breast cancers.

Primary triple-negative breast cancers (TNBCs), a tumour type defined by lack of oestrogen receptor, progesterone receptor and ERBB2 gene amplification, represent approximately 16% of all breast cancers. Here we show in 104 TNBC cases that at the time of diagnosis these cancers exhibit a wide and continuous spectrum of genomic evolution, with some having only a handful of coding somatic aberrations in a few pathways, whereas others contain hundreds of coding somatic mutations. High-throughput RNA sequencing (RNA-seq) revealed that only approximately 36% of mutations are expressed. Using deep re-sequencing measurements of allelic abundance for 2,414 somatic mutations, we determine for the first time-to our knowledge-in an epithelial tumour subtype, the relative abundance of clonal frequencies among cases representative of the population. We show that TNBCs vary widely in their clonal frequencies at the time of diagnosis, with the basal subtype of TNBC showing more variation than non-basal TNBC. Although p53 (also known as TP53), PIK3CA and PTEN somatic mutations seem to be clonally dominant compared to other genes, in some tumours their clonal frequencies are incompatible with founder status. Mutations in cytoskeletal, cell shape and motility proteins occurred at lower clonal frequencies, suggesting that they occurred later during tumour progression. Taken together, our results show that understanding the biology and therapeutic responses of patients with TNBC will require the determination of individual tumour clonal genotypes.

Human mammary cancer progression model recapitulates methylation events associated with breast premalignancy.

INTRODUCTION: We have previously identified a rare subpopulation of variant human mammary epithelial cells (vHMEC) with repressed p16INK4A that exist in disease-free women yet display premalignant properties, suggesting that they have engaged the process of malignant transformation. In order to gain insight into the molecular alterations required for vHMEC to progress to malignancy, and to characterize the epigenetic events associated with early progression, we examined the effect of oncogenic stress on the behavior of these cells. METHODS: HMEC that express p16INK4A and vHMEC that do not, were transduced with constitutively active Ha-rasV12 and subsequently exposed to serum to determine whether signals from the cellular microenvironment could cooperate with ras to promote the malignant transformation of vHMEC. Epigenetic alterations were assessed using methylation-specific polymerase chain reaction (PCR). RESULTS: vHMEC expressing Ha-rasV12 (vHMEC-ras) bypassed the classic proliferative arrest that has been previously documented in normal fibroblasts following oncogenic stress, and that we also observe here in normal HMEC. Moreover, vHMEC-ras cells exhibited many additional alterations that are observed during progression to malignancy such as the generation of chromosomal abnormalities, upregulation of telomerase activity, immortalization following exposure to serum, and anchorage-independent growth, but they did not form tumors following orthotopic injection in vivo. Associated with their early progression to malignancy was an increase in the number of genes methylated, two of which (RASSF1A and SFRP1) were also methylated in other immortalized mammary cell lines as well as in breast cancer cells and tissues. CONCLUSIONS: We have characterized a mammary progression model that recapitulates molecular and methylation alterations observed in many breast cancers. Our data suggest that concomitant methylation of RASSF1A and SFRP1 marks an early event in mammary transformation and may thus have prognostic potential.

Sustained induction of epithelial to mesenchymal transition activates DNA methylation of genes silenced in basal-like breast cancers.

The active acquisition of epigenetic changes is a poorly understood but important process in development, differentiation, and disease. Our work has shown that repression of the p16/pRb pathway in human epithelial cells, a condition common to stem cells and many tumor cells, induces dynamic epigenetic remodeling resulting in the targeted methylation of a selected group of CpG islands. We hypothesized that cells in this epigenetically plastic state could be programmed by the microenvironment to acquire epigenetic changes associated with tumorigenesis. Here, we describe an in vitro model system where epigenetically plastic cells were placed in an environment that induced epithelial to mesenchymal transition (EMT) and led to a program of acquired de novo DNA methylation at targeted sites. In this model, we found that repression of E-cadherin transcription preceded the subsequent acquisition of methylated CpG sites. Furthermore, the induction of EMT was accompanied by de novo methylation of several other gene promoters, including those of the estrogen receptor and Twist. These data demonstrate that signals from the microenvironment can induce phenotypic and gene expression changes associated with targeted de novo epigenetic alterations important in tumor progression, and that these alterations occur through a deterministic, rather than stochastic, mechanism. Given the dynamic epigenetic reprogramming that occurs in these cells, DNA methylation profiles observed in human tumors may reflect the history of environmental exposures during the genesis of a tumor.

Tumor suppressor p16INK4A regulates polycomb-mediated DNA hypermethylation in human mammary epithelial cells.

Alterations in DNA methylation are important in cancer, but the acquisition of these alterations is poorly understood. Using an unbiased global screen for CpG island methylation events, we have identified a non-random pattern of DNA hypermethylation acquired in p16-repressed cells. Interestingly, this pattern included loci located upstream of a number of homeobox genes. Upon removal of p16(INK4A) activity in primary human mammary epithelial cells, polycomb repressors, EZH2 and SUZ12, are up-regulated and recruited to HOXA9, a locus expressed during normal breast development and epigenetically silenced in breast cancer. We demonstrate that at this targeted locus, the up-regulation of polycomb repressors is accompanied by the recruitment of DNA methyltransferases and the hypermethylation of DNA, an endpoint, which we show to be dependent on SUZ12 expression. These results demonstrate a causal role of p16(INK4A) disruption in modulating DNA hypermethylation, and identify a dynamic and active process whereby epigenetic modulation of gene expression is activated as an early event in breast tumor progression.