DNA methylation and Mycoplasma genomes

DNA methylation is one of the many hypotheses proposed to explain the observed deficiency in CpG dinucleotides in a variety of genomes covering a wide taxonomic distribution. Recent studies challenged the methylation hypothesis on empirical grounds. First, it cannot explain why the Mycoplasma genitalium genome exhibits strong CpG deficiency without DNA methylation. Second, it cannot explain the great variation in CpG deficiency between M. genitalium and M. pneumoniae that also does not have CpG-specific methyltransferase genes. I analyzed the genomic sequences of these Mycoplasma species together with the recently sequenced genomes of M. pulmonis, Ureaplasma urealyticum, and Staphylococcus aureus, and found the results fully compatible with the methylation hypothesis. In particular, I present compelling empirical evidence to support the following scenario. The common ancestor of the three Mycoplasma species has CpG-specific methyltransferases, and has evolved strong CpG deficiency as a result of the specific DNA methylation. Subsequently, this ancestral genome diverged into M. pulmonis and the common ancestor of M. pneumoniae and M. genitalium. M. pulmonis has retained methyltransferases and exhibits the strongest CpG deficiency. The common ancestor lost the methyltransferase gene and then diverged into M. genitalium and M. pneumoniae. M. genitalium and M. pneumoniae, after losing methylation activities, began to regain CpG dinucleotides through random mutation. M. genitalium evolved more slowly than M. pneumoniae, gained relatively fewer CpG dinucleotides, and is more CpG-deficient.     

Genome-wide analysis of DNA methylation patterns

Cytosine methylation is the most common covalent modification of DNA in eukaryotes. DNA methylation has an important role in many aspects of biology, including development and disease. Methylation can be detected using bisulfite conversion, methylation-sensitive restriction enzymes, methyl-binding proteins and anti-methylcytosine antibodies. Combining these techniques with DNA microarrays and high-throughput sequencing has made the mapping of DNA methylation feasible on a genome-wide scale. Here we discuss recent developments and future directions for identifying and mapping methylation, in an effort to help colleagues to identify the approaches that best serve their research interests.     

Genome-wide DNA methylation analysis in precursor B-cells

DNA methylation is responsible for regulating gene expression and cellular differentiation and for maintaining genomic stability during normal human development. Furthermore, it plays a significant role in the regulation of hematopoiesis. In order to elucidate the influence of DNA methylation during B-cell development, genome-wide DNA methylation status of pro-B, pre-BI, pre-BII, and naïve-B-cells isolated from human umbilical cord blood was determined using the methylated CpG island recovery assay followed by next generation sequencing. On average, 182–200 million sequences were generated for each precursor B-cell subset in 10 biological replicates. An overall decrease in methylation was observed during the transition from pro-B to pre-BI, whereas no differential methylation was observed in the pre-BI to pre-BII transition or in the pre-BII to naïve B-cell transition. Most of the methylated regions were located within intergenic and intronic regions not present in a CpG island context. Putative novel enhancers were identified in these regions that were differentially methylated between pro-B and pre-BI cells. The genome-wide methylation profiles are publically available and may be used to gain a better understanding of the involvement of atypical DNA methylation in the pathogenesis of malignancies associated with precursor B-cells.     

Genome-wide DNA methylation analysis in precursor B-cells

DNA methylation is responsible for regulating gene expression and cellular differentiation and for maintaining genomic stability during normal human development. Furthermore, it plays a significant role in the regulation of hematopoiesis. In order to elucidate the influence of DNA methylation during B-cell development, genome-wide DNA methylation status of pro-B, pre-BI, pre-BII, and naïve-B-cells isolated from human umbilical cord blood was determined using the methylated CpG island recovery assay followed by next generation sequencing. On average, 182–200 million sequences were generated for each precursor B-cell subset in 10 biological replicates. An overall decrease in methylation was observed during the transition from pro-B to pre-BI, whereas no differential methylation was observed in the pre-BI to pre-BII transition or in the pre-BII to naïve B-cell transition. Most of the methylated regions were located within intergenic and intronic regions not present in a CpG island context. Putative novel enhancers were identified in these regions that were differentially methylated between pro-B and pre-BI cells. The genome-wide methylation profiles are publically available and may be used to gain a better understanding of the involvement of atypical DNA methylation in the pathogenesis of malignancies associated with precursor B-cells.     

