Significance of selenium supplementation in root- shoot reactions under manganese stress in wheat seedlings – biochemical and cytological studies

Agronomic practices are one of the reasons for the increasing accumulation of elements in the soil, including manganese (Mn). Our previous studies have shown that selenium (Se) ions can reduce the toxic actions of metal stress. Those, we studied the effects of Mn—treated as a stressor and Se – as a potential defense in plants. Mn ions (10 mM) or/and Se (15 μM) were added into hydroponic nutrients of two wheat cultivares. The evaluation of the stress-generating and protective actions were analyzed by biochemical methods and microscopic observations in leaves and roots. Moreover the level of DNA methylation for these tissues was determined. Mn application caused an increase of lipid peroxidation and hydrogen peroxide content in both leaves and roots and was accompanied with a greater absorption of this element by the roots. For other elements (K, Fe, S, P), with the exception of Ca, the reduced their uptake was registered, especially in roots. For roots, Mn stimulated greater, microscopically observed, desorganization in cell structure as compared to leaves, which was accompanied by a quantitative increase in 5-methylcytosine (5-metC) in root meristem. Se application diminished the effects of Mn-stress. These studies is the first in which indicated that global 5-metC level in roots enhancing from dividing meristematic cells to elongating cells of the axial cylinder and cortex. It was suggested that the rise in Ca level can lead to modification of root cells differentiations what may be one of the steps in defense mechanisms.


Introduction
Understanding the mechanism of action of biotic and abiotic stresses in plants still requires extensive research in order to unveil biochemical properties of cells that lead to minimization of toxic effects. The intensity of a stressor action, as well as specific response of species and even cultivars to the stress, are some of the reasons that make a detailed description of stages in stress adaptation process difficult. In most research papers overproduction of reactive oxygen species (ROS), is accepted as an important step in the mechanism the stress action (Nita and Grzybowski 2016). One of those ROS are the superoxide anion radical and H 2 O 2 responsible, among others, for inducing changes in the structure of cell membranes by increasing saturation of the fatty acids (recorded as an increase in MDA content) (Ayala et al. 2014). Thus, MDA level fluctuation is demonstrated as one of the factors informing about the stress intensity. The modification of the lipid structure may disturb the activity of protein channels and transporters located in the membranes and affect ion transportation. Among the ions especially important for proper functioning of cells are the cations of calcium (Ca), magnesium (Mg), and potassium (K), and the anions of sulfur (S) and phosphorus (P) (analyzed in mg concentrations), and also the cations of manganese (Mn), zinc (Zn), copper (Cu), and iron (Fe) (present in µg amounts). The role of Ca in plant cell homeostasis is particularly significant, as it mediates signals from the roots to the shoots in the process of adaptation to the environmental conditions and ROS overproduction (Sharma et al. 2017). Magnesium presence in activity centers of organic components is essential to stimulate many metabolic reactions that depend on protein synthesis, phosphorylation, photooxidation or enzyme activation (Cakmak and Yazici 2010). Similarly, K is involved in protein synthesis, enzyme activation and energy transfer, but also in osmoregulation and cation-anion balance, which were analyzed in aspects of biotic and abiotic stresses . Sulfur presence in amino acids of numerous enzymes increases plant tolerance to stress factors (Liang et al. 2016). The importance of P is mainly considered in terms of its contribution to the structure of nucleic acids. This macronutrient stabilizes the membrane structures and mediates the activation of enzymes and carbohydrate metabolism (Razaq et al. 2017;Malhotra et al. 2018). Membrane proteins containing S and/or P may serve as transporters in the translocation of various ions from and to cells (Sasaki et al. 2012). Zinc, copper and iron serve as co-factors of many enzymes, including those participating in antioxidant protection (Andresen et al. 2018). Zinc participates in the regulation of gene expression shaping plant tolerance under stress and is involved in protein synthesis and carbohydrate metabolism (Hafeez et al. 2013). Iron is a fundamental element of biomolecules controlling the structure and function of chloroplasts and synthesis of chlorophyll (Rout and Sahoo 2011). Uptake of Fe into cells may be supported by Cu ions (Scheers 2013).
In recent years, high demand for Mn has significantly increased emission of this element into environment as an effect of intensification of Mn ore mining mainly for Mn alloy production, mainly in welding, coke ovens and to production of dry alkaline batteries. Additionally, many countries use Mncontaining agrochemicals such as fungicides and fertilizers (Röllin and Nogueira 2011). The ion is a cofactor and activator of numerous plant enzymes, such as oxalate oxidase, oxalate decarboxylase, enzymes involved in fatty acid synthesis, or synthesis of deoxyribonucleic acid and ribonucleic acid (Broadley et al. 2012;Schmidt and Husted 2019). Manganese presence is indispensable to maintain the redox balance and electron transfer reactions, which is why it plays a valid role in the regulation of excess ROS generation under stress (Maathius and Diatloff 2013). It is suggested that exogenous application of Mn ions limited damage to plant cells by improving antioxidant defense systems, ion homeostasis and diminishing toxic changes in lipid peroxidation, particularly under metal stress conditions (Rahmann et al. 2016;Sebastian and Prasad 2015;Gzyl-Malcher et al. 2019). However, too high concentration of Mn, above the optimal content in cells (specific for plant cultivars), may induce destructive effects in plant tissues by initiating stress reactions. Similarly, as under other stresses, ROS generation and activation of oxidation processes were demonstrated under Mn excess in plant cells (Führs et al. 2009(Führs et al. , 2012Li et al. 2010;Bradley et al. 2012). It was suggested that Mn toxicity is related to oxidation of Mn in the apoplast, and that Mn may serve as an oxidizer of lipids and proteins (Fecht-Christoffers et al. 2003). The ability of Mn to bind with sulfhydryl groups of proteins disturbs their chemical structure and in consequence disrupts the movement of basic components in plant cells Sharma and Dietz 2009). A decrease in the absorption of elements such as Fe, Mg, Ca and P, characterized by a similar size of ions or the strength of ligand binding, was shown by Shi and Zhu (2008) and Millaleo et al. 2013). Chen et al. (2016) suggested that the modification of protein structure in the root system was the primary reason leading to growth inhibition in these organs. This suggests that the anatomical and chemical structure of the root system is essential to ensure proper uptake of substances from the soil and their transport to the leaves. Mn can be transferred via the plasma membranes of the root cells in a symplastic transport, or via the space between the root cells (apoplastic transport). Following absorption in the roots, Mn can be accumulated in the vacuole (as its main storage place) and/or the cell wall, and may also be transported to the leaves (via xylem). Rajpoot et al. (2021) showed that excessive Mn content (1-2 mM) in the growth media drastically changed the physicochemical conditions in the apical part of the roots.
One of the methods of protecting plants against stress induced by an excess of metals in the soil is the application of micro-amounts of selenium (Se) ions (Sieprawska et al. 2015;Mechora 2019). Such studies in the aspect of stress caused by Mn have not been conducted so far. Selenium is not required for plant physiological processes, but its presence is valued, especially in plants used as food products for humans and animals (Gupta and Gupta 2017). It is suggested that the mechanism of Se-based protection involves its blocking the uptake and transport of toxic substances (Koselski et al. 2019) and/or reducing the overproduction of ROS (Sieprawska et al. 2015;. Another mechanism may be a Se-mediated modification of the methylation pattern disturbed by stress factors (Filek et al. 2008;Kleiber et al. 2018).
DNA methylation regulates many biological processes in plants, including their growth and development (Grafi et al. 2011;Chen and Zhou 2013;Guzy-Wróbelska et al. 2013). Changes in DNA methylation were also documented in plants exposed to biotic and abiotic stresses, however, its direction (demethylation and de-novo methylation) depended on factors such as plant species, developmental stage and stress intensity (Karan et al. 2012;Greco et al. 2012). DNA methylation may enhance crop adaptation to stress by serving as a protective system that maintains genome integrity, and as the main epigenetic mechanism that controls gene expression in cells (Tan 2010;Shan et al. 2013). Moreover, it was shown that stress resistant genotypes responded to stress with a decrease in the percentage of 5-methylcytosine (5-metC), in the DNA pattern, while sensitive ones with an increase in methylation level (Gayacharan 2013). However, other studies demonstrated the same pattern of changes in cytosine methylation irrespective of the level of stress tolerance of the genotypes (Karan et al. 2012). Therefore, further research in this area should be carried out to learn more about the mechanism of action of both stress and adaptation processes.
The presented study aims to demonstrate a possible protective effect of Se ions in wheat subjected to toxic concentration of Mn and to estimate the influence of Mn and Se treatment on DNA methylation. Two cultivars of wheat, differing in visible response to Mn (including the length of leaves and roots and alterations in the size and number of chlorotic spots), were selected for analysis. The aim of the research was achieved by: i) determining the amount of hydrogen peroxide and the degree of lipid peroxidation as indicators of stress intensity; ii) assessing changes in the content of macro-and microelements, including Mn and Se, and modification of the membrane permeability as effect Mn and Se mediation in the transport of nutrients from the roots to the leaves; iii) measuring the content of soluble sugars and starch, important metabolites that may indicate disturbances in metabolic processes; iv) visualizing 5-metC content in the leaf and root cells; v) estimating the proportion of demethylated and methylated cytosine residues within CCGG sequence in DNA from the whole leaves and roots using Methylation-Sensitive Amplified Polymorphism (MSAP) method. Special attention was paid to the reactions occurring in the root system, the place of direct contact of Mn / Mn + Se with the plant tissues.

