Epigenetic inactivation of the HOXA gene cluster in breast cancer

Analytical and therapeutic profiles of DNA methylation alterations in cancer; an overview of changes in chromatin arrangement and alterations in histone surfaces

DNA methylation is the most important epigenetic element that activates the inhibition of gene transcription and is included in the pathogenesis of all types of malignancies. Remarkably, the effectors of DNA methylation are DNMTs (DNA methyltransferases) that catalyze de novo or keep methylation of hemimethylated DNA after the DNA replication process. DNA methylation structures in cancer are altered, with three procedures by which DNA methylation helps cancer development which are including direct mutagenesis, hypomethylation of the cancer genome, and also focal hypermethylation of the promoters of TSGs (tumor suppressor genes). Conspicuously, DNA methylation, nucleosome remodeling, RNA-mediated targeting, and histone modification balance modulate many biological activities that are essential and indispensable to the genesis of cancer and also can impact many epigenetic changes including DNA methylation and histone modifications as well as adjusting of non-coding miRNAs expression in prevention and treatment of many cancers. Epigenetics points to heritable modifications in gene expression that do not comprise alterations in the DNA sequence. The nucleosome is the basic unit of chromatin, consisting of 147 base pairs (bp) of DNA bound around a histone octamer comprised of one H3/H4 tetramer and two H2A/H2B dimers. DNA methylation is preferentially distributed over nucleosome regions and is less increased over flanking nucleosome-depleted DNA, implying a connection between nucleosome positioning and DNA methylation. In carcinogenesis, aberrations in the epigenome may also include in the progression of drug resistance. In this report, we report the rudimentary notes behind these epigenetic signaling pathways and emphasize the proofs recommending that their misregulation can conclude in cancer. These findings in conjunction with the promising preclinical and clinical consequences observed with epigenetic drugs against chromatin regulators, confirm the important role of epigenetics in cancer therapy.

DNA methylation fingerprint of hepatocellular carcinoma from tissue and liquid biopsies

Hepatocellular carcinoma (HCC) is amongst the cancers with highest mortality rates and is the most common malignancy of the liver. Early detection is vital to provide the best treatment possible and liquid biopsies combined with analysis of circulating tumour DNA methylation show great promise as a non-invasive approach for early cancer diagnosis and monitoring with low false negative rates. To identify reliable diagnostic biomarkers of early HCC, we performed a systematic analysis of multiple hepatocellular studies and datasets comprising > 1500 genome-wide DNA methylation arrays, to define a methylation signature predictive of HCC in both tissue and cell-free DNA liquid biopsy samples. Our machine learning pipeline identified differentially methylated regions in HCC, some associated with transcriptional repression of genes related with cancer progression, that benchmarked positively against independent methylation signatures. Combining our signature of 38 DNA methylation regions, we derived a HCC detection score which confirmed the utility of our approach by identifying in an independent dataset 96% of HCC tissue samples with a precision of 98%, and most importantly successfully separated cfDNA of tumour samples from healthy controls. Notably, our risk score could identify cell-free DNA samples from patients with other tumours, including colorectal cancer. Taken together, we propose a comprehensive HCC DNA methylation fingerprint and an associated risk score for detection of HCC from tissue and liquid biopsies.

Comprehensive Landscape of HOXA2, HOXA9, and HOXA10 as Potential Biomarkers for Predicting Progression and Prognosis in Prostate Cancer

Prostate cancer (PCa) is recognized as a common malignancy in male patients. The homeobox A cluster (HOXA) family members have been confirmed to be implicated in the development of several types of tumors. However, the expression pattern and prognostic values of HOXA genes in PCa have not been investigated. In this study, we analyzed TCGA datasets and identified six HOXA family members which showed a dysregulated expression in PCa specimens compared with nontumor specimens. We also explored the potential mechanisms involved in the dysregulation of HOXA family members in PCa, and the results of Pearson's correlation revealed that most HOXA members were negatively related to the methylation degree. Moreover, we explored the prognostic values of HOXA family members and identified six survival-related HOXA members. Importantly, HOXA2, HOXA9, and HOXA10 were identified as critical PCa-related genes which were abnormally expressed in PCa and associated with clinical outcomes of PCa patients. Then, we explored the association between the above three genes and immune cell infiltration. We observed that the levels of HOXA2, HOXA9, and HOXA10 were associated with the levels of immune infiltration of several kinds of immune cells. Overall, our findings identified the potential values of the HOXA family for outcome prediction in PCa, which might facilitate personalized counselling and treatment in PCa.

Identifying and characterizing lincRNA genomic clusters reveals its cooperative functions in human cancer

Background

Emerging evidence has revealed that some long intergenic non-coding RNAs (lincRNAs) are likely to form clusters on the same chromosome, and lincRNA genomic clusters might play critical roles in the pathophysiological mechanism. However, the comprehensive investigation of lincRNA clustering is rarely studied, particularly the characterization of their functional significance across different cancer types.

Methods

In this study, we firstly constructed a computational method basing a sliding window approach for systematically identifying lincRNA genomic clusters. We then dissected these lincRNA genomic clusters to identify common characteristics in cooperative expression, conservation among divergent species, targeted miRNAs, and CNV frequency. Next, we performed comprehensive analyses in differentially-expressed patterns and overall survival outcomes for patients from The Cancer Genome Atlas (TCGA) and The Genotype-Tissue Expression (GTEx) across multiple cancer types. Finally, we explored the underlying mechanisms of lincRNA genomic clusters by functional enrichment analysis, pathway analysis, and drug-target interaction.

Results

We identified lincRNA genomic clusters according to the algorithm. Clustering lincRNAs tended to be co-expressed, highly conserved, targeted by more miRNAs, and with similar deletion and duplication frequency, suggesting that lincRNA genomic clusters may exert their effects by acting in combination. We further systematically explored conserved and cancer-specific lincRNA genomic clusters, indicating they were involved in some important mechanisms of disease occurrence through diverse approaches. Furthermore, lincRNA genomic clusters can serve as biomarkers with potential clinical significance and involve in specific pathological processes in the development of cancer. Moreover, a lincRNA genomic cluster named Cluster127 in DLK1-DIO3 imprinted locus was discovered, which contained MEG3, MEG8, MEG9, MIR381HG, LINC02285, AL132709.5, and AL132709.1. Further analysis indicated that Cluster127 may have the potential for predicting prognosis in cancer and could play their roles by participating in the regulation of PI3K-AKT signaling pathway.

Conclusions

Clarification of the lincRNA genomic clusters specific roles in human cancers could be beneficial for understanding the molecular pathogenesis of different cancer types.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-021-03179-5.

Nup93 and CTCF modulate spatiotemporal dynamics and function of the HOXA gene locus during differentiation

Nucleoporins regulate nuclear transport and are also involved in DNA damage, repair, cell cycle, chromatin organization, and gene expression. Here, we studied the role of nucleoporin Nup93 and the chromatin organizer CTCF in regulating HOXA expression during differentiation. ChIP sequencing revealed a significant overlap between Nup93 and CTCF peaks. Interestingly, Nup93 and CTCF are associated with the 3' and 5'HOXA genes respectively. Depletions of Nup93 and CTCF antagonistically modulate expression levels of 3'and 5'HOXA genes in undifferentiated NT2/D1 cells. Nup93 also regulates the localization of the HOXA gene locus, which disengages from the nuclear periphery upon Nup93 but not CTCF depletion, consistent with its upregulation. The dynamic association of Nup93 and CTCF with the HOXA locus during differentiation correlates with its spatial positioning and expression. While Nup93 tethers the HOXA locus to the nuclear periphery, CTCF potentially regulates looping of the HOXA gene cluster in a temporal manner. In summary, Nup93 and CTCF complement one another in modulating the spatiotemporal dynamics and function of the HOXA gene locus during differentiation.

