Phyton-International Journal of Experimental Botany

DOI: 10.32604/phyton.2022.016397

ARTICLE

Effect of UV Radiation and Other Abiotic Stress Factors on DNA of Different Wild Plant Species Grown in Three Successive Seasons in Alpine and Subalpine Regions

Svetla P. Gateva1,*, Gabriele Jovtchev1, Tsveta V. Angelova1, Tzvetana P. Nonova2, Nikolay Tyutyundzhiev3, Elena G. Geleva2, Kostadin Katrandzhiev1, Nina A. Nikolova2, Dimitar Dimitrov4 and Christo V. Angelov2

1Department of Ecosystem Research, Environmental Risk Assessment and Conservation Biology, Institute of Biodiversity and Ecosystem Research, BAS, 2 Gagarin Str., 1113, Sofia, Bulgaria
2Nuclear Scientific Experimental Base, Institute for Nuclear Researches and Nuclear Energy, BAS, 72 Tzarigradsko Chaussee Blvd, 1784, Sofia, Bulgaria
3Central Laboratory of Solar Energy and New Energy Sources, BAS, 72 Tzarigradsko Chaussee Blvd, 1784, Sofia, Bulgaria
4National Museum of Natural History, BAS, 1 Tsar Osvoboditel Blvd, 1000, Sofia, Bulgaria
*Corresponding Author: Svetla P. Gateva. Email: spetkova2002@yahoo.co.uk
Received: 03 March 2021; Accepted: 10 June 2021

1  Introduction

As a result of climate change, along with growing anthropogenic pressure, plants in the natural ecosystems are forced to survive in extreme environmental conditions, with rapid changes in UV radiation, ionizing radiation (IR), temperature, humidity, precipitation, etc. How do plants respond to the simultaneous changes of various environmental factors and do they differ in their tolerance to these factors? Gudkov et al. [1] commented that plants show modified responses under the combined effect of radiation and other abiotic stress factors. At the ecosystem level, the effects are less well understood than at the molecular and organismal levels. Many of the most important, yet less predictable consequences will be indirect effects of elevated UV-B rays acting through changes in the chemical composition and form of plants and through changes in the abiotic environment [2]. DNA is the main carrier of hereditary information in cells, which can undergo damage as a result of various environmental stress factors.

UV radiation plays a major role in normal plant photosynthesis, growth and development. The UV spectrum consists of UV rays with different wavelength. UV-B (280–315 nm) is about 1% of the total spectrum of UV radiation, whereas about 95% of UV radiation is UV-A (315–400 nm) [3]. Although, the exposure of plants to UV-A photons is 10–100 times more than to UV-B photons, the energy of UV-B radiation is high enough to induce numerous deleterious effects [4]. It is widely accepted that UV-C (100–280 nm) is fully absorbed by the Earth’s stratospheric ozone layer. Around 50 percent of global radiation is photosynthetically active radiation (PAR) which is in the range of 400 nm to 700 nm [5]. PAR is defined as the electromagnetic radiation in the range of the light spectrum that phototrophic organisms mainly use in photosynthesis.

In recent decades, the thickness of the stratospheric Earth’s ozone layer has undergone changes mainly because of the anthropogenic factor [6]. Its thinning leads to a higher amount of UV radiation, especially of UV-B, reaching the Earth’s surface. Increase in the level of UV-B radiation can lead to excess elevation of reactive oxygen species (ROS) and oxidative stress induction in plant cells [7]. As a result, cells can suffer DNA damage, damage of their photosynthetic apparatus, alteration of membrane properties, cell cycle arrest, morphological, physiological, biochemical changes, growth inhibition, and cell death [8–14]. To reduce ROS and oxidative stress induced by high UV-B radiation as well as to avoid damaging effects in the main cellular structures, many plants have developed several defense responses [6,10,14], increase in the levels of UV-absorbing compounds [9,15–19] and development of defensive morphological alteration [19]. Plants can activate a variety of DNA repair mechanisms, some specialized in the repair of UV damage [20,21]. The degree of damage in plant cells depends on the effectiveness of these responses [22].

UV-A light is less reactive than UV-B, and is less effective in damaging cellular structures [23,24]. In addition to the damaging effect in plants, UV-A light can stimulate induction of UV- absorbing pigments and antioxidants, enhance the stomatal conductance, etc. [25–28].

In addition to UV radiation, plants in natural ecosystems are exposed to a constantly changing natural background of ionizing radiation (IR) of geological as well as of cosmic origin. The main sources of natural beta-radioactivity are the natural long-lived isotope 40K and the radionuclide 238U and 232Th decay series [29,30]. In some areas, natural radioactivity is significantly increased due to the presence of radioactive elements in the rocks and/or soil [31,32]. The enhanced levels of natural radioactivity in the soil and water in the areas rich in natural radionuclides could affect plants, because radionuclides can easily pass into the plants from the soil through the roots and then to the other parts of the food chain [33].

The biosphere in the northern hemisphere has been contaminated by radionuclides of artificial origin such as 137Cs and 90Sr after nuclear incidents such as the Chernobyl nuclear power station accident and the Fukushima disaster [34–37] and the detrimental effect of IR on plants is well documented [34–40]. It affects DNA by direct energy deposition to DNA, leading to a single ionization in a DNA molecule, as well as indirectly-by ionization of water molecules generating multiple hydroxyl radicals and clustered DNA damage within a few base pairs [41,42]. Irradiation with higher doses leads to DNA SSBs and DSBs, which, if not repaired or if repaired incorrectly, can result in chromosomal aberrations and genomic alterations. The injuries induced by IR can be repaired and recovered but it depends on the type of irradiation (acute and chronic), dose rate, species radio sensitivity and the developmental stage [20,43,44]. In rice irradiated with low-dose radiation (4 μSv/h) in the area near Fukushima for three days, there was rapid activation of DNA reparation, antioxidant defense, photosynthesis, secondary metabolism, and cell death in the leaves [45].

