Epigenetic response of imprinted domains during carcinogenesis
© The Author(s). 2017
Received: 8 May 2017
Accepted: 16 August 2017
Published: 25 August 2017
Imprinted domains have been identified as targets for aberrant DNA methylation during carcinogenesis, but it remains unclear when these epigenetic alterations occur and how they contribute to tumor progression. Epigenetic instability at key cis-regulatory elements within imprinted domains can concomitantly activate proto-oncogenes and turn off tumor suppressor genes. Thus, to further characterize the epigenetic response of imprinted domains during carcinogenesis, we compared the stability of DNA methylation at a variety of cis-regulatory elements within imprinted domains in two fundamentally different mouse tumors, benign and malignant, induced by the KrasG12D mutation.
We report that imprinted domains remain stable in benign processes but are highly susceptible to epigenetic alterations in infiltrative lesions. The preservation of DNA methylation within imprinted domains in benign tumors throughout their duration suggests that imprinted genes are not involved with the initiation of carcinogenesis or the growth of tumors. However, the frequent detection of DNA methylation changes at imprinting control regions in infiltrative lesions suggest that imprinted genes are associated with tumor cells gaining the ability to defy tissue boundaries.
Overall, this study demonstrates that imprinted domains are targeted for DNA hypermethylation when benign tumor cells transition to malignant. Thus, monitoring DNA methylation within imprinted domains may be useful in evaluating the progression of neoplasms.
1Imprinted genes are mainly expressed from a single allele based upon the allele’s parental origin . Thus far, there have been just over 100 imprinted genes identified by conventional methodology; however, a far greater number of imprinted genes may exist due to recent data obtained by implementing modern whole-transcriptome sequencing technologies [2, 3]. Nevertheless, the well-known imprinted genes are clustered in discrete chromosomal domains, which often contain a differentially methylated master cis-regulatory element termed an imprinting control region (ICR) [1, 4, 5]. ICRs are accompanied by other cis-regulatory elements such as enhancers and differentially methylated promoters that act synergistically within their respective imprinted domain to elegantly control imprinted gene dosage [1, 6]. While DNA methylation is transient among the other cis-regulatory elements depending on developmental contexts and cell types, DNA methylation at ICRs is exceedingly stable [1, 6]. The stability of this epigenetic mark at ICRs is paramount due to their control over several imprinted genes within a domain that have key roles in cell growth, division, or death [7–9]. Consistent with this, imprinted genes have gained much attention as both tumor suppressor genes and oncogenes [10–14]. Indeed, epigenetic perturbations within imprinted domains have been reported in a variety of both human and mouse malignant neoplasms [15, 16]. However, the timing in which epigenetic change occurs within imprinted domains, especially at ICRs, is still unclear.
DNA methylation is an important epigenetic regulator of gene expression. Misregulation of DNA methylation in tumor cells is well recognized as an epigenetic alteration that results in significant expression level differences of genes that contribute to the carcinogenic process [13, 14, 17–20]. The genome of any given normal cell by default is virtually stifled with DNA methylation, with the exception of specific gene promoters and cis-regulatory elements such as enhancers that contribute to the expression of genes necessary for the proper function of the respective cell. However, during carcinogenesis, the levels of DNA methylation are drastically reduced genome wide . There are two distinct drastic reductions in 5-methylcytosine that occur throughout the genome during the carcinogenic stages: first, during the transition of normal cells to immortalized tumor cells, and second, as benign tumor cells acquire infiltrative capacity . Interestingly, hypermethylation at specific regulatory elements such as tumor suppressor gene promoters occurs concomitantly with genome-wide hypomethylation during these critical transitional stages. Although ICRs have been well identified as targets for hypermethylation, the functional role that this epigenetic aberration plays during the carcinogenesis process remains largely uncharacterized [15, 16]. Does hypermethylation at ICRs contribute to the initiation of tumor formation, the proliferation of tumor cells, or the transition of a benign process to that of an infiltrative one?
In the current study, we exploited the KrasG12D genomic mutation to initiate carcinogenesis in mice and surveyed the epigenetic stability at ICRs in two fundamentally different neoplasms: squamous papilloma that remains a benign process but exhibits continual growth, and T cell lymphoblastic lymphoma that rapidly acquires infiltrative capacity. This comparison allowed further characterization of the functional role that aberrant DNA methylation at ICRs plays during the carcinogenic process. According to the results, ICRs remain epigenetically stable during the initiation of carcinogenesis and throughout the growth of tumors. However, DNA hypermethylation was observed across virtually all ICRs in infiltrative tumor cell populations regardless of their duration. These data suggest ICRs are targeted during the second phase of DNA methylation alterations to the genome and epigenetic instability among ICRs contributes to the infiltrative capacity of tumor cells. Detailed results and the relevant discussion are presented below.
The KrasG12D mutation is sufficient to initiate both an infiltrative T cell neoplasm and benign squamous cell papilloma [15, 24, 25]. Furthermore, the KrasG12D mutation has also been exploited to study non-small cell lung cancer, pancreatic cancer, and colorectal tumors. Previously, we reported the epigenetic instability at ICRs within T cell lymphoma driven by the KrasG12D mutation, which suggested that ICRs are targeted for hypermethylation in infiltrative T cell populations . With the current study, we sought to further support the previous report by comparing the T cell neoplasm to a benign or non-infiltrative tumor type. To accomplish this, we bred mice that express Cre under the Mouse Mammary Tumor Virus long terminal repeat promoter (MMTV-Cre) with the LSL-Kras G12D line to simultaneously target oncogenic KrasG12D to T cells and epithelial cells of the mucous membranes (Fig. 1b) [26, 27]. Using this breeding scheme, we obtained 80 mice from 11 litters with an average litter size of 7. Of the 80 total mice, 28 were positive for the recombined KrasG12D allele, which deviated significantly from the expected Mendelian ratio (P = 0.007, chi-square). All of the recombinant progeny developed squamous papilloma and nearly 40% developed thymic lymphoma as well as squamous papilloma (Fig. 1c, d) . Unlike the thymic lymphoma, which infiltrated and extensively replaced the surrounding tissues (Fig. 1e), the squamous papillomas were benign. The squamous papilloma lesions retained an exophytic growth pattern with an overall thickening of the epithelial layer into finger-like projections (Fig. 1f). In sum, a mouse breeding scheme utilizing KrasG12D/MMTV-Cre mice was successful in obtaining mice that developed both benign and infiltrative lesions.
Many efforts have been focused on determining the potential role that imprinted genes may play during carcinogenesis [10, 11, 15, 16, 21, 33, 34]. Analyses of human malignant neoplasms have revealed that many ICRs and putative enhancers within imprinted domains show significant epigenetic change pertaining to DNA methylation and that these changes often correlate with the expression level differences of nearby imprinted genes within the domain . However, the question whether this epigenetic instability within imprinted domains contributes to the infiltration potential of tumor cells remains elusive . In the current study, we initiated tumorigenesis in mice with the KrasG12D mutation and subsequently compared the epigenetic stability of various regulatory sequence elements within imprinted domains between two fundamentally different neoplasms, benign squamous papilloma, and malignant T cell lymphoma. Although there is some overlap with our previous report on the T cell lymphoma, we further show with this report that there is a stark contrast in the epigenetic stability at ICRs between the benign and infiltrative lesions studied: while DNA methylation among ICRs was stable in the benign lesions, including in the thymic lymphoid hyperplastic lesions, it was highly unstable in the infiltrative thymic T cell lesions. The epigenetic instability at imprinting control regions was also seen in our previous report on the infiltrative thymic lymphoma, which were isolated from mice using a different breeding scheme involving a conditional knockout construct for Peg3 that was bred into the KrasG12D model. Thus, it should be noted that even though the tumors from each study were isolated from mice of two different breeding schemes, the thymic tumors in both studies were of the same genotype (MMTV-Cre/KrasG12D). We also expanded our analysis of imprinted domains in the current study and show that DNA methylation is unstable at certain regulatory sequence elements specific to their imprinted domain in both the benign and malignant lesions. Overall, it appears that imprinted domains are most vulnerable to epigenetic change in infiltrative tumor cells.
DNA methylation settings at ICRs are critical for maintaining parent-of-origin specific expression of many genes within an imprinted domain . ICRs contain various sequence elements that commonly manifest aberrant DNA methylation trends in tumor cells ; thus, the stability of DNA methylation at ICRs is concerning in the context of carcinogenesis. Moreover, ICRs are susceptible to the Knudson two-hit hypothesis because one allele is already methylated and silenced in the native state. Therefore, ICRs can potentially respond to either hypomethylation events or hypermethylation events. Indeed, the Knudson two-hit hypothesis has been well supported with hypermethylation of imprinted genes in malignant tumors [13, 15, 16]. However, the timing when ICRs manifest epigenetic change, particularly DNA methylation changes, is debatable. First, do epigenetic changes at ICRs coincide with the initiation of carcinogenesis or, second, do epigenetic changes at ICRs coincide with neoplastic transformation of tumor cells?
Epigenetic stability at ICRs is challenged during the second reconstruction of the DNA methylome (Fig. 5). This wave of DNA methylation changes is associated with increase in the infiltrative potential of tumor cells. This is demonstrated by the results from the T cell lymphoblastic lymphoma panel (Fig. 2). In addition to neoplastic thymic lymphoid lesions, this panel also included hyperplastic and preneoplastic thymic lymphoid lesions . ICRs remained epigenetically stable in the hyperplastic and preneoplastic specimens regardless of the duration of the lesions. In contrast, malignant neoplastic specimens displayed marked epigenetic instability with nearly all the ICRs becoming hypermethylated, regardless of the duration of the lesion. ICRs acquiring hypermethylation do not necessarily confer repression of imprinted genes as seen with hypermethylation in tumor suppressor gene promoters . For instance, hypermethylation at the ICR of the H19/Igf2 domain can potentially result in complete repression of H19 and bi-allelic expression of Igf2, a known oncogene, whereas hypermethylation at the ICR of the Peg3 domain can potentially result in repression of Peg3, a putative tumor suppressor . Thus, concomitant hits to multiple ICRs can simultaneously activate oncogenes while suppressing tumor suppressor genes. Overall, it appears ICRs are targeted to be hypermethylated in the midst of the second drastic genome-wide reduction in 5-methylcytosine, a hallmark of tumor progression to malignant disease. This raises the notion that epigenetic instability among ICRs could enhance the infiltrative capacity of tumor cells.
The epigenetic stability among the ICRs in the benign tumors is not a common theme among other regulatory sequences within their respective imprinted domain. This is demonstrated by our analysis of the H19/Igf2 and Peg3 imprinted domains (Fig. 4). For instance, the Zfp264 promoter within the Peg3 imprinted domain became hypermethylated in several benign specimens while the ICR retained normal DNA methylation patterns. This observation demonstrates the specificity of DNA methylation changes among sequence elements within imprinted domains during tumor progression. The observed specificity of DNA methylation changes within a single imprinted domain further supports the notion that epigenetic changes at ICRs could potentially be involved in enhancing the infiltrative potential of tumor cells as ICRs are only targeted during the second reconstruction of the DNA methylome.
There are limitations to consider with the methodology used to measure DNA methylation for the current study. First, COBRA is not a very sensitive method and it only measures the methylation at one CpG site at a time. Although densitometry can be performed to quantify COBRA data as we did in the current study, COBRA is also inherently qualitative. To address these issues and to validate the results from COBRA, we performed bisulfite sequencing of the Peg3 bisulfite PCR products. The results from sequencing reinforced the results from COBRA; however, there are also limitations to consider with traditional bisulfite sequencing due to the low number of reads compared to Mass Array or Next Generation sequencing approaches. Again, caution should be used when interpreting the results from bisulfite sequencing, as they are more qualitative than quantitative. For studies that require more quantitative methodologies such as in the case with human cancer studies, it is advised to employ more sensitive methods that utilize mass array and/or next generation-based sequencing technology.
To our knowledge, this is the first study to compare the epigenetic stability among ICRs between benign and malignant tumor cells driven by the same genomic mutation, KrasG12D. The results from the DNA methylation analyses of this study demonstrate that ICRs remain stable during the first reconstruction of the DNA methylation landscape but are targeted specifically by DNA hypermethylation during the second reconstruction of the DNA methylation landscape—a hallmark that defines the transition of benign tumor cells to malignant. Further investigation is required to determine whether the observed changes in DNA methylation at ICRs directly enhance the infiltrative potential of tumor cells. However, considering the effects that aberrant DNA methylation at an ICR can have on an entire imprinted domain, it is likely that concomitant hits to several ICRs can result in an aberrant protein environment within a tumor cell whereby several tumor suppressors are shut down and several oncogenes are activated. These dynamic changes in critical imprinted gene products may very well enhance a tumor cell’s ability to infiltrate tissue boundaries. Lastly, monitoring DNA methylation at ICRs may be a useful diagnostic in the clinical setting in determining the progression of a neoplasm.
Mouse strains and breeding
Two mouse strains were purchased from Jackson laboratories: B6.129-Krastm4Tyj/Nci (LSL- KrasG12D) and STOCK Tg(MMTV-Cre)4Mam/J (MMTV-Cre) [22, 23]. The LSL-KrasG12D mice were maintained as heterozygotes and bred with homozygous MMTV-Cre mice.
Mice weaned at 21 days postpartum were separated by sex and marked by a hole punch with varying positions in the left ear. Ear clips were taken from the right ear for genomic DNA isolation and subsequent PCR analysis. Sample preparation and PCR protocols have been previously described . Primer sets for genotyping can be found in Additional file 5: Table S2.
Necropsy and histopathology
Mice were monitored daily for signs of distress or weight loss equaling 15% of total body weight and sacrificed according to the guidelines set forth in the IACUC protocol #16–060. Upon euthanasia, mice were submitted to full necropsy. Samples of squamous papilloma and thymic lesions were placed directly in tail lysis buffer with Proteinase K for subsequent DNA isolation. Detailed buffer compositions and DNA extraction protocols were previously reported . The remaining portion of these lesions and representative samples of all other tissues were fixed for at least 48 h in 10% neutral buffered formalin (Thermo Scientific, Cat. # 5725) and then transferred to 70% ethanol. Tissues were further processed according to standard protocols for hematoxylin and eosin (H&E) staining and examined by a board-certified pathologist (IML).
DNA methylation analysis
Combined Bisulfite Restriction Analysis (COBRA) and bisulfite sequencing was performed for DNA methylation analyses . Genomic DNA from the thymic lymphoma and squamous papilloma lesions was isolated using a commercial kit (Genomic DNA clean and concentrator-25, Zymo Research, Cat. No. D4065), and subsequently, bisulfite-converted using another commercial kit (EZ DNA methylation kit, Zymo Research, Cat. No. D5002). Bisulfite-converted DNA was then used as template for PCR. Additional file 5: Table S2 contains detailed information pertaining to the genomic locations analyzed and their respective oligonucleotide bisulfite PCR primer sets. Each bisulfite PCR product was digested with at least one restriction enzyme. The resulting fragments were visualized using gel electrophoresis, and band densities were calculates as described below. Images were exported as tiff files into ImageJ software. Briefly, data was inverted, the background was removed, the brightness/contrast was adjusted, the bands were selected using the rectangular tool to generate density plots, density peaks were gated using the line tool, and the area of each peak was then calculated using the wand tool. The results were then copied into a Microsoft Excel spreadsheet for further processing and statistics. DNA methylation values were calculated in the following manner: 100*((area of peak from digested DNA/s)/(area of peak from digested band/s + area of peak from undigested DNA)). These % methylation values were then used for statistical analyses. First, a single factor ANOVA test was performed to identify significant differences among the sample. Upon a significant P value from ANOVA, the means of the tumor samples were then compared to the normal sample using the Student t test (two sample assuming equal variance). Three independent technical replicates were performed for each sample and each primer set. For the squamous papilloma panel, 15 samples were analyzed initially; however, only seven samples were used for statistical comparisons based on the consistent results from successful bisulfite PCRs across all primer sets (Additional file 6: Table S3). For bisulfite sequencing, PCR products were individually cloned using a commercial vector system (pGEM®-T Easy Vector System I, Promega, A1360). Ten clones were selected for sequencing. In brief, plasmids were purified using a commercial kit (DNA-spin plasmid DNA purification kit, iNtRON Biotechnology, 17097), restriction digest was performed to check for insert, BigDye Terminator v3.1 reactions were performed, products were precipitated and re-suspended in ABI Hi-Di formamide and sequences using a 3130XL Genetic Analyzer. Sequencing results were then manually processed to remove vector sequence, placed in the same orientation, and aligned using ClustalW. The methylation of each CpG site in the sequencing reads was then manually scored in an Excel spreadsheet. Sequencing reads with less than 95% bisulfite conversion efficiency were removed.
We would like to thank the members of the histology core facility at the Louisiana Animal Disease Diagnostic Laboratory (LADDL) for their help with processing and embedding the tissues.
This research was supported by the National Institutes of Health (J.K. R01-GM066225 and R01-GM 097074).
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its supplementary information files.
CLB designed and executed the breeding experiments for the study, conducted all necropsies, performed all molecular and statistical analyses, and wrote the paper. IML performed the histopathological analyses and contributed to writing the paper. JK provided the original concept of the study, supervised the study, and contributed to writing the paper. All authors have read and approved the final manuscript.
Ethics approval and consent to participate
All the experiments related to mice were performed in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals and also approved by the Louisiana State University Institutional Animal Care and Use Committee (IACUC), protocol #16–060.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
- Bartolomei MS, Ferguson-Smith AC. (2011) Mammalian genomic imprinting. Cold Spring Harb Perspect Biol 3: PMID: 21576252. doi:https://doi.org/10.1101/cshperspect.a002592.
- Gregg C, Zhang J, Weissbourd B, Luo S, Schroth GP, et al. (2010) High-resolution analysis of parent-of-origin allelic expression in the mouse brain. Science 329: 643–648. PMID: 20616232. doi:https://doi.org/10.1126/science.1190830.
- Kelsey G, Bartolomei MS. Imprinted genes … and the number is? PLoS Genet. 2012;8(3):e1002601. https://doi.org/10.1371/journal.pgen.1002601.View ArticlePubMedPubMed CentralGoogle Scholar
- Bartolomei MS. (2009) Genomic imprinting: employing and avoiding epigenetic processes. Genes Dev 23(18): 2124–2133. PMID: 19759261. doi:https://doi.org/10.1101/gad.1841409.
- Spahn L, Barlow DP. An ICE pattern crystallizes. Nat Genet. 2003;35:11–12. PMCID: PMC2847178. DOI:10.1038/ng0903-11.
- Barlow DP, Bartolomei MS. (2014) Genomic imprinting in mammals. Cold Spring Harb Perspect Biol 6. PMID: 9442905. doi:https://doi.org/10.1101/cshperspect.a018382.
- Relaix F, Wei XJ, Li W, Pan J, Lin Y, Bowtell DD, et al. (2000) Pw1/Peg3 is a potential cell death mediator and cooperates with Siah1a in p53-mediated apoptosis. Proc Natl Acad Sci U S A 97(5): 2105–2110. PMID: 10681424. doi:https://doi.org/10.1073/pnas.040378897.
- Reik W, Davies K, Dean W, Kelsey G, Constância M. (2001) Imprinted genes and the coordination of fetal and postnatal growth in mammals. Novartis Found Symp 237: 19-31; discussion 31-42. PMID: 11444044. DOI:https://doi.org/10.1002/0470846666.ch3.
- Rezvani G, Lui JC, Barnes KM, Baron J. (2012) A set of imprinted genes required for normal body growth also promotes growth of Rhabdomyosarcoma cells. Pediatr Res 71(1): 32-38. PMID: 22289848. doi:https://doi.org/10.1038/pr.2011.6.
- Lim DH, Maher ER. (2010) Genomic imprinting syndromes and cancer. Adv Genet 70:145-175. PMID: 20920748. doi:https://doi.org/10.1016/B978-0-12-380866-0.60006-X.
- Uribe-Lewis S, Woodfine K, Stojic L, Murrell A. (2011) Molecular mechanisms of genomic imprinting and clinical implications for cancer. Expert Rev Mol Med 13: e2. PMID: 21262060. doi:https://doi.org/10.1017/S1462399410001717.
- Monk D. (2010) Deciphering the cancer imprintome. Brief Funct Genomics 9(4): 329-339. PMID: 20551091. doi:https://doi.org/10.1093/bfgp/elq013.
- Plass C, Soloway PD. (2002) DNA methylation, imprinting and cancer. Eur J Hum Genet 10(1): 6-16. PMID: 11896451. doi:https://doi.org/10.1038/sj/ejhg/5200768.
- Jelinic P, Shaw P. (2007) Loss of imprinting and cancer. J Pathol 211(3): 261-268. PMID: 17177177. 96451. doi:https://doi.org/10.1038/sj/ejhg/5200768 doi:https://doi.org/10.1101/gad.943001.
- Bretz CL, Langohr IM, Lee S, Kim J. (2015) Epigenetic instability at imprinting control regions in a Kras(G12D)-induced T-cell neoplasm. Epigenetics 10(12): 1111-1120. PMID 26507119. doi:https://doi.org/10.1080/15592294.2015.1110672.
- Kim J, Bretz CL, Lee S. (2015) Epigenetic instability of imprinted genes in human cancers. Nucleic Acids Res 43(22): 10689-10699. PMID 26338779. doi:https://doi.org/10.1093/nar/gkv867.
- Aran D, Sabato S, Hellman A. (2013) DNA methylation of distal regulatory sites characterizes dysregulation of cancer genes. Genome Biol 14(3): R21. PMID: 23497655. doi:https://doi.org/10.1186/gb-2013-14-3-r21.
- Ehrlich M. (2002) DNA methylation in cancer: too much, but also too little. Oncogene 21(35): 5400-5413. PMID: 12154403. doi:https://doi.org/10.1038/sj.onc.1205651.
- Sharma S, Kelly TK, Jones PA. Epigenetics in cancer. Carcinogenesis. 31(1): 27–36. PMID: 19752007. doi:https://doi.org/10.1093/carcin/bgp220.
- Gokul, G, and Sanjeev K. (2013) DNA methylation and cancer. Epigenetics: Development and Disease. Springer Netherlands, 597-625. doi:https://doi.org/10.1007/978-94-007-4525-4_26.
- Fraga MF, Herranz M, Espada J, Ballestar E, Paz MF et al. (2004) A mouse skin multistage carcinogenesis model reflects the aberrant DNA methylation patterns of human tumors. Cancer Res 64(16): 5527-5534. PMID: 15313885. doi:https://doi.org/10.1158/0008-5472.CAN-03-4061.
- Jackson EL, Willis N, Mercer K, Bronson RT, Crowley D, Montoya R, et al. (2001) Analysis of lung tumor initiation and progression using conditional expression of oncogenic K-ras. Genes Dev 15: 3243-3248. PMID: 11751630. doi:https://doi.org/10.1101/gad.943001.
- Johnson, L., Mercer, K., Greenbaum, D., Bronson, R.T., Crowley, D., Tuveson, D.A., et al. (2001) Somatic activation of the K-ras oncogene causes early onset lung cancer in mice. Nature 410: 1111– 1116. PMID: 11323676. doi:https://doi.org/10.1038/35074129.
- Chan IT, Kutok JL, Williams IR, Cohen S, Kelly L, Shigematsu H, et al. (2004) Conditional expression of oncogenic K-ras from its endogenous promoter induces a myeloproliferative disease. J Clin Invest 113: 528-538. PMID: 14966562. doi:https://doi.org/10.1172/JCI200420476.
- Kong G, Du J, Liu Y, Meline B, Chang YI, Ranheim EA, et al. (2013) Notch1 gene mutations target KRAS G12D-expressing CD8+ cells and contribute to their leukemogenic transformation. J Biol Chem 288(25):18219-18227. PMID: 23673656. doi:https://doi.org/10.1074/jbc.M113.475376.
- Wagner K, Wall RJ, St-Onge L, Gruss P, Wynshaw-Boris A, Garrett L, Li M, et al. (1997) Cre-mediated gene deletion in the mammary gland. Nuc Acids Res 25: 4323–4330. PMID: 9336464. doi:https://doi.org/10.1172/JCI200420476.
- Wagner KU, McAllister K, Ward T, Davis B, Wiseman R, Hennighausen L. (2001) Spatial and temporal expression of the Cre gene under the control of the MMTV-LTR in different lines of transgenic mice. Transgenic Res 10(6) : 545-553. PMID: 11817542. doi:https://doi.org/10.1023/A:1013063514007.
- Xiong Z, Laird PW. (1997) COBRA: a sensitive and quantitative DNA methylation assay. Nucleic Acids Res 25: 2532-2534. PMID: 9171110. doi:https://doi.org/10.1093/nar/25.12.2532.
- Shiura H, Nakamura K, Hikichi T, Hino T, Oda K, Suzuki-Migishima R, Kohda T et al. (2009) Paternal deletion of Meg1/Grb10 DMR causes maternalization of the Meg1/Grb10 cluster in mouse proximal chromosome 11 leading to severe pre- and postnatal growth retardation. Hum Mol Genet 18(8): 1424-1438. PMID: 19174477. doi:https://doi.org/10.1093/hmg/ddp049.
- Quenneville S, Verde G, Corsinotti A, Kapopoulou A, Jakobsson J, Offner S, Baglivo I, et al. (2011) In Embryonic stem cells, ZFP57/KAP1 recognize a Methylated Hexanucleotide to affect chromatin and DNA Methylation of imprinting control regions. Mol Cell. 44(3): 361-372. PMID: 22055183. https://doi.org/10.1016/j.molcel.2011.08.032
- Pearse G. (2006) Histopathology of the thymus. Toxicol Pathol 34(5):515-547. PMID: 17067942. doi:https://doi.org/10.1080/01926230600978458.
- Ward JM. Lymphomas and leukemias in mice. Exp Toxicol Pathol 57(5–6):377-381. PMID: 16713211. doi:https://doi.org/10.1016/j.etp.2006.01.007.
- Barrow TM, Barault L, Ellsworth RE, Harris HR, Binder AM, Valente AL, et al. (2015) Aberrant methylation of imprinted genes is associated with negative hormone receptor status in invasive breast cancer. Int J Cancer 137(3): 537-547. PMID: 25560175. doi:https://doi.org/10.1002/ijc.29419.
- Aran D, Hellman A. (2013) DNA methylation of transcriptional enhancers and cancer predisposition. Cell 154(1):11-13. PMID: 23827668. doi:https://doi.org/10.1016/j.cell.2013.06.018.