Genome-wide mapping of DNA methylation in Nile Tilapia

Cytosine DNA methylation is crucial for gene regulation and maintenance of genome stability. However, the detailed nile tilapia methylome remains uncharacterized. In this study, we present the first high-resolution methylome of tilapia gonad generated using methylated DNA immunoprecipitation (MeDIP) and high-throughput sequencing. In the ovary, 265 and 56 methylation peaks were identified in the genebody and promoter region of 145 genes, respectively. In the testis, 293 and 80 methylation peaks were identified in the genebody and promoter region of 144 genes. Furthermore, 8 and 49 genes showed differentially higher and lower promoter-region methylation rates, respectively, in the ovary relative to those of the testis. Quantitative PCR results revealed that the expression level of fibroblast growth factor 16 (fgf16), sialidase-3-like, fibroblast growth factor 20, aromatase (cyp19a), estrogen receptor, and gonadotropin receptor II precursor were negatively correlated to their methylation levels in the ovary and testis. The methylated levels of cyp19a and fgf16 were validated by bisulfite sequencing PCR technology, and the results were consistent with the MeDIP results. Thus, apart from generating the first methylation map, this study produced a candidate gene repository that provides additional options to explore the relationship between DNA methylation and sex differentiation or maintenance.     

Genome-wide analysis of retroviral DNA integration

Retroviral vectors are often used to introduce therapeutic sequences into patients' cells. In recent years, gene therapy with retroviral vectors has had impressive therapeutic successes, but has also resulted in three cases of leukaemia caused by insertional mutagenesis, which has focused attention on the molecular determinants of retroviral-integration target-site selection. Here, we review retroviral DNA integration, with emphasis on recent genome-wide studies of targeting and on the status of efforts to modulate target-site selection.     

Genome-wide DNA methylation studies suggest distinct DNA methylation patterns in pediatric embryonal and alveolar rhabdomyosarcomas

Rhabdomyosarcoma is the most common soft-tissue sarcoma in children. While cytogenetic abnormalities have been well characterized in this disease, aberrant epigenetic events such as DNA hypermethylation have not been described in genome-wide studies. We have analyzed the methylation status of 25,500 promoters in normal skeletal muscle, and in cell lines and tumor samples of embryonal and alveolar rhabdomyosarcoma from pediatric patients. We identified over 1,900 CpG islands that are hypermethylated in rhabdomyosarcomas relative to skeletal muscle. Genes involved in tissue development, differentiation, and oncogenesis such as DNAJA4, HES5, IRX1, BMP8A, GATA4, GATA6, ALX3, and P4HTM were hypermethylated in both RMS cell lines and primary samples, implicating aberrant DNA methylation in the pathogenesis of rhabdomyosarcoma. Furthermore, cluster analysis revealed embryonal and alveolar subtypes had distinct DNA methylation patterns, with the alveolar subtype being enriched in DNA hypermethylation of polycomb target genes. These results suggest that DNA methylation signatures may aid in the diagnosis and risk stratification of pediatric rhabdomyosarcoma and help identify new targets for therapy.     

Genome-wide DNA methylation studies suggest distinct DNA methylation patterns in pediatric embryonal and alveolar rhabdomyosarcomas

Rhabdomyosarcoma is the most common soft-tissue sarcoma in children. While cytogenetic abnormalities have been well characterized in this disease, aberrant epigenetic events such as DNA hypermethylation have not been described in genome-wide studies. We have analyzed the methylation status of 25,500 promoters in normal skeletal muscle, and in cell lines and tumor samples of embryonal and alveolar rhabdomyosarcoma from pediatric patients. We identified over 1,900 CpG islands that are hypermethylated in rhabdomyosarcomas relative to skeletal muscle. Genes involved in tissue development, differentiation, and oncogenesis such as DNAJA4, HES5, IRX1, BMP8A, GATA4, GATA6, ALX3, and P4HTM were hypermethylated in both RMS cell lines and primary samples, implicating aberrant DNA methylation in the pathogenesis of rhabdomyosarcoma. Furthermore, cluster analysis revealed embryonal and alveolar subtypes had distinct DNA methylation patterns, with the alveolar subtype being enriched in DNA hypermethylation of polycomb target genes. These results suggest that DNA methylation signatures may aid in the diagnosis and risk stratification of pediatric rhabdomyosarcoma and help identify new targets for therapy.     

