A genome-wide study of DNA methylation patterns and gene expression levels in multiple human and chimpanzee tissues

Tissue-specific DNA methylation is conserved across human, mouse, and rat, and driven by primary sequence conservation

Background

Uncovering mechanisms of epigenome evolution is an essential step towards understanding the evolution of different cellular phenotypes. While studies have confirmed DNA methylation as a conserved epigenetic mechanism in mammalian development, little is known about the conservation of tissue-specific genome-wide DNA methylation patterns.

Results

Using a comparative epigenomics approach, we identified and compared the tissue-specific DNA methylation patterns of rat against those of mouse and human across three shared tissue types. We confirmed that tissue-specific differentially methylated regions are strongly associated with tissue-specific regulatory elements. Comparisons between species revealed that at a minimum 11-37% of tissue-specific DNA methylation patterns are conserved, a phenomenon that we define as epigenetic conservation. Conserved DNA methylation is accompanied by conservation of other epigenetic marks including histone modifications. Although a significant amount of locus-specific methylation is epigenetically conserved, the majority of tissue-specific DNA methylation is not conserved across the species and tissue types that we investigated. Examination of the genetic underpinning of epigenetic conservation suggests that primary sequence conservation is a driving force behind epigenetic conservation. In contrast, evolutionary dynamics of tissue-specific DNA methylation are best explained by the maintenance or turnover of binding sites for important transcription factors.

Conclusions

Our study extends the limited literature of comparative epigenomics and suggests a new paradigm for epigenetic conservation without genetic conservation through analysis of transcription factor binding sites.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4115-6) contains supplementary material, which is available to authorized users.

Epigenome-Wide Association Study Identifies Cardiac Gene Patterning and a Novel Class of Biomarkers for Heart Failure

BACKGROUND

Biochemical DNA modification resembles a crucial regulatory layer among genetic information, environmental factors, and the transcriptome. To identify epigenetic susceptibility regions and novel biomarkers linked to myocardial dysfunction and heart failure, we performed the first multi-omics study in myocardial tissue and blood of patients with dilated cardiomyopathy and controls.

METHODS

Infinium human methylation 450 was used for high-density epigenome-wide mapping of DNA methylation in left-ventricular biopsies and whole peripheral blood of living probands. RNA deep sequencing was performed on the same samples in parallel. Whole-genome sequencing of all patients allowed exclusion of promiscuous genotype-induced methylation calls.

RESULTS

In the screening stage, we detected 59 epigenetic loci that are significantly associated with dilated cardiomyopathy (false discovery corrected P≤0.05), with 3 of them reaching epigenome-wide significance at P≤5×10^-8. Twenty-seven (46%) of these loci could be replicated in independent cohorts, underlining the role of epigenetic regulation of key cardiac transcription regulators. Using a staged multi-omics study design, we link a subset of 517 epigenetic loci with dilated cardiomyopathy and cardiac gene expression. Furthermore, we identified distinct epigenetic methylation patterns that are conserved across tissues, rendering these CpGs novel epigenetic biomarkers for heart failure.

CONCLUSIONS

The present study provides to our knowledge the first epigenome-wide association study in living patients with heart failure using a multi-omics approach.

Measuring Animal Age with DNA Methylation: From Humans to Wild Animals

DNA methylation (DNAm) is a key mechanism for regulating gene expression in animals and levels are known to change with age. Recent studies have used DNAm changes as a biomarker to estimate chronological age in humans and these techniques are now also being applied to domestic and wild animals. Animal age is widely used to track ongoing changes in ecosystems, however chronological age information is often unavailable for wild animals. An ability to estimate age would lead to improved monitoring of (i) population trends and status and (ii) demographic properties such as age structure and reproductive performance. Recent studies have revealed new examples of DNAm age association in several new species increasing the potential for developing DNAm age biomarkers for a broad range of wild animals. Emerging technologies for measuring DNAm will also enhance our ability to study age-related DNAm changes and to develop new molecular age biomarkers.

Evolution of the Human Nervous System Function, Structure, and Development

The nervous system-in particular, the brain and its cognitive abilities-is among humans' most distinctive and impressive attributes. How the nervous system has changed in the human lineage and how it differs from that of closely related primates is not well understood. Here, we consider recent comparative analyses of extant species that are uncovering new evidence for evolutionary changes in the size and the number of neurons in the human nervous system, as well as the cellular and molecular reorganization of its neural circuits. We also discuss the developmental mechanisms and underlying genetic and molecular changes that generate these structural and functional differences. As relevant new information and tools materialize at an unprecedented pace, the field is now ripe for systematic and functionally relevant studies of the development and evolution of human nervous system specializations.

A complex association between DNA methylation and gene expression in human placenta at first and third trimesters

The human placenta is a maternal-fetal organ essential for normal fetal development and maternal health. During pregnancy, the placenta undergoes many structural and functional changes in response to fetal needs and environmental exposures. Previous studies have demonstrated widespread epigenetic and gene expression changes from early to late pregnancy. However, on the global level, how DNA methylation changes impact on gene expression in human placenta is not yet well understood. We performed DNA methylome analysis by reduced representation bisulfite sequencing (RRBS) and gene expression analysis by RNA-Seq for both first and third trimester human placenta tissues. From first to third trimester, 199 promoters (corresponding to 189 genes) and 2,297 gene bodies were differentially methylated, with a clear dominance of hypermethylation (96.8% and 93.0% for promoters and gene bodies, respectively). A total of 2,447 genes were differentially expressed, of which 77.2% were down-regulated. Gene ontology analysis using differentially expressed genes were enriched for cell cycle and immune response functions. The correlation between DNA methylation and gene expression was non-linear and complex, depending on the genomic context (promoter or gene body) and gene expression levels. A wide range of DNA methylation and gene expression changes were observed at different gestational ages. The non-linear association between DNA methylation and gene expression indicates that epigenetic regulation of placenta development is more complex than previously envisioned.

Regulation of DNA methylation on EEF1D and RPL8 expression in cattle

Dynamic changes to the epigenome play a critical role in a variety of biology processes and complex traits. Many important candidate genes have been identified through our previous genome wide association study (GWAS) on milk production traits in dairy cattle. However, the underlying mechanism of candidate genes have not yet been clearly understood. In this study, we analyzed the methylation variation of the candidate genes, EEF1D and RPL8, which were identified to be strongly associated with milk production traits in dairy cattle in our previous studies, and its effect on protein and mRNA expression. We compared DNA methylation profiles and gene expression levels of EEF1D and RPL8 in five different tissues (heart, liver, mammary gland, ovary and muscle) of three cows. Both genes showed the highest expression level in mammary gland. For RPL8, there was no difference in the DNA methylation pattern in the five tissues, suggesting no effect of DNA methylation on gene expression. For EEF1D, the DNA methylation levels of its first CpG island differed in the five tissues and were negatively correlated with the gene expression levels. To further investigate the function of DNA methylation on the expression of EEF1D, we collected blood samples of three cows at early stage of lactation and in dry period and analyzed its expression and the methylation status of the first CpG island in blood. As a result, the mRNA expression of EEF1D in the dry period was higher than that at the early stage of lactation, while the DNA methylation level in the dry period was lower than that at the early stage of lactation. Our result suggests that the DNA methylation of EEF1D plays an important role in the spatial and temporal regulation of its expression and possibly have an effect on the milk production traits.

