Plastome phylogenomics and historical biogeography of aquatic plant genus Hydrocharis (Hydrocharitaceae)
BMC Plant Biology volume 22, Article number: 106 (2022)
Hydrocharis L. and Limnobium Rich. are small aquatic genera, including three and two species, respectively. The taxonomic status, phylogenetic relationships and biogeographical history of these genera have remained unclear, owing to the lack of Central African endemic H. chevalieri from all previous studies. We sequenced and assembled plastomes of all three Hydrocharis species and Limnobium laevigatum to explore the phylogenetic and biogeographical history of these aquatic plants.
All four newly generated plastomes were conserved in genome structure, gene content, and gene order. However, they differed in size, the number of repeat sequences, and inverted repeat borders. Our phylogenomic analyses recovered non-monophyletic Hydrocharis. The African species H. chevalieri was fully supported as sister to the rest of the species, and L. laevigatum was nested in Hydrocharis as a sister to H. dubia. Hydrocharis-Limnobium initially diverged from the remaining genera at ca. 53.3 Ma, then began to diversify at ca. 30.9 Ma. The biogeographic analysis suggested that Hydrocharis probably originated in Europe and Central Africa.
Based on the phylogenetic results, morphological similarity and small size of the genera, the most reasonable taxonomic solution to the non-monophyly of Hydrocharis is to treat Limnobium as its synonym. The African endemic H. chevalieri is fully supported as a sister to the remaining species. Hydrocharis mainly diversified in the Miocene, during which rapid climate change may have contributed to the speciation and extinctions. The American species of former Limnobium probably dispersed to America through the Bering Land Bridge during the Miocene.
Hydrocharis L., an aquatic monocot genus, has a broad native distribution in tropical and temperate regions of the Old World . As traditionally circumscribed, the genus comprises only three allopatric species: Hydrocharis chevalieri (De Wild.) Dandy, H. morsus-ranae L. and H. dubia (Blume) Backer. Among them, H. chevalieri is only found in small sedge swamps or wetlands in central Africa [1, 2]. It is easily distinguished from the other two species by unique morphological features, e.g., long and erect petioles and the number of primary veins (> 4) in leaf lamina . In comparison, H. morsus-ranae and H. dubia are more widely distributed. The native range of the former covers West and North Eurasia, whereas the latter is natively distributed in Southeast Asia [1, 3]. Beyond their natural range, H. morsus-ranae and H. dubia have invaded Northeast America and Australia, respectively , and cause serious damage to the local environment . However, due to the deterioration of the aquatic habitats resulting from human activities, H. morsus-ranae and H. dubia have been considered threatened in part of their native distributions . Hydrocharis and the closely related Limnobium are unusual in Hydrocharitaceae by having aerial leaves only, and therefore are important for understanding the evolution of aquatic adaptations in the family [1, 4]. Additionally, H. morsus-ranae has been shown to have the ability to accumulate heavy metal elements making it a possible bioindicator of tracing element pollution in freshwaters .
Since the taxonomic revision of the genus Hydrocharis by Cook and Lüönd , only a few studies have focused on the systematics and taxonomy of the genus. Chen et al.  reappraised the generic relationships of Hydrocharitaceae using eight genes representing nuclear, plastid, and mitochondria genomes. The genus Hydrocharis was resolved as monophyletic and highly supported sister of Limnobium, a genus with two currently accepted species . In that study, both Limnobium species were sampled, but one of the Hydrocharis species remained unsampled. Similarly, Ross et al.  resolved Hydrocharis and Limnobium as sisters in a phylogenomic analysis of 83 plastid genes, but they only sampled one species of each genus. Bernardini and Lucchese  further studied the phylogeny of Hydrocharitaceae using a dataset of ITS and five plastid genes. In their study, the taxon sampling and results concerning Limnobium and Hydrocharis were comparable with Chen et al. . They suggested that these two genera could be merged into a single genus. This has also been supported by morphologic evidence [6, 9, 10] and recently Plant Gateway  concluded – without further justification – that Limnobium is part of Hydrocharis. However, the systematic relationships of Hydrocharis and Limnobium remain unclear due to the lack of morphologically intermediate Central African endemic H. chevalieri from all previous studies.
Previous biogeographic analyses suggested that the most recent common ancestor (MRCA) of the Hydrocharis-Limnobium clade originated in the Oriental region approximately 15.9 Ma, and long distance dispersal (LDD) might contribute to its current distributions . However, numerous fossils, including more than 10 Hydrocharis species, which mainly occurred during the Oligocene and Miocene epochs, were discovered in Europe and the Far East [3, 12,13,14,15]. No fossils have been found in Africa. The biogeographic history of the African species H. chevalieri has been unsolved, limiting a comprehensive understanding of the biogeographic history of the genus.
As the plastome has become more easily retrievable, it has dramatically enhanced the resolution of phylogenetic relationships at various taxonomic levels [7, 16,17,18]. Here, we sequenced and assembled plastomes of all three Hydrocharis species and Limnobium laevigatum. The objectives of this study were to 1) investigate the plastome evolution in this group; 2) clarify the phylogenetic relationships of Hydrocharis and Limnobium; and 3) infer the historical biogeography of the group.
After de novo assembly, we obtained the complete plastome of three Hydrocharis species and L. laevigatum from genome skimming data. The coverage of each plastome ranged from ~ 807× to ~ 878× (Table 1). The size of newly assembled plastomes ranged from 153,373 bp to 159,698 bp and exhibited the typical quadripartite structure (Fig. 1), including a pair of IRs (22,292–30,605 bp) separated by the Large Single Copy (LSC, 85,578–89,581 bp) and Small Single Copy (SSC, 11,278–21,476 bp) regions (Table 1). Except for rps16, lost from the plastome of H. chevalieri, the other plastomes encoded 113 unique genes, which is similar to most of the plastomes in angiosperms. Unique genes comprised 79 PCGs, 30 tRNA genes, and four rRNA genes. Both the average gene density (~ 0.8) and the overall GC content (~ 37%) of plastomes were conserved across the species. Furthermore, a small inversion of the trnN-GUU gene was detected in H. dubia plastome (Fig. 1).
Comparison of border regions and sequence identity
The IR/LSC and IR/SSC junctions were compared to assess the IR expansion or contraction of the plastome. Here, the LSC/IRa (JLA) and LSC/IRb (JLB) junctions were identical for all four species; IR expanded into the rpl23 gene at the JLB, and the JLA were located at the trnI-trnH intergenic spacer (Fig. 2). substantial divergence was detected at the SSC/IRa (JSA) and SSC/IRb (JSB) junctions. In H. dubia IR expanded into the ycf1 gene at the JSA, whereas in H. chevalieri the ndhA gene overlapped with the JSA and in H. morsus-ranae and L. laevigatum JSA was placed in ycf1-trnN intergenic spacer. In comparison, at the JSB the IR regions expanded into genes in H. dubia (ycf1), H. morsus-ranae (ndhF) and L. laevigatum (ndhF). The JSB of H. chevalieri lay in the ndhH-ndhF intergenic spacer.
Both the sequence identity analysis and Pi values revealed a high differentiation among the newly assembled plastomes. The divergent regions were mostly derived from IGS regions (Fig. S1,S2), e.g., ycf1-trnN-GUU, rps16-trnQ-UUG, trnI-CAU-ycf2, ycf2-trnI-CAU, trnK-UUU-rps16. Additionally, the top five PCGs, ycf1, rps12, rps18, rps3, accD, showed relatively high nucleotide diversity and had the potential to be developed into molecular markers for further phylogenetic research.
Sequence repeat and SSR analysis
The SSR analysis revealed a similar number of SSRs (Table S1) , ranged from 64 (H. morsus-ranae) to 69 (H. dubia), within Hydrocharis species. In contrast, a clearly higher number of SSRs was detected in L. laevigatum (85). Moreover, species of Hydrocharis had the largest number of tri-, tetra-, and hexa-nucleotide repeats, whereas L. laevigatum was rich in mono-, and di- nucleotide repeats were more common in L. laevigatum (Table S1). Four types of repeat elements, forward, palindromic, reverse, and complementary, were identified for Hydrocharis and L. laevigatum plastomes. The total number varied from 383 (H. morsus-ranae) to 689 (H. dubia) repeats of 30–317 bp length for the newly generated plastomes (Fig. S3). Forward repeats were the most common and complementary repeats were the rarest repeat type. Based on the repeat size, we further divided all repeats into the following categories: 35–45 bp, 46–60 bp, 61–75 bp, 76–90 bp, and > 90 bp (Fig. 3). At least half of the repeats belonged to the size class of 35–45 bp, and the majority of long repeats (> 90 bp) were located in the plastome of H. dubia (107 repeats).
Both ML and BI methods, regardless of partitioning strategy, revealed completely identical topologies and robust support for most of the nodes (Fig. 3). Two major clades were identified with high support (BS = 99/100, PP = 1/1) and the genus Stratiotes was fully resolved as the sister to these two clades. Hydrocharis and Limnobium were clustered with full support (BS = 100/100, PP = 1/1) and resolved as the sister group of the remaining genera within Clade I. However, Hydrocharis did not form a monophyletic group, but L. laevigatum was nested within it as a fully supported sister to H. dubia. Additionally, African species H. chevalieri was fully supported as sister to the rest of the species within (Hydrocharis + Limnobium) group, and H. morsus-ranae shared a common ancestor with (L. laevigatum + H. dubia).
Divergence time estimation and biogeographical reconstruction
The MRCA of Hydrocharitaceae was inferred at ca. 61.7 Ma [Fig. 3, 95% highest posterior densities (HPD): 54.6–71.8 Ma] and the Clade I and Clade II split from each other at 57.6 Ma (95% HPD: 49.5–67.4 Ma). In Clade I, Hydrocharis-Limnobium initially diverged from the remaining genera at ca. 53.3 Ma (95% HPD: 43.7–63.1 Ma), and then began to diversify at ca. 30.9 Ma (HPD: 13.8–49.9 Ma). The MRCA of H. dubia and L. laevigatum was dated at ca. 12.84 Ma (95% HPD: 3.8–24.1 Ma) and diverged of H. morsus-ranae and the remaining species at ca. 17.5 Ma (95% HPD: 6.6–30.5 Ma).
Biogeographical reconstruction was implemented using RASP  with DIVALIKE model supported as the best-fit model (LnL = − 8.23, AICc = 24.46, AICc_wt = 0.56; Table 2) from BioGeoBEARS analysis. However, the ancestral area remained uncertain. The scenario with the highest probability (39%) suggested that the genus originated in Europe and Central Africa, followed by diversification due to three vicariance and two dispersal events (Fig. 3).
Similar to the most reported plastomes of angiosperms [17, 20], all the plastomes in our study exhibited relatively conservative genome structure, gene content, and gene order [21, 22]. However, when compared to H. morsus-ranae and L. laevigatum, H. dubia (159,698 bp) and H. chevalieri (158,066 bp) had larger plastomes. One reason for plastome size variation is the expansion or contraction of IRs [23, 24]. Indeed, H. dubia and H. chevalieri had one (ycf1) and three more genes (ycf1, rps15, and ndhH) in IRa, respectively, than the other two species. In addition, typical patterns of the plastome evolution associated with the fluctuation of plastome length include the gain/loss of genes, pseudogenization, and variations in intergenic regions [24,25,26]. Our results revealed the loss of gene rps16 in African endemic H. chevalieri. This gene has been reported missing from plastomes in many plant groups, e.g., Oxalidaceae [27, 28], Podostemaceae , and Violaceae . It seems that rps16 has been frequently transferred to the nuclear genome [20, 30]. Additionally, a small inversion of the trnN-GUU gene near the IR boundary was detected in H. dubia. The inversion was likely related to the shift in the IR boundary (Fig. 2), as identified and demonstrated by Zhu et al. .
Previous studies [4, 8] have applied only five plastid genes (rbcL, matK, rpoB, rpoC1, trnK) to resolve the phylogeny with the genus Hydrocharis. However, there was only low nucleotide diversity(Pi < 0.06) within these five genes, which might have affected the accuracy and resolution of phylogenetic reconstruction . In our study, a number of regions (Fig. S2), including five plastid genes (ycf1, rps12, rps18, rps3, accD) and five intergenic spacers (ycf1-trnN-GUU, rps16-trnQ-UUG, trnI-CAU-ycf2, ycf2-trnI-CAU, trnK-UUU-rps16), were detected as hotpots and may be helpful even in population genetic studies in the genus.
The high frequency of repeat sequences has been demonstrated to be one of the leading causes of plastome rearrangement and divergence . Our analysis identified many repetitions (> 350) with more than 30 bp length in all newly assembled plastomes. Among these, short repeat sequences (30–45 bp) were dominant, similar to other plant plastomes that have not undergone large-scale structural variation. SSRs are widely distributed in plant plastomes and exhibit relatively high polymorphism, which can be used in population genetics [34,35,36]. The majority of SSRs were mono- and di- nucleotides in four newly assembled plastomes, which have been reported in other angiosperm plastomes, e.g., Primula , Dendrosenecio , and Oxalis . The cpSSRs reported here could be used as genetic markers for future studies into the genetic diversity of Hydrocharis.
Phylogeny and biogeographical reconstruction
Prior to this study, no molecular information was available for the African endemic H. chevalieri. Here we have clarified the relationships within Hydrocharis by assembling and reporting the plastomes of all three Hydrocharis species and Limnobium laevigatum. Unlike previous studies [4, 8], H. morsus-ranae was fully supported (Fig. 3, BS = 100/100; PP = 1/1) as sister to the (L. laevigatum + H. dubia), instead of getting together with H. dubia, which was likely due to the inadequacy of L. spongia in our current study. Moreover, we recovered non-monophyletic Hydrocharis. The genus Limnobium represented here by L. laevigatum was nested in Hydrocharis as sister to H. dubia with robust support (Fig. 3, BS = 100/100; PP = 1/1). This contradicting result might be because previous studies did not include H. chevalieri and used only a limited number of molecular markers [4, 8]. Furthermore, a series of morphological features support the current phylogenomic relationships (Fig. 3). Vegetatively, Hydrocharis and Limnobium are indistinguishable, but H. chevalieri has stout, erect petioles, and laminas with a large number of primary veins, some of which originate from the lower half of the midrib . The remaining species have slender petioles and primary veins originating from the point of petiolar attachment . Limnobium differs from Hydrocharis by its rudimentary petals . Given the morphological similarity and small size of the genera, only two species in Limnobium and three species in Hydrocharis, the most reasonable taxonomic solution to the non-monophyly of Hydrocharis is to treat Limnobium as the synonym of Hydrocharis. Nomenclatural combinations already exist, i.e., Hydrocharis spongia Bosc and Hydrocharis laevigata (Willd.) Byng & Christenh.
Our time-calibrated tree indicated stem node age of about 53 Ma for Hydrocharis (Fig. 3), which is similar to the previously reported median age of 54.7 Ma  and in agreement with the oldest known Hydrocharis fossils from the Eocene . The crown age of Hydrocharis was estimated to be approximately 31 Ma, much older than the 15.9 Ma reported in Chen et al.  based on analysis that lacked H. chevalieri, the sister species of the remaining genus. Based on our results, Hydrocharis mainly diversified in the Miocene, which is consistent with the time inferred for many other extant plant species [39, 40]. Rapid climate change during the Miocene may have contributed to the speciation of Hydrocharis as well as extinctions, given that the fossil diversity is high in comparison to the limited extant diversity [41, 42]. Additionally, paleoclimatic changes may have induced such changes in water bodies that some groups were driven to extinction . This might explain why there are only three extant species, and a majority of fossil groups only occurred in Europe during the Miocene .
The area of origin remained uncertain for the genus because several alternative areas with low probabilities were recovered for the deeper nodes (Fig. 3). This may, of course, also reflect a widespread ancestral distribution. Although, L. spongia from America was not sampled here, our model with the highest probabilities for the ancestor of the genus in Europe, Central Asia, and Central African regions (Fig. 3) contrasted the biogeographical model of Chen et al. , which suggested origination in the Oriental area. Additionally, our results indicated that at least three vicariance and two dispersal events had shaped the current distribution of the genus (Fig. 3). One dispersal event was from Central Africa to Europe and Central Asia, probably via the last direct connection between Eurasia and Africa prior to the Miocene. The rapid climate change during the Miocene likely resulted in the isolation of H. chevalieri, as well as contributed to the extinction of many Hydrocharis species in Europe, as also proposed for Pyrularia , Allioideae  and some Compositae groups . And other one dispersal route from Europe and Central Asia to America, East and Southeast Asia. The Bering Land Bridge (BLB) has been proved to play an important role in the dispersal of different plant lineages between Eurasia and America [47,48,49]. The former genus Limnobium, naturally recorded only in America, may have dispersed to America through the BLB during the Miocene [49, 50]. The LDD event might be interpreted by the transport of seeds or stolons in mud on the feet of aquatic birds, and It has been considered a credible explanation for the LDD of aquatic groups [4, 17]. The geographic isolation between Eurasia and America likely contributed to the isolation of Limnobium. Although the phylogeny and biogeographic history of Hydrocharis has been greatly improved here, it should be mentioned that the incomplete sampling of Limnobium in America and only one individual was used for each species. There is still a need to include L. spongia and more populations to solve a more detailed biogeographic history of Hydrocharis.
Taxon sampling and DNA extraction
All the samples used in this study were listed in Table 1. The fresh leaves of two Hydrocharis species, H. chevalieri and H. dubia, were collected in 2019 from Cameroon and China, respectively. Limnobium laevigatum was sampled from a small swamp in Brazil. The field collection followed the ethics and legality of the local government and was permitted by the government. Total genomic DNA of these three species was extracted from silica-dried leaves using the MagicMag Genomic DNA Micro Kit (Sangon Biotech, Shanghai, China) following the manufacturer’s protocol. The genomic DNA of H. morsus-ranae was retrieved from the Kew DNA bank (ID: 14513) and used for the subsequent analysis.
Genome skimming, plastome assembly, and annotation
The construction of all sequencing libraries and genome skimming was carried out following Li et al. . Approximately 8.5G paired-end reads (150 bp) were yielded for each sample. After trimming and filtering using Fastp v. 0.20.0  with default parameters, at least 7G clean reads were obtained. The complete plastomes were assembled de novo using GetOrganelle v. 1.7.5 with the recommended settings . Web applications Geseq  and PGA  were applied to annotate genes in the newly generated plastomes using default settings, with manual adjustment of the start/stop codons through comparison with a reference plastome (Ottelia alismoides; NC_057494). Circular maps of plastome were created using OGDRAW v1.3.1 . All newly generated sequencing data and annotated plastomes were submitted to the GenBank (Table 1).
Comparative analysis and divergence hotspot identification
The web-based software mVISTA  was used to identify the sequence and structural variations in Hydrocharis and Limnobium, using the plastome of H. dubia as reference. The comparison of expansions/contraction of IR in Hydrocharis and Limnobium was conducted using Geneious v 5.6.3 .
All the protein coding genes (PCGs) and intergenic spacers (IGSs) were extracted separately, and aligned in MAFFT v. 7.221  with default settings. The program DnaSP v.6  was used to assess the nucleotide diversity (Pi) for all PCGs and IGS of the studied plastomes. Additionally, the online application MISA  was employed to predict the simple sequence repeats (SSRs) for each plastome. The minimum number of repetitions was set to 10, 5, 4, 3, 3, and 3, for mono-, di-, tri-, tetra-, penta-, and hexanucleotides, respectively. Furthermore, we explored the forward, reverse, palindromic, and complementary repeats for each plastome using REPuter web-based tools , with a minimum repeat size of 30 bp and Hamming distance equal to 3.
Six plastomes of Hydrocharitaceae (Elodea canadensis NC_018541; Enhalus acoroides NC_048519; Halophila beccarii NC_051970; Najas flexilis NC_021936; Ottelia cordata NC_046891; Thalassia hemprichii NC_043774) were downloaded from GenBank. The PCGs were extracted from each of the plastomes and used in the phylogenetic analysis. A total of 78 PCGs in 15 Hydrocharitaceae species, three Alismataceae species (Alisma triviale, Baldellia ranunculoides, and Hydrocleys nymphoides), and Butomus umbellatus representing Butomaceae were retrieved from Ross et al. . Finally, we compiled all 78 PCGs from 29 samples for our phylogenetic analysis (Table S2). Three Alismataceae species and B. umbellatus were used as outgroups.
All PCGs were aligned using MAFFT v. 7.221  and ambiguously aligned regions were removed using trimAl v. 1.2 . Unpartitioned and partitioned strategies of both Maximum likelihood (ML) and Bayesian inference (BI) methods were utilized for phylogenetic inference. For the unpartitioned strategy, we concatenated all 78 PCGs sequences as a supermatrix to infer the phylogenetic relationship with the best-fit model of nucleotide substitution estimated by ModelFinder . For partitioned strategy, PartitionFinder 2  was used to select the best-fit partitioning schemes and models for 78 PCGs via the rcluster algorithm. All phylogenetic analyses were performed with ML using IQ-Tree v. 1.6.12  and BI using MrBayes v. 3.2.7 . An ultrafast bootstrap approximation was used to estimate ML branch support values with 1000 replicates. BI was run using two independent runs of 8 million generations, and four Markov chains were set for each run, sampling trees every 5000 generations. After verifying that an average standard deviation of split frequencies was less than 0.001, the initial 25% of the trees were discarded as burn-in. A consensus tree with Bayesian posterior probabilities (PP) was constructed from the remaining trees.
Molecular dating and ancestral area reconstruction
Three credible fossils from Hydrocharitaceae with Lognormal prior distributions were employed to calibrate the divergence times of the lineages within Hydrocharitaceae, following previous studies [4, 17, 67]. The first calibration point was assigned between Stratiotes and the rest of the Hydrocharitaceae based on the fossil seeds of Stratiotes from the Paleocene-Eocene boundary . This node was constrained to a minimum of 54.5 Ma (offset = 54.5, mean = 1.0, SD = 1.0). The fossil Thalassites was applied to constrain the stem node of Enhalus + Halophila + Thalassia to the Middle Eocene. Thus, the split between (Enhalus + Halophila + Thalassia) and (Nechamandra + Vallisneria) was constrained to a minimum of 38.0 Ma (offset = 38.0, mean = 0.8, SD = 0.9). The oldest fossil from Ottelia recorded in the Upper Eocene was used to constrain the split between Blyxa and Ottelia with a minimum age of 33.7 Ma (offset = 33.7, mean = 1.1, SD = 1.2). An uncorrelated Lognormal relaxed clock model with Yule tree prior was applied for molecular dating implemented in BEAST v.1.10.4 . We ran 8.0 × 108 iterations of Markov chain Monte Carlo (MCMC) and sampled every 2000 iterations. The program Tracer v. 1.7.1  was used to check the effective sample size (ESS) for the convergence of each parameter discarding the initial 25% generations as burn-in.
According to the existing distribution patterns of Hydrocharis [1, 3], four major areas were identified for biogeographic analysis: A) East and Southeast Asia; B) Europe and Central Asia; C) Central Africa; D) America. The ancestral states were reconstructed using the BioGeoBEARS package  implemented in RASP v. 4.0 . The best fit biogeographic model was selected by BioGeoBEARS  based on the Akaike Information Criterion cumulative weight (AICc_wt). The condensed tree and 100,000 sampled trees from BEAST analysis were used as input.
In this study, we firstly assembled plastomes of three Hydrocharis species and Limnobium laevigatum. Our phylogenomic analysis recovered non-monophyletic Hydrocharis and fully supported H. chevalieri as a sister to the rest of the species within the (Hydrocharis + Limnobium) group. Based on the morphological and phylogenomic evidence, we suggested treating Limnobium as the synonym of Hydrocharis. Moreover, we reappraised the divergence time and historical biogeography of Hydrocharis. Hydrocharis mainly diversified in the Miocene, rapid climate change during the Miocene may have contributed to the speciation of Hydrocharis and extinctions. The American members of the genus, i.e. former Limnobium, probably have dispersed to America through the Bering Land Bridge during the Miocene. In summary, our study provides new insights into the plastome evolution, phylogeny, and biogeography of the genus Hydrocharis.
Availability of data and materials
All newly annotated plastomes in this study are available from the National Center for Biotechnology Information (NCBI) (accession numbers: OK326868- OK326871). The associated BioProject, BioSample numbers, and SRA are PRJNA771243, SAMN22263148-SAMN22263151, and SRR16469678- SRR16469681.
Akaike Information Criterion cumulative weight
Effective sample size
Internal transcribed spacer
Long distance dispersal
Large single copy
Markov chain Monte Carlo
Most recent common ancestor
Protein coding genes
Small single copy
Simple sequence repeat
Cook CDK, Lüönd R. A revision of the genus Hydrocharis (Hydrocharitaceae). Aquat Bot. 1982;14:177–204.
Bean AR. Hydrocharis dubia (Blume) backer (Hydrocharitaceae) is an alien species in Australia. Austrobaileya. 2011;8:435–7.
Efremov AN, Grishina VS, Kislov DE, Mesterházy A, Toma C. The genus Hydrocharis L. (Hydrocharitaceae): distribution features and conservation status. Bot Pac. 2020;9:83–94.
Chen LY, Chen JM, Gituru R, Wang QF. Generic phylogeny, historical biogeography and character evolution of the cosmopolitan aquatic plant family Hydrocharitaceae. BMC Evol Biol. 2012;12:30.
Polechońska L, Klink A. Validation of Hydrocharis morsus-ranae as a possible bioindicator of trace element pollution in freshwaters using Ceratophyllum demersum as a reference species. Environ Pollut. 2021;269:116145.
Cook CDK, Urmi-König K. A revision of the genus Limnobium including Hydromistria (Hydrocharitaceae). Aquat Bot. 1983;17:1–27.
Ross TG, Barrett CF, Soto-Gomez M, Lam VKY, Henriquez CL, Les DH, et al. Plastid phylogenomics and molecular evolution of Alismatales. Cladistics. 2016;32:160–78.
Bernardini B, Lucchese F. New phylogenetic insights into Hydrocharitaceae. Ann Bot. 2018;8:45–58.
Shaffer-Feher M. The endotegmen tuberculae: an account of little-known structures from the seed coat of the Hydrocharitoideae (Hydrocharitaceae) and Najas (Najadaceae). Bot J Linn Soc. 1991;107:169–88.
Les DH, Moody ML, Soros CL. A reappraisal of phylogenetic relationships in the monocotyledon family Hydrocharitaceae (Alismatidae). Aliso. 2006;22:211–30.
Gateway P. A practical Flora to vascular plant species of the world. Vol. 4 of The Global Flora. Special ed. Bradford: GLOVAP, Nomenclature, Part 1; 2018.
Mai DH, Walther H. Die pliozänen Floren von Thü ringen. Quartärpaläontologie. 1988;7:55–295.
Carrión JS, Dupré M. Late Quaternary vegetational history at Navarrés, eastern Spain. A two-core approach. New Phytol. 1996;134:177–91.
Velichkevich FY, Zastawniak E. The Pliocene flora of Kholmech, South-Eastern Belarus and its correlation with other Pliocene floras of Europe. Acta Palaeobot. 2003;43:137–259.
Yao YF, Bruch AA, Mosbrugger V, Li CS. Quantitative reconstruction of Miocene climate patterns and evolution in southern China based on plant fossils. Palaeogeogr Palaeoclimatol Palaeoecol. 2011;304:291–307.
Li P, Lu RS, Xu WQ, Ohi-Toma T, Cai MQ, Qiu YX, et al. Comparative genomics and phylogenomics of east Asian tulips (Amana, Liliaceae). Front Plant Sci. 2017;8:1–12.
Li ZZ, Lehtonen S, Martins K, Gichira AW, Wu S, Li W, et al. Phylogenomics of the aquatic plant genus Ottelia ( Hydrocharitaceae ): implications for historical biogeography. Mol Phylogenet Evol. 2020;152:106939.
Gichira AW, Chen LY, Li ZZ, Hu GW, Saina JK, Gituru RW, et al. Plastid phylogenomics and insights into the inter-mountain dispersal of the eastern African giant senecios (Dendrosenecio, Asteraceae). Mol Phylogenet Evol. 2021;164(11):107271.
Yu Y, Blair C, He XJ. RASP 4: ancestral state reconstruction tool for multiple genes and characters. Mol Biol Evol. 2020;37:604–6.
Mwanzia VM, Nzei JM, Yan DY, Kamau PW, Chen JM, Li ZZ. The complete chloroplast genomes of two species in threatened monocot genus Caldesia in China. Genetica. 2019;147:381–90.
Ravi V, Khurana JP, Tyagi AK, Khurana P. An update on chloroplast genomes. Plant Syst Evol. 2008;271:101–22.
Wicke S, Schneeweiss GM, dePamphilis CW, Müller KF, Quandt D. The evolution of the plastid chromosome in land plants: gene content, gene order, gene function. Plant Mol Biol. 2011;76:273–97.
Jansen RK, Ruhlman TA. Plastid Genomes of Seed Plants. In: In Genomics of Chloroplasts and Mitochondria. Dordrecht: Springer; 2012. p. 103–26.
Mower JP, Vickrey TL. Structural diversity among plastid genomes of land plants. Adv Bot Res. 2018;85:263–92.
Jansen RK, Cai Z, Raubeson LA, Daniell H, dePamphilis CW, Leebens-Mack J, et al. Analysis of 81 genes from 64 plastid genomes resolves relationships in angiosperms and identifies genome-scale evolutionary patterns. Proc Natl Acad Sci U S A. 2007;104:19369–74.
Ruhlman TA, Jansen RK. The plastid genomes of flowering plants. In: Maliga P, editor. Chloroplast biotechnology, methods in molecular biology (methods and protocols), vol. 1132. Totowa, NJ: Humana Press; 2014. p. 3e38.
Schmickl R, Liston A, Zeisek V, Oberlander K, Weitemier K, Straub SCK, et al. Phylogenetic marker development for target enrichment from transcriptome and genome skim data: the pipeline and its application in southern African Oxalis (Oxalidaceae). Mol Ecol Resour. 2015;16:1124–35.
Li X, Zhao Y, Tu X, Li C, Zhu Y, Zhong H, et al. Comparative analysis of plastomes in Oxalidaceae: phylogenetic relationships and potential molecular markers. Plant Divers. 2021;43:281–91.
Bedoya AM, Ruhfel BR, Philbrick CT, Madriñán S, Bove CP, Mesterházy A, et al. Plastid genomes of five species of riverweeds (Podostemaceae): structural organization and comparative analysis in Malpighiales. Front Plant Sci. 2019;10:1–14.
Sudianto E, Chaw S. Two independent plastid accD transfers to the nuclear genome of Gnetum and other insights on acetyl-CoA carboxylase evolution in gymnosperms. Genome Biol Evol. 2019;11:1691–705.
Zhu AD, Guo WH, Gupta S, Fan WS, Mower JP. Evolutionary dynamics of the plastid inverted repeat: the effects of expansion, contraction, and loss on substitution rates. New Phytol. 2016;209:1747–56.
Som A. Causes, consequences and solutions of phylogenetic incongruence. Brief Bioinform. 2015;16:536–48.
Weng ML, Blazier JC, Govindu M, Jansen RK. Reconstruction of the ancestral plastid genome in Geraniaceae reveals a correlation between genome rearrangements, repeats, and nucleotide substitution rates. Mol Biol Evol. 2014;31:645–59.
Provan J, Powell W, Hollingsworth PM. Chloroplast microsatellites: new tools for studies in plant ecology and evolution. Trends Ecol Evol. 2001;16:142–7.
Liu L, Wang Z, Huang L, Wang T, Su Y. Chloroplast population genetics reveals low levels of genetic variation and conformation to the central–marginal hypothesis in Taxus wallichiana var. mairei, an endangered conifer endemic to China. Ecol Evol. 2019;9(20):11944–56.
Kosiński P, Sękiewicz K, Walas L, Boratyński A, Dering M. Spatial genetic structure of the endemic alpine plant Salix serpillifolia: genetic swamping on nunataks due to secondary colonization? Alp Bot. 2019;129(2):107–21.
Ren T, Yang Y, Zhou T, Liu ZL. Comparative plastid genomes of Primula species: sequence divergence and phylogenetic relationships. Int J Mol Sci. 2018;19:1050.
Stockey RA. The fossil record of basal monocots. Aliso. 2006;22:91–106.
Zachos J, Pagani H, Sloan L, Thomas E, Billups K. Trends, rhythms, and aberrations in global climate 65 ma to present. Science. 2001;292:686–93.
Kong H, Condamine FL, Harris AJ, Chen J, Pan B, Möller M, et al. Both temperature fluctuations and east Asian monsoons have driven plant diversification in the karst ecosystems from southern China. Mol Ecol. 2017;26:6414–29.
Crowley TJ, North GR. Abrupt climate change and extinction events in earth history. Science. 1988;240:996–1002.
Shukla A, Mehrotra RC, Guleria JS. Emergence and extinction of Dipterocarpaceae in western India with reference to climate change: fossil wood evidences. J Earth Syst Sci. 2013;122:1373–86.
Alçiçek H, Jiménez-Moreno G. Late Miocene to Plio-Pleistocene fluvio-lacustrine system in the Karacasu Basin (SW Anatolia, Turkey): depositional, paleogeographic and paleoclimatic implications. Sediment Geol. 2013;291:62–83.
Zhou Z, Hu JJ, Wen J, Sun H. Morphometric, phylogenetic and biogeographic analyses of Pyrularia (Santalales), a parasitic disjunct lineage between eastern Asia and eastern North America. Taxon. 2019;68:47–71.
Namgung J, Do HDK, Kim C, Choi HJ, Kim JH. Complete chloroplast genomes shed light on phylogenetic relationships, divergence time, and biogeography of Allioideae (Amaryllidaceae). Sci Rep. 2021;11:1–13.
Keeley SC, Cantley JT, Gallaher TJ. The “evil tribe” spreads across the land : a dated molecular phylogeny provides insight into dispersal, expansion , and biogeographic relationships within one of the largest tribes of the sunflower family ( Vernonieae : Compositae ). Am J Bot. 2021;108:505–19.
Tiffney BH, Manchester SR. The use of geological and paleontological ev- idence in evaluating plant phylogeographic hypotheses in the northern hemisphere tertiary. Int J Plant Sci. 2001;162:S3–S17.
Donoghue MJ, Smith SA. Patterns in the assembly of temperate forests around the northern hemisphere. Philos Trans R Soc Lond Ser B Biol Sci. 2004;2004(359):1633–44.
Jia LB, Manchester SR, Su T, Xing YW, Chen WY, Huang YJ, et al. First occurrence of Cedrelospermum (Ulmaceae) in Asia and its biogeographic implications. J Plant Res. 2015;128:747–61.
Marincovich L, Gladenkov AY. Evidence for an early opening of the Bering Strait. Nature. 1999;397:149–51.
Chen SF, Zhou YQ, Chen YR, Gu J. Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34:884–90.
Jin JJ, Yu WB, Yang JB, Song Y, dePamphilis CW, Yi TS, et al. GetOrganelle: a fast and versatile toolkit for accurate de novo assembly of organelle genomes. Genome Biol. 2020;21:241.
Tillich M, Lehwark P, Pellizzer T, Ulbricht-Jones ES, Fischer A, Bock R, et al. GeSeq-versatile and accurate annotation of organelle genomes. Nucleic Acids Res. 2017;45:6–11.
Qu XJ, Moore MJ, Li DZ, Yi TS. PGA: a software package for rapid, accurate, and flexible batch annotation of plastomes. Plant Methods. 2019;15:50.
Lohse M, Drechsel O, Kahlau S, Bock R. OrganellarGenomeDRAW—a suite of tools for generating physical maps of plastid and mitochondrial genomes and visualizing expression data sets. Nucleic Acids Res. 2013;41:W575–81.
Mayor C, Brudno M, Schwartz JR, Poliakov A, Rubin EM, Frazer KA, et al. VISTA: visualizing global DNA sequence alignments of arbitrary length. Bioinformatics. 2000;16:1046–7.
Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28:1647–9.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.
Rozas J, Ferrer-Mata A, Sanchez-DelBarrio JC, Guirao-Rico S, Librado P, Ramos-Onsins SE, et al. DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol Biol Evol. 2017;34:3299–302.
Beier S, Thiel T, Münch T, Scholz U, Mascher M. MISA-web: a web server for microsatellite prediction. Bioinformatics. 2017;33:2583–5.
Kurtz S, Choudhuri JV, Ohlebusch E, Schleiermacher C, Stoye J, Giegerich R. REPuter: the manifold applications of repeat analysis on a genomic scale. Nucleic Acids Res. 2001;29:4633–42.
Capella-Gutierrez S, Silla-Martinez JM, Gabaldon T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009;25:1972–3.
Kalyaanamoorthy S, Minh BQ, Wong TKF, von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14:587–9.
Lanfear R, Frandsen PB, Wright AM, Senfeld T, Calcott B. Partitionfinder 2: new methods for selecting partitioned models of evolution for molecular and morphological phylogenetic analyses. Mol Biol Evol. 2017;34:772–3.
Nguyen LT, Schmidt HA, von Haeseler A, Minh BQ. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015;32:268–74.
Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61:539–42.
Iles WJ, Smith SY, Gandolfo MA, Graham SW. Monocot fossils suitable for molecular dating analyses. Bot J Linn Soc. 2015;178:346–74.
Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.
Rambaut A, Suchard MA, Xie D, Drummond AJ. Tracer v1.6. 2014. Available at: http://beast.bio.ed.ac.uk/Tracer. Accessed April 20, 2021.
Matzke NJ. Model selection in historical biogeography reveals that founder-event speciation is a crucial process in island clades. Syst Biol. 2014;63:951–70.
We are grateful to Rosie Woods for help in the DNA collection of Hydrocharis morsus-ranae. We appreciate three anonymous reviewers for their valuable suggestions to improve the manuscript. This research was conducted under Permit No. 0000072/MINRESI/B00/C00/C10/C12 in Cameroon. Permission for sampling Limnobium laevigata was issued by the Brazilian Ministry of Environment (ICMBIO/SISBIO 67605-2) and the shipment to Wuhan Botanical Garden was registered in Brazilian SISGEN platform (No. RF546AF).
This work was supported by grants from the Strategic Priority Research Program of the Chinese Academy of Sciences (No. XDB31000000) and the National Natural Science Foundation of China (Nos. 32100186, 32070231).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Genome alignment of plastomes of Hydrocharis species and Limnobium laevigatum. H. dubia was used as a reference.
Nucleotide diversity hotspot regions in plastomes of Hydrocharis species and Limnobium laevigatum.
Analysis of repeat sequences in plastomes of Hydrocharis species and Limnobium laevigatum. (A) Total number of four repeat types. (B) Number of repeats divided by size.
Distribution of SSRs in plastomes of Hydrocharis species and Limnobium laevigatum.
List of the species downloaded from GenBank and Ross et al. (2016) used for phylogenomic and dating analyses in this study.
About this article
Cite this article
Li, ZZ., Lehtonen, S., Gichira, A.W. et al. Plastome phylogenomics and historical biogeography of aquatic plant genus Hydrocharis (Hydrocharitaceae). BMC Plant Biol 22, 106 (2022). https://doi.org/10.1186/s12870-022-03483-2
- Divergence time
- Hydrocharis chevalieri