Skip to main content

A relative quantitative Methylation-Sensitive Amplified Polymorphism (MSAP) method for the analysis of abiotic stress



We present a new methylation-sensitive amplified polymorphism (MSAP) approach for the evaluation of relative quantitative characteristics such as demethylation, de novo methylation, and preservation of methylation status of CCGG sequences, which are recognized by the isoschizomers HpaII and MspI. We applied the technique to analyze aluminum (Al)-tolerant and non-tolerant control and Al-stressed inbred triticale lines. The approach is based on detailed analysis of events affecting HpaII and MspI restriction sites in control and stressed samples, and takes advantage of molecular marker profiles generated by EcoRI/HpaII and EcoRI/MspI MSAP platforms.


Five Al-tolerant and five non-tolerant triticale lines were exposed to aluminum stress using the physiologicaltest. Total genomic DNA was isolated from root tips of all tolerant and non-tolerant lines before and after Al stress following metAFLP and MSAP approaches. Based on codes reflecting events affecting cytosines within a given restriction site recognized by HpaII and MspI in control and stressed samples demethylation (DM), de novo methylation (DNM), preservation of methylated sites (MSP), and preservation of nonmethylatedsites (NMSP) were evaluated. MSAP profiles were used for Agglomerative hierarchicalclustering (AHC) based on Squared Euclidean distance and Ward’s Agglomeration method whereas MSAP characteristics for ANOVA.


Relative quantitative MSAP analysis revealed that both Al-tolerant and non-tolerant triticale lines subjected to Al stress underwent demethylation, with demethylation of CG predominating over CHG. The rate of de novo methylation in the CG context was ~3-fold lower than demethylation, whereas de novo methylation of CHG was observed only in Al-tolerant lines.


Our relative quantitative MSAP approach, based on methylation events affecting cytosines within HpaII–MspI recognition sequences, was capable of quantifying de novo methylation, demethylation, methylation, and non-methylated status in control and stressed Al-tolerant and non-tolerant triticale inbred lines. The method could also be used to analyze methylation events affecting CG and CHG contexts, which were differentially methylated under Al stress. We cannot exclude that the methylation changes revealed among lines as well as between Al-tolerant and non-tolerant groups of lines were due to some experimental errors or that the number of lines was too small for ANOVA to prove the influence of Al stress. Nevertheless, we suspect that Al tolerance in triticale could be partly regulated by epigenetic factors acting at the level of DNA methylation. This method provides a valuable tool for studies of abiotic stresses in plants.


When plants are stressed, their DNA methylation patterns may change, as indicated by numerous morphological [13], biochemical [4, 5], and molecular [6, 7] studies. Alterations of epigenetic states play essential roles in protecting organisms from environmental stresses [7] and are believed to be involved in plant immunity [8]. Abiotic stresses such as drought [9], cold [10], salt [11], or heavy metals in the soil [12] can influence methylation patterns. Consistent with this, DNA methylation in Poa annua populations [13] is influenced by climatic factors, and similar evidence was obtained for Deschampsia antarctica [14]. Even short-term stresses affect DNA methylation, e.g., in barley regenerants derived via in vitro tissue culture from anthers [15]. In both barley [15] and Gentiana pannonica [16], global DNA methylation in regenerants increased relative to donor levels; in triticale, the DNA methylation of regenerants and their progeny decreased in comparison to those of donor plants and did not return to the initial level even after several reproductive cycles [17]. Furthermore, DNA methylation changes induced by abiotic stresses can be stably passed to progeny [18].

In plants, DNA methylation on cytosine occurs in all sequence contexts: symmetric CG and CHG (where H equals either A, C, or T) and asymmetric CHH. Both CG and CHG methylation can be copied during DNA replication, whereas CHH methylation must be reestablished in each generation via mechanisms involving DRM1, DRM2, and iRNAs [8]. DNA methylation patterns may change in response to stressful conditions [6], and such changes in symmetric CG sites may be involved in the regulation of gene expression [19, 20]. By contrast, methylation of CHG and CHH is linked to regulation of transposon activity via chromatin remodeling [2123] in response to environmental factors.

Many approaches are available for quantitative analysis of global genome DNA methylation. Probably the simplest and the most robust of these is based on RP-HPLC [24]. Recently, the methylation-sensitive amplified fragment length polymorphism (metAFLP) method was introduced [15] and then extended [25] to enable quantification of changes in sequence and DNA methylation sites (i.e., de novo methylation and demethylation) in a single experiment. Another option is methylation-sensitive amplified polymorphism (MSAP), which is based on AFLP technology and takes advantage of two isoschizomers that recognize the same restriction site [26]. This approach was adapted to the analysis of DNA methylation pattern in plants [27] and successfully used to study changes in cytosine methylation under abiotic stresses [11, 2830]. Using this technique, changes in DNA methylation sites could be evaluated by comparing molecular patterns [9]. Consensus has not been reached regarding how best to interpret MSAP outputs [31, 32], and at least five alternatives have been proposed [33]. However, none of the options takes into account the multiple events (reflecting various methylation pattern states or changes) that must take place simultaneously to explain the differences in individual digestion patterns between control and stressed materials. In addition, because mutations are hardly difficult to identify using the MSAP approach, they are excluded from the underlying models [33]. Further confusing interpretation of the results, no systematic studies have investigated the dependence of HpaII and MspI activities on the cytosine methylation state at all dCs in dsDNA [31]. Thus, a detailed analysis of the effects of methylation changes at restriction sites recognized by HpaII and MspI endonucleases would allow the potential of MSAP to be more fully realized. Such an analysis would need to incorporate an assessment of the background changes at restriction sites caused by abiotic stresses.

