The genetic diversity of Asplenium viride (Aspleniaceae) fern colonizing heavy metal-polluted sites

DOI: https://doi.org/10.21203/rs.3.rs-1165755/v1

Abstract

Heavy metals can affect the morphology, physiology and evolution of plants. Asplenium viride is a diploid species, belonging to the largest genus of the cosmopolitan fern family Aspleniaceae, and occurring on various types of alkaline rocks. It is known to colonize sites with high concentrations of heavy metals, exhibiting changes in frond morphology. Microevolutionary processes, manifesting as changes in genome size and new genotype formation, can ultimately lead to the formation of new subspecies and speciation. This study aimed to evaluate the morphological and genetic diversity of A. viride, and to test for a potential correlation between variability and heavy metal concentration. Analysis of A. viride specimens, from one metalliferous site and five non-metalliferous localities, showed no apparent variation in genome size between plants from affected and non-affected sites. There was no significant correlation between genetic variability and heavy metal concentration. This was possibly due to intragametophytic selfing, caused by patchy habitats and subsequent founder effects, resulting from long-distance colonization by single spores.

Introduction

Contamination with heavy metals can have very significant effects on the metabolism, ecology and microevolution of plants. An excess of heavy metals can exert harmful effects on plants based on the disturbance of important metabolic pathways, by interfering with the functions of various enzymes (DalCorso 2012). This results in significant differences between affected and control plants, in morphological and anatomical features, such as growth, habit, size and thickness of leaves and the density of root hairs (Meyer et al 2010; Przedpełska and Wierzbicka 2007).

High concentrations of heavy metals also have significant effects on plant reproductive processes and fertility, which are manifested as disorders of male and female gametophyte development (Czapik 2002; Koivunen et al. 2004). Microevolutionary processes induced by the presence of heavy metals often result in changes in the structure and number of chromosomes, the formation of new cytotypes (Coulaud et al. 1999), changes in genome size (Temsch et al. 2010), and the generation of new genotypes (Hildebrandt et al. 2006; Słomka et al. 2011).

Asplenium L. is the largest genus belonging to the cosmopolitan family of Aspleniaceae Newm., which comprises about 720 species, and is found in humid tropical and temperate zones across all continents (Kramer and Viane 1990; Smith et al. 2006). Some species have a sub-cosmopolitan range, or occur as endemics on islands and oceanic archipelagos such as Madagascar, New Zealand, Juan Fernandez and Tristan da Cunha (Brownsey 1977; Kramer and Viane 1990; Roux 1993). The mountainous regions of Asia and the Americas are secondary centers of Aspleniaceae species diversity (Kramer and Viane 1990; Werth et al. 1985). The diversity of genome sizes within the Asplenium genus range from 8.22 pg/2C (A. ceterach L.) to 30.17 pg/2C (A. monanthes L., Bennett and Leitch 2012; Dyer et al. 2013). This variation can be explained by a high variability of chromosome number and multiple polyploidization events during the speciation of the Asplenium genus. The ploidy of Asplenium species ranges from diploid to octoploid (base chromosome number x = 36; Brownsey 1977; Smith et al. 2006).

Reticulate evolution plays an important role in diversity, as many species are of hybrid and allopolyploid origin, and frequently exhibit parallel and recurring polyploid speciation (Dyer et al. 2012; Perrie et al. 2010; Schneider et al. 2013). Asplenium viride Huds. is a diploid species (2n = 2x = 72), which occurs mainly in the mountainous regions of Europe (Ivanova and Piȩkoś-Mirkowa 2003; Valentine and Moore 1993). In Poland, the species’ range includes the areas of the Sudetes and the Carpathians, the Kraków-Częstochowa Upland, the Lublin Upland and the Świętokrzyskie Mountains (Zając and Zając 2001). It occurs in rock crevices, in damp and shaded places. A. viride prefers alkaline rocks, mostly limestone (Valentine and Moore 1993), and serpentine rocks (Johnston and Proctor 1979; Marin and Tatić 2001; Proctor and Woodell 1971), being an indicator of ultrabasic soils in some regions, for example in Scandinavia (Lepp 2001).

The aim of this study was to (1) estimate ploidy level and genome size, and evaluate morphological and genetic variability of selected populations of Asplenium viride, as well as to (2) ascertain whether there is a correlation between heavy metal concentration and specimen variability.

Materials And Methods

Study areas and sampling