Environmental factors change with altitude. The intensity of UV radiation, especially the UV-B fraction, increases because of the increased atmospheric transparency with the elevation. In this way, the plants that naturally grow at higher altitudes are exposed to a higher extent of both UV-A and UV-B radiation. On the other hand, many other factors are relevant to the amount of UV radiation reaching the surface, such as local meteorological conditions, cloud cover, season, time of day, solar zenith angle, latitude and air pollution [46]. With increasing height, the temperature shows extreme levels, the humidity and the moisture vary in a specific way, a different cloud mode is observed, the atmospheric pressure and partial pressure of O2 and CO2 decrease, the winds increase, the growth season duration decreases, the nutrients in the soil decrease. So in the areas with high altitudes, the plants are under combined abiotic stress [47]. This suggests a negative effect on the development and growth of the plants [48,49]. Barnes et al. [50] conducted UV screening of native and non-native species along an elevation gradient spanning 2600–3800 m to determine if certain species differ in their abilities to adjust their levels of UV protection in response to natural variation in UV-B exposure. Evidence was found that high levels of UV screening are not restricted to plant species native to the high UV-B conditions of the tropical alpine and that plasticity in epidermal UV transmittance is a mechanism employed by certain species to cope with varying solar UV exposures.

Plants have developed various adaptation strategies to overcome the environmental stress and to survive in extreme environmental conditions at high altitudes, such as synthesis of photoprotective compounds [16], activation of the plants’ antioxidant defense system, stress related hormones, osmotic adjustment, morphological alterations in some structures, etc. [17,19,48,51–54].

Many studies on the effect of UV light and altitude on plants exist assessing biochemical, physiological and morphological alterations but some of them are short-term and on a small number of species [28], others were conducted under controlled conditions in growth chambers and greenhouses [22,55], or the effect of only one stress factor was investigated. Ultimately, the question is whether the UV-induced reactions have been developed solely as a reaction to UV light, or is there a correlation between the high UV exposure and other loads [12]? Long-term monitoring studies at field conditions can provide more realistic assessment of the combined effect of UV radiation and other environmental stress factors (temperature, availability of water, soil type and composition, soil salinity) on plants from natural ecosystems at different altitudes. Consideration should be given to constantly varying environmental parameters such as temperature and humidity. Studies including the IR impact accompanied with others abiotic stress factors would be also useful to throw some light on plant adaptation to overcome environmental stress in areas of high altitude. Assessment of the effect on the main hereditary molecule DNA would give valuable information. In our previous study, we examined the effect of UV radiation, IR, temperature, humidity on DNA of wild plants growing at different altitudes in Rila Mountain as well as the extent of their adaptation for only one growing season [56].

The aim of the present study is to investigate DNA alterations of some wild plants from natural ecosystems at different altitudes in Rila Mountain, Bulgaria (1500, 1782, and 2925 m a.s.l.) exposed to UV radiation, IR and other changing abiotic stress factors such as temperature, deficient or excessive water in three successive growth seasons and to assess the extent of plant tolerance to the constantly changing environmental conditions. For this purpose, physicochemical, cytogenetic, and molecular methods were applied. Sensitive tests-induction of micronuclei and comet assay were used for assessment of the genotoxic effect of the combined impact of UV radiation, IR and other abiotic stress factors.

2  Material and Methods

2.1 Characteristics of the Sampling Sites and Plant Samples

Three sampling sites located at different altitudes in Rila Mountain, Bulgaria were chosen for the experiments.

The locality named “Skakavcite” is at 1500 m a.s.l. (42° 12′ 08.52″ N, 23° 33′ 06.38″ E). The average annual temperature at this experimental site is 8.6°C, and the average annual precipitation 667 mm. The soil type is brown forest with low clay content and high content of organic matter. Its pH is in the range of 4.5 to 6.0 (Rila National Park, Management Plan 2001–2010). In this area during the growth seasons of three years (2017–2019) we detected four tree species of four families, one shrub species, six Poaceae species, six Fabaceae species, two Juncaceae species, and 46 grass species of 23 families.

The locality named “Third window” is at 1782 m a.s.l. (42° 10′ 27.70″ N, 23° 33′33.20″ E). The soil type at this experimental site is dark mountain-forest soils, with acidic pH ranging between 4.0 and 5.0 [57]. For the three years’ period of our study, we detected there seven tree species from four families, three shrub species from three families, six Poaceae species, five Fabaceae species, two Juncaceae species, and 42 grass species of 25 families.

The third experimental site is situated at alpine altitude. It is the highest peak of Rila Mountain, Moussala Peak at 2925 m a.s.l. (42° 10′ 45.13″ N, 23° 35′ 06.81″ E). It is characterized with alpine type climate with average annual temperature of –3.0°C, and annual average precipitation of 1193 mm [57]. The soil type is mountain-meadow alpine with a pH value of approximately 4.5. Azonal scree and rock vegetation, calcareous and acid grasslands and sedge heaths, snow beds, soligenous mires, and spring vegetation is typical here. Above the tree line, the most important habitats are grasslands associated with late snow cover > 210 days/year. The species that grow here can survive at low temperature, low precipitation, higher atmospheric pressure, wind, poor soil, and short growing season. Only a few species of several plant families were obtained for the three years’ period of our study. The most species are from the Poaceae family, followed by acid grasses of Juncaceae, and sedges Cyperaceae, some species of Saxifragaceae, Asteraceae, Scrophulariaceae, and Boraginaceae.

2.2 Sample Collection

The following herbaceous plant species characteristic of the habitats at the three altitudes in Rila Mountain were collected:

- Epilobium angustifolium L. (Onagraceae), Fragaria vesca L. (Rosaceae), Myosotis sylvatica Ehrh. (Boraginaceae), Achillea clusiana Tausch (Asteraceae) and Dactylis glomerata L. (Poaceae) from habitats at 1500 m a.s.l. and at 1782 m a.s.l. in growth seasons June/August 2017, 2018 and 2019.