The aim of this study was to develop a theoretical model for quantification of cytosine methylation status at restriction sites of the isoschizomers HpaII and MspI, using only the documented activities of these enzymes, based on the MSAP profiles of triticale plants grown under Al stress.


Physiological tests revealed that root regrowth of five Al-tolerant triticale lines (T1, T2, T3, T4, and T5) ranged from 0.8 to 2.8 cm in the presence of Al3+ ions, whereas no regrowth was observed in five non-tolerant lines (NT1, NT2, NT3, NT4, and NT5) (Fig. 1).

Fig. 1
figure 1

Aluminum tolerance test according to Anioł A [54]. Roots were examined 48 h after staining with Eriochrome Cyanine R

Analysis of documented methylation patterns of the restriction site recognized by HpaII and MspI showed that HpaII can cut at sites that are either non-methylated or contain one methylated external cytosine, whereas MspI cuts non-methylated sites and those with one or two methylated internal cytosines [26]. Comparison of the digestion patterns of these isoschizomers could be used to evaluate the extent of restriction site DNA methylation. The original MSAP approach can be used to analyze abiotic stresses when identical samples are analyzed under different conditions (i.e., control and stressed) [11, 28, 29]. In such cases, the presence or absence of molecular markers generated by EcoRI/HpaII and EcoRI/MspI digestion in control and stressed samples can be represented by 16 four-bit codes. Each code, reflecting a specific banding pattern, has its own rationalization, as depicted in Table 1. Four molecular patterns could be easily explained by comparison of control and stressed restriction sites: code (1111) reflects the situation in which all non-methylated cytosines of a given restriction site present in controls remain non-methylated in stressed samples; (1110) indicates de novo methylation of the external cytosine in stressed samples; (1010) reflects preservation of the external mC and lack of methylation of the other C; and finally (1011) indicates demethylation of the external C and preservation of non-methylated status of the remaining Cs in stressed samples. Similarly, code (0000) pertains to the situation in which a site that is fully methylated in controls is not altered in the stressed samples. However, such a pattern cannot be easily recognized unless several replicates of stressed samples are used, and the pattern is identified in at least one of them. The other codes have more complicated explanations, and could result from differences in methylation patterns of a given restriction site between control and stressed samples (Table 1). Nevertheless, all 16 codes can be scored (Table 2). Assuming that all of the events (de novo methylation (DNM), demethylation (DM), and preservation of non-methylated (NMSP) and methylated (MSP) sites) are equally probable, it is possible to estimate the number of events of a given type, corresponding to each four-bit code, considering only the known properties of the isoschizomers. By multiplying the number of individual events participating in the explanation of a given code by the number of MSAP profiles depicted by that code, followed by summation of events of the same kind and normalization of the data (expressed as percentages), it is possible to evaluate the relative quantitative characteristics of de novo methylation, demethylation, and preservation of either methylated and non-methylated cytosines at a restriction site. Additionally, one can also determine the total number of methylated and non-methylated cytosines in stressed samples (Tables 1 and 2). Because both external and internal cytosine methylation states are being analyzed, similar reasoning enables quantification of de novo methylation and demethylation events, reflecting the CHG and CG symmetric methylation states of the restriction sites recognized by HpaII and MspI.

Table 1 Theoretical background used for the explanation of de novo methylation (DNM), demethylation (DM), preservation of non-methylated sites (NMSP) and preservation of methylated sites (MSP) events based on known properties of HpaII and MspI isoschizomers
Table 2 The arrangement of all possible events explaining the sixteen four-bit binary MSAP codes

Quantifying methylation events in Al-stressed triticale inbred lines

Of the 16 available MSAP profiles, we identified only seven in ten triticale lines that were maintained under control and Al-stressed conditions (Table 3). The most frequent were the cases corresponding to codes 1010, 0101, and 1111, whereas code 1100 was detected only twice. The patterns of the seven codes were more or less uniform among the samples analyzed.

Table 3 The arrangement of the profiles reflecting given MSAP four-bit binary code evaluated among control and Al stressed triticale inbred lines

When the codes were converted into events of the corresponding type (Table 4), and then expressed as normalized relative quantitative characteristics (Table 5), we observed slight differences among the analyzed lines. When stressed, two tolerant lines (T1 and T2) increased the percent de novo methylation (DNM%) of some cytosines (average, 4.19%) compared with that of other cytosines (average, 4.04%). The demethylation level of one tolerant line (T5) (5.89%) was lower than the average score (6.28 ± 0.25). By contrast, two non-tolerant lines (NT2 and NT5) exhibited reduced DM% relative to the three others (NT1, NT3, and NT4), whereas DNM% was comparable (average 4.04%) in all five lines.

Table 4 The arrangement of events recalculated based on data presented in Tables 3 and 2
Table 5 Quantitative characteristics evaluated based on MSAP data expressed in percentages

The other characteristics of the analyzed materials were comparable. Detailed analysis of DM-CG% (percentages of demethylated CG sequences within restriction sites recognized by HpaII and MspI), DM-CHG%, DNM-CG%, and DNM-CHG% revealed that DM-CHG% and DNM-CHG% were lower than their CG counterparts. Moreover, DNM-CG% was slightly lower than DM-CG%, whereas DNM-CHG% was always many times lower than DM-CHG%. It should be emphasized, however, that DNM-CHG% was determined exclusively in two tolerant lines (NT1 and NT2).

Agglomerative hierarchical clustering performed on the data presented in Table 5 demonstrated that NT1, NT3, NT4, T1, T2, T3, and T4 formed one cluster, with minor differences between T1 and T2 as well as among NT1, NT3, and T3, whereas NT2, NT5, and T5 clustered separately (Fig. 2).