Plant material was collected from both heavy metalliferous (HM) and non-metalliferous (NM) sites in the summer of 2013. The polluted sites (heaps of waste Triassic dolomite rocks left by historical Zn-Pb ore mining) are located in western Olkusz Upland (Southern Poland), in the vicinity of the former mining village of Galman (G), approx. 5 km north-east of the town of Trzebinia. Detailed characteristics of the geology and climate of the region, the history of mining, the mine waste heaps as well as the vegetation covering them are described in Woch (2015) and Woch et al. (2017). Plants at the metalliferous site were collected from five heaps, scattered over approx. 2 km2. Three fern specimens were taken from each heap and pooled to obtain one sample per heap (n = 6). Additionally, three soil samples (0–15 cm in depth) were collected from the same heaps, and pooled into one composite sample per heap. Control ferns were taken from five non-metalliferous sites located in Southern Poland: 1) the Czupel (C) mountain (n = 4), located in the Beskidy Małe mountains, 6 km south-east of the city of Bielsko-Biała, 480 m ASL, the Carpathian flysch; 2) the Cisowe Skały (CS) mountain (n = 2) located in western Olkusz Upland, 5 km north-east from the town of Trzebinia, 400 m ASL, the Triassic dolomite; 3) the Murońka (Mu) mountain (n = 2) located in the Beskid Śląski mountains, 10 km south-west from the town of Żywiec, 950 m ASL, the Carpathian flysch; 4) the Kletno (K) village (n = 2) located in the Masyw Śnieżnika mountains, 7 km south-west from the town of Stronie Śląskie, 830 m ASL, the Devonian marble; 5) the Marzysz (M) village (n = 5) located in the Świętokrzyskie Mountains, 9 km south-east from the city of Kielce, 240 m ASL, the Devonian limestone.

 

Analysis of element concentrations in the soil and plant tissues

Soil samples were sieved using a 2 mm mesh and dried at 105°C. Total As, Cd, Pb, Tl and Zn were extracted by digestion of the ground soil with hot concentrated HClO4 (Foss Tecator Digestor Auto). Asplenium viride specimens from the metalliferous site were washed carefully in running tap water, followed by double-distilled and deionized water. The plants were divided into shoots and roots and dried at 80°C for 48 h. To analyze As, Cd, Pb, Tl and Zn content, the plant material was ground and digested in a hot concentrated mixture of HNO3 and HClO4 (4:1; Foss Tecator Digestor Auto). The elemental contents of the soil and plant extracts were analyzed using flame or graphite furnace atomic absorption spectrometry (Varian AA280FS; Varian AA280Z, GTA 120). Certified reference materials were used to estimate the quality of the metal analyses of soil (CRM048-050; RTC) and plants (Oriental Basma Tobacco Leaves, INCT-OBTL-5; The Institute of Nuclear Chemistry and Technology and moss Pleurozium schreberi M2; The Finnish Forest Research Institute).

 

Plant material

Ploidy level and genome size were estimated in fresh and young leaves of A. viride, collected from metalliferous (HM) and non-metalliferous (NM) sites. Ploidy was determined using the diploid plants of A. viride as a reference standard. Leaves of Vicia faba ‘Inovec’ (2C=26.90 pg; Doležel et al. 1992) were used as an internal standard for nuclear DNA content estimation. Molecular analyses were performed on 24 specimens originating from one metalliferous and five non-metalliferous sites.

 

Flow cytometry measurements

Samples for cytometry analysis were prepared according to Galbraith et al. (1983), with some modifications. Plant tissues of each sample were chopped simultaneously with a sharp razor blade in a plastic Petri dish, with 1 ml of TrisMgCl(200 mM Tris, 4mM MgCl2∙6H2O, 0.5% (v/v) Triton X-100, pH 7.5) nucleus-isolation buffer, supplemented with fluorochrome 4,6'-diamidino-2-phenylindole (DAPI, 2μg/mL) and 2 % (v/v) polyvinylpyrrolidone (PVP-10) for ploidy level, or propidium iodide (PI, 50 μg/mL), ribonuclease A (50 μg/mL) and 2 % PVP-10 for genome size estimation. Ploidy level was estimated by comparison of the position of the G0/G1 peak of the target sample on a histogram with that of the diploid reference standard. Analyses were performed for 5000-7000 nuclei using a Partec CCA flow cytometer (Partec GmbH, Münster, Germany), equipped with a mercury UV lamp. The obtained histograms were analyzed using DPAC v.2.2 software (Partec, GmbH). For nuclear DNA content, 7000 – 10 000 nuclei were measured using a Partec CyFlow SL Green flow cytometer (Partec GmbH, Münster, Germany), equipped with a high-grade solid-state laser with green light emission at 532 nm and side (SSC) and forward (FSC) scatters. Analyses were performed on 15 individual samples from the HM and NM sites. The obtained histograms were evaluated manually using a FloMax software (Partec, GmbH). Genome size was estimated using the linear relationship between the ratios of target species and the internal standard 2C peak positions on the histograms. The coefficient of variation (CV) of the G0/G1 peak of Asplenium samples ranged between 4.45 and 4.49%. The 2C DNA contents (pg) were transformed to megabase pairs of nucleotides using the following conversion: 1 pg = 978 Mbp (Doležel and Bartoš 2005).

 

DNA extraction

