Comparative DNA methylome analysis of endometrial carcinoma reveals complex and distinct deregulation of cancer promoters and enhancers

Background

Aberrant DNA methylation is a hallmark of many cancers. Classically there are two types of endometrial cancer, endometrioid adenocarcinoma (EAC), or Type I, and uterine papillary serous carcinoma (UPSC), or Type II. However, the whole genome DNA methylation changes in these two classical types of endometrial cancer is still unknown.

Results

Here we described complete genome-wide DNA methylome maps of EAC, UPSC, and normal endometrium by applying a combined strategy of methylated DNA immunoprecipitation sequencing (MeDIP-seq) and methylation-sensitive restriction enzyme digestion sequencing (MRE-seq). We discovered distinct genome-wide DNA methylation patterns in EAC and UPSC: 27,009 and 15,676 recurrent differentially methylated regions (DMRs) were identified respectively, compared with normal endometrium. Over 80% of DMRs were in intergenic and intronic regions. The majority of these DMRs were not interrogated on the commonly used Infinium 450K array platform. Large-scale demethylation of chromosome X was detected in UPSC, accompanied by decreased XIST expression. Importantly, we discovered that the majority of the DMRs harbored promoter or enhancer functions and are specifically associated with genes related to uterine development and disease. Among these, abnormal methylation of transposable elements (TEs) may provide a novel mechanism to deregulate normal endometrium-specific enhancers derived from specific TEs.

Conclusions

DNA methylation changes are an important signature of endometrial cancer and regulate gene expression by affecting not only proximal promoters but also distal enhancers.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-868) contains supplementary material, which is available to authorized users.     

Genome-wide DNA methylation changes in skeletal muscle between young and middle-aged pigs

Background

Age-related physiological, biochemical and functional changes in mammalian skeletal muscle have been shown to begin at the mid-point of the lifespan. However, the underlying changes in DNA methylation that occur during this turning point of the muscle aging process have not been clarified. To explore age-related genomic methylation changes in skeletal muscle, we employed young (0.5 years old) and middle-aged (7 years old) pigs as models to survey genome-wide DNA methylation in the longissimus dorsi muscle using a methylated DNA immunoprecipitation sequencing approach.

Results

We observed a tendency toward a global loss of DNA methylation in the gene-body region of the skeletal muscle of the middle-aged pigs compared with the young group. We determined the genome-wide gene expression pattern in the longissimus dorsi muscle using microarray analysis and performed a correlation analysis using DMR (differentially methylated region)-mRNA pairs, and we found a significant negative correlation between the changes in methylation levels within gene bodies and gene expression. Furthermore, we identified numerous genes that show age-related methylation changes that are potentially involved in the aging process. The methylation status of these genes was confirmed using bisulfite sequencing PCR. The genes that exhibited a hypomethylated gene body in middle-aged pigs were over-represented in various proteolysis and protein catabolic processes, suggesting an important role for these genes in age-related muscle atrophy. In addition, genes associated with tumorigenesis exhibited aged-related differences in methylation and expression levels, suggesting an increased risk of disease associated with increased age.

Conclusions

This study provides a comprehensive analysis of genome-wide DNA methylation patterns in aging pig skeletal muscle. Our findings will serve as a valuable resource in aging studies, promoting the pig as a model organism for human aging research and accelerating the development of comparative animal models in aging research.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-653) contains supplementary material, which is available to authorized users.     

Genome-wide and single-base resolution DNA methylomes of the Pacific oyster Crassostrea gigas provide insight into the evolution of invertebrate CpG methylation

Background

Studies of DNA methylomes in a wide range of eukaryotes have revealed both conserved and divergent characteristics of DNA methylation among phylogenetic groups. However, data on invertebrates particularly molluscs are limited, which hinders our understanding of the evolution of DNA methylation in metazoa. The sequencing of the Pacific oyster Crassostrea gigas genome provides an opportunity for genome-wide profiling of DNA methylation in this model mollusc.

Results