Fig. 2
figure 2

Agglomerative hierarchical clustering of T and NT lines. The dotted line represents the automatic truncation allowing discrimination among groups of lines

Analysis of variance based on all MSAP characteristics together showed no differences among the ten lines, demonstrating that the 1% variability of the dependent variable could be explained by the explanatory variable. Similarly, ANOVA failed to detect differences between tolerant and non-tolerant lines. When individual relative quantitative MSAP characteristics were considered, they explained from 0 to 100% (in the case of DNM%) of the variance in all lines, but the differences were not significant. ANOVA performed on DNM% data between tolerant and non-tolerant lines explained 10% of the variability of the dependent variable and was negligible. When DN-CG%, DN-CHG%, DM-CG%, and DM-CHG% were used as MSAP characteristics, ANOVA also failed to identify differences among lines as well as between Al tolerant and non-tolerant groups of lines. However, up to 39% of the variability among all lines was due to differences in DN-CHG% and DM-CHG% data whereas up to 25% of the differences between T and NT lines could be explained by differences in DN-CHG%. Analysis of variance demonstrated that the differences among DN%, DM%, MP%, NMP%, M% and NM% were significant (F = 5.706, p-value <0.0001). Similar data were evaluated for DN-CG%, DN-CHG%, DM-CG% and DM-CHG% characteristics (F = 12.84, p-value <0.0001).


AFLP technologies (e.g., MSAP and metAFLP), which use isoschizomers with different specificities for DNA methylation at the restriction site, have been extensively exploited to analyze variations in DNA methylation patterns in plants under stress [9, 12, 15, 34]. MSAP was introduced by Reyna-Lopez et al. [26] in studies on fungi and adapted for used in plants by Xiong et al. [27], whereas metAFLP was developed for analyses of variation induced during in vitro tissue culture plant regeneration [15]. MetAFLP, which is focused on quantification of ‘dynamic’ changes that can be assessed between control and stressed materials, can accommodate both sequence and methylation changes at the same time. By contrast, MSAP was developed to study site-specific DNA methylation status. By default, in MSAP, mutations within restriction sites are hardly difficult to identify because both endonucleases are methylation sensitive [31]. Thus, such mutations are not considered in the model being assumed to be rare and not affecting the results. This assumption is supported by experiments showing that methylation at mCmCG sequence motifs is frequent [35], whereas sequence changes are comparatively rare [36]. However, data obtained using metAFLP in triticale [18, 25] revealed that sequence changes in in vitro regenerated plants were not rare, in some cases exceeding the frequency of alterations in methylation patterns. On the other hand, similar data in barley [15] confirmed that the rate of sequence changes was significantly lower than that of changes in DNA methylation. Therefore, further studies are required to conclusively determine whether MSAP analysis is influenced by sequence changes, as well as whether such changes are species-specific. In this study, however, mutations were not incorporated into the model, following the reasoning presented by the others [29, 3436].

In addition to the unsolved problem of point mutations, which may influence the MSAP results, it remains unclear how best to interpret the molecular profiles generated by the approach; currently, at least five alternatives exist [37]. All of them are strongly correlated, indicating that they reflect similar phenomena [32, 37]. However, the proposed interpretations do not take into account the fact that all MSAP molecular profiles reflect multiple linked events responsible for changes to, or preservation of, various methylation patterns. For example, some profiles could be explained by de novo methylation of one cytosine and demethylation of another, while the methylation status of the remaining cytosines is preserved. Moreover, the same molecular profiles could be explained by multiple types of changes in the methylation pattern. Nevertheless, by assuming that the documented activities of HpaII and MspI are accurate, and by analyzing the events underlying the molecular profiles generated by the MSAP approach from control and stressed samples, it is possible to quantify events that were previously invisible to established methods for MSAP data analysis. Given that our proposed method is based on an equal probability of events affecting restriction sites rather than numbers of analyzed bands, it is capable of evaluating relative quantitative characteristics such as DNM%, DM%, MSP%, and NMSP%. In that sense, our method is similar to metAFLP, where the same strategy is applied [15, 25]. It should be emphasized, however, that the method only allows strictly defined types of restriction site methylation patterns, described in the literature as recognized by HpaII and MspI. This assumption is likely valid because certain methylation patterns are extremely rare (5′-CmCGG-3′/5′-GGCmC-3′; 5′-CCGG-3′/5′-GGCC-3′), whereas others (i.e. 5′-mCCGG-3′/5′-GGCmC-3′; 5′-mCCGG-3′/5′-GGCmCm-3′; 5′-CmCGG-3′/5′-GGCmCm-3′), to the best of our knowledge, have never been described [31]. Thus, it would not be reasonable to incorporate them into the model. Nevertheless, we cannot exclude the possibility that, in some cases, even such rare patterns may appear. Studies of the properties of HpaII and MspI, including all putative methylation patterns of the restriction site, could help resolve this issue.

As stressed earlier, the proposed method assumes that there is an equal probability of events affecting restriction sites. In our earlier studies using the metAFLP approach [15, 25], we demonstrated, however, that de novo demethylation and sequence changes varied in barley [15] and triticale [25]. Thus, the equal probability assumption should be validated e.g., by a whole-genome bisulfite approach [38, 39] conducted on a given genome and restricted to restriction sites recognized by HpaII and MspI isoschizomers. However, despite the availability of triticale sequencing data [40], this data cannot be adopted for evaluation in a more advanced MSAP model because information about methylation changes has not been evaluated yet. If available, the information could be easily incorporated into the relative quantitative MSAP approach by using relevant probability factors that most probably will be species specific.