Genomic DNA was extracted from silica gel-dried leaf tissue of HM and NM samples using a CTAB extraction protocol (Rogers and Bendich 1994) modified by Trewick et al. (2002) and Kwiatkowska et al. (2019). Leaves with the sori removed underwent homogenization in 500 μl 2% CTAB isolation buffer (2% CTAB, 1.4 M NaCl, 100 mM Tris, 20 mM EDTA), following the addition of 50 μl 10 % sodium dodecyl sulphate buffer (10% SDS, 100 mM Tris-HCl pH 8.0, 20 mM EDTA) and 5 μl 98% β-mercaptoethanol, after which the whole mixture was incubated at 60° C for 1 h. An equal volume of chloroform:isoamyl alcohol mixture 24:1 (v/v) was added, mixed for 10 minutes and centrifuged for 3 minutes at 13000 rpm at room temperature. The collected aqueous phase was mixed with 12.5 μl sodium acetate (pH 5.2) per 100 μl of the mixture. Precipitation was performed by adding a 2/3 volume of isopropanol and incubating at -20° C overnight. DNA pellets were centrifuged at 13000 rpm for 3 minutes at 4° C and rinsed in 500 µl 70% ethanol, centrifuging at 13000 rpm for 1 minute. Vacuum-dried pellets were dissolved in 50 μl nuclease-free water (EURx, Gdańsk, Poland). The quality and quantity of extracted DNA was established using a UV/Vis Q5000 spectrophotometer (Quawell Technology, Inc, CA, USA) and gel electrophoresis on 1% agarose gel. Only samples found to be high quality were used for ISSR-PCR.

 

ISSR analysis

After an initial screening of 15 ISSR primers, five were chosen for further analysis (Tab. 1). ISSR-PCR amplifications were performed in 25 µl reaction volumes containing 2.5 µl 10X DreamTaq™ Green Buffer (Thermo Scientific), 2 µl 10 mM dNTPs (Thermo Scientific), 0.5 µl bovine serum albumin (BSA), 0.25 µl of each primer, 1 µl template DNA and 0.25 µl DreamTaq™ DNA polymerase. PCR reactions were performed as follows: initial denaturation at 94° C for 2 minutes and 35 cycles of denaturation at 94° C for 30 seconds, annealing at 44° C for 45 seconds and elongation at 72° C for 90 seconds. Final extension was performed at 72° C for 20 minutes. Obtained ISSR products were separated through gel electrophoresis on 1% (v/v) agarose gel. The bands were visualized using a GelDoc®-it2 imager (UVP, Jena, Germany), sized using GeneRuler 1kb Plus DNA ladder (Thermo Scientific), and scored by the presence/absence of each feature using PyElph v. 1.4 (Pavel and Vasile 2012).

 

Statistical analysis

Student’s t-test for paired samples was used to compare element concentrations between A. viride shoots and roots. The nuclear DNA content results were analyzed using one-way analysis of variance (ANOVA) and Student’s t-test (STATISTICA v. 10, StatSoft, Poland) to determine possible differences in nuclear DNA content among A. viride plants collected from HM and NM sites. We estimated genetic diversity by calculating parameters both at the population and the whole study area levels. These included: number (P) and proportion (%poly) of polymorphic markers, number of private and discriminating markers (Nprt, Nd; Schlüter 2013), Nei’s gene diversity (Hj; Nei 1973), total gene diversity (HT), mean gene diversity within populations (HS; Nei 1973), gene diversity among populations (GST) and estimated gene flow (Nm; Slatkin and Barton, 1989). All calculations were performed in FAMD v. 1.31 (Schlüter and Harris 2006) and POPGENE v. 1.31 (Yeh et al. 1999). 

A dendrogram, based on Unweighted Pair Group Method with Arithmetic Mean (UPGMA) and NeigborNet, was constructed using Treecon v. 1.3b (Van de Peer and De Wachter 1994) and SplitsTree v. 4.6 (Huson and Bryant 2006), based on a matrix of Nei-Li coefficients (Nei and Li 1979). Bootstrap analysis was performed using the Neighbor-Joining method, with 2000 replicates. Population genetic structure analysis was completed using Bayesian clustering inference in STRUCTURE v. 2.3 (Falush et al. 2007; Pritchard et al. 2000), following the protocol described in Migdałek et al. (2017) and Kwiatkowska et al. (2019). For these calculations, a recessive allele model for dominant markers, admixture and independent allele frequencies between clusters were assumed. Ten independent runs were performed for each K value ranging from 2 to 5, with a burn-in of 200000, followed by 1000000 Markov Chain Monte Carlo replicates. The estimated mean logarithmic likelihoods of K values and ∆K values were calculated to determine the optimal K value (Evanno et al. 2005) using Structure Harvester v. 0.6 (Earl and VonHoldt 2012). STRUCTURE results were summarized in CLUMPAK (Kopelman et al. 2015), using the LargeKGreedy search method and 2000 random input repeats. A three-level hierarchical analysis of molecular variance (AMOVA) was performed in ARLEQUIN v. 3.5 (Excoffier et al. 2005), aiming on testing the statistical significance of inferred groups. Calculations were run on all individuals and population groups suggested by STRUCTURE using a pairwise difference distance matrix at P = 0.05.

Results And Discussion

Soil contamination measurements