- Poa alpina L. (Poaceae), Sesleria coerulans Friv. (Poaceae), Saxifraga cymosa Waldst and Kit (Saxifragaceae), Myosotis sylvatica Ehrh. (Boraginaceae) and Achillea clusiana Tausch (Asteraceae) from habitats at 2925 m a.s.l., Moussala Peak in growth season August/September 2017. Festuca valida (R.Uechtr.) Pénzes (Poaceae), Achillea multifida (DC.) Boiss. (Asteraceae), Saxifraga pedemontana All. (Saxifragaceae) were collected in growth season August/September 2019.

A model plant, Hordeum vulgare L. (Poaceae), standard genotype (2n = 14), non-irradiated and unaffected by other abiotic stress factors was used as a control. The barley seeds were germinated for 19 h and grown to mature plants at standard laboratory conditions (12-h light/photoperiod/:12-h dark) at a constant temperature of 20–22°C in soil (alluvial-meadow) without any pest control techniques, watering during the growing stage and were used further in the experimental procedures given below. The experiments were conducted at the Institute of Biodiversity and Ecosystem Research, Sofia. This plant is widely used as a model test system in cytogenetics and molecular analyses.

Samplings with all wild plants as well as with the model one were performed on fully expanded leaves in mature plants after flowering, so all plants were at the same phenological stage.

The wild plant material for cytogenetic and molecular investigations was collected in moisture paper and after cooling at 4°C was put for prompt further laboratory analysis. The plants for gross beta activity determination were dried and subjected to the procedure given below.

2.3 Climatological Data

Climatological data (temperature and relative humidity) about Moussala Peak were provided by the Basic Environmental Observatory (BEO) Moussala located at the Moussala Peak (2925 m a.s.l.), which is in the frame of Institute for Nuclear Research and Nuclear Energy at Bulgarian Academy of Sciences (INRNE-BAS). Data were obtained by Vaisala QML201C Data Logger, equipped with sensor HMP45C – Temperature and Relative Humidity Probe mounted on a tripod placed 1 m above the ground. It measures relative humidity over the range of 0.8% to 100% RH and temperature over the range of −40 to +60oC. Data about the other two experimental sites were obtained using the database of the European Centre for Medium-Range Weather Forecasts (ECMWF) (https://apps.ecmwf.int/datasets/data/interim-full-daily/levtype=sfc/).

2.4 UV Parameters

Available data for average daily value of UV-A radiation and average daily UV-B radiation for each experimental site were investigated in various sources. The results are presented as Wh/m2. Information for these parameters was obtained from the SoDa database (provider: MINES ParisTech/Armines/Transvalor (France) (http://www.soda-pro.com) derived from broadband HelioClim-3 v5 Solar Irradiance (satellite data). Average daily values for five successive months (May, June, July, August and September) are calculated.

2.5 Gross Beta Activity Determination

The gross beta activity was determined for soils and plant tissues. The samples were measured after cleaning from mechanical impurities and drying. The plant samples were milled and the soil samples fractionated, homogenized and then ashed in a muffle furnace at 500°C [58]. Gross beta activity was measured by universal measuring device UVJ-01 with a Geiger detector with 0.23 cps Bq−1 efficiency for 90Sr/90Y (MK-30 Measuring Chamber) and background ≤2 cps (VF s.r.o. Žilina, http://www.vf.cz). The gross α/β activity in aerosol filters was measured by Multi-Detector Systems PIC-MDS-8 (Protean Instr. http://www.proteaninstrument.com), at background of 0.03–0.07 cpm for α and 0.4–0.7 cpm for β activity. The efficiency varies between 35% and 55% depending on the energy and geometry of sources. The values of gross beta activity are presented as Bq/g ash mean ± SD. All plant tissues were used for its determination.

2.6 Gamma Background Detection

At present the gamma background at BEO Moussala-INRNE-BAS is measured by a IGS-421 B1 gamma probe. The probe consists of two low-energy and one high-energy Geiger Müller tubes (ENVINET, Munich, Germany, http://www.ccsbg.com/userfiles/files/Data_Sheet_IGS421.pdf) with a sensitivity range of 10 nGy/h–10 Gy/h, with an accuracy of 15% with respect to 137Cs [59]. Data for gamma background is given in μSv/h.

2.7 Assessment of Plant DNA Susceptibility

To assess the genotoxic effect of the combined impact of UV radiation, IR and other abiotic stress factors and plant DNA susceptibility, respectively, cytogenetic and molecular methods were applied. Due to the specificity of different plant species both methods were modified by us for each of them.

2.7.1 Micronucleus Test

This sensitive and widely used technique for analysis of genotoxicity was adapted and applied for the following plant species: E. angustifolium, F. vesca, M. sylvatica, D. glomerata, P. alpina, S. coerulans, S. cymosa, A. clusiana, A. multifida, F. valida as well as S. cymosa and S. pedemontana.

The experiments were conducted as described previously by [56]. The time for hydrolysis varied depending on the plant species from 9 to 20 min. The time for maceration of Poaceae species was extended to 1 h (60°C). A total of 3000 micronuclei per plant species and experimental area were examined. The experiments were duplicated. The values of the micronuclei are presented as MN (%) mean ± SD.

2.7.2 Comet Assay (SCGE)

In our experimental setup we used the alkaline denaturation/neutral gel electrophoresis (A/N) variant of the comet assay, modified for barley by Jovtchev et al. [60]. This molecular method for detection of DNA damage in single cells was modified and adapted here for E. angustifolium, D. glomerata, A. clusiana, A. multifida, F. valida and S. pedemontana.

In brief, healthy leaves of each plant were carefully sliced in 400 μl of 1× PBS on ice in dim light and passed through a filter (0.45 μm) to isolate nuclei [56]. Thirty microliters of cell suspension from each plant sample were mixed with 60 μl of 1% low gelling agarose at 42°C and dropped on warm slides previously coated with 0.5% normal melting point agarose. Each drop was covered with a coverslip, and slides were stored on ice. After removal of the coverslips, leaf nuclei were denatured for 10 min at room temperature. The slides were washed triple in ice-cold 1× TBE for 5 min and electrophoresed for 6 min in1500 ml 1× TBE at 15–17 mA and 1 V/cm (31 V) at room temperature. After that the slides were placed in 70% and 96% ethanol for 5 min each, dried on air and stained with acridine orange (6 μg/ml). The comet analyses were done on a fluorescence microscope (Jenalumar, Carl Zeiss, Germany) at 480–490 nm and a digital camera (Canon, Power Shot A95, China). A total of 200–400 cells per each plant species and experimental area were evaluated. COMET image analysis system CASP (Comet Assay Software Project, https://www.casp.of.pl) software was applied to measure DNA content in the head and tail of each comet. The experiments were duplicated. The values of DNA in the tail are presented as (%) mean ± SD.

2.8 Statistical Data Analysis

Two-tailed Fisher’s exact test was used for group comparison of different plant species based on data obtained by comet assay, micronucleus test, and gross beta activity.

3  Results

3.1 Climatological Data

Climatological data for air temperature and relative humidity were collected in each growing season for a period of three years (2017, 2018 and 2019) for the three experimental areas. Clearly expressed temperature dynamics were observed for the growth seasons in the three years depending on altitudes (Fig. 1). The highest annual temperatures with low relative humidity were measured in 2017.

Figure 1: Average monthly values of temperature and relative humidity for the growth season of 2017, 2018, 2019 at “Skakavcite” locality A); B), and Moussala Peak C); D)

3.2 UV Radiation

The average daily values of UV-A and UV-B radiation in the period of five months (May, June, July, August and September over a period of three years) are given in Fig. 2.

Figure 2: Average daily values of UV-A and UV-B radiation detected at locality Beli Iskar “Skakavcite” 1500 m a.s.l.; locality “Third window” 1782 m a.s.l.; Moussala Peak 2925 m a.s.l. for the growth season of: A) 2017, B) 2018, C) 2019