Moreover, the relative quantification model assumes that DNA digestion is complete and that co-migrating markers can be neglected. Thus, information obtained from a quantitative treatment of the marker band intensities) could be lost. Nevertheless, quantification of band intensities could be incorporated into the model assuming that the denser bands originate from the same event types. Quantification could be accomplished by multiplication of the events represented by more intense bands by the degree to which their intensity exceeds the average level of the given marker. However, without sequencing data, it will be difficult to be sure that putative differences in band intensities are not related to co-migrating markers of distinct origin and are not necessarily related to the same methylation change as the weaker markers. Assuming that differences in band density are relatively rare and the fact they may not necessary reflect identical events quantification of such cases was omitted by us.

The proposed MSAP approach allowed us to evaluate quantitative characteristic of ten triticale inbred lines that differed in regard to Al tolerance. Al stress resulted in elevated demethylation in both non-tolerant and tolerant plants, with de novo methylation occurring less frequently than demethylation. Total demethylation due to Al stress was about 2% in both cases. These observations are in close agreement with data obtained from salt-tolerant and non-tolerant forms of rice [37], consistent with the notion that abiotic stress results in demethylation and de novo methylation in sensitive and tolerant variants, respectively. However, this is not always the case: the opposite changes were detected in sorghum under Al stress [41], indicating that the trend could depend on the species or the specific type of stress.

The advantage of the relative quantitative MSAP approach is that it can evaluate subtle effects related to methylation events affecting CG and CHG contexts, using reasoning similar to that applied to the general characteristics. Such sequences play crucial roles in plant genomes [8, 42] because they contribute to the symmetric methylation of DNA. Symmetric and asymmetric methylation may predominate in different regions of the genome, such as genes and transposable elements (TEs), as well as in their up and downstream regions [43]. Moreover, symmetric methylation is associated with either fixed or recent changes in methylation patterns [7, 44]. Changes in DNA methylation affecting CG, CHG, and CHH contexts are involved in adaptive responses to abiotic stresses [43, 45, 46]. Such alterations could be of value when determining stress levels in analyzed materials [19, 46, 47]. Therefore, the development of a relatively inexpensive marker-based approach would be a valuable alternative to whole-genome bisulfite sequencing [4850].

Although the results of relative quantitative MSAP, which reveals differences in methylation status of restriction sites between control and stressed samples, are generally in good agreement with metAFLP and HPLC data [17], they do not reveal what happens in the CG and CHG contexts. In the triticale lines we studied, CG sequences underwent most of the changes in methylation status in response to Al stress, with DM% slightly exceeding DNM%. Similarly, DM% was higher than DNM% in CHG, but the overall level of changes in this context was usually 3-fold lower than in CG. In addition, in both tolerant and non-tolerant lines, demethylation of CG and CHG exceeded de novo methylation. Only in two Al-tolerant lines, we observed an increase in de novo methylation of CHG sequences. In general, our data are congruent with those for tobacco [19], in which Al stress results in demethylation of CG sequences within coding regions. However, we observed no effects in the CHG or CHH context, regardless of the stress [19]. Demethylation of CG and CHG sequences was also observed in cotton under alkali stress [28]. Similar demethylation of symmetric and asymmetric sequences within promoter regions was observed in wheat [46], whereas in Chloris virgata, de novo CHG methylation under salt and alkaline stresses occurred in roots without an accompanying increase in CG methylation [45]. On the other hand, in sorghum, Al stress induces de novo methylation of CG or CHG sequences, depending on the line [41]. Thus, it is possible for methylation changes to affect symmetric and asymmetric contexts in different directions, and detailed studies are required to understand this phenomenon and its role in abiotic stresses.

The results of relative quantitative MSAP revealed that the analyzed triticale lines could differ from each other with respect to DNA methylation patterns induced by Al stress. Two of the non-tolerant lines and one tolerant line were distinct from the other samples (see also dendrogram Fig. 2), and in particular exhibited lower levels of demethylation. It has been suggested that changes in methylation pattern are to some extent stochastic [47, 51, 52]. On the other hand, the high degree of similarity among the Al-tolerant lines suggest that in any given case methylation could be directed towards specific genomic regions. This notion is supported by the observed increase in CHG de novo methylation in Al-tolerant lines. However, these changes were rather small and limited to a few lines. Moreover, hierarchical clustering used for the visualization of putative differences among triticale lines could be the result of chance, which could lead to misinterpretation of the data. Putative misinterpretation of the results of cluster analysis seems to have been confirmed by ANOVA, which failed to reveal differences between Al tolerant and non-tolerant triticale lines based on any relative quantitative MSAP characteristic. Alternatively, analysis of variance may suggest that the number of lines was too small to discriminate between the analyzed materials at the level of methylation changes affecting CCGG restriction sites recognized by the HpaII and MspI endonucleases used in the MSAP approach. Thus, it remains unclear whether the methylation differences reflect biologically relevant phenomenon not univocally supported by statistics or experimental errors. Interestingly, however, ANOVA demonstrated significant differences among all MSAP characteristics, suggesting that they may reflect biologically important phenomenon. Evidently, further studies are needed to determine whether effects on DNA methylation patterns as subtle as those observed for CHG could be associated with the Al tolerance. On the other hand, known Al-tolerance QTLs explain up to 30% of the phenotypic variance in triticale [53], whereas the remaining 70% is unexplained. Thus, our data may suggest that epigenetic factors may also participate in Al tolerance in triticale.