Evolution of microRNA in primates

MicroRNA play an important role in post-transcriptional regulation of most transcripts in the human genome, but their evolution across the primate lineage is largely uncharacterized. A particular miRNA can have one to thousands of messenger RNA targets, establishing the potential for a small change in sequence or overall miRNA structure to have profound phenotypic effects. However, the majority of non-human primate miRNA is predicted solely by homology to the human genome and lacks experimental validation. In the present study, we sequenced thirteen species representing a wide range of the primate phylogeny. Hundreds of miRNA were validated, and the number of species with experimentally validated miRNA was tripled. These species include a sister taxon to humans (bonobo) and basal primates (aye-aye, mouse lemur, galago). Consistent with previous studies, we found the seed region and mature miRNA to be highly conserved across primates, with overall structural conservation of the pre-miRNA hairpin. However, there were a number of interesting exceptions, including a seed shift due to structural changes in miR-501. We also identified an increase in the number of miR-320 paralogs throughout primate evolution. Many of these non-conserved miRNA appear to regulate neuronal processes, illustrating the importance of investigating miRNA to learn more about human evolution.

Gene expression and adaptive noncoding changes during human evolution

Background

Despite evidence for adaptive changes in both gene expression and non-protein-coding, putatively regulatory regions of the genome during human evolution, the relationship between gene expression and adaptive changes in cis-regulatory regions remains unclear.

Results

Here we present new measurements of gene expression in five tissues of humans and chimpanzees, and use them to assess this relationship. We then compare our results with previous studies of adaptive noncoding changes, analyzing correlations at the level of gene ontology groups, in order to gain statistical power to detect correlations.

Conclusions

Consistent with previous studies, we find little correlation between gene expression and adaptive noncoding changes at the level of individual genes; however, we do find significant correlations at the level of biological function ontology groups. The types of function include processes regulated by specific transcription factors, responses to genetic or chemical perturbations, and differentiation of cell types within the immune system. Among functional categories co-enriched with both differential expression and noncoding adaptation, prominent themes include cancer, particularly epithelial cancers, and neural development and function.

Electronic supplementary material

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

Genome-wide DNA methylation at birth in relation to in utero arsenic exposure and the associated health in later life

BACKGROUND

In utero arsenic exposure may alter fetal developmental programming by altering DNA methylation, which may result in a higher risk of disease in later life. We evaluated the association between in utero arsenic exposure and DNA methylation (DNAm) in cord blood and its influence in later life.

METHODS

Genome-wide DNA methylation in cord blood from 64 subjects in the Taiwanese maternal infant and birth cohort was analyzed. Robust regressions were applied to assess the association of DNA methylation with in utero arsenic exposure. Multiple testing was adjusted by controlling false discovery rate (FDR) of 0.05. The DAVID bioinformatics tool was implemented for functional annotation analyses on the detected CpGs. The identified CpGs were further tested in an independent cohort. For the CpGs replicated in the independent cohort, linear mixed models were applied to assess the association of DNA methylation with low-density lipoprotein (LDL) at different ages (2, 5, 8, 11 and 14 years).

RESULTS

In total, 579 out of 385,183 CpGs were identified after adjusting for multiple testing (FDR = 0.05), of which ~60% were positively associated with arsenic exposure. Functional annotation analysis on these CpGs detected 17 KEGG pathways (FDR = 0.05) including pathways for cardiovascular diseases (CVD) and diabetes mellitus. In the independent cohort, about 46% (252 out of 553 CpGs) of the identified CpGs showed associations consistent with those in the study cohort. In total, 11 CpGs replicated in the independent cohort were in the pathways related to CVD and diabetes mellitus. Via longitudinal analyses, we found at 5 out of the 11 CpGs methylation was associated with LDL over time and interactions between DNA methylation and time were observed at 4 of the 5 CpGs, cg25189764 (coeff = 0.157, p-value = 0.047), cg04986899 (coeff. For interaction [coeff.int] = 0.030, p-value = 0.024), cg04903360 (coeff.int = 0.026, p-value = 0.032), cg08198265 (coeff.int = -0.063, p-value = 0.0021), cg10473311 (coeff.int = -0.021, p-value = 0.027).

CONCLUSION

In utero arsenic exposure was associated with cord blood DNA methylation at various CpGs. The identified CpGs may help determine pathological epigenetic mechanisms linked to in utero arsenic exposure. Five CpGs (cg25189764, cg04986899, cg04903360, cg08198265 and cg10473311) may serve as epigenetic markers for changes in LDL later in life.

DNA Methylation Changes in Intron 1 of Triggering Receptor Expressed on Myeloid Cell 2 in Japanese Schizophrenia Subjects

A hypothesis for schizophrenia (SCZ) called the “microglia hypothesis” has been suggested. In SCZ, expression of triggering receptor expressed on myeloid cell 2 (TREM2) mRNA is higher in leukocytes than in healthy individuals. Here, the methylation rates of four CpG sites in TREM2 intron 1 that may bind important transcription factors and the correlation between the methylation rate and mRNA expression were determined. We compared the methylation rates in SCZ patients and age-matched controls (n = 50 each). SCZ patients had significantly lower methylation rates of CpG 2 (17.0 ± 6.7 vs. 20.2 ± 5.0; p = 0.02) and CpG 3 (23.8 ± 8.2 vs. 28.1 ± 6.2; p = 0.01). The average methylation rate (15.3 ± 5.2 vs. 17.6 ± 3.9; p = 0.009) was also lower. A significant negative correlation was found between TREM2 mRNA expression and the methylation rate of CpG 2 (r = −0.252, p = 0.012). SCZ susceptibility markers may include low methylation at TREM2 intron 1 and increased TREM2 mRNA levels. Our pilot study requires validation with higher numbers of participants and with other myeloid cell types.

Prenatal Maternal Stress Predicts Methylation of Genes Regulating the Hypothalamic-Pituitary-Adrenocortical System in Mothers and Newborns in the Democratic Republic of Congo

