DNA methylation of the LIN28 pseudogene family
© Davis et al.; licensee BioMed Central. 2015
Received: 15 August 2014
Accepted: 25 March 2015
Published: 11 April 2015
DNA methylation directs the epigenetic silencing of selected regions of DNA, including the regulation of pseudogenes, and is widespread throughout the genome. Pseudogenes are decayed copies of duplicated genes that have spread throughout the genome by transposition. Pseudogenes are transcriptionally silenced by DNA methylation, but little is known about how pseudogenes are targeted for methylation or how methylation levels are maintained in different tissues.
We employed bisulfite next generation sequencing to examine the methylation status of the LIN28 gene and four processed pseudogenes derived from LIN28. The objective was to determine whether LIN28 pseudogenes maintain the same pattern of methylation as the parental gene or acquire a methylation pattern independent of the gene of origin. In this study, we determined that the methylation status of LIN28 pseudogenes does not resemble the pattern evident for the LIN28 gene, but rather these pseudogenes appear to acquire methylation patterns independent of the parental gene. Furthermore, we observed that methylation levels of the examined pseudogenes correlate to the location of insertion within the genome. LIN28 pseudogenes inserted into gene bodies were highly methylated in all tissues examined. In contrast, pseudogenes inserted into genomic regions that are not proximal to genes were differentially methylated in various tissue types.
Our analysis suggests that Lin28 pseudogenes do not aquire patterns of tissue-specific methylation as for the parental gene, but rather are methylated in patterns specific to the local genomic environment into which they were inserted.
DNA methylation controls diverse aspects of genome regulation and transcriptional activity. Methylation of mammalian DNA involves the addition of a methyl group to the 5’-carbon of the cytosine in a cytosine-guanine (CpG) dinucleotide. This system of methylation likely evolved from a genomic defense system responsible for preventing the spread of parasitic genetic elements. DNA methylation has since evolved to play an active role in maintaining genetic structure and genome regulation . Methylation is involved in X-chromosome inactivation [2,3], silencing of transposable elements [4-7], tissue-specific gene expression [8-11], and gene imprinting [12-15]. DNA methylation is widespread throughout the genome, and the maintenance of methylation patterns is highly regulated and tissue-specific [16,17]. DNA methyltransferases (DNMTs) are responsible for de novo and maintenance methylation of the genome [18-20]. The proper establishment and maintenance of methylation patterns is critical for early development and the absence of DNA methylation results in embryonic lethality [21-23].
DNA methylation also regulates pseudogenes within the genome [16,24]. Pseudogenes are decayed copies of active genes that have arisen from either a duplication event, in which the entire gene or portion of a gene is duplicated (non-processed pseudogenes), or from the retrotransposition of an RNA transcript into the genome (processed pseudogenes). An analysis of the human genome estimates that as many as 19,000 pseudogenes are evenly distributed throughout the genome, and approximately 70% of these are processed pseudogenes [20,25]. Ten percent of genes within the human genome have at least one corresponding pseudogene [20,26], and pseudogenes primarily arise from parental genes that are transcriptionally active within the germ line.
The same regulatory network that inhibits transposable element movement likely induces DNA methylation on pseudogenes. Methylation of pseudogenes is elevated in embryos, likely as a mechanism for preventing the spread of transposable elements during embryogenesis [16,27]. In plants, the inactivation of methyltransferases resulted in the widespread activation of transposable elements and pseudogenes , demonstrating that DNA methylation is sufficient to prevent the activation of pseudogenes. In humans, pseudogenes are highly methylated, presumably to prevent transcription and further transposition .
Characterizing methylation patterns of pseudogenes is critical, as pseudogenes with high sequence identity to parental genes can lead to misinterpretation of results in methylation studies [28,29]. The characterization of pseudogene methylation signatures also reveals how DNA segments are regulated by methylation networks once inserted into the genome. In order to better understand how methylation patterns are established and maintained on pseudogenes, we examined the methylation status of four pseudogenes derived from the translational enhancer LIN28. This gene is important in early embryo development and can also act as a reprogramming factor in the production of induced pluripotent stem cells [30,31].
LIN28 has given rise to at least ten processed pseudogenes within the bovine genome that vary in length between 100 to 4000 bp (Btau_4.6.1, released Nov 2, 2011). The protein-coding region of LIN28 contains a high concentration of CpG sites, making the gene a potential target for DNA methylation once inserted elsewhere in the genome as a pseudogene. By measuring the methylation levels of selected pseudogenes and the LIN28 gene, we sought to determine whether the same regulatory mechanism that directs and maintains methylation of the LIN28 gene also controls the methylation status of LIN28 pseudogenes. Additionally we examined the expression of genes near the insertion site to determine whether pseudogene methylation is involved in transcriptional control of adjacent genes. This study is the first to characterize the methylation status of LIN28 and its associated pseudogenes, and is the first such research to characterize the methylation status of a pseudogene family within the bovine genome.
Methylation status of the LIN28 family
Our results indicate that LIN28 and all examined pseudogenes were methylated in all bovine cell and tissue types examined. Due to the sequence similarity between LIN28 and its derived pseudogenes we were only able to obtain coverage of seven CpG dinucleotides within the LIN28 gene. This region was selected because of primer-specific amplification adjacent to the region examined in the pseudogenes. These dinucleotides are located in the third exon of LIN28. All of the CpG sites were methylated in each tissue to varying degrees (Figure 3; Additional file 1: Figures S5-S9). In addition, all pseudogenes examined were methylated in all tissue samples. However, each pseudogene demonstrated a distinct methylation pattern that differed sharply from the LIN28 gene, and that correlates with the pseudogene insertion location.
The pseudogenes LIN28P-Ch:3 and LIN28P-Ch:28 were inserted into a location isolated from any gene, and both of these pseudogenes share a similar overall pattern of methylation (Figure 3B, E; Additional file 1: Figures S6, S9). We analyzed methylation at nine CpG dinucleotides within LIN28P-Ch:3 and eight CpG dinucleotides in LIN28P-Ch:28. In both pseudogenes, oocytes had the lowest levels of overall methylation and blastocysts had the highest level of methylation. The methylation patterns of both pseudogenes apparently deviated from patterns observed in the LIN28 gene. Within the LIN28 gene, oocytes had the highest frequency of methylation. Alternatively, the methylation of the parental gene was lowest in blastocyst embryos. In LIN28P-Ch:28 we observed three overlapping CpG dinucleotides with the LIN28 gene. Although limited, within these three overlapping sites, the inverse pattern of methylation between the LIN28 gene and pseudogenes LIN28P-Ch:3 and LIN28P-Ch:28 is maintained.
In contrast to the tissue-specific methylation patterns, pseudogenes LIN28P-Ch:7 and LIN28P-Ch:26 both are highly methylated (Figure 3C-D; Additional file 1: Figures S7-S8). Both LIN28P-Ch:7 and LIN28P-Ch:26 were inserted into gene introns. LIN28P-Ch:7 is inserted into the fourth intron of the gene MAN2A1. A majority of the 27 CpG sites inspected within this pseudogene were methylated in all tissue samples (Figure 3C; Additional file 1: Figure S7). The LIN28P-Ch:26 pseudogene lies within the first intron of the gene ACADSB and is 20 kbp from the gene IKZF5. Similar to LIN28P-Ch:7, all 26 CpG sites examined for the LIN28P-Ch:26 pseudogene were highly methylated at a high frequency (Figure 3D; Additional file 1: Figure S8). The high levels of methylation observed in both of these pseudogenes is in sharp contrast to the moderate levels of methylation of the LIN28 gene and pseudogenes not inserted into gene bodies, none of which are as highly methylated for all tissue samples.
These apparent differences in methylation correlate to the location of insertion of the pseudogene. LIN28 pseudogenes inserted into gene bodies were highly methylated. Alternatively, LIN28 pseudogenes with an insertion location distant from a gene varied in the pattern of methylation in a tissue-specific manner, although these tissue-specific patterns were inverse with respect to the LIN28 gene (Figure 3; Additional file 1: Figures S5-S9).
We observed two cytosine nucleotides that were not contained within CpG dinucleotides that were methylated to some degree in all tissues examined. Both cytosine nucleotides were located within LIN28P-Ch:7. Both non-CpG methylation sites were located within close proximity to CpG sites. One non-CpG methylated cytosine was located between CpG sites three and four (CGCCG) on LIN28P-Ch:7. The percentage of sequences that were methylated at this location in the tissues examined were: 9% in testes, 13% in liver, 20% in brain, 46% in fibroblast, 46% in IVF blastocyst, and 4% in oocyte. Both CpG sites flanking the specified cytosine are highly methylated in all samples (Figure 3C; Additional file 1: Figure S7). The second non-CpG methylated cytosine observed is upstream from CpG site 13 in LIN28P-Ch:7 (CCG). The percentage of methylated sequences at this cytosine were observed to be: 8% in testes, 10% in liver, 16% in brain, 38% in fibroblast, 40% in IVF blastocyst, and 100% in oocyte. Methylation of CpG site 13 on LIN28P-Ch:7 is high in all tissue samples (Figure 3C; Additional file 1: Figure S7). Additionally, the methylation status of either non-CpG cytosine nucleotides is independent of the other.
Single CpG demethylation
LIN28P-Ch:26 contains a single CpG dinucleotide that is characterized by a low frequency of methylation in a tissue-specific manner. This CpG dinucleotide is under-methylated in a tissue-specific manner at CpG site 16 in IVF blastocyst and oocyte samples (Figure 3D; Additional file 1: Figure S8). Both samples are also under-methylated relative to adjacent CpG sites with a 13% methylation frequency in IVF blastocyst and 5% frequency in oocytes, yet highly methylated for other tissue types.
Expression of genes adjacent to LIN28 pseudogenes
The objective of our study was to determine whether the methylation of pseudogenes followed the same patterns as the functional parental gene by examining LIN28 as a case study. An observation that the LIN28 pseudogenes maintained an identical or highly similar methylation pattern as the LIN28 gene would indicate that regulation of LIN28 methylation is intrinsic to the gene sequence and that methylation of pseudogenes with high sequence identity is likely controlled by the same mechanism(s) that maintain methylation of the LIN28 gene. Alternatively, an observation that pseudogene methylation patterns deviated from the parental gene would indicate that pseudogenes are subject to local regulation of methylation patterns. In this study, we observed that LIN28 pseudogenes do not recapitulate the same methylation status as LIN28, but rather appear to acquire methylation patterns independent of the parental gene. Furthermore, we observed that methylation levels of the examined pseudogenes correlate to the location of insertion. To our knowledge, this study is the first to characterize the methylation signatures of the LIN28 pseudogene family and to identify an effect of genome location on pseudogene methylation. Previous studies have shown that methylation patterns of pseudogenes deviate from those of the parental gene, and our findings are consistent with these observations [16,27,29,32].
Our examination of an entire pseudogene family highlights the diversity of methylation patterns apparent for highly similar sequences. LIN28 pseudogenes inserted into gene bodies were highly methylated in all tissues examined. In contrast, pseudogenes inserted into genomic regions that are not proximal to genes had reduced overall methylation and were differentially methylated in unique tissue types. The measurement of methylation patterns in six distinct tissue types showed that methylation of pseudogenes can be highly variable in different tissues (Figure 3; Additional file 1: Figures S5-S9) and suggested that methylation of CpG sites may be differentially regulated in these tissues. Pseudogenes not associated with genes had less methylation in tissue samples that were highly methylated in the parental gene.
We measured methylation of seven CpG dinucleotides located within the third exon of the LIN28 gene. Methylation of CpG dinucleotides within the gene body is generally associated with transcriptionally active genes [25,26]. This pattern is in contrast to methylation of the 5’ upstream and promoter regions of genes, which are typically associated with transcriptional silencing. Methylation of these seven CpG dinucleotides was dependent upon tissue type. Notably, oocytes had the highest levels of methylation within the LIN28 gene. Oocytes generally maintain low levels of global methylation relative to somatic cells [33,34]. Following fertilization, global methylation levels decline further and are then reestablished during embryonic and somatic cell development [35-38]. The methylation level observed in the LIN28 gene is counter to this pattern. However, methylation levels measured in the LIN28P-Ch:3 and LIN28P-Ch:28 pseudogenes are consistent with changes in methylation levels on a global level. Within both these pseudogenes, oocytes maintained low levels of methylation, whereas somatic cells maintained high levels of methylation. It is possible that the methylation patterns observed for LIN28P-Ch:3 and LIN28P-Ch:28 are maintained by the same mechanism that maintains global levels of methylation, while methylation of the LIN28 gene is maintained by a separate mechanism.
We also observed that methylation of pseudogenes depends on the genomic context into which the pseudogene was inserted. Both pseudogenes inserted into the intron of a gene, LIN28P-Ch:7 and LIN28P-Ch:26, were highly methylated in all tissues. The absence of variability among tissues is revealing, as insertion of a pseudogene into a gene body appears to induce high levels of methylation. On the other hand, both pseudogenes not associated with a gene, LIN28P-Ch:3 and LIN28P-Ch:28, vary in the pattern of methylation within different tissue samples, but tissue methylation levels were similar between the two pseudogenes. The similarity of methylation levels among all tissue types between similar pseudogenes indicates that pseudogenes share a common regulatory mechanism that establishes and maintains the methylation signature. When considering the location of the LIN28 pseudogenes in context of chromosomal structure, we originally postulated that LIN28P-Ch:7 and LIN28P-Ch:26 exist in a region that is likely euchromatin, as both pseudogenes were inserted into gene bodies (Figure 2). Furthermore, expression data obtained for MAN2A1 and IKZF5 (Figure 4) suggest that these regions are transcriptionally active, providing evidence that the region maintains euchromatin structure in the tissues examined. Because LIN28P-Ch:3 and LIN28P-Ch:28 are isolated from any local gene, it is reasonable to expect that the local sequence environment for these genes is heterochromatin in structure. However, our expression analysis did not directly test this hypothesis, since these two pseudogenes do not lie within a gene body.
Other researchers have shown that a limited number of pseudogenes have evolved to play a regulatory function within the genome [39-44]. It is possible that the introduction of a pseudogene within a the promoter of a gene may provide a CpG-rich sequence that can be utilized as a site for gene silencing via DNA methylation, although this observation was not evident in our study. As indicated by the results of our analysis of LIN28, these pseudogenes may undergo tissue-specific methylation unique to the region of pseudogene insertion, which lead us to investigate the relationship between methylated pseudogenes and expression of adjacent genes. However, evidence from PCR analysis of genes associated with the four pseudogenes, including DHDDS, MAN2A1, ACADSB, and IKZF5, suggested that LIN28 pseudogene methylation status was not correlated with decreased gene expression. This observation suggests that in the case of the LIN28 pseudogene family, methylation status of the pseudogene does not play a major role in regulation of expression of the local genes. However, methylation of a single CpG dinucleotide in the LIN28 pseudogene was correlated with gene expression changes for the gene IKZF5. The LIN28P-Ch:26 pseudogene is located within the first intron of ACADSB and 20 Kbp upstream from IKZF5. Within the LIN28P-Ch:26 pseudogene, a single CpG dinucleotide at site 16 is hypomethylated in the oocyte and blastocyst samples, a sharp contrast to the hypermethylated status of all surrounding CpG dinucleotides as well as the same CpG site in the four remaining tissues. Interestingly, IKZF5 was not expressed in either oocytes or blastocyst embryos. Although the expression pattern of IKZF5 was correlated with the tissue-specific methylation patterns observed for the pseudogene, it is doubtful that this single CpG site is involved in the transcriptional regulation of IKZF5.
Within LIN28P-Ch:7, we identified two cytosine nucleotides that are methylated outside of the CpG dinucleotide context. In both cases, the methylation has occurred on a cytosine immediately upstream from a CpG dinucleotide. Non-CpG methylation has been observed to occur with higher frequency in non-dividing cells and gametes, although its function remains unknown . Interestingly, the same cytosine nucleotides are observed within the LIN28P-Ch:26 and are unmethylated in all samples. This observation further suggests that upon integration into the genome, methylation of the LIN28 pseudogene family is location specific and the methylation occurs independent of their sequence.
In conclusion, the four LIN28 pseudogenes examined in this study maintained methylation patterns that deviate from those found in the parental gene as well as from one another according to the genomic location into which the pseudogene was inserted. To our knowledge, this is a unique phenomena identified in our study, and examination of more pseudogene families will be required to determine whether the same observation is consistent for other integrated pseudogenes. Pseudogenes derived from LIN28 have undergone mutations and no longer maintain an exact sequence identity to the LIN28 gene. Sequence identity ranges from 84% to 97% (Figure 1), and it is possible that the differences in methylation pattern are a result of changes in the pseudogene sequence. Additionally, our conclusions are based on methylation of CpG sites within the pseudogenes only, as well as a limited number of CpG sites within the LIN28 gene. Analysis of CpG dinucleotides throughout the LIN28 parent gene as well as CpG dinucleotides flanking the pseudogene point on insertion would further add to our findings. Future work should focus on these CpG sites and would further help determine how DNA methylation is targeted to specific genomic regions. Although methylation status of the LIN28 pseudogenes was not associated with changes gene expression of proximal genes, this observation does not rule out the possibility that CpG-rich pseudogenes could serve as sites for regulation of gene expression by methylation, a hypothesis that may also be addressed by survey of other pseudogene families. New knowledge on the regulation of pseudogenes via DNA methylation could contribute to greater understanding of the maintenance of global and/or regional patterns of methylation. Future work on this topic should focus on characterizing methylation patterns for other pseudogene families to determine whether all pseudogenes are maintained in a similar manner or whether sequence specific patterns can be identified through analysis of pseudogenes.
Fibroblast cell culture
Bovine fibroblasts isolated from skin were cultured in DMEM F12 (Thermo Scientific HyClone Laboratories, Logan, UT) supplemented with 15% fetal bovine serum (FBS) (Thermo Scientific HyClone Laboratories), 100 U/ml penicillin, and 100 mg/ml streptomycin. Cells were cultured at 37°C with 5% CO2. For cell collection fibroblasts were treated with 0.25% trypsin prior to collection for RNA and DNA isolation.
The Utah State University Institutional Animal Care and Use Committee approved all procedures for the use of animals in this study (protocol #1506). Bovine ovaries were collected at a local abattoir (E.A. Miller, Hyrum, UT) and used for collecting oocytes. Oocytes with 3 to 8 mm follicles were aspirated along with cumulus complexes. Following aspiration, cumulus oocyte complexes were cultured at 37°C with 5% CO2 for 18 to 22 hr in TCM 199 maturation medium containing 10% FBS (Thermo Scientific HyClone Laboratories), 0.05 mg/ml follicle stimulating hormone, 5 mg/ml luteinizing hormone, 100 U/ml penicillin, and 100 mg/ml streptomycin.
In vitro fertilization
Following 18 to 22 hr of oocyte maturation, cryopreserved bovine semen (Hoffman AI, Logan, UT) was thawed in a 37°C water bath. Live sperm were isolated by centrifugation through a 45%/90% Percoll® gradient. Sperm were suspended in Tyrode’s albumin lactate pyruvate (TALP) and used for oocyte fertilization. Twenty-four hr post fertilization, cumulus cells were removed by vortexing the cumulus oocyte complex in phosphate buffered saline (PBS) containing 0.32 mM sodium pyruvate, 5.55 mM glucose, 3 mg/ml BSA, and 10 mg/ml hyaluronidase. Oocytes were washed through six drops of PBS and placed in co-culture dishes plated with cultured cumulus cells and cultured in CR2 medium.
Tissue collection and RNA isolation
Twenty-five pooled blastocyst embryos were collected after 8 days of culture in CR2 medium, snap frozen with liquid nitrogen, and stored at -80°C. Twenty-five pooled oocytes were collected after 22 hr of maturation, vortexed for 5 min in PBS to remove cumulus cells, snap frozen with liquid nitrogen, and stored at -80°C. RNA was isolated from oocyte and blastocyst embryo samples using the RNeasy Mini Kit (Qiagen, Germantown, MD) following manufacturer’s instructions. Bovine brain, liver, testes tissue samples were collected immediately after slaughter and suspended in RNALater (Ambion, Austin, TX). Samples were stored overnight at 4°C. Fibroblasts were collected as described above. RNA was isolated from brain, liver, testes, and fibroblasts using TRIzol reagent (Life Technologies, Carlsbad, CA). Tissue samples were homogenized in 3 ml TRIzol with a tissue homogenizer. Cells were incubated for 5 min at room temperature, combined with 0.6 ml chloroform, and mixed by inversion. Samples were centrifuged at 12,000 × g for 15 min at 4°C, the upper aqueous phase was removed and combined with 1.5 ml isopropyl alcohol, then centrifuged at 12,000 × g for 10 min at 4°C. Supernatant was removed and the RNA was washed with 75% ethanol and centrifuged again, then dried and resuspended in H2O. Isolated RNA was immediately converted to cDNA using the Superscript III Reverse Transcriptase Kit (Life Technologies) following manufacturer’s protocol. Samples were stored at -20°C until use.
DNA isolation and bisulfite conversion
Twenty-five pooled oocytes and blastocyst embryos were snap frozen with liquid nitrogen and stored at -80°C until direct bisulfite conversion. DNA from brain, liver, testes, and fibroblast samples were collected from the interphase of the TRIzol treatment used for RNA collection. The collected interphase was combined with 0.9 ml ethanol, mixed by inversion, and centrifuged at 2000 × g for 5 min at 4°C. The pellet was suspended with 75% ethanol and washed by centrifugation three times. The DNA pellet was resuspended in 0.1 μM sodium citrate. The isolated DNA underwent bisulfite conversion using the EZ DNA Methylation Kit (Zymo Research, Irvine, CA) according to manufacturer’s recommendation. DNA was stored at -20°C until use. Bisulfite PCR and 454 sequencing
Primers for bisulfite sequencing
Primers for end-point RT-PCR
Sequencing data was analyzed to generate the frequency of methylation using BISMA analysis software . Any sequence with a conversion rate less than 98% was excluded from sequence analysis. The average conversion rate of all sequences was 99.2%. Data are presented as the frequency of methylation on a per CpG site basis (number of transcripts with methylated CpG site ÷ total number of transcripts sequenced × 100). In order to account for the potential of a C-to-T transition causing misinterpretation of our results, known C/T SNPs were identified in the Single Nucleotide Polymorphism database (http://www.ncbi.nlm.nih.gov/snp). The database identified two CpG sites that exist as either a cytosine or thymine in normal populations. CpG site 23 (rs43296056) in LIN28P-Ch:7 and CpG site 9 (rs208890804) in LIN28P-Ch:26 (Additional file 1: Figures S2-S3). The identified sites were excluded from all methylation analysis.
Financial support for this research was obtained from U.S. Department of Agriculture, National Institute of Food and Agriculture grant numbers 08-34526-19199 and 09-34526-19808 (to KLW). This research was also supported by the Utah Agricultural Experiment Station, Utah State University, and was approved as journal paper number 8711.
- Chen RZ, Pettersson U, Beard C, Jackson-Grusby L, Jaenisch R. DNA hypomethylation leads to elevated mutation rates. Nature. 1998;395(6697):89–93.View ArticlePubMedGoogle Scholar
- Csankovszki G, Nagy A, Jaenisch R. Synergism of Xist RNA, DNA methylation, and histone hypoacetylation in maintaining X chromosome inactivation. J Cell Biol. 2001;153(4):773–84.View ArticlePubMed CentralPubMedGoogle Scholar
- Mohandas T, Sparkes RS, Shapiro LJ. Reactivation of an inactive human X chromosome: evidence for X inactivation by DNA methylation. Science. 1981;211(4480):393–6.View ArticlePubMedGoogle Scholar
- Liu WM, Maraia RJ, Rubin CM, Schmid CW. Alu transcripts: cytoplasmic localisation and regulation by DNA methylation. Nucleic Acids Res. 1994;22(6):1087–95.View ArticlePubMed CentralPubMedGoogle Scholar
- Walsh CP, Chaillet JR, Bestor TH. Transcription of IAP endogenous retroviruses is constrained by cytosine methylation. Nat Genet. 1998;20(2):116–7.View ArticlePubMedGoogle Scholar
- Woodcock DM, Lawler CB, Linsenmeyer ME, Doherty JP, Warren WD. Asymmetric methylation in the hypermethylated CpG promoter region of the human L1 retrotransposon. J Biol Chem. 1997;272(12):7810–6.View ArticlePubMedGoogle Scholar
- Yoder JA, Walsh CP, Bestor TH. Cytosine methylation and the ecology of intragenomic parasites. Trends Genet. 1997;13(8):335–40.View ArticlePubMedGoogle Scholar
- De Smet C, De Backer O, Faraoni I, Lurquin C, Brasseur F, Boon T. The activation of human gene MAGE-1 in tumor cells is correlated with genome-wide demethylation. Proc Natl Acad Sci USA. 1996;93(14):7149–53.View ArticlePubMed CentralPubMedGoogle Scholar
- De Smet C, Lurquin C, Lethe B, Martelange V, Boon T. DNA methylation is the primary silencing mechanism for a set of germ line- and tumor-specific genes with a CpG-rich promoter. Mol Cell Biol. 1999;19(11):7327–35.PubMed CentralPubMedGoogle Scholar
- Jaenisch R, Bird A. Epigenetic regulation of gene expression: how the genome integrates intrinsic and environmental signals. Nat Genet. 2003;33(Suppl):245–54.View ArticlePubMedGoogle Scholar
- Shiota K. DNA methylation profiles of CpG islands for cellular differentiation and development in mammals. Cytogenet Genome Res. 2004;105(2–4):325–34.View ArticlePubMedGoogle Scholar
- Bartolomei MS, Webber AL, Brunkow ME, Tilghman SM. Epigenetic mechanisms underlying the imprinting of the mouse H19 gene. Genes Dev. 1993;7(9):1663–73.View ArticlePubMedGoogle Scholar
- Bird A. DNA methylation patterns and epigenetic memory. Genes Dev. 2002;16(1):6–21.View ArticlePubMedGoogle Scholar
- Li E, Beard C, Jaenisch R. Role for DNA methylation in genomic imprinting. Nature. 1993;366(6453):362–5.View ArticlePubMedGoogle Scholar
- Stoger R, Kubicka P, Liu CG, Kafri T, Razin A, Cedar H, et al. Maternal-specific methylation of the imprinted mouse Igf2r locus identifies the expressed locus as carrying the imprinting signal. Cell. 1993;73(1):61–71.View ArticlePubMedGoogle Scholar
- Eckhardt F, Lewin J, Cortese R, Rakyan VK, Attwood J, Burger M, et al. DNA methylation profiling of human chromosomes 6, 20 and 22. Nat Genet. 2006;38(12):1378–85.View ArticlePubMed CentralPubMedGoogle Scholar
- Illingworth R, Kerr A, Desousa D, Jorgensen H, Ellis P, Stalker J, et al. A novel CpG island set identifies tissue-specific methylation at developmental gene loci. PLoS Biol. 2008;6(1):e22.View ArticlePubMed CentralPubMedGoogle Scholar
- Holliday R, Pugh JE. DNA modification mechanisms and gene activity during development. Science. 1975;187(4173):226–32.View ArticlePubMedGoogle Scholar
- Jones PA, Liang G. Rethinking how DNA methylation patterns are maintained. Nat Rev Genet. 2009;10(11):805–11.View ArticlePubMed CentralPubMedGoogle Scholar
- Riggs AD. X chromosome inactivation, differentiation, and DNA methylation revisited, with a tribute to Susumu Ohno. Cytogenet Genome Res. 2002;99(1–4):17–24.View ArticlePubMedGoogle Scholar
- Goll MG, Bestor TH. Eukaryotic cytosine methyltransferases. Annu Rev Biochem. 2005;74:481–514.View ArticlePubMedGoogle Scholar
- Li E, Bestor TH, Jaenisch R. Targeted mutation of the DNA methyltransferase gene results in embryonic lethality. Cell. 1992;69(6):915–26.View ArticlePubMedGoogle Scholar
- Okano M, Bell DW, Haber DA, Li E. DNA methyltransferases Dnmt3a and Dnmt3b are essential for de novo methylation and mammalian development. Cell. 1999;99(3):247–57.View ArticlePubMedGoogle Scholar
- Grunau C, Hindermann W, Rosenthal A. Large-scale methylation analysis of human genomic DNA reveals tissue-specific differences between the methylation profiles of genes and pseudogenes. Hum Mol Genet. 2000;9(18):2651–63.View ArticlePubMedGoogle Scholar
- Rauch TA, Wu X, Zhong X, Riggs AD, Pfeifer GP. A human B cell methylome at 100-base pair resolution. Proc Natl Acad Sci USA. 2009;106(3):671–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Laurent L, Wong E, Li G, Huynh T, Tsirigos A, Ong CT, et al. Dynamic changes in the human methylome during differentiation. Genome Res. 2010;20(3):320–31.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhang X, Yazaki J, Sundaresan A, Cokus S, Chan SW, Chen H, et al. Genome-wide high-resolution mapping and functional analysis of DNA methylation in arabidopsis. Cell. 2006;126(6):1189–201.View ArticlePubMedGoogle Scholar
- Hesson LB, Packham D, Pontzer E, Funchain P, Eng C, Ward RL. A reinvestigation of somatic hypermethylation at the PTEN CpG island in cancer cell lines. Biol Proced Online. 2012;14(1):5.View ArticlePubMed CentralPubMedGoogle Scholar
- Zysman MA, Chapman WB, Bapat B. Considerations when analyzing the methylation status of PTEN tumor suppressor gene. Am J Pathol. 2002;160(3):795–800.View ArticlePubMed CentralPubMedGoogle Scholar
- Darr H, Benvenisty N. Genetic analysis of the role of the reprogramming gene LIN-28 in human embryonic stem cells. Stem Cells. 2009;27(2):352–62.View ArticlePubMedGoogle Scholar
- Yu J, Vodyanik MA, Smuga-Otto K, Antosiewicz-Bourget J, Frane JL, Tian S, et al. Induced pluripotent stem cell lines derived from human somatic cells. Science. 2007;318(5858):1917–20.View ArticlePubMedGoogle Scholar
- Cortese R, Krispin M, Weiss G, Berlin K, Eckhardt F. DNA methylation profiling of pseudogene-parental gene pairs and two gene families. Genomics. 2008;91(6):492–502.View ArticlePubMedGoogle Scholar
- Monk M, Boubelik M, Lehnert S. Temporal and regional changes in DNA methylation in the embryonic, extraembryonic and germ cell lineages during mouse embryo development. Development. 1987;99(3):371–82.PubMedGoogle Scholar
- Razin A, Webb C, Szyf M, Yisraeli J, Rosenthal A, Naveh-Many T, et al. Variations in DNA methylation during mouse cell differentiation in vivo and in vitro. Proc Natl Acad Sci USA. 1984;81(8):2275–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Hajkova P, Erhardt S, Lane N, Haaf T, El-Maarri O, Reik W, et al. Epigenetic reprogramming in mouse primordial germ cells. Mech Dev. 2002;117(1–2):15–23.View ArticlePubMedGoogle Scholar
- Mayer W, Niveleau A, Walter J, Fundele R, Haaf T. Demethylation of the zygotic paternal genome. Nature. 2000;403(6769):501–2.View ArticlePubMedGoogle Scholar
- Morgan HD, Santos F, Green K, Dean W, Reik W. Epigenetic reprogramming in mammals. Hum Mol Genet. 2005;14(Spec No 1):R47–58.View ArticlePubMedGoogle Scholar
- Oswald J, Engemann S, Lane N, Mayer W, Olek A, Fundele R, et al. Active demethylation of the paternal genome in the mouse zygote. Curr Biol. 2000;10(8):475–8.View ArticlePubMedGoogle Scholar
- Duret L, Chureau C, Samain S, Weissenbach J, Avner P. The Xist RNA gene evolved in eutherians by pseudogenization of a protein-coding gene. Science. 2006;312(5780):1653–5.View ArticlePubMedGoogle Scholar
- Guo X, Lin M, Rockowitz S, Lachman HM, Zheng D. Characterization of human pseudogene-derived non-coding RNAs for functional potential. PLoS One. 2014;9(4):e93972.View ArticlePubMed CentralPubMedGoogle Scholar
- Hawkins PG, Morris KV. Transcriptional regulation of Oct4 by a long non-coding RNA antisense to Oct4-pseudogene 5. Transcription. 2010;1(3):165–75.View ArticlePubMed CentralPubMedGoogle Scholar
- Hirotsune S, Yoshida N, Chen A, Garrett L, Sugiyama F, Takahashi S, et al. An expressed pseudogene regulates the messenger-RNA stability of its homologous coding gene. Nature. 2003;423(6935):91–6.View ArticlePubMedGoogle Scholar
- Korneev SA, Park JH, O'Shea M. Neuronal expression of neural nitric oxide synthase (nNOS) protein is suppressed by an antisense RNA transcribed from an NOS pseudogene. J Neurosci. 1999;19(18):7711–20.PubMedGoogle Scholar
- Poliseno L, Salmena L, Zhang J, Carver B, Haveman WJ, Pandolfi PP. A coding-independent function of gene and pseudogene mRNAs regulates tumour biology. Nature. 2010;465(7301):1033–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Patil V, Ward RL, Hesson LB. The evidence for functional non-CpG methylation in mammalian cells. Epigenetics. 2014;9(6):823–8.View ArticlePubMedGoogle Scholar
- Rohde C, Zhang Y, Reinhardt R, Jeltsch A. BISMA–fast and accurate bisulfite sequencing data analysis of individual clones from unique and repetitive sequences. BMC Bioinformatics. 2010;11:230.View ArticlePubMed CentralPubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.