Finally, it is good experimental practice to validate any approach using alternative methods [32, 39]. The results of the relative quantitative MSAP approach could be validated by comparing them to those obtained by RP-HPLC [32]. This type of analysis was performed previously by us in the case of the metAFLP analysis [17, 18], and some correlation between methylation parameters was detected by the two methods. In general, the metAFLP approach underestimated global changes. The incongruence was explained by the fact that metAFLP analyzes site DNA methylation changes affecting restriction sequences recognized by Acc65I and KpnI isoschizomers whereas the RP-HPLC approach detects global methylation changes. Moreover, RP-HPLC is robust whereas characteristics evaluated in the case of Al tolerant and non-tolerant lines were relatively small. Thus, RP-HPLC seems to be inappropriate for the validation of most MSAP characteristics (at least in the case of our experimental data) because it delivers only limited and averaged information on DNA methylation changes. Alternatively, MSAP results could be validated by whole-genome bisulfite sequencing data [38]. Such an analysis would be beneficial for the evaluation of additional factors that could be incorporated into the MSAP approach, and would allow testing of the assumption of equal probability of methylation changes assumed in the present study. Future studies should allow validation of the MSAP approach as well as its further development via the incorporation of whole-genome sequencing data. Thus, in the case of preliminary studies where e.g., abiotic factors are suspected of affecting stressed plants the relative quantitative MSAP approach could be the method of choice.


Our proposed MSAP approach is capable of quantifying events such as de novo methylation, demethylation, methylation, and preservation of methylated status within HpaII–MspI restriction sites. Moreover, it can quantify events affecting CG and CHG sequences, providing important insight into the role of epigenetics of abiotic stresses in plants.

The quantitative results clearly demonstrated that in triticale, Al stress resulted in demethylation and de novo methylation in both tolerant and non-tolerant triticale lines; de novo methylation of CHG sequence was affected in tolerant lines but not in non-tolerant lines. MSAP analysis was partially successful in differentiating the triticale lines. However, ANOVA showed that differences among lines as well as differences between Al tolerant and non-tolerant groups of lines based on global and individual MSAP characteristics (including CHG% and CG%, DN%, and DM%) were insignificant. This result could have been due to experimental error or to the fact that methylation changes between Al tolerant and non-tolerant lines evaluated by MSAP are not necessarily related to the abiotic stress analyzed in this study. Alternatively, the number of lines used was insufficient to perform statistical analysis on the relatively weak differences revealed by this MSAP method.


Plant material and growing conditions

Five Al-tolerant (T1, T2, T3, T4, and T5) and five non-tolerant (NT1, NT2, NT3, NT4, and NT5) triticale lines were exposed to aluminum stress using the physiological test described previously by Anioł [54]. Seeds were surface-sterilized by soaking in 10% sodium hypochlorite for 10 min, followed by three washes with deionized water and germination on moist filter paper in Petri dishes for 24 h. Seeds representative of each triticale line were sown on two individual polyethylene nets floated in a tray filled with basic medium containing 2.0 mM CaCl2, 3.25 mM KNO3, 1.25 mM MgCl2, 0.5 mM (NH4)2SO4, and 0.2 mM NH4NO3 (pH 4.5). Seedlings were maintained in a controlled-environment growth cabinet (POL-EKO-APARATURA, ST500 B40 FOT10) at 25 °C, with a photoperiod of 12 h/12 h day/night, light intensity of 40 W m−2, and aeration. After 3 days, one of two nets representing each line was transferred for 24 h onto the same medium supplemented with 20 ppm AlCl3. Then, stressed seedlings were washed for 2–3 min in running water and transferred to nutrient solution without Al for 48 h. Finally, 7-day old seedlings were removed from control and stressful nutrient solutions and washed. Root tips (3–4 mm long) of control and stressed plants were cut off and kept at −70 °C prior to DNA isolation.

DNA isolation

Total genomic DNA was isolated from root tips (3–4 mm long) using the Plant DNeasy MiniKit (Qiagen), followed by spectrophotometric quantification and integrity testing by 1.0% agarose gel electrophoresis.

Methylation-sensitive amplification polymorphism (MSAP) profiling

The MSAP procedure is based on the metAFLP method [15]. However, instead of the Acc65I and KpnI endonucleases, HpaII and MspI were used according to Xiong et al. [27]. Briefly, genomic DNA samples were digested with HpaII/EcoRI or MspI/EcoRI (New England Biolabs, Beverly, MA, USA). Digestion mixtures contained 5.0 U HpaII or MspI, 5.0 U EcoRI, 1× ligation buffer, 50 nM NaCl, and 0.5 mg/ml BSA. Digestion, as well as adapter ligation and pre-selective and selective PCR, were performed as described for metAFLP [15]. The sequences of EcoRI and HpaII/MspI adapters and pre-selective and selective primers were as described in Xiong et al. [27]. In order to reduce the number of fragments to be amplified by each selective primer, the primer contained selective nucleotides at its 3′-terminal end. Fourteen primer combinations with three selective nucleotides for the EcoRI ends and three or four selective nucleotides for the HpaII-MspI ends were used (HpaII-MspI/EcoRI: TTG/ACG, TAC/AAG, TTC/AAC, TGA/AGG, TCG/ACT, TGC/AGT, TGT/ATT, TCA/ATC, TCG/AGA, TGC/ACA, TTGC/AGC, TCAA/ACC, TCAA/ATG, TCAA/ATT). Each HpaII–MspI selective primer was end-labeled with 10–30 μCi γ-[32P] ATP using T4 polynucleotide kinase prior to performing PCR. The denatured PCR products were separated on 7% denaturing polyacrylamide gels and exposed to X-ray film overnight at −70 °C. MSAP analysis was performed twice for each primer combination.

Methylation pattern analysis