Exposure to stress early in life permanently shapes activity of the hypothalamic-pituitary-adrenocortical (HPA) axis and the brain. Prenatally, glucocorticoids pass through the placenta to the fetus with postnatal impacts on brain development, birth weight (BW), and HPA axis functioning. Little is known about the biological mechanisms by which prenatal stress affects postnatal functioning. This study addresses this gap by examining the effect of chronic stress and traumatic war-related stress on epigenetic changes in four key genes regulating the HPA axis in neonatal cord blood, placenta, and maternal blood: CRH, CRHBP, NR3C1, and FKBP5. Participants were 24 mother-newborn dyads in the conflict-ridden region of the eastern Democratic Republic of Congo. BW data were collected at delivery and maternal interviews were conducted to assess culturally relevant chronic and war-related stressors. Chronic stress and war trauma had widespread effects on HPA axis gene methylation, with significant effects observed at transcription factor binding (TFB) sites in all target genes tested. Some changes in methylation were unique to chronic or war stress, whereas others were observed across both stressor types. Moreover, stress exposures impacted maternal and fetal tissues differently, supporting theoretical models that stress impacts vary according to life phase. Methylation in several NR3C1 and CRH CpG sites, all located at TFB sites, was associated with BW. These findings suggest that prenatal stress exposure impacts development via epigenetic changes in HPA axis genes.

Low methylation rates of dopamine receptor D2 gene promoter sites in Japanese schizophrenia subjects

OBJECTIVES

According to the dopamine hypothesis, several studies on the gene for the dopamine receptor D2 (DRD2) have been conducted. However, no trait biomarkers on DRD2 are available. We examined whether the methylation rates in the upstream region of DRD2 in leukocytes are different in schizophrenia (SCZ) subjects compared to control subjects.

METHODS

We selected seven CpG sites in the upstream region of DRD2 that may theoretically bind major transcription factors. The methylation rates in these regions of 50 medicated and 18 drug-naïve SCZ subjects were compared with those of age-matched control subjects.

RESULTS

The methylation rates were significantly lower in medicated (CpG2, P < 0.0001; CpG4, P = 0.013; CpG7, P < 0.0001; and average: 12.9 ± 1.8 vs. 14.1 ± 2.2, P = 0.005) and drug-naïve SCZ subjects (CpG1, P = 0.006; CpG2, P = 0.001; CpG3, P = 0.001; CpG5, P = 0.02; CpG6, P = 0.015; CpG7, P = 0.027; and average: 9.86 ± 0.9 vs. 11.2 ± 1.3, P = 0.002).

CONCLUSIONS

We confirmed low methylation rates in the upstream region of DRD2 in both medicated and drug-naïve SCZ subjects. Low methylation rates of DRD2 in leukocytes may be a trait biomarker for SCZ.

Comparative Methylome Analyses Identify Epigenetic Regulatory Loci of Human Brain Evolution

How do epigenetic modifications change across species and how do these modifications affect evolution? These are fundamental questions at the forefront of our evolutionary epigenomic understanding. Our previous work investigated human and chimpanzee brain methylomes, but it was limited by the lack of outgroup data which is critical for comparative (epi)genomic studies. Here, we compared whole genome DNA methylation maps from brains of humans, chimpanzees and also rhesus macaques (outgroup) to elucidate DNA methylation changes during human brain evolution. Moreover, we validated that our approach is highly robust by further examining 38 human-specific DMRs using targeted deep genomic and bisulfite sequencing in an independent panel of 37 individuals from five primate species. Our unbiased genome-scan identified human brain differentially methylated regions (DMRs), irrespective of their associations with annotated genes. Remarkably, over half of the newly identified DMRs locate in intergenic regions or gene bodies. Nevertheless, their regulatory potential is on par with those of promoter DMRs. An intriguing observation is that DMRs are enriched in active chromatin loops, suggesting human-specific evolutionary remodeling at a higher-order chromatin structure. These findings indicate that there is substantial reprogramming of epigenomic landscapes during human brain evolution involving noncoding regions.

Identification of Multiple Forms of RNA Transcripts Associated with Human-Specific Retrotransposed Gene Copies

The human genome contains thousands of retrocopies, mostly as processed pseudogenes, which were recently shown to be prevalently transcribed. In particular, those specifically acquired in the human lineage are able to modulate gene expression in a manner that contributed to the evolution of human-specific traits. Therefore, knowledge of the human-specific retrocopies that are transcribed or their full-length transcript structure contributes to better understand human genome evolution. In this study, we identified 16 human-specific retrocopies that harbor 5' CpG islands by in silico analysis and showed that 12 were transcribed in normal tissues and cancer cell lines with a variety of expression patterns, including cancer-specific expression. Determination of the structure of the transcripts associated with the retrocopies revealed that none were transcribed from their 5' CpG islands, but rather, from inside the 3' UTR and the nearby 5' flanking region of the retrocopies as well as the promoter of neighboring genes. The multiple forms of the transcripts, such as chimeric and individual transcripts in both the sense and antisense orientation, might have introduced novel post-transcriptional regulation into the genome during human evolution. These results shed light on the potential role of human-specific retrocopies in the evolution of gene regulation and genomic disorders.

Molecular characterization and analysis of the porcine NURR1 gene

Orphan receptor NURR1 (also termed NR4A2) belongs to the nuclear receptor superfamily and functions as a regulatory factor of differentiation, migration, maturation and maintenance of mesencephalic dopaminergic neurons. NURR1 plays an important role in nigrostriatal dopamine neuron development and is therefore implicated in the pathogenesis of neurodegenerative diseases linked to the dopamine system of the midbrain.

Here we report the isolation and characterization of porcine NURR1 cDNA. The NURR1 cDNA was RT-PCR cloned using NURR1-specific oligonucleotide primers derived from in silico sequences. The porcine NURR1 cDNA encodes a polypeptide of 598 amino acids, displaying a very high similarity with bovine, human and mouse (99%) NURR1 protein. Expression analysis revealed a differential NURR1 mRNA expression in various organs and tissues. NURR1 transcripts could be detected as early as at 60 days of embryo development in different brain tissues. A significant increase in NURR1 transcript in the cerebellum and a decrease in NURR1 transcript in the basal ganglia was observed during embryo development. The porcine NURR1 gene was mapped to chromosome 15. Two missense mutations were found in exon 3, the first coding exon of NURR1. Methylation analysis of the porcine NURR1 gene body revealed a high methylation degree in brain tissue, whereas methylation of the promoter was very low. A decrease in DNA methylation in a discrete region of the NURR1 promoter was observed in pig frontal cortex during pig embryo development. This observation correlated with an increase in NURR1 transcripts. Therefore, methylation might be a determinant of NURR1 expression at certain time points in embryo development.

•

The porcine NURR1 gene was cloned and characterized.

•

NURR1 transcript was detected early in pig embryo brain development.

•

Methylation status of NURR1 may be a determinant for its expression.

Human age estimation from blood using mRNA, DNA methylation, DNA rearrangement, and telomere length

