Evolutionary unraveling: new insights into the Persicaria amphibia complex

The Persicaria amphibia complex exhibits significant morphological variation depending on its habitat, existing in either aquatic or terrestrial forms. Traditionally, four distinct elements have been recognized based on morphological features along with their distinct geographical distributions. Recent studies suggest that the Asian element may be genetically distinct from the European and American elements. However, a comprehensive study on the genetic differentiation among all four elements remains lacking. This study aimed to leverage whole plastid genome sequences and ITS2 haplotypes to comprehensively assess the genomic diversity within the P. amphibia complex. Notably, we included multiple individuals from New York State to resolve the ongoing debate regarding the taxonomic status of two American elements – whether they represent a single species or distinct entities. Our analysis revealed a well-supported monophyletic clade encompassing all four elements, endorsing their own section, Amphibia. Notably, the terrestrial form of the American element is sister to all other elements, suggesting it deserves its own species status. This reinstates its historical name, P. coccinea, separating it from the broader P. amphibia. Furthermore, distinct compositions of the ITS2 haplotypes differentiated the four elements, although the European element should be further investigated with more sampling. The most intriguing discovery is the identification of putative hybrids between the two American elements. In one population out of four putative hybrid populations, all three entities – the two parent species and their hybrid offspring – thrive together, showcasing a fascinating microcosm of ongoing evolutionary processes. Unraveling the intricate genetic tapestry within each American species and their hybrid populations remains a compelling next step. By delving deeper into their genetic makeup, we can gain a richer understanding of their evolutionary trajectories and the intricacies of their interactions. Finally, it is estimated that the two species of sect. Amphibia diverged approximately 4.02 million years ago during the Pliocene epoch, when there was a significant global cooling and drying trend.

1 Introduction

Aquatic plants are on the front lines of climate change, such as temperature increases or low precipitation, exhibiting shifting in species composition, a shrinking of their range and distribution, as well as a decline in species richness (Short et al., 2016). Submerged freshwater plants are vital players in the delicate ecosystems they inhabit. They offer refuge, habitat, and a food source for diverse organisms. Their loss would ripple through the entire wetland community, making their proper management and protection a priority for conservation efforts. Identifying distinct species is crucial for effective conservation programs. This knowledge allows us to develop our conservation strategies to meet the specific needs of each vulnerable plant, ensuring their continued survival and the health of the ecosystems they sustain.

Persicaria amphibia (L.) Delarbre is an illustrative case of a species whose identity lacks clear consensus within the scientific community. Nonetheless, it holds significant ecological importance within wetland ecosystems. Serving as a primary producer, it provides habitat for various aquatic and semi-aquatic organisms. It contributes to the improvement of water quality by absorbing excess nutrients and filtering pollutants, like many other wetland plants.

The P. amphibia complex is the most widespread amphibious species within Polygonaceae. It exhibits dynamic morphological transformations depending on where it grows. In aquatic forms, the plants bloom in water with short-cylindric inflorescences (< 4 cm long) (Figure 1A). They are glabrous and have long, floating stems with elliptic, obtuse, or rounded leaves and adventitious roots. When they grow on land, they often develop flared ocreae, but this character has only been reported from North American plants (Figure 1B) (Mitchell, 1968, 1971). In terrestrial forms, the plants bloom on soil with elongate-cylindric inflorescences (> 4 cm long) (Figure 1C). They are often pubescent and have erect or decumbent stems with relatively large linear or ovate, acuminate, strigose leaves and erect roots. In contrast to aquatic forms, they have ocreae that are entire (Mitchell, 1968, 1971; Partridge, 2001; Reveal and Atha, 2010). This terrestrial form can persist under water, but it does not develop floating leaves (Figure 1D).

Figure 1

Figure 1 Two North American elements of the P. amphibia complex. (A) P. amphibia var. stipulacea (aquatic form), (B) P. amphibia var. stipulacea (terrestrial form), (C) P. coccinea (= P. amphibia var. emersa) (terrestrial form), (D) P. coccinea (aquatic form). Photo credits: (A) Nick Belliveau, (B) Mi-Jeong Yoo, (C) Charley Amos, (D) Justin Chappelle.

Being native to the northern hemisphere and spanning across Europe, Asia, and North America, this species complex has led to the recognition of four distinct elements based on their morphological features and geographic distribution patterns: 1) The Asian element (var. amurensis (Korsh.) H. Hara) characterized by white flowers and cordate leaves; 2) The European element (var. amphibia) with pink flowers, which resembles the Asian element; and 3) two different American elements (var. emersa (Michx.) J.C. Hickman and var. stipulacea (N. Coleman) H. Hara).

In North America, two elements were recognized based on their morphological characteristics; the terrestrial form, P. amphibia var. emersa, which occurs on soil but does not develop floating leaves when submerged (Figures 1C, D), and the aquatic form, P. amphibia var. stipulacea, which has flared ocreae when stranded on soil (Figures 1A, B). However, highly variable morphological transformations were often observed depending on their habitats among these varieties, so their taxonomic entities have been treated differently. For example, Hinds and Freeman (Hinds and Freeman, 2005) integrated these two types into one species without variety ranks. However, Reveal and Atha (Reveal and Atha, 2012) differentiated P. amphibia var. emersa into a distinct taxon called P. coccinea (Muhl. ex Willd.) Greene and maintained P. amphibia var. stipulacea as it is. In this study, we will use the Reveal and Atha (2012) nomenclature and refer to P. amphibia var. emersa as P. coccinea throughout the paper.