MSAP markers were evaluated for all triticale inbred lines (including controls and their Al-stressed counterparts). MSAP profiles were recorded as 1/0 binary matrices, where 1 indicates the presence and 0 the absence of a given fragment. The resultant code, expressed as four binary digits (bits), describes the presence/absence of each fragment in the EcoRI/HpaII and EcoRI/MspI digests of the control materials and their Al-stressed counterparts. Theoretically, 16 permutations are possible for a binary code with four positions. The number of identical codes was summarized and presented in a form of a table (Table 1). Each code reflects a set of events affecting cytosines within a given restriction site recognized by HpaII and MspI in control and stressed samples. Demethylation (DM), de novo methylation (DNM), preservation of methylated sites (MSP), and preservation of non-methylated sites (NMSP) are the events that must be considered simultaneously.

The two isoschizomers digest methylated sites differently [31, 55, 56] depending on which cytosines are or are not methylated. For example, code (1111) represents a fragment unaffected by Al stress: the four cytosines in questions are non-methylated in both the control and stressed samples from that line. Similarly, code (1110) could reflect the preservation of methylation at an external position, or non-methylation of an internal position and de novo methylation of the other external cytosine. In both instances, to interpret the code, it is necessary to take into account four independent events. The other codes can be interpreted according to similar logic (Table 1). The total number of events [related to de novo methylation (DN), demethylation (DM), and preservation of non-methylated (NMSP) and methylated (MSP) sites] was calculated as the sum of all codes reflecting the given event type, where each code was multiplied by the number of corresponding events (Table 2). The sum of all events related to all 16 codes was calculated and used as the denominator for quantitative calculations. Each event type was then expressed as a percentage.

Eventually, it may become possible to quantify de novo methylation and demethylation of cytosines at symmetric CHG and CG sites. For this type of analysis, the methylation of cytosines in the CHG context was assumed when the external nucleotide in the 5′-CCGG-3′ sequence underwent change. When the internal cytosines of the restriction sites were affected the CG methylation was assumed. To quantify both cases, the denominator described above was used, and the results were expressed as percentages. A full interpretation of all codes is given in Table 1.

Statistical analysis

Agglomerative hierarchical clustering (AHC) based on Squared Euclidean distance and Ward’s Agglomeration method was performed in XlStat version 2015.6. ANOVA was performed using the Best model (adjusted R2) variable selection method, with a fixed intercept, a confidence interval equal to 95%, 0.0001 tolerance, and random validation settings in XlStat version 2015.6.







De novo methylation


Methylation sensitive amplification polymorphism


Preservation of methylated site


Preservation of non-methylated site


Non-tolerant triticale line