Establishing the age of unknown persons, or persons with unknown age, can provide important leads in police investigations, disaster victim identification, fraud cases, and in other legal affairs. Previous methods mostly relied on morphological features available from teeth or skeletal parts. The development of molecular methods for age estimation allowing to use human specimens that possess no morphological age information, such as bloodstains, is extremely valuable as this type of samples is commonly found at crime scenes. Recently, we introduced a DNA-based approach for human age estimation from blood based on the quantification of T-cell specific DNA rearrangements (sjTRECs), which achieves accurate assignment of blood DNA samples to one of four 20-year-interval age categories. Aiming at improving the accuracy of molecular age estimation from blood, we investigated different types of biomarkers. We started out by systematic genome-wide surveys for new age-informative mRNA and DNA methylation markers in blood from the same young and old individuals using microarray technologies. The obtained candidate markers were validated in independent samples covering a wide age range using alternative technologies together with previously proposed DNA methylation, sjTREC, and telomere length markers. Cross-validated multiple regression analysis was applied for estimating and validating the age predictive power of various sets of biomarkers within and across different marker types. We found that DNA methylation markers outperformed mRNA, sjTREC, and telomere length in age predictive power. The best performing model included 8 DNA methylation markers derived from 3 CpG islands reaching a high level of accuracy (cross-validated R(2)=0.88, SE±6.97 years, mean absolute deviation 5.07 years). However, our data also suggest that mRNA markers can provide independent age information: a model using a combined set of 5 DNA methylation markers and one mRNA marker could provide similarly high accuracy (cross-validated R(2)=0.86, SE±7.62 years, mean absolute deviation 4.60 years). Overall, our study provides new and confirms previously suggested molecular biomarkers for age estimation from blood. Moreover, our comparative study design revealed that DNA methylation markers are superior for this purpose over other types of molecular biomarkers tested. While the new and some previous findings are highly promising, before molecular age estimation can eventually meet forensic practice, the proposed biomarkers should be tested further in larger sets of blood samples from both healthy and unhealthy individuals, and markers and genotyping methods shall be validated to meet forensic standards.

Evolution of Epigenetic Regulation in Vertebrate Genomes

Empirical models of sequence evolution have spurred progress in the field of evolutionary genetics for decades. We are now realizing the importance and complexity of the eukaryotic epigenome. While epigenome analysis has been applied to genomes from single-cell eukaryotes to human, comparative analyses are still relatively few and computational algorithms to quantify epigenome evolution remain scarce. Accordingly, a quantitative model of epigenome evolution remains to be established. We review here the comparative epigenomics literature and synthesize its overarching themes. We also suggest one mechanism, transcription factor binding site (TFBS) turnover, which relates sequence evolution to epigenetic conservation or divergence. Lastly, we propose a framework for how the field can move forward to build a coherent quantitative model of epigenome evolution.

Linking short tandem repeat polymorphisms with cytosine modifications in human lymphoblastoid cell lines

Inter-individual variation in cytosine modifications has been linked to complex traits in humans. Cytosine modification variation is partially controlled by single nucleotide polymorphisms (SNPs), known as modified cytosine quantitative trait loci (mQTL). However, little is known about the role of short tandem repeat polymorphisms (STRPs), a class of structural genetic variants, in regulating cytosine modifications. Utilizing the published data on the International HapMap Project lymphoblastoid cell lines (LCLs), we assessed the relationships between 721 STRPs and the modification levels of 283,540 autosomal CpG sites. Our findings suggest that, in contrast to the predominant cis-acting mode for SNP-based mQTL, STRPs are associated with cytosine modification levels in both cis-acting (local) and trans-acting (distant) modes. In local scans within the ±1 Mb windows of target CpGs, 21, 9, and 21 cis-acting STRP-based mQTL were detected in CEU (Caucasian residents from Utah, USA), YRI (Yoruba people from Ibadan, Nigeria), and the combined samples, respectively. In contrast, 139,420, 76,817, and 121,866 trans-acting STRP-based mQTL were identified in CEU, YRI, and the combined samples, respectively. A substantial proportion of CpG sites detected with local STRP-based mQTL were not associated with SNP-based mQTL, suggesting that STRPs represent an independent class of mQTL. Functionally, genetic variants neighboring CpG-associated STRPs are enriched with genome-wide association study (GWAS) loci for a variety of complex traits and diseases, including cancers, based on the National Human Genome Research Institute (NHGRI) GWAS Catalog. Therefore, elucidating these STRP-based mQTL in addition to SNP-based mQTL can provide novel insights into the genetic architectures of complex traits.

A Study of Alterations in DNA Epigenetic Modifications (5mC and 5hmC) and Gene Expression Influenced by Simulated Microgravity in Human Lymphoblastoid Cells

Cells alter their gene expression in response to exposure to various environmental changes. Epigenetic mechanisms such as DNA methylation are believed to regulate the alterations in gene expression patterns. In vitro and in vivo studies have documented changes in cellular proliferation, cytoskeletal remodeling, signal transduction, bone mineralization and immune deficiency under the influence of microgravity conditions experienced in space. However microgravity induced changes in the epigenome have not been well characterized. In this study we have used Next-generation Sequencing (NGS) to profile ground-based “simulated” microgravity induced changes on DNA methylation (5-methylcytosine or 5mC), hydroxymethylation (5-hydroxymethylcytosine or 5hmC), and simultaneous gene expression in cultured human lymphoblastoid cells. Our results indicate that simulated microgravity induced alterations in the methylome (~60% of the differentially methylated regions or DMRs are hypomethylated and ~92% of the differentially hydroxymethylated regions or DHMRs are hyperhydroxymethylated). Simulated microgravity also induced differential expression in 370 transcripts that were associated with crucial biological processes such as oxidative stress response, carbohydrate metabolism and regulation of transcription. While we were not able to obtain any global trend correlating the changes of methylation/ hydroxylation with gene expression, we have been able to profile the simulated microgravity induced changes of 5mC over some of the differentially expressed genes that includes five genes undergoing differential methylation over their promoters and twenty five genes undergoing differential methylation over their gene-bodies. To the best of our knowledge, this is the first NGS-based study to profile epigenomic patterns induced by short time exposure of simulated microgravity and we believe that our findings can be a valuable resource for future explorations.

DNA Methylation: Insights into Human Evolution

A fundamental initiative for evolutionary biologists is to understand the molecular basis underlying phenotypic diversity. A long-standing hypothesis states that species-specific traits may be explained by differences in gene regulation rather than differences at the protein level. Over the past few years, evolutionary studies have shifted from mere sequence comparisons to integrative analyses in which gene regulation is key to understanding species evolution. DNA methylation is an important epigenetic modification involved in the regulation of numerous biological processes. Nevertheless, the evolution of the human methylome and the processes driving such changes are poorly understood. Here, we review the close interplay between Cytosine-phosphate-Guanine (CpG) methylation and the underlying genome sequence, as well as its evolutionary impact. We also summarize the latest advances in the field, revisiting the main literature on human and nonhuman primates. We hope to encourage the scientific community to address the many challenges posed by the field of comparative epigenomics.