To overcome the limitation of morphological traits in differentiation of this complex, molecular approaches have been taken, which have revealed a clear distinction between P. amphibia and other Persicaria species (Kim et al., 2008; Kim and Donoghue, 2008a, 2008b). Another study of morphological data of 107 individuals from the P. amphibia complex from Asia and North America showed that aquatic forms can be discerned from terrestrial forms by leaf size, shape, and petiole length (Gao et al., 2013). Despite their morphological differences, sequences from four regions of the chloroplast genome failed to distinguish between the two forms, indicating a low level of genetic differentiation between them. In Gao et al. (2013)’s study, 13 individuals from two populations in the United States (California and New York) were included, but only one individual represented the aquatic form. Among the 13 individuals examined, two haplotypes were identified, with a difference of only 16 of 2,760 bp (0.6%). This suggests that the variation in the chloroplast genome within this complex may be very minimal. However, this study focused on a small portion of the chloroplast genome. Therefore, examining both the whole plastome sequences and nuclear regions could provide a more comprehensive understanding of the genetic divergence within this complex. Furthermore, such analysis would help elucidate whether the proposed four elements demonstrate genetic differentiation.

As a polyploid species (2n = 66, 88, or 96 for P. amphibia, 2n = 66 for P. coccinea) (Jaretzky, 1928; Sultan et al., 1998), the P. amphibia complex typically harbors multiple copies of genes within its nuclear genome, often complicating the analysis of nuclear regions in phylogenetic studies. This applied to the internal transcribed spacer (ITS1 and ITS2) of a segment of nuclear ribosomal DNA (nrDNA). ITS regions have been extensively sequenced and widely employed as a molecular marker in plant systematics due to their valuable features. These include conserved regions suitable for designing universal primers, ease of amplification, and sufficient variability to differentiate closely related species (Coleman, 2003; Schultz and Wolf, 2009; Gao et al., 2010; Pang et al., 2011). Although ITS regions are well known for their homogeneous sequences resulting from concerted evolution, intra-individual polymorphisms in nrDNA have been documented across various plant lineages, including non-hybrid diploids (Xiao et al., 2010; Huang et al., 2016; Xu et al., 2017; Asanuma et al., 2019; Li et al., 2019; Zhang et al., 2022). Thus, the exploration of nuclear markers like ITS regions or single-copy gene areas (e.g., LEAFY second intron) requires a laborious cloning process (Kim et al., 2008), often limited to detecting copies occurring in high abundance. Recently, Next Generation Sequencing (NGS) technology has facilitated the direct sequencing of ITS regions at a sufficient depth. This advancement has unveiled significantly higher percentages of intragenomic polymorphisms within eukaryote genomes, including those of plants (Wang et al., 2023).

In this study, we investigated genomic diversity and population differentiation of the P. amphibia complex using the whole plastome sequences and the ITS2 region. This study aimed to elucidate the species status of the P. amphibia complex in New York State (NY) and tried to answer the following questions: 1) How much genomic variation exists among 12 populations in NY? 2) Are two varieties of this species genetically distinct? 3) If yes, when had these two diverged? 4) Can ITS2 effectively identify these two varieties and their hybrids, if present? 5) Can whole plastome and ITS2 sequences distinguish each element from three regions (America, Asia, and Europe)? The findings of this study will inform updates to the species information in the Flora of North America and provide enhanced insights into the taxonomic status and evolutionary history of this complex.

2 Materials and methods

2.1 Taxon sampling

To investigate the taxonomic status of the P. amphibia complex, leaf tissues were collected from at least ten individuals in each of the 12 populations found in New York (Table 1). These twelve populations encompass the two American elements and their putative hybrid populations based on their morphological features (Table 2). One individual from each population was selected for DNA library construction, except the Sandy Creek population, where four individuals were chosen. Additionally, one individual from the United Kingdom and two individuals from Asia were obtained from the Royal Botanic Gardens, Kew, and the National Institute of Biological Resources (NIBR; South Korea), respectively (Table 1).

Table 1

Table 1 The information of P. amphibia complex samples and the population acronyms used in this study.

Table 2

Table 2 Morphological features of samples used in plastome sequencing.

2.2 DNA library preparation, sequencing, and assembly

