Attila I. Engloner
Kitti Németh1
Emese Meglécz- 1Institute of Aquatic Ecology, Centre for Ecological Research, Budapest, Hungary
- 2Institute of Plant Biology, Biological Research Centre, Eötvös Loránd Research Network, Szeged, Hungary
- 3Department of Biotechnology, Faculty of Science and Informatics, Szeged University, Szeged, Hungary
- 4Aix Marseille University, Avignon University, French National Center for Scientific Research (CNRS), French National Research Institute for Sustainable Development (IRD), Mediterranean Institute of Marine and Terrestrial Biodiversity and Ecology (IMBE), Marseille, France
- 5Molecular Taxonomy Laboratory, Hungarian Natural History Museum, Budapest, Hungary
The adaptability of plant populations to a changing environment depends on their genetic diversity, which in turn is influenced by the degree of sexual reproduction and gene flow from distant areas. Aquatic macrophytes can reproduce both sexually and asexually, and their reproductive fragments are spread in various ways (e.g. by water). Although these plants are obviously exposed to hydrological changes, the degree of vulnerability may depend on the types of their reproduction and distribution, as well as the hydrological differences of habitats. The aim of this study was to investigate the genetic diversity of the cosmopolitan macrophyte Ceratophyllum demersum in hydrologically different aquatic habitats, i.e. rivers and backwaters separated from the main river bed to a different extent. For this purpose, the first microsatellite primer set was developed for this species. Using 10 developed primer pairs, a high level of genetic variation was explored in C. demersum populations. Overall, more than 80% of the loci were found to be polymorphic, a total of 46 different multilocus genotypes and 18 private alleles were detected in the 63 individuals examined. The results demonstrated that microsatellite polymorphism in this species depends on habitat hydrology. The greatest genetic variability was revealed in populations of rivers, where flowing water provides constant longitudinal connections with distant habitats. The populations of the hydrologically isolated backwaters showed the lowest microsatellite polymorphism, while plants from an oxbow occasionally flooded by the main river had medium genetic diversity. The results highlight that in contrast to species that spread independently of water flow or among hydrologically isolated water bodies, macrophytes with exclusive or dominant hydrochory may be most severely affected by habitat fragmentation, for example due to climate change.
1 Introduction
Aquatic macrophytes are widely distributed and of great importance in various freshwater and marine ecosystems. They play a key role in aquatic food webs and biochemical cycles, influence the hydrology and sediment dynamics of aquatic habitats by altering current velocity, and provide habitat and refuge for many living organisms from microbiota to vertebrates and multiple benefits for humans (Bornette and Puijalon, 2011; Hossain et al., 2017; Thomaz, 2023). It is widely recognized that climate change threatens aquatic plants in several ways, including rising temperature, CO2 concentration and dissolved organic carbon, as well as changes in nutrient availability, light conditions and salinity. Climate change has been proven to greatly affect the growth and physiology of these plants, as well as the species composition and ecosystem functioning (Hossain et al., 2017; Li et al., 2017; Reitsema et al., 2018; Reitsema et al., 2020). First, aquatic macrophytes are completely dependent on the water supply of their habitat and are obviously exposed to hydrological changes caused by global climate change. It is clear that where less precipitation falls and the water supply of the habitats decreases, the living conditions of macrophytes deteriorate, or, if the climate becomes wetter, their potential habitats may increase. In river ecosystems, altered flow regimes have significant effects on macrophytes (Reitsema et al., 2020; Goldenberg-Vilar et al., 2022; Rivaes et al., 2022). Droughts and floods associated with climate change have both been shown to simplify physical habitats by leading to shallow flow conditions and reduced hydraulic diversity, especially in ecologically important pool and riffle features(O’Briain, 2019).
Not only the amount, presence or absence of water in a given habitat can affect the colonization of aquatic macrophytes, but also the degree of connectivity between different habitats can be of great importance. These plants can reproduce both sexually and asexually, and their seeds and vegetative shoot fragments can spread in several ways; by wind, birds or water (Santamaría, 2002; Capers, 2003; Jones et al., 2020). It is easy to understand, that dispersion by water flow (hydrochory) is not possible between hydrologically isolated habitats. In addition to sexual reproduction, gene inflow from distant areas can significantly increase the genetic diversity of populations, which can improve their ability to adapt to changing environments (Cao et al., 2020; Cao et al., 2021). Conversely, if a plant species reproduces mainly vegetatively and has only limited ‘external’ propagating material introduced into its population (or none at all), it is expected to be more vulnerable to habitat changes.
Changes in the genetic diversity of populations can be caused or enhanced by various environmental factors, among which habitat heterogeneity and dispersal can play decisive role (Orsini et al., 2013; Davis et al., 2018). The environmental factors have also been suggested as important drivers of genetic variation in aquatic plants, however, only a few studies focus on these organisms (Foust et al., 2016; Robertson et al., 2017; Li et al., 2022 and references therein). Temperature, salinity-gradient, heavy metals and bird-mediated dispersal have been shown to affect macrophyte genetic diversity (Gupta and Sarin, 2009; Triest et al., 2010; Cao et al., 2017; Li et al., 2022).
Ceratophyllum demersum is a cosmopolitan submersed macrophyte with wide climatic tolerance and all the ecosystem functions detailed above (Les, 1986; Qadri et al., 2022). This species is known for low sexual reproduction and predominance of clonal growth (Les, 1988; Les, 1991). Although the main reproductive mechanism of C. demersum is reported to be shoot fragmentation (Arber, 1920; Les, 1991; Fukuhara et al., 1997), and dispersion of this species relies on these vegetative fragments flowing in water rather than pollen or seeds (Wade, 1993; Capers, 2003), the effect of hydrological conditions on its genetic variability has never been studied.
Most of the molecular studies on C. demersum are phylogenetic (Moore et al., 2007; Szalontai et al., 2018; Albert and Renner, 2020; Yang et al., 2020), and further investigations concern the genetic alteration in response to heavy metal tolerance (Gupta and Sarin, 2009; Khaleel et al., 2022) and the isolation of genes encoding enzymes potentially required for heavy metal accumulation (Shukla et al., 2012). A variety of molecular markers and techniques have been used for the genetic investigations of Ceratophyllum, from karyology to AFLP, MSAP and ISSR markers (Triest et al., 2010; Cao et al., 2017; Gargiulo et al., 2022; Li et al., 2022), however, microsatellite primers and studies based on them have not yet been published for this species. Microsatellites (i.e. short tandem repeats) have been the most widely used markers for genotyping plants over the past 20 years (Vieira et al., 2016), as they are multi-allelic markers with a mutation rate of 103-106 per cell generation, which can be 10 orders of magnitude greater than point mutations (Gemayel et al., 2012). Their codominant nature, biparental mode of inheritance and elevated levels of polymorphism (Goldstein and Schlötterer, 1999; Schlötterer, 2000; Ellegren, 2004) have made them particularly informative and powerful for investigating genetic diversity and structure as well as demographic processes (Kim and Sappington, 2013). The advantages of using microsatellite markers, especially in population genetic studies, were even demonstrated for macrophyte species (Ouborg et al., 1999; Triest et al., 2010; Kong et al., 2019).
For the above reasons, the aim of the present paper is to explore whether the microsatellite polymorphism of the worldwide macrophyte C. demersum depends on the hydrology of the habitat and how it can affect the genetic diversity of its populations to changing environments. For this purpose, we have developed and provide the first microsatellite primers set for this species.
2 Materials and methods
2.1 Marker development
To produce high-quality DNA extract, C. demersum samples were first rinsed and impurities were removed with toothbrushes. The stems were then discarded and the leaves were used for subsequent processes. About 80 g leaf material were grinded in 200 ml of “grinding buffer” (450 mM Sucrose, 1.5 mM EGTA, 0.2% [w/v] bovine serum albumin [BSA], 0.6% [w/v] polyvinylpyrrolidone 40, 10 mM dithiothreitol [DTT], 0.2 mM phenylmethylsulfonyl fluoride [PMSF], and 15 mM MOPS [3-(N-morpholino)-propanesulfonic acid]/KOH, pH 7.4). The cells were disrupted by homogenizing for three periods of 15 s using a Waring blender. Cell debris was removed by centrifugation at 3,000 g for 5 min. According to the recommendations, the sample DNA was isolated using Thermo Scientific™ Genomic DNA Purification Kit (Thermo Fisher Scientific, Waltham, MA, USA) and sequenced on Illumina Miseq sequencer using v3 kit for 2x300 bp PE protocol. The sequences were aligned using Mira4 (Chevreux et al., 1999) in accurate mode and the contigs exceeding 1 kb length were used for microsatellite prediction.
The potential SSR loci were mined from the assembled, total of 289 contigs, by QDD v. 3.1.2 (Meglécz et al., 2014) using default parameters. The following criteria were used to select primer pairs for laboratory testing: (i) the sequence contained only pure microsatellites in the target region with at least six repeats; (ii) did not contain repeats of (AT)n or (CG)n; (iii) the primer alignment score to the amplified sequence was lower than six; (iv) primers were at least five bases away from the microsatellite motif; (v) no BLAST hit to non-Viridiplantae sequences in GenBank; and only one microsatellite was selected from each contig to avoid linkage disequilibrium. After the above selection, a total of 50 primer pairs were selected for initial screening.
Among the potential primers tested according to the protocol detailed in the next section, 16 consistently amplified the target sequences and showed signs of variability between the tested individuals based on visual inspection of the agarose gel. Of the 16 loci, those with at least two alleles were retained for the final marker set.
Finally, our marker development resulted in 10 recommendable primer pairs to explore genetic polymorphism of C. demersum (Table 1). The use of these primers is facilitated by the same optimal annealing temperature and a wide range of the amplicon size, which allows multiplexing of PCR products in a single reaction.
Table 1 Characteristics of 10 microsatellite loci of Ceratophyllum demersum.
2.2 Sampling and laboratory work to explore microsatellite polymorphism of Ceratophyllum populations
Plant material was collected from five aquatic habitats belonging to two main different hydrological types (Figure 1). Two habitats represent rivers (hereinafter referred to as ‘river habitats’), namely (i) the main channel of the River Danube (RD) and (ii) the Soroksári Danube (RS), the second largest side arm of River Danube in Hungary between the 1642 and 1586 river kms. The latter is closed at both ends by sluices regulating its water flow. Due to the hydrological regulation, only a small amount of water enters the side arm from the main channel, it is the fortieth to sixtieth part of the Danube’s water at mean water-flow (Hungarian Hydrological Forecasting Service, 2023). In addition to the two ‘river habitats’, three oxbows were also selected; (iii) Schisler oxbow (OS) and (iv) Zátonyi Danube (OZ) in the Szigetköz, the watery plain with many branches, islands and backwaters between the main Danube channel and the Moson-arm, and (v) the oxbow Mocskos Danube (OM), a former Danube side branch lying on the floodplain of Béda-Karapancsa Landscape Protection Area. Of the five habitats, only the main Danube channel and the Soroksári Danube have water flowing continuously, and only these habitats have steady longitudinal connections with distant river sections. Among the backwaters, the Mocskos Danube can temporarily connect to the Danube during the river’s floods, however, these events have only occurred on a third of the days in the last fifty years (Hungarian Hydrological Forecasting Service, 2023).
Figure 1 Sampling locations along the Danube River. RD – the main Danube channel; RS, Soroksári Danube; OS, Schisler oxbow; OM, Mocskos Danube; OZ, Zátonyi Danube. The numbers after the abbreviations indicate the serial numbers of samples, while numbers in the rectangles show the river kilometers in the main river channel.
The two selected ‘river habitats’ allowed both large-scale and fine-scale sampling: 13 plants were collected between the 1680-425 river kms of the main Danube channel, and 16 samples were taken along the 1600 m long stretch of the Soroksári Danube, at river km 31. Sampling in the oxbows was determined by the size of the water bodies and the amount of plant cover in them. Furthermore, in order to avoid sampling from the same plant (i.e. collecting multiple leaves of the same individual), a distance of at least 30 m was kept between two samples. In this way, the C. demersum stands in the 500, 400 and 900 meter long sections of OS, OZ and OM provided 16, 8 and 10 plant samples for the genetic investigations. Plant material was stored at 4°C during transport to the laboratory, and then at -80°C until DNA isolation.
DNA was extracted by homogenizing 200 mg frozen leaves in 800 μl CTAB isolation buffer (2%) following the protocol described in Bousquet et al. (1990). The amplification procedure from 5 μl of DNA extracts was carried out in 15 μl final reaction volumes containing 10X PCR buffer, 3 mM MgCl2, 0.2 mM dNTPs, 0.05 units/μl of Taq DNA polymerase (Taq DNA Polymerase, recombinant, Fermentas) and 0.5 μM of each primer. Forward primers were labelled fluorescently at their 5’ end. The following cycling conditions were used: initial denaturation for 5 min at 94°C, 40 cycles of 45 s at 94°C, 45 s at 55°C, 45 s at 72°C; final elongation of 10 min at 72°C. The success of PCR was checked by running 2 µl of product on 1.4% agarose gels stained with GelRed Nucleic Acid Stain (Biotium Inc., Fremont, CA, USA). After amplification, PCR products were multiplexed in a single reaction and fragment analysis was carried out on an ABI 3130 Genetic Analyser. Allele sizes were estimated using Peak Scanner software (Thermo Fisher Scientific, Waltham, MA, USA).
2.3 Statistical analyses
Checking for the presence of null alleles was performed using Micro-Checker 2.2.3 (Van Oosterhout et al., 2004) by Monte Carlo simulation of expected homozygote frequencies and heterozygote allele size differences. Since null alleles were only observed in a single case, the entire dataset was used for further analyses. Parameters of polymorphism were determined using GenAlEx v. 6.5 (Peakall and Smouse, 2006; Peakall and Smouse, 2012) and Fstat v. 2.9.4 (Goudet, 2003). Based on the genetic distance matrix calculated by Identix v.1.1 (Belkhir et al., 2002), UPGMA tree was constructed using PAST v. 4.12. A minimum spanning network was generated by the bruvo.msn() function of the ‘poppr’ package v. 2.9.3 in R (Kamvar et al., 2014; Kamvar et al., 2015) using the Bruvo method (Bruvo et al., 2004) to calculate genetic distances. To analyse the correlation between the genetic and geographic distance matrices, a Mantel test was carried out using GenAlEx v. 6.5 (Peakall and Smouse, 2006; Peakall and Smouse, 2012) with 9999 permutations. Genetic distances among individuals were calculated by GenAlEx v. 6.5, while the geographic distance matrix was produced by PAST v. 4.12 (Hammer et al., 2001).
The genetic structure of the populations has been analysed using Bayesian-clustering method (Pritchard et al., 2000). The most probable number of genetically differentiated groups (K) in the populations was estimated and the individuals were assigned to these groups. Structure 2.3.4 was run to carry out these analyses using default settings with an initial burn in of 100,000 steps and running length of 500,000 steps. In the evaluation of the results, ΔK was computed which indicates the change in log probability between successive K values (Evanno et al., 2005). Structure Harvester Web 0.6.94 (Earl and vonHoldt, 2012) was used to compute the ΔK values. The package ‘pophelper’ in R (Francis, 2017) was applied to average the ten runs of the most probable K value given by Structure and correct for label switching.
To reveal the genetic differentiation among plant samples from the five selected habitats, the microsatellite allele frequency data were evaluated by standardized principal component analysis (PCA) using SYN-TAX 2000 computer program package (Podani, 2001).
3 Results
3.1 Level of microsatellite polymorphism
All indices of polymorphism indicated a high level of variation in C. demersum populations (Table 2). Overall, more than 80% of the loci were found to be polymorphic. The average number of alleles per locus was slightly more than 3 and the average frequency of heterozygotes was over 57%. A total of 46 different multilocus genotypes and 18 private alleles were detected in the 63 individuals examined. The variability parameters corrected for sample size, i.e. the number of effective alleles and the allelic richness in total were 2.222 and 2.905, respectively.
Table 2 The parameters of variability based on 10 microsatellite loci.
The data clearly show greater genetic variability in river populations than in backwaters (e.g., Ne 2.567 and 1.992; AR 3.606 and 2.437, for the former and the latter, respectively). The average number of alleles per locus, the number of effective alleles, allelic richness and Shannon’s information index all reached the highest values in the Danube River and the Soroksári side arm, while among the backwaters, most of these values were the highest in the Mocskos oxbow (Table 2).
3.2 Structure of the genetic variation
According to the hierarchical classification, none of the populations formed completely separate clusters, i.e. samples from the different habitats have mixed with each other to some extent (Figure 2). However, with the exception of three samples of Mocskos Danube oxbow, plants of backwater habitats were included in one main clade (marked with a dashed box at the hierarchical level of 0.450 on the UPGMA tree). The largest distance between all oxbow samples was 0.500 (marked with black arrow in Figure 2). Contrary, plants from the two river habitats showed higher mixing and dissimilarity (the latter indicated by an empty arrow at the distance level of 0.650 on the UPGMA tree). Geographically closer samples were not clearly placed in closer clades, which is consistent with the result of Mantel test showing a statistically significant positive but weak correlation between geographic and genetic distance (R=0.249, p=0.0001). It can also be observed in Figure 2, that multilocus genotypes consisting of several plant samples occurred in both river habitats and backwaters (up to six plants per genotype, in both habitat types). Minimum spanning network also shows that plants belonging to the river habitats can appear quite far from each other, and larger circles, i.e. multilocus genotypes with more than one sample, are colored differently, indicating that they were observed in all habitats (Figure 3).
Figure 2 UPGMA tree based on genetic distance matrix among individuals. RD, the main Danube channel; RS, Soroksári Danube; OS, Schisler oxbow; OM, Mocskos Danube; OZ, Zátonyi Danube. Samples from river habitats and oxbows are indicated by bluish and yellowish-red colors, respectively. The dashed box shows the main clade that contains the majority of oxbow plants. Black and empty arrows indicate clades that contain all backwater and river samples, respectively.
Figure 3 Minimum spanning network based on Bruvo’s distance. Each node on the graph represents a different multilocus genotype, while the edges show the genetic distances between them. Areas of circles are proportional to the number of multilocus genotypes displayed. RD, the main Danube channel; RS, Soroksári Danube; OS, Schisler oxbow; OM, Mocskos Danube; OZ, Zátonyi Danube.
According to the Structure analysis, the most probable number of genetically differentiated groups (K) proved to be three (Figure 4). The plants of the main Danube channel consisted predominantly of a single genetic cluster, in addition to which, another genetic cluster appeared in high proportion in the Soroksári Danube. The samples of these two river habitats were clearly separated from those of Schisler and Zátonyi oxbows, which were mainly characterized by the third cluster. In the Mocskos Danube, all three clusters were mixed.
Figure 4 Bayesian assignment of individuals based on 10 microsatellite loci. RD, the main Danube channel; RS, Soroksári Danube; OS, Schisler oxbow; OM, Mocskos Danube; OZ, Zátonyi Danube.
The ordination of microsatellite data resulted in two large, overlapped groups of samples taken from the main Danube River and the Soroksári side arm, and three smaller and better-separated groups representing the oxbow samples (Figure 5). This means that the populations of the three oxbows were separated from each other while the plants of two river habitats were mixed based on microsatellites. As the areas of the polygons enclosing objects in the scattergram are proportional to the variability between the objects (i.e. the larger the area the higher the variability), the result demonstrated that genetic variability of C. demersum in river habitats was higher than in backwaters. Comparing the backwater samples, the plants of the Mocskos oxbow showed the highest variability.
Figure 5 PCA ordination of Ceratophyllum demersum samples based on microsatellite allele frequency data. Convex polygons enclose plant samples collected from the same habitats. RD, the main Danube channel; RS, Soroksári Danube; OS, Schisler oxbow; OM, Mocskos Danube; OZ, Zátonyi Danube.
4 Discussion
Although microsatellites have become the most widely used markers for genotyping plants, and the benefits of their use were demonstrated even in population genetic studies of macrophytes, microsatellite primers have not yet been published for C. demersum. Based on the newly developed loci, this study revealed differences in the microsatellite polymorphism of C. demersum populations related to the hydrological types of the habitats. Although the number of samples obtained from the populations was not exactly the same due to differences in size and vegetation cover of the water bodies, the number of effective alleles and allelic richness make it possible to compare the genetic variability of populations with different sample sizes. Both calculated variability parameters showed greater genetic variability in riverine habitats than in backwaters. In addition, the applied multivariate analyses revealed low genetic variability from backwaters based on both the smaller (Zátonyi and Mocskos oxbows) and larger sample sizes (Schisler oxbow). The results also demonstrated that a high number of multilocus genotypes could be detected from both river and backwater populations, however, the observed genotypes were genetically closer to each other in the latter than in the former habitat.
Genetic variation is highly dependent on the ratio of sexual and asexual reproduction, as well as the short- and long-distance dispersal of sexual and vegetative propagules (Eckert et al., 2016). Asexual spread is often assumed to be important in local, short-distance dispersal while regional, long-distance dispersal is thought to be based on sexual reproduction (Santamaría, 2002). However, the low sexual reproduction and the predominance of clonal growth and vegetative dispersion of C. demersum (Les, 1988; Les, 1991; Wade, 1993; Capers, 2003) suggest that the revealed genetic variability is the result of differences in vegetative propagation. In addition, the clear relationship observed between the genetic diversity of populations and the hydrological types of the habitats indicates that the mode of propagule spread is mainly hydrochory. The greatest genetic variability was revealed in river populations, where continuously flowing water provides a constant longitudinal connection with distant habitats. Contrary, C. demersum stands in the hydrologically isolated backwaters showed the lowest genetic diversity. The population of Mocskos oxbow, occasionally flooded by the Danube, which enables gene flow, had medium genetic diversity. If the majority of C. demersum propagules were not spread by water but by wind or birds, they could easily reach the backwaters that are hydrologically isolated from the main river but geographically very close to it.
Topography and geographic distance affect population genetic structure of plants, especially in aquatic habitats (Kong et al., 2019). However, unlike ponds and isolated backwaters that are discrete in terrestrial landscapes, water flows continuously in rivers and interconnected aquatic habitats, and the longitudinal connection of distant sections is ensured in these habitats. The importance of water flow in the establishment of plants was also demonstrated in the case of other submerged macrophytes, for which, however, seeds and the wind and animal-mediated dispersal were of great importance (Chen et al., 2009; Pollux et al., 2009; Triest et al., 2010). The spread of these macrophytes can be significant even against the water flow (Chen et al., 2009; Jones et al., 2020) or between water bodies with no hydrological connection, potentially providing their populations with high genetic diversity. Conversely, when hydrochory is exclusive or dominant, hydrological connections between habitats are crucial for gene flow, as the results of the present study suggest for C. demersum.
Of course, genetic diversity detected at marker loci in clonally spreading populations may be attributed to somatic mutations, which are enhanced by rapid clonal growth and a large amount of vegetative propagules (de Witte and Stöcklin, 2010; Eckert et al., 2016; Zhu et al., 2017; Kong et al., 2019). However, there are no reports on somatic mutation related to C. demersum.
For the above reasons, although habitat deterioration affects all macrophytes, the loss of hydrological connections between habitats is expected to be more severe for species, such as C. demersum, whose population genetic diversity is highly dependent on flow conditions and hydrological connectivity. The harmful effects of this kind of habitat fragmentation may appear as a result of any local anthropogenic intervention (cf. Cao et al., 2021) or global processes, such as climate change. The effect of the latter is also enhanced by temperature changes, which further reduce the genetic diversity of macrophytes by altering reproductive strategies, resulting in fewer flowers and thus less sexual reproduction (Li et al., 2017). That is, the main processes of climate change can reinforce each other, impairing the adaptability of macrophyte populations and endangering their ecosystem functions.
Data availability statement
The sequence data has been deposited in the SRA database of NCBI with the code PRJNA1012867.
Author contributions
AE: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Visualization, Writing – original draft. KN: Investigation, Visualization, Writing – original draft. PK: Investigation, Methodology, Writing – original draft. EM: Investigation, Methodology, Writing – original draft. JB: Formal Analysis, Investigation, Methodology, Visualization, Writing – original draft.
Funding
This work was partially funded by the National Research, Development and Innovation Office (OTKA K106177).
Acknowledgments
The authors thank Mária Tuschek for the technical assistance.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Albert, V. A., Renner, T. (2020). Aquatic angiosperm ambiguities answered. Nat. Plants 6, 181–183. doi: 10.1038/s41477-020-0607-5
Belkhir, K., Castric, V., Bonhomme, F. (2002). IDENTIX, a software to test for relatedness in a population using permutation methods. Mol. Ecol. Notes 2, 611–614. doi: 10.1046/j.1471-8286.2002.00273.x
Bornette, G., Puijalon, S. (2011). Response of aquatic plants to abiotic factors: a review. Aquat. Sci. 73, 1–14. doi: 10.1007/s00027-010-0162-7
Bousquet, J., Simon, L., Lalonde, M. (1990). DNA amplification from vegetative and sexual tissues of trees using polymerase chain reaction. Can. J. For. Res. 20, 254–257. doi: 10.1139/x90-037
Bruvo, R., Michiels, N. K., D’Souza, T. G., Schulenburg, H. (2004). A simple method for the calculation of microsatellite genotype distances irrespective of ploidy level. Mol. Ecol. 13, 2101–2106. doi: 10.1111/j.1365-294X.2004.02209.x
Cao, Q., Hu, F., Liu, N. (2020). Local-scale patterns of genetic variation in coexisting floating-leaved Nymphoides peltata and submerged Myriophyllum spicatum in Donghu Lake. J. Oceanol. Limnol. 38, 1825–1834. doi: 10.1007/s00343-019-9068-y
Cao, Q., Liu, B., Hu, F. (2021). Effects of hydrological connection and human disturbance on genetic variation of submerged Vallisneria natans populations in four lakes in China. J. Oceanol. Limnol. 39, 1403–1416. doi: 10.1007/s00343-020-0277-1
Cao, Q. J., Mei, F. F., Wang, L. (2017). Population genetic structure in six sympatric and widespread aquatic plants inhabiting diverse lake environments in China. Ecol. Evol. 7, 5713–5723. doi: 10.1002/ece3.3141
Capers, R. S. (2003). Macrophyte colonization in a freshwater tidal wetland (Lyme, CT, USA). Aquat. Bot. 77, 325–338. doi: 10.1016/j.aquabot.2003.08.001
Chen, Y., Li, X., Yin, L., Cheng, Y., Li, W. (2009). Genetic diversity and migration patterns of the aquatic macrophyte Potamogeton malaianus in a potamolacustrine system. Freshw. Biol. 54, 1178–1188. doi: 10.1111/j.1365-2427.2008.02165.x
Chevreux, B., Wetter, T., Suhai, S. (1999). “Genome sequence assembly using trace signals and additional sequence information,” in Computer Science and Biology: Proceedings of the German Conference on Bioinformatics (GCB), vol. 99 . Ed. Wingender, E. (GBF-Braunschweig: Department of Bioinfrormatics), 45–56.
Davis, C. D., Epps, C. W., Flitcroft, R. L., Banks, M. A. (2018). Refining and defining riverscape genetics: How rivers influence population genetic structure. Wiley Interdiscip. Rev. Water 5, e1269. doi: 10.1002/wat2.1269
de Witte, L. C., Stöcklin, J. (2010). Longevity of clonal plants: why it matters and how to measure it. Ann. Bot. 106, 859–870. doi: 10.1093/aob/mcq191
Earl, D. A., vonHoldt, B. M. (2012). STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour. 4, 359–361. doi: 10.1007/s12686-011-9548-7
Eckert, G. C., Dorken, M. E., Barrett, S. C. H. (2016). Ecological and evolutionary consequences of sexual and clonal reproduction in aquatic plants. Aquat. Bot. 135, 46–61. doi: 10.1016/j.aquabot.2016.03.006
Ellegren, H. (2004). Microsatellites: sample sequences with complex evolution. Nat. Rev. Genet. 5, 435–445. doi: 10.1038/nrg1348
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. doi: 10.1111/j.1365-294X.2005.02553.x
Foust, C., Preite, V., Schrey, A. W., Alvarez, M., Robertson, M., Verhoeven, K., et al. (2016). Genetic and epigenetic differences associated with environmental gradients in replicate populations of two salt marsh perennials. Mol. Ecol. 25, 1639–1652. doi: 10.1111/mec.13522
Francis, R. M. (2017). POPHELPER: An R package and web app to analyse and visualize population structure. Mol. Ecol. Resour. 17, 27–32. doi: 10.1111/1755-0998.12509
Fukuhara, H., Tanaka, T., Izumi, M. (1997). Growth and turion formation of Ceratophyllum demersum in a shallow lake in Japan. Jpn. J. Limnol. 58, 335–347. doi: 10.3739/RIKUSUI.58.335
Gargiulo, G. M., El Bakkouri, B., Crisafulli, A., Donato, M., Picone, R. (2022). Polysomaty and chromosome number variation in a population of Ceratophyllum demersum L. from Aquila Lake (Aspromonte Mountains, Calabria, Italy). Aquat. Bot. 180, 103530. doi: 10.1016/j.aquabot.2022.103530
Gemayel, R., Cho, J., Boeynaems, S., Verstrepen, K. J. (2012). Beyond junk-variable tandem repeats as facilitators of rapid evolution of regulatory and coding sequences. Genes 3, 461–480. doi: 10.3390/genes3030461
Goldenberg-Vilar, A., Delgado, C., Penas, F. J., Barquín, J. (2022). The effect of altered flow regimes on aquatic primary producer communities: Diatoms and macrophytes. Ecohydrology 15, e2353. doi: 10.1002/eco.2353
Goldstein, D. B., Schlötterer, C. (1999). Microsatellites: Evolution and Applications (Oxford: Oxford University Press).
Goudet, J. (2003) Fstat (ver. 2.9.4), a program to estimate and test population genetics parameters. Available at: http://www2.unil.ch/popgen/softwares/fstat.htm (Accessed August 11, 2023).
Gupta, M., Sarin, N. B. (2009). Heavy metal induced DNA changes in aquatic macrophytes: Random amplified polymorphic DNA analysis and identification of sequence characterized amplified region marker. J. Environ. Sci. 21, 686–690. doi: 10.1016/S1001-0742(08)62324-4
Hammer, Ø., Harper, D. A. T., Ryan, P. D. (2001). PAST: Paleontological statistics software package for education and data analysis. Palaeontol. Electron. 4, 1–9.
Hossain, K., Yadav, S., Quaik, S., Pant, G., Ismail, M. N. (2017). Vulnerabilities of macrophytes distribution due to climate change. Theor. Appl. Climatol. 129, 1123–1132. doi: 10.1007/s00704-016-1837-3
Hungarian Hydrological Forecasting Service (2023). Available at: https://www.hydroinfo.hu/en/index.html (Accessed August 11, 2023).
Jones, P. E., Consuegra, S., Börger, L., Jones, J., Garcia de Leaniz, C. (2020). Impacts of artificial barriers on the connectivity and dispersal of vascular macrophytes in rivers: A critical review. Freshw. Biol. 65, 1165–1180. doi: 10.1111/fwb.13493
Kamvar, Z. N., Brooks, J. C., Grünwald, N. J. (2015). Novel R tools for analysis of genome-wide population genetic data with emphasis on clonality. Front. Genet. 6. doi: 10.3389/fgene.2015.00208
Kamvar, Z. N., Tabima, J. F., Grünwald, N. J. (2014). Poppr: an R package for genetic analysis of populations with clonal, partially clonal, and/or sexual reproduction. PeerJ 2, e281. doi: 10.7717/peerj.281
Khaleel, H. A., Jawad, H. J., Alrufaye, Z. T. A. (2022). Genotoxic effect of heavy metals on Ceratophyllum demersum L. using RAPD markers. Iran. J. Ichthyol. 9, 46–53.
Kim, K. S., Sappington, T. W. (2013). “Microsatellite data analysis for population genetics,” in Microsatellites: Methods and Protocols. Ed. Kantartzi, S. K. (New York: Humana Press, New York), 271–295.
Kong, F., Wu, Z., Wang, H., Chen, J., Xu, X. (2019). Population genetic structure of the whorl-leaf watermilfoil Myriophyllum verticillatum shaped by topography and geographic distance. Hydrobiologia 838, 55–64. doi: 10.1007/s10750-019-03977-5
Les, D. H. (1986). The phytogeography of Ceratophyllum demersum and C. eChinatum (Ceratophyllaceae) in glaciated North America. Can. J. Bot. 64, 498–509. doi: 10.1139/b86-064
Les, D. H. (1988). The origin and affinities of the Ceratophyllaceae. Taxon 37, 326–345. doi: 10.2307/1222142
Les, D. H. (1991). Genetic diversity in the monoecious hydrophile Ceratophyllum (Ceratophyllaceae). Am. J. Bot. 78, 1070–1082. doi: 10.2307/2444896
Li, Z., He, L., Zhang, H., Urrutia-Cordero, P., Ekvall, M. K., Hollander, J., et al. (2017). Climate warming and heat waves affect reproductive strategies and interactions between submerged macrophytes. Glob. Change Biol. 23, 108–116. doi: 10.1111/gcb.13405
Li, Y., Xia, M., Zhao, X., Yang, J., Li, G., Sun, Z., et al. (2022). Temperature is a cryptic factor to shape the geographical pattern of genetic variation in Ceratophyllum demersum across a subtropical freshwater lake. Authorea October 06, pp. 1–12. doi: 10.22541/au.166505256.63665199/v1
Meglécz, E., Pech, N., Gilles, A., Dubut, V., Hingamp, P., Trilles, A., et al. (2014). QDD version 3.1: A user-friendly computer program for microsatellite selection and primer design revisited: Experimental validation of variables determining genotyping success rate. Mol. Ecol. Resour. 14, 1302–1313. doi: 10.1111/1755-0998.12271
Moore, M. J., Bell, C. D., Soltis, P. S., Soltis, D. E. (2007). Using plastid genome-scale data to resolve enigmatic relationships among basal angiosperms. Proc. Natl. Acad. Sci. U.S.A. 104, 19363–19368. doi: 10.1073/pnas.0708072104
O’Briain, R. (2019). Climate change and European rivers: An eco-hydromorphological perspective. Ecohydrology 12, e2099. doi: 10.1002/eco.2099
Orsini, L., Vanoverbeke, J., Swillen, I., Mergeay, J., De Meester, L. (2013). Drivers of population genetic differentiation in the wild: isolation by dispersal limitation, isolation by adaptation and isolation by colonization. Mol. Ecol. 22, 5983–5999. doi: 10.1111/mec.12561
Ouborg, N. J., Piquot, Y., Van Groenendael, J. M. (1999). Population genetics, molecular markers and the study of dispersal in plants. J. Ecol. 87, 551–568. doi: 10.1046/j.1365-2745.1999.00389.x
Peakall, R., Smouse, P. E. (2006). GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol. Ecol. Notes 6, 288–295. doi: 10.1111/j.1471-8286.2005.01155.x
Peakall, R., Smouse, P. E. (2012). GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics 28, 2537–2539. doi: 10.1093/bioinformatics/bts460
Podani, J. (2001). “SYN-TAX computer programs for data analysis in ecology and systematics,” in User’s Manual (Budapest: Scientia).
Pollux, B. J. A., Luteijn, A., Van Groenendael, J. M., Ouborg, N. J. (2009). Gene flow and genetic structure of the aquatic macrophyte Sparganium emersum in a linear unidirectional river. Freshw. Biol. 54, 64–76. doi: 10.1111/j.1365-2427.2008.02100.x
Pritchard, J. K., Stephens, M., Donnelly, P. (2000). Inference of population structure using multilocus genotype data. Genetics 155, 945–959. doi: 10.1093/genetics/155.2.945
Qadri, H., Uqab, B., Javeed, O., Dar, G. H., Bhat, R. A. (2022). Ceratophyllum demersum-An accretion biotool for heavy metal remediation. Sci. Total. Environ. 806, 150548. doi: 10.1016/j.scitotenv.2021.150548
Reitsema, R. E., Meire, P., Schoelynck, J. (2018). The future of freshwater macrophytes in a changing world: Dissolved organic carbon quantity and quality and its interactions with macrophytes. Front. Plant Sci. 9. doi: 10.3389/fpls.2018.00629
Reitsema, R. E., Wolters, J.-W., Preiner, S., Meire, P., Hein, T., De Boeck, G., et al. (2020). Response of submerged macrophyte growth, morphology, chlorophyll content and nutrient stoichiometry to increased flow velocity and elevated CO2 and dissolved organic carbon concentrations. Front. Environ. Sci. 11. doi: 10.3389/fenvs.2020.527801
Rivaes, R. P., Feio, M. J., Almeida, S. F. P., Calapez, A. R., Sales, M., Gebler, D., et al. (2022). River ecosystem endangerment from climate change-driven regulated flow regimes. Sci. Total Environ. 818, 151857. doi: 10.1016/j.scitotenv.2021.151857
Robertson, M., Schrey, A., Shayter, A., Moss, C. J., Richards, C. (2017). Genetic and epigenetic variation in Spartina alterniflora following the Deepwater Horizon oil spill. Evol. Appl. 10, 792–801. doi: 10.1111/eva.12482
Santamaría, L. (2002). Why are most aquatic plants widely distributed? Dispersal, clonal growth and small-scale heterogeneity in a stressful environment. Acta Oecol. 23, 137–154. doi: 10.1016/S1146-609X(02)01146-3
Schlötterer, C. (2000). Evolutionary dynamics of microsatellite DNA. Chromosoma 109, 365–371. doi: 10.1007/s004120000089
Shukla, D., Kesari, R., Mishra, S., Dwivedi, S., Tripathi, R. D., Nath, P., et al. (2012). Expression of phytochelatin synthase from aquatic macrophyte Ceratophyllum demersum L. enhances cadmium and arsenic accumulation in tobacco. Plant Cell Rep. 31, 1687–1699. doi: 10.1007/s00299-012-1283-3
Szalontai, B., Stranczinger, S., Mesterházy, A., Scribailo, R. W., Les, D. H., Efremov, A. N., et al. (2018). Molecular phylogenetic analysis of Ceratophyllum L. taxa: a new perspective. Bot. J. Linn. Soc 188, 161–172. doi: 10.1093/botlinnean/boy057
Thomaz, S. M. (2023). Ecosystem services provided by freshwater macrophytes. Hydrobiologia 850, 2757–2777. doi: 10.1007/s10750-021-04739-y
Triest, L., Tran Thi, V., Le Thi, D., Sierens, T., Van Geert, A. (2010). Genetic differentiation of submerged plant populations and taxa between habitats. Hydrobiologia 656, 15–27. doi: 10.1007/s10750-010-0434-1
Van Oosterhout, C., Hutchinson, W. F., Wills, D. P. M., Shipley, P. (2004). MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol. Ecol. Notes 4, 535–538. doi: 10.1111/j.1471-8286.2004.00684.x
Vieira, M. L. C., Santini, L., Diniz, A. L., Munhoz, C. D. F. (2016). Microsatellite markers: what they mean and why they are so useful. Genet. Mol. Biol. 39, 312–328. doi: 10.1590/1678-4685-GMB-2016-0027
Wade, P. M. (1993). The influence of vegetation pre-dredging on the post-dredging community. J. Aquat. Plant Manage. 31, 141–144.
Yang, Y., Sun, P., LV, L., Wang, D., Ru, D., Li, Y., et al. (2020). Prickly waterlily and rigid hornwort genomes shed light on early angiosperm evolution. Nat. Plants 6, 215–222. doi: 10.1038/s41477-020-0594-6
Keywords: microsatellites, multilocus genotypes, vegetative dispersion, Danube, river channel, backwaters
Citation: Engloner AI, Németh K, Kós PB, Meglécz E and Bereczki J (2023) Genetic diversity of the submerged macrophyte Ceratophyllum demersum depends on habitat hydrology and habitat fragmentation. Front. Plant Sci. 14:1277916. doi: 10.3389/fpls.2023.1277916
Received: 15 August 2023; Accepted: 23 October 2023;
Published: 09 November 2023.
Edited by:
Guixiang Yuan, Hunan Agricultural University, ChinaReviewed by:
Wei Li, Nanchang Institute of Technology, ChinaBeibei Hao, Guangdong Academy of Sciences, China
Qian Hu, Lushan Botanical Garden (CAS), China
Copyright © 2023 Engloner, Németh, Kós, Meglécz and Bereczki. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Attila I. Engloner, engloner.attila@ecolres.hu