The values of UV-A and UV-B radiation could vary throughout each season and each day. The information obtained about the average daily values both for UV-A and UV-B at all three experimental sites showed that the values of these parameters for the growing season of 2017 were close to that of growing season 2019 at all three altitudes, as in 2019 these average daily values were slightly lower. The lowest values were detected in May/September 2018.

Clear dependence on these values at the altitude was observed (Fig. 2), as the highest value both for UV-A and UV-B radiation was obtained at the highest experimental site Moussala Peak 2925 m a.s.l. for all the three successive active seasons (Fig. 2). As it is known, the maximum of UV radiation in the Northern hemisphere is in July/August. In Fig. 2 it can be seen that at the Moussala Peak 2925 m a.s.l. the average daily value of UV-A for July 2017 was 423.06 Wh/m2, and UV-B 9.94 Wh/m2, whereas at the lowest experimental site “Skakavcite”, UV-A was 388.55 Wh/m2 and UV-B 9.26 Wh/m2, respectively. For July 2018 these parameters at the highest experimental site were 376.65 Wh/m2 (UV-A), and 9.03 Wh/m2 (UV-B), whereas at the lowest site UV-A was 339.57 Wh/m2 and UV-B 8.25 Wh/m2, respectively. The same tendency was observed for 2019, namely UV-A value was 410.74 Wh/m2 and UV-B 9.68 Wh/m2 at 2925 m a.s.l., whereas at 1500 m a.s.l. UV-A value was 378.48 Wh/m2 and UV-B 9.06 Wh/m2.

3.3 Gamma Background and Gross Beta Activity

The gamma background was monitored during the experimental period of three years (2017–2019) and was measured on Moussala Peak and in Sofia. Its average value (0.152 ± 0.001 μSv/h) reported on Moussala Peak in June/September 2017 varied in very narrow ranges and was very close to that in 2019 (0.157 ± 0.002 μSv/h). The average values of gamma background in Sofia remained lower: 0.085 ± 0.001μSv/h in 2017 and 0.120 ± 0.002 μSv/h in 2019 (Fig. 3).

Figure 3: Average monthly values of gamma background for the growth season of 2017, 2018, 2019 at Sofia A), and Moussala Peak B)

The value of natural gross beta activity depends on the radionuclides in the rocks and soils. We detected the value of gross beta activity in the soil samples collected from the three experimental sites. The highest value was obtained in the soil samples from “Skakavcite” locality (0.950 ± 0.140 Bq/g dry matter), followed by that in the soils taken at the altitude of alpine type 2925 m a.s.l.-Moussala Peak (0.488 ± 0.069 Bq/g dry matter), and the lowest one was detected in the sample from the locality at an altitude of 1782 m a.s.l.-“Third window” (0.354 ± 0.068 Bq/g dry matter). Beta activity detected in the soil sample from Sofia at 595 m a.s.l. was relatively high, namely 0.796 ± 0.037 Bq/g dry matter.

We detected the value of gross beta activity in the plant samples collected from the three experimental sites for a period of three years (2017, 2018, 2019). As a control we used the model plant H. vulgare L. of Poaceae grown under standard laboratory conditions. The value of beta-activity detected in this control plant was 0.221 ± 0.008 Bq/g dry matter. The data obtained for gross beta activity in the wild plants are presented in Tab. 1.

The gross beta-activity was measured in the biennial herb M. sylvatica in plants grown at 1500 m a.s.l. (2018), 1782 m a.s.l. (2018 and 2019) and at 2925 m a.s.l. (2019). As can be seen from Tab. 1, this plant species of Boraginaceae had the highest beta-activity (p ≤ 0.001) compared to the other plants at all three experimental sites on different altitudes. The value of gross beta activity in this plant species varied within a narrow range, with no dependence on the annual growing season. A slight increase was observed with the elevation. The value of gross beta-activity in perennial plants E. angustifolium, F. vesca and A. clusiana measured at the two altitudes 1500 m a.s.l. and 1782 m a.s.l. did not show statistically significant differences in the samples from different altitudes and growth seasons (Tab. 1). Contrary to this, the species D. glomerata (Poaceaе) showed a lower value (p ≤ 0.001) of beta activity than the other grass species growing at 1500 m a.s.l. and 1782 m a.s.l., namely 0.324 ± 0.025 Bq/g at “Skakavcite” (2019) and 0.342 ± 0.051 Bq/g at “Third window”(2019), respectively (Tab. 1).

