Email updates

Keep up to date with the latest news and content from BMC Genetics and BioMed Central.

Open Access Highly Accessed Methodology article

How to interpret Methylation Sensitive Amplified Polymorphism (MSAP) profiles?

Jaroslav Fulneček* and Aleš Kovařík

Author Affiliations

Institute of Biophysics, Academy of Sciences of the Czech Republic, v.v.i., Královopolská 135, Brno CZ-612 65, Czech Republic

For all author emails, please log on.

BMC Genetics 2014, 15:2  doi:10.1186/1471-2156-15-2

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2156/15/2


Received:30 July 2013
Accepted:10 December 2013
Published:6 January 2014

© 2014 Fulneček and Kovařík; licensee BioMed Central Ltd.

This is an open access article distributed under the terms of the Creative Commons Attribution License ( http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

DNA methylation plays a key role in development, contributes to genome stability, and may also respond to external factors supporting adaptation and evolution. To connect different types of stimuli with particular biological processes, identifying genome regions with altered 5-methylcytosine distribution at a genome-wide scale is important. Many researchers are using the simple, reliable, and relatively inexpensive Methylation Sensitive Amplified Polymorphism (MSAP) method that is particularly useful in studies of epigenetic variation. However, electrophoretic patterns produced by the method are rather difficult to interpret, particularly when MspI and HpaII isoschizomers are used because these enzymes are methylation-sensitive, and any C within the CCGG recognition motif can be methylated in plant DNA.

Results

Here, we evaluate MSAP patterns with respect to current knowledge of the enzyme activities and the level and distribution of 5-methylcytosine in plant and vertebrate genomes. We discuss potential caveats related to complex MSAP patterns and provide clues regarding how to interpret them. We further show that addition of combined HpaII + MspI digestion would assist in the interpretation of the most controversial MSAP pattern represented by the signal in the HpaII but not in the MspI profile.

Conclusions

We recommend modification of the MSAP protocol that definitely discerns between putative hemimethylated mCCGG and internal CmCGG sites. We believe that our view and the simple improvement will assist in correct MSAP data interpretation.

Keywords:
MSAP; DNA methylation; Methylcytosine; MspI; HpaII; Hemimethylated site; Data interpretation

Background

The DNA of most eukaryotic organisms contains 5-methylcytosine (mC) residues, which represent important epigenetic information involved in the regulation of gene expression during various developmental processes such as cell differentiation, imprinting, or X chromosome inactivation [1-3]. DNA methylation also contributes to genome stability by repressing harmful genetic elements [4]. The DNA methylation pattern may be changed by diet or stressful external conditions [5,6]. In plants, epialleles occurring naturally or after genetically/chemically induced epimutations may produce heritable phenotypic diversity supporting adaptation and evolution [7].

In a genome, only a fraction of cytosine residues is postreplicatively methylated by DNA (cytosine-5) methyltransferases [8]. The level of mC varies among individual species, and mC distribution along DNA is not uniform. Within an individual, the tissues or cell types differ in mC distribution. Repetitive sequences usually contain more mC than genic sequences [9]. Methylation of promoters leads to transcriptional inactivation of linked genes in most cases; the function of evolutionary conserved gene body methylation is still unknown [10,11]. Vertebrate DNA is predominantly methylated in CG dinucleotides [12-14]; however, in plant DNA, a cytosine in any sequence context can be methylated [15] with decreased frequency for CG, CHG, and CHH motifs (H = A or T or C; [16]).

To identify changes in the methylation of genomic DNA connected with biological processes or different types of treatments, a convenient method for mC detection should be used. Recently, several genome-wide methods were developed to analyze mC [17] based on bisulfite modification [16] or immunoprecipitation [18]. Because these methods are entirely dependent on detailed knowledge of the genome sequence, many scientists use Methylation Sensitive Amplified Polymorphism (MSAP), particularly in non-model organisms. The MSAP method based on digestion with methylation-sensitive restriction endonucleases followed by amplification of restriction fragments is independent on the availability of genome sequence information and has been used frequently to analyze DNA methylation changes in plants and animals [19-21]. The MSAP procedure was first introduced by Reyna-Lopez and co-workers [22]. In their modification of the original amplified fragment length polymorphism protocol [23], frequent cutter MseI was replaced by methylation-sensitive MspI and HpaII restriction enzymes. However, the choice of these particular enzymes may lead to ambiguous interpretation of MSAP data, and some findings may be inconsistent with our current knowledge of mammalian and plant DNA methylation. Here, we reviewed 16 years of experience with the MSAP technique and attempted to explain the ambiguous characteristics of some restriction patterns considering the methylation sensitivities of restriction enzymes and known methylation frequencies of various sequence motifs in vertebrate and plant genomes. We present a modification of the MSAP protocol that would help to interpret the patterns of the most commonly used enzymes, MspI and HpaII.

Results and discussion

Principle of MSAP

In the original MSAP protocol [22], genomic DNA is divided into two parts, and each part is digested with EcoRI, which recognizes the GAATTC target site and is thought to be negligibly influenced by DNA cytosine methylation, and one of the methylation-sensitive MspI or HpaII isoschizomers, which can both cleave CCGG sequences. MspI can cleave non-methylated CCGG sequences and hemi (mC in one DNA strand only) or fully methylated CmCGG sequences but not hemi and fully methylated mCCGG and mCmCGG sequences [22,24]. HpaII is presumed to digest only non-methylated CCGG sequences and hemimethylated mCCGG sequences from all possible methylated CCGG variants [22,25]. The DNA samples digested with EcoRI and MspI or EcoRI and HpaII are ligated to two dsDNA adapters compatible with EcoRI- and MspI/HpaII-generated ends. Ligated fragments are preamplified using non-selective or pre-selective primers complementary to the adapters followed by amplification with a pair of selective primers (these are one- to three-base extended variants of non-selective or pre-selective primers at 3′ ends). Such amplification produces a reduced population of fragments that are separated in denaturing polyacrylamide gels. Visualization of fragments is usually accomplished by the introduction of a radioactively or fluorescently-labeled primer during a selective amplification step. In the gel, each DNA sample is represented by two neighboring lanes of fragments resulting from EcoRI, MspI- and EcoRI, HpaII-digested DNA (M and H lane, respectively). Thus, four MH fragment pattern variants "+ +, - -, + -, - +", referring to the presence (+) or absence (-) of a fragment, can be observed for each position in the gel. The (+, +) pattern (a fragment of definite length visualized in both the M and H lanes) is attributed to digestion by both enzymes at a non-methylated CCGG site. The (-, -) pattern indicates inhibition of digestion with both enzymes at a fully methylated mCmCGG site when another treated sample shows the presence of a fragment at that position. The (-, -) pattern may also represent a mutated site when genetically distinct samples are compared. The (+, -) pattern representing a fragment of definite length visualized in the M but not in the H lane corresponds to digestion with MspI but not HpaII and refers to the presence of a CmCGG site. Finally, the (-, +) pattern seems to be quite difficult to interpret, particularly in plant genomes where methylation occurs in both the CCG and CG motifs.

Methylation sensitivity of restriction enzymes

The interpretation of MSAP data is predominantly based on known restriction enzyme (RE) activities at recognition sequences modified by methylation. Data concerning the methylation sensitivity of REs and corresponding literature can be found at the website of The Restriction Enzyme Database (REBASE; http://rebase.neb.com webcite; [26]). At REBASE, 18 different references were identified concerning MspI and/or HpaII enzyme activities at the recognition sequence containing 5-mC residues (13 for MspI and 14 for HpaII; http://rebase.neb.com webcite).

The CCGG recognition site can theoretically occur in 10 differentially methylated forms (Table  1), and not all of them were tested in MspI/HpaII activity assays. In addition, different methylated forms may provide identical MspI and HpaII digestion profiles, suggesting only limited ability of the MSAP method to detect all methylation changes at CCGG sites. As follows from the function of bacterial HpaII and MspI restriction-modification systems, CCGG sites should be resistant to digestion when the inner cytosines (for HpaII) or outer cytosines (for MspI) are methylated by cognate methyltransferases [27,28]. Simultaneously, postreplicatively generated sites methylated only at one DNA strand (termed hemimethylated) should also be resistant to digestion. It is known that MspI can digest HpaII-methylated sites (CmCGG) [29,30]; however, HpaII cannot digest MspI-methylated sites (mCCGG) [24,31]. Inconsistency has been found regarding the sensitivity of HpaII to methylation of the outer C (mCCGG) in the hemimethylated state. While earlier work reported cleavage of hemimethylated mCCGG sites with HpaII [25], other studies did not completely support such activity [24,31-33]. This discrepancy may be explained by the known slow HpaII nicking activity of the non-methylated DNA strand [24] possibly combining with impurity-driven nonspecific nuclease activity in some samples. Such DNA degradation might be supported by the observed lower intensity of HpaII digestion products compared with non-digested or HpaII-methylated DNA templates [25]. Indeed, Korch and Hagblom found smearing of HpaII-digested products after longer incubation or a higher enzyme concentration and suggested contaminating unspecific nuclease(s) in all HpaII enzyme preparations they used [32]. The current view supports HpaII nicking activity of the non-methylated but not the methylated strand in hemimethylated mCCGG sequences [31]. However, there is continuing controversy that the specificity of the HpaII enzyme and its ability to digest hemimethylated mCCGG sequences may be dependent on the enzyme concentration and incubation time.

Table 1. Frequency, effect on cleavage and expected MSAP profile of ten different methylated forms of CCGG site

Although the EcoRI enzyme recognizing GAATTC is used as an enzyme that is not sensitive to methylation in MSAP, EcoRI did not digest GAATTmC sequences when cytosines at both strands are methylated [34,35]. Thus, EcoRI digestion may be inhibited by overlapping C methylation, indicating that not all changes in MH patterns always reflect changes in the methylation status of MspI/HpaII sites.

Frequency of mC at CCGG sites in plant and vertebrate genomes

Despite earlier reports regarding the high frequency of mCCG methylation in plant genomes, [36] the more recent data support the following frequencies: CmCGG > mCmCGG> > mCCGG [16,37] which is in accord with DNA methylation and demethylation studies in plants [38,39]. Although the expected frequency of MspI/HpaII sites is one per 256 bp, the exact spacing is highly variable depending mostly on the GC content of a particular chromosomal region. In regions with closely spaced CCGG sequences, the supposed internal CmCGG site(s) of the HpaII-EcoRI fragment would be digested by MspI, and the resulting shorter MspI-EcoRI fragment may be lost during selective amplification or may be too short for electrophoretic detection. This finding is consistent with the observation that most (-, +) MH type fragments are long, while most (+, -) MH type fragments are typically short [40]. If so, all isolated HpaII fragments of the (-, +) MH type should contain at least one additional internal CCGG site; indeed, this was reported for all such sequenced fragments [41]. Internal CCGG sequences were also identified in HpaII-specific fragments obtained from restriction landmark genome scanning of A. thaliana and rice, and their CmCGG methylation status was confirmed by Southern blot hybridization and/or by polymerase chain reaction (PCR) [42,43]. More evidence for the presence of internal CmCGG site(s) in the fragments of the (-, +) MH type comes from MSAP analysis of seven swine tissues where approximately 26.2% and 25.7% of MH patterns were of the (+, -) and (-, +) type, respectively [21]. Considering that vertebrate genomes are methylated almost exclusively at CG sites and a nearly equal number of (+, -) and (-, +) MH profiles, we may expect an identical origin of these two patterns; essentially, the occurrence of CmCGG sequences only. This presumption is certainly valid also for the other two vertebrate species analyzed by MSAP [44,45].

Hyper- and hypomethylation may not be easily deciphered from MSAP patterns

Particularly for sequences with closely spaced MspI/HpaII sites, up to half of the different MSAP pattern changes are ambiguous and potentially misinterpreted either at the level of methylation change or as a change in nucleotide sequence (Table  2). Considering a model DNA molecule with closely spaced CCGG sites occurring in five different methylation variants (Figure  1), variants 1–3 are found in both vertebrate and plant DNA while variants 4 and 5 are plant specific. First, we would expect amplification of the long fragment (L) due to the specific combination of preselective and selective primers. Variant 2 will show a (-, +) MH pattern (due to an internal CmCGG site) that will be changed either by hypomethylation (variant 1) or hypermethylation (variant 3) into a (-, -) MH pattern. A shift from the (-, +) to the (+, +) MH pattern is not expected in vertebrate DNA samples; however, it may be quite frequent in plant samples (variant 4) and is caused by methylation of the outer cytosine in the internal CmCGG site. However, the most commonly used interpretation of the shift from the (-, +) to (+, +) MH pattern is hypomethylation of a putative mCCGG site.

Table 2. Possible changes in MSAP MH profiles from Figure1and their interpretation

thumbnailFigure 1. DNA with two closely spaced MspI/HpaII sites in five methylation states and the expected MSAP signal. A model DNA molecule with two closely spaced CCGG sites is depicted in five different methylated variants. To simplify, only one DNA strand is shown, and a symmetrical methylation pattern is expected. Variants 1–3 are found in both vertebrate and plant DNA, whereas variants 4 and 5 are plant specific. The methylation state of the two CCGG sites (Site s and l) influences digestion by MspI and HpaII. The expected MSAP MH patterns of the generated short (S) and long (L) fragments (Frag.) are drawn in boxes for the five methylated variants. Note that digested S and L DNA fragments may be amplified by different selective primers, and, if so, they are not visible together on one MSAP profile.

A modified MSAP protocol improves problematic interpretation of the (-, +) MH pattern

As stated above, in plant genomes, the (-, +) MH pattern seems to be particularly problematic because it may be interpreted as two different situations: (i) the cutting of the hemimethylated mCCGG sites with HpaII but not MspI and (ii) the presence of internal CmCGG site(s) between the cleaved distal CCGG and the EcoRI site (Figure  2). Correct interpretation of the (-, +) MH pattern is important for classification of methylation changes that may involve both hypermethylation and hypomethylation. Here we show that introduction of a "third" lane containing DNA digested with both MspI and HpaII enzymes (M + H lane) may be highly informative. The persistence of the "H only" band leading to the (-, +, +) (M, H, M + H, in this order) pattern would support hemimethylation of the mCCGG. By contrast, loss of the "H only" band (-, +, -) indicates the presence of an internal CmCGG site. We decided to test this theoretical construct by analysis of control and hypomethylated tobacco DNA (Figure  3, Additional file 1: Figure S1). A typical (-, +) MH pattern is shown in Figure  3 containing two HpaII fragments (1, 3) without corresponding peaks in the MspI lane. Neither fragment 1 nor fragment 3 was resistant to combined HpaII + MspI cleavage (-, +, -), indicating that fragments 1 and 3 contained two consecutive differentially methylated CCGG sites: one proximal to the EcoRI site methylated at CmCGG and one distal non-methylated CCGG site. The products of MspI digestion are not visualized probably because they are too short and/or are not amplified by the selective primers. In hypomethylated DNA, we observed a loss of fragments 1 and 3 after HpaII digestion (Figure  3, panel H_h), further supporting the presence of internal CmCGG site(s) in the fragments. The only truly non-methylated fragment is fragment 2 that is visualized in the H, M and H + M lanes and indicates proper DNA digestion.

thumbnailFigure 2. Verification of the internal CmCGG site in (-, +) MH EcoRI-HpaII fragments by the addition of the M + H profile. Two different DNA molecules producing a (-, +) MH pattern for the long (L) fragment can contain either an internal CmCGG site (1) or hemimethylated mCCGG site (2). Combined digestion by both MspI and HpaII enzymes (M + H lanes) distinguish between the two molecules. Note that digested short (S) and long (L) DNA fragments may be amplified by different selective primers; if so, they are not visible together on one MSAP profile.

thumbnailFigure 3. Example of MSAP electropherograms documenting assessment of (-, +) MH patterns using MspI +HpaII double digestion. Fragments 1 and 3 represent the (-, +) MH pattern, fragment 2 represents the (+, +) MH pattern, and fragment 4 represents the (+, -) MH pattern. Note that fragments 1 and 3 disappear after application of the combination of the HpaII, MspI and EcoRI enzymes (M + H panel) or after DNA hypomethylation (H_h panel) by 100 μM DHPA [48]. Each restriction fragment is represented by two peaks (green and blue) derived from two complementary DNA strands (differently labeled selective primers were used). The presence of both peaks at the same position or in close proximity (note that the mobility of complementary strands of a DNA fragment may not be equal; [23]) indicates that the fragment had both EcoRI and MspI/HpaII ends.

Additional file 1: Figure S1. MSAP analysis of DNA samples isolated from tobacco seedlings treated with 0 μM (DHPA 0), 10 μM (DHPA 10) and 100 μM (DHPA 100) 9-(S)-(2,3-dihydroxypropyl)-adenine (DHPA).

Format: PDF Size: 367KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Dinucleotide "mCG only" methylation is localized in coding regions known to possess high GC content and a relatively high frequency of CCGG sites compared with the rest of the plant genome. Thus, the (-, +, -) patterns are likely formed by fragments having their origin in coding regions. By contrast, the (-, +, +) patterns, indicating mCCG methylation, would originate prevalently from repetitive non-coding regions because trinucleotide mCHG methylation is primarily localized to satellite repeats and transposons [46]. Thus, the proposed MSAP modification utilizing combined M + H digestion may help to determine the genomic origin of methylated fragments. The mCmCGG sites, which are common in plant heterochromatin, will not be digested with either MspI or HpaII (-, -, - pattern). This may explain why heterochromatic sequences are usually underrepresented in eluted MSAP fragments.

Regarding vertebrate DNA, which is methylated almost exclusively in CG sequences, we infer that the (-, +) MH pattern may not be interpreted as mCCGG methylation [21] but rather as internal CmCGG sites.

Because there is no entirely methylation-insensitive isoschizomer cutting at the CCGG motif, discrimination between methylation (epigenetic change) and mutation (genetic change) may be difficult, particularly among individuals of different genotypes. Therefore, employing other restriction sites may be an interesting option. Among others, Sau3AI (sensitive to methylation at GATmC sites), EcoT22I (ATGmCAT, [47]), NlaIII (mCATG), and ScrFI (CmCNGG) may be considered. For Sau3AI and NlaIII, methylation-insensitive MboI and CviAII izoschizomers are available, respectively.

Conclusions

Clearly, the MSAP method is a reliable, inexpensive, and relatively simple genome-wide method for the identification of genome regions with putative changes in DNA cytosine methylation in response to environmental and developmental stimuli. While the interpretation of the (+, -) and (+, +) MspI/HpaII patterns is straightforward, the presence of an HpaII signal without a concomitant MspI signal (-, +) is more ambiguous. In the latter case, we recommend modifying the MSAP protocol by including an additional lane in which the DNA sample is digested with a combination of both enzymes. Such improvement helps to distinguish between di- and trinucleotide methylation and between methylation in coding and noncoding areas of the plant genome.

Methods

DNA samples

We used DNA samples isolated from tobacco seedlings treated with 0, 10, and 100 μM 9-(S)-(2,3-dihydroxypropyl)-adenine (DHPA) and analyzed previously by MspI/HpaII digestion and Southern blot hybridization [48]. DHPA is a reversible competitive inhibitor of S-adenosylhomocysteine-hydrolase that elevates the intracellular SAH concentration [49]. SAH competes with S-adenosylmethionine (SAM), leading to inhibition of most SAM-dependent methyltransferases. In plants, DHPA treatment induces hypomethylation of CHG sites preferentially (CHG DNA methylation represents a metabolic bottleneck due to its dependence on histone methylation) and also hypomethylation of CG sites to some extent at elevated DHPA concentrations [48].

MSAP procedure

The MSAP procedure was performed in two technical replicates starting from the same DNA sample. DNA (1 μg) was digested in NEB buffer 1 with 10 U of HpaII (NEB) in a 20-μl reaction volume for 3 hours at 37°C. Ten microliters of the reaction mixture was transferred to 10 μl of NEB 2 buffer with 5 U of MspI (NEB) and incubated for 3 hours at 37°C. Simultaneously, 0.5 μg of the same original DNA sample was incubated in NEB 2 buffer with 5 U of MspI (NEB) in a 10-μl reaction volume for 3 hours at 37°C. Next, 10 μl of EcoRI buffer (NEB) with 5 U of EcoRI (NEB) were added to 10-μl reactions (HpaII and MspI digestions) and 5 μl of 2.5× EcoRI buffer (NEB) with 5 U of EcoRI (NEB) were added to a 20-μl HpaII + MspI reaction volume followed by incubation for an additional 1.5 hours at 37°C. EcoRI_ADAPTER1(F) 5′-CTCGTAGACTGCGTACC-3′ (10 μM)/ EcoRI_ADAPTER2(R) 5′-AATTGGTACGCAGTCTAC-3′ (10 μM) and HpaII/MspI_ADAPTER1(F) 5′-GATCATGAGTCCTGCT-3′ (100 μM)/ HpaII/MspI_ADAPTER2(R) 5′-CGAGCAGGACTCATGA-3′ (100 μM) were used to prepare 5 μM of EcoRI and 50 μM of MspI/HpaII adapters, respectively, by mixing appropriate oligonucleotides and incubating at 98°C for 5 min followed by slow cooling in a polystyrene box (for approximately 6 hours). Six microliters of digested DNA samples was added to 3.75 μl of mix comprising 1 μl of 5 μM EcoRI adapter (5 pmol), 1 μl of 50 μM MspI/HpaII adapter (50 pmol), 1 μl of fresh 10× T4 DNA ligase buffer, and 0.75 μl of sterile redistilled water. After mixing and heating to 37°C, 0.25 μl of T4 DNA ligase (100 U, NEB) was added, and the mix was incubated at 37°C for 3 hours followed by incubation at 20°C overnight. Non-ligated adapters were removed, and ligated DNA was purified using the NucleoSpin Gel and PCR Clean-up Kit (Macherey-Nagel). DNAs were eluted in 2 × 10 μl of 5 mM Tris-Cl (pH 8), and 3 μl was used in 10-μl preselective amplification reactions with EcoRI_A 5′-GACTGCGTACCAATTCA-3′ and HpaII/MspI_T 5′-ATGAGTCCTGCTCGGT-3′ primers at 200 nM concentration each. In addition, the preselective reactions contain DynazymeII buffer, 125 μM dNTPs, and 0.5 U of DynazymeII DNA Polymerase (Thermo Scientific). The conditions of the preselective amplification were as follows: 72°C for 5 min, 94°C for 3 min followed by 30 cycles of 94°C for 30 s, 56°C for 30 s, and 72°C for 2 min. A final step at 60°C for 10 min was also added. Five microliters of preselective products was checked by electrophoresis in 2% agarose gels (visible as a smear from 100 to 1000 bp) and the remaining 5 μl was diluted by adding 50 μl of redistilled water. The diluted preselective product (2.5 μl) was amplified using fluorescently labeled EcoRI_ACT[HEX] (green peaks in Figure  3 and Additional file 1: Figure S1) and labeled HpaII/MspI_TTC[6-FAM] (blue peaks in Figure  3) or HpaII/MspI_TAG (Additional file 1: Figure S1) selective primers (200 nM each) in a 12.5-μl total reaction volume containing KAPA buffer B, 125 μM dNTPs, and 0.5 U of KAPA Taq DNA Polymerase (KAPA Biosystems). The conditions of the selective amplification were as follows: 94°C for 3 min, 13 cycles of 94°C for 30 s, 65°C for 30 s reduced by 0.7°C per cycle, and 72°C for 2 min followed by 24 cycles of 94°C for 30 s, 56°C for 30 s, and 72°C for 2 min. A final step at 72°C for 10 min was also added. The product of selective amplification (5 μl) was checked by 2% agarose gel electrophoresis. Next, 1.5 μl of the selective amplification product was added to a mix of 10 μl of deionized formamide for capillary electrophoresis (Genomac) and 0.3 μl of GMC GT500-L DNA standard (Genomac) that is fluorescently labeled with LIZ dye (orange peaks in Figure  3 and Additional file 1: Figure S1). DNA was denaturated by heating at 98°C for 3 min followed by cooling on ice. Denatured DNA was separated using the fragment analysis option in the ABI Prism 3100 Genetic Analyzer. The obtained electropherograms were aligned according the signal of the DNA standard, analyzed and visualized using GeneMarker version 1.80 (SoftGenetics).

Authors’ contributions

JF and AK designed the study, JF performed the experiments. JF and AK analyzed and interpreted the data as well as drafted the manuscript. Both authors read and approved the final manuscript.

Acknowledgements

This work was supported by the Czech Science Foundation (P501/12/G090, 13-10057S, 206/09/1751) and the Academy of Sciences of the Czech Republic (RVO68081707). The authors thank Dr. Bohuslav Janoušek and Ing. Vojtěch Hudzieczek for providing the EcoRI adapters, preselective primer, and fluorescently labeled primer, as well as for the valuable suggestions concerning fragment analysis. The authors also thank Dr. Jiří Široký for technical advice concerning fragmentation analysis (ABI Prism Genetic Analyzer) and providing a trial version of GeneMarker software.

References

  1. Rai K, Jafri IF, Chidester S, James SR, Karpf AR, Cairns BR, Jones DA: Dnmt3 and G9a cooperate for tissue-specific development in zebrafish.

    J Biol Chem 2010, 285(6):4110-4121. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Renfree MB, Suzuki S, Kaneko-Ishino T: The origin and evolution of genomic imprinting and viviparity in mammals.

    Philos Trans R Soc Lond B Biol Sci 2013, 368(1609):20120151. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  3. Sharp AJ, Stathaki E, Migliavacca E, Brahmachary M, Montgomery SB, Dupre Y, Antonarakis SE: DNA methylation profiles of human active and inactive X chromosomes.

    Genome Res 2011, 21(10):1592-1600. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Yoder JA, Walsh CP, Bestor TH: Cytosine methylation and the ecology of intragenomic parasites.

    Trends Genet 1997, 13(8):335-340. PubMed Abstract | Publisher Full Text OpenURL

  5. Lim U, Song MA: Dietary and lifestyle factors of DNA methylation.

    Methods Mol Biol 2012, 863:359-376. PubMed Abstract | Publisher Full Text OpenURL

  6. Solis MT, Rodriguez-Serrano M, Meijon M, Canal MJ, Cifuentes A, Risueno MC, Testillano PS: DNA methylation dynamics and MET1a-like gene expression changes during stress-induced pollen reprogramming to embryogenesis.

    J Exp Bot 2012, 63(18):6431-6444. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  7. Durand S, Bouche N, Perez Strand E, Loudet O, Camilleri C: Rapid establishment of genetic incompatibility through natural epigenetic variation.

    Curr Biol 2012, 22(4):326-331. PubMed Abstract | Publisher Full Text OpenURL

  8. Law JA, Jacobsen SE: Establishing, maintaining and modifying DNA methylation patterns in plants and animals.

    Nature rev 2010, 11(3):204-220. Publisher Full Text OpenURL

  9. Saze H, Kakutani T: Differentiation of epigenetic modifications between transposons and genes.

    Curr Opin Plant Biol 2011, 14(1):81-87. PubMed Abstract | Publisher Full Text OpenURL

  10. Lunerova-Bedrichova J, Bleys A, Fojtova M, Khaitova L, Depicker A, Kovarik A: Trans-generation inheritance of methylation patterns in a tobacco transgene following a post-transcriptional silencing event.

    Plant J 2008, 54(6):1049-1062. PubMed Abstract | Publisher Full Text OpenURL

  11. Zemach A, McDaniel IE, Silva P, Zilberman D: Genome-wide evolutionary analysis of eukaryotic DNA methylation.

    Science 2010, 328(5980):916-919. PubMed Abstract | Publisher Full Text OpenURL

  12. Su Z, Han L, Zhao Z: Conservation and divergence of DNA methylation in eukaryotes: new insights from single base-resolution DNA methylomes.

    Epigenetics 2011, 6(2):134-140. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Ramsahoye BH, Biniszkiewicz D, Lyko F, Clark V, Bird AP, Jaenisch R: Non-CpG methylation is prevalent in embryonic stem cells and may be mediated by DNA methyltransferase 3a.

    Proc Natl Acad Sci USA 2000, 97(10):5237-5242. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Ichiyanagi T, Ichiyanagi K, Miyake M, Sasaki H: Accumulation and loss of asymmetric non-CpG methylation during male germ-cell development.

    Nucleic Acids Res 2013, 41(2):738-745. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Fulnecek J, Matyasek R, Kovarik A, Bezdek M: Mapping of 5-methylcytosine residues in Nicotiana tabacum 5S rRNA genes by genomic sequencing.

    Mol Gen Genet 1998, 259(2):133-141. PubMed Abstract | Publisher Full Text OpenURL

  16. Cokus SJ, Feng S, Zhang X, Chen Z, Merriman B, Haudenschild CD, Pradhan S, Nelson SF, Pellegrini M, Jacobsen SE: Shotgun bisulphite sequencing of the Arabidopsis genome reveals DNA methylation patterning.

    Nature 2008, 452(7184):215-219. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Lister R, Ecker JR: Finding the fifth base: Genome-wide sequencing of cytosine methylation.

    Genome Res 2009, 19(6):959-966. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Down TA, Rakyan VK, Turner DJ, Flicek P, Li H, Kulesha E, Graf S, Johnson N, Herrero J, Tomazou EM, et al.: A Bayesian deconvolution strategy for immunoprecipitation-based DNA methylome analysis.

    Nat Biotechnol 2008, 26(7):779-785. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Paun O, Bateman RM, Fay MF, Hedren M, Civeyrel L, Chase MW: Stable epigenetic effects impact adaptation in allopolyploid orchids (Dactylorhiza: Orchidaceae).

    Mol Biol Evol 2010, 27(11):2465-2473. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Xu Y, Zhao Q, Mei S, Wang J: Genomic and transcriptomic alterations following hybridisation and genome doubling in trigenomic allohexaploid Brassica carinata x Brassica rapa.

    Plant biol (Stuttg) 2012, 14(5):734-744. Publisher Full Text OpenURL

  21. Yang C, Zhang M, Niu W, Yang R, Zhang Y, Qiu Z, Sun B, Zhao Z: Analysis of DNA methylation in various swine tissues.

    PLoS One 2011, 6(1):e16229. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Reyna-Lopez GE, Simpson J, Ruiz-Herrera J: Differences in DNA methylation patterns are detectable during the dimorphic transition of fungi by amplification of restriction polymorphisms.

    Mol Gen Genet 1997, 253(6):703-710. PubMed Abstract | Publisher Full Text OpenURL

  23. Vos P, Hogers R, Bleeker M, Reijans M, van de Lee T, Hornes M, Frijters A, Pot J, Peleman J, Kuiper M, et al.: AFLP: a new technique for DNA fingerprinting.

    Nucleic Acids Res 1995, 23(21):4407-4414. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Walder RY, Langtimm CJ, Chatterjee R, Walder JA: Cloning of the MspI modification enzyme. The site of modification and its effects on cleavage by MspI and HpaII.

    J Biol Chem 1983, 258(2):1235-1241. PubMed Abstract | Publisher Full Text OpenURL

  25. Mann MB, Smith HO: Specificity of Hpa II and Hae III DNA methylases.

    Nucleic Acids Res 1977, 4(12):4211-4221. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Roberts RJ, Vincze T, Posfai J, Macelis D: REBASE--a database for DNA restriction and modification: enzymes, genes and genomes.

    Nucleic Acids Res 2010, 38(Database issue):D234-D236. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Lin PM, Lee CH, Roberts RJ: Cloning and characterization of the genes encoding the MspI restriction modification system.

    Nucleic Acids Res 1989, 17(8):3001-3011. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Card CO, Wilson GG, Weule K, Hasapes J, Kiss A, Roberts RJ: Cloning and characterization of the HpaII methylase gene.

    Nucleic Acids Res 1990, 18(6):1377-1383. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Waalwijk C, Flavell RA: MspI, an isoschizomer of hpaII which cleaves both unmethylated and methylated hpaII sites.

    Nucleic Acids Res 1978, 5(9):3231-3236. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Ben-Hattar J, Jiricny J: Effect of cytosine methylation on the cleavage of oligonucleotide duplexes with restriction endonucleases HpaII and MspI.

    Nucleic Acids Res 1988, 16(9):4160. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Tardy-Planechaud S, Fujimoto J, Lin SS, Sowers LC: Solid phase synthesis and restriction endonuclease cleavage of oligodeoxynucleotides containing 5-(hydroxymethyl)-cytosine.

    Nucleic Acids Res 1997, 25(3):553-559. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Korch C, Hagblom P: In-vivo-modified gonococcal plasmid pJD1. A model system for analysis of restriction enzyme sensitivity to DNA modifications.

    Eur J Biochem 1986, 161(3):519-524. PubMed Abstract | Publisher Full Text OpenURL

  33. Butkus V, Petrauskiene L, Maneliene Z, Klimasauskas S, Laucys V, Janulaitis A: Cleavage of methylated CCCGGG sequences containing either N4-methylcytosine or 5-methylcytosine with MspI, HpaII, SmaI, XmaI and Cfr9I restriction endonucleases.

    Nucleic Acids Res 1987, 15(17):7091-7102. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Huang LH, Farnet CM, Ehrlich KC, Ehrlich M: Digestion of highly modified bacteriophage DNA by restriction endonucleases.

    Nucleic Acids Res 1982, 10(5):1579-1591. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Brennan CA, Van Cleve MD, Gumport RI: The effects of base analogue substitutions on the cleavage by the EcoRI restriction endonuclease of octadeoxyribonucleotides containing modified EcoRI recognition sequences.

    J Biol Chem 1986, 261(16):7270-7278. PubMed Abstract | Publisher Full Text OpenURL

  36. Jeddeloh JA, Richards EJ: mCCG methylation in angiosperms.

    Plant J 1996, 9(5):579-586. PubMed Abstract | Publisher Full Text OpenURL

  37. Fulnecek J, Matyasek R, Kovarik A: Distribution of 5-methylcytosine residues in 5S rRNA genes in Arabidopsis thaliana and Secale cereale.

    Mol Gen Genomics 2002, 268(4):510-517. Publisher Full Text OpenURL

  38. Meyer P: DNA methylation systems and targets in plants.

    FEBS Lett 2011, 585(13):2008-2015. PubMed Abstract | Publisher Full Text OpenURL

  39. Mathieu O, Reinders J, Caikovski M, Smathajitt C, Paszkowski J: Transgenerational stability of the Arabidopsis epigenome is coordinated by CG methylation.

    Cell 2007, 130(5):851-862. PubMed Abstract | Publisher Full Text OpenURL

  40. Sasaki T, Fujimoto R, Kishitani S, Nishio T: Analysis of target sequences of DDM1s in Brassica rapa by MSAP.

    Plant Cell Rep 2011, 30(1):81-88. PubMed Abstract | Publisher Full Text OpenURL

  41. Xu M, Li X, Korban SS: DNA-methylation alterations and exchanges during in vitro cellular differentiation in rose (Rosa hybrida L.).

    Theor Appl Genet 2004, 109(5):899-910. PubMed Abstract | Publisher Full Text OpenURL

  42. Takamiya T, Hosobuchi S, Asai K, Nakamura E, Tomioka K, Kawase M, Kakutani T, Paterson AH, Murakami Y, Okuizumi H: Restriction landmark genome scanning method using isoschizomers (MspI/HpaII) for DNA methylation analysis.

    Electrophoresis 2006, 27(14):2846-2856. PubMed Abstract | Publisher Full Text OpenURL

  43. Takamiya T, Hosobuchi S, Noguchi T, Asai K, Nakamura E, Habu Y, Paterson AH, Iijima H, Murakami Y, Okuizumi H: Inheritance and alteration of genome methylation in F1 hybrid rice.

    Electrophoresis 2008, 29(19):4088-4095. PubMed Abstract | Publisher Full Text OpenURL

  44. Koroma AP, Jones R, Michalak P: Snapshot of DNA methylation changes associated with hybridization in Xenopus.

    Physiol genomics 2011, 43(22):1276-1280. PubMed Abstract | Publisher Full Text OpenURL

  45. Xiao J, Song C, Liu S, Tao M, Hu J, Wang J, Liu W, Zeng M, Liu Y: DNA methylation analysis of allotetraploid hybrids of red crucian carp (Carassius auratus red var.) and common carp (Cyprinus carpio L.).

    PLoS One 2013, 8(2):e56409. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  46. Zakrzewski F, Weisshaar B, Fuchs J, Bannack E, Minoche AE, Dohm JC, Himmelbauer H, Schmidt T: Epigenetic profiling of heterochromatic satellite DNA.

    Chromosoma 2011, 120(4):409-422. PubMed Abstract | Publisher Full Text OpenURL

  47. You W, Tyczewska A, Spencer M, Daxinger L, Schmid MW, Grossniklaus U, Simon SA, Meyers BC, Matzke AJ, Matzke M: Atypical DNA methylation of genes encoding cysteine-rich peptides in Arabidopsis thaliana.

    BMC Plant Biol 2012, 12:51. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  48. Fulnecek J, Matyasek R, Votruba I, Holy A, Krizova K, Kovarik A: Inhibition of SAH-hydrolase activity during seed germination leads to deregulation of flowering genes and altered flower morphology in tobacco.

    Mol Genet Genomics 2011, 285(3):225-236. PubMed Abstract | Publisher Full Text OpenURL

  49. Fojtova M, Kovarik A, Votruba I, Holy A: Evaluation of the impact of S-adenosylhomocysteine metabolic pools on cytosine methylation of the tobacco genome.

    Eur J Biochem 1998, 252(3):347-352. PubMed Abstract | Publisher Full Text OpenURL