Homologous searches against the C. gigas genome identified functional orthologs for key genes involved in DNA methylation: DNMT1, DNMT2, DNMT3, MBD2/3 and UHRF1. Whole-genome bisulfite sequencing (BS-seq) of the oyster’s mantle tissues revealed that more than 99% methylation modification was restricted to cytosines in CpG context and methylated CpGs accumulated in the bodies of genes that were moderately expressed. Young repeat elements were another major targets of CpG methylation in oysters. Comparison with other invertebrate methylomes suggested that the 5’-end bias of gene body methylation and the negative correlation between gene body methylation and gene length were the derived features probably limited to the insect lineage. Interestingly, phylostratigraphic analysis showed that CpG methylation preferentially targeted genes originating in the common ancestor of eukaryotes rather than the oldest genes originating in the common ancestor of cellular organisms.

Conclusions

Comparative analysis of the oyster DNA methylomes and that of other animal species revealed that the characteristics of DNA methylation were generally conserved during invertebrate evolution, while some unique features were derived in the insect lineage. The preference of methylation modification on genes originating in the eukaryotic ancestor rather than the oldest genes is unexpected, probably implying that the emergence of methylation regulation in these ''relatively young’ genes was critical for the origin and radiation of eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1119) contains supplementary material, which is available to authorized users.     

MAGEB2 is Activated by Promoter Demethylation in Head and Neck Squamous Cell Carcinoma

Purpose

Although promoter hypermethylation has been an accepted means of tumor suppressor gene inactivation, activation of otherwise normally repressed proto-oncogenes by promoter demethylation has been infrequently documented.

Experimental Design

In this study we performed an integrative, whole-genome analysis for discovery of epigenetically activated proto-oncogenes in head and neck cancer tumors. We used the 47K GeneChip U133 Plus 2.0 Affymetrix expression microarray platform to obtain re-expression data from 5-aza treated normal cell line and expression data from primary head and neck squamous cell carcinoma (HNSCC) tumor tissues and normal mucosa tissues. We then investigated candidate genes by screening promoter regions for CpG islands and bisulfite sequencing followed by QUMSP and RT PCR for the best candidate genes. Finally, functional studies were performed on the top candidate gene.

Results

From the top 178 screened candidates 96 had CpG islands in their promoter region. Seven candidate genes showed promoter region methylation in normal mucosa samples and promoter demethylation in a small cohort of primary HNSCC tissues. We then studied the demethylation of the top 3 candidate genes in an expanded cohort of 76 HNSCC tissue samples and 17 normal mucosa samples. We identified MAGEB2 as having significant promoter demethylation in primary head and neck squamous cell carcinoma tissues. We then found significantly higher expression of MAGEB2 in tumors in a separate cohort of 73 primary HNSCC tissues and 31 normal tissues. Finally, we found that MAGEB2 has growth promoting effects on minimally transformed oral keratinocyte cell lines but not a definite effect on HNSCC cell lines.

Conclusion

In conclusion, we identified MAGEB2 as activated by promoter demethylation in HNSCCand demonstrates growth promoting effects in a minimally transformed oral keratinocyte cell line. More studies are needed to evaluate MAGBE2''s exact role in HNSCC.     

OPCML is a broad tumor suppressor for multiple carcinomas and lymphomas with frequently epigenetic inactivation

Background

Identification of tumor suppressor genes (TSGs) silenced by CpG methylation uncovers the molecular mechanism of tumorigenesis and potential tumor biomarkers. Loss of heterozygosity at 11q25 is common in multiple tumors including nasopharyngeal carcinoma (NPC). OPCML, located at 11q25, is one of the downregulated genes we identified through digital expression subtraction.

Methodology/Principal Findings