Advances in mapping the epigenetic modifications of 5-methylcytosine (5mC), N6-methyladenine (6mA), and N4-methylcytosine (4mC)

DNA modification plays a pivotal role in regulating gene expression in cell development. As prevalent markers on DNA, 5-methylcytosine (5mC), N6-methyladenine (6mA), and N4-methylcytosine (4mC) can be recognized by specific methyltransferases, facilitating cellular defense and the versatile regulation of gene expression in eukaryotes and prokaryotes. Recent advances in DNA sequencing technology have permitted the positions of different modifications to be resolved at the genome-wide scale, which has led to the discovery of several novel insights into the complexity and functions of multiple methylations. In this review, we summarize differences in the various mapping approaches and discuss their pros and cons with respect to their relative read depths, speeds, and costs. We also discuss the development of future sequencing technologies and strategies for improving the detection resolution of current sequencing technologies. Lastly, we speculate on the potentially instrumental role that these sequencing technologies might play in future research.

Luminal A Breast Cancer Co-expression Network: Structural and Functional Alterations

Luminal A is the most common breast cancer molecular subtype in women worldwide. These tumors have characteristic yet heterogeneous alterations at the genomic and transcriptomic level. Gene co-expression networks (GCNs) have contributed to better characterize the cancerous phenotype. We have previously shown an imbalance in the proportion of intra-chromosomal (cis-) over inter-chromosomal (trans-) interactions when comparing cancer and healthy tissue GCNs. In particular, for breast cancer molecular subtypes (Luminal A included), the majority of high co-expression interactions connect gene-pairs in the same chromosome, a phenomenon that we have called loss of trans- co-expression. Despite this phenomenon has been described, the functional implication of this specific network topology has not been studied yet. To understand the biological role that communities of co-expressed genes may have, we constructed GCNs for healthy and Luminal A phenotypes. Network modules were obtained based on their connectivity patterns and they were classified according to their chromosomal homophily (proportion of cis-/trans- interactions). A functional overrepresentation analysis was performed on communities in both networks to observe the significantly enriched processes for each community. We also investigated possible mechanisms for which the loss of trans- co-expression emerges in cancer GCN. To this end we evaluated transcription factor binding sites, CTCF binding sites, differential gene expression and copy number alterations (CNAs) in the cancer GCN. We found that trans- communities in Luminal A present more significantly enriched categories than cis- ones. Processes, such as angiogenesis, cell proliferation, or cell adhesion were found in trans- modules. The differential expression analysis showed that FOXM1, CENPA, and CIITA transcription factors, exert a major regulatory role on their communities by regulating expression of their target genes in other chromosomes. Finally, identification of CNAs, displayed a high enrichment of deletion peaks in cis- communities. With this approach, we demonstrate that network topology determine, to at certain extent, the function in Luminal A breast cancer network. Furthermore, several mechanisms seem to be acting together to avoid trans- co-expression. Since this phenomenon has been observed in other cancer tissues, a remaining question is whether the loss of long distance co-expression is a novel hallmark of cancer.

Roles of the HOX Proteins in Cancer Invasion and Metastasis

Invasion and metastasis correspond to the foremost cause of cancer-related death, and the molecular networks behind these two processes are extremely complex and dependent on the intra- and extracellular conditions along with the prime of the premetastatic niche. Currently, several studies suggest an association between the levels of HOX genes expression and cancer cell invasion and metastasis, which favour the formation of novel tumour masses. The deregulation of HOX genes by HMGA2/TET1 signalling and the regulatory effect of noncoding RNAs generated by the HOX loci can also promote invasion and metastasis, interfering with the expression of HOX genes or other genes relevant to these processes. In this review, we present five molecular mechanisms of HOX deregulation by which the HOX clusters products may affect invasion and metastatic processes in solid tumours.

Profiling of immune related genes silenced in EBV-positive gastric carcinoma identified novel restriction factors of human gammaherpesviruses

EBV-associated gastric cancer (EBVaGC) is characterized by high frequency of DNA methylation. In this study, we investigated how epigenetic alteration of host genome contributes to pathogenesis of EBVaGC through the analysis of transcriptomic and epigenomic datasets from NIH TCGA (The Cancer Genome Atlas) consortium. We identified that immune related genes (IRGs) is a group of host genes preferentially silenced in EBV-positive gastric cancers through DNA hypermethylation. Further functional characterizations of selected IRGs reveal their novel antiviral activity against not only EBV but also KSHV. In particular, we showed that metallothionein-1 (MT1) and homeobox A (HOXA) gene clusters are down-regulated via EBV-driven DNA hypermethylation. Several MT1 isoforms suppress EBV lytic replication and release of progeny virions as well as KSHV lytic reactivation, suggesting functional redundancy of these genes. In addition, single HOXA10 isoform exerts antiviral activity against both EBV and KSHV. We also confirmed the antiviral effect of other dysregulated IRGs, such as IRAK2 and MAL, in scenario of EBV and KSHV lytic reactivation. Collectively, our results demonstrated that epigenetic silencing of IRGs is a viral strategy to escape immune surveillance and promote viral propagation, which is overall beneficial to viral oncogenesis of human gamma-herpesviruses (EBV and KSHV), considering that these IRGs possess antiviral activities against these oncoviruses.

HOTAIRM1 lncRNA is downregulated in clear cell renal cell carcinoma and inhibits the hypoxia pathway

HOXA Transcript Antisense RNA, Myeloid-Specific 1 (HOTAIRM1) is a conserved long non-coding RNA (lncRNA) involved in myeloid and neural differentiation that is deregulated in acute myeloid leukemia and other cancers. Previous studies focused on the nuclear unspliced HOTAIRM1 transcript, however cytoplasmic splice variants exist whose roles have remained unknown. Here, we report novel functions of HOTAIRM1 in the kidney. HOTAIRM1 transcripts are induced during renal lineage differentiation of embryonic stem cells and required for expression of specific renal differentiation genes. We show that the major HOTAIRM1 transcript in differentiated cells is the spliced cytoplasmic HM1-3 isoform and that HM1-3 is downregulated in >90% of clear cell renal cell carcinomas (ccRCCs). Knockdown of HM1-3 in renal cells deregulates hypoxia-responsive and angiogenic genes, including ANGPTL4. Furthermore, HOTAIRM1 transcripts are downregulated by hypoxia-mimetic stress and knockdown of the cytoplasmic HM1-3 isoform in normoxic cells post-transcriptionally induces Hypoxia-Inducible Factor 1α (HIF1α) protein, a key activator of ANGPTL4. Our results demonstrate the pervasive downregulation of the specific HOTAIRM1 cytoplasmic isoform HM1-3 in ccRCC and suggest possible roles of HOTAIRM1 in kidney differentiation and suppression of HIF1-dependent angiogenic pathways.

E3 ligase RCHY1 negatively regulates HDAC2

HDAC2, one of the class I histone deacetylase regulates epigenetic landscape through histone modification. Because HDAC2 is overexpressed in many cancers, cancer therapeutics against HDAC2 have been developed. Here we show novel mechanism of HDAC2 regulation by E3 ligase RCHY1. We found inverse correlation RCHY1 and HDAC2 levels in tumor tissue from six independent dataset using meta-analysis. Ectopic expression of RCHY1 decreased the level of HDAC2 from cancer cells including p53 wildtype, mutant and null cells. In addition, HDAC2 was increased by RCHY1 knockdown. RCHY1 directly interacts with HDAC2. Ectopic expression of wild type but not RING mutant RCHY1 increased HDAC2 levels. These data provide an evidence that RCHY1 negatively regulates HDAC2.