Total DNAs were extracted from silica-gel dried leaves with a DNeasy Plant Mini kit (Qiagen, Hilden, Germany) and quantified using a NanoDrop™ One (Thermo, Waltham, MA, USA). Seventeen DNA libraries (two Asian individuals, one European individual, and 14 American individuals) were constructed with 100 ng of DNA using the NEBNext^® Ultra™ II DNA Library Prep Kit for Illumina^® (New England Biolabs, Ipswich, MA, USA). Each library was indexed with multiplex oligos and quantified with the NEBNext Library Quant Kit (New England Biolabs, Ipswich, MA, USA). After quantification, equimolar libraries (5 nM each) were pooled together and sequenced with 2 ×150 bp on an Illumina NovaSeq 6000 (Illumina, San Diego, CA, USA) at the Interdisciplinary Center for Biotechnology Research (ICBR) of the University of Florida. The obtained sequences were cleaned with Trimmomatic v. 0.36 (Bolger et al., 2014) and assembled into the plastome sequences and nuclear ribosomal DNA (nrDNA) sequences using GetOrganelle 1.7.5.0 (Jin et al., 2020). The assembled sequences were then confirmed with Geneious v.2024.0.1 (http://www.geneious.com).

2.3 Sequence alignment and phylogenetic analysis

Assemble plastome sequences were annotated using GeSeq (Tillich et al., 2017) and visualized with OGDRAW (Greiner et al., 2019). These sequences were aligned with MAFFT v.7.520 (Katoh and Standley, 2013) and manually adjusted. Phylogeny was inferred with the maximum likelihood (ML) method using IQ-TREE v.2.2.2.7 (Nguyen et al., 2015). For the ML tree, molecular evolution models were evaluated using ModelFinder (Kalyaanamoorthy et al., 2017), and the standard nonparametric bootstrap method (-b 1000) was applied to evaluate the support of retrieved clades. Plastome sequences from members of sections Persicaria, Tovara, Echinocaulon, and Cephalophilon were retrieved from NCBI GenBank. The sequences of sect. Cephalophilon were utilized as outgroups in the analysis.

2.4 Genotyping population using two chloroplast regions

To evaluate the origin of putative hybrid populations and genetic differentiation within and between populations, two variable regions of the plastome between P. amphibia var. stipulacea and P. coccinea were investigated. Out of seven intergenic spacers (IGS) and intron regions analyzed, we chose two regions that showed the most substantial size difference (Supplementary Table 1). Two target areas were amplified in a 50 µL PCR reaction that included 5 – 50 ng of DNA, 0.5 µM of each primer, and 25 µL of OneTaq^® 2× Master Mix (New England Biolabs, Ipswich, MA, USA). The PCR was performed as follows: 3 minutes of initial denaturation at 95°C, followed by 35 cycles of 30 seconds at 95°C, 30 seconds at 55°C, 45 seconds at 72°C, and the final extension of 5 minutes at 72°C. The amplicons were assessed using a 1.2% agarose gel, running at 100 V for 30 min.

2.5 ITS2 haplotype identification

To identify ITS2 haplotypes, we initially aligned short reads to the ITS regions and assessed sequence variations. Given that ITS2 showed more variability than ITS1, we extracted representative haplotypes from each sequencing dataset. Subsequently, 133 bp of ITS2 haplotypes were tallied from the short reads of each DNA library using a custom script. During this process, short reads that fully encompassed the entire 133 bp of ITS2 were counted, and the results were presented as a pie chart.

2.6 Divergence time estimates

The divergence times of the P. amphibia complex were estimated using BEAST v.2.6.2 (Bouckaert et al., 2019) based on a plastid ML tree. First, BEAST input XML file was generated using BEAUti v2.7.5. A GTRGAMMA evolutionary model was used with 10 gamma rate categories. A log-normal relaxed clock (Drummond et al., 2006) was the clock model, using the ML tree estimated above as the starting tree, and with speciation modeled by the Yule process. Age calibrations were implemented as follows: a log-normal distribution with a mean of 1.0, a standard deviation of 1.0, and an offset of 31.14 for the crown group of sections Persicaria, Tovara, and Echinocaulon, and an offset of 14.29 and 17.39 for sections Cephalophilon and Echinocaulon, respectively (Cao et al., 2022). We ran five separate runs for a total of 50 million generations, sampling every 1000. Runs were combined using LogCombiner v.2.7.5 with a 10% burnin, and a maximum clade credibility (MCC) tree with mean node heights was created using a sample of 18,002 trees in TreeAnnotator v.2.7.5. The chronogram was drawn using FigTree v.1.4.3 (http://tree.bio.ed.ac.uk/software/figtree/).

3 Results

3.1 Assembled nrDNA and plastome sequences

Assembled nrDNA is about 5,868 bp, but there was little sequence variation (0 – 10 bp differences among the P. amphibia complex; data not shown). In the ITS region of nrDNA, the Asian and European elements have identical sequences with some gaps, while the American elements showed a difference of seven out of 750 bp compared to the Asian and European elements. Thus, due to their low sequence divergence and heterogeneity, nrDNA sequences were not employed in phylogenetic analysis. Instead, a thorough investigation was conducted on the compositions of ITS2 haplotypes (refer to the information below).

Assembled plastome sequences range from 159,320 bp (P. amphibia var. stipulacea from Sandy Creek) to 159,539 bp (P. coccinea from Vischer Ferry) and have a typical structure of a plastid genome with two inverted repeats (IRA, IRB), which separate a long single copy section (LSC) from a short single copy section (SSC). However, the assembled plastomes showed no structural variations among the individuals investigated here (Supplementary Figure 1). There are 79 protein-coding genes, 4 rRNA genes, and 30 tRNA genes arranged in the same gene order (Supplementary Figure 1).

3.2 Phylogenetic analysis of plastome sequences

The analysis of plastome sequences revealed a strongly supported monophyletic clade that includes all four elements, endorsing their own section, Amphibia, as well as confirming the monophyly of each section of Persicaria (Figure 2). However, the American elements failed to form a single clade: the American element I (P. amphibia var. stipulacea) was closely placed in the clade containing the Asian and European elements, while the American element II (P. coccinea) was sister to the clade that includes the other three elements (Figure 2). The genetic distance between American element II and the other three elements was much greater compared to those observed among the members of sect. Persicaria and Tovara. Thus, American element II should be considered a distinct species, P. coccinea.

Figure 2

Figure 2 A maximum-likelihood tree inferred from the entire chloroplast sequences with the best fit model of K3Pu+F+R2. Four elements were indicated by different colors. Red: Asian element, purple: European element, green: American element I, blue: American element II. The aligned length was 168,017 bp, of which 90.8% (152,514 bp) were constant across taxa and 6.8% sites (11,365) were parsimony informative. The navy circle on the branch shows 100% bootstrap support value. Asterisk (*) indicates putative hybrids. For sequences downloaded from NCBI, their accession numbers were presented next to the species name.

Interestingly, all putative hybrid species were included in the clade of P. amphibia var. stipulacea, except one individual (SC2) from Sandy Creek (Figure 2). This individual is morphologically similar to that of P. coccinea, while the other three individuals (SC1, SC3, and SC4), containing the plastome of P. amphibia var. stipulacea, exhibit similarities to P. coccinea and P. amphibia var. stipulacea, respectively (Tables 2, 3; Figure 3). Thus, there might be all three taxa, two parental species and their putative hybrid species, occurring together in Sandy Creek (see below).

Table 3

Table 3 Comparison of morphological and molecular characteristics of putative hybrids.

Figure 3

Figure 3 Putative hybrids from four populations. (A) P. amphibia var. stipulacea-like aquatic form and (B) P. coccinea-like terrestrial form from Guppy Falls (GF); (C) P. amphibia var. stipulacea-like aquatic form and (D) P. coccinea-like terrestrial form from Lake Alice (LA); (E) P. amphibia var. stipulacea-like terrestrial form and (F) P. coccinea-like terrestrial form from Sandy Creek (SC); (G) P. amphibia var. stipulacea-like aquatic form and (H) P. coccinea-like terrestrial form from Saranac Lake (SC). Photo credits: (A, B) Randy A. Nonenmacher, (C–H) Mi-Jeong Yoo.

In the meantime, the European element was placed as a sister to the Asian element in the clade of P. amphibia (Figure 2), but its branch length was too short, and its sequence divergence was too low (0.06%) to differentiate it from the Asian element. Therefore, we analyzed the partial chloroplast sequences of a Greek individual from the previous work (Gitsopoulos et al., 2013) alongside the current data to determine the differentiation between the European and Asian elements. While the monophyletic clades of two species, P. amphibia and P. coccinea, were retrieved in agreement with the ML tree inferred from the entire plastome data, two individuals from Europe (UK, Greece) failed to form their own clade (Supplementary Figure 2).

3.3 Genotyping populations using two chloroplast regions

Based on the chloroplast sequences, we designed primers targeting the rpl20-rpl12 IGS and rpl33-rps18 IGS, which yielded different-sized fragments between two American elements. Utilizing the regions, we conducted screenings of individuals from 12 populations. Based on morphological features, four hybrid populations were identified. The plants in four hybrid populations exhibited a morphological resemblance to P. coccinea. However, noteworthy variations were observed; for example, certain individuals developed floating leaves when grown in water and exhibited flared ocreae when stranded on dry land (Table 2; Figure 3).

The screening results unequivocally confirmed the identities of P. amphibia var. stipulacea, P. coccinea, and their hybrids (Supplementary Figure 3). Of particular interest, the Sandy Creek population displayed the co-occurrence of two plastomes, while other hybrid populations contained only the plastome of P. amphibia var. stipulacea (Supplementary Figure 3). Additionally, it is worth noting that most hybrids possess a plastome that is similar to the American element I (P. amphibia var. stipulacea). This suggests that P. amphibia var. stipulacea may have been the maternal parent during the hybridization events.

3.4 ITS2 haplotype analysis

To find out how many different ITS2 haplotypes exist, short reads were assembled into nrDNA and examined manually. The study revealed the presence of three ITS2 haplotypes in the Asian element, five in the European element, and seven in the American elements (Figure 4A). Then, the frequencies of each haplotype were determined based on around 1000 short reads, that provided complete coverage of the 133 bp region of the ITS2 haplotype. As a result, the Type2a haplotype accounted for 88% and 86% of the ITS2 haplotypes in individuals from Korea and Russia, respectively. In the European element, the Type1, Type2a, and Type2d haplotypes together represented 92% of the ITS2 haplotypes (Figure 4B). The composition of ITS2 haplotypes in the American elements exhibited variations. Specifically, four individuals of P. amphibia var. stipulacea (GL, HY, ML, and WT) and one individual from Sandy Creek (SC4) were found to have just two haplotypes (Type1 and US1). On the other hand, P. coccinea and putative hybrids had a higher number of haplotypes, ranging from five to seven (Figure 4C). In addition, the geographical differences were evident in the composition of ITS2 haplotypes in P. coccinea. Among the P. coccinea individuals studied, two individuals from the southern region (RL, VF) and one individual from Sandy Creek (SC2) were found to lack Type1 and US1, which are prominent haplotypes observed in the remaining individuals. Furthermore, three distinct compositions of ITS2 haplotypes were identified in the Sandy Creek populations: 1) SC1 and SC3 exhibit similarities to northern P. coccinea (PD and BL), 2) SC2 shows similarity to southern P. coccinea (RL and VF), and 3) SC4 displays characteristics comparable to P. amphibia var. stipulacea (Figure 4C).