Genome-wide DNA methylation scan in major depressive disorder

While genome-wide association studies are ongoing to identify sequence variation influencing susceptibility to major depressive disorder (MDD), epigenetic marks, such as DNA methylation, which can be influenced by environment, might also play a role. Here we present the first genome-wide DNA methylation (DNAm) scan in MDD. We compared 39 postmortem frontal cortex MDD samples to 26 controls. DNA was hybridized to our Comprehensive High-throughput Arrays for Relative Methylation (CHARM) platform, covering 3.5 million CpGs. CHARM identified 224 candidate regions with DNAm differences >10%. These regions are highly enriched for neuronal growth and development genes. Ten of 17 regions for which validation was attempted showed true DNAm differences; the greatest were in PRIMA1, with 12-15% increased DNAm in MDD (p = 0.0002-0.0003), and a concomitant decrease in gene expression. These results must be considered pilot data, however, as we could only test replication in a small number of additional brain samples (n = 16), which showed no significant difference in PRIMA1. Because PRIMA1 anchors acetylcholinesterase in neuronal membranes, decreased expression could result in decreased enzyme function and increased cholinergic transmission, consistent with a role in MDD. We observed decreased immunoreactivity for acetylcholinesterase in MDD brain with increased PRIMA1 DNAm, non-significant at p = 0.08.While we cannot draw firm conclusions about PRIMA1 DNAm in MDD, the involvement of neuronal development genes across the set showing differential methylation suggests a role for epigenetics in the illness. Further studies using limbic system brain regions might shed additional light on this role.     

Genome-wide high throughput analysis of DNA methylation in eukaryotes

Cytosine methylation is the quintessential epigenetic mark. Two well-established methods, bisulfite sequencing and methyl-DNA immunoprecipitation (MeDIP) lend themselves to the genome-wide analysis of DNA methylation by high throughput sequencing. Here we provide an overview and brief review of these methods. We summarize our experience with MeDIP followed by high throughput Illumina/Solexa sequencing, exemplified by the analysis of the methylated fraction of the Neurospora crassa genome ("methylome"). We provide detailed methods for DNA isolation, processing and the generation of in vitro libraries for Illumina/Solexa sequencing. We discuss potential problems in the generation of sequencing libraries. Finally, we provide an overview of software that is appropriate for the analysis of high throughput sequencing data generated by Illumina/Solexa-type sequencing by synthesis, with a special emphasis on approaches and applications that can generate more accurate depictions of sequence reads that fall in repeated regions of a chosen reference genome.     

Genome-wide analysis of differential DNA methylation in Silver-Russell syndrome

Silver-Russell Syndrome(SRS) is clinically heterogeneous disorder characterized by low birth weight, postnatal growth restriction, and variable dysmorphic features. Current evidence strongly implicates imprinted genes as an important etiology of SRS. Although almost half of the patients showed DNA hypomethylation at the H19/IGF2 imprinted domain, and approximately7%–10% of SRS patients have maternal uniparental disomy of chromosome 7(UPD(7) mat); the rest of the SRS patients shows unknown etiology. In this study, we investigate whether there are further DNA methylation defects in SRS patients. We measured DNA methylation in seven SRS patients and five controls at more than 485,000 CpG sites using DNA methylation microarrays. We analyzed methylation changes genome-wide and identified the differentially methylated regions(DMRs) using bisulfite sequencing and digital PCR. Our analysis identifies epimutations at the previously characterized domains of H19/IGF2,providing proof of principle that our methodology can detect the changes in DNA methylation at imprinted loci. In addition,our results showed a novel SRS associated imprinted gene OSBPL5 located on chromosome 11p14 with the probe cg25963939,which is hypomethylated in 4/7 patients(P=0.023, β=.0.243). We also report DMRs in other genes including TGFβ3, HSF1,GAP43, NOTCH4 and MYH14. These DMRs were found to be associated with SRS using GO pathway analysis. In this study,we identified the probe cg25963939, located at the 5′UTR of imprinted gene OSBPL5, as a novel DMR that is associated with SRS. This finding provides new insights into the mechanism of SRS etiology and aid the further stratification of SRS patients by molecular phenotypes.     

Genome-wide DNA methylation profiling of recurrent and non-recurrent chordomas