Genome-wide association study of germline variants and breast cancer-specific mortality

BACKGROUND

We examined the associations between germline variants and breast cancer mortality using a large meta-analysis of women of European ancestry.

METHODS

Meta-analyses included summary estimates based on Cox models of twelve datasets using ~10.4 million variants for 96,661 women with breast cancer and 7697 events (breast cancer-specific deaths). Oestrogen receptor (ER)-specific analyses were based on 64,171 ER-positive (4116) and 16,172 ER-negative (2125) patients. We evaluated the probability of a signal to be a true positive using the Bayesian false discovery probability (BFDP).

RESULTS

We did not find any variant associated with breast cancer-specific mortality at P < 5 × 10-8. For ER-positive disease, the most significantly associated variant was chr7:rs4717568 (BFDP = 7%, P = 1.28 × 10-7, hazard ratio [HR] = 0.88, 95% confidence interval [CI] = 0.84-0.92); the closest gene is AUTS2. For ER-negative disease, the most significant variant was chr7:rs67918676 (BFDP = 11%, P = 1.38 × 10-7, HR = 1.27, 95% CI = 1.16-1.39); located within a long intergenic non-coding RNA gene (AC004009.3), close to the HOXA gene cluster.

CONCLUSIONS

We uncovered germline variants on chromosome 7 at BFDP < 15% close to genes for which there is biological evidence related to breast cancer outcome. However, the paucity of variants associated with mortality at genome-wide significance underpins the challenge in providing genetic-based individualised prognostic information for breast cancer patients.

Dynamic chromosomal tuning of a novel GAU1 lncing driver at chr12p13.32 accelerates tumorigenesis

Aberrant chromatin transformation dysregulates gene expression and may be an important driver of tumorigenesis. However, the functional role of chromosomal dynamics in tumorigenesis remains to be elucidated. Here, using in vitro and in vivo experiments, we reveal a novel long noncoding (lncing) driver at chr12p13.3, in which a novel lncRNA GALNT8 Antisense Upstream 1 (GAU1) is initially activated by an open chromatin status, triggering recruitment of the transcription elongation factor TCEA1 at the oncogene GALNT8 promoter and cis-activates the expression of GALNT8. Analysis of The Cancer Genome Atlas (TCGA) clinical database revealed that the GAU1/GALNT8 driver serves as an important indicative biomarker, and targeted silencing of GAU1 via the HKP-encapsulated method exhibited therapeutic efficacy in orthotopic xenografts. Our study presents a novel oncogenetic mechanism in which aberrant tuning of the chromatin state at specific chromosomal loci exposes factor-binding sites, leading to recruitment of trans-factor and activation of oncogenetic driver, thereby provide a novel alternative concept of chromatin dynamics in tumorigenesis.

Homeobox A11 hypermethylation indicates unfavorable prognosis in breast cancer

Homeobox A11 (HOXA11) is one of the hypermethylated genes in breast cancer and its function in breast tumorigenesis remains elusive. In this study, we analyzed the methylation status of HOXA11 in 264 paired breast cancer and normal tissue as well as in matched serum samples by MethyLight assay. Further, the function of HOXA11 in breast tumorigenesis was analyzed by cell proliferation and migration assays. We found that HOXA11 was hypermethylated in cancer tissues (45.08%), especially in invasive ductal carcinomas (P<0.001), patients with a family history of cancer (P=0.033), cases with metastatic lymph nodes (P=0.004) and P53 positive group (P=0.017). Kaplan-Meier survival analysis and Cox regression analysis revealed that HOXA11 hypermethylation is an independent predictor of poor outcomes. The over expression of HOXA11 suppressed cell growth in MDA-MB-231, MCF7, SKBR3 and BT474 cells. In conclusion, the hypermethylation of HOXA11 is an independent prognostic biomarker in breast cancer. Additionally, HOXA11 can be a potential tumor suppressor.

Upregulated HOXA9 expression is associated with lymph node metastasis in colorectal cancer

Homeobox A (HOXA) cluster genes, members of the HOX family, perform an important role in normal organ development. It has previously been reported that HOXA gene expression in various types of cancer is associated with poor patient outcomes. However, the role of HOXA genes, as well as their expression, in colorectal cancers (CRC) remains unknown. Therefore, the present study investigated HOXA gene expression in patients with CRC and revealed that HOXA9 expression was significantly increased in tumor tissues compared with non-tumor tissues. Additionally, the functional role of HOXA9 was assessed by knocking down the HOXA9 gene in CRC cells and by evaluating cell growth. Regarding gene expression, cases with positive HOXA9 expression (as detected by immunohistochemical staining) were significantly associated with higher TNM stage and positive lymph node metastasis, although no association was observed between increased HOXA9 levels and the rate of overall survival in the present cohort. Regarding the functional role, HOXA9 expression was demonstrated to be upregulated in patients with CRC and was associated with lymph node metastasis.

Effect of 2,3',4,4',5-Pentachlorobiphenyl Exposure on Endometrial Receptivity and the Methylation of HOXA10

Polychlorinated biphenyls (PCBs) are one of the most common endocrine-disrupting chemicals and have obvious toxicity on human reproductive development. The aim of our study was to investigate the toxicity of chronic 2,3',4,4',5-pentachlorobiphenyl (PCB 118) exposure on embryo implantation and endometrial receptivity, with the possible mechanism of DNA methylation involved. Virgin CD-1 female mice (3 weeks old) were housed and orally treated with PCB 118 (0, 1, 10, 100 μg/kg) for a month. After mating with fertile males, the pregnant mice were killed on gestation day 4.5. Compared with the control group, implantation failures were observed in 1 μg/kg PCB 118- and 100 μg/kg PCB 118-treated groups. Abnormal endometrial morphology with open uterine lumens and densely compact stromal cells and poorly developed pinopodes were substantially in response to PCB 118 doses above, as well as the significant downregulation of implantation-associated genes (estrogen receptor 1, homeobox A10 [HOXA10], integrin subunit beta 3) and hypermethylation in the promoter region of HOXA10 further. It was confirmed that chronic exposure to PCB 118 produced an increased number of implantation failures in association with a defective uterine morphology during the implantation period. Alterations in methylation of HOXA10 could explain, at least in part, the mechanism of effects of PCB 118 exposure on the implantation process.

DNA methylation of extracellular matrix remodeling genes in children exposed to arsenic

Several novel mechanistic findings regarding to arsenic's pathogenesis has been reported and some of them suggest that the etiology of some arsenic induced diseases are due in part to heritable changes to the genome via epigenetic processes such as DNA methylation, histone maintenance, and mRNA expression. Recently, we reported that arsenic exposure during in utero and early life was associated with impairment in the lung function and abnormal receptor for advanced glycation endproducts (RAGE), matrix metalloproteinase-9 (MMP-9) and tissue inhibitor of matrix metalloproteinase-1 (TIMP-1) sputum levels. Based on our results and the reported arsenic impacts on DNA methylation, we designed this study in our cohort of children exposed in utero and early childhood to arsenic with the aim to associate DNA methylation of MMP9, TIMP1 and RAGE genes with its protein sputum levels and with urinary and toenail arsenic levels. The results disclosed hypermethylation in MMP9 promotor region in the most exposed children; and an increase in the RAGE sputum levels among children with the mid methylation level; there were also positive associations between MMP9 DNA methylation with arsenic toenail concentrations; RAGE DNA methylation with iAs, and %DMA; and finally between TIMP1 DNA methylation with the first arsenic methylation. A negative correlation between MMP9 sputum levels with its DNA methylation was registered. In conclusion, arsenic levels were positive associated with the DNA methylation of extracellular matrix remodeling genes;, which in turn could modifies the biological process in which they are involved causing or predisposing to lung diseases.