Figure 4

Figure 4 ITS2 haplotypes sequenced in the P. amphibia complex. (A) Sequences of ITS2 haplotypes, (B) pie charts showing the compositions of ITS2 haplotypes in P. amphibia var. amurensis (Korea, Russia) and P. amphibia var. amphibia (UK), (C) Pie charts showing the compositions of ITS2 haplotypes in P. coccinea (blue circles), P. amphibia var. stipulacea (red circles), and their putative hybrids (yellow circles) collected from New York, USA.

3.5 Divergence time estimates

To estimate divergence times for the P. amphibia complex, we utilized the ML tree obtained above and three calibration points that were reported previously (Cao et al., 2022). The estimated divergence times for lineages exhibited variation in comparison to the previous work based on 80 plastome sequences of Polygonaceae (Cao et al., 2022). For example, the genus Persicaria diverged 36.10 million years ago (Ma) (95% HPD: 33.18–39.01 Ma) in this study, which was younger than 42.78 Ma (95% HPD: 34.72–45.36 Ma) in the previous work of Cao et al. (2022) (Figure 5). However, the 95% HPD for the major clades showed overlap with the 95% HPD of Cao et al. (2022)’s estimates.

Figure 5

Figure 5 Divergence time estimates for Persicaria species inferred from the whole plastome sequences. Blue bars represent 95% highest posterior density (HPD) intervals. Orange stars show the calibration points used in this study.