Multivariate Boosting for Integrative Analysis of High-Dimensional Cancer Genomic Data

In this paper, we propose a novel multivariate component-wise boosting method for fitting multivariate response regression models under the high-dimension, low sample size setting. Our method is motivated by modeling the association among different biological molecules based on multiple types of high-dimensional genomic data. Particularly, we are interested in two applications: studying the influence of DNA copy number alterations on RNA transcript levels and investigating the association between DNA methylation and gene expression. For this purpose, we model the dependence of the RNA expression levels on DNA copy number alterations and the dependence of gene expression on DNA methylation through multivariate regression models and utilize boosting-type method to handle the high dimensionality as well as model the possible nonlinear associations. The performance of the proposed method is demonstrated through simulation studies. Finally, our multivariate boosting method is applied to two breast cancer studies.

Derivation of induced pluripotent stem cells from orangutan skin fibroblasts

Background

Orangutans are an endangered species whose natural habitats are restricted to the Southeast Asian islands of Borneo and Sumatra. Along with the African great apes, orangutans are among the closest living relatives to humans. For potential species conservation and functional genomics studies, we derived induced pluripotent stem cells (iPSCs) from cryopreserved somatic cells obtained from captive orangutans.

Results

Primary skin fibroblasts from two Sumatran orangutans were transduced with retroviral vectors expressing the human OCT4, SOX2, KLF4, and c-MYC factors. Candidate orangutan iPSCs were characterized by global gene expression and DNA copy number analysis. All were consistent with pluripotency and provided no evidence of large genomic insertions or deletions. In addition, orangutan iPSCs were capable of producing cells derived from all three germ layers in vitro through embryoid body differentiation assays and in vivo through teratoma formation in immune-compromised mice.

Conclusions

We demonstrate that orangutan skin fibroblasts are capable of being reprogrammed into iPSCs with hallmark molecular signatures and differentiation potential. We suggest that reprogramming orangutan somatic cells in genome resource banks could provide new opportunities for advancing assisted reproductive technologies relevant for species conservation efforts. Furthermore, orangutan iPSCs could have applications for investigating the phenotypic relevance of genomic changes that occurred in the human, African great ape, and/or orangutan lineages. This provides opportunities for orangutan cell culture models that would otherwise be impossible to develop from living donors due to the invasive nature of the procedures required for obtaining primary cells.

Electronic supplementary material

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

Methadone induces hypermethylation of human DNA

BACKGROUND

Increased global DNA methylation in the blood of patients chronically exposed to opioids had been interpreted as an indication of an epigenetic action of this drug class.

MATERIALS & METHODS

To strengthen the causality, human MCF7 cells were cultured in media with the addition of several known or potential modulators of DNA methylation including methadone.

RESULTS

Following 3 days of incubation with several different known or potential epigenetic modulators, global DNA methylation, quantified at LINE-1 CpG islands, showed a large variability across all treatments ranging from 27.8 to 63%. Based on distribution analysis of the global methylation of human DNA exposed to various potential modulators, present in vitro experiments showed that treatment with the opioid methadone was associated with an increased probability of hypermethylation.

CONCLUSION

This strengthens the evidence that opioids interfere with mechanisms of classical epigenetics.

Characterizing 5-hydroxymethylcytosine in human prefrontal cortex at single base resolution

Background

The recent discovery that methylated cytosines are converted to 5-hydroxymethylated cytosines (5hmC) by the family of ten-eleven translocation enzymes has sparked significant interest on the genomic location, the abundance in different tissues, the putative functions, and the stability of this epigenetic mark. 5hmC plays a key role in the brain, where it is particularly abundant and dynamic during development.

Results

Here, we comprehensively characterize 5hmC in the prefrontal cortices of 24 subjects. We show that, although there is inter-individual variability in 5hmC content among unrelated individuals, approximately 8 % of all CpGs on autosomal chromosomes contain 5hmC, while sex chromosomes contain far less. Our data also provide evidence suggesting that 5hmC has transcriptional regulatory properties, as the density of 5hmC was highest in enhancer regions and within exons. Furthermore, we link increased 5hmC density to histone modification binding sites, to the gene bodies of actively transcribed genes, and to exon-intron boundaries. Finally, we provide several genomic regions of interest that contain gender-specific 5hmC.

Conclusions

Collectively, these results present an important reference for the growing number of studies that are interested in the investigation of the role of 5hmC in brain and mental disorders.

Electronic supplementary material

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

Epigenetic modifications are associated with inter-species gene expression variation in primates

BACKGROUND

Changes in gene regulation have long been thought to play an important role in evolution and speciation, especially in primates. Over the past decade, comparative genomic studies have revealed extensive inter-species differences in gene expression levels, yet we know much less about the extent to which regulatory mechanisms differ between species.

RESULTS

To begin addressing this gap, we perform a comparative epigenetic study in primate lymphoblastoid cell lines, to query the contribution of RNA polymerase II and four histone modifications, H3K4me1, H3K4me3, H3K27ac, and H3K27me3, to inter-species variation in gene expression levels. We find that inter-species differences in mark enrichment near transcription start sites are significantly more often associated with inter-species differences in the corresponding gene expression level than expected by chance alone. Interestingly, we also find that first-order interactions among the five marks, as well as chromatin states, do not markedly contribute to the degree of association between the marks and inter-species variation in gene expression levels, suggesting that the marginal effects of the five marks dominate this contribution.

CONCLUSIONS

Our observations suggest that epigenetic modifications are substantially associated with changes in gene expression levels among primates and may represent important molecular mechanisms in primate evolution.

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.

Layered genetic control of DNA methylation and gene expression: a locus of multiple sclerosis in healthy individuals

DNA methylation may contribute to the etiology of complex genetic disorders through its impact on genome integrity and gene expression; it is modulated by DNA-sequence variants, named methylation quantitative trait loci (meQTLs). Most meQTLs influence methylation of a few CpG dinucleotides within short genomic regions (<3 kb). Here, we identified a layered genetic control of DNA methylation at numerous CpGs across a long 300 kb genomic region. This control involved a single long-range meQTL and multiple local meQTLs. The long-range meQTL explained up to 75% of variance in methylation of CpGs located over extended areas of the 300 kb region. The meQTL was identified in four samples (P = 2.8 × 10(-17), 3.1 × 10(-31), 4.0 × 10(-71) and 5.2 × 10(-199)), comprising a total of 2796 individuals. The long-range meQTL was strongly associated not only with DNA methylation but also with mRNA expression of several genes within the 300 kb region (P = 7.1 × 10(-18)-1.0 × 10(-123)). The associations of the meQTL with gene expression became attenuated when adjusted for DNA methylation (causal inference test: P = 2.4 × 10(-13)-7.1 × 10(-20)), indicating coordinated regulation of DNA methylation and gene expression. Further, the long-range meQTL was found to be in linkage disequilibrium with the most replicated locus of multiple sclerosis, a disease affecting primarily the brain white matter. In middle-aged adults free of the disease, we observed that the risk allele was associated with subtle structural properties of the brain white matter found in multiple sclerosis (P = 0.02). In summary, we identified a long-range meQTL that controls methylation and expression of several genes and may be involved in increasing brain vulnerability to multiple sclerosis.