The genus Asplenium comprises about 50 species in Europe, half of which are diploid species. Most of Europe, including Scandinavia and the British Isles, is occupied by polyploid species, while the range of diploid species extends mainly in the south of Europe and along the Mediterranean (Vogel et al. 1999a). The genus has a significant potential for tolerance of trace elements, including heavy metals. Our study showed that concentrations of As, Cd, Pb and Zn at the metalliferous site were very high. Concentrations in soils collected from heaps ranged from 11 to 237 mg As kg-1, from 25 to 213 mg Cd kg-1, from 776 to 5437 mg Pb kg-1, from 0.40 to 2.36 mg Tl kg-1 and from 3548 to 44659 mg Zn kg-1. The maximum concentrations of these elements in soil were approximately 5 (As), 20 (Cd), 10 (Pb) and 45 (Zn) times higher than allowable limits established by Polish law for forest soils (Minister of the Environment, 2016). In contrast, the Tl concentrations were closer to the values measured typically in unpolluted soils (Kabata-Pendias and Pendias 2011).

The concentrations of most elements in A. viride samples from the metalliferous site were elevated in both shoots and roots, with quite extreme values in the latter, when compared to the “reference plant” introduced by Markert (1992, Tab. 2). The maximum concentrations in the roots reached 0.5 mg kg-1 As, 56.5 mg kg-1 Cd, 3241 mg kg-1 Pb, 2.3 mg kg-1 Tl and 5124 mg kg-1 Zn (Tab. 2). Maximum element concentrations measured in A. viride roots were 5 (As), 1130 (Cd), 3241 (Pb), 45 (Tl) and 102 (Zn) times higher than in the “reference plant” (Markert 1992). 

Our results indicate that A. viride represents the “excluder” strategy for survival on soils contaminated with heavy metals, which has also been observed by other authors (Punz et al. 1994; Punz and Sieghardt 1993). If the concentration of heavy metals is increased in the above-ground organs, this is typically due to air pollution (Mróz and Rudecki 2009). Some A. viride variants can exhibit high bioaccumulation coefficients during in vitro treatment with Pb (Soare et al. 2015), while sub-tropical species (e.g., A. achilleifolium and A. nidus L.) manifest significant tolerance to arsenic (As) (Chang et al. 2009; Liu et al. 2006). In addition to tolerance, the leaves of A. nidus can be used as a biosorbent for the removal of Pb and Ni from aqueous environments (Dissanayake et al. 2016). 

 

Ploidy level study

All analyzed samples of A. viride collected from HM and NM sites had the same ploidy level as the diploid reference standard (Fig. 1). The genome size of the studied plants of A. viride was very similar across both (HM and NM) sites and ranged from 8.67 pg/2C to 8.69 pg/2C, which correspond to 8479 and 8499 Mbp, respectively. Statistical analysis indicated no significant differences in 2C-values between the two studied habitats. The obtained histograms of genome size, along with the internal standard (V. faba ‘Inovec’), showed two distinct G0/G1 peaks (Fig. 2). Flow cytometry measurements revealed that all studied plants possessed intermediate genome sizes according to Soltis et al.’s (2003) classification. The mean 2C DNA content estimated for this species was 8.68 pg/2C. To the best of our knowledge, this is the first report indicating the genome size of A. viride

 

Genetic diversity studies

Samples from the HM and NM groups (24 individuals), were screened using five ISSR primers, which resulted in the amplification of 48 clearly resolved bands (6-14 per primer), with 37 polymorphic loci. The percentage of polymorphic loci was in excess of 77%. There was a low overall total genetic diversity (HT = 0.25). The within-population diversity (HS = 0.09) was also much lower than diversity among populations (GST = 0.62). The highest within population diversity was present in the M population (H = 0.17), while the lowest was found in the K population (H = 0.02). The metalliferous group showed similar or lower diversity (H = 0.15) than some of the natural (NM) populations. 

The properties of the substrates on which certain populations grow are known to cause morphological variability within Asplenium species, which is maintained in later generations. In case of A. adiantum-nigrum, populations growing on calciferous substrate are morphologically distinct from serpentine populations, and often fail to germinate on serpentine soil (Sleep, 1985). A. viride serpentine populations have been proved to be genetically distinct, and are probably more primitive than those colonizing limestone soils (James et al. 2008).