Based on our analysis, we determined that sections Persicaria and Amphibia diverged approximately 13.46 Ma (95% HPD: 11.01–16.21 Ma), while two members of sect. Amphibia split apart around 4.02 Ma (95% HPD: 2.95–5.42 Ma) (Figure 5). Finally, P. amphibia var. stipulacea diverged approximately 1.21 Ma from other elements.

4 Discussion

4.1 Support of section Amphibia and two distinct species in American elements

The genus Persicaria has about 100 species and can be further classified into four sections: Persicaria, Tovara, Echinocaulon, and Cephalophilon (Haraldson, 1978). Among them, sect. Persicaria stands out with the highest number of taxa. Notably, the P. amphibia complex has long resided within this section. However, recent molecular evidence (Kwak et al., 2006; Kim and Donoghue, 2008b; Cao et al., 2022) and their unique morphological features, such as pollen typology (Kong et al., 2021) and the presence of the rhizome, placed this complex as its own section of Amphibia. Nevertheless, the previous investigations failed to examine all four elements of the P. amphibia complex. In this study, we incorporated all four elements for the first time to elucidate the species entities and their relationships.

The P. amphibia complex has been well known for its high polymorphisms, depending on its environment. This variability had led to many synonyms; for example, 187 names were listed in Plants of the World Online (POWO). However, it has long been recognized that the American elements were clearly differentiated from Asian or European elements as they possessed flared ocreas when stranded. A recent analysis of nearly 10,000 observations on iNaturalist shed further light on this taxonomic tangle by confirming a clear morphological distinction between the two American elements, although there were a handful of exceptions (less than 0.1% of all observations) (Atha, 2024). Therefore, Atha (2024) proposed the existence of two species, P. amphibia var. stipulacea and P. coccinea, within the American elements. However, no molecular data has been shown to support this claim.

Here, we delved deep into the plastome sequences of the American elements, seeking to establish their species entity. Furthermore, we included Asian and European elements to evaluate whether they could be genetically differentiated from the American elements. We assembled 17 plastome sequences, and upon analysis, it was seen that four elements were clustered together, receiving 100% bootstrap support (Figure 2).

In the Amphibia clade, there were remarkable genetic similarities among four individuals of the Asian element, while the European individual exhibited slight genetic differentiation from the Asian individuals (Figure 2). Even though the European element was sister to the clade of the Asian element (Figure 2), it is important to examine this element more thoroughly because of its wide distribution in Europe and the failure to form a monophyletic group in the analysis of two European individuals (Supplementary Figure 2). In addition, there were some shared features in the ITS2 haplotypes between Asian and European elements, such as Type2a and Type2d (Figure 4B). These results indicate that further investigation with more samples is required to confirm the differentiation of the European element from the Asian element and its taxonomic status. In contrast, the American elements formed two separate clades. The first one included all the aquatic forms, along with some terrestrial forms, and is closely related to a clade consisting of the Asian and European elements (Figure 2). All the aquatic forms (HY, GL, ML, and WT) have typical morphological features of P. amphibia var. stipulacea, such as ovoid and short inflorescence, oblong leaves, short plant height, and flared ocreae when stranded. In addition, they had only two types of ITS2 haplotypes, Type1 of which was shared with the European element (Figures 4B, C), supporting their close relationship as shown in the plastome phylogeny (Figure 2). Therefore, based on their consistent morphology and plastome sequences, the members of this clade are classified as P. amphibia var. stipulacea.

The other clade contains terrestrial forms of P. coccinea-like plants only and is sister to the other elements (Figure 2). This group, genetically separated from the other elements by 0.32 – 0.34% in the whole plastome sequences (twice as high as that between P. orientalis and P. lapathifolia), deserves its own species recognition. Also, the members of this clade (BL, PD, RL, and VF) had five to seven different types of ITS2 haplotypes (Figure 4C). Thus, the American element II should be recognized as a distinct species, P. coccinea, based on its unique morphology and significant genetic differentiation from the other elements.

As a result, the phylogenetic analysis of plastome sequences strongly corroborates sect. Amphibia, which comprises two species, P. amphibia and P. coccinea (Figure 2). Notably, the morphological and molecular characteristics of the full plastome sequences and ITS2 haplotype composition provide evidence for the existence of two distinct species in American elements: P. amphibia var. stipulacea and P. coccinea. This suggests independent evolutionary trajectories for these distinct species (Figure 4C). Specifically, we estimated that their divergence occurred 4.02 Ma (Figure 5), during the Pliocene epoch when there was a significant global cooling and drying trend. This suggests that their distinct environmental preferences may have influenced their evolutionary process over a period. For instance, despite both species being amphibious, these species exhibit differences in population size depending on their habitat. Persicaria amphibia forms significantly larger populations in aquatic environments compared to terrestrial habitats, suggesting a better adaptation to such habitats. Conversely, P. coccinea prefers terrestrial habitats but can tolerate temporary flooding. These adaptations are reflected in their distribution patterns: while both species are found throughout North America, P. amphibia var. stipulacea is absent south of the Laurentide Ice Sheet, where aquatic habitats would differ from those to the north of it (Atha, 2024).