Tolerant triticale line


  1. Kooke R, Johannes F, Wardenaar R, Becker F, Etcheverry M, Colot V, Vreugdenhil D, Keurentjes JJB. Epigenetic basis of morphological variation and phenotypic plasticity in Arabidopsis thaliana. Plant Cell. 2015;27:337–48.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  2. Suter L, Widmer A. Phenotypic effects of salt and heat stress over three generations in Arabidopsis thaliana. PLoS One. 2013;8:e80819.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Zhang Y-Y, Fischer M, Colot V, Bossdorf O. Epigenetic variation creates potential for evolution of plant phenotypic plasticity. New Phytol. 2013;197:314–22.

    CAS  Article  PubMed  Google Scholar 

  4. Loenen WAM, Raleigh EA. The other face of restriction: modification-dependent enzymes. Nucl Acids Res. 2014;42:56–69.

    CAS  Article  PubMed  Google Scholar 

  5. Golldack D, Lüking I, Yang O. Plant tolerance to drought and salinity: stress regulating transcription factors and their functional significance in the cellular transcriptional network. Plant Cell Rep. 2011;30:1383–91.

    CAS  Article  PubMed  Google Scholar 

  6. Mirouze M, Paszkowski J. Epigenetic contribution to stress adaptation in plants. Curr Opin Plant Biol. 2011;14:267–74.

    CAS  Article  PubMed  Google Scholar 

  7. Iwasaki M, Paszkowski J. Identification of genes preventing transgenerational transmission of stress-induced epigenetic states. Proc Natl Acad Sci U S A. 2014;111:8547–52.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. Law JA, Jacobsen SE. Establishing, maintaining and modifying DNA methylation patterns in plants and animals. Nat Rev Genet. 2010;11:204–20.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  9. Wang W-S, Pan Y-J, Zhao X-Q, Dwivedi D, Zhu L-H, Ali J, Fu B-Y, Li Z-K. Drought-induced site-specific DNA methylation and its association with drought tolerance in rice (Oryza sativa L.). J Exp Botany. 2011;62:1951–60.

    CAS  Article  Google Scholar 

  10. Pan Y, Wang W, Zhao X, Zhu L, Fu B, Li Z. DNA methylation alterations of rice in response to cold stress. POJ. 2011;4:364–9. ISSN:1836–3644

    CAS  Google Scholar 

  11. Marconi G, Pace R, Traini A, Raggi L, Lutts S, Chiusano M, et al. Use of MSAP markers to analyse the effects of salt stress on DNA methylation in rapeseed (Brassica napus var. oleifera). PLoS One. 2013;8:e75597.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  12. Abratowska A, Wąsowicz P, Bednarek PT, Telka J, Wierzbicka M. Morphological and genetic distinctiveness of metallicolous and non-metallicolous populations of Armeria maritima s.L. (Plumbaginaceae) in Poland. Plant Biol. 2012;14:586–95.

    CAS  Article  PubMed  Google Scholar 

  13. Chwedorzewska KJ, Bednarek PT. Genetic and epigenetic variation in a cosmopolitan grass Poa annua from Antarctic and polish populations. Pol Polar Res. 2012;33:63–80.

    Google Scholar 

  14. Chwedorzewska KJ, Bednarek PT. Genetic and epigenetic studies on populations of Deschampsia antarctica Desv. From contrasting environments on king George Island. Pol Polar Res. 2011;32:15–26.

    Google Scholar 

  15. Bednarek PT, Orłowska R, Koebner RMD, Zimny J. Quantification of the tissue-culture induced variation in barley (Hordeum vulgare L.). BMC Plant Biol. 2007;7:10.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Fiuk A, Bednarek PT, Rybczyński JJ. Flow cytometry, HPLC-RP, and metAFLP analyses to assess genetic variability in somatic embryo-derived plantlets of Gentiana pannonica Scop. Plant Mol Biol Rep. 2010;28:413–20.

    CAS  Article  Google Scholar 

  17. Machczyńska J, Orłowska R, Mańkowski DR, Zimny J, Bednarek PT. DNA methylation changes in triticale due to in vitro culture plant regeneration and consecutive reproduction. Plant Cell Tissue Organ Cult. 2014a;119:289–99.

    Article  Google Scholar 

  18. Machczyńska J, Zimny J, Bednarek PT. Tissue culture-induced genetic and epigenetic variation in triticale (xTriticosecale spp. Wittmack ex A. Camus 1927) regenerants. Plant Mol Biol. 2015;89:279–92.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Choi C-S, Sano H. Abiotic-stress induces demethylation and transcriptional activation of a gene encoding a glycerophosphodiesterase-like protein in tobacco plants. Mol Gen Genomics. 2007;277:589–600.

    CAS  Article  Google Scholar 

  20. Boyko A, Kathiria P, Zemp FJ, Yao Y, Pogribny I, Kovalchuk I. Transgenerational changes in the genome stability and methylation in pathogen-infected plants (virus-induced plant genome instability). Nucleic Acids Res. 2007;35:1714–25.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  21. Hashida SN, Uchiyama T, Martin C, Kishima Y, Sano Y, Mikami T. The temperature-dependent change in methylation of the Antirrhinum transposon Tam3 is controlled by the activity of its transposase. Plant Cell. 2006;18:104–18.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. Zemach A, McDaniel IE, Silva P, Zilberman D. Genome-wide evolutionary analysis of eucaryotic DNA methylation. Science. 2010;14:916–9.

    Article  Google Scholar 

  23. Gent JI, Ellis NA, Guo L, Harkess AE, Yao Y, Zhang X, Dawe RK. CHH islands: de novo DNA methylation in near-gene chromatin regulation in maize. Genome Res. 2013;23:628–37.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  24. Johnston JW, Harding K, Bremner DH, Souch G, Green J, Lynch PT, Grout B, Benson EE. HPLC analysis of plant DNA methylation: a study of critical methodological factors. Plant Physiol Biochem. 2005;43:844–53.

    CAS  Article  PubMed  Google Scholar 

  25. Machczyńska J, Orłowska R, Zimny J, Bednarek PT. Extended metAFLP approach in studies of tissue culture induced variation (TCIV) in triticale. Mol Breed. 2014b;34:845–54.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Reyna-López 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:703–10.

    Article  PubMed  Google Scholar 

  27. Xiong LZ, Xu CG, Saghai Maroof MA, Zhang Q. Patterns of cytosine methylation in an elite rice hybrid and its parental lines, detected by a methylation-sensitive amplification polymorphism technique. Mol Gen Genet. 1999;261:439–46.

    CAS  Article  PubMed  Google Scholar 

  28. Cao DH, Gao X, Liu J, Kimatu JN, Geng SJ, et al. Methylation sensitive amplified polymorphism (MSAP) reveals that alkali stress triggers more DNA hypomethylation levels in cotton (Gossypium hirsutum L.) roots than salt stress. Afr J Biotechnol. 2011;10:18971–80.

    CAS  Google Scholar 

  29. Karan R, DeLeon T, Biradar H, Subudhi PK. Salt stress induced variation in DNA methylation pattern and its influence on gene expression in contrasting rice genotypes. PLoS One. 2012;7:e40203.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  30. Gayacharan, Joel AJ. Epigenetic responses to drought stress in rice (Oryza sativa L.). Physiol Mol Biol Plants. 2013;19:379–87.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  31. Fulneček J, Kovařik A. How to interpret Methylation sensitive amplified polymorphism (MSAP) profiles? BMC Genet. 2014;15:2.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Alonso C, Perez R, Bazaga P, Medrano M, Herrera CM. MSAP markers and global cytosine methylation in plants: a literature survey and comparative analysis for a wild-growing species. Mol Ecol Resour. 2016;16:80–90.

    CAS  Article  PubMed  Google Scholar 

  33. Schulz B, Eckstein RL, Durka W. Scoring and analysis of methylation-sensitive amplification polymorphisms for epigenetic population studies. Mol Ecol Resour. 2013;13:642–53.

    CAS  Article  PubMed  Google Scholar 

  34. Zhong L, Xu Y, Wang J. DNA-methylation changes induced by salt stress in wheat Triticum aestivum. Afr J Biotechnol. 2009;8:6201–7.

    CAS  Article  Google Scholar 

  35. Fulneček J, Matyášek R, Kovařik A. Distribution of 5-methylcytosine residues in 5S rRNA genes in Arabidopsis thaliana and Secale cereale. Mol Gen Genomics. 2002;268:510–7.

    Article  Google Scholar 

  36. Zhang MS, Yan HY, Zhao N, Lin XY, Pang JS, Xu KZ, Liu LX, Liu B. Endosperm-specific hypomethylation, and meiotic inheritance and variation of DNA methylation level and pattern in sorghum (Sorghum bicolor L.) inter-strain hybrids. Theor Appl Genet. 2007;115:195–207.

    CAS  Article  PubMed  Google Scholar 

  37. Feng Q, Yang C, Lin X, Wang J, Ou X, Zhang C, Chen Y, Liu B. Salt and alkaline stress induced transgenerational alteration in DNA methylation of rice (Oryza sativa). AJSC. 2012;6:877–83.

    CAS  Google Scholar 

  38. Hardcastle TJ. High-throughput sequencing of cytosine methylation in plant DNA. Plant Methods. 2013;9:16.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  39. Zhang Y, Jeltsch A. The application of next generation sequencing in DNA methylation analysis. Genes. 2010;1:85–101.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Tyrka M, Bednarek PT, Kilian A, Wędzony M, Hura T, Bauer E. Genetic map of triticale compiling DArT, SSR, and AFLP markers. Genome. 2011;54:391–401.

    CAS  Article  PubMed  Google Scholar 

  41. Kimatu JN, Diarso M, Song C, Agboola RS, Pang J, Qi X, Liu B. DNA cytosine methylation alterations associated with aluminium toxicity and low pH in Sorghum bicolor. Afr J Agric Res. 2011;6:4579–93.

    Google Scholar 

  42. Stroud H, Ding B, Simon SA, Feng S, Bellizzi M, Pellegrini M, Wang G-L, Meyers BC, Jacobsen SE. Plants regenerated from tissue culture contain stable epigenome changes in rice. elife. 2015;2:e00354.

    Google Scholar 

  43. Garg R, Chevala VVSN, Shankar R, Jain M. Divergent DNA methylation patterns associated with gene expression in rice cultivars with contrasting drought and salinity stress response. Sci Rep. 2015;5:14922.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  44. Sahu PP, Pandey G, Sharma N, Puranik S, Muthamilarasan M, Prasad M. Epigenetic mechanisms of plant stress responses and adaptation. Plant Cell Rep. 2013;32:1151–9.

    CAS  Article  PubMed  Google Scholar 

  45. Cao D, Gao X, Liu J, Wang X, Geng S, Yang C, Liu B, Shi D. Root-specific DNA methylation in Chloris virgata, a natural alkaline-resistant halophyte, in response to salt and alkaline stresses. Plant Mol Biol Rep. 2012;30:1102–9.

    CAS  Article  Google Scholar 

  46. Wang M, Qin L, Xie C, Li W, Yuan J, Kong L, Yu W, Xia G, Liu S. Induced and constitutive DNA methylation in a salinity-tolerant wheat introgression line. Plant Cell Physiol. 2014;55:1354–65.

    CAS  Article  PubMed  Google Scholar 

  47. Peng H, Zhang J. Plant genomic DNA methylation in response to stresses. Potential applications and challenges in plant breeding. Prog Nat Sci. 2009;19:1037–45.

    CAS  Article  Google Scholar 

  48. Hayatsu H, Wataya Y, Kazushige K. The addition of sodium bisulfite to uracil and to cytosine. JACS. 1970;92:724–6.

    CAS  Article  Google Scholar 

  49. Cokus SJ, Feng S, Zhang X, Chen Z, Merriman B, Haudenschild CD, Pradhan S, Nelson SF, Pellegrini M, Jacobsen SE. Shotgun bisulfite sequencing of the Arabidopsis genome reveals DNA methylation patterning. Nature. 2008;452:215–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  50. Lister R, O'Malley RC, Tonti-Filippini J, Gregory BD, Berry CC, Millar AH, Ecker JR. Highly integrated single-base resolution maps of the epigenome in Arabidopsis. Cell. 2008;133:523–36.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  51. Iyer LM, Kumpatla SP, Chandrasekharan MB, Hall TC. Transgene silencing in monocots. Plant Mol Biol. 2000;43:323–46.

    CAS  Article  PubMed  Google Scholar 

  52. Eichten SR, Springer NM. Minimal evidence for consistent changes in maize DNA methylation patterns following environmental stress. PMCID. 2015, PMC4422006.

  53. Niedziela A, Bednarek PT, Labudda M, Mańkowski DR, Anioł A. Genetic mapping of a 7R al tolerance QTL in triticale (X Triticosecale Wittmack). J Appl Genet. 2014;55:1–14.

    CAS  Article  PubMed  Google Scholar 

  54. Anioł A. Induction of aluminium tolerance in wheat seedlings by low doses of aluminium in the nutrient solution. Plant Physiol. 1984;75:551–5.

    Google Scholar 

  55. McClelland M, Nelson M, Raschke E. Effect of site-specific modification on restriction endonucleases and DNA modification methyltransferases. Nucl Acids Res. 1994;22:3640–59.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  56. Aversano R, Caruso I, Aronne G, De Micco V, Scognamiglio N, Carputo D. Stochastic changes affect Solanum wild species following autopolyploidization. J Exp Bot. 2013;64:625–35.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

Download references


This research was financed by the Polish Ministry of Agriculture and Rural Development.

Availability of data and materials

The data sets supporting the results of this article are included within the article.

Authors’ contributions

PB developed the concept of the study and the theoretical background of quantitative MSAP, and wrote the manuscript; RO did ANOVA and participated in describing the results; AN performed experiments, verified the theoretical data, and participated in writing the manuscript. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Piotr T. Bednarek.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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 ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Bednarek, P.T., Orłowska, R. & Niedziela, A. A relative quantitative Methylation-Sensitive Amplified Polymorphism (MSAP) method for the analysis of abiotic stress. BMC Plant Biol 17, 79 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Methylation Pattern
  • Methylation Change
  • Stressed Sample
  • Tolerant Line
  • Triticale Line