The interplay between DNA methylation and sequence divergence in recent human evolution

Despite the increasing knowledge about DNA methylation, the understanding of human epigenome evolution is in its infancy. Using whole genome bisulfite sequencing we identified hundreds of differentially methylated regions (DMRs) in humans compared to non-human primates and estimated that ∼25% of these regions were detectable throughout several human tissues. Human DMRs were enriched for specific histone modifications and the majority were located distal to transcription start sites, highlighting the importance of regions outside the direct regulatory context. We also found a significant excess of endogenous retrovirus elements in human-specific hypomethylated.

We reported for the first time a close interplay between inter-species genetic and epigenetic variation in regions of incomplete lineage sorting, transcription factor binding sites and human differentially hypermethylated regions. Specifically, we observed an excess of human-specific substitutions in transcription factor binding sites located within human DMRs, suggesting that alteration of regulatory motifs underlies some human-specific methylation patterns. We also found that the acquisition of DNA hypermethylation in the human lineage is frequently coupled with a rapid evolution at nucleotide level in the neighborhood of these CpG sites. Taken together, our results reveal new insights into the mechanistic basis of human-specific DNA methylation patterns and the interpretation of inter-species non-coding variation.

Genome-wide mapping reveals conservation of promoter DNA methylation following chicken domestication

It is well-known that environment influences DNA methylation, however, the extent of heritable DNA methylation variation following animal domestication remains largely unknown. Using meDIP-chip we mapped the promoter methylomes for 23,316 genes in muscle tissues of ancestral and domestic chickens. We systematically examined the variation of promoter DNA methylation in terms of different breeds, differentially expressed genes, SNPs and genes undergo genetic selection sweeps. While considerable changes in DNA sequence and gene expression programs were prevalent, we found that the inter-strain DNA methylation patterns were highly conserved in promoter region between the wild and domestic chicken breeds. Our data suggests a global preservation of DNA methylation between the wild and domestic chicken breeds in either a genome-wide or locus-specific scale in chick muscle tissues.

Extrapolating histone marks across developmental stages, tissues, and species: an enhancer prediction case study

Background

Dynamic activation and inactivation of gene regulatory DNA produce the expression changes that drive the differentiation of cellular lineages. Identifying regulatory regions active during developmental transitions is necessary to understand how the genome specifies complex developmental programs and how these processes are disrupted in disease. Gene regulatory dynamics are mediated by many factors, including the binding of transcription factors (TFs) and the methylation and acetylation of DNA and histones. Genome-wide maps of TF binding and DNA and histone modifications have been generated for many cellular contexts; however, given the diversity and complexity of animal development, these data cover only a small fraction of the cellular and developmental contexts of interest. Thus, there is a need for methods that use existing epigenetic and functional genomics data to analyze the thousands of contexts that remain uncharacterized.

Results

To investigate the utility of histone modification data in the analysis of cellular contexts without such data, I evaluated how well genome-wide H3K27ac and H3K4me1 data collected in different developmental stages, tissues, and species were able to predict experimentally validated heart enhancers active at embryonic day 11.5 (E11.5) in mouse. Using a machine-learning approach to integrate the data from different contexts, I found that E11.5 heart enhancers can often be predicted accurately from data from other contexts, and I quantified the contribution of each data source to the predictions. The utility of each dataset correlated with nearness in developmental time and tissue to the target context: data from late developmental stages and adult heart tissues were most informative for predicting E11.5 enhancers, while marks from stem cells and early developmental stages were less informative. Predictions based on data collected in non-heart tissues and in human hearts were better than random, but worse than using data from mouse hearts.

Conclusions

The ability of these algorithms to accurately predict developmental enhancers based on data from related, but distinct, cellular contexts suggests that combining computational models with epigenetic data sampled from relevant contexts may be sufficient to enable functional characterization of many cellular contexts of interest.

Electronic supplementary material

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

The genetic and mechanistic basis for variation in gene regulation

It is now well established that noncoding regulatory variants play a central role in the genetics of common diseases and in evolution. However, until recently, we have known little about the mechanisms by which most regulatory variants act. For instance, what types of functional elements in DNA, RNA, or proteins are most often affected by regulatory variants? Which stages of gene regulation are typically altered? How can we predict which variants are most likely to impact regulation in a given cell type? Recent studies, in many cases using quantitative trait loci (QTL)-mapping approaches in cell lines or tissue samples, have provided us with considerable insight into the properties of genetic loci that have regulatory roles. Such studies have uncovered novel biochemical regulatory interactions and led to the identification of previously unrecognized regulatory mechanisms. We have learned that genetic variation is often directly associated with variation in regulatory activities (namely, we can map regulatory QTLs, not just expression QTLs [eQTLs]), and we have taken the first steps towards understanding the causal order of regulatory events (for example, the role of pioneer transcription factors). Yet, in most cases, we still do not know how to interpret overlapping combinations of regulatory interactions, and we are still far from being able to predict how variation in regulatory mechanisms is propagated through a chain of interactions to eventually result in changes in gene expression profiles.

Evolutionary dynamics of coding and non-coding transcriptomes

Gene expression changes may underlie much of phenotypic evolution. The development of high-throughput RNA sequencing protocols has opened the door to unprecedented large-scale and cross-species transcriptome comparisons by allowing accurate and sensitive assessments of transcript sequences and expression levels. Here, we review the initial wave of the new generation of comparative transcriptomic studies in mammals and vertebrate outgroup species in the context of earlier work. Together with various large-scale genomic and epigenomic data, these studies have unveiled commonalities and differences in the dynamics of gene expression evolution for various types of coding and non-coding genes across mammalian lineages, organs, developmental stages, chromosomes and sexes. They have also provided intriguing new clues to the regulatory basis and phenotypic implications of evolutionary gene expression changes.

Dysregulation of genome-wide gene expression and DNA methylation in abnormal cloned piglets

Background

Epigenetic modifications (especially altered DNA methylation) resulting in altered gene expression may be one reason for development failure or abnormalities in cloned animals, but the underlying mechanism of the abnormal phenotype in cloned piglets remains unknown. Some cloned piglets in our study showed abnormal phenotypes such as large tongue (longer and thicker), weak muscles, and exomphalos. Here we conducted DNA methylation (DNAm) immunoprecipitation and high throughput sequencing (MeDIP-seq) and RNA sequencing (RNA-seq) of muscle tissues of cloned piglets to investigate the relationship of abnormal DNAm with gene dysregulation and the unusual phenotypes in cloned piglets.