4.2 Hybridization between Two American Elements

As Atha (2024) noted, there are several reports on American elements that have intermediate phenotypes between P. amphibia var. stipulacea and P. coccinea. We also identified four populations that display intermediate morphological characteristics (GF, LA, SC, and SL; Table 2).

Interestingly, in a plastome phylogenetic tree, the clade of P. amphibia var. stipulacea includes some individuals which have long, cylindrical inflorescence, ovate leaves, and tall plant height that are the hallmarks of the other American element, P. coccinea (Figure 2). This unexpected merging raises intriguing inquiries regarding the delineations between these two American entities. When two chloroplast regions were used to look at multiple individuals from 12 populations, all putative hybrid populations except Sandy Creek had the fragment types of P. amphibia var. stipulacea (Supplementary Figure 3). These populations seem to be of hybrid origin based on their phenotypes and contradicting plastome sequences (Tables 2, 3; Figure 3). Although multiple individuals with P. coccinea-like phenotypes from those populations were sequenced for the complete plastome sequences, they had P. amphibia var. stipulacea-like genomes (Table 3; Figures 2, 3). In addition, individuals that phenotypically resemble P. amphibia var. stipulacea also exist in these populations with flared ocreae when stranded (Figure 3E). However, unlike typical P. amphibia var. stipulacea, the aquatic forms did not grow far from the water’s edge, preferring to root at the water-soil border (Figures 3A, C, G). This behavior further strengthens the case for their hybrid origin. Thus, all four hybrid populations might have experienced hybridization between the two elements, with P. amphibia var. stipulacea likely playing as the maternal parent.

The clade of P. coccinea contains one putative hybrid individual (SC2), which morphologically resembles P. coccinea, from the Sandy Creek population. As mentioned above, in contrast to the other hybrid populations with one type of plastome, the Sandy Creek population held two different types of plastomes (Supplementary Figure 3). The complete plastome sequences from four individuals revealed an intriguing story: SC1, which looked like P. coccinea, had a P. amphibia var. stipulacea-like plastid genome. SC3 and SC4, on the other hand, showed a typical morphology of P. amphibia var. stipulacea with a P. amphibia var. stipulacea-like plastid genome. However, SC2 resembled P. coccinea morphologically and had a P. coccinea-like plastome (Figure 2; Table 3).

When the compositions of ITS2 haplotypes were examined, it was found that two individuals (SC1, SC3) had ITS2 haplotypes that closely resembled those of P. coccinea obtained from northern populations, PD and BL (Table 3). Interestingly, SC2 had ITS2 haplotypes that were similar to those of P. coccinea from the southern regions, RL and VF. Nevertheless, SC4 contained two distinct ITS2 haplotypes of P. amphibia var. stipulacea. These findings indicate that all three forms – the two parent species and their hybrid offspring – thrive together in the Sandy Creek population. However, a P. coccinea-like individual (SC2) may not be considered a parent species due to its different ITS2 haplotypes compared to putative hybrids (SC1, SC3). Thus, comprehensive study is warranted to investigate the species composition and population structure within this “hybrid swarm”, which serves as a tangible example of ongoing evolutionary processes. Furthermore, given its specific geographical location adjacent to Lake Ontario, it is critical to thoroughly investigate the extent of these hybrid swarms near Lake Ontario and along the St. Lawrence River.

5 Conclusions

Recognition of accurate species identity is crucial for establishing conservation programs aimed at preserving unique genetic resources. While morphological variation often presents challenges in species determination due to its continuous nature, molecular tools offer a more reliable approach. In this study, we employed whole plastome sequences and ITS2 haplotype composition to differentiate species within the P. amphibia complex, focusing on two American elements. Our findings support the recognition of P. amphibia and P. coccinea as distinct species within sect. Amphibia, while also highlighting significant genetic divergence among four geographic elements. Despite the challenges posed by the polyploid nature of the species in traditional sequencing methods, NGS enabled us to successfully determine both ITS2 haplotype composition and whole plastome sequences of members of sect. Amphibia. However, it is worth noting that this study utilized only one representative individual per population and focused on one nuclear region. Therefore, further investigation is necessary to explore population structure, utilizing a broader sample size and examining additional genomic features.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, PRJNA1093078 and https://figshare.com/articles/dataset/Phylogenetic_analysis_using_plastid_genome_sequences/23668770.

Author contributions

GB: Formal analysis, Investigation, Writing – review & editing. AN: Formal analysis, Investigation, Writing – review & editing. YH: Formal analysis, Investigation, Writing – review & editing. MK: Investigation, Writing – review & editing. M-JY: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was supported by the 2023 Huyck Research Grant.

Acknowledgments

We express our gratitude to the Royal Botanic Garden, Kew, for providing the DNA sample of the European element. We are also thankful to Clarkson University and the Office of Information Technology for providing computational resources and support. Our appreciation extends to the iNaturalist users, whose observations were crucial in identifying the location of the plant materials. Finally, we would like to thank the Huyck Preserve for financial support for this endeavor.

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.