Tethered Oligonucleotide-Primed Sequencing, TOP-Seq: A High-Resolution Economical Approach for DNA Epigenome Profiling

Modification of CG dinucleotides in DNA is part of epigenetic regulation of gene function in vertebrates and is associated with complex human disease. Bisulfite sequencing permits high-resolution analysis of cytosine modification in mammalian genomes; however, its utility is often limited due to substantial cost. Here, we describe an alternative epigenome profiling approach, named TOP-seq, which is based on covalent tagging of individual unmodified CG sites followed by non-homologous priming of the DNA polymerase action at these sites to directly produce adjoining regions for their sequencing and precise genomic mapping. Pilot TOP-seq analyses of bacterial and human genomes showed a better agreement of TOP-seq with published bisulfite sequencing maps as compared to widely used MBD-seq and MRE-seq and permitted identification of long-range and gene-level differential methylation among human tissues and neuroblastoma cell types. Altogether, we propose an affordable single CG-resolution technique well suited for large-scale epigenome studies.

Mutual regulation of microRNAs and DNA methylation in human cancers

microRNAs (miRNAs) and DNA methylation are the 2 epigenetic modifications that have emerged in recent years as the most critical players in the regulation of gene expression. Compelling evidence has indicated the roles of miRNAs and DNA methylation in modulating cellular transformation and tumorigenesis. miRNAs act as negative regulators of gene expression and are involved in the regulation of both physiologic conditions and during diseases, such as cancer, inflammatory diseases, and psychiatric disorders, among others. Meanwhile, aberrant DNA methylation manifests in both global genome changes and in localized gene promoter changes, which influences the transcription of cancer genes. In this review, we described the mutual regulation of miRNAs and DNA methylation in human cancers. miRNAs regulate DNA methylation by targeting DNA methyltransferases or methylation-related proteins. On the other hand, both hyper- and hypo-methylation of miRNAs occur frequently in human cancers and represent a new level of complexity in gene regulation. Therefore, understanding the mechanisms underlying the mutual regulation of miRNAs and DNA methylation may provide helpful insights in the development of efficient therapeutic approaches.

HOXA repression is mediated by nucleoporin Nup93 assisted by its interactors Nup188 and Nup205

Background

The nuclear pore complex (NPC) mediates nuclear transport of RNA and proteins into and out of the nucleus. Certain nucleoporins have additional functions in chromatin organization and transcription regulation. Nup93 is a scaffold nucleoporin at the nuclear pore complex which is associated with human chromosomes 5, 7 and 16 and with the promoters of the HOXA gene as revealed by ChIP-on-chip studies using tiling microarrays for these chromosomes. However, the functional consequences of the association of Nup93 with HOXA is unknown.

Results

Here, we examined the association of Nup93 with the HOXA gene cluster and its consequences on HOXA gene expression in diploid colorectal cancer cells (DLD1). Nup93 showed a specific enrichment ~1 Kb upstream of the transcription start site of each of the HOXA1, HOXA3 and HOXA5 promoters, respectively. Furthermore, the association of Nup93 with HOXA was assisted by its interacting partners Nup188 and Nup205. The depletion of the Nup93 sub-complex significantly upregulated HOXA gene expression levels. However, expression levels of a control gene locus (GLCCI1) on human chromosome 7 were unaffected. Three-dimensional fluorescence in situ hybridization (3D-FISH) analyses revealed that the depletion of the Nup93 sub-complex (but not Nup98) disengages the HOXA gene locus from the nuclear periphery, suggesting a potential role for Nup93 in tethering and repressing the HOXA gene cluster. Consistently, Nup93 knockdown increased active histone marks (H3K9ac), decreased repressive histone marks (H3K27me3) on the HOXA1 promoter and increased transcription elongation marks (H3K36me3) within the HOXA1 gene. Moreover, the combined depletion of Nup93 and CTCF (a known organizer of HOXA gene cluster) but not Nup93 alone, significantly increased GLCCI1 gene expression levels. Taken together, this suggests a novel role for Nup93 and its interactors in repressing the HOXA gene cluster.

Conclusions

This study reveals that the nucleoporin Nup93 assisted by its interactors Nup188 and Nup205 mediates the repression of HOXA gene expression.

Electronic supplementary material

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

Genome-wide association studies and epigenome-wide association studies go together in cancer control

Completion of the human genome a decade ago laid the foundation for: using genetic information in assessing risk to identify individuals and populations that are likely to develop cancer, and designing treatments based on a person's genetic profiling (precision medicine). Genome-wide association studies (GWAS) completed during the past few years have identified risk-associated single nucleotide polymorphisms that can be used as screening tools in epidemiologic studies of a variety of tumor types. This led to the conduct of epigenome-wide association studies (EWAS). This article discusses the current status, challenges and research opportunities in GWAS and EWAS. Information gained from GWAS and EWAS has potential applications in cancer control and treatment.

Systematically Prioritizing Functional Differentially Methylated Regions (fDMRs) by Integrating Multi-omics Data in Colorectal Cancer

While genome-wide differential DNA methylation regions (DMRs) have been extensively identified, the comprehensive prioritization of their functional importance is still poorly explored. Here, we aggregated multiple data resources rooted in the genome, epigenome and transcriptome to systematically prioritize functional DMRs (fDMRs) in colorectal cancer (CRC). As demonstrated, the top-ranked fDMRs from all of the data resources showed a strong enrichment for known methylated genes. Additionally, we analyzed those top 5% DMR-coupled coding genes using functional enrichment, which resulted in significant disease-related biological functions in contrast to the tail 5% genes. To further confirm the functional importance of the top-ranked fDMRs, we applied chromatin modification alterations of CRC cell lines to characterize their functional regulation. Specifically, we extended the utility of the top-ranked DMR-coupled genes to serve as classification and survival biomarkers, which showed a robust performance across diverse independent data sets. Collectively, our results established an integrative framework to prioritize fDMRs, which could help characterize aberrant DNA methylation-induced potential mechanisms underlying tumorigenesis and uncover epigenome-based biomarkers for clinical diagnosis and prognosis.

Estrogen-induced chromatin decondensation and nuclear re-organization linked to regional epigenetic regulation in breast cancer

BACKGROUND

Epigenetic changes are being increasingly recognized as a prominent feature of cancer. This occurs not only at individual genes, but also over larger chromosomal domains. To investigate this, we set out to identify large chromosomal domains of epigenetic dysregulation in breast cancers.

RESULTS

We identify large regions of coordinate down-regulation of gene expression, and other regions of coordinate activation, in breast cancers and show that these regions are linked to tumor subtype. In particular we show that a group of coordinately regulated regions are expressed in luminal, estrogen-receptor positive breast tumors and cell lines. For one of these regions of coordinate gene activation, we show that regional epigenetic regulation is accompanied by visible unfolding of large-scale chromatin structure and a repositioning of the region within the nucleus. In MCF7 cells, we show that this depends on the presence of estrogen.

CONCLUSIONS

Our data suggest that the liganded estrogen receptor is linked to long-range changes in higher-order chromatin organization and epigenetic dysregulation in cancer. This may suggest that as well as drugs targeting histone modifications, it will be valuable to investigate the inhibition of protein complexes involved in chromatin folding in cancer cells.

Genome-wide DNA methylation analysis reveals estrogen-mediated epigenetic repression of metallothionein-1 gene cluster in breast cancer

Background