Results

Analysis of the methylomes revealed that abnormal cloned piglets suffered more hypomethylation than hypermethylation compared to the normal cloned piglets, although the DNAm level in the CpG Island was higher in the abnormal cloned piglets. Some repetitive elements, such as SINE/tRNA-Glu Satellite/centr also showed differences. We detected 1,711 differentially expressed genes (DEGs) between the two groups, of which 243 genes also changed methylation level in the abnormal cloned piglets. The altered DNA methylation mainly affected the low and silently expressed genes. There were differences in both pathways and genes, such as the MAPK signalling pathway, the hypertrophic cardiomyopathy pathway, and the imprinted gene PLAGL1; all of which may play important roles in development of the abnormal phenotype.

Conclusions

The abnormal cloned piglets showed substantial changes both in the DNAm and the gene expression. Our data may provide new insights into understanding the molecular mechanisms of the reprogramming of genetic information in cloned animals.

Electronic supplementary material

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

Comparative studies of gene regulatory mechanisms

It has become increasingly clear that changes in gene regulation have played an important role in adaptive evolution both between and within species. Over the past five years, comparative studies have moved beyond simple characterizations of differences in gene expression levels within and between species to studying variation in regulatory mechanisms. We still know relatively little about the precise chain of events that lead to most regulatory adaptations, but we have taken significant steps towards understanding the relative importance of changes in different mechanisms of gene regulatory evolution. In this review, we first discuss insights from comparative studies in model organisms, where the available experimental toolkit is extensive. We then focus on a few recent comparative studies in primates, where the limited feasibility of experimental manipulation dictates the approaches that can be used to study gene regulatory evolution.

DNA methylation age of human tissues and cell types

Background

It is not yet known whether DNA methylation levels can be used to accurately predict age across a broad spectrum of human tissues and cell types, nor whether the resulting age prediction is a biologically meaningful measure.

Results

I developed a multi-tissue predictor of age that allows one to estimate the DNA methylation age of most tissues and cell types. The predictor, which is freely available, was developed using 8,000 samples from 82 Illumina DNA methylation array datasets, encompassing 51 healthy tissues and cell types. I found that DNA methylation age has the following properties: first, it is close to zero for embryonic and induced pluripotent stem cells; second, it correlates with cell passage number; third, it gives rise to a highly heritable measure of age acceleration; and, fourth, it is applicable to chimpanzee tissues. Analysis of 6,000 cancer samples from 32 datasets showed that all of the considered 20 cancer types exhibit significant age acceleration, with an average of 36 years. Low age-acceleration of cancer tissue is associated with a high number of somatic mutations and TP53 mutations, while mutations in steroid receptors greatly accelerate DNA methylation age in breast cancer. Finally, I characterize the 353 CpG sites that together form an aging clock in terms of chromatin states and tissue variance.

Conclusions

I propose that DNA methylation age measures the cumulative effect of an epigenetic maintenance system. This novel epigenetic clock can be used to address a host of questions in developmental biology, cancer and aging research.

HOXA11 gene is hypermethylation and aberrant expression in gastric cancer

Background

Aberrant DNA methylation is an acquired epigenetic alteration that serves as an alternative to genetic defects in the inactivation of tumor suppressor genes and other genes in diverse human cancers. Gastric carcinoma is one of the tumors with a high frequency of aberrant methylation in promoter region. Hence we investigated the promoter methylation status and expression level of HOXA11 gene which may involve in GC development.

Methods

Thirty-two surgical excised gastric cancer specimens, twelve paired adjacent non-cancerous specimens and seven normal gastric mucosas were examined. The methylation status and expression level of HOXA11 gene were determined by bisulfite sequencing polymerase chain reaction (BSP), real-time polymerase chain reaction (RT-PCR) and immunohistochemistry (IHC) respectively. HOXA11 expression was knocked-down with siRNA to mimic HOXA11 gene hypermethylation and ability of cell proliferation and migration was determinate. In addition, we analyzed and correlated the findings with clinicopathological features.

Results

The methylation level of HOXA11 gene in gastric cancer tissues and adjacent non-cancerous tissues were higher than those in normal gastric mucosa (P < 0.05). The methylation level was higher in TNM III and IV patients of GC than those in TNM I and II patients (P < 0.05). The expression of HOXA11 mRNA and protein decreased in normal gastric mucosa, peri-cancer tissue and GC (P < 0.05). HOXA11 expression was inversely correlated with DNA methylation (P < 0.05). Knocked-down of HOXA11 expression with siRNA in BGC-823 cells enhanced cell proliferation compared with control, but no significant different was observed in migration ability.

Conclusion

Hypermethylation and decreased expression of HOXA11 gene may be involved in the carcinogenesis and development of GC and may provide useful information for the prediction of the malignant behaviors of GC. And the expression of HOXA11 is impaired by DNA methylation. However, repression of HOXA11 expression promoted BGC-823 cell proliferation.

Whole-genome bisulfite sequencing of multiple individuals reveals complementary roles of promoter and gene body methylation in transcriptional regulation

BACKGROUND

DNA methylation is an important type of epigenetic modification involved in gene regulation. Although strong DNA methylation at promoters is widely recognized to be associated with transcriptional repression, many aspects of DNA methylation remain not fully understood, including the quantitative relationships between DNA methylation and expression levels, and the individual roles of promoter and gene body methylation.

RESULTS

Here we present an integrated analysis of whole-genome bisulfite sequencing and RNA sequencing data from human samples and cell lines. We find that while promoter methylation inversely correlates with gene expression as generally observed, the repressive effect is clear only on genes with a very high DNA methylation level. By means of statistical modeling, we find that DNA methylation is indicative of the expression class of a gene in general, but gene body methylation is a better indicator than promoter methylation. These findings are general in that a model constructed from a sample or cell line could accurately fit the unseen data from another. We further find that promoter and gene body methylation have minimal redundancy, and either one is sufficient to signify low expression. Finally, we obtain increased modeling power by integrating histone modification data with the DNA methylation data, showing that neither type of information fully subsumes the other.

CONCLUSION

Our results suggest that DNA methylation outside promoters also plays critical roles in gene regulation. Future studies on gene regulatory mechanisms and disease-associated differential methylation should pay more attention to DNA methylation at gene bodies and other non-promoter regions.

Whole-genome bisulfite sequencing of multiple individuals reveals complementary roles of promoter and gene body methylation in transcriptional regulation

Background