The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2024.1408810/full#supplementary-material

References

Asanuma, M., Zhu, S., Okura, N., Cai, S. Q., Yoshimatsu, K., Komatsu, K. (2019). Genetic polymorphism of Japanese cultivated Rheum species in the internal transcribed spacer region of nuclear ribosomal DNA. J. Nat. Med. 73, 541–554. doi: 10.1007/s11418-019-01298-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Atha, D. E. (2024) The Case for Recognizing Persicaria amphibia and Persicaria coccinea as Distinct Species (iNaturalist). Available at: https://www.inaturalist.org/journal/danielatha/43123-the-case-for-recognizing-persicaria-amphibia-and-persicaria-coccinea-as-distinct-species (Accessed December 11 2023).

Google Scholar

Bolger, A. M., Lohse, M., Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170

PubMed Abstract | CrossRef Full Text | Google Scholar

Bouckaert, R., Vaughan, T. G., Barido-Sottani, J., Duchene, S., Fourment, M., Gavryushkina, A., et al. (2019). BEAST 2.5: An advanced software platform for Bayesian evolutionary analysis. PloS Comput. Biol. 15, e1006650. doi: 10.1371/journal.pcbi.1006650

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, D. L., Zhang, X. J., Qu, X. J., Fan, S. J. (2022). Plastid phylogenomics sheds light on divergence time and ecological adaptations of the tribe Persicarieae (Polygonaceae). Front. Plant Sci. 13. doi: 10.3389/fpls.2022.1046253

CrossRef Full Text | Google Scholar

Coleman, A. W. (2003). ITS2 is a double-edged tool for eukaryote evolutionary comparisons. Trends Genet. 19, 370–375. doi: 10.1016/S0168-9525(03)00118-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Drummond, A. J., Ho, S. Y., Phillips, M. J., Rambaut, A. (2006). Relaxed phylogenetics and dating with confidence. PloS Biol. 4, e88. doi: 10.1371/journal.pbio.0040088

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, T., Yao, H., Song, J., Liu, C., Zhu, Y., Ma, X., et al. (2010). Identification of medicinal plants in the family Fabaceae using a potential DNA barcode ITS2. J. Ethnopharmacol. 130, 116–121. doi: 10.1016/j.jep.2010.04.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, Y., Bhandari, G. S., Park, J. H., Park, C.-W. (2013). A systematic study of the Polygonum amphibium L. complex (Polygonaceae) based on chloroplast DNA sequences. Korean J. Pl. Taxon. 43, 34–45. doi: 10.11110/kjpt.2013.43.1.34

CrossRef Full Text | Google Scholar

Gitsopoulos, T. K., Vasilakoglou, I., Tsoktouridis, G. (2013). Persicaria amphibia, a serious terrestrial weed in northern Greece: A combined molecular and morphological approach to identification and taxonomy. Biotechnol. Biotechnol. Equip. 27, 4236–4242. doi: 10.5504/bbeq.2013.0088

CrossRef Full Text | Google Scholar

Greiner, S., Lehwark, P., Bock, R. (2019). OrganellarGenomeDRAW (OGDRAW) version 1.3.1: expanded toolkit for the graphical visualization of organellar genomes. Nucleic Acids Res. 47, W59–W64. doi: 10.1093/nar/gkz238

PubMed Abstract | CrossRef Full Text | Google Scholar

Haraldson, K. (1978). Anatomy and taxonomy in Polygonaceae subfam. Polygonoideae Meissn. emend. Jaretzky (Univ ; Almqvist & Wiksell international (distr.), Uppsala : Stockholm ; New York), 1–95.

Google Scholar

Hinds, H. R., Freeman, C. C. (2005). “Persicaria,” in Flora of North America North of Mexico, ed. F.O.N.a.E. Committee (Oxford University, New York), 574–594.

Google Scholar

Huang, J. F., Li, L., van der Werff, H., Li, H. W., Rohwer, J. G., Crayn, D. M., et al. (2016). Origins and evolution of cinnamon and camphor: A phylogenetic and historical biogeographical analysis of the Cinnamomum group (Lauraceae). Mol. Phylogenet. Evol. 96, 33–44. doi: 10.1016/j.ympev.2015.12.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Jaretzky, R. (1928). Histologische und karyologische studien an Polygonaceen. Jahrb. f. wiss. Botanik 69, 357–490.

Google Scholar

Jin, J. J., Yu, W. B., Yang, J. B., Song, Y., Depamphilis, C. W., Yi, T. S., et al. (2020). GetOrganelle: a fast and versatile toolkit for accurate de novo assembly of organelle genomes. Genome Biol. 21, 241. doi: 10.1186/s13059-020-02154-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Kalyaanamoorthy, S., Minh, B. Q., Wong, T. K. F., Von Haeseler, A., Jermiin, L. S. (2017). ModelFinder: fast model selection for accurate phylogenetic estimates. Nat. Methods 14, 587–589. doi: 10.1038/nmeth.4285

PubMed Abstract | CrossRef Full Text | Google Scholar

Katoh, K., Standley, D. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772–780. doi: 10.1093/molbev/mst010

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, S. T., Donoghue, M. J. (2008a). Incongruence between cpDNA and nrITS trees indicates extensive hybridization within Eupersicaria (Polygonaceae). Am. J. Bot. 95, 1122–1135. doi: 10.3732/ajb.0700008

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, S. T., Donoghue, M. J. (2008b). Molecular phylogeny of persicaria (Persicarieae, polygonaceae). Syst. Bot. 33, 77–86. doi: 10.1600/036364408783887302