The results obtained for the plant material from Moussala Peak (2925 m a.s.l.) in 2017 and 2019, showed a lower beta activity for the species of Poaceae (P. alpina–0.413 ± 0.011 Bq/g; 0.227 ± 0.031 Bq/g, S. coerulans–0.394 ± 0.012; F. valida–0.401 ± 0.019 Bq/g) than other grass species. In comparison, A. clusiana 1.253 ± 0.031 Bq/g (2017) and A. multifida 0.819 ± 0.037 Bq/g (2019) of Asteraceae family demonstrated a significantly higher total beta activity (p ≤ 0.001) (Tab. 1).

The values of beta activity detected in the wild representatives of Poaceae grown in mountain conditions were close to that of the model plant H. vulgare.

3.4 DNA Susceptibility

To test the plant DNA susceptibility of wild plants to UV irradiation, IR and other abiotic stresses (temperature, availability of water, etc.) at different altitudes, micronucleus test and comet assay were used.

Data about the frequency of MN formation assessed by the micronucleus test in 11 wild plant species collected from three different altitudes in Rila Mountain in three consecutive annual growth seasons are given in Figs. 4 and 6A.

Figure 4: An overview of all observed micronuclei induced by UV radiation, IR and other abiotic stress factors in wild plant species from three altitudes in Rila Mountain, Bulgaria–1500 m a.s.l. A), 1782 m a.s.l. B), 2925 m a.s.l. C) in three successive years. The dashed line represents data from unaffected control plant H. vulgare aDifferences between species at the same altitude/cereals vs. other grasses/ bDifferences between the same species at different altitudes/cereals versus other grasses/***p ≤ 0.001; **p ≤ 0.01;-not significant.

Data obtained in 2017 and 2018 for the plants growing at the lowest altitude-1500 m a.s.l. showed that the most sensitive plant species was F. vesca (Rosaceae) where the value of detected MN was 0.330 ± 0.033% in 2017 and 0.210 ± 0.039% in 2018, followed by representative of Boraginaceae M. sylvatica (0.267 ± 0.088% in 2017 and 0.210 ± 0.043% in 2018), and E. angustifolium (0.200 ± 0.058% in 2017 and 0.160 ± 0.040% in 2018) of Onagraceae. D. glomerata (Poaceae) had the lowest susceptibility. The frequency of MN obtained in that plant samples was 0.066 ± 0.033% in 2017 and 0.100 ± 0.020% in 2018. The results in 2019 showed minor differences that are likely due to the climate change and the sampling time. E. angustifolium (0.221 ± 0.034%) and M. sylvatica (0.205 ± 0.069%) were found to be the most sensitive species, followed by D. glomerata (0.167 ± 0.024%), F. vesca (0.133 ± 0.023%) and A. clusiana (0.099 ± 0.041%).

Similar results were obtained in relation to the DNA susceptibility to UV radiation, IR, temperature, water availability in plants growing at 1782 m a.s.l. As can be seen in Fig. 4, F. vesca was again the most sensitive plant species (0.367 ± 0.067% in 2017 and 0.320 ± 0.030% in 2018), followed by E. angustifolium (0.267 ± 0.033% in 2017 and 0.330 ± 0.060% in 2018), M. sylvatica (0.300 ± 0.057% in 2017 and 0.200 ± 0.030% in 2018). The Poaceae representative D. glomerata (0.100 ± 0.058% for 2017 and 0.100 ± 0.030% 2018) had the lowest value of induced MN.

The sensitivity of the investigated plant species in 2019 showed the following gradation: M. sylvatica (0.398 ± 0.086%) > E. angustifolium (0.391 ± 0.058%) > F. vesca (0.128 ± 0.054%) > D. glomerata (0.100 ± 0.034%) > A. clusiana (0.098 ± 0.056%). The plants of Poaceae and Asteraceae possessed the lowest value of DNA damage.

The most sensitive species F. vesca and M. sylvatica had the highest values in growing season of 2017 compared to these detected in 2018 and 2019. Not significant differences in MN values were observed in E. angustifolium for each growing season at both lower altitudes (Fig. 4).

In the samples from the highest altitude Moussala Peak at 2925 m a.s.l., the Poaceaе species had lower susceptibility than the other plants. The frequency of detected MN both in the samples of P. alpina and S. coerulans were (0.033 ± 0.033%) in 2017 and in F. valida-0.034 ± 0.033% in 2019 (Fig. 4). A. clusiana (Asteraceae) did not show significant differences in the frequency of induced micronuclei (0.233 ± 0.033%) for 2017 compared to that in other subspecies from the same species investigated in 2019 (A. multifida – 0.223 ± 0.030%) (Fig. 4). The two representatives of the Saxifragaceae family, namely S. cymosa and S. pedemontana, show results close to the Asteraceae family (A. clusiana and A. multifida, grown at 2925 m a.s.l.) in the MN test.

It should be noted that the representatives of the same plant family show similar DNA sensitivity to injuries.

The frequency of MN observed in the non-irradiated model plant H. vulgare (Poaceae) grown under standard laboratory conditions was 0.089 ± 0.085%. Most of the wild species from all three altitudes in the Rila Mountains had higher value of MN than that detected in the control plant (Fig. 4), but D. glomerata (Poaceae) and A. clusiana (Asteraceae) collected at 1500 and 1782 m a.s.l. showed close MN to that of the control plant. It is interesting to note, that MN observed in the wild plant species from Poaceae grown at 2925 m a.s.l. had lower frequencies than that of the model.

The data obtained by molecular comet assay for strand DNA damage in the wild plants are presented in Figs. 5 and 6B. In 2017, we succeeded in adapting and applying comet analysis for two different species: E. angustifolium and D. glomerata, and in the following two years for another four species: F. valida, S. coerulans, A. clusiana and A. multifida.

Figure 5: An overview of all observed DNA damage induced by UV radiation, IR and other abiotic stress factors in wild plant species from three altitudes in Rila Mountain, Bulgaria–1500 m a.s.l. A), 1782 m a.s.l. B), 2925 m a.s.l. C) in three successive years. The dashed line represents data from unaffected control plant H. vulgare aDifferences between species at the same altitude/cereals vs. other grasses/bDifferences between the same species at different altitudes/cereals vs. other grasses/ cDifferences between species at different altitudes/cereals versus other grasses/ **p ≤ 0.01;-not significant

Figure 6: Images of micronuclei (A) and DNA damage expressed as DNA migration–comets (B) observed in various wild plant species exposed to UV radiation and other stressors at different altitudes in Rila Mountain

The results from three years of investigations on plant samples collected from 1500 m a.s.l. showed that DNA of E. angustifolium was the most susceptible to damage compared to the other tested species (Figs. 5 and 6B). The level of DNA damage assessed by the amount of DNA migrated in the comet tail was 38.737 ± 2.223% in 2017, 26.784 ± 1.998% in 2018 and 30.249 ± 2.643% in 2019, respectively. The DNA strand breaks in A. clusiana were found to be 18.009 ± 2.443% in the samples from 2019. D. glomerata, the representative of Poaceae family, had the lowest level of DNA in tail measured in the samples for all three successive seasons.

The molecular analysis of the samples collected from 1782 m a.s.l. showed that DNA of E. angustifolium is again the most sensitive compared to that of the other plant species (Fig. 5). The level of DNA in the tail was highest in 2017 (39.652 ± 1.231%), whereas in 2018 it was 30.462 ± 3.367% and 34.715 ± 2.554% in 2019. The amount of migrated DNA in the tail of the samples of D. glomerata was lower than that of E. angustifolium. The highest values were obtained in 2017 (23.450 ± 0.492%), whereas in 2018 and 2019, the DNA damage in D. glomerata had close levels 15.490 ± 2.270% and 13.840 ± 1.770%, respectively. A. clusiana had strand DNA damage of 21.899 ± 2.651% in 2017.

From Fig. 5 it can be seen that damage detected by comet assay in DNA of E. angustifolium and D. glomerata had the highest values in growing season of 2017 compared to these observed in 2018 and 2019 at both lower altitudes.

Analysis of DNA damage induced by UV and other abiotic factors in plants growing at 2925 m a.s.l. using comet assay was carried out on three species. The most susceptible to injuries was the DNA of A. multifida. The level of DNA in the tail was 20.009 ± 2.194%. The obtained results by this molecular technique showed that the DNA of S. pedemontana (level of DNA in tail was 10.944 ± 1.568%) and the cereal plant F. valida (11.500 ± 1.156% tail DNA) showed a significantly lower sensitivity to damage than the other plant species from this altitude.

The level of strand DNA damage in control unaffected plant H. vulgare (Poaceae) was 7.520 ± 2.159%. The obtained here results show that DNA damages in plants grown in mountain conditions are slightly, but still significantly higher than that in unaffected control plant H. vulgare. Plants from the highest altitude of 2925 m a.s.l. from Poaceae and Saxifragaceae families had close values of DNA damage to that of the control model plants.

Overall, based on the results obtained over the past three years by both cytogenetic and molecular test, it should be noted that the DNA of representatives from Poaceae is less sensitive to damage than other grasses at all three altitudes. These results correspond to the low values of the gross beta activity detected in these plant species.

4  Discussion

Alpine and subalpine plant species are of special interest in ecology and ecophysiology because they represent life at climate limit where various abiotic factors change particularly sharply with elevation. This could provoke diverse responses and/or adaptation of plants to changing climate and radiation conditions.

It is essential to study the combined effect of different environmental factors on the biota over time. DNA is the main target of UV radiation and IR radiation, potentially resulting in impairment of genomic integrity. Therefore, it is of particular importance to study the effect of interacting environmental factors on the plant’s hereditary material. Our study presents data based on three years of monitoring the effect of UV irradiation, ionizing irradiation of natural origin and other abiotic factors on DNA of wild plants at different altitudes in Rila Mountain, and assesses their tolerance to constantly changing environmental conditions.

The monitoring of natural background radiation levels, due to the presence of radionuclides in the rocks and soil gives information about the radiation load of the environment. The values of the natural radioactive elements 40K, 232Th, 238U and 226Ra and their decay product series in the environment depend on the local geological conditions, the type of rock from which the soil is formed [61]. On the other hand, anthropogenic waste generated from weapons and incidents in nuclear power plants [40,62,63] can increase the level of heavier long lived radioactive elements, e.g., 137Cs, 134Cs, 90Sr, and 238Pu, 239Pu, and hence increase the level of environmental contamination. Data obtained by us about the gross beta activity give information about the beta-emitting isotopes and dynamics of beta activity values in the soil and plants at three localities in Rila Mountain for the three-year study. The soil samples from the locality “Skakavcite” showed the highest value of gross beta activity, followed by that of Moussala Peak. The lowest one was detected in the soil from “Third window”. Gamma background radiation at the peak had a higher level for the three years’ study (0.152 ± 0.001–0.157 ± 0.002 μSv/h) compared to that of Sofia. The persistently high values found for three successive seasons both at Moussala Peak and “Skakavcite” is probably due to the peculiarities of the rocks in the Rila Mountains, as well as to the secondary ionization from the cosmic radiation typical of high altitude alpine areas. The rocks in the Rila Mountain and especially in the alpine part where Moussala Peak is located are predominantly granite gneiss [57]. Xinwei et al. [61] detected that igneous rock such as granite gneiss has a higher average value of radionuclide 40K than the typical world average value. 40K contributes also to the gamma radiation exposure because it is a gamma-ray emitter in addition to beta decays. The soil from the locality “Skakavcite” near Beli Iskar is known as a site with a high natural activity concentration of radionuclides such as 238U and 226Ra [64].

The plants at different altitudes in Rila Mountain showed a variation of gross beta activity and a species specificity independent of the altitude of the habitats. The biennial herb M. sylvatica (Boraginaceae) was highly susceptible compared to the other plant species and had the highest level of gross beta activity. This was manifested both in 2018 and 2019 with a trend to increase with the elevation. The species A. clusiana (Asteraceae) possessed a higher value (1.253 ± 0.031 Bq/g) of gross beta activity than the other plants growing at 2925 m a.s.l. in 2017. All grasses from Poaceae family–D. glomerata, P. alpina, F. valida, S. coerulans-had lower ability to accumulate radionuclides, where the values of gross beta activity were low regardless of the altitude. The variation of gross beta activity observed in the tested plant species was in agreement with the studies of Iovtchev et al. [65,66]. The authors obtained a broad variation in beta activity between phytomonitor plant species at different altitudes on Rila Mountain, which occur independently of the altitude of the habitats. Gupta et al. [67] reported in their review that plants differ in their capability to accumulate non-essential radionuclides from the soil. It depends on the plant life cycle, stage of growth, plant root structure, root-shoot ratio as well as on the reactivity of the radioisotope, which may affect the water solubility of the isotope and transport through the soil pore within the roots.

To detect DNA damage and compare DNA susceptibility of some wild plants from Rila Mountain to UV radiation, IR and other abiotic factors, we adapted and applied cytogenetic and molecular methods. Evaluation of DSBs is an essential step in the study of deleterious effects induced in DNA by the exogenous sources. Using the molecular test comet assay it was obtained that plant DNA of different species possessed different sensitivity to the combined effect of UV radiation and ionizing radiation. The highest value of induced DSBs was detected in E. angustifolium (Onagraceae), whereas the species from Poaceae family were more resistant. In the plants of Poaceae, the amount of DNA damage assessed by value of DNA migrated in the comet tail was less than in the other plants at both lower altitudes (1500 m a.s.l. and 1782 m a.s.l.). A similar low level of DNA damage was detected in the species F. valida (Poaceae) at the highest altitude of 2925 m a.s.l. It is interesting to note that the DNA injuries had the highest values in some plants in 2017, when the levels of some abiotic environmental parameters (temperature, UV radiation) were reportedly the highest, and the relative humidity was lower.

Micronucleus test was applied to most of the plant species from the three altitudes. The results showed different DNA species-specific sensitivity to the combined effect of UV radiation, IR and other abiotic stress factors (temperature, deficient or excessive water) in the natural environment. DNA of perennial herb F. vesca showed the highest sensitivity compared with other plants in the first two years of the study (2017, 2018) at both lower altitudes, whereas in 2019 these were E. angustifolium and M. sylvatica. The studied plant species of Poaceae demonstrated the lowest sensitivity to the stress factors with the lowest induction of MN. This was most expressed at the highest altitude 2925 m a.s.l. Moreover, the lower levels of DNA damage of Poaceae species observed both by cytogenetic and molecular tests corresponded with the lower ability of these plants to accumulate radionuclides, which is an indication of their higher ability to adapt to these environmental conditions.

It is interesting to note that a correlation was observed between the levels of induced DNA damage in plants and the change of UV-A and UV-B values and other environmental factors such as temperature and humidity following their dynamics of over a period of three years of our investigation. The plants grown in 2017 have shown the highest levels of DNA damage in most of the plants, which corresponded with the highest levels of the parameters for UV radiation (UV-A, UV-B) and temperature, and the lowest levels of humidity, suggesting water deficit. These values indicate that the effect was better pronounced in F. vesca, M. sylvatica and E. angustifolium, whose DNA was more susceptible to damage. Lower levels of DNA damage were detected in plants grown in 2018 and 2019, where the levels of the environment parameters such as temperature and UV radiation intensity were lower and humidity was higher. It is necessary to note that F. vesca and E. angustifolium do not occur at the highest altitude of 2925 m a.s.l. which is an indicator of their high susceptibility to the local environmental conditions and a lack of adaptation to survive in alpine regions.

The low level of DNA damage in some plants is probably due to some defense strategy developed by plants to become more resistant to constant variations of the abiotic factors and successfully adapt to the alpine environment.

The adaptation strategies to the elevation can be triggered in the plants as a result of combined effects of multiple environmental factors such as increase in UV radiation, increase in the atmospheric pressure, decrease in the temperatures, alteration in precipitation, alteration in CO2 and O2 concentrations with elevation [54]. Herbs that have a cushion-like habit, such as representative of the Saxifragaceae family, can also be significant to local species richness. Cushion plants are low-statured with small microphyllus leaves that are typically upright in orientation and often within a few cm of the soil surface on windy ridges. Munné-Bosch et al. [19] reported that the alpine plant species S. longifolia showed an increasing level of a-tocopherol, reduction in lipid peroxidation, and/or changes in the reproductive strategy. S. hostii exhibited a high photochemical efficiency as well as stomatal conductance under high UV radiation. These properties make the plant well acclimated with a better light-use performance at the alpine altitude [28]. Tusevski et al. [52] reported that the methanolic extracts from leaves and flowers of wild plants M. sylvatica, E. angustifolium and A. millefolium growing at a high altitude (1500–2000 m a.s.l.) in Jablanica Mountain possess high antioxidant activity. Additionally, the content of non-enzymatic antioxidants correlated with that of phenolic compounds. Cui et al. [54] obtained that the increase in proline and soluble sugar contents with the elevation may play a key role in osmotic adjustment in plants of alpine regions. The reduction in chlorophyll content and increase in the carotenе can help plants to reduce light absorption and to avoid photo-damage of the chloroplast.

A variation in physiological traits related to plant size is also observed to overcome the photo-oxidative stress during the periods of high light and low precipitation [68]. Deckmyn and Impens [69] investigated the effects of changes in solar UV-B on the growth and pigmentation of six grass species (Poaceae) from cold-temperate grasslands-among others F. arundinacea, F. rubra and D. glomerata, in spring and summer. They came to the conclusion that the sensitivity of the species is related to their ability to increase their leaf thickness in response to UV-B, whereas the change in concentration of protective pigments was species specific.

Spectral balance of PAR, UV-A and UV-B has also been shown to be important in determining plant sensitivity in field studies. A decrease in total biomass and plant height is usually associated with decreasing PAR and increasing UV-B. The protective effects of high PAR against elevated UV-B can be indirectly realized by increasing leaf thickness and concentration of flavonoids and other phenolic compounds [70].

Some plants are more adapted than others to the local environment at highest altitudes. We detected variability in DNA sensitivity in response to UV radiation, IR radiation and other abiotic factors (temperature, deficient or excessive water) among plant species possessing different genotype. Our study shows that plants from Poaceae and Saxifragaceae families growing at altitudes higher than 2000 m a.s.l. were more tolerant to combined stresses than other plants. These plants are UV adapted, which corresponds to the other high light related traits. The levels of damage detected in their DNA were lower and close to that in the control unaffected plant. In our previous laboratory investigations with UV-simulator using the same model test-system we obtained that the representative of Poaceae H. vulgare can adapt relatively quickly to UV radiation. Results from the comet assay showed that after an initial increase in DNA damage after 10 days of UV exposure in a dose close to that at the real conditions on Moussala Peak, the value of DNA migration (DNA in tail) after 20 days of irradiation was close to that of the unaffected control plants. Our data with four mature wild plant species from Poaceae grown in mountain conditions showed a lower frequency of induced micronuclei compared to that in H. vulgare after prolonged UV-A, UV-B irradiation (43 days, 27.408 ± 0.348 Wh/m2) in laboratory conditions. The induction of micronuclei in H. vulgare, on the other hand, showed a concentration-dependent effect, but there too was a significantly slower increase in the number of induced MN at the 20th day after the start of irradiation [71].

Monocots appear to be generally more resistant to UV-B enhancement than dicots [72]. Native species show pronounced resistance to increased UV-B radiation in relation to changes in biomass reduction and show minor changes in shoot elongation and leaf size. These morphological changes may have important consequences for natural ecosystems by altering the competitive balance in border areas [70]. Liu et al. [73] pointed out some advantages of Poaceae such as stomatal size and densities (affecting the closure speed), width of leaves, which make them more adaptive to changing light intensity, cold and drought. Teramura et al. [74] found that all species collected by them at an altitude higher than 2000 m a.s.l. were tolerant to UV–B radiation, where different plants vary in the degree of adaptation. Different species sensitivity has also been observed by other authors and in other plant species [21,75]. Rau et al. [76] showed that the content of UV-B absorbing compounds induced by UV-B irradiation was increased in the alpine species from Brassicaceae that became less sensitive to UV-B. The authors reported that the amount of the induced compounds depends on the genotype.

In the present study, A. multifida (Asteraceae) showed more susceptibility than other plant species at the same high altitude, but the fact that the plant survives in these alpine conditions can reveal an adaptation strategy. Sardrodi et al. [77] obtained that A. aucheri, another species of the genus Achillea, family (Asteraceae), growing at an altitude of over 3900 m a.s.l., had an increased essential oil content by 0.68%-0.95% based on dry weight, with increasing the height, and decreasing biomass at higher altitudes. This shows evidence of the role of advanced antioxidant systems in plants to overcome the stress and survive under adverse conditions. The altitude-dependent changes in the environment affect not only the plant physiology, but also the plant morphology [19]. Fitness analysis on Poahiemata (Poaceae) suggested that high altitude selection favored plants with short leaves and larger circumferences, whereas those at low altitude selected in favor of the opposite traits [78].

The ever-changing UV radiation and other abiotic factors may act antagonistically or synergistically in the plant adaptation. Štroch et al. [27] found that UV-A can mitigate the damaging effects of UV-B on the photosynthetic apparatus under photosynthetically active radiation (PAR) conditions in barley. The enhanced UV-B reduced plant growth, chlorophyll content and net photosynthesis rate, but these effects could be alleviated by higher temperature [79].

In this way, a complex mechanism of adaptation in plants was included.

5  Conclusion

Wild plant species’ specificity was detected to damage induced in the DNA by combined environmental stress factors in different altitudes in Rila Mountain, applying the adapted by us cytogenetic and molecular methods. The DNA from Poaceae plants was less sensitive to damage, suggesting a higher degree of adaptation to the environmental conditions than the other grasses at all three altitudes. The lower levels of DNA damage of Poaceae species corresponded to the lower ability of these plants to accumulate radionuclides. A particularly pronounced low level of DNA injuries in the plant species at the highest altitude indicated their successful adaptation to the alpine environment.

Following the study of DNA damage in plants over the three successive seasons it can be said that the effect of natural UV radiation and IR is probably impacted by other abiotic stress factors (temperature, deficient or excessive water). Our initial data are promising and further longitudinal studies are needed in order to understand more deeply the mechanisms related to plants’ responses to UV radiation and other abiotic factors as well as how the different factors interact and induce plant response/adaptation to the changing environmental conditions.

Author Contributions: Conceptualization and writing of the original draft, Svetla Gateva; Writing, review and editing, Gabriele Jovtchev; contributed to the analysis using MN test and comet assay, Tsveta Angelova; performed the analyses of gross beta-activity data and statistics, Tzvetana Nonova; sampling of plant material and analysis of the gross beta-activity data, Nikolay Tyutyundzhiev: contributed to the analysis of the UV data and approval of the final version, Elena Geleva; assembled the climate data, made the statistical analysis, Kostadin Katrandzhiev; assembled the plant material, edited the manuscript, Nina Nikolova; sampled the plant material and took part in data curation, Dimitar Dimitrov; performed project administration and contributed to the writing and approval of the final version Christo Angelov.

Funding Statement: This work was supported by a grant of the National Science Fund of the Republic of Bulgaria under Project No. DN 04/1, 13.12.2016 entitled: “Study of the combined effect of the natural radioactivity background, the UV radiation, the climate changes and the cosmic rays on model groups of plant and animal organisms in mountain ecosystems”.

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.

This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.