The UPGMA clustering and NeighborNet analysis showed a lack of diverging population groups (using bootstrap > 50) (Figs. 3, 4 The Bayesian STRUCTURE analysis revealed K = 2 as the optimal number of groups. One main, larger group comprised most of the populations (CS, C, G and Mu), while the remaining two populations (K and M) clustered in a second smaller group. Within the M population, some of the individuals showed admixture of both genetic groups (Fig. 5). Analysis of molecular variance (AMOVA) showed the highest percentage of variation within populations (55.22%) and the lowest percentage of variation between groups (14.76%). The genetic variation between populations relative to the total variance was also the highest among the studied populations (FST = 0.45), while the variation between groups was the lowest (FCT = 0.15) and was not statistically significant (P = 0.072; Tab. 3). 

Despite some apparent morphological traits of plant populations in heavy metal contaminated areas, notably dichotomic fronds, no evident genetic differences were found between metalliferous and non-metalliferous plants in this study; Heavy metal populations did not form a separate, highly supported clade. The studied populations did not exhibit significant differences in genetic structure when studied with Bayesian analysis. Populations from ore-bearing sites clustered together with non-metalliferous populations, which does not support the presence of a heavy metal ecotype. 

All populations exhibited rather high between-population gene diversity, which could suggest some influence of intragametophytic selfing on genetic variability. Selfing as a cause of low within-population genetic variability is particularly apparent in polyploid fern species, which are less vulnerable to the effects of inbreeding depression (Schneller and Holderegger 1996; Vogel et al. 1999b). Diploid species generally favor an outcrossing mating system, however selfing can occur due to the patchy character of their habitats (de Groot et al. 2012a; Vogel et al. 1999a). Long-distance colonization from single spores can also cause founder effects, that can increase variation between populations (de Groot et al. 2012a; b). The consequences of a prior founder effect might be responsible for the lower variability observed within the heavy metal population group.

Conclusions

Our studies have shown that A. viride displays high overall tolerance to heavy metals, exhibiting an ‘excluder’ strategy. However, no evidence was found that increased heavy metal concentration has a significant impact on genetic diversity and on the formation of new ecotypes.

Declarations

Conflict of interest

The authors declare no competing of interests.

 

Acknowledgements 

The research was funded by the Institute of Biology of the Nicolaus Copernicus University in Toruń, the Institute of Biology of the Pedagogical University of Kraków, the Laboratory of Molecular Biology and Cytometry of UTP University of Science and Technology in Bydgoszcz and the Władysław Szafer Institute of Botany of the Polish Academy of Sciences. 

References

  1. Bennett MD, Leitch IJ (2012) Plant DNA C-values Database (release 6.0, December 2012). https://data.kew.org/cvalues/
  2. Brownsey P (1977) A taxonomic revision of the New Zealand species of Asplenium. N. Z. J. Bot., 15: 39–86. https//doi:10.1080/0028825X.1977.10429618
  3. Chang JS, Yoon IH, Kim KW (2009) Heavy metal and arsenic accumulating fern species as potential ecological indicators in As-contaminated abandoned mines. Ecol. Indic., 9: 1275–1279. https//doi:10.1016/j.ecolind.2009.03.011
  4. Coulaud J, Barghi N, Lefèbvre C, Siljak-Yakovlev S (1999) Cytogenetic variation in populations of Armeria maritima (Mill.) Willd. in relation to geographical distribution and soil stress tolerances. Can. J. Bot., 77: 673–685. https//doi:10.1139/cjb-77-5-673
  5. Czapik R (2002) Embryological and cytological variability of plants in polluted environment. Polish Botanical Studies (Vol. 15). Kraków: Szafer Institute of Botany Polish Academy of Sciences.
  6. DalCorso G (2012) Heavy Metal Toxicity in Plants. Pp. 1–25 in Furini A. (ed.), Plants and Heavy Metals. Dordrecht: Springer Netherlands.
  7. de Groot GA, During HJ, Ansell SW, Schneider H, Bremer P, Wubs ERJ (2012a) Diverse spore rains and limited local exchange shape fern genetic diversity in a recently created habitat colonized by long-distance dispersal. Ann. Bot., 109: 965–978. https//doi:10.1093/aob/mcs013
  8. de Groot G Arjen, Verduyn B, Wubs EJ, Erkens RH, During HJ (2012b) Inter-and intraspecific variation in fern mating systems after long-distance colonization: the importance of selfing. BMC Plant. Biol., 12: 3. https//doi:10.1186/1471-2229-12-3
  9. Dissanayake DMREA., Wijesinghe WMKEH, Iqbal SS, Priyantha N, Iqbal MCM (2016) Isotherm and kinetic study on Ni(II) and Pb(II) biosorption by the fern Asplenium nidus L. Ecol. Eng., 88: 237–241. https//doi:10.1016/j.ecoleng.2015.12.028
  10. Doležel J, Bartoš J (2005) Plant DNA flow cytometry and estimation of nuclear genome size. Ann. Bot., 95: 99–110. https//doi:10.1093/aob/mci005
  11. Doležel J, Sgorbati S, Lucretti S (1992) Comparison of three DNA fluorochromes for flow cytometric estimation of nuclear DNA content in plants. Physiol. Plant., 85: 625–631. https//doi:10.1111/j.1399-3054.1992.tb04764.x
  12. Dyer RJ, Savolainen V, Schneider H (2012) Apomixis and reticulate evolution in the Asplenium monanthes fern complex. Ann. Bot., 110: 1515–1529. https//doi:10.1093/aob/mcs202
  13. Earl DA, VonHoldt BM (2012) STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour., 4: 359–361. https//doi:10.1007/s12686-011-9548-7
  14. Evanno G, Regnaut S, Goudet J (2005) Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol. Ecol., 14: 2611–2620. https//doi:10.1111/j.1365-294X.2005.02553.x
  15. Excoffier L, Laval G, Schneider S (2005) Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evolutionary Bioinformatics Online, 1: 47–50.
  16. Falush D, Stephens M, Pritchard JK (2007) Inference of population structure using multilocus genotype data: dominant markers and null alleles. Mol. Ecol. Notes, 7: 574–578. https//doi:10.1111/j.1471-8286.2007.01758.x
  17. Galbraith DW, Harkins KR, Maddox JM, Ayres NM, Sharma DP, Firoozabady E (1983) Rapid flow cytometric analysis of the cell cycle in intact plant tissues. Science, 220: 1049–1051. https//doi:10.1126/science.220.4601.1049
  18. Hildebrandt U, Hoef-Emden K, Backhausen S, Bothe H, Bożek M, Siuta, Kuta E (2006) The rare, endemic zinc violets of Central Europe originate from Viola lutea Huds. Plant. Sys. Evol., 257: 205–222. https//doi:10.1007/s00606-005-0387-4
  19. Huson DH, Bryant D (2006) Application of phylogenetic networks in evolutionary studies. Mol. Biol. Evol., 23: 254–267. https//doi:10.1093/molbev/msj030
  20. Ivanova D, Piȩkoś-Mirkowa H (2003) Chromosome numbers of Polish ferns. Acta. Biol. Crac. Ser. Bot., 45: 93–99.
  21. James KE, Schneider H Ansell SW, Evers M, Robba L, Uszynski G (2008) Diversity Arrays Technology (DArT) for pan-genomic evolutionary studies of non-model organisms. PLoS ONE, 3: e1682. https//doi:10.1371/journal.pone.0001682
  22. Johnston WR, Proctor J (1979) Ecological studies on the Lime Hill serpentine, Scotland. Transactions of the Botanical Society of Edinburgh, 43: 145–150. https//doi:10.1080/03746607908685347
  23. Kabata-Pendias A, Pendias H (2011) Trace elements in soils and plants (4th ed.). Boca Raton London New York: Taylor and Francis Group.
  24. Koivunen S, Saikkonen K, Vuorisalo T, Mutikainen P (2004) Heavy metals modify costs of reproduction and clonal growth in the stoloniferous herb Potentilla anserina. Evol. Ecol., 18: 541–561. https//doi:10.1007/s10682-004-5143-7
  25. Kopelman NM, Mayzel J, Jakobsson M, Rosenberg NA, Mayrose I (2015) Clumpak: A program for identifying clustering modes and packaging population structure inferences across K. Mol. Ecol. Resour., 15: 1179–1191. https//doi:10.1111/1755-0998.12387
  26. Kramer KU, Viane RLL (1990) Aspleniaceae. Pp 52–57 in Kramer, K.U. and Green, P.S. (eds.), Pteridophytes and Gymnosperms. Berlin: Springer Berlin Heidelberg. https//doi:10.1007/978-3-662-02604-5_14
  27. Kwiatkowska M, Zabicka J, Migdalek G, Zabicki P, Cubała M, Bohdanowicz J (2019) Comprehensive characteristics and genetic diversity of the endemic Australian Viola banksii (section Erpetion, Violaceae). Aust. J. Bot., 67: 81–98. https//doi:10.1071/BT18233
  28. Lepp NW (2001) Bryophytes and Pteridophytes. In MNV Prasad (Ed.), Metals in the Environment: Analysis by Biodiversity (pp. 159–170). New York, Basel: Marcel Dekker, Inc.
  29. Liu WJ, Kim KW, Zhu YG, Lee SW, Chang PC Kwak JH (2006) A survey of arsenic and other heavy metals in vegetation from markets or mine tailings. J. Environ. Sci. (China), 18: 287–291.
  30. Marin PD, Tatić B (2001) Serpentine soil and plant diversity, with emphasis balkan Peninsula. Bocconea, 13, 145–150.
  31. Markert B (1992) Establishing of “Reference Plant” for inorganic characterization of different plant species by chemical fingerprinting. Water Air Soil Pollut., 64: 533–538. https//doi:10.1007/BF00483363
  32. Meyer CL, Kostecka AA, Saumitou-Laprade P, Créach A, Castric V, Pauwels M, Frérot H (2010) Variability of zinc tolerance among and within populations of the pseudometallophyte species Arabidopsis halleri and possible role of directional selection. New Phytol., 185(1), 130–142. https//doi:10.1111/j.1469-8137.2009.03062.x
  33. Migdałek G, Nowak J, Saługa M, Cieślak E, Szczepaniak M, Ronikier M (2017) No evidence of contemporary interploidy gene flow between the closely related European woodland violets Viola reichenbachiana and V. riviniana (sect. Viola, Violaceae). Plant Biol., 19: 542–551. https//doi:10.1111/plb.12571
  34. Minister of the Environment, 2016. Regulation of the Minister of the Environment of 1 September 2016 on the Method of the Contamination Assessment of the Earth Surface. In Journal of Laws of 2016; Item 1395, Wydawnictwo Sejmowe: Warsaw, Poland (In Polish).
  35. Mróz L, Rudecki AL (2009) Comparison of chemical composition of two co-occurring chasmophytes, Asplenium ruta-muraria L. (Pteridophyta) and Cymbalaria muralis Gaer., Mey. and Scher. (Spermatophyta) from habitats differing in air pollution. Pol. J. Environ. Stud., 18: 853–863.
  36. Nei M (1973) Analysis of gene diversity in subdivided populations. Proceedings of the National Academy of Sciences, 70: 3321–3323.
  37. Nei M, Li WH (1979) Mathematical model for studying genetic variation in terms of restriction endonucleases. Proceedings of the National Academy of Sciences, 76: 5269–5273.
  38. Pavel AB, Vasile CI (2012) PyElph - a software tool for gel images analysis and phylogenetics. BMC Bioinformatics, 13: 9. https//doi:10.1186/1471-2105-13-9
  39. Perrie LR, Shepherd LD, De Lange PJ, Brownsey PJ (2010) Parallel polyploid speciation: Distinct sympatric gene-pools of recurrently derived allo-octoploid Asplenium ferns. Molecular Ecology, 19: 2916–2932. https//doi:10.1111/j.1365-294X.2010.04705.x
  40. Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics, 155: 945–959.
  41. Proctor J, Woodell SRJ (1971) The Plant Ecology of Serpentine: I. Serpentine Vegetation of England and Scotland. Journal of Ecology, 59: 375–395.
  42. Przedpełska E, Wierzbicka M (2007) Arabidopsis arenosa (Brassicaceae) from a lead-zinc waste heap in southern Poland - A plant with high tolerance to heavy metals. Plant and Soil, 299: 43–53. https//doi:10.1007/s11104-007-9359-5
  43. Punz WF, Sieghardt H (1993) The response of roots of herbaceous plant species to heavy metals. Environmental and Experimental Botany, 33: 85–98. https//doi:10.1016/0098-8472(93)90058-N
  44. Punz W, Kovacs G, Körber-Ulrich SM, Thonke A, Wieländer B, Wieshofer I (1994) Schwermetallstandorte im mittleren Alpenraum und ihre Vegetation - neue Befunde. Verhandlungen der Zoologisch-Botanischen Gesellschaft Österreich, 131: 1–26.
  45. Rogers SO, Bendich AJ (1994) Extraction of total cellular DNA from plants, algae and fungi. Plant Molecular Biology: 183–190. https//doi:10.1007/978-94-011-0511-8_12
  46. Roux JP (1993) The genus Asplenium L. (Aspleniaceae: Pteridophyta) in the Tristan-Gough Island Group. Kew Bulletin, 48: 79–97. https//doi:10.2307/4115750
  47. Schlüter PM (2013) FAMD - Fingerprint Analysis with Missing Data 1.31 Manual. Zurich: University of Zurich.
  48. Schlüter PM, Harris SA (2006) Analysis of multilocus fingerprinting data sets containing missing data. Molecular Ecology Notes, 6: 569–572. https//doi:10.1111/j.1471-8286.2006.01225.x
  49. Schneider H, Navarro-Gomez A, Russell SJ, Ansell S, Grundmann M, Vogel J (2013) Exploring the utility of three nuclear regions to reconstruct reticulate evolution in the fern genus Asplenium. Journal of Systematics and Evolution, 51: 142–153. https//doi:10.1111/j.1759-6831.2012.00226.x
  50. Schneller JJ, Holderegger R (1996) Genetic variation in small, isolated fern populations. Journal of Vegetation Science, 7: 113–120. https//doi:10.2307/3236423
  51. Slatkin M, Barton NH (1989) A comparison of three indirect methods for estimating average levels of gene flow. Evolution, 43(7), 1349–1368. https//doi:10.2307/2409452
  52. Sleep A (1985) Speciation in relation to edaphic factors in the Asplenium adiantum-nigrum group. Proceedings of the Royal Society of Edinburgh. Section B. Biological Sciences, 86: 325–334. https//doi:10.1017/S0269727000008290
  53. Słomka A, Sutkowska A, Szczepaniak M, Malec P, Mitka J, Kuta E (2011) Increased genetic diversity of Viola tricolor L. (Violaceae) in metal-polluted environments. Chemosphere, 83: 435–442. https//doi:10.1016/j.chemosphere.2010.12.081
  54. Smith AR, Pryer KM, Schuettpelz E, Korall P, Schneider H, Wolf PG (2006) A classification for extant ferns. Taxon, 55: 705–731. https//doi:10.2307/25065646
  55. Soare LC, Visoiu E, Bejan C, Dobrescu CM, Fierascu I, Iosub I, Paunescu A (2015) Research on the in vitro bioaccumulation capacity of lead in some pteridophyte species of the Romanian flora. Revista de Chimie, 66: 2017–2020.
  56. Soltis DE, Soltis PS, Bennett MD, Leitch IJ (2003) Evolution of genome size in the angiosperms. American Journal of Botany, 90: 1596–1603. https//doi:10.3732/ajb.90.11.1596
  57. Temsch EM, Temsch W, Ehrendorfer-Schratt L, Greilhuber J (2010) Heavy metal pollution, selection, and genome size: the species of the žerjav study revisited with flow cytometry. Journal of Botany, 2010: 1–11. https//doi:10.1155/2010/596542
  58. Trewick S, Morgan-Richards M, Russell SJ, Henderson S, Rumsey FJ, Pintér I (2002) Polyploidy, phylogeography and Pleistocene refugia of the rockfern Asplenium ceterach: Evidence from chloroplast DNA. Molecular Ecology, 11: 2003–2012. https//doi:10.1046/j.1365-294X.2002.01583.x
  59. Valentine DH, Moore DM (1993) Aspleniaceae. Pp. 18–23 in Tutin, T.G.N.A., Burges, A.O., Chater, J.R., Edmondson, V.H., Heywood, Moore, D.M., Valentine, D.H., Walters, S.M. and Webb, D.A. (eds.), Flora Europaea 1. Cambridge: Cambridge University Press.
  60. Van de Peer Y, De Wachter R (1994) TREECON for Windows: a software package for the construction and drawing of evolutionary trees for the Microsoft Windows environment. Computer applications in the biosciences: CABIOS, 10: 569–570.
  61. Vogel JC, Rumsey FJ, Schneller JJ, Barret JA, Gibby M (1999a) Where are the glacial refugia in Europe? Evidence from pteridophytes. Biological Journal of the Linnean Society, 66: 23–37. https//doi:10.1111/j.1095-8312.1999.tb01915.x
  62. Vogel JC, Rumsey FJ, Russell SJ, Cox CJ, Holmes JS, Bujnoch W (1999b) Genetic structure, reproductive biology and ecology of isolated populations of Asplenium csikii (Aspleniaceae, pteridophyta). Heredity, 83: 604–612. https//doi:10.1038/sj.hdy.6886120
  63. Werth CR, Guttman SI, Eshbaugh WH (1985) Recurring origins of allopolyploid species in Asplenium. Science, 228: 731–733. https//doi:10.1126/science.228.4700.731
  64. Woch M, Stefanowicz AM, Stanek M (2017) Waste heaps left by historical Zn-Pb ore mining are hotspots of species diversity of beech forest understory vegetation. Science of the Total Environment, 599–600: 32–41. https//doi:10.1016/j.scitotenv.2017.04.197
  65. Woch MW (2015) Characteristics of landscape features related tomining and metallurgy in the Olkusz region. Pp. 43–53, in Godzik, B. (ed.), Natural and historical values of the Olkusz Ore-bearing region Kraków: W. Szafer Institute of Botany, Polish Academy of Sciences.
  66. Yeh FC, Yang RC, Boyle T (1999) PopGene Version 131: Microsoft Windows-based freeware for population genetic analysis. University of Alberta and Centre for International Forestry Research, 11–23.
  67. Zając A, Zając M (2001) Distribution atlas of vascular plant in Poland. Kraków: Laboratory of Computer Corology Institute of Botany Jagiellonian University.

Tables

Table 1. Sequences of primers used in ISSR analysis, along with total number of bands (Nbands) and number of polymorphic bands (Npoly) generated by each primer.

Name 

Primer sequence (5’-3’)

Nbands

 

Npoly

 

 

 

HB10

GAGAGAGAGAGACC

10

 

8

 

 

 

HB11

GTGTGTGTGTGTCC

11

 

6

 

 

 

HB12

CACCACCACGC

14

 

12

 

 

 

HB14

CTCCTCCTCGC

6

 

6

 

 

 

HB15

GTGGTGGTGGC

7

 

5

 

 

 

 

 

Table 2. Element concentrations (mg kg-1) in the tissues of Asplenium viride from metalliferrous site (N=5).

 

Shoot

 

 

Root

 

 

t-test

Minimum

Mean (SD)

Maximum

Minimum

Mean (SD)

Maximum

p

As 

0.004

0.023 (0.018)

0.045

0.098

0.244 (0.171)

0.504

0.005

Cd

1.1

2.1 (0.9)

3.5

10.2

34.7 (20.8)

56.5

<0.001

Pb

7

33 (50)

122

297

1474 (1341)

3241

<0.001

Tl

0.04

0.09 (0.05)

0.18

0.97

1.50 (0.49)

2.27

<0.001

Zn

117

200 (110)

376

1063

2992 (1641)

5124

<0.001

Student’s t-test for paired samples was used to compare element concentrations between shoots and roots.

 

Table 3. Result parameters of the AMOVA performed on 6 populations of Asplenium viride grouped according to STRUCTURE results.

Source of variation

df

Sum of squares

Variance components

Percentage of variation

Fixation Indices

Among groups

1

23.24

1.08

14.76

FCT = 0.15 P = 0.072

Among populations within groups

4

43.62

2.19

30.02

FSC = 0.35*

Within populations

15

60.38

4.03

55.22

FST = 0.45*

Total

20

127.24

7.29

-









d.f. – degrees of freedom, FCT – variance among groups relative to total variance, FSC – variance among populations within groups,  FST – variance among populations relative to total variance, P – significance level (if F values marked by *, then < 0.001).