Recent genome-wide analysis has shown that DNA methylation spans long stretches of chromosome regions consisting of clusters of contiguous CpG islands or gene families. Hypermethylation of various gene clusters has been reported in many types of cancer. In this study, we conducted methyl-binding domain capture (MBDCap) sequencing (MBD-seq) analysis on a breast cancer cohort consisting of 77 patients and 10 normal controls, as well as a panel of 38 breast cancer cell lines.

Results

Bioinformatics analysis determined seven gene clusters with a significant difference in overall survival (OS) and further revealed a distinct feature that the conservation of a large gene cluster (approximately 70 kb) metallothionein-1 (MT1) among 45 species is much lower than the average of all RefSeq genes. Furthermore, we found that DNA methylation is an important epigenetic regulator contributing to gene repression of MT1 gene cluster in both ERα positive (ERα+) and ERα negative (ERα−) breast tumors. In silico analysis revealed much lower gene expression of this cluster in The Cancer Genome Atlas (TCGA) cohort for ERα + tumors. To further investigate the role of estrogen, we conducted 17β-estradiol (E2) and demethylating agent 5-aza-2′-deoxycytidine (DAC) treatment in various breast cancer cell types. Cell proliferation and invasion assays suggested MT1F and MT1M may play an anti-oncogenic role in breast cancer.

Conclusions

Our data suggests that DNA methylation in large contiguous gene clusters can be potential prognostic markers of breast cancer. Further investigation of these clusters revealed that estrogen mediates epigenetic repression of MT1 cluster in ERα + breast cancer cell lines. In all, our studies identify thousands of breast tumor hypermethylated regions for the first time, in particular, discovering seven large contiguous hypermethylated gene clusters.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-015-0045-9) contains supplementary material, which is available to authorized users.

SUZ12 promotes gastric cancer cell proliferation and metastasis by regulating KLF2 and E-cadherin

SUZ12 is a core component of the polycomb repressive complex 2 (PRC2), which could silence gene transcription by generating trimethylation on lysine 27 residue of histone H3 (H3K27Me3). Meanwhile, SUZ12 has been found to be overexpressed in multiple cancers; however, the clinical significance and molecular mechanisms of SUZ12 controlling gastric cancer cell proliferation and metastasis are unclear. In this study, we found that SUZ12 expression was significantly increased in 64 gastric tumor tissues compared with normal tissues. Additionally, SUZ12 expression was associated with pathological stage, metastasis distance, and shorter overall survival of gastric cancer patients. Knockdown of SUZ12 expression impaired cell proliferation and invasion in vitro, leading to the inhibition of metastasis in vivo. Upregulation of SUZ12 was found to play a key role in gastric cancer cell proliferation and metastasis through the regulation of EMT and KLF2 expression.

PRC2-independent chromatin compaction and transcriptional repression in cancer

The silencing of large chromosomal regions by epigenetic mechanisms has been reported to occur frequently in cancer. Epigenetic marks, such as histone methylation and acetylation, are altered at these loci. However, the mechanisms of formation of such aberrant gene clusters remain largely unknown. Here, we show that, in cancer cells, the epigenetic remodeling of chromatin into hypoacetylated domains covered with histone H3K27 trimethylation is paralleled by changes in higher-order chromatin structures. Using fluorescence in situ hybridization, we demonstrate that regional epigenetic silencing corresponds to the establishment of compact chromatin domains. We show that gene repression is tightly correlated to the state of chromatin compaction and not to the levels of H3K27me3-its removal through the knockdown of EZH2 does not induce significant gene expression nor chromatin decompaction. Moreover, transcription can occur with intact high-H3K27me3 levels; treatment with histone deacetylase inhibitors can relieve chromatin compaction and gene repression, without altering H3K27me3 levels. Our findings imply that compaction and subsequent repression of large chromatin domains are not direct consequences of PRC2 deregulation in cancer cells. By challenging the role of EZH2 in aberrant gene silencing in cancer, these findings have therapeutical implications, notably for the choice of epigenetic drugs for tumors with multiple regional epigenetic alterations.

BAZ2A (TIP5) is involved in epigenetic alterations in prostate cancer and its overexpression predicts disease recurrence

Prostate cancer is driven by a combination of genetic and/or epigenetic alterations. Epigenetic alterations are frequently observed in all human cancers, yet how aberrant epigenetic signatures are established is poorly understood. Here we show that the gene encoding BAZ2A (TIP5), a factor previously implicated in epigenetic rRNA gene silencing, is overexpressed in prostate cancer and is paradoxically involved in maintaining prostate cancer cell growth, a feature specific to cancer cells. BAZ2A regulates numerous protein-coding genes and directly interacts with EZH2 to maintain epigenetic silencing at genes repressed in metastasis. BAZ2A overexpression is tightly associated with a molecular subtype displaying a CpG island methylator phenotype (CIMP). Finally, high BAZ2A levels serve as an independent predictor of biochemical recurrence in a cohort of 7,682 individuals with prostate cancer. This work identifies a new aberrant role for the epigenetic regulator BAZ2A, which can also serve as a useful marker for metastatic potential in prostate cancer.

Exploring breast carcinogenesis through integrative genomics and epigenomics analyses

There have been many DNA methylation studies on breast cancer which showed various methylation patterns involving tumour suppressor genes and oncogenes but only a few of those studies link the methylation data with gene expression. More data are required especially from the Asian region and to analyse how the epigenome data correlate with the transcriptome. DNA methylation profiling was carried out on 76 fresh frozen primary breast tumour tissues and 25 adjacent non-cancerous breast tissues using the Illumina Infinium(®) HumanMethylation27 BeadChip. Validation of methylation results was performed on 7 genes using either MS-MLPA or MS-qPCR. Gene expression profiling was done on 15 breast tumours and 5 adjacent non-cancerous breast tissues using the Affymetrix GeneChip(®) Human Gene 1.0 ST array. The overlapping genes between DNA methylation and gene expression datasets were further mapped to the KEGG database to identify the molecular pathways that linked these genes together. Supervised hierarchical cluster analysis revealed 1,389 hypermethylated CpG sites and 22 hypomethylated CpG sites in cancer compared to the normal samples. Gene expression microarray analysis using a fold-change of at least 1.5 and a false discovery rate (FDR) at p>0.05 identified 404 upregulated and 463 downregulated genes in cancer samples. Integration of both datasets identified 51 genes with hypermethylation with low expression (negative association) and 13 genes with hypermethylation with high expression (positive association). Most of the overlapping genes belong to the focal adhesion and extracellular matrix-receptor interaction that play important roles in breast carcinogenesis. The present study displayed the value of using multiple datasets in the same set of tissues and how the integrative analysis can create a list of well-focused genes as well as to show the correlation between epigenetic changes and gene expression. These gene signatures can help us understand the epigenetic regulation of gene expression and could be potential targets for therapeutic intervention in the future.

Aldehyde dehydrogenase 1A3 influences breast cancer progression via differential retinoic acid signaling

