Heterogeneous pattern of DNA methylation in developmentally important genes correlates with its chromatin conformation
© The Author(s) 2017
Received: 8 October 2016
Accepted: 13 December 2016
Published: 11 January 2017
DNA methylation is a major epigenetic modification, playing a crucial role in the development and differentiation of higher organisms. DNA methylation is also known to regulate transcription by gene repression. Various developmental genes such as c-mos, HoxB5, Sox11, and Sry show tissue-specific gene expression that was shown to be regulated by promoter DNA methylation. The aim of the present study is to investigate the establishment of chromatin marks (active or repressive) in relation to heterogeneous methylation in the promoter regions of these developmentally important genes.
Chromatin-immunoprecipitation (ChIP) assays were performed to immuno-precipitate chromatin by antibodies against both active (H3K4me3) and repressive (H3K9me3) chromatin regions. The analysis of ChIP results showed that both the percentage input and fold enrichment of activated chromatin was higher in tissues expressing the respective genes as compared to the tissues not expressing the same set of genes. This was true for all the genes selected for the study (c-mos, HoxB5, Sox11, and Sry). These findings illustrate that inconsistent DNA methylation patterns (sporadic, mosaic and heterogeneous) may also influence gene regulation, thereby resulting in the modulation of chromatin conformation.
These findings illustrate that various patterns of DNA methylation (asynchronous, mosaic and heterogeneous) correlates with chromatin modification, resulting in the gene regulation.
KeywordsDNA methylation Chromatin conformation Chromatin-immunoprecipitation (ChIP) Sox11 C-mos HoxB5 Tissue-specific gene expression
DNA methylation is a major component of epigenetic mechanism which is stably inherited and playing important roles during transcriptional regulation [1, 2]. This event implicates the cytosine residues to be directly modified which is immediately followed by guanine residues (CpGs). These CpG dinucleotides are extremely underrepresented in mammalian genomes and are usually present in small stretches known as CpG islands. These CpG islands are found in approximately 70% of annotated gene promoters and are normally unmethylated in normal somatic tissues. Methylation of these CpG targets leads to loss of gene expression during embryonic development . DNA methylation patterns are well regulated, non-random and tissue-specific which is consistent with its functional importance . DNA methylation is also involved in other fundamental processes like X chromosome inactivation, genomic imprinting, suppression of retrotransposon elements, and is essential for normal development . The principal importance of DNA methylation has been studied in myriad of biological context and the mechanism involves the influence of overall chromatin structure. DNA methylation causes the attraction of repressive complexes and the recruitment of methyl CpG binding proteins (including MeCP2) thereby making the chromatin structure inaccessible for binding of transcription factors leading to the formation of closed chromatin structure . Lucifero et al. [7, 8] reported differences in DNA methylation between maternal and paternal alleles of many imprinted genes (Snrpn, Mest and Peg3). Acquisition of methylation was asynchoronous and heterogeneous at these different genes however the varied patterns of methylation on this subset of imprinted genes have not been analyzed for their role in regulating gene expression. Snrpn methylation data showed unique mosaic pattern of specific cytosine methylation during post natal oocyte development. E-cadherin, a tumor suppressor gene, showed de novo methylation of upstream and downstream regions in neoplastic tissues as predominantly methylated/unmethylated CpG islands .
Alternatively DNA methylation plays a pivotal role in embryonic development and differentiation. Developmental genes belonging to different gene families such as Hox and Sox family elicits distinct developmental programs to regulate developmental plasticity during evolution. Various developmentally important genes have been selected for the present study (HoxB5, c-mos, Sox11 and Sry) showing varied temporal and spatial tissue-specific gene expression profile and this variation has been shown to be influenced by DNA methylation patterns directly or indirectly during ontogenesis [11–14].
Therefore the present study focuses to understand the correlation of inconsistency and heterogeneity in DNA methylation patterns during developmental time with altered chromatin conformation and thereby gene regulation. However no reports were found to support this evidence enlightening its relevance in regulation of gene expression and chromatin modifications. Distribution of active (open) or inactive (closed) chromatin marks are studied using ChIP (chromatin-immunoprecipitation) assay for gaining insight towards a gross correlation between the heterogeneous methylation pattern and the distribution of the two chromatin states leading to transcriptional regulation.
Genomic DNA extraction
The Parkes strain of mouse was used in this study. Work was approved by the Institutional Ethics Committee (IEC, Ref No. 176/R&C/13-14) of Motilal Nehru National Institute of Technology Allahabad and the Project Ref. No. is BT/PR13317/GBD/27/252/2009. Mice were randomly bred and maintained under laboratory conditions. Animals were killed by cervical dislocation and different somatic and germinal tissues of both fetal and adult were excised for DNA isolation. Mouse genomic DNA was isolated in three independent sets from adult tissues while 6–7 fetuses were pooled in fetal/neonatal stage especially whole mesonephros gonadal complex from 12.5dpc embryos. DNA was isolated with the help of standard protocol using Proteinase K digestion (50 μg/ml) at 37 °C for 12–14 h. Phenol: chloroform: isoamylalcohol/chloroform: isoamyl alcohol extraction was done at 25 °C. Finally, DNA precipitation was done by adding 1/30th volume of 3 M sodium acetate (pH 5.0) along with two volumes of absolute chilled ethanol.
Sodium bisulfite treatment
Sodium bisulfite treatment converts cytosine residues to uracil except the methylated cytosine residues which imparts unique distinguishing pattern between methylated and non-methylated cytosines in DNA. The bisulfite conversion was done using the BisulFlash™ DNA Modification Kit (Cat. P-1026, USA) according to manufacturer’s instructions. 400–500 ng of input DNA was used as starting template for bisulfite conversion reaction. The bisulfite converted DNA approx. 1.0–2.0 μl was used for each PCR reaction and was amplified in duplicates which was then pooled together during purification. The primers (forward and reverse each with 20 µM) for c-mos gene were as follows: Region (1970–2189 bp): FP 5′-TGTTTTATGTGATTGTTTTATTTG-3′, RP 5′-CACAAAAACACCATAATAAATAAC-3′. Region (2242 bp-2537 bp): FP 5′-ATTTATTATGGTGTTTTTGTGGTTA-3′, RP 5′-ATTCACTAACTTCAAATCCAAATAC-3′. PCR conditions include denaturation @ 94 °C, annealing @ 58 °C and 54 °C for region I and II respectively, extension @ 72 °C repeated for 30 cycles.
Cloning and sequencing of positive clones
Amplified PCR products were gel excised and purified from low melting agarose using a DNA purification kit (nucleopore) as per instructions manual. The ligation of purified PCR product was done in T-vector using an InsTAclone™ PCR cloning kit (Thermo Scientific USA) and the ligated product was finally transformed and cloned in E.coli DH5α. Positive transformed clones were selected on the basis of blue-white screening and plasmid isolation was done using an alkaline lysis method (Thermo Scientific Plasmid Miniprep Kit). Sequencing of 8-10 independent clones were done by automated DNA sequencer (ABI 3130 genetic analyzer). Plasmid DNA 50–100 ng was used for sequencing with BigDye® Terminator v3.1 Cycle Sequencing Kit as per manufacturer’s instructions.
RNA isolation and quantitative real time PCR
Total RNA was isolated from mouse tissues using TRI reagent (sigma) according to the manufacturer’s protocol. Concentration of RNA was determined spectrophotometrically at OD 260 nm both before and after DNAse treatment. RNA was reverse transcribed in equal amount using oligo-dT(18) primers and 200 U of MMLV Reverse transcriptase (NEB). cDNA 0.5 µl was used in PCR amplification using SsoFast™ EvaGreen Supermix with Low ROX (2X) along with 5 µM each of forward and reverse primers in real time step one plus PCR machine (ABI). PCR reactions were carried out in triplicates for both GAPDH and c-mos. Expression of mRNA level was estimated and normalized relative to the mRNA level of GAPDH. The average of the cycle threshold (Ct) values and standard deviation were determined. Fold change in gene expression was calculated according to 2− ∆∆Ct method . The primer sequences of GAPDH are FP: 5′GGAGCCAAACGGGTCATCATCTC3′ and RP-5′GAGGGGCCATCCACAGTCTTCT 3′; c-mos gene FP 5′-TACGCCACGACAACATAGTTCG-3′ RP 5′-CTTGCTCACTGATCAAAATGTTGG-3′.
ChIP assay was performed according to the instructions manual (Diagenode ChIP kit Cat. No. kch-orgHIS-012). Chromatin was isolated from different somatic (brain, spleen and kidney) and germinal tissues (testis) of adult, fetal and neonatal stages of mouse. The excised tissues were homogenized and subjected to collagenase treatment (50–200 U/ml) followed by incubation for 2–3 h at 37 °C. Single cell suspension was made by pipetting during the incubation time and cell counting was performed using haemocytometer. The minimum number of cells required to perform ChIP experiments is 1 × 106 cells. Cell cross-linking was done by adding 37% formaldehyde (w/v, final concentration 1%) kept for 10 min at 25 °C on a rotating wheel followed by quenching with 1.25 M glycine (final concentration 125 mM) for 5 min at 25 °C, centrifuged at 4 °C for 5–8 min. Supernatant was discarded and the cell pellet was resuspended in lysis buffer (containing protease inhibitors). The cell suspension was subjected to sonication using a sonicator (SKN-IIDN) at the rate of 3 s ON/1 s OFF for 3–4 cycles for obtaining the desired chromatin range from 200–800 bp. The sheared chromatin was then processed for pre-clearing by adding an IP-incubation mix and pre-blocked beads. Antibodies specific for capturing the desired protein and interacting DNA were used (H3K4me3, Diagenode MAb-152-050 and H3K9me3, Diagenode, MAb-146-050, concentration 1 µg/µl). Negative control IgG antibody (Diagenode, C15400001 (C15200001) was used which binds with non-specific target and the associated DNA fragments were immuno-precipitated. The addition of specific antibodies was followed by incubation on a rotating wheel at 4 °C for overnight. Bead washing with wash buffer—1, 2 and 3 removes non-associated DNA fragments and Protein/DNA complexes were found to get eluted from pre-blocked beads by the addition of elution buffer. The eluted complex was reversibly cross-linked and purified using phenol: chloroform: iso-amyl alcohol/chloroform: iso-amyl alcohol. DNA fragments were precipitated by adding DNA precipitant, DNA co-precipitant and absolute chilled ethanol. The DNA pellet was resuspended in 30 µl of milliQ water and the relative amount of specifically immunoprecipitated DNA was analyzed through PCR amplification using quantitative real-time PCR (ABI step one plus) with 1.0 µl of DNA, SsoFast™ EvaGreen Supermix (2X) with Low ROX (Biorad) and gene specific primers forward and reverse 5 µM each. Control primers GAPDH (c17021045, Diagenode used as positive control against activated chromatin regions) and TSH2β (c17021042, Diagenode used as positive control against repressed chromatin regions) were used. The percentage input and fold enrichment was calculated which represents the enrichment of certain histone modifications on specific region using the ChIP reactions performed in triplicate.
Shows the primer sequence of different genes used for ChIP-qPCR reactions
We have selected three more developmentally important genes whose methylation pattern has been already studied in different tissues include kidney, brain, spleen, testis and mesonephros gonadal cells (MGCs) [11–13]. In HoxB5 methylation analysis was performed in fetal and adult stages of mouse kidney and spleen tissues whereas the same analysis was carried out in fetal and adult stages of mouse brain and kidney tissues in Sox11. Promoter methylation pattern of Sry was determined in mouse 12.5 dpc mesonephros gonadal complex (MGCs) and in adult testis. Comparing the methylation status of HoxB5 and Sox11, higher methylation was found in adult as compared to fetus. In Sry, less methylation was present in fetal gonads as compared to fetal liver. The work done in the above cited references showed a direct inverse correlation between tissue-specific developmental methylation pattern and gene expression.
Methylation and expression pattern of c-mos gene
ChIP (Chromatin-Immunoprecipitation) results
The results of Chromatin-Immunoprecipitation demonstrate the occupancy of the fractionated DNA fragments precipitated with a particular antibody against activated (H3K4me3) or repressed (H3K9me3) chromatin domains for specific gene in adult, fetal and neonatal stages of various somatic and germinal tissues of mice. The data were represented as percentage input and in terms of fold enrichment (FE).
The chromatin interaction results of c-mos gene in adult testis illustrates that the percentage input of activated chromatin (H3K4me3) was higher than adult kidney (Fig. 1f, h). In a similar way fold enrichment was also higher (19 fold) in activated chromatin regions of adult testis as compared to adult kidney (Fig. 1g, i).
The heterogeneous methylation pattern present in the somatic tissues (adult kidney shows the average percentage of all CpGs around 55) and almost complete absence of methylation in the germinal tissues (testis and ovary) of c-mos gene correlates inversely with its tissue specific expression. c-mos displays higher expression in adult testis and ovary and almost no expression in somatic tissues. It further strengthens the crucial role of c-mos gene in meiotic maturation of germinal tissues [14, 15] and in spermatogenesis  and oogenesis . These finding correlates with the chromatin–immunoprecipitation results that displays the fractionation of chromatin fragments having the higher percentage input of activated chromatin regions (H3K4me3) in the expressing tissue (germinal tissues) in contrast to repressed chromatin regions (H3K9me3) in non-expressing tissues (somatic). Similarly HoxB5, a member of Hox family genes playing important roles in antero-posterior axial body patterning and the RNA transcripts were found in a wide range of fetal tissues such as lung, gut, kidney, liver and spleen [18–20]. Examination of the methylation pattern of HoxB5 shows higher methylation density in adult as compared to fetal stages of mouse tissues (Kidney and spleen) . ChIP data also suggests higher input fraction of activated chromatin during fetal/neonatal stages of tissues in comparison to their adult counter parts.
Another developmental gene, Sox11, whose expression was prominent in developing central nervous system of mouse embryos, suggests its important role in neuronal growth, maturation and survival [21–23]. Promoter methylation mapping of Sox11 reveals higher methylation in both adult brain and kidney (average percentage of all CpGs around 48 and 68% in brain and kidney respectively) while the fetal and neonatal stages of the same tissues showed lower methylation density (14 and 31% in fetal brain and kidney respectively) . The ChIP results of Sox11 in the present study are highly consistent with its tissue specific gene expression which correlates with changes in chromatin conformation.
A study by Nishino et al.  demonstrate that the mechanistic role of Sry gene regulation occurs by DNA methylation during mouse gonadal development and the Sry transcripts were expressed only in pre Sertoli cells of developing male gonad during 10.5–12.5dpc (average percentage of all CpGs around 62% and <3% in 11.5dpc liver and gonads respectively) [24, 25]. Hypermethylated status of the two regions analyzed was maintained in the tissues that did not express Sry however the methylation pattern was not purely consistent. Our ChIP results coincides with this finding as higher percentage input and fold enrichment values of activated chromatin fragments in fetal mesonephros gonads are in order with tissue specific gene expression  and inversely correlates with the promoter methylation of Sry gene.
Curradi et al.  analyzed the inhibitory effects of promoter DNA methylation through in vitro experiments and the results showed that a definite number of methylated CpG residues are required to organize a stable, diffusible chromatin structure. Histone deacetylation significantly imparts to transcriptional repression only when the number of modified cytosine residues is adequate to suppress the gene over a long range. These data suggest the importance of methylated cytosine residues for its ability to create a repressive effect and this repression can propagate in both directions for several hundred base pairs of the altered DNA only when enough number of methylated CpGs are present .
The potential role of DNA methylation is conserved from yeast to mammals. Some recent reports exhibited tissue and life-stage specific gene expression and its regulation by DNA methylation in species other than mammals suggesting functional similarity between transcriptional regulation [27, 28]. However, other studies also show impact of heterogeneous DNA methylation on gene expression including Oct-4 whose expression is confined to the germ line and pluripotent stem cells, which is critical for normal mouse development. The epigenetic regulation of Oct-4 show heterogeneous DNA methylation in its regulatory region among the population of adult somatic cells . Similarly the heterogeneity in DNA methylation patterns across the distinct promoter and 5′UTR regions of differentiating stem cells and during reprogramming events suggest new insight into the methylation dynamics pertaining to embryonic stem cells [30, 31]. It has also been reported that the methylation patterns present in several imprinted genes and repetitive elements of embryonic germ (EG) cells may exhibit DNA methylation changes occurring in the germ cell lineage during the differentiation process . Induction of mosaic DNA methylation patterns also leads to tumorigenesis by participating in important cancer related pathways and thus causing loss of co-expression connectivity in colorectal cancer . The presence of epigenetic heterogeneity of developmentally important genes has important implications for assisted reproduction outcomes . Another significant role was observed in the epigenetic regulation of stress reactivity in humans  and the formation and maintenance of memory in plasticity genes .
These findings elucidate that the tissue and stage specific gene regulation of several developmentally important genes of mouse is mediated through promoter/regulatory regions DNA methylation. The ChIP results are also consistently indicating a direct correlation between the percentage of input chromatin fraction and the level of gene expression.
Transcriptional regulation mediated by the establishment of tissue specific methylation patterns in the promoter or regulatory regions of genes examined in this study might speculate to be a crucial mechanism which may be different from non-developmental genes. This might reflect the specialized functions of various developmental genes during the process of embryonic development. A significant level of heterogeneous/asynchronous methylation exists as the development progresses from embryo to adult and in most of the cases this heterogeneity is associated with transcriptional gene regulation and chromatin modifications. Further, higher density of methylation present in the promoter regions might render the DNA inaccessible (either fully or partially) for binding of transcription factors, leading to altered chromatin assembly and gene expression.
days post coitum
days post partum
methyl cytosine binding protein 2
PS carried out the experimental work and drafted the manuscript. MS and KS supervised the study. All authors read and approved the final manuscript.
This work was supported by Department of Biotechnology (DBT), Government of India. We thank Cytogenetic Lab, Department of Zoology, Banaras Hindu University (BHU) for providing animal house facility.
The authors declare that they have no competing interests.
Availability of data and materials
All data supporting the conclusion of this article are contained within the manuscript.
Ethics approval and consent to participate
All animal experiments were approved by the Institutional Ethics Committee (IEC, Ref No. 176/R&C/13-14) of Motilal Nehru National Institute of Technology Allahabad.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Kulis M, Esteller M. DNA methylation and cancer. Adv Genet. 2010;70:27–56.PubMedGoogle Scholar
- Suzuki MM, Bird A. DNA methylation landscapes: provocative insights from epigenomics. Nature Rev Genet. 2008;9(6):465–76.View ArticlePubMedGoogle Scholar
- Deaton AM, Bird A. CpG islands and the regulation of transcription. Genes Dev. 2011;25:1010–22.View ArticlePubMedPubMed CentralGoogle Scholar
- Illingworth R, Kerr A, Desousa D, Jørgensen H, Ellis P, Stalker J, Jackson D, Clee C, Plumb R, Rogers J, Humphray S, Cox T, Langford C, Bird A. A novel CpG island set identifies tissue-specific methylation at developmental gene loci. PLoS Biol. 2008;6(1):37–51.View ArticleGoogle Scholar
- Chen ZX, Riggs AD. DNA methylation and demethylation in mammals. J Biol Chem. 2011;286(21):18347–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Singal R, Ginder GD. DNA methylation. Blood. 1999;93(12):4059–70.PubMedGoogle Scholar
- Lucifero D, Mertineit C, Clarke HJ, Bestor TH, Trasler JM. Methylation dynamics of imprinted genes in mouse germ cells. Genomics. 2002;79(4):530–8.View ArticlePubMedGoogle Scholar
- Lucifero D, Mann MR, Bartolomei MS, Trasler JM. Gene-specific timing and epigenetic memory in oocyte imprinting. Hum Mol Genet. 2004;13(8):839–49.View ArticlePubMedGoogle Scholar
- Graff JR, Herman JG, Lapidus RG, Chopra H, Xu R, Jarrard DF, Isaacs WB, Pitha PM, Davidson NE, Baylin SB. E-cadherin expression is silenced by DNA hypermethylation in human breast and prostate carcinomas. Cancer Res. 1995;55(22):5195–9.PubMedGoogle Scholar
- Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25:402–8.View ArticlePubMedGoogle Scholar
- Sachan M, Raman R. Developmental methylation of the regulatory region of HoxB5 gene in mouse correlates with its tissue-specific expression. Gene. 2006;380(2):151–8.View ArticlePubMedGoogle Scholar
- Pamnani M, Sinha P, Nara S, Sachan M. Study of promoter DNA methylation of Sox11 and its correlation with tissue-specific expression in the laboratory mouse. Gene. 2014;552(1):133–9.View ArticlePubMedGoogle Scholar
- Nishino K, Hattori N, Tanaka S, Shiota K. DNA methylation-mediated control of sry gene expression in mouse gonadal development. J Biol Chem. 2004;279(21):22306–13.View ArticlePubMedGoogle Scholar
- Kanduri C, Raman R. Tissue-specific characterization of DNA methylation in the gonad-specific proto-oncogene, c-mos, in the male laboratory mouse. Int J Dev Biol. 1999;43:91–4.PubMedGoogle Scholar
- Mutter GL, Grills GS, Wolgemuth DJ. Evidence for the involvement of the proto-oncogene c-mos in mammalian meiotic maturation and possibly very early embryogenesis. EMBO J. 1988;7(3):683–9.PubMedPubMed CentralGoogle Scholar
- Cao SF, Li D, Yuan Q, Guan X, Xu C. Spatial and temporal expression of c-mos in mouse testis during postnatal development. Asian J Androl. 2008;10(2):277–85.View ArticlePubMedGoogle Scholar
- Jałocha I, Gabryś MS, Bal J. The crucial role of the proto-oncogene c-mos in regulation of oocyte maturation. Postepy Hig Med Dosw. 2010;64:636–41.Google Scholar
- Hershko AY, Kafri T, Fainsod A, Razin A. Methylation of HoxA5 and HoxB5 and its relevance to expression during mouse development. Gene. 2003;302(1–2):65–72.View ArticlePubMedGoogle Scholar
- Krumlauf RO, Holland PW, McVEY JH, Hogan BL. Developmental and spatial patterns of expression of the mouse homeobox gene, Hox 2.1. Development. 1987;99(4):603–17.PubMedGoogle Scholar
- Fu M, Lui VC, Sham MH, Cheung AN, Tam PK. HOXB5 expression is spatially and temporarily regulated in human embryonic gut during neural crest cell colonization and differentiation of enteric neuroblasts. Dev Dyn. 2003;228(1):1–10.View ArticlePubMedGoogle Scholar
- Hargrave M, Wright E, Kun J, Emery J, Cooper L, Koopman P. Expression of the Sox11 gene in mouse embryos suggests roles in neuronal maturation and epithelio-mesenchymal induction. Dev Dyn. 1997;210(2):79–86.View ArticlePubMedGoogle Scholar
- Jankowski MP, Cornuet PK, McIlwrath S, Koerber HR, Albers KM. SRY-box containing gene 11 (Sox11) transcription factor is required for neuron survival and neurite growth. Neuroscience. 2006;143(2):501–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang Z, Reynolds A, Kirry A, Nienhaus C, Blackmore MG. Overexpression of Sox11 promotes corticospinal tract regeneration after spinal injury while interfering with functional recovery. J Neurosci. 2015;35(7):3139–45.View ArticlePubMedPubMed CentralGoogle Scholar
- Rossi P, Dolci S, Albanesi C, Grimaldi P, Geremia R. Direct evidence that the mouse sex-determining gene Sry is expressed in the somatic cells of male fetal gonads and in the germ cell line in the adult testis. Mol Reprod Dev. 1993;34(4):369–73.View ArticlePubMedGoogle Scholar
- Hacker A, Capel B, Goodfellow P, Lovell-Badge R. Expression of Sry, the mouse sex determining gene. Development. 1995;121:1603–14.PubMedGoogle Scholar
- Curradi M, Izzo A, Badaracco G, Landsberger N. Molecular mechanisms of gene silencing mediated by DNA methylation. Mol Cell Biol. 2002;22(9):3157–73.View ArticlePubMedPubMed CentralGoogle Scholar
- Li S, Zhu Y, Zhi L, Han X, Shen J, Liu Y, Yao J, Yang X. DNA methylation variation trends during the embryonic development of chicken. PLoS ONE. 2016;11(7):e0159230.View ArticlePubMedPubMed CentralGoogle Scholar
- Mayasich SA, Bemis LT, Clarke BL. DNA methylation in a sea lamprey vasotocin receptor gene promoter correlates with tissue-and life-stage-specific mRNA expression. Comp Biochem Physiol. 2016;202:56–66.View ArticleGoogle Scholar
- Marikawa Y, Fujita TC, Alarcón VB. Heterogeneous DNA methylation status of the regulatory element of the mouse Oct4 gene in adult somatic cell population. Cloning Stem Cells. 2005;7(1):8–16.View ArticlePubMedGoogle Scholar
- Shao X, Zhang C, Sun MA, Lu X, Xie H. Deciphering the heterogeneity in DNA methylation patterns during stem cell differentiation and reprogramming. BMC Genom. 2014;15:978–88.View ArticleGoogle Scholar
- Singer ZS, Yong J, Tischler J, Hackett JA, Altinok A, Surani MA, Cai L, Elowitz MB. Dynamic heterogeneity and DNA methylation in embryonic stem cells. Mol Cell. 2014;55(2):319–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Shovlin TC, Durcova-Hills G, Surani A, McLaren A. Heterogeneity in imprinted methylation patterns of pluripotent embryonic germ cells derived from pre-migratory mouse germ cells. Dev Biol. 2008;313(2):674–81.View ArticlePubMedGoogle Scholar
- Wang Q, Jia P, Cheng F, Zhao Z. Heterogeneous DNA methylation contributes to tumorigenesis through inducing the loss of coexpression connectivity in colorectal cancer. Genes Chromosomes Cancer. 2015;54(2):110–21.View ArticlePubMedGoogle Scholar
- Kuhtz J, Schneider E, Hajj NE, Zimmermann L, Fust O, Linek B, Seufert R, Hahn T, Schorsch M, Haaf T. Epigenetic heterogeneity of developmentally important genes in human sperm: implications for assisted reproduction outcome. Epigenetics. 2014;9(12):1648–58.View ArticlePubMedGoogle Scholar
- Houtepen LC, Vinkers CH, Carrillo-Roa T, Hiemstra M, van Lier PA, Meeus W, Branje S, Heim CM, Nemeroff CB, Mill J, Schalkwyk LC, Creyghton MP, Kahn RS, Joëls M, Binder EB, Boks MP. Genome-wide DNA methylation levels and altered cortisol stress reactivity following childhood trauma in humans. Nat Commun. 2016;7:10967.View ArticlePubMedPubMed CentralGoogle Scholar
- Halder R, Hennion M, Vidal RO, Shomroni O, Rahman RU, Rajput A, Centeno TP, Bebber FV, Capece V, Vizcaino JCG, Schuetz AL, Burkhardt S, Benito E, Sala MN, Javan SB, Haass C, Schmid B, Fischer A, Bonn S. DNA methylation changes in plasticity genes accompany the formation and maintenance of memory. Nature Neurosci. 2016;19:102–17.View ArticlePubMedGoogle Scholar