Genome-wide analysis reveals DNA methylation markers that vary with both age and obesity

The combination of the obesity epidemic and an aging population presents growing challenges for the healthcare system. Obesity and aging are major risk factors for a diverse number of diseases and it is of importance to understand their interaction and the underlying molecular mechanisms. Herein the authors examined the methylation levels of 27578 CpG sites in 46 samples from adult peripheral blood. The effect of obesity and aging was ascertained with general linear models. More than one hundred probes were correlated to aging, nine of which belonged to the KEGG group map04080. Additionally, 10 CpG sites had diverse methylation profiles in obese and lean individuals, one of which was the telomerase catalytic subunit (TERT). In eight of ten cases the methylation change was reverted between obese and lean individuals. One region proved to be differentially methylated with obesity (LINC00304) independent of age. This study provides evidence that obesity influences age driven epigenetic changes, which provides a molecular link between aging and obesity. This link and the identified markers may prove to be valuable biomarkers for the understanding of the molecular basis of aging, obesity and associated diseases.     

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.     

DNA methylation errors in cloned mice disappear with advancement of aging

Cloned animals have various health problems. Aberrant DNA methylation is a possible cause of the problems. Restriction landmark genomic scanning (RLGS) that enabled us to analyze more than 1,000 CpG islands simultaneously demonstrated that all cloned newborns had aberrant DNA methylation. To study whether this aberration persists throughout the life of cloned individuals, we examined genome-wide DNA methylation status of newborn (19.5 dpc, n=2), adult (8-11 months old, n=3), and aged (23-27 months old, n=4) cloned mice using kidney cells as representatives. In the adult and aged groups, cloning was repeated using cumulus cells of the adult founder clone of each group as nucleus donor. Two newborn clones had three with aberrantly methylated loci, which is consistent with previous reports that all cloned newborns had DNA methylation aberrations. Interestingly, we could detect only one aberrantly methylated locus in two of the three adult clones in mid-age and none of four senescent clones, indicating that errors in DNA methylation disappear with advancement of animals' aging.     

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 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 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.     

Proteomic characterization in the hippocampus of prenatally stressed rats

Rats exposed to early life stress are considered as a valuable model for the study of epigenetic programming leading to mood disorders and anxiety in the adult life. Rats submitted to prenatal restraint stress (PRS) are characterized by an anxious/depressive phenotype associated with neuroadaptive changes in the hippocampus. We used the model of PRS to identify proteins that are specifically affected by early life stress. We therefore performed a proteomic analysis in the hippocampus of adult male PRS rats. We found that PRS induced changes in the expression profile of a number of proteins, involved in the regulation of signal transduction, synaptic vesicles, protein synthesis, cytoskeleton dynamics, and energetic metabolism. Immunoblot analysis showed significant changes in the expression of proteins, such as LASP-1, fascin, and prohibitin, which may lie at the core of the developmental programming triggered by early life stress.     

Neonatal DNA methylation patterns associate with gestational age

Risk for adverse neonatal outcome increases with declining gestational age (GA), and changes in DNA methylation may contribute to the relationship between GA and adverse health outcomes in offspring. To test this hypothesis, we evaluated the association between GA and more than 27,000 CpG sites in neonatal DNA extracted from umbilical cord blood from two prospectively-characterized cohorts: (1) a discovery cohort consisting of 259 neonates from women with a history of neuropsychiatric disorders and (2) a replication cohort consisting of 194 neonates of uncomplicated mothers. GA was determined by obstetrician report and maternal last menstrual period. The associations between proportion of DNA methylated and GA were evaluated by fitting a separate linear mixed effects model for each CpG site, adjusting for relevant covariates including neonatal sex, race, parity, birth weight percentile and chip effects. CpG sites in 39 genes were associated with GA (false discovery rate < 0.05) in the discovery cohort. The same CpG sites in 25 of these genes replicated in the replication cohort, with each association replicating in the same direction. Notably, these CpG sites were located in genes previously implicated in labor and delivery (e.g., AVP, OXT, CRHBP and ESR1) or that may influence the risk for adverse health outcomes later in life (e.g., DUOX2, TMEM176A and CASP8). All associations were independent of method of delivery or induction of labor. These results suggest neonatal DNA methylation varies with GA even in term deliveries. The potential contribution of these changes to clinically significant postnatal outcomes warrants further investigation.     

Neonatal DNA methylation patterns associate with gestational age

Risk for adverse neonatal outcome increases with declining gestational age (GA), and changes in DNA methylation may contribute to the relationship between GA and adverse health outcomes in offspring. To test this hypothesis, we evaluated the association between GA and more than 27,000 CpG sites in neonatal DNA extracted from umbilical cord blood from two prospectively-characterized cohorts: (1) a discovery cohort consisting of 259 neonates from women with a history of neuropsychiatric disorders and (2) a replication cohort consisting of 194 neonates of uncomplicated mothers. GA was determined by obstetrician report and maternal last menstrual period. The associations between proportion of DNA methylated and GA were evaluated by fitting a separate linear mixed effects model for each CpG site, adjusting for relevant covariates including neonatal sex, race, parity, birth weight percentile and chip effects. CpG sites in 39 genes were associated with GA (false discovery rate &lt; 0.05) in the discovery cohort. The same CpG sites in 25 of these genes replicated in the replication cohort, with each association replicating in the same direction. Notably, these CpG sites were located in genes previously implicated in labor and delivery (e.g., AVP, OXT, CRHBP and ESR1) or that may influence the risk for adverse health outcomes later in life (e.g., DUOX2, TMEM176A and CASP8). All associations were independent of method of delivery or induction of labor. These results suggest neonatal DNA methylation varies with GA even in term deliveries. The potential contribution of these changes to clinically significant postnatal outcomes warrants further investigation.