Chordomas are an aggressive rare type of malignant bone tumors arising from the remnant of the notochord. Chordomas occur mainly in vertebral bones and account for 1–4% of malignant bone tumors. Management and treatment of chordomas are difficult as they are resistant to conventional chemotherapy; therefore, they are mainly treated with surgery and radiation therapy. In this study, we performed DNA methylation profiling of 26 chordomas and normal nucleus pulposus samples plus UCH-1 chordoma cell line using the Illumina Infinium HumanMethylation450 BeadChips. Combined bisulfite restriction analysis and bisulfite sequencing was used to confirm the methylation data. Gene expression was analyzed using RT-PCR before and after 5-aza-2’-deoxycytidine (5-azaDC) treatment of chordoma cell lines. Analysis of the HumanMethylation450 BeadChip data led to the identification of 8,819 loci (2.9%) that were significantly differentially methylated (>0.2 average β-value difference) between chordomas and nucleus pulposus samples (adjusted P < 0.05). Among these, 5,868 probes (66.5%) were hypomethylated, compared to 2,951 (33.5%) loci that were hypermethylated in chordomas compared to controls. From the 2,951 differentially hypermethylated probes, 33.3% were localized in the promoter region (982 probes) and, among these, 104 probes showed cancer-specific hypermethylation. Ingenuity Pathway Analysis indicates that the cancer-specific differentially methylated loci are involved in various networks including cancer disease, nervous system development and function, cell death and survival, cellular growth, cellular development, and proliferation. Furthermore, we identified a subset of probes that were differentially methylated between recurrent and non-recurrent chordomas. BeadChip methylation data was confirmed for these genes and gene expression was shown to be upregulated in methylated chordoma cell lines after treatment with 5-azaDC. Understanding epigenetic changes in chordomas may provide insights into chordoma tumorigenesis and development of epigenetic biomarkers.     

Selection againstdam methylation sites in the genomes of DNA of enterobacteriophages

Summary Postreplicative methylation of adenine inEscherichia coli DNA to produce G^6m ATC (where^6mA is 6-methyladenine) has been associated with preferential daughter-strand repair and possibly regulation of replication. An analysis was undertaken to determine if these, or other, as yet unknown roles of GATC, have had an effect on the frequency of GATC inE. coli or bacteriophage DNA. It was first ascertained that the most accurate predictions of GATC frequency were based on the observed frequencies of GAT and ATC, which would be expected since these predictors take into account preferences in codon usage. The predicted frequencies were compared with observed GATC frequencies in all available bacterial and phage nucleotide sequences. The frequency of GATC was close to the predicted frequency in most genes ofE. coli and its RNA bacteriophages and in the genes of nonenteric bacteria and their bacteriophages. However, for DNA enterobacteriophages the observed frequency of DNA enterobacteriophages the observed frequency of GATC was generally significantly lower than predicted when assessed by the chi square test. No elevation in the rate of mutation of^6mA in GATC relative to other bases was found when pairs of DNA sequences from closely related phages or pairs of homologous genes from enterobacteria were compared, nor was any preferred pathway for mutation of^6mA evident in theE. coli DNA bacteriophages. This situation contrasts with that of 5-methylcytosine, which is hypermutable, with a preferred pathway to thymine. Thus, the low level of GATC in enterobacteriophages is probably due not to^6mA hypermutability, but to selection against GATC in order to bypass a GATC-mediated host function.     

An overview of the analysis of DNA methylation in mammalian genomes

DNA methylation at position C5 of the pyrimidine ring of cytosine in mammalian genomes has received a great deal of research interest due to its importance in many biological phenomena. It is associated with events such as epigenetic gene silencing and the maintenance of genome integrity. Aberrant DNA methylation, particularly that of chromosomal regions called CpG islands, is an important step in carcinogenesis. In order to elucidate methylation profiling of complex genomes, various methods have been developed. Many of these methods are based on the differential reactivity of cytosine and 5-methylcytosine to various chemicals. The combined use of these chemical reactions and other preexisting methods has enabled the discrimination of cytosine and 5-methylcytosine in complex genomes. The use of proteins that preferentially bind to methylated DNA has also successfully been used to discriminate between methylated and unmethylated sites. The chemical and structural dissection of the in vivo processes of enzymatic methylation and the binding of methyl-CpG binding proteins provides evidence for the complex mechanisms that nature has acquired. In this review we summarize the methods available for the discrimination between cytosine and 5-methylcytosine in complex genomes.