Aldehyde dehydrogenase (ALDH) 1A enzymes produce retinoic acid (RA), a transcription induction molecule. To investigate if ALDH1A1 or ALDH1A3-mediated RA signaling has an active role in breast cancer tumorigenesis, we performed gene expression and tumor xenograft studies. Analysis of breast patient tumors revealed that high levels of ALDH1A3 correlated with expression of RA-inducible genes with retinoic acid response elements (RAREs), poorer patient survival and triple-negative breast cancers. This suggests a potential link between ALDH1A3 expression and RA signaling especially in aggressive and/or triple-negative breast cancers. In MDA-MB-231, MDA-MB-468 and MDA-MB-435 cells, ALDH1A3 and RA increased expression of RA-inducible genes. Interestingly, ALDH1A3 had opposing effects in tumor xenografts, increasing tumor growth and metastasis of MDA-MB-231 and MDA-MB-435 cells, but decreasing tumor growth of MDA-MB-468 cells. Exogenous RA replaced ALDH1A3 in inducing the same opposing tumor growth and metastasis effects, suggesting that ALDH1A3 mediates these effects by promoting RA signaling. Genome expression analysis revealed that ALDH1A3 induced largely divergent gene expression in MDA-MB-231 and MDA-MB-468 cells which likely resulted in the opposing tumor growth effects. Treatment with DNA methylation inhibitor 5-aza-2'deoxycytidine restored uniform RA-inducibility of RARE-containing HOXA1 and MUC4 in MDA-MB-231 and MDA-MB-468 cells, suggesting that differences in epigenetic modifications contribute to differential ALDH1A3/RA-induced gene expression in breast cancer. In summary, ALDH1A3 induces differential RA signaling in breast cancer cells which affects the rate of breast cancer progression.

Genomic and transcriptomic analyses match medulloblastoma mouse models to their human counterparts

Medulloblastoma is a malignant embryonal brain tumor with highly variable outcome. In order to study the biology of this tumor and to perform preclinical treatment studies, a lot of effort has been put into the generation of appropriate mouse models. The usage of these models, however, has become debatable with the advances in human medulloblastoma subgrouping. This study brings together multiple relevant mouse models and matches genetic alterations and gene expression data of 140 murine tumors with 423 human medulloblastomas in a global way. Using AGDEX analysis and k-means clustering, we show that the Blbp-cre::Ctnnb1(ex3)(Fl/+)Trp53 (Fl/Fl) mouse model fits well to human WNT medulloblastoma, and that, among various Myc- or Mycn-based mouse medulloblastomas, tumors in Glt1-tTA::TRE-MYCN/Luc mice proved to be most specific for human group 3 medulloblastoma. None of the analyzed models displayed a significant match to group 4 tumors. Intriguingly, mice with Ptch1 or Smo mutations selectively modeled SHH medulloblastomas of adulthood, although such mutations occur in all human age groups. We therefore suggest that the infantile or adult gene expression pattern of SHH MBs are not solely determined by specific mutations. This is supported by the observation that human medulloblastomas with PTCH1 mutations displayed more similarities to PTCH1 wild-type tumors of the same age group than to PTCH1-mutated tumors of the other age group. Together, we provide novel insights into previously unrecognized specificity of distinct models and suggest these findings as a solid basis to choose the appropriate model for preclinical studies on medulloblastoma.

A signature predicting poor prognosis in gastric and ovarian cancer represents a coordinated macrophage and stromal response

PURPOSE

Gene-expression profiling has revolutionized the way we think about cancer and confers the ability to observe the synchronous expression of thousands of genes. The use of putative genome-level expression profiles has allowed biologists to observe the complex interactions of genes that constitute recognized biologic pathways. We used gastric and ovarian datasets to identify gene-expression signatures and determine any functional significance.

EXPERIMENTAL DESIGN

Microarray data of 94-tumor and 45-benign samples derived from patients with gastric cancer were interrogated using Hierarchical Ordered Partitioning and Collapsing Hybrid analysis identifying clusters of coexpressed genes. Clusters were further characterized with respect to biologic significance, gene ontology, and ability to discriminate between normal and tumor tissue. Tumor tissues were separated into epithelial and stromal compartments and immunohistochemical analysis performed to further elucidate specific cell lineages expressing genes contained in the signature.

RESULTS

We identified a "stromal-response" expression signature, highly enriched for inflammatory, extracellular matrix, cytokine, and growth factor proteins. The majority of genes in the signature are expressed in the tumor-associated stroma but were absent in associated premalignant conditions. In gastric cancer, this module almost perfectly differentiates tumor from nonmalignant gastric tissue and hence can be regarded as a highly tumor-specific gene-expression signature.

CONCLUSIONS

We show that these genes are consistently coexpressed across a range of independent gastric datasets as well as other cancer types suggesting a conserved functional role in cancer. In addition, we show that this signature can be a surrogate marker for M2 macrophage activity and has significant prognostic implications in gastric and ovarian high-grade serous cancer.

A key role for EZH2 in epigenetic silencing of HOX genes in mantle cell lymphoma

The chromatin modifier EZH2 is overexpressed and associated with inferior outcome in mantle cell lymphoma (MCL). Recently, we demonstrated preferential DNA methylation of HOX genes in MCL compared with chronic lymphocytic leukemia (CLL), despite these genes not being expressed in either entity. Since EZH2 has been shown to regulate HOX gene expression, to gain further insight into its possible role in differential silencing of HOX genes in MCL vs. CLL, we performed detailed epigenetic characterization using representative cell lines and primary samples. We observed significant overexpression of EZH2 in MCL vs. CLL. Chromatin immune precipitation (ChIP) assays revealed that EZH2 catalyzed repressive H3 lysine 27 trimethylation (H3K27me3), which was sufficient to silence HOX genes in CLL, whereas in MCL H3K27me3 is accompanied by DNA methylation for a more stable repression. More importantly, hypermethylation of the HOX genes in MCL resulted from EZH2 overexpression and subsequent recruitment of the DNA methylation machinery onto HOX gene promoters. The importance of EZH2 upregulation in this process was further underscored by siRNA transfection and EZH2 inhibitor experiments. Altogether, these observations implicate EZH2 in the long-term silencing of HOX genes in MCL, and allude to its potential as a therapeutic target with clinical impact.

Real-time methylomic aberrations during initiation and progression of induced human mammary epithelial cell tumorigenesis

AIM

Neoplastic transformation provides one of the few existing opportunities to analyze molecular changes in real time during the initiation and progression of breast cancer.

MATERIALS & METHODS

Human mammary epithelial cells underwent neoplastic reprogramming, generating one line of semitransformed, premalignant cells and two separate, temporal lines of fully transformed human mammary epithelial cells (THMECs). An Illumina Infinium HumanMethylation27 BeadChip was used to analyze DNA methylation alterations in 27,578 CpG loci at three consecutive time points over an 80-day (d) transformation period.

RESULTS

The mean β value for semitransformed human mammary epithelial cells CpG loci (0.245) was much greater than for either THMEC-40d (0.055) or THMEC-80d (0.066), indicating a large loss of methylation after neoplastic induction. In addition, 54% of CpG loci were hypermethylated during the THMEC-40d to THMEC-80d transition. We observed that the CpG loci exhibiting DNA methylation changes during early oncogenesis were enriched for biological functions like cellular movement; this was distinctly different than in the later, more progressive stages of the transformation process enriched for processes involving differentiation.

CONCLUSION

The timing of major methylomic changes may be important in directing the cell toward a more cancerous phenotype. In addition, gene-specific hypermethylation appears to silence developmentally related genes, leading to dedifferentiation.

HMGA2/TET1/HOXA9 signaling pathway regulates breast cancer growth and metastasis

The ten-eleven translocation (TET) family of methylcytosine dioxygenases initiates demethylation of DNA and is associated with tumorigenesis in many cancers; however, the mechanism is mostly unknown. Here we identify upstream activators and downstream effectors of TET1 in breast cancer using human breast cancer cells and a genetically engineered mouse model. We show that depleting the architectural transcription factor high mobility group AT-hook 2 (HMGA2) induces TET1. TET1 binds and demethylates its own promoter and the promoter of homeobox A (HOXA) genes, enhancing its own expression and stimulating expression of HOXA genes including HOXA7 and HOXA9. Both TET1 and HOXA9 suppress breast tumor growth and metastasis in mouse xenografts. The genes comprising the HMGA2-TET1-HOXA9 pathway are coordinately regulated in breast cancer and together encompass a prognostic signature for patient survival. These results implicate the HMGA2-TET1-HOX signaling pathway in the epigenetic regulation of human breast cancer and highlight the importance of targeting methylation in specific subpopulations as a potential therapeutic strategy.