Plant material
Grains of spring wheat (cv. 'Nimfa', cv. 'Alibi'), obtained from a Polish Plant Breeding Station in Strzelce, were sterilized with 80% ethanol and 10% perhydrol, and rinsed four times with distilled water. After sterilization the seeds were placed in Petri dishes and germinated for two days, at 20 °C in darkness. Next, the seedlings were placed individually in the holes of polystyrene discs put on top of the vessels (10 l) filled with 50% Hoagland medium (1:1; v:v) as a control, or supplemented with 10 mM MnCl 2 (Mn treatment), 10 mM MnCl 2 + 15 µM SeCl 2 (Mn + Se treatment), or 15 µM SeCl 2 (Se). The seedlings were grown in the hydroponic conditions with a relative humidity 45-50%, light intensity 800 µmol (photon) m −2 s −1 (SQS, Hansatech Ltd, Kings Lynn, United Kingdom), and 16 h photoperiod (17/20 °C night/ day) until the phase of the third leaf (about 12 days). Then the plants were collected, shoots were separated from roots, tissues were weighed and frozen in liquid nitrogen for biochemical and MSAP analysis. Microscopic analyzes were carried out in fresh leaves and roots.

Determination of hydrogen peroxide
One gram of fresh weight leaves or roots was homogenized with 2 ml of a solution containing 0.25 ml trichloroacetic acid (TCA) (0.1% (w:v) (Sergiev et al. 1997). The homogenate was centrifuged (MPW 251R) at 12,000 × g for 15 min at 4˚C. Then 0.5 ml of the supernatant was taken and mixed with 0.5 ml KI (1 M) and 0.25 ml potassium phosphate buffer (10 mM, pH 7.2) at 4˚C. The quantitative determination was performed based on calibration curves for H 2 O 2 (0.0 -1.0 μmol . g −1 ) prepared in 0.1% TCA at λ = 395 nm.

Determination of soluble sugars and starch
For analyses of soluble sugars, a method described by Janeczko et al. (2010) with same modifications was used. About 1 g of fresh weight leaves were homogenized with 80% (v:v) ethanol and the extract was centrifuged at 5,000 × g for 10 min, at 4 °C. Then, the supernatant was heated at 95 °C with anthrone (0.2 g anthrone in 1000 ml 72% H 2 SO 4 ) for 10 min, and fast cooled to 20 °C. Spectrophotometric measurements (Thermo Scientific Evolution 201/220 UV-Visible) at λ = 625 nm allowed us to determine the content of carbohydrates. For quantitative calculations, a curve with glucose (1 mg . ml −1 ) as a standard was prepared.
Starch content was detected by the combined methods of Huber and Israel (1982) and Janeczko et al. (2010), in leaf samples (1 g) homogenized with 0.2 M KOH. After heating for 30 min, the samples were cooled and 1 M acetic acid was added to neutralize the mixtures. Starch depolymerization into glucose was carried out in a medium with 0.5 unit of α-amylase (EC 3211) for 15 h at 50 °C. After centrifuging at 6,000 × g for 5 min, glucose content was determined by the anthrone method, as described above.

Malonedialdehyde (MDA) analysis
One gram of fresh weight leaves and roots were homogenized in 0.5% (w/v) TCA (trichloroacetic acid) in a cold room according to Dhindsa et al. (1981) and Tobiasz et al. (2014). After homogenization, the samples were centrifuged at 5000 × g for 20 min, at 4 °C (MPW-351R). Then, 0.25% (w/v) TBA (thiobarbituric acid) in 10% (w/v) TCA was added into the supernatant and the mixture was heated at 95 °C for 30 min, quickly cooled on ice, and again centrifuged at 5,000 × g for 20 min. Absorbance of the supernatant samples was measured at λ = 532 nm and λ = 600 nm using a Thermo Scientific Evolution 201/220 UV-Visible spectrophotometer.

Macro-and micronutrient analysis
Portions (1.2 g) of lyophilized plant material were weighed (with an accuracy of 10 -5 g) and mineralized with 6 ml of concentrated nitric acid. Macronutrient content was detected by an inductively coupled plasma optic emission spectrometer ICP OES Optima 2100 (Perkin Elmer) at the following wavelengths: Mn 259.372 nm, Mg 285.213 nm, Fe 238.204 nm, Ca 317.933 nm, S 181.975 nm, P 213.617 nm in axial plasma observation mode; K 766.490 nm in radial observation mode, and with high argon flow through the monochromator system. Micronutrient content determination was based on Cu63, Zn68 and Fe57 isotopes as standards and inductively coupled plasma mass spectrometry (ICPMS; Elan DRC-e, Perkin Elmer, Shelton, USA) (Tobiasz et al. 2014).
Microscopic observations and determination of 5-metC localization (5-metC) The samples dissected from the leaves and roots were fixed with 4% paraformaldehyde (Polysciences) and 0.25% glutaraldehyde (Sigma) in phosphate-buffered saline (PBS), pH 7.2, for 24 h at 4 °C and prepared for embedding in BMM resin (butyl methacrylate, methyl methacrylate, 0.5% benzoyl ethyl ether (Sigma) with 10 mM DDT (Thermo Fisher Scientific)) according to Niedojadło et al. (2015). The samples were cut on Leica UCT ultramicrotome into serial semithin cross Sects. (1.5 µm), which were collected on Thermo Sci-entificTM Polysine adhesion microscope slides. The sections were hot stained with methylene blue for 10 min. They were flushed several times with distilled water, and then sealed in DPX Mountant for histology. The slides were assessed and documented with a BX50 microscope (Olympus) with NIS Elements AR 3.00 NIKON software.
For 5-metC localization the resin was removed with two changes of acetone and the slides were washed in distilled water and PBS pH 7.2. Before immunocytochemical reactions the sections were incubated with 0.1% RNAse for 30 min at 37˚C and washed in PBS pH 7.2 (3 × 5 min., 30 min at 37˚C). After neutralization of the medium with 0.1 M Na 2 B 4 O 7 × 10H 2 O (10 min at room temperature), the sections were incubated with anti-5-metC mouse (Abcam) primary antibody (1:200 in 1% BSA in PBS pH 7.2 for 12 h at 4˚C) and secondary goat anti-mouse antibody Alexa Fluor 488 (Invitrogen) diluted 1:1000 in 1% BSA -1 h at 37˚C. DNA was stained with 4,6-diamidino-2-phenylindole (DAPI; Fluka). Finally, the sections were washed in PBS pH 7.2, dried at room temperature and covered with Pro-LongTMGold antifade reagent (Thermo Fisher Scientific). The control reactions omitted the incubation with primary antibodies. Semithin sections were analyzed with an Olympus BX50 fluorescence microscope, with a UPlanFI 1009 (N.A. 1.3) oil immersion lens and narrow band filters (U-MNU, U-MNG). The results were recorded with an Olympus XC50 digital color camera and CellB software (Olympus Soft Imaging Solutions GmbH, Germany).

Quantitative evaluation of 5-metC content
The quantitative measurements of 5-metC levels were performed in the root apex cells. For this study, each experiment was carried out at the same temperature, incubation times, and antibody concentrations (Kozłowska et al. 2016). ImageJ (1.47v) software was used for image processing and analysis. Before quantification of the fluorescence intensity, the background was eliminated by adjusting a threshold according to autofluorescence based on the negative control. The fluorescence intensity was measured for five semithin sections per experimental variant at the same magnification and constant exposure time to ensure comparable results. For each zone of the root apex 15 to 30 cells were analyzed depending on the experiment. The level of fluorescence was expressed in arbitrary units (as mean intensity per μm 2 ). Figures were prepared in Microsoft Excel 2010.

DNA extraction and Methylation-Sensitive Amplified Polymorphism (MSAP) procedure
Tissues for the MSAP procedure were collected from three biological replicates. Each replicate consisted of leaves (second leaf) and the whole root system from three individual plants per treatment combination (control, Mn only, a combination of Mn and Se, and Se only). The leaves and roots were homogenized in liquid nitrogen and subjected to DNA extraction using micro-CTAB method (Doyle and Doyle 1987). Then, 250 μg of DNA were used for restriction with two combinations of enzymes in 10 μl of the reaction mixture: HpaII / EcoRI and MspI / EcoRI (New England Biolabs) in 1 × CutSmart Buffer with 1.25 U of each of the restriction enzymes. The restriction reaction was run at 37 °C for 3 h with subsequent inactivation of the enzymes at 70 °C for 10 min. Two technical replications of both restriction combinations were used to cut each DNA sample. The restriction fragments were then subjected to 16 h of adapter ligation at 37 °C, by adding 6 μl of a mix containing 30 pmol of HpaII/MspI adapter, 3 pmol of EcoRI adapter (Table A1), 0.6 U of T4 DNA ligase (Thermo Fischer Scientific) and 1 × T4 DNA ligase buffer. Then, 1 ml of DNA after ligation was used for preamplification reaction (10 ml of total volume) and mixed with 15 ng of HpaII/MspI + T and Eco + A primers (Table A1), 200 mM of each dNTPs (Promega), 0.25 U Taq polymerase (DreamTaq, ThermoFischer Scientific) and 1 × PCR buffer. The PCR conditions were as follows: 20 cycles of 30 s at 94 °C, 40 s at 56 °C and 50 s at 72 °C. The preamplified DNA was diluted 10 times and 2.5 μl were taken as a template for the selective amplification. In this step, 0.5 pmol of EcoRI selective primer labeled with a fluorescent dye IRD 800 (MWG), 12.5 ng of MseI primer, 200 μM of each dNTPs (Amersham Pharmacia), 0.3 U of Taq polymerase (DreamTaq, ThermoFischer Scientific), 1 mM MgCl 2 and 1 × PCR buffer supplemented with 2 mM MgCl 2 were used in a total volume of 10 μl reaction. The PCR was performed in the following conditions: 12 cycles of "touchdown" PCR of 30 s at 94 °C, 30 s at 68 °C (decreasing 0.7 °C per cycle) and 1 min at 72, followed by 23 cycles of 30 s at 94 °C, 30 s at 59 °C, and 1 min at 72 °C. All together 11 selective primer combinations were used in the analysis (Table A2). After the PCR, 4 μl of formamide gel loading buffer were added into the samples and, after denaturation at 95 °C for 3 min, 0.8 μl samples were loaded into polyacrylamide gels composed of 6% acrylamide/bis-acrylamide 19:1 solution (Sigma), 7 M urea (Amersham Pharmacia), and 1 × TBE buffer. Li-Cor DNA sequencer was used for the electrophoresis at the following parameters: 1400 V, 30 mA, 30 W. Medium speed of a laser scanning was used.

MSAP data analysis
To describe the changes in DNA methylation all MSAP bands that met the following two conditions were analyzed: (i) the band pattern was replicable in two technical replications of each DNA sample and (ii) the same band pattern was present in all three biological replicates representing particular treatment combination. Several profiles in MSAP loci were possible: (i) monomorphic loci, characterized by the presence of a band in all samples and in both restriction enzymes combinations; (ii) polymorphic loci of type I, characterized by the presence of a band in one restriction enzyme combination and the lack of a band in the second combination, while this presence/absence of a band was identical for all DNA samples. These loci reflected particular methylation status of cytosines in CCGG sequence recognized by HpaII and MspI restriction enzymes, but were unrelated to any methylation changes triggered by the treatment; (iii) polymorphic loci of type II, characterized by differences in the band pattern between treatment combinations, reflecting methylation changes in CCGG sequence related to the growth conditions. Only the last category of MSAP loci was taken into consideration.
The analysis of methylation changes was based on differences in sensitivity to cytosine methylation in CCGG sequence displayed by HpaII and MspI isoschizomers. Briefly, both enzymes can cleave non-methylated CCGG sequence. HpaII can also cut this sequence when external cytosine is hemimethylated (Mann and Smith 1977;Reyna-Lopez et al. 1997), and MspI is capable of cleaving this site when one or two internal cytosines are methylated (Reyna-Lopez et al. 1997). Such a profile does not allow for identification of the methylation status of each cytosine in CCGG sequence with a perfect precision. It allows, however, for quantification of the most probable changes in methylation using one of the existing models of MSAP data interpretation. In our analysis, we followed the model proposed by Bednarek et al. (2017) which takes into consideration all possible methylation statuses of both external and internal cytosines in CCGG sequence that may underline the presence or the lack of MSAP band after HpaII and MspI digestion. Based on this model, the percentage of demethylation (later encoded as de-Met) and de novo methylation (de novo-Met) events that may occur after a particular treatment vs. the control conditions was calculated. Additionally, this model also allows for quantifying the probable share of demethylation and de novo methylation of cytosines at symmetric CHG and CG sites (later encoded as de-Met-CHG and de-Met-CG or de novo-Met-CHG and de novo-Met-CG). A more detailed description of the data quantification method is provided in the Supplemental Tables.

Statistical analysis
The results of biological experiments were reported as an average of triplicate repetitions (± SD) and analyzed statistically using one-way ANOVA, followed by post-hoc (Duncan test) comparisons with Statistica Software version 9.1 (SAS Institute Inc., Cary, NC, USA). The significance level was set at p < 0.05.

Biochemical analysis
No significant visual differences between the tested cultivars were found in the seedlings grown under the control conditions. This was also demonstrated in the similarities of the fresh weight of their leaves and roots. Only a slightly higher weight (statistically insignificant) was recorded in cv. Nimfa's leaves as compared with cv. Alibi (Table 1). Growth in the presence of Mn decreased the weight of both leaves and roots by about three times (leaves) or about seven times (roots), independently of the cultivar. Moreover, both cultivars showed signs of leaf edge yellowing. Damage was also visible in the root system, mainly in cv. Alibi. In this cultivar, shortening and deformation of the root system and its stronger browning were more significant than in cv. Nimfa. In Mn + Se media, a weight increase of the examined organs was noted as compared with Mn only treatment. For plants grown on the media containing only Se, fresh weight of both leaves and roots was comparable to that of the control.
The control seedlings accumulated a greater amount of Mn in the roots than in the leaves ( Table 1). The cultivation on Mn-containing media resulted in a significant increase in the content of this element in both organs, especially in the root systems of the investigated cultivars. The changes were more pronounced in cv. Nimfa than in cv. Alibi. Following joint application of Mn and Se, Mn level slightly increased in the leaves and decreased in the roots of cv. Nimfa, while in cv. Alibi the opposite effect was observed (as compared with Mn treatment). The presence of Se alone did not affect the content of Mn in An analysis of Se content in the tested cultivars showed only trace amounts of this element in the control and Mn variants. In the roots, Se administration, both alone and with Mn, increased Se content by about six to nine times for cv. Nimfa and cv. Alibi, respectively. In the leaves, the presence of Mn more than doubled decreased the Se uptake, as compared with Se added individually, especially in cv. Alibi.
The studied cultivars differed in their content of macro-and micronutrients (Table 1). In the control leaves of cv. Nimfa, higher concentrations of Ca, Mg, Zn and Fe were recorded than in cv. Alibi, while the root system of cv. Nimfa was richer in Mg, S, Cu and Fe than that of cv. Alibi. Mn presence restricted the accumulation of all tested elements in the leaves of both cultivars. In the root system, an increase in the content of P, S and Zn (both cultivars) and Ca (cv. Nimfa) was recorded. Simultaneous treatment with Mn and Se reduced the content of P, Cu and Fe in the leaves of the tested cultivars and additionally of K in cv. Nimfa, in comparison with the plants exposed to Mn alone. Moreover, in the leaves of cv. Alibi a drop in the content of Ca, Mg and S, was noticed. In the root system, there was an increase in the level of P and Fe for both cultivars, while the concentration of the other elements decreased in cv. Nimfa and rose in cv. Alibi (as compared with Mn alone). In the presence of Se alone, the content of all elements in both leaves and roots was close to that recorded for the controls.
An analysis of H 2 O 2 level indicated its much more abundant presence in the leaves than in the roots of all studied plants (Fig. 1A). In control plants, no significant differences in the content of this compound between the tested plants were found. After treatment with Mn both cultivars showed increased concentration of H 2 O 2 in the leaves, and in cv. Alibi it rose by about 1.5 times. In the roots of both cultivars the content of H 2 O 2 grew by about three to five times (more in cv. Nimfa). When Mn and Se were administered simultaneously, the concentration of H 2 O 2 in the leaves of both cultivars decreased in relation to that observed for Mn (especially in cv. Alibi). In the roots, hydrogen peroxide level changed only to a small extent as compared with Mn-only treatment. The presence of Se added individually did not significantly affect the level of H 2 O 2 of the studied cultivars and its concentration in both leaves and roots was close to that in the control.
The change in MDA level, which is considered an indicator of the level of lipid oxidation by ROS, was slightly higher in cv. Alibi both in the roots and leaves (Fig. 1B). After growing in Mn media, MDA increased by about three times in the leaves and by about 1.5-2 times in the roots of both wheat cultivars. Cultivation on Mn + Se media diminished MDA values in all tested cultivars, as compared with Mnonly conditions, except for cv. Nimfa leaves, where no significant differences between Mn + Se and Mn treatment were found. Cv. Nimfa plants grown on Se-only media had similar MDA levels to controls in Data are averages from three independent biological replications ± SE. The same letters indicate not significant differences between treatments (p ≤ 0.5) both leaves and roots, while in cv. Alibi a decrease in MDA (vs. control) was registered.
The contents of soluble sugars and starch are presented in Table 2. Control leaves of both cultivars contained similar levels of starch, while significantly higher values of soluble sugars were found in cv. Alibi. In Mn-only media, an increase in the level of both simple and compound carbohydrates was indicated. It was especially visible for soluble sugars, the amount of which spiked five times following the metal treatment. Joint treatment with Mn and Se caused a drop in starch concentration (vs. Mn-only variant) in both cultivars. The presence of Se, as the only supplement to the media, did not significantly affect the concentration of any tested carbohydrates.

Microscopic observations and analysis of DNA methylation
Microscopic observations after staining with tolouidine blue revealed numerous chloroplasts in the mesophyll cells of both cultivars grown under control conditions (a representative example is shown in Fig. 2A). Figure 2A' shows chloroplasts localized under the membrane uniformly stained with toluidine blue. The unstained areas of the plastids indicate the presence of starch grains. Mn induced visible changes in the chloroplasts, as it boosted starch accumulation and the formation of amyloplasts in mesophyll cells (Fig. 2B, C, D). In cv. Alibi the response varied widely, ranging from slight changes in the plastids where starch accumulated (Fig. 2B') to complete degeneration of the leaves with numerous amyloplasts in their cells (Fig. 2C, C'). In cv. Nimfa, the vast majority of plastids were changed and accumulated starch (Fig. 2D, D').
In Mn + Se variants, plastids in cv. Alibi were mainly amyloplasts (Fig. 2E') or small, nearly uniformly colored structures (Fig. 2E''). In cv. Nimfa, only amyloplasts were found in the cells of the assimilation parenchyma (Fig. 2F, F'). The presence of Se did not change the microscopic image of the leaves in any of the cultivars, which was the same as that recorded for the control (data not shown).
In the leaves of both wheat cultivars, 5-metC was localized in the cell nuclei of all tissues, i.e., in the epidermis, assimilation mesophyll and vascular bundle (Fig. 3). There were no significant differences in the nucleus signal intensity in the control leaves (Fig. 3A, B) and those grown in the presence of Mn (Fig. 3C, D) and Mn + Se (Fig. 3E, F). Unexpectedly, microscopic analysis of these sections revealed that Mn uptake resulted not only in starch accumulation in mesophyll plastids. The membranes of the transformed plastids showed autofluorescence, which may indicate changes in their structure (Fig. 3C', D'). In either cultivar, chloroplasts observed in the control mesophyll cells did not show autofluorescence (Fig. 3A, B). These changes in plastids were also seen in the leaves of plants grown in the presence of Mn + Se (Fig. 3E', F'). In both cultivars, amyloplasts ( Fig. 3E and E', F and F') were mainly localized in the cells of the assimilation parenchyma.
Microscopic observations of the control roots ( Fig. 4A) revealed similar anatomy for both cultivars, with the growth apex surrounded by cap cells, and the meristem composed of rows of cuboid-shaped cells (Fig. 4A', A''). Above those, there was a zone of cell extension and differentiation with procambium and procortex layers. In Mn treated plants, the root anatomical structures showed significant differences between the cultivars. In cv. Alibi roots (Fig. 4B), the meristematic cells did not have a regular shape (Fig. 4B'') and the apex was surrounded by a smaller number of cap cells, often with irregular shapes and deformed nuclei. Moreover, the area between the meristem and the zone of differentiation comprised groups of oval-shaped cells that were larger than the surrounding cells and absent in the control (Fig. 4B'). Mn uptake led to drastic anatomical changes in cv. Nimfa roots (Fig. 4C). These changes were often accompanied by degeneration of major growth apexes. In such roots, numerous cap cells were damaged and the meristematic cells were of abnormal shapes without visible nuclei (Fig. 4C''). Additionally, a characteristic group of compact, small cells with large nuclei at the meristem differentiation border was found (Fig. 4C'), which was not observed in the control roots. When Mn was provided together with Se, the anatomy of the root growth apex in cv. Alibi (Fig. 4D) and cv. Nimfa (Fig. 4E) was more similar to the control, and there were no distinguished groups of cells above the meristem. In cv. Alibi, the meristematic cells were of correct shape (Fig. 4D'), whereas irregularly shaped cells but with visible nuclei were still present in cv. Nimfa (Fig. 4E').
In the roots a similar pattern of 5-metC localization (Fig. 5.I.A-C) and levels ( Fig. 6A and B) was observed in the plants grown under control conditions. A very weak signal was localized in the meristematic and cap cells (Fig. 5.I.A). A successive increase in 5-metC was found in the nuclei of cells that entered the differentiation and elongation pathway (Figs. 5.I.B and C, 6A and B). In the root elongation zone, the strongest signal was detected in the nuclei of the cells of the differentiating vascular cylinder (Fig. 5.I.C). No signal was localized in control material, incubated without anti-5-metC antibodies (data not shown). In the roots of cv. Alibi grown in the presence of Mn the increase in DNA methylation was observed from the differentiation zone to cell elongation zone, similarly as in control (Figs. 5.II.A-C and 6A). In the cap and meristematic cells, DNA was only slightly more methylated than in control (Fig. 6A). However, the tissues exhibiting anatomical abnormalities (Fig. 4) showed also changes in their DNA methylation pattern. In the cap zone the enhanced signal was visible only in single cells (Fig. 5.II.A), while the group of oval-shaped cells, not present in control (Fig. 4B'), showed lower level of 5-metC than the surrounding cells (Figs. 5.II.B, 6A -"different cells"). The elevated level of 5-metC was detected only in the cells on the procambium-procortex border (Fig. 5.II.B). In the elongation zone, the cells were irregular, some with no visible nuclei, and some with nuclei with a strong 5-metC signal (Fig. 5.II.C).
Mn uptake leading to anatomical changes in cv. Nimfa roots (Fig. 4) was also reflected in 5-metC signal intensity in the nuclei (Figs. 5.III.A-C, 6B). A significantly stronger signal was detected in the damaged cells of the root apex, i.e., the cap cells with irregular shapes and deformed nuclei and single meristematic cells with nuclei still present (Fig. 5.III.A). Additionally, higher level of 5-metC was found in the specific group of small cells with highly condensed DNA (Fig. 4C') on the meristem -differentiation zone border (Figs. 5.III.B, 6B -"different cells"). In cv. Nimfa, in which the anatomy of these root zones was slightly impaired in some roots (Fig. 4.E), the level of 5-metC was also slightly higher than in control (Figs. 5.V.A, 6B). In both cultivars, no significant disturbances in the differentiation and root elongation zones were found (Fig. 5.IV.B-C, V.B-C) and, similarly to the control, a progressive increase in DNA methylation was observed in the differentiating cells (Fig. 6).
A quantitative assessment of 5-metC level is shown in Fig. 6. In both cultivars, regardless of the treatment, we found an increase in methylation level from the root cap to the differentiation zone. Mn presence generally enhanced methylation intensity. In cv. Alibi, significant changes were restricted to the differentiation zone, while in cv. Nimfa they occurred throughout the root apex. A particular increase in 5-metC was recorded in the root cap (ca. three-fold) and in the apical meristem (ca. 2.5-fold). Under Mn + Se conditions, methylation levels were similar to those found in the control.
The analysis of DNA methylation profiles using MSAP method showed no changes in the band pattern in the leaves collected from control plants and Fig. 4 Microscopic observation of wheat roots grown under control conditions and after treatment with 10 mM MnCl 2 and 10 mM MnCl 2 + 15 µM SeCl 2 . cc -cap cells, m -meristem, pc -procambium, pcx -procortex from those grown in the presence of Mn, Mn + Se or Se (data not shown). In this organ, only two types of loci were found: monomorphic and polymorphic of type I, which reflected a particular methylation status of cytosines in CCGG sequence recognized by HpaII and MspI restriction enzymes, but unrelated to methylation changes triggered by growth conditions. Contrary to that, in the roots we found 101 polymorphic loci of type II, representing differences in methylation of CCGG sites. One of these loci was polymorphic between cvs. Alibi and Nimfa, but did not show any methylation changes between the treatments: the bands were present in all cv. Alibi samples and absent in all cv. Nimfa samples, and thus this locus was not taken into account in further analysis. The other two loci also distinguished the cultivars of Alibi and Nimfa. In the first locus, absent in cv. Nimfa, the bands were present in both restriction combinations in cv. Alibi roots treated with a combination of Mn and Se. The second locus was absent in cv. Alibi, whereas in all cv. Nimfa samples the bands were present after HpaII/EcoRI digestion only. The bands in the remaining loci were present in both cultivars and represented various patterns of MSAP profiles. Altogether 100 polymorphic loci were used to calculate the proportion of demethylation (de-Met) and de novo methylation (de novo-Met) of cytosines at CCGG sites. Additionally, the proportion of symmetric methylation changes of cytosines in CG and CHG context (where H represents internal cytosine in CCGG sequence) was estimated taking into account: (i) demethylation of CHG context (de-Met-CHG), (ii) demethylation of CG context (de-Met-CG), (iii) de novo methylation of CHG context (de novo-Met-CHG), and (iv) de novo methylation of CG context (de novo-Met-CG).
The methylation status of CCGG sites in cv. Alibi and cv. Nimfa roots was comparable and changed in a similar manner in all treatments in both cultivars (Table 3). A higher proportion of Nimfa (B) root growth apex (control, Mn -manganese, Mn + Se -manganese + selenium). Data are averages from three independent biological replications ± SE. The same letters indicate not significant differences between treatments (p ≤ 0.5) demethylation than de novo methylation events was detected at CCGG sites after Mn and Mn + Se application regardless of the genotype. At the same time, more changes in methylation status were noticed for the joint treatment with Mn and Se than for Mn-only treatment. Plant cultivation on the media supplemented with Se alone resulted in the lowest proportion of changes with similar number of demethylation and de novo methylation events. An analysis of the proportion of methylation changes at symmetric cytosines also showed a greater percentage of de-methylation than de novo methylation events at CHG and CG sites following Mn treatment. Under Mn + Se conditions, the proportion of all tested methylation events increased, especially in the context of de novo methylation of CHG sites. However, the presence of Se only reduced the amount of methylation changes in comparison with the controls, particularly of demethylation or de novo methylation of CHG context (Table 3). When the analysis was limited to the 20 polymorphic loci that represented the genotype specific changes in the methylation events, some minor differences in the proportion of methylation changes were noticed as compared with those described above (Table 4). The combined Mn + Se treatment resulted in a higher proportion of de novo methylations of CCGG sequences in cv. Nimfa roots. Also, Se-only treatment increased the proportion of methylation changes in cvs. Nimfa and Alibi, including the higher number of demethylation or de novo methylation events at CHG and CG context, except for demethylations at CHG in cv. Alibi roots.

Biochemical analysis
The effect of a stressful stimulus on a plant triggers a cascade of signals, the effect of which is the modification of metabolism enabling survival in stressful conditions. Plenty of stress factors trigger similar morphological alteration in plants such as: weight loss, yellowing or browning, as well as tissue deformation and shortening demonstrated i.e., under metal stress by Mishra et al. (2006), Selvaraj (2018), Huang et al. (2019).
In our study, the tissue browning indicated that the applied dose of Mn was stressful for both tested cultivars, and slightly more intense in cv. Alibi. Greater differences between the cultivars involved the response of their root systems. The larger loss of weight (in relation to the control) noted in cv. Alibi, may suggest that for this cultivar Mn treatment was more stressful than for cv. Nimfa. This conclusion is consistent with the generally accepted assumption that the decrease of tissue mass is an indicator of the intensity of toxic effects in plants (Jaleel et al. 2009). Selenium used at the selected concentration did not affect the appearance of tissues (macroscopic and microscopic observations) or their fresh weight. However, increased weight of both leaves and roots (vs. Mn-only treatment) and reduction of visual effects of Mn, indicate possible Se-induced protection against Mn stress. Thus, our data suggest the possibility of using Se to improve plant resistance to Mn excess. Yao et al. (2012), who analyzed Mn level in grapes, indicated that depending on the cultivar, either more of this element was retained in the roots or was Table 4 The comparison of methylation status in cvs. Nimfa and Alibi root samples expressed as the percentage of methylation changes between control samples (C) and each of the treatments (Mn only, Mn + Se and Se only). The data are based on the 20 polymorphic loci that represent the genotype specific changes in the methylation events distributed to the leaves, even in the resistant plants. It is suggested that the accumulation of greater amounts of stressors (metal ions, toxins) in tissues of susceptible cultivars, as compared with resistant ones, is the main cause of their greater damage (Sieprawska et al. 2017;Kim et al. 2009). Also, in the studies of Millaleo et al. (2010), the correlation between Mn content in selected organs and their destruction occurred even in the same genotype. In our experiments, greater amount of Mn accumulated in the root system than in the leaves of the tested plants may indicate that the roots were more sensitive to Mn stress. However, higher level of Mn (in both leaves and roots) was found in cv. Nimfa than in cv. Alibi, even though cv.
Alibi responded with a greater weight loss. Moreover, the addition of Se (suggested above as a potential protector) further increased the accumulation of Mn. However, the distribution of this element between leaves and roots depended on the cultivar, which is in agreement with the suggestions of Yao et al. (2012). Mn translocation is also influenced by the chemical form of Se introduced into the culture medium. As shown by Guerrero et al. (2014), the exposure of wheat plants to selenate enhanced Mn root uptake, whereas selenite treatment reduced Mn levels in both roots and shoots, with a particularly strong effect in the shoots. Thus, it seems that it cannot be unequivocally concluded that greater accumulation of Mn is the only indicator that differentiates sensitivity of the studied cultivars to this stressor. Also, specific properties of the membranes of the investigated cultivars should be taken into account (especially permeability, depending on the content of non-saturated lipids and the presence of ion transporters). Yang et al. (2014) researched specific transporters participating in Mn translocation. Independently of the uptake by Mntransporters, Mn can be also absorbed by transporters of other ions of similar value and size. The change in MDA is considered an indicator of the membrane permeability and fluidity, which depend on the content of non-saturated fatty acids (Ayala et al. 2014). A decrease in lipid unsaturation (recorded as an increase in MDA value) may increase the membrane stiffness and, consequently, modify the activity of membrane transporters and ion channels. The decrease in the nutrient content in the presence of Mn correlated with a decrease in membrane permeability, indirectly shown as an increase in MDA concentration. Significant negative correlations (r = -0.802 to r = -0.991, p < 0.05) between accumulation of the investigated elements in the leaves and Mn uptake (both alone and in combination with Se) suggested "blocking the transport" of the other elements under Mn stress. However, in the roots of both cultivars, also positive correlations of Mn with Zn, P and S (r = 0.844 to r = 0. 999; p < 0.05), and additionally with Ca (r = 0.951, p < 0.05) for Mn + Se variant were shown. Moreover, in the roots of cv. Nimfa, a correlation between Mn and Ca (r = 0.997) was noted when Mn was administered alone. The reduced uptake of the studied elements under Mn stress was demonstrated earlier in the leaves of various plant species (Juice et al. 2006;Kleiber 2014). Also, in chloroplasts of wheat seedlings, diminished accumulation of Zn, Cu, Fe, Mg and Ca was shown under Mn stress (Sieprawska et al. 2017). Altered concentrations of these nutrients may disturb metabolic processes and consequently curb plant growth as manifested in a decrease of fresh weight in the presence of Mn.
Interestingly, the increased uptake of Ca by the roots of both cultivars under Mn + Se treatment, and in cv. Nimfa roots under Mn stress was noted. The importance of Ca ion as a signal molecule (so-called "secondary messenger"), under stress was described in numerous papers (Thor 2019;Wilkins et al. 2016). As revealed by White et al. (2002), this ion is uptaken from a growth media through plasma membrane channels located in the roots. This may explain why the increase in its amount was evident in the roots. The presence of Se ions can also activate calcium channels partly blocked by Mn, as shown for another stressor (Koselski et al. 2019). Additionally, in Mndamaged roots, extra Ca uptake stimulated by Se may contribute to "rebuilding" of the tissues. The increase in the concentration of Ca ions during the initiation of developmental processes was demonstrated in many studies (Hepler 2005;Lamport and Varnai 2013;Thor 2019). This role of Ca may explain the appearance of meristematic centers in the roots of cv. Nimfa treated with Mn (Fig. 4). Ca influx may be also regulated by the generation of ROS during stress reactions (Demidchik et al. 2007). Among ROS H 2 O 2 , similarly to Ca, is the molecule recognized as an important element of cell signaling (Demidchik et al. 2007), independently of its ability to oxidize biomolecules and disrupt their functioning. Thus, a significant increase in H 2 O 2 concentration in the presence of Mn, especially visible in the roots (as compared with control), may indicate both greater intensity of stress in these organs, but also, as in the case of Ca, its involvement in "chemical signaling process" leading to limited tissue destruction (roots).

Microscopic observations and analysis of DNA methylation
Special function of the roots was indicated by the fact that microscopic observations of the leaf structures showed smaller damage of these organs in the presence of Mn than the roots. Main changes in the leaf tissues involved chloroplast formation and the amount and distribution of starch in these organelles. The quantitative changes and redistribution of starch in chloroplasts were found under various environmental stressors . Our earlier studies also showed a genotype resistance-dependent increase in the accumulation of this polysaccharide (Sieprawska et al. 2017). Small differences regarding the chloroplasts modification of the investigated plants treated with Mn, were mainly registered as changes in distribution of these organelles in the cells. Starch distribution within the chloroplasts (Fig. 2) correlated with the increase in its concentration, registered biochemically (Table 2). It may indicate a comparable susceptibility of both cultivars to Mn. However, in cv. Alibi we also found partly damaged leaves with changes in the chloroplasts typical of less tolerant plants. The stimulation of protective mechanisms under Mn stress was evidenced by an almost twofold increase in the concentration of soluble sugars in both cultivars. These molecules are regarded as osmoprotective agents, the accumulation of which may reduce negative effects of many stressors (Rosa et al. 2009). Relatively high level of soluble sugars accompanied by reduced starch content and smaller microscopically visible damage in the leaves may confirm the involvement of Se in the protective reactions against Mn stress, as suggested based on the experiments described above.
Larger anatomical differences between the examined cultivars exposed to Mn were observed in the roots. In cv. Alibi, changes in the shape of meristematic cells were noted, and a reduction in the number and extent of structural damage of the cap cells, as well as the appearance, in the area between the meristem and the differentiation zone, of additional (compared to the control) groups of oval cells, larger than the surrounding cells (Fig. 4B'). Their presence may indicate a potential change in the direction of cell differentiation in this area of the root. Moreover, in the elongation zone, cells were often irregular in shape, some without visible nuclei. Such modifications may suggest that Mn induced faster apoptosis of the cap cells and impaired differentiation of root tissues. In cv. Nimfa roots treated with Mn we found a decrease in the number of cap cells in comparison with control (Fig. 4C). A characteristic change, observed only in this cultivar, was the presence of a group of small cells, with large nuclei and condensed chromatin, visible at the border of the meristem and differentiation zone (Fig. 4C'). Presumably, these are procortex cells that differentiate into a new root apex. Creation of inappropriate root branching as a consequence of damage to the apical meristem were found by Clarkson (1969) and Llugany et al. (2003). We suggest that the formation of the new root apexes could be an element of plant defence mechanisms allowing for plant regeneration even after death of the primary roots. It may be related to the increase in Ca content in these cultivars observed in the presence of Mn, as suggested above. Independently of the stress reactions induced by Mn, the protective role of Se was also revealed in anatomical observations of the roots, as a reduction of changes caused by the presence of Mn.
The anatomical modifications induced by Mn in the leaves of the studied cultivars, mainly in the form of starch accumulation in plastids, were not accompanied by changes in DNA methylation determined based on the global level of 5-metC detected with immunocytochemical method and the analysis of methylation status of CCGG sequences with MSAP method.
Contrary to that, our analysis of wheat roots indicated that the global 5-metC level varied in different root zones and tissues. Kumpf and Nowack (2015) reported a very low level of DNA methylation in dividing, undifferentiated meristematic cells and short-lived cap cells, which soon enter the path of apoptosis. However, we showed for the first time a progressive increase in DNA methylation from dividing meristematic cells to differentiating and elongating cells of the vascular cylinder and cortex, as measured by the global 5-metC level. The role of DNA methylation in the processes of plant tissue differentiation is well known (Brown 1989;Ikeuchi et al. 2015). Elevated levels of 5-metC in the meristem and cap cells in the presence of Mn likely reflect the process of their apoptosis. The higher 5-metC content in cv. Alibi was observed only in a few cells of the cap and meristem, while in cv. Nimfa it was found in all visible cell nuclei. We may then conclude that in this cultivar Mn absorption could lead to cell apoptosis and death of the main growth apex. The phenomenon of DNA methylation during apoptosis was described in Danio rerio cells subjected to malnutrition stress (Fan et al. 2020).
Moreover, microscopic observations indicated that in both cultivars, cells of the differentiation zone had the potential ability to change their differentiation pattern in the presence of Mn, probably towards formation of an additional growth apex. This process seemed to be faster in the cultivar Nimfa, in which Mn stress strongly damaged meristematic root cells. Mn uptake resulted in faster apoptosis of the cap cells and impaired differentiation of root tissues. These processes were accompanied by changes in DNA methylation. Our study is the first to report a progressive increase in DNA methylation from dividing meristematic cells to elongating cells of the axial cylinder and cortex, as measured by global 5-metC level. The increase in 5-metC methylation signals after the treatment with Mn may result from the activation of pathways that counteract the transposon movement. A characteristic feature observed in both cultivars under Mn was the presence of groups of distinct cells in the meristem-differentiation zone. In cv. Alibi they were larger cells (Fig. 4B') with slightly reduced 5-metC level (Figs. 5IIB, 6A). In cv. Nimfa they were a group of small cells with highly condensed DNA (Fig. 4C') and increased level of DNA methylation (Figs. 5IIIB, 6B). It can be assumed that these are cells whose differentiation was altered. Mn induced process of changing the direction of differentiation was likely accompanied by changes in DNA methylation. In cv. Nimfa, with the death of the main grows apex, these cells are probably differentiating into new root apexes. The formation of a new root apex could serve as a defense mechanism for the plant in the event of the main root apex dieback.
Contrary to the increase in cytosine methylation reflected by global level of 5-metC in the root nuclei, our MSAP analysis showed the predominance of demethylation events at CCGG sites after Mn and Mn + Se treatments in the roots of both cultivars.
These discrepancies may be explained in two ways. The MSAP analysis was based on the homogenized material from the entire root system, thus the results were averaged across all root zones and all tissue types, whereas the quantification of global level of 5-metC was directed to specific zones and tissues. The second possible explanation is related to the resolution of the analysis and the type of sequence context where the methylated cytosine may be present. Here, the MSAP method gives much higher resolution than in situ hybridization, as it allows us to detect the methylation status of sequences located in very close proximity to each other. On the other hand, the MSAP method, which utilizes HpaII and MspI restriction enzymes, is directed for the analysis of CCGG sequences only and in situ hybridization allows for detecting also other sequences with methylated cytosine residues. This kind of explanation seems even more plausible, as it indicates that the stress related to high doses of Mn may preferentially induce de novo methylation of cytosines in other sequence contexts, for example in asymmetric CHH sites (where H indicates any nucleotide except for G). In many species, the asymmetric CHH sites are often present in heterochromatin and within or nearby transposons and retrotransposons, which need to be silenced to prevent the instability of the genome (Zhang et al. 2018). We may therefore hypothesize that the increase in 5-metC methylation signals after the treatment with Mn may result from the activation of pathways that counteract the transposon movement. At the same time our in situ hybridization analysis showed that the addition of Se to the medium with Mn resulted in similar methylation levels as found in control. This may support the hypothesis of the protective role of Se ions in wheat subjected to high concentration of Mn, as in their presence the activation of the mechanisms of genome stability protection may not be necessary.
Following up on this hypothesis, we may assume that at least some CCGG sequences targeted by the MSAP method may be located in the gene coding regions. It was showed that methylated CG clusters are present in the gene body of many constitutively expressed genes (Tran et al. 2005), and that methylation of CHG sequences can also correlate with gene expression activation (Bewick and Schmitz 2017). The demethylation events detected in our MSAP analysis after Mn and Mn + Se treatments may reflect the necessity of differential regulation of gene expression at particular loci in the presence of these ions. It is also known, however, that the link between methylation level of CG and CHG sites and the level of gene expression is highly gene-specific (Bewick and Schmitz 2017), which mitigates the straightforward interpretation of such results. Nonetheless, the complex changes in DNA methylation patterns that emerge from our study encourage a much deeper analysis based on DNA bisulfide sequencing in order to identify specific genes that may be regulated by differential methylation status in response to Mn and Se ions. This may support the hypothesis of protective role of Se ions in wheat subjected to high concentration of Mn, as in their presence the activation of the mechanisms of genome stability protection may not be necessary.

Conclusions
The presented research demonstrates that the root system of wheat seedlings is where most plant reactions to toxic doses of Mn take place. The responses were indicated as a considerable reduction of root weight and disturbances of their macroscopic and microscopic structure, especially within the apical parts. Se ions may have a protective effect on plants also under Mn stress. The mechanism of action of Se ions was not related to a direct blockage of Mn influx to the cells. Se-stimulated uptake of Ca ions seems to be an important step in this process, enabling regeneration of root tissues.
Microscopic observations indicated the potential ability of the meristem zone of roots to modify their differentiation pattern under Mn stress, more significantly in the case of strong damage to the meristematic and cap cells. These changes were related to a progressive increase in global DNA methylation from dividing meristematic cells to elongating cells of the axial vascular cylinder and cortex, what was indicated for the first time in our studies. Based on the DNA methylation profiles detected by MSAP, we hypothesize that different types of methylation sites may undergo different changes in the methylation pattern. We presume that asymmetric CHH sequences may come under preferential methylation in the presence of high doses of Mn, which is blocked or slowed down by Se addition, whereas other types of sequences, i.e., CGs and CHGs may undergo preferential demethylation both under Mn and Mn + Se treatment.
In summary, it may be assumed that Mn ions accumulate mainly in the cells of the "main" root apex. After exceeding its maximum concentration, the meristematic cells die (apoptosis), but new root buds differentiate rapidly and compensate for the loss of the dead cells. Thus, the plant can still take up Mn that is transported to the leaves. This may explain why in cv. Nimfa roots, where we observed the formation of new meristematic centers, the toxic processes were not as pronounced as in cv. Alibi, despite greater Mn uptake. Another explanation may be the possibility of Mn storage in degenerated fragments of roots. Such effects were not seen in cv. Alibi, which absorbed less Mn but transported it to the leaves. Data availability Not applicable.

Conflicts of interest
The author declares no conflict of interest.

Informed consent statement Not applicable.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.