CrossRef Full Text | Google Scholar

Kim, S. T., Sultan, S. E., Donoghue, M. J. (2008). Allopolyploid speciation in Persicaria (Polygonaceae): insights from a low-copy nuclear region. Proc. Natl. Acad. Sci. U.S.A. 105, 12370–12375. doi: 10.1073/pnas.0805141105

PubMed Abstract | CrossRef Full Text | Google Scholar

Kong, M. J., Song, J. H., Hong, S. P. (2021). Pollen morphology of Persicaria section Cephalophilon (Polygonaceae). Pl. Syst. Evol. 307, 16. doi: 10.1007/s00606-020-01738-9

CrossRef Full Text | Google Scholar

Kwak, M., Kim, M.-H., Won, H., Park, C.-W. (2006). Variation of nuclear ribosomal ITS sequences of Polygonum section Persicaria (Polygonaceae) in Korea. Korean J. Pl. Taxon. 36, 21–39. doi: 10.11110/kjpt.2006.36.1.021

CrossRef Full Text | Google Scholar

Li, M., Zhao, H., Zhao, F., Jiang, L., Peng, H., Zhang, W., et al. (2019). Alternative analyses of compensatory base changes in an ITS2 phylogeny of Corydalis (Papaveraceae). Ann. Bot. 124, 233–243. doi: 10.1093/aob/mcz062

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitchell, R. S. (1968). Variation in the Polygonum amphibium complex and its taxonomic significance Vol. 45 (Oakland, CA: University of California Press), 1–65.

Google Scholar

Mitchell, R. S. (1971). A guide to aquatic smartweeds (Polygonum) of the United States Vol. 41 (Blacksburg, VA: Virginia Polytechnic Institute and State University), 1–52.

Google Scholar

Nguyen, L. T., Schmidt, H. A., Von Haeseler, A., Minh, B. Q. (2015). IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 32, 268–274. doi: 10.1093/molbev/msu300

PubMed Abstract | CrossRef Full Text | Google Scholar

Pang, X., Song, J., Zhu, Y., Xu, H., Huang, L., Chen, S. (2011). Applying plant DNA barcodes for Rosaceae species identification. Cladistics 27, 165–170. doi: 10.1111/cla.2011.27.issue-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Partridge, J. W. (2001). Persicaria amphibia (L.) Gray (Polygonum amphibium L.). J. Ecol. 89, 487–501. doi: 10.1046/j.1365-2745.2001.00571.x

CrossRef Full Text | Google Scholar

Reveal, J. L., Atha, D. E. (2010). New combinations and typifications in Bistorta, Persicaria, Polygonum, and Rumex (Polygonaceae). Brittonia 62, 243–263. doi: 10.1007/s12228-010-9141-5

CrossRef Full Text | Google Scholar

Reveal, J. L., Atha, D. E. (2012). “Persicaria (L.) Mill. Smartweed,” in Intermountain Flora. Ed. Al, C. E. (Bronx, NY: The New York Botanical Garden Press), 236–250.

Google Scholar

Schultz, J., Wolf, M. (2009). ITS2 sequence-structure analysis in phylogenetics: a how-to manual for molecular systematics. Mol. Phylogenet. Evol. 52, 520–523. doi: 10.1016/j.ympev.2009.01.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Short, F. T., Kosten, S., Morgan, P. A., Malone, S., Moore, G. E. (2016). Impacts of climate change on submerged and emergent wetland plants. Aquat. Bot. 135, 3–17. doi: 10.1016/j.aquabot.2016.06.006

CrossRef Full Text | Google Scholar

Sultan, S. E., Wilczek, A. M., Bell, D. L., Hand, G. (1998). Physiological response to complex environments in annual Polygonum species of contrasting ecological breadth. Oecologia 115, 564–578. doi: 10.1007/s004420050554

PubMed Abstract | CrossRef Full Text | Google Scholar

Tillich, M., Lehwark, P., Pellizzer, T., Ulbricht-Jones, E. S., Fischer, A., Bock, R., et al. (2017). GeSeq - versatile and accurate annotation of organelle genomes. Nucleic Acids Res. 45, W6–W11. doi: 10.1093/nar/gkx391

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, W., Zhang, X., Garcia, S., Leitch, A. R., Kovarik, A. (2023). Intragenomic rDNA variation - the product of concerted evolution, mutation, or something in between? Heredity (Edinb) 131, 179–188. doi: 10.1038/s41437–023-00634–5

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiao, L. Q., Moller, M., Zhu, H. (2010). High nrDNA ITS polymorphism in the ancient extant seed plant Cycas: incomplete concerted evolution and the origin of pseudogenes. Mol. Phylogenet. Evol. 55, 168–177. doi: 10.1016/j.ympev.2009.11.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, B., Zeng, X. M., Gao, X. F., Jin, D. P., Zhang, L. B. (2017). ITS non-concerted evolution and rampant hybridization in the legume genus Lespedeza (Fabaceae). Sci. Rep. 7, 40057. doi: 10.1038/srep40057

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Chi, X., Zhong, J., Fernie, A., Alseekh, S., Huang, L., et al. (2022). Extensive nrDNA ITS polymorphism in Lycium: Non-concerted evolution and the identification of pseudogenes. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.984579

CrossRef Full Text | Google Scholar