miRNA gene promoters are frequent targets of aberrant DNA methylation in human breast cancer

miRNAs are important regulators of gene expression that are frequently deregulated in cancer, with aberrant DNA methylation being an epigenetic mechanism involved in this process. We previously identified miRNA promoter regions active in normal mammary cell types and here we analyzed which of these promoters are targets of aberrant DNA methylation in human breast cancer cell lines and breast tumor specimens. Using 5-methylcytosine immunoprecipitation coupled to miRNA tiling microarray hybridization, we performed comprehensive evaluation of DNA methylation of miRNA gene promoters in breast cancer. We found almost one third (55/167) of miRNA promoters were targets for aberrant methylation in breast cancer cell lines. Breast tumor specimens displayed DNA methylation of majority of these miRNA promoters, indicating that these changes in DNA methylation might be clinically relevant. Aberrantly methylated miRNA promoters were, similar to protein coding genes, enriched for promoters targeted by polycomb in normal cells. Detailed analysis of selected miRNA promoters revealed decreased expression of miRNA linked to increased promoter methylation for mir-31, mir-130a, let-7a-3/let-7b, mir-155, mir-137 and mir-34b/mir-34c genes. The proportion of miRNA promoters we found aberrantly methylated in breast cancer is several fold larger than that observed for protein coding genes, indicating an important role of DNA methylation in miRNA deregulation in cancer.

Clustered protocadherins and neuronal diversity

Neuronal diversity is a fundamental requirement for complex neuronal networks and brain function. The clustered protocadherin (Pcdh) family possesses several characteristic features that are important for the molecular basis of neuronal diversity. Clustered Pcdhs are expressed predominantly in the central nervous system, in neurites, growth cones, and synapses. They consist of about 60 isoforms, and their expression is stochastically and combinatorially regulated in individual neurons. The multiple clustered Pcdhs expressed in individual neurons form heteromultimeric protein complexes that exhibit homophilic adhesion properties. Theoretically, the clustered Pcdhs could generate more than 3×10(10) possible variations in each neuron and 12,720 types of cis-tetramers per neuron. The clustered Pcdhs are important for normal neuronal development. The clustered Pcdh genes have also attracted attention as a target for epigenetic regulation.

Epigenetic changes during cell transformation

Malignant cancer emerges from normal healthy cells in a multistep -process that involves both genetic and epigenetic lesions. Both genetic and environmental inputs participate in driving the epigenetic changes that occur during human carcinogenesis. The pathologic changes seen in DNA methylation and histone posttranslational modifications are complex, deeply intertwined, and act in concert to produce malignant transformation. To better understand the causes and consequences of the pathoepigenetic changes in cancer formation, a variety of experimentally tractable human cell line model systems that accurately reflect the molecular alterations seen in the clinical disease have been developed. Results from studies using these cell line model systems suggest that early critical epigenetic events occur in a stepwise fashion prior to cell immortalization. These epigenetic steps coincide with the cell's transition through well-defined cell proliferation barriers of stasis and telomere dysfunction. Following cell immortalization, stressors, such as environmental toxicants, can induce malignant transformation in a process in which the epigenetic changes occur in a smoother progressive fashion, in contrast to the stark stepwise epigenetic changes seen prior to cell immortalization. It is hoped that developing a clearer understanding of the identity, timing, and consequences of these epigenetic lesions will prove useful in future clinical applications that range from early disease detection to therapeutic intervention in malignant cancer.

Cell-type specific DNA methylation patterns define human breast cellular identity

DNA methylation plays a role in a variety of biological processes including embryonic development, imprinting, X-chromosome inactivation, and stem cell differentiation. Tissue specific differential methylation has also been well characterized. We sought to extend these studies to create a map of differential DNA methylation between different cell types derived from a single tissue. Using three pairs of isogenic human mammary epithelial and fibroblast cells, promoter region DNA methylation was characterized using MeDIP coupled to microarray analysis. Comparison of DNA methylation between these cell types revealed nearly three thousand cell-type specific differentially methylated regions (ctDMRs). MassARRAY was performed upon 87 ctDMRs to confirm and quantify differential DNA methylation. Each of the examined regions exhibited statistically significant differences ranging from 10-70%. Gene ontology analysis revealed the overrepresentation of many transcription factors involved in developmental processes. Additionally, we have shown that ctDMRs are associated with histone related epigenetic marks and are often aberrantly methylated in breast cancer. Overall, our data suggest that there are thousands of ctDMRs which consistently exhibit differential DNA methylation and may underlie cell type specificity in human breast tissue. In addition, we describe the pathways affected by these differences and provide insight into the molecular mechanisms and physiological overlap between normal cellular differentiation and breast carcinogenesis.

HOX gene methylation status analysis in patients with hereditary breast cancer

Cancer development is related not only to genetic alterations but also to aberrant epigenetic changes that could lead to heritable gene patterns critical for neoplastic initiation and progression. Knowledge of epigenetic regulation in cancer cells is useful for both the understanding of carcinogenesis and for the possibility of using epigenetic drugs. HOX genes deregulation have a crucial role in oncogenesis process and tumor suppression. In this report, the methylation of HOXA1, HOXA9, HOXA10, HOXB13, HNF1B, OTX1, TLX1 genes have been analyzed in patients with hereditary breast cancer. This is the first study analyzing BRCA mutational status of patients with respect to methylation of HOX genes. HOXA10 has been found to be methylated in all patients analyzed but never in healthy subjects. With respect to clinical pathological information, hypermethylation of all studied genes, with the exception of OTX1, was significantly associated with absence of HER2 neu expression (P<0.05). Moreover, hypermethylation of HOXB13, HOXA10 and HOXA1 was associated with a high proliferation index (Mib1≥10%, P<0.05) and hypermethylation of HOXB13 and HOXA10 also with high expression of estrogen and progesterone receptors. These preliminary data suggest a possible involvement of HOX genes in familial breast cancer as marker helpful to identify high-risk patients.

Agglomerates of aberrant DNA methylation are associated with toxicant-induced malignant transformation

Epigenetic dysfunction is a known contributor in carcinogenesis, and is emerging as a mechanism involved in toxicant-induced malignant transformation for environmental carcinogens such as arsenicals or cadmium. In addition to aberrant DNA methylation of single genes, another manifestation of epigenetic dysfunction in cancer is agglomerative DNA methylation, which can participate in long-range epigenetic silencing that targets many neighboring genes and has been shown to occur in several types of clinical cancers. Using in vitro model systems of toxicant-induced malignant transformation, we found hundreds of aberrant DNA methylation events that emerge during malignant transformation, some of which occur in an agglomerative fashion. In an arsenite-transformed prostate epithelial cell line, the protocadherin (PCDH), HOXC and HOXD gene family clusters are targeted for agglomerative DNA methylation. The agglomerative DNA methylation changes induced by arsenicals appear to be common and clinically relevant events, since they occur in other human cancer cell lines and models of malignant transformation, as well as clinical cancer specimens. Aberrant DNA methylation in general occurred more often within histone H3 lysine-27 trimethylation stem cell domains. We found a striking association between enrichment of histone H3 lysine-9 trimethylation stem cell domains and toxicant-induced agglomerative DNA methylation, suggesting these epigenetic modifications may become aberrantly linked during malignant transformation. In summary, we found an association between toxicant-induced malignant transformation and agglomerative DNA methylation, which lends further support to the hypothesis that epigenetic dysfunction plays an important role in toxicant-induced malignant transformation.

FAM83B mediates EGFR- and RAS-driven oncogenic transformation

Aberrant regulation of growth signaling is a hallmark of cancer development that often occurs through the constitutive activation of growth factor receptors or their downstream effectors. Using validation-based insertional mutagenesis (VBIM), we identified family with sequence similarity 83, member B (FAM83B), based on its ability to substitute for RAS in the transformation of immortalized human mammary epithelial cells (HMECs). We found that FAM83B coprecipitated with a downstream effector of RAS, CRAF. Binding of FAM83B with CRAF disrupted CRAF/14-3-3 interactions and increased CRAF membrane localization, resulting in elevated MAPK and mammalian target of rapamycin (mTOR) signaling. Ablation of FAM83B inhibited the proliferation and malignant phenotype of tumor-derived cells or RAS-transformed HMECs, implicating FAM83B as a key intermediary in EGFR/RAS/MAPK signaling. Analysis of human tumor specimens revealed that FAM83B expression was significantly elevated in cancer and was associated with specific cancer subtypes, increased tumor grade, and decreased overall survival. Cumulatively, these results suggest that FAM83B is an oncogene and potentially represents a new target for therapeutic intervention.

Hydroxylation of 5-methylcytosine by TET2 maintains the active state of the mammalian HOXA cluster

Differentiation is accompanied by extensive epigenomic reprogramming, leading to the repression of stemness factors and the transcriptional maintenance of activated lineage-specific genes. Here we use the mammalian Hoxa cluster of developmental genes as a model system to follow changes in DNA modification patterns during retinoic acid-induced differentiation. We find the inactive cluster to be marked by defined patterns of 5-methylcytosine (5mC). Upon the induction of differentiation, the active anterior part of the cluster becomes increasingly enriched in 5-hydroxymethylcytosine (5hmC), following closely the colinear activation pattern of the gene array, which is paralleled by the reduction of 5mC. Depletion of the 5hmC generating dioxygenase Tet2 impairs the maintenance of Hoxa activity and partially restores 5mC levels. Our results indicate that gene-specific 5mC-5hmC conversion by Tet2 is crucial for the maintenance of active chromatin states at lineage-specific loci.

Evidence of epigenetic regulation of the tumor suppressor gene cluster flanking RASSF1 in breast cancer cell lines

Epigenetic mechanisms are frequently deregulated in cancer cells and can lead to the silencing of genes with tumor suppressor activities. The isoform A of the Ras-association domain family member 1 (RASSF1A) gene is one of the most frequently silenced transcripts in human tumors, however, few studies have simultaneously investigated epigenetic abnormalities associated with the 3p21.3 tumor suppressor gene cluster flanking RASSF1 (i.e., SEMA3B, HYAL3, HYAL2, HYAL1, TUSC2, RASSF1, ZMYND10, NPRL2, TMEM115, and CACNA2D2). This study aimed to investigate the role of epigenetic changes to these genes in seventeen breast cancer cell lines and in three non-tumorigenic epithelial breast cell lines (184A1, 184B5, and MCF 10A) and to evaluate the effect on gene expression of treatment with the demethylating agent 5-Aza-2'-deoxycytidine and/or Trichostatin A (TSA), a histone deacetylase inhibitor. We report that, although the RASSF1A isoform was determined to be epigenetically silenced in 15 of the 17 breast cancer cell lines, all the cell lines expressed the RASSF1C isoform. Five breast cancer cell lines overexpressed RASSF1C, when compared to the normal epithelial cell line 184A1. Furthermore, the genes HYAL1 and CACNA2D2 were significantly overexpressed after the treatments. After the combinated treatment, RASSF1A re-expression was accompanied by an increase in expression levels of the flanking genes. The Spearman's correlation coefficient indicated a positive co-regulation of the following gene pairs: RASSF1 and TUSC2 (r=0.64, p=0.002), RASSF1 and ZMYND10 (r=0.58, p=0.07), RASSF1 and NPRL2 (r=0.48, p=0.03), ZMYND10 and NPRL2 (r=0.71; p=0,0004), and NPRL2 and TMEM115 (r=0.66, p=0.001). Interestingly, the genes TUSC2, NPRL2 and TMEM115 were found to be unmethylated in each of the untreated cell lines. Chromatin immunoprecipitation using antibodies against the acetylated and trimethylated lysine 9 of histone H3 demonstrated low levels of histone methylation in these genes, which are located closest to RASSF1. These results provide evidence that epigenetic repression is involved in the down-regulation of multiple genes at 3p21.3 in breast cancer cells.

Rare and frequent promoter methylation, respectively, of TSHZ2 and 3 genes that are both downregulated in expression in breast and prostate cancers

Background

Neoplastic cells harbor both hypomethylated and hypermethylated regions of DNA. Whereas hypomethylation is found mainly in repeat sequences, regional hypermethylation has been linked to the transcriptional silencing of certain tumor suppressor genes. We attempted to search for candidate genes involved in breast/prostate carcinogenesis, using the criteria that they should be expressed in primary cultures of normal breast/prostate epithelial cells but are frequently downregulated in breast/prostate cancer cell lines and that their promoters are hypermethylated.

Methodology/Principal Findings

We identified several dozens of candidates among 194 homeobox and related genes using Systematic Multiplex RT-PCR and among 23,000 known genes and 23,000 other expressed sequences in the human genome by DNA microarray hybridization. An additional examination, by real-time qRT-PCR of clinical specimens of breast cancer, further narrowed the list of the candidates. Among them, the most frequently downregulated genes in tumors were NP_775756 and ZNF537, from the homeobox gene search and the genome-wide search, respectively. To our surprise, we later discovered that these genes belong to the same gene family, the 3-member Teashirt family, bearing the new names of TSHZ2 and TSHZ3. We subsequently determined the methylation status of their gene promoters. The TSHZ3 gene promoter was found to be methylated in all the breast/prostate cancer cell lines and some of the breast cancer clinical specimens analyzed. The TSHZ2 gene promoter, on the other hand, was unmethylated except for the MDA-MB-231 breast cancer cell line. The TSHZ1 gene was always expressed, and its promoter was unmethylated in all cases.

Conclusions/Significance

TSHZ2 and TSHZ3 genes turned out to be the most interesting candidates for novel tumor suppressor genes. Expression of both genes is downregulated. However, differential promoter methylation suggests the existence of distinctive mechanisms of transcriptional inactivation for these genes.

DNA methylation changes in murine breast adenocarcinomas allow the identification of candidate genes for human breast carcinogenesis

Epigenetic inactivation due to aberrant promoter methylation is a key process in breast tumorigenesis. Murine models for human breast cancer have been established for nearly every important human oncogene or tumor suppressor gene. Mouse-to-human comparative gene expression and cytogenetic profiling have been widely investigated for these models; however, little is known about the conservation of epigenetic alterations during tumorigenesis. To determine if this key process in human breast tumorigenesis is also mirrored in a murine breast cancer model, we mapped cytosine methylation changes in primary adenocarcinomas and paired lung metastases derived from the polyomavirus middle T antigen mouse model. Global changes in methylcytosine levels were observed in all tumors when compared to the normal mammary gland. Aberrant methylation and associated gene silencing was observed for Hoxa7, a gene that is differentially methylated in human breast tumors, and Gata2, a novel candidate gene. Analysis of HOXA7 and GATA2 expression in a bank of human primary tumors confirms that the expression of these genes is also reduced in human breast cancer. In addition, HOXA7 hypermethylation is observed in breast cancer tissues when compared to adjacent tumor-free tissue. Based on these studies, we present a model in which comparative epigenetic techniques can be used to identify novel candidate genes important for human breast tumorigenesis, in both primary and metastatic tumors.