DNA methylation is an important type of epigenetic modification involved in gene regulation. Although strong DNA methylation at promoters is widely recognized to be associated with transcriptional repression, many aspects of DNA methylation remain not fully understood, including the quantitative relationships between DNA methylation and expression levels, and the individual roles of promoter and gene body methylation.

Results

Here we present an integrated analysis of whole-genome bisulfite sequencing and RNA sequencing data from human samples and cell lines. We find that while promoter methylation inversely correlates with gene expression as generally observed, the repressive effect is clear only on genes with a very high DNA methylation level. By means of statistical modeling, we find that DNA methylation is indicative of the expression class of a gene in general, but gene body methylation is a better indicator than promoter methylation. These findings are general in that a model constructed from a sample or cell line could accurately fit the unseen data from another. We further find that promoter and gene body methylation have minimal redundancy, and either one is sufficient to signify low expression. Finally, we obtain increased modeling power by integrating histone modification data with the DNA methylation data, showing that neither type of information fully subsumes the other.

Conclusion

Our results suggest that DNA methylation outside promoters also plays critical roles in gene regulation. Future studies on gene regulatory mechanisms and disease-associated differential methylation should pay more attention to DNA methylation at gene bodies and other non-promoter regions.

Electronic supplementary material

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

Stability of gene expression and epigenetic profiles highlights the utility of patient-derived paediatric acute lymphoblastic leukaemia xenografts for investigating molecular mechanisms of drug resistance

BACKGROUND

Patient-derived tumour xenografts are an attractive model for preclinical testing of anti-cancer drugs. Insights into tumour biology and biomarkers predictive of responses to chemotherapeutic drugs can also be gained from investigating xenograft models. As a first step towards examining the equivalence of epigenetic profiles between xenografts and primary tumours in paediatric leukaemia, we performed genome-scale DNA methylation and gene expression profiling on a panel of 10 paediatric B-cell precursor acute lymphoblastic leukaemia (BCP-ALL) tumours that were stratified by prednisolone response.

RESULTS

We found high correlations in DNA methylation and gene expression profiles between matching primary and xenograft tumour samples with Pearson's correlation coefficients ranging between 0.85 and 0.98. In order to demonstrate the potential utility of epigenetic analyses in BCP-ALL xenografts, we identified DNA methylation biomarkers that correlated with prednisolone responsiveness of the original tumour samples. Differential methylation of CAPS2, ARHGAP21, ARX and HOXB6 were confirmed by locus specific analysis. We identified 20 genes showing an inverse relationship between DNA methylation and gene expression in association with prednisolone response. Pathway analysis of these genes implicated apoptosis, cell signalling and cell structure networks in prednisolone responsiveness.

CONCLUSIONS

The findings of this study confirm the stability of epigenetic and gene expression profiles of paediatric BCP-ALL propagated in mouse xenograft models. Further, our preliminary investigation of prednisolone sensitivity highlights the utility of mouse xenograft models for preclinical development of novel drug regimens with parallel investigation of underlying gene expression and epigenetic responses associated with novel drug responses.

Epigenetics and evolution

Epigenetic mechanisms traditionally have been studied in the domains of development and disease, but they may also play important roles in ecological and evolutionary processes. In this article, we revisit historical as well as recent studies that indicate significant impacts of epigenetic processes on evolution. Our main focus is DNA methylation, which is a prevalent chemical modification of genomic DNA. First, it has been long known that DNA methylation acts as a major mutational facilitator in animal genomes and influences nucleotide compositions of genomes. More recently, genome-wide analyses have demonstrated that the current levels of DNA methylation can be predicted from the evolutionary signatures of DNA methylation, indicating that these two processes are intimately correlated. Indeed, the recent explosive growth in the knowledge of genomic DNA methylation in wide-ranging taxa has revealed that patterns of DNA methylation are surprisingly conserved across deep phylogenies. Interestingly, comparative analyses of humans and closely related primate species show that genomic regions that do show evolutionary divergence of DNA methylation are enriched for developmental and tissue specializations. A key question is how epigenetic patterns transmit between generations and impact evolutionary dynamics. On the one hand, some studies report direct transmissions of epigenetic features to the next generation. On the other hand, it is becoming clear that genomic sequence variants exist that encode and presumably regulate distinctive epigenetic patterns. For instance, numerous single-nucleotide polymorphisms that affect DNA-methylation patterns have been discovered in human populations. These studies begin to unveil a dynamic interplay between genomic and epigenomic factors across long and short evolutionary timescales.

A novel method for identification and quantification of consistently differentially methylated regions

Advances in biotechnology have resulted in large-scale studies of DNA methylation. A differentially methylated region (DMR) is a genomic region with multiple adjacent CpG sites that exhibit different methylation statuses among multiple samples. Many so-called “supervised” methods have been established to identify DMRs between two or more comparison groups. Methods for the identification of DMRs without reference to phenotypic information are, however, less well studied. An alternative “unsupervised” approach was proposed, in which DMRs in studied samples were identified with consideration of nature dependence structure of methylation measurements between neighboring probes from tiling arrays. Through simulation study, we investigated effects of dependencies between neighboring probes on determining DMRs where a lot of spurious signals would be produced if the methylation data were analyzed independently of the probe. In contrast, our newly proposed method could successfully correct for this effect with a well-controlled false positive rate and a comparable sensitivity. By applying to two real datasets, we demonstrated that our method could provide a global picture of methylation variation in studied samples. R source codes to implement the proposed method were freely available at http://www.csjfann.ibms.sinica.edu.tw/eag/programlist/ICDMR/ICDMR.html.

Insights from a chimpanzee adipose stromal cell population: opportunities for adult stem cells to expand primate functional genomics

Comparisons between humans and chimpanzees are essential for understanding traits unique to each species. However, linking important phenotypic differences to underlying molecular changes is often challenging. The ability to generate, differentiate, and profile adult stem cells provides a powerful but underutilized opportunity to investigate the molecular basis for trait differences between species within specific cell types and in a controlled environment. Here, we characterize adipose stromal cells (ASCs) from Clint, the chimpanzee whose genome was first sequenced. Using imaging and RNA-Seq, we compare the chimpanzee ASCs with three comparable human cell lines. Consistent with previous studies on ASCs in humans, the chimpanzee cells have fibroblast-like morphology and express genes encoding components of the extracellular matrix at high levels. Differentially expressed genes are enriched for distinct functional classes between species: immunity and protein processing are higher in chimpanzees, whereas cell cycle and DNA processing are higher in humans. Although hesitant to draw definitive conclusions from these data given the limited sample size, we wish to stress the opportunities that adult stem cells offer for studying primate evolution. In particular, adult stem cells provide a powerful means to investigate the profound disease susceptibilities unique to humans and a promising tool for conservation efforts with nonhuman primates. By allowing for experimental perturbations in relevant cell types, adult stem cells promise to complement classic comparative primate genomics based on in vivo sampling.