Semi-quantitative RT-PCR showed frequent OPCML silencing in NPC and other common tumors, with no homozygous deletion detected by multiplex differential DNA-PCR. Instead, promoter methylation of OPCML was frequently detected in multiple carcinoma cell lines (nasopharyngeal, esophageal, lung, gastric, colon, liver, breast, cervix, prostate), lymphoma cell lines (non-Hodgkin and Hodgkin lymphoma, nasal NK/T-cell lymphoma) and primary tumors, but not in any non-tumor cell line and seldom weakly methylated in normal epithelial tissues. Pharmacological and genetic demethylation restored OPCML expression, indicating a direct epigenetic silencing. We further found that OPCML is stress-responsive, but this response is epigenetically impaired when its promoter becomes methylated. Ecotopic expression of OPCML led to significant inhibition of both anchorage-dependent and -independent growth of carcinoma cells with endogenous silencing.

Conclusions/Significance

Thus, through functional epigenetics, we identified OPCML as a broad tumor suppressor, which is frequently inactivated by methylation in multiple malignancies.     

Genes suppressed by DNA methylation in non-small cell lung cancer reveal the epigenetics of epithelial–mesenchymal transition

Background

DNA methylation is associated with aberrant gene expression in cancer, and has been shown to correlate with therapeutic response and disease prognosis in some types of cancer. We sought to investigate the biological significance of DNA methylation in lung cancer.

Results

We integrated the gene expression profiles and data of gene promoter methylation for a large panel of non-small cell lung cancer cell lines, and identified 578 candidate genes with expression levels that were inversely correlated to the degree of DNA methylation. We found these candidate genes to be differentially methylated in normal lung tissue versus non-small cell lung cancer tumors, and segregated by histologic and tumor subtypes. We used gene set enrichment analysis of the genes ranked by the degree of correlation between gene expression and DNA methylation to identify gene sets involved in cellular migration and metastasis. Our unsupervised hierarchical clustering of the candidate genes segregated cell lines according to the epithelial-to-mesenchymal transition phenotype. Genes related to the epithelial-to-mesenchymal transition, such as AXL, ESRP1, HoxB4, and SPINT1/2, were among the nearly 20% of the candidate genes that were differentially methylated between epithelial and mesenchymal cells. Greater numbers of genes were methylated in the mesenchymal cells and their expressions were upregulated by 5-azacytidine treatment. Methylation of the candidate genes was associated with erlotinib resistance in wild-type EGFR cell lines. The expression profiles of the candidate genes were associated with 8-week disease control in patients with wild-type EGFR who had unresectable non-small cell lung cancer treated with erlotinib, but not in patients treated with sorafenib.

Conclusions

Our results demonstrate that the underlying biology of genes regulated by DNA methylation may have predictive value in lung cancer that can be exploited therapeutically.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1079) contains supplementary material, which is available to authorized users.     

GAP-Seq: a method for identification of DNA palindromes

Background

Closely spaced long inverted repeats, also known as DNA palindromes, can undergo intrastrand annealing to form DNA hairpins. The ability to form these hairpins results in genome instability, difficulties in maintaining clones in Escherichia coli and major problems for most DNA sequencing approaches. Because of their role in genomic instability and gene amplification in some human cancers, it is important to develop systematic approaches to detect and characterize DNA palindromes.

Results

We developed a new protocol to identify palindromes that couples the S1 nuclease treated Cot0 DNA (GAPF) with high-throughput sequencing (GAP-Seq). Unlike earlier protocols, it does not involve restriction enzymatic digestion prior to DNA snap-back thereby preserving longer DNA sequences. It also indicates the location of the novel junction, which can then be recovered. Using MCF-7 breast cancer cell line as the proof-of-principle analysis, we have identified 35 palindrome candidates and physically characterized the top 5 candidates and their junctions. Because this protocol eliminates many of the false positives that plague earlier techniques, we have improved palindrome identification.

Conclusions

The GAP-Seq approach underscores the importance of developing new tools for identifying and characterizing palindromes, and provides a new strategy to systematically assess palindromes in genomes. It will be useful for studying human cancers and other diseases associated with palindromes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-394) contains supplementary material, which is available to authorized users.     

Characterization of the MLO gene family in Rosaceae and gene expression analysis in Malus domestica

Background

Powdery mildew (PM) is a major fungal disease of thousands of plant species, including many cultivated Rosaceae. PM pathogenesis is associated with up-regulation of MLO genes during early stages of infection, causing down-regulation of plant defense pathways. Specific members of the MLO gene family act as PM-susceptibility genes, as their loss-of-function mutations grant durable and broad-spectrum resistance.

Results

We carried out a genome-wide characterization of the MLO gene family in apple, peach and strawberry, and we isolated apricot MLO homologs through a PCR-approach. Evolutionary relationships between MLO homologs were studied and syntenic blocks constructed. Homologs that are candidates for being PM susceptibility genes were inferred by phylogenetic relationships with functionally characterized MLO genes and, in apple, by monitoring their expression following inoculation with the PM causal pathogen Podosphaera leucotricha.

Conclusions

Genomic tools available for Rosaceae were exploited in order to characterize the MLO gene family. Candidate MLO susceptibility genes were identified. In follow-up studies it can be investigated whether silencing or a loss-of-function mutations in one or more of these candidate genes leads to PM resistance.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-618) contains supplementary material, which is available to authorized users.     

The UBC-40 Urothelial Bladder Cancer cell line index: a genomic resource for functional studies

Background

Urothelial bladder cancer is a highly heterogeneous disease. Cancer cell lines are useful tools for its study. This is a comprehensive genomic characterization of 40 urothelial bladder carcinoma (UBC) cell lines including information on origin, mutation status of genes implicated in bladder cancer (FGFR3, PIK3CA, TP53, and RAS), copy number alterations assessed using high density SNP arrays, uniparental disomy (UPD) events, and gene expression.

Results

Based on gene mutation patterns and genomic changes we identify lines representative of the FGFR3-driven tumor pathway and of the TP53/RB tumor suppressor-driven pathway. High-density array copy number analysis identified significant focal gains (1q32, 5p13.1-12, 7q11, and 7q33) and losses (i.e. 6p22.1) in regions altered in tumors but not previously described as affected in bladder cell lines. We also identify new evidence for frequent regions of UPD, often coinciding with regions reported to be lost in tumors. Previously undescribed chromosome X losses found in UBC lines also point to potential tumor suppressor genes. Cell lines representative of the FGFR3-driven pathway showed a lower number of UPD events.

Conclusions

Overall, there is a predominance of more aggressive tumor subtypes among the cell lines. We provide a cell line classification that establishes their relatedness to the major molecularly-defined bladder tumor subtypes. The compiled information should serve as a useful reference to the bladder cancer research community and should help to select cell lines appropriate for the functional analysis of bladder cancer genes, for example those being identified through massive parallel sequencing.

Electronic supplementary material

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

The complete methylome of Helicobacter pylori UM032

Background

The genome of the human gastric pathogen Helicobacter pylori encodes a large number of DNA methyltransferases (MTases), some of which are shared among many strains, and others of which are unique to a given strain. The MTases have potential roles in the survival of the bacterium. In this study, we sequenced a Malaysian H. pylori clinical strain, designated UM032, by using a combination of PacBio Single Molecule, Real-Time (SMRT) and Illumina MiSeq next generation sequencing platforms, and used the SMRT data to characterize the set of methylated bases (the methylome).

Results

The N4-methylcytosine and N6-methyladenine modifications detected at single-base resolution using SMRT technology revealed 17 methylated sequence motifs corresponding to one Type I and 16 Type II restriction-modification (R-M) systems. Previously unassigned methylation motifs were now assigned to their respective MTases-coding genes. Furthermore, one gene that appears to be inactive in the H. pylori UM032 genome during normal growth was characterized by cloning.

Conclusion

Consistent with previously-studied H. pylori strains, we show that strain UM032 contains a relatively large number of R-M systems, including some MTase activities with novel specificities. Additional studies are underway to further elucidating the biological significance of the R-M systems in the physiology and pathogenesis of H. pylori.

Electronic supplementary material

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

Secretome analysis reveals effector candidates associated with broad host range necrotrophy in the fungal plant pathogen Sclerotinia sclerotiorum

Background

The white mold fungus Sclerotinia sclerotiorum is a devastating necrotrophic plant pathogen with a remarkably broad host range. The interaction of necrotrophs with their hosts is more complex than initially thought, and still poorly understood.

Results

We combined bioinformatics approaches to determine the repertoire of S. sclerotiorum effector candidates and conducted detailed sequence and expression analyses on selected candidates. We identified 486 S. sclerotiorum secreted protein genes expressed in planta, many of which have no predicted enzymatic activity and may be involved in the interaction between the fungus and its hosts. We focused on those showing (i) protein domains and motifs found in known fungal effectors, (ii) signatures of positive selection, (iii) recent gene duplication, or (iv) being S. sclerotiorum-specific. We identified 78 effector candidates based on these properties. We analyzed the expression pattern of 16 representative effector candidate genes on four host plants and revealed diverse expression patterns.

Conclusions

These results reveal diverse predicted functions and expression patterns in the repertoire of S. sclerotiorum effector candidates. They will facilitate the functional analysis of fungal pathogenicity determinants and should prove useful in the search for plant quantitative disease resistance components active against the white mold.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-336) contains supplementary material, which is available to authorized users.     

Ontogeny,distribution and potential roles of 5-hydroxymethylcytosine in human liver function

Background

Interindividual differences in liver functions such as protein synthesis, lipid and carbohydrate metabolism and drug metabolism are influenced by epigenetic factors. The role of the epigenetic machinery in such processes has, however, been barely investigated. 5-hydroxymethylcytosine (5hmC) is a recently re-discovered epigenetic DNA modification that plays an important role in the control of gene expression.

Results

In this study, we investigate 5hmC occurrence and genomic distribution in 8 fetal and 7 adult human liver samples in relation to ontogeny and function. LC-MS analysis shows that in the adult liver samples 5hmC comprises up to 1% of the total cytosine content, whereas in all fetal livers it is below 0.125%. Immunohistostaining of liver sections with a polyclonal anti-5hmC antibody shows that 5hmC is detected in most of the hepatocytes. Genome-wide mapping of the distribution of 5hmC in human liver samples by next-generation sequencing shows significant differences between fetal and adult livers. In adult livers, 5hmC occupancy is overrepresented in genes involved in active catabolic and metabolic processes, whereas 5hmC elements which are found in genes exclusively in fetal livers and disappear in the adult state, are more specific to pathways for differentiation and development.

Conclusions

Our findings suggest that 5-hydroxymethylcytosine plays an important role in the development and function of the human liver and might be an important determinant for development of liver diseases as well as of the interindividual differences in drug metabolism and toxicity.     

Bimodal signatures of germline methylation are linked with gene expression plasticity in the coral Acropora millepora

Background

In invertebrates, genes belonging to dynamically regulated functional categories appear to be less methylated than “housekeeping” genes, suggesting that DNA methylation may modulate gene expression plasticity. To date, however, experimental evidence to support this hypothesis across different natural habitats has been lacking.

Results

Gene expression profiles were generated from 30 pairs of genetically identical fragments of coral Acropora millepora reciprocally transplanted between distinct natural habitats for 3 months. Gene expression was analyzed in the context of normalized CpG content, a well-established signature of historical germline DNA methylation. Genes with weak methylation signatures were more likely to demonstrate differential expression based on both transplant environment and population of origin than genes with strong methylation signatures. Moreover, the magnitude of expression differences due to environment and population were greater for genes with weak methylation signatures.

Conclusions

Our results support a connection between differential germline methylation and gene expression flexibility across environments and populations. Studies of phylogenetically basal invertebrates such as corals will further elucidate the fundamental functional aspects of gene body methylation in Metazoa.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1109) contains supplementary material, which is available to authorized users.