Generic phylogeny, historical biogeography and character evolution of the cosmopolitan aquatic plant family Hydrocharitaceae
© Chen et al; licensee BioMed Central Ltd. 2012
Received: 4 January 2012
Accepted: 10 March 2012
Published: 10 March 2012
Hydrocharitaceae is a fully aquatic monocot family, consists of 18 genera with approximately 120 species. The family includes both fresh and marine aquatics and exhibits great diversity in form and habit including annual and perennial life histories; submersed, partially submersed and floating leaf habits and linear to orbicular leaf shapes. The family has a cosmopolitan distribution and is well represented in the Tertiary fossil record in Europe. At present, the historical biogeography of the family is not well understood and the generic relationships remain controversial. In this study we investigated the phylogeny and biogeography of Hydrocharitaceae by integrating fossils and DNA sequences from eight genes. We also conducted ancestral state reconstruction for three morphological characters.
Phylogenetic analyses produced a phylogeny with most branches strongly supported by bootstrap values greater than 95 and Bayesian posterior probability values of 1.0. Stratiotes is the first diverging lineage with the remaining genera in two clades, one clade consists of Lagarosiphon, Ottelia, Blyxa, Apalanthe, Elodea and Egeria; and the other consists of Hydrocharis-Limnobium, Thalassia, Enhalus, Halophila, Najas, Hydrilla, Vallisneria, Nechamandra and Maidenia. Biogeographic analyses (DIVA, Mesquite) and divergence time estimates (BEAST) resolved the most recent common ancestor of Hydrocharitaceae as being in Asia during the Late Cretaceous and Palaeocene (54.7-72.6 Ma). Dispersals (including long-distance dispersal and migrations through Tethys seaway and land bridges) probably played major roles in the intercontinental distribution of this family. Ancestral state reconstruction suggested that in Hydrocharitaceae evolution of dioecy is bidirectional, viz., from dioecy to hermaphroditism, and from hermaphroditism to dioecy, and that the aerial-submerged leaf habit and short-linear leaf shape are the ancestral states.
Our study has shed light on the previously controversial generic phylogeny of Hydrocharitaceae. The study has resolved the historical biogeography of this family and supported dispersal as the most likely explanation for the intercontinental distribution. We have also provided valuable information for understanding the evolution of breeding system and leaf phenotype in aquatic monocots.
KeywordsHydrocharitaceae Phylogeny Historical biogeography Dispersal Vicariance Morphological character evolution
Hydrocharitaceae is a fully aquatic monocot family, consists of 18 genera and approximately 120 species [1, 2] with a cosmopolitan distribution. The group is an important component of aquatic ecosystems and occurs in both freshwater and marine habitats. It also exhibits great diversity in form and habit including annual and perennial life histories; submersed, partially submersed and floating leaf habits, and linear to orbicular leaf shapes . It includes a number of valuable aquarium plants and some species (e.g., Hydrilla verticillata and Elodea canadensis) provide fish breeding sites and fodder for the poultry raising industry . Owing to habitat destruction and unprecedented increase in commercial shipping activities, several species (e.g., Ottelia alismoides and Blyxa japonica) are threatened, while some other species (e.g., Hydrilla verticillata) have become invasive weeds of great concern . Similar to many aquatic plants, the group displays considerable convergent adaptations and morphological reduction, which make the study of phylogenetics and taxonomy of the group difficult [1, 5–7]. There exists little consensus in classification systems of the group based on morphological characters [1, 8].
Although molecular phylogenetic studies of Hydrocharitaceae have created consensus on the relationships between certain genera, disagreements on the relationships of other genera still exist. The seagrass subclade which includes Halophila, Enhalus and Thalassia was resolved as sister to the subclade comprising Najas, Hydrilla and Vallisneria by analyses using rbcL+ matK  and rbcL+ cob + atp1 . However, the seagrasses was resolved as sister to the subclade comprising Nechamandra, Vallisneria and Maidenia by analyses using morphological + rbcL+ matK+ trnK intron + ITS . Similarly, Hydrilla, which was reported to be closely related to Vallisneria by analyses using rbcL, matK [2, 8], rbcL, cob and atp1 , was reported to be closely related to Najas by analyses using morphological character + rbcL+ matK+ trnK intron + ITS . Furthermore, Stratiotes has been resolved as sister to the rest of this family by analyses using rbcL  and mitochondrial genes (nad5, cob, ccmB, mtt2, atp1) [9, 10]. The position is not in agreement with the results of two other studies using chloroplast and nuclear sequences [1, 2]. In addition, a single species was selected for each genus in previous studies. The phylogeny of this family could be improved by a better sampling of taxon and DNA sequence, and perhaps a careful outgroup selection .
The divergence time of Hydrocharitaceae is still a subject of debate, and two competing ages (one much more recent than the other) have been proposed. Janssen and Bremer (2004)  placed the crown node age of this family in the Late Cretaceous (75 Ma) by analyses using rbcL and fossil calibrations. Notably, in that study the family was represented by only 16 terminals and there was no internal calibration point. This could be improved by better sampling and by incorporating internal fossil calibration points . Kato et al (2003)  dated the seagrasses within Hydrocharitaceae at 119 ± 11 Ma by analyses using the substitution rates of rbcL and matK. However, this time overlaps with the generally accepted age of the order Alismatales thus putting the validity of the results of that study into doubt . He et al. (1991)  proposed that Ottelia had originated no later than the Cretaceous based on the distribution of the genus and the predictions of the continental drift theory. Fossils of Hydrocharitaceae have been found in Europe including those of the extant genera Vallisneria, Hydrilla, Ottelia, Thalassia, Stratiotes, Hydrocharis and Najas from the Eocene, Oligocene and Miocene [16–18]. The oldest fossil of Hydrocharitacae (genus Stratiotes) is from the Late Paleocene . These fossils have prior to the present study not been incorporated in divergence time estimates.
The geographic origin of Hydrocharitaceae remains unresolved. The diversity centre of the family has been suggested to be in tropical Asia . However, the diversity centre of a taxon does not necessarily correspond to its centre of origin. Numerous fossils of Hydrocharitaceae have been found in Europe , implying a possible European origin of the family. A biogeographic analysis is required to elucidate the geographical origin of the family.
The transoceanic distribution of angiosperms has long intrigued biologists. Two competing hypotheses have been proposed to explain this phenomenon: the first attributing it to dispersal [21–23] and the second to vicariance (continental drift) . Les et al. (2003)  proposed that dispersal is the major factor accounting for the disjunctive distribution of aquatic plants. This is contrary to the traditional viewpoint which considered vicariance  as the major cause of the disjunctive distribution in aquatic taxa such as Limnocharitaceae , Ottelia (Hydrocharitaceae)  and Sagittaria (Alismataceae) . Hydrocharitaceae exhibits a wide transoceanic distribution at genera and species levels. The family can serve as a suitable model to investigate transoceanic distribution in aquatic monocots.
About 10% of all flowering plants have unisexual flowers which have traditionally been regarded as a derived state in angiosperms . Most species of Hydrocharitaceae are unisexual while some are hermaphrodite. Hermaphroditism is regarded as the ancestral condition which gave rise to unisexual flowers . However, hermaphroditic flowers also occur in more recently evolved genera such as Ottelia . This suggests that the view that hermaphroditism is the ancestral state needs to be re-examined. Leaf habit in Hydrocharitaceae varies from aerial, aerial-submerged to fully submerged, and leaf shape varies from circular, linear, to ribbon like . Ancestral state reconstruction is useful in understanding the evolutionary history of reproductive system and leaf morphology and their significance in adapting the plants to the aquatic environment.
In this study, we sampled 17 genera of Hydrocharitaceae, using DNA sequences from 8 genes. The aims of our study were: 1) to reappraise the generic relationships of Hydrocharitaceae; 2) to estimate the divergence times; 3) to investigate the area of origin and the major factor(s) shaping the global distribution; and 4) to investigate the evolution of reproductive system and leaf morphology in this family.
Details of the voucher and DNA sequence information are provided in Additional file 1. We successfully generated 117 sequences. Other sequences used in this study were downloaded from GenBank. Seventeen of the 18 genera in Hydrocharitaceae were sampled. Genus Appertiella was not included because of unavailability of molecular data. The supermatrix dataset, which resulted from assembling the DNA sequences of the 8 genes, was 8,086 nt in length (2,035 mt, 4,406 chl, 1,645 nu). The dataset consisted of 38 terminals of Hydrocharitaceae and was submitted to TreeBASE (accession number S12110). The dataset comprised about 29% missing characters mainly due to unavailability of some sequences.
Dating analysis and ancestral area reconstruction
All BEAST MCMC runs yielded sufficient effective sample sizes (> 200) for all relevant parameters and converged on topologies identical to the tree in Figure 1. The crown node age of Hydrocharitaceae was dated at 65.2 Ma (95% HPD: 54.6-79.6 Ma). The mean divergence between clade A and clade B, estimated to be 63.1 Ma. For the three calibration nodes, mean posterior estimates were consistent with prior node ages, suggesting that the calibration points were sufficiently concordant [29, 30].
Results of the biogeographic analyses using species as terminal taxa suggested that Hydrocharitaceae originated in Orient (Figure 2c). A minimum of 76 dispersal events was inferred from DIVA to explain the current distribution of Hydrocharitaceae. The ancestor of Stratiotes was suggested to have been in Orient and dispersed to Europe during the Late Cretaceous and Paleocene (Figure 2c, d, arrow 1). This route is similar to that which has been reported for Alangiopollis (Alangiaceae) . The analyses using genera or species as terminal taxa yielded comparable results on the ancestral area of Hydrocharitaceae. This indicates that incomplete sampling may have little effect on the accuracy of investigation into the geographical origin of Hydrocharitaceae. However, the analyses using species as terminal successfully resolved the ancestral areas for more genera than those using genera as terminal.
The MRCAs of clade A, the seagrass genera and the seagrass subclade were shown to be of Oriental origin. The taxa then diversified in Oriental region (DIVA, Mesquite; Figure 2c). DIVA suggested that Vallisneria originated in Oriental and Australasian regions. Mesquite suggested an Oriental origin for Vallisneria, followed by diversification and dispersal to other continents (Figure 2c). Long distance dispersal (LDD) for Vallisneria from Oriental area to Africa and Europe was inferred from V. spinulosa, V. spiralis and V. denseserrulata (Figure 2c, e, arrow 1 & 2). LDD from Australasia to North America for this genus was inferred from V. neotropicalis and V. americana (Figure 2c, e, arrow 7). This route is similar to that envisaged for the taxon Scaevola (Goodeniaceae) . LDD from Australasia to Asia similar to that which has been recorded for the plant family Cucurbitaceae  was inferred from taxa including V. asiatica and V. natans (Figure 2c, e, arrow 4). The genus Najas was inferred to have originated in Oriental area during the Oligocene. LDD from Asia to North America (inferred from N. gracillima), to Europe (inferred from N. minor), to Africa (inferred from N. minor) and to Australia (inferred from subclade N. browniana + N. tenuifolia) was suggested (Figure 2c, e, arrow 6, 1, 2 & 5).
The ancestor of clade B was inferred to have dispersed from Orient to Afrotropical region (Mesquite) or Southern hemisphere (DIVA) during the Eocene (Figure 2c, d, arrow 2), followed by diversification in Southern hemisphere during the Tertiary. MRCA of Ottelia was suggested to have lived in Oriental and Afrotropical regions (Figure 2c). LDD from West Africa into South America for this genus was inferred from O. brasiliensis (Figure 2c, e, arrow 9). This route is similar to the one suggested for the dispersal of Gossypium (Malvaceae)  and is further supported by the fact that Ottelia in South America is confined to the southeastern area . LDD from S.E. Asia to Australasia was inferred from O. ovalifolia (Figure 2c, e, arrow 5). Dispersal from Africa into Asia was also inferred from O. cordata and O. mesenterium (Figure 2c, e, arrow 3). This route is similar to that envisaged for the two genera Coccinia and Momordica in family Cucurbitaceae . The MRCA of the subclade comprising Apalanthe, Egeria and Elodea was suggested to have arisen in Neotropical region, while the ancestor of this subclade may have come from Afrotropical region (Figure 2c, e, arrow 9). Dispersal from South America to North America during the Miocene was inferred for the ancestor of Elodea, probably via the island chains of Central America (Figure 2c, e, arrow 10).
Ancestral character state reconstructions
By merging diverse sequences into a supermatix data set, we obtained a well-resolved phylogeny with most branches strongly supported by BT values greater than 95% and PP values of 1.0. This indicates that increasing the number of taxa and the number of molecular markers improved the resultant phylogeny, and it further supports the notion that a supermatrix can be used to obtain a well-resolved and strongly supported phylogeny in cases where some data are missing [37–39]. The phylogenetic relationships of Egeria, Elodea, Ottelia, Blyxa, Apalanthe and Lagarosiphon have remained unchanged in all the earlier molecular phylogenetic studies [1, 2, 8, 9, 40, 41]. Our analyses resolved the generic relationships that are largely similar to those reported in those studies. However, incongruences existed for the other genera, which we briefly address here below.
The genus Hydrilla comprises only one species, H. verticillata. Based on rbcL, matK, trnK intron, ITS and morphological data, Les et al. (2006)  suggested that Hydrilla was most closely related to Najas, despite their being quite divergent at the phenotypic level. Our phylogenetic analyses suggested that Hydrilla is most closely related to the subclade comprising Nechamandra, Vallisneria and Maidenia (Figure 1). This position is consistent with all previous phylogenetic studies (except Les et al., 2006) based on molecular and morphological data [2, 8, 9, 40, 41]. A close relationship between the subclade (Hydrilla (Nechamandra (Vallisneria + Maidenia)) and Najas was strongly supported (BS = 96, PP = 1.0, Figure 1). This is in agreement with earlier results from rbcL+ matK+ trnK intron analyses with ML approach , rbcL , and rbcL+ matK . However, the results did not support the close affinity between the subclade and seagrasses which had been inferred from rbcL  and rbcL+ matK+ trnK intron + ITS analysis .
Stratiotes was resolved as the first diverging lineage within Hydrocharitaceae (BS = 100, PP = 1.0; Figure 1). This is in agreement with analyses based on rbcL , mtt2 and nad5  and the fossil records of Stratiotes which include the most abundant and the oldest fossils of' Hydrocharitaceae [19, 42]. However, the phylogenetic position of Stratiotes seems to be mainly derived from the mitochondrial sequences (cob, atp1) which are prone to flaws in plant phylogenetic analysis . Therefore, further studies are required to confirm the position of the genus obtained in this study.
Results of divergence time estimates are in agreement with the fossil records of Hydrocharitaceae. The 95% HPD of Najas was 11.9-34.3 Ma, consistent with the oldest fossil of this genus in the Oligocene  (Figure 1). The stem node age of Hydrocharis-Limnobium was dated around 54.7 Ma. However, the crown node age of this subclade was dated around 15.9 Ma, younger than the oldest fossil of Hydrocharis from the Upper Eocene [17, 44, 45]. This could be interpreted as an indication that Limnobium had split from the relatively ancient Hydrocharis in the Miocene (Figure 1), and the great morphological similarity between the two genera is probably due to the short evolutionary history of Limnobium. Although the present study has yielded improved divergence time estimates, it is possible that the estimates of the time of origin for some genera such as Ottelia, Vallisneria, Najas and Blyxa may have been affected by under-representation in sampling.
The age of Hydrocharitaceae estimated in this study (mean: 65.5 Ma, 95% HPD: 54.6-79.6 Ma) is in agreement with that based on rbcL analysis and external fossil calibration points (crown node age = 75 Ma) . However, the stem node age of seagrasses estimated in this study (15.9-41.3 Ma) (Figure 1), is more recent than the 119 ± 11 Ma suggested from analysis using the substitution rates of rbcL and matK . Similarly, our estimates of the stem node age of Ottelia (8.1-33.3 Ma) is more recent than the Cretaceous origin suggested by He et al. (1991) . The split between Zosteraceae and Potamogetonaceae has been dated at 47 Ma by rbcL and fossil calibration . The time is also more recent than the 100 Ma inferred from analyses using the substitution rates of rbcL and matK . These discrepancies indicate that for estimating divergence times in aquatic plants, incorporating fossil calibration point would be more reliable.
The Oriental origin of Hydrocharitaceae inferred from our analysis is supported by the known existence of regions with humid and warm conditions in southeastern Asia during the late Cretaceous and Palaeocene [46, 47] and the fact that the genetic diversity centre of this family is in tropical Asia . The ancestor of clade A was inferred to have originated and diversified in the Orient, while that of clade B dispersed from the Orient to the Southern Hemisphere during the Late Cretaceous and Paleocene (Figure 2c, d, arrow 2). Different environments and oceanic barriers among the major continents (vicariance mechanism) during the Tertiary probably contributed to the diversification of this family resulting in taxa such as the African endemic Lagarosiphon.
Most fossils of Hydrocharitaceae and its close relatives Butomaceae and Alismataceae have been found in Europe (Butomaceae in the Neogene of south Aral region, Miocene of northwest and east Caucasus [18, 48]; Alismataceae in the Tertiary of Europe, a few in North America). The fossil records seem to be inconsistent with the Oriental origin of this family. However, the absence of reports of fossils from Asia most likely reflects a bias in paleobotany, rather than an indication of the origin and past distribution of Hydrocharitaceae. A similar situation exists in Rhinolophus (Rhinolophidae), for which, although the genus is thought to have originated in Asia, fossils have only been reported from Europe and Africa but not from Asia [49, 50].
Does the origin fit with dispersal?
The modern continents viz. South America, Africa, Eurasia, Australia and North America have been separated by oceans since at least ca. 90 Ma [46, 51, 52], earlier than the origin of Hydrocharitaceae. Therefore, dispersal must have played a dominant role in the transoceanic distribution of this family. This contradicts the view that the transoceanic distribution of Ottelia mainly resulted from vicariance . The role of dispersal in transoceanic distribution has been supported by evidences from the studies of geological events and land plant families. Ocean currents are a viable means of dispersal of plants , and a tropical westward-flowing ocean current had spanned the world from the Cretaceous to Paleocene [54, 55]. Island chains existed in the Tethys from Cretaceous to Eocene, which served as a stepping-stone in biotic dispersal between S.E. Asia, Africa and southern Europe [32, 56, 57]. The Malay Archipelago probably facilitated biotic dispersal between S.E. Asia and Australia during the Miocene . The North Atlantic Land Bridge (NALB) aided plant migration between North America and Europe during the late Cretaceous and early Tertiary [59–61]. The Bering Land Bridge (BLB) was open from at least the early Paleocene until its closure ca. 7.4-4.8 Ma . Several recent studies of angiosperms based on molecular and fossil data have supported dispersal as the dominant factor responsible for transoceanic distribution, e.g., in Cucurbitaceae , Sapindaceae , Chrysophylloideae (Sapotaceae) , Burseraceae  and Malphigiaceae . It is probable that Hydrocharitaceae have dispersed to all continents of the world via island chains, land bridges and ocean currents.
Biogeographic studies have suggested that the sub-cosmopolitan distribution of the aquatic plant family Alismataceae has mainly resulted from dispersals (the work will be reported in a separate paper). It is probable that dispersal is the dominant factor, accounting for transoceanic distribution of aquatic angiosperms. However, more studies on aquatic angiosperms are required to investigate this idea further.
Historical biogeography of some genera of Hydrocharitaceae
The ancestor of Stratiotes was suggested to have dispersed from Orient into Europe during the late Cretaceous and Palaeocene (Figure 2b & 2c, 2d, arrow 1), which coincided with the existence of the Tethys seaway (TESW) . Alternatively, the ancestor may have migrated from Orient to Europe across Eurasia. Abundant fossils (15 fossil species) of this genus in Europe  suggested that the genus had diversified widely in this region adapting to wet swamps in the Late Cretaceous .
The genus Hydrilla is native to Eurasia and Australia , and introduced to Americas  and parts of Africa . The centre of differentiation of the genus was thought by Cook and Luond (1982) to lie in tropical Asia . This idea got support from genetic diversity analysis which revealed that the highest diversity is located in China and with lower albeit similar genetic types occurring in Africa, India and USA . Hydrilla might have arisen in the Orient dispersing to Europe and Australia (Figure 2, e, arrow 1 & 5).
The MRCA of the seagrasses within Hydrocharitaceae were suggested to have lived in Oriental area during the Oligocene and Miocene. The result is in agreement with the view that seagrasses possibly originated in the S.E. Asia [72, 73]. The result is supported by the environments of S.E. Asia which was characterized by abundant islands, spacious shallow-seas, warm temperature and plenty of isolated seas . However, the result denied the Cretaceous origin of the group which has been suggested in previous studies [13, 75, 76]. The seagrasses were suggested to have been dispersed from Oriental to other regions (Figure 2c), probably by ocean currents [73, 77]. For example, the warm northward Kuroshio Current carried seagrasses from the equatorial region to the Nansei Islands . Seagrasses are capable of surviving during the LDD between major ocean systems .
Vallisneria has a world-wide distribution, with the highest number of species in Australia [79, 80]. Les et al. (2008)  resolved the phylogeny of this genus, but they conceded that the geographical origin is difficult to pinpoint. In this study by DIVA analysis, Oriental and Australasian areas were suggested as the co-existed ancestral areas of Vallisneria. However, Oriental area is more likely the centre of origin considering the following facts: the closest relative of Vallisneria namely Nechamandra is confined to Asia ; the ancestral species in Vallisneria namely Vallisneria spinulosa, V. spiralis and V. denseserrulata are confined to the Old World [67, 80].
Evolution of morphological characters
Ancestral state reconstruction of reproductive system in Hydrocharitaceae provides empirical evidence that evolution of dioecy in plants has been a bidirectional, viz. from dioecy to hermaphroditism, and from hermaphroditism to dioecy (Figure 3a). This view is supported by Delph (2009)  and Canovas et al. (2011) , but rejects the view that hermaphroditism is the ancestral state in Hydrocharitaceae .
The evolution of leaf habit and leaf shape in Hydrocharitaceae provides several cases of evolutionary adaptation to diverse habitats. The evolution from aerial-submerged leaf to submerged leaf is probably due to change in habitat from shallow to deep waters . The reverse evolution from submerged leaf to aerial-submerged leaf in Ottelia is probably an adaptation to change in habitat from deep to shallow water or some other disadvantageous habitat(s). Taxa with broad-circular leaves (e.g., Ottelia and Hydrocharis) usually occur in still water, while those with ribbon like leaves such as Enhalus and Thalassia occur in coastal waters with strong waves .
In summary, this study has reconstructed the phylogeny of Hydrocharitaceae. The family was suggested to originate in Oriental area during the Late Cretaceous and Paleocene (54.7-72.6 Ma). Dispersal is the most likely factor shaping the transoceanic distribution of this family. Ancestral character state reconstruction of gender and leaf morphology offered valuable information for understanding adaptive evolution in aquatic monocots. However, the historical biogeography for some genera (e.g., Ottelia, Vallisneria) may suffer from under-representation in sampling, and require re-evaluation in future studies.
Sampling and molecular protocols
Most materials used for DNA sequencing was collected from Wuhan Botanical Garden. Some were collected from natural populations in China. Eight genes were used, among which 18S is from nuclear; rbcL, matK, trnK5' intron, rpoB and rpoC1 are from chloroplasts; cob and atp1 are from mitochondria. A detailed list including the voucher information and GenBank accession numbers is provided in Additional file 1.
Genomic DNA was extracted from silica-dried leaves using the Plant Genomic DNA Isolation Kit (Dingguo Biotech, Beijing, China). All polymerase chain reactions (PCR) were conducted in the ABI 2720 Thermal Cycler (Applied Biosystems) in 40 μl volume containing 4 μl of 10 × amplification buffer (200 mM Tris-HCl (pH 8.4), 200 mM KCl, 100 mM (NH4)2SO4, 20 mM MgSO4), 0.8 μl of each primer (10 μM), 0.8 μl of dNTPs (10 mM), 2 U of Taq DNA polymerase (TransGen Biotech Co., Beijing, China) and 60 ng of DNA template. For cob, atp1, 18S, trnK 5' intron, rpoB and rpoC1, the following PCR profile was adopted: 94°C for 3 min, 35 cycles of 30 s at 94°C, 30 s at 50°C, 1 min at 72°C and a final step for 10 min at 70°C. For matK, 55°C Tm and 1 min 30 s extension times were used. Primer sequences were obtained from previous studies: cob (COB1F, COB1R) and atp1 (atpAF1.5, atpAl137r) ; rbcL (1 F, 1204R) ; matK ; trnK 5' intron (3914-F, TRANK2-R) ; 18S (N-NS1, C-18 L) ; rpoB (1f, 4r) and rpoC1 (2f, 4r) . Purified PCR products were double direction sequenced using an automated DNA sequencer (ABI 3730, Applied Biosystems). All newly generated sequences were deposited in GenBank (Additional file 1).
All sequences were aligned individually using Clustal X v2.0 . The output was manually inspected, and ambiguously aligned parts were excluded. Inspired by other studies (e.g. [37, 39, 89, 90]), we assembled all the aligned sequences into a supermatrix data set (combined data set), which was used in phylogenetic analyses. Butomus (Butomaceae) was used as outgroup according to Les et al. (2006) . We also selected (Butomus + Alisma + Cymodocea + Hydrocleys + Potamogeton) as outgroup in order to investigate the influence of outgroup in topology of Hydrocharitaceae.
ML analysis was conducted using RAxML v7.2.5  via the Cyberinfrastructure for Phylogenetic Research (CIPRES) Portal http://www.phylo.org. Two strategies were employed, one involved partitioning the supermatrix data set into eight genes, while another did not partition the data set. GTRCAT and GTRGAMMA options were used, 1000 rapid bootstrap replicates were conducted to assess bootstrap values.
Bayesian analysis was conducted in MrBayes v3.1.2 . The best-fit model of nucleotide substitution was chosen by MrModeltest v2.3  according to the Akaike Information Criterion (combined data set: GTR + I + G model). Two separate runs of four concurrent runs (one cold, three heated each) of 16,000,000 generations were employed with sampling at every 1,000 generations. The stationarity of the likelihood scores of sampled trees was evaluated in Tracer v1.5 , and the first 10% generations were discarded as burn-in.
Divergence time estimates
Divergence time estimates were conducted in BEAST v1.5.4  using the supermatrix data set. Butomus, Alisma, Cymodocea, Hydrocleys and Potamogeton were selected as outgroup. To prevent the negative effects from heterogeneity of substitution rates and uncertainty of fossil data, we used a relaxed clock and Uncorrelated Lognormal (UCLN) model [96, 97]. GTR + I + G model with Gamma Categories set to 6 was adopted. The starting tree was randomly generated with a Yule process prior. More than 90,000,000 generations of MCMC were implemented of which every 1,000 generations were sampled. Tracer v1.5  was used to check the parameters and the first 10% generations were discarded as burn-in.
Lognormal distribution was selected for each calibration point according to Adamson et al. (2010) . This distribution defined the minimum ages for calibrated nodes but allowed the maximum ages to be estimated following a lognormal distribution without hard limit . Three calibration points were incorporated. The oldest reliable fossils of Hydrilla and Vallisneria were reported from the Upper Eocene (33.7-55.8 Ma) [17, 44]. Therefore, the split between Hydilla and (Vallisneria + Nechamandra) was constrained to a minimum of 33.7 Ma (offset = 33.7, mean = 1.1, SD = 1.2). In addition, the oldest reliable fossil of Ottelia was from the Upper Eocene (33.7-55.8 Ma) [17, 44]. Therefore, the split between Blyxa and Ottelia was constrained to a minimum of 33.7 Ma (offset = 33.7, mean = 1.1, SD = 1.2). Lastly, the oldest fossil of this family (genus Stratiotes) was 0.1 Ma younger than the Paleocene-Eocene boundary (54.6 Ma) ). Therefore, the fossil was used to set the split between Stratiotes and the remaining genera of this family not later than 54.5 Ma (offset = 54.5, mean = 1.0, SD = 1.0).
Seven biogeographic areas were recognized according to Dr Morse  (Figure 2a). Biogeographic distribution of Hydrocharitaceae was mainly compiled from literature [36, 67, 72, 78, 80, 100–106]. Fossil was not considered in the area coding as no fossil has been found outside the natural distribution of any genus. Distribution that is known to have been caused by human activities was not included in the analyses. Two methods were used in the analyses: dispersal-vicariance analysis implemented in DIVA v1.2  with the maximum number of ancestral areas at each node constrained to four, and parsimony ancestral state reconstruction implemented in Mesquite .
Two strategies were applied in the biogeographic analyses. One strategy used genera as terminal taxa in the analyses. A tree (Figure 2b) that represented the generic topologies inferred from the phylogenetic analysis using the supermatrix (Figure 1) was constructed, and it was used in the analyses. Each genus was coded based on the current distribution. Details of the distribution are provided in Figure 2b. This strategy followed the methods in the biogeography study of Ranunculaceae  and the suggestion from Yan Yu (one of the authors of S-DIVA). Another strategy used species as terminal taxa and a tree including 72 species of Hydrocharitaceae was used (Figure 2c). Species belonging to Ottelia , Halophila  and Vallisneria  with known phylogenetic relationships were manually added to the tree which resulted from the phylogenetic analysis (Figure 1). In addition, seven species of Najas with the topology from a ML analysis based on rbcL were also added to the tree. Each species was coded based on the current distribution. Details of the distribution are provided in Figure 2c. The purpose of using species as terminal taxa was to reconstruct the ancestral areas at the family and genus levels.
Ancestral character state reconstructions
Information on the reproductive system, leaf habit and leaf shape was mainly compiled from literature [1, 36, 67, 103]. Details of the phenotypic data are provided in Figure 3. Parsimony ancestral state reconstruction was performed using the Mesquite  and the tree inferred from the phylogenetic analyses using the supermatrix dataset.
Additional file 1 Taxa included in this study with voucher information and GenBank accession numbers (DOC 99 kb)
We thank Mao-Xian He, Chun-Feng Yang, Shu-Ying Zhao, Fan Liu, Zhi-Yuan Du, Kuo Liao, Zhong-Ming Ye and Wen-Kui Dai for help in the field work. Ji Yang, Yu-Guo Wang and Yong-Qing Zhu assisted in data analyses. We thank Can Dai and Jun Wen for revising the manuscript. We thank Yang-Zhou Wang and Zhao-Bo Li for comments on this manuscript. We are grateful to Prof Miguel Vences and the two anonymous reviewers for valuable suggestions on this manuscript. This work received financial support from the Foundation of the Chinese Academy of Sciences granted to Qing-Feng Wang (KSCX2-YW-Z-0805) and the National Natural Science Foundation of China (No. 30970202 and 30970195).
- Les DH, Moody ML, Soros CL: A reappraisal of phylogenetic relationships in the monocotyledon family Hydrocharitaceae (Alismatidae). Aliso. 2006, 22: 211-230.Google Scholar
- Les DH, Cleland MA, Waycott M: Phylogenetic studies in Alismatidae, II: Evolution of marine angiosperms (seagrasses) and hydrophily. Syst Bot. 1997, 22: 443-463. 10.2307/2419820.Google Scholar
- He JB: Systematic Botanical and Biosystematic Studies on Ottelia in China. 1991, Wuhan: Wuhan University Press, 1Google Scholar
- Chambers PA, Lacoul P, Murphy KJ, Thomaz SM: Global diversity of aquatic macrophytes in freshwater. Hydrobiologia. 2008, 595: 9-26. 10.1007/s10750-007-9154-6.Google Scholar
- Les DH, Moody ML, Jacobs SWL: Phylogeny and systematics of Aponogeton (Aponogetonaceae): The Australian species. Syst Bot. 2005, 30: 503-519. 10.1600/0363644054782215.Google Scholar
- Les DH, Haynes RR: Systematics of subclass Alismatidae: a synthesis of approaches. Monocotyledons: Systematics and Evolution. Edited by: Rudall PJ, Cribb PJ, Cutler DF, Humphries CJ. 1995, Kew: Royal Botanic Gardens, 353-377.Google Scholar
- Lehtonen S, Myllys L: Cladistic analysis of Echinodorus (Alismataceae): simultaneous analysis of molecular and morphological data. Cladistics. 2008, 24: 218-239. 10.1111/j.1096-0031.2007.00177.x.Google Scholar
- Tanaka N, Setoguchi H, Murata J: Phylogeny of the family Hydrocharitaceae inferred from rbcL and matK gene sequence data. J Plant Res. 1997, 110: 329-337. 10.1007/BF02524931.Google Scholar
- Petersen G, Seberg O, Davis JI, Stevenson DW: RNA editing and phylogenetic reconstruction in two monocot mitochondrial genes. Taxon. 2006, 55: 871-886. 10.2307/25065682.Google Scholar
- Cuenca A, Petersen G, Seberg O, Davis JI, Stevenson DW: Are substitution rates and RNA editing correlated?. BMC Evol Biol. 2010, 10: 349-PubMedPubMed CentralGoogle Scholar
- Janssen T, Bremer K: The age of major monocot groups inferred from 800 + rbcL sequences. Bot J Linn Soc. 2004, 146: 385-398. 10.1111/j.1095-8339.2004.00345.x.Google Scholar
- Ramirez SR, Gravendeel B, Singer RB, Marshall CR, Pierce NE: Dating the origin of the Orchidaceae from a fossil orchid with its pollinator. Nature. 2007, 448: 1042-1045. 10.1038/nature06039.PubMedGoogle Scholar
- Kato Y, Aioi K, Omori Y, Takahata N, Satta Y: Phylogenetic analyses of Zostera species based on rbcL and matK nucleotide sequences: Implications for the origin and diversification of seagrasses in Japanese waters. Genes Genet Syst. 2003, 78: 329-342. 10.1266/ggs.78.329.PubMedGoogle Scholar
- Magallon S, Castillo A: Angiosperm Diversification through Time. Am J Bot. 2009, 96: 349-365. 10.3732/ajb.0800060.PubMedGoogle Scholar
- He JB, Sun XZ, Zhong Y, Huang DS: Cladistic studies on the genus Ottelia (Hydrocharitaceae). J Wuhan Bot Res. 1991, 9: 121-129.Google Scholar
- Kvacek Z: Bilina: a window on Early Miocene marshland environments. Rev Palaeobot Palynol. 1998, 101: 111-123. 10.1016/S0034-6667(97)00072-9.Google Scholar
- Stockey RA: The fossil record of basal monocots. Aliso. 2006, 22: 91-106.Google Scholar
- Daghlian CP: A review of the fossil record of monocotyledons. Bot Rev. 1981, 47: 517-555. 10.1007/BF02860540.Google Scholar
- Sille NP, Collinson ME, Kucera M, Hooker JJ: Morphological evolution of Stratiotes through the Paleogene in England: An example of microevolution in flowering plants. Palaios. 2006, 21: 272-288. 10.2110/palo.2005.P05-21e.Google Scholar
- Cook CDK: Hydrocharitaceae. Flowering Plants, Monocotyledons: Alismatanae and Commelinanae (except Gramineae). Edited by: Kubitzki K, Huber H. 1998, Berlin: Springer, 234-247.Google Scholar
- Les DH, Crawford DJ, Kimball RT, Moody ML, Landolt E: Biogeography of discontinuously distributed hydrophytes: A molecular appraisal of intercontinental disjunctions. Int J Plant Sci. 2003, 164: 917-932. 10.1086/378650.Google Scholar
- Emadzade K, Horandl E: Northern Hemisphere origin, transoceanic dispersal, and diversification of Ranunculeae DC. (Ranunculaceae) in the Cenozoic. J Biogeogr. 2011, 38: 517-530. 10.1111/j.1365-2699.2010.02404.x.Google Scholar
- Emadzade K, Gehrke B, Linder HP, Horandl E: The biogeographical history of the cosmopolitan genus Ranunculus L. (Ranunculaceae) in the temperate to meridional zones. Mol Phylogenet Evol. 2011, 58: 4-21. 10.1016/j.ympev.2010.11.002.PubMedGoogle Scholar
- Raven PH, Axelrod DI: Angiosperm biogeography and past continental movements. Ann Mo Bot Gard. 1974, 61: 539-673. 10.2307/2395021.Google Scholar
- Haynes RR, Holm-Nielsen LB: The Limnocharitaceae. Flora Neotropica. 1992, New York: New York Botanical Garden Press, 56: 1-32.Google Scholar
- Chen JK: Systematic and Evolutionary Botanical Studies on Chinese Sagittaria. 1989, Wuhan: Wuhan University Press, 1Google Scholar
- Barrett SCH: The evolution of plant sexual diversity. Nat Rev Genet. 2002, 3: 274-284. 10.1038/nrg776.PubMedGoogle Scholar
- Sculthorpe CD: The Biology of Aquatic Vascular Plants. 1967, London: Edward ArnoldGoogle Scholar
- Sanders KL, Mumpuni , Hamidy A, Head JJ, Gower DJ: Phylogeny and divergence times of filesnakes (Acrochordus): Inferences from morphology, fossils and three molecular loci. Mol Phylogenet Evol. 2010, 56: 857-867. 10.1016/j.ympev.2010.04.031.PubMedGoogle Scholar
- Sanders KL, Lee MSY: Evaluating molecular clock calibrations using Bayesian analyses with soft and hard bounds. Biol Lett. 2007, 3: 275-279. 10.1098/rsbl.2007.0063.PubMedPubMed CentralGoogle Scholar
- Morse JC: Trichoptera World Checklist. 2011, [http://entweb.clemson.edu/database/trichopt]Google Scholar
- Feng CM, Manchester SR, Xiang QY: Phylogeny and biogeography of Alangiaceae (Cornales) inferred from DNA sequences, morphology, and fossils. Mol Phylogenet Evol. 2009, 51: 201-214. 10.1016/j.ympev.2009.01.017.PubMedGoogle Scholar
- Howarth DG, Gustafsson MHG, Baum DA, Motley TJ: Phylogenetics of the genus Scaevola (Goodeniaceae): implication for dispersal patterns across the Pacific Basin and colonization of the Hawaiian Islands. Am J Bot. 2003, 90: 915-923. 10.3732/ajb.90.6.915.PubMedGoogle Scholar
- Schaefer H, Heibl C, Renner SS: Gourds afloat: a dated phylogeny reveals an Asian origin of the gourd family (Cucurbitaceae) and numerous oversea dispersal events. Proc R Soc B Biol Sci. 2009, 276: 843-851. 10.1098/rspb.2008.1447.Google Scholar
- Cronn R, Wendel JF: Cryptic trysts, genomic mergers, and plant speciation. New Phytol. 2004, 161: 133-142.Google Scholar
- Cook CDK, Urmikonig K: A revision of the genus Ottelia (Hydrocharitaceae). 2. The species of Eurasia, Australasia and America. Aquat Bot. 1984, 20: 131-177. 10.1016/0304-3770(84)90033-0.Google Scholar
- McGowen MR, Spaulding M, Gatesy J: Divergence date estimation and a comprehensive molecular tree of extant cetaceans. Mol Phylogenet Evol. 2009, 53: 891-906. 10.1016/j.ympev.2009.08.018.PubMedGoogle Scholar
- Wiens JJ, Moen DS: Missing data and the accuracy of Bayesian phylogenetics. J Syst Evol. 2008, 46: 307-314.Google Scholar
- Thomson RC, Shaffer HB: Sparse supermatrices for phylogenetic inference: taxonomy, alignment, rogue taxa, and the phylogeny of living turtles. Syst Biol. 2010, 59: 42-58. 10.1093/sysbio/syp075.PubMedGoogle Scholar
- Li XX, Zhou ZK: Phylogenetic studies of the core Alismatales inferred from morphology and rbcL sequences. Prog Nat Sci. 2009, 19: 931-945. 10.1016/j.pnsc.2008.09.008.Google Scholar
- Chen JM, Chen D, Gituru WR, Wang QF, Guo YH: Evolution of apocarpy in Alismatidae using phylogenetic evidence from chloroplast rbcL gene sequence data. Bot Bull Acad Sin. 2004, 45: 33-40.Google Scholar
- Cook CDK, Urmikonig K: A Revision of the Genus Stratiotes (Hydrocharitaceae). Aquat Bot. 1983, 16: 213-249. 10.1016/0304-3770(83)90035-9.Google Scholar
- Hiesel R, von Haeseler A, Brennicke A: Plant mitochondrial nucleic acid sequences as a tool for phylogenetic analysis. Proc Natl Acad Sci USA. 1994, 91: 634-638. 10.1073/pnas.91.2.634.PubMedPubMed CentralGoogle Scholar
- Mai DH, Walther H: Die Floren der Haselbacher Serie im Weisselster-Becken (Bezirk Leipzig, DDR). Abhand-lungen des Staatlichen Museums für Mineralogie und Geologie zu Dresden. 1978, 28: 1-101.Google Scholar
- Mai DH, Walther H: Die obereozänen Floren des Weisselsterbeckens und seiner Randgebiete. Abhandlungen des Staatlichen Museums für Mineralogie und Geologie zu Dresden. 1985, 33: 1-260.Google Scholar
- Markwick PJ: Paul's Palaeo Pages. 2011, [http://www.palaeogeography.net]Google Scholar
- Okada H: Mateer NJ: Cretaceous Environments of Asia. 2000, Netherlands: Elsevier Science B.VGoogle Scholar
- Takhtajan A, Vakrameev V: Radchenko GP: Fundamentals of paleontology: Gymnosperms and angiosperm. 1963, Moscow, Russia: Akadamiia Nauk S.S.S.RGoogle Scholar
- Stoffberg S, Jacobs DS, Mackie IJ, Matthee CA: Molecular phylogenetics and historical biogeography of Rhinolophus bats. Mol Phylogenet Evol. 2010, 54: 1-9. 10.1016/j.ympev.2009.09.021.PubMedGoogle Scholar
- Stoffberg S: Rhinolophus capensis (Chiroptera: Rhinolophidae). Mamm Species. 2008, 810: 1-4.Google Scholar
- Lomolino MV, Riddle BR, Brown JH: Biogeography. 2006, Sunderland: Sinauer Associates, 3Google Scholar
- Davis CC, Bell CD, Mathews S, Donoghue MJ: Laurasian migration explains Gondwanan disjunctions: Evidence from Malpighiaceae. Proc Natl Acad Sci USA. 2002, 99: 6833-6837. 10.1073/pnas.102175899.PubMedPubMed CentralGoogle Scholar
- Renner S: Plant dispersal across the tropical Atlantic by wind and sea currents. Int J Plant Sci. 2004, 165: 23-33. 10.1086/383334.Google Scholar
- Bush ABG: Numerical simulation of the Cretaceous Tethys circumglobal current. Science. 1997, 275: 807-810. 10.1126/science.275.5301.807.PubMedGoogle Scholar
- Gordon WA: Marine life and ocean surface currents in the Cretaceous. J Geol. 1973, 81: 269-284. 10.1086/627870.Google Scholar
- Tiffney BH, Manchester SR: The use of geological and paleontological evidence in evaluating plant phylogeographic hypotheses in the Northern Hemisphere tertiary. Int J Plant Sci. 2001, 162: S3-S17. 10.1086/323880.Google Scholar
- Heaney LR: Conservation biogeography in oceanic archipelagoes. Frontiers of Biogeography: New directions in the Geography of Nature. Edited by: Lomolino MV. 2004, Sunderland: Sinauer Associates, 345-360.Google Scholar
- Van Welzen PC, Slik JWF, Alahuhta J: Plant distribution patterns and plate tectonics in Malesia. Biol Skr. 2005, 55: 199-217.Google Scholar
- Tiffney BH: The Eocene North Atlantic land bridge: its importance in Tertiary and modern phytogeography of the Northern Hemisphere. J Arnold Arbor. 1985, 66: 243-273.Google Scholar
- Denk T, Grimsson F, Zetter R: Episodic migration of oaks to Iceland: Evidence for a North Atlantic "land bridge" in the latest Miocene. Am J Bot. 2010, 97: 276-287. 10.3732/ajb.0900195.PubMedGoogle Scholar
- Mao K, Hao G, Liu J, Adams RP, Milne RI: Diversification and biogeography of Juniperus (Cupressaceae): variable diversification rates and multiple intercontinental dispersals. New Phytol. 2010, 188: 254-272. 10.1111/j.1469-8137.2010.03351.x.PubMedGoogle Scholar
- Marincovich L, Gladenkov AY: Evidence for an early opening of the Bering Strait. Nature. 1999, 397: 149-151. 10.1038/16446.Google Scholar
- Buerki S, Forest F, Alvarez N, Nylander JAA, Arrigo N, Sanmartin I: An evaluation of new parsimony-based versus parametric inference methods in biogeography: a case study using the globally distributed plant family Sapindaceae. J Biogeogr. 2011, 38: 531-550. 10.1111/j.1365-2699.2010.02432.x.Google Scholar
- Bartish IV, Antonelli A, Richardson JE, Swenson U: Vicariance or long-distance dispersal: historical biogeography of the pantropical subfamily Chrysophylloideae (Sapotaceae). J Biogeogr. 2011, 38: 177-190. 10.1111/j.1365-2699.2010.02389.x.Google Scholar
- Weeks A, Daly DC, Simpson BB: The phylogenetic history and biogeography of the frankincense and myrrh family (Burseraceae) based on nuclear and chloroplast sequence data. Mol Phylogenet Evol. 2005, 35: 85-101. 10.1016/j.ympev.2004.12.021.PubMedGoogle Scholar
- Col J: Enchanted Learning. 2011, [http://www.enchantedlearning.com/subjects/dinosaurs/mesozoic/cretaceous/lc.shtml]Google Scholar
- Wang QF, Guo YH, Haynes RR, Hellquist CB: Hydrocharitaceae. Flora of China. Edited by: Wu ZY, Peter HR. 2010, Beijing and St. Louis: Science Press and Missouri Botanical Garden Press, 23: 91-102.Google Scholar
- Haller WT, Sutton DL: Community structure and competition between Hydrilla and Vallisneria. Hyacinth Control J. 1975, 13: 48-50.Google Scholar
- Madeira PT, Coetzee JA, Center TD, White EE, Tipping PW: The origin of Hydrilla verticillata recently discovered at a South African dam. Aquat Bot. 2007, 87: 176-180. 10.1016/j.aquabot.2007.04.008.Google Scholar
- Cook CDK, Luond R: A revision of the genus Hydrilla (Hydrocharitaceae). Aquat Bot. 1982, 13: 485-504.Google Scholar
- Overholt WA, Copeland R, Williams D, Cuda J, Nzigidahera B, Nkubaye E, Wanda F, Gidudu B: Exploration for natural enemies of Hydrilla verticillata in East/Central Africa and genetic characterization of worldwide populations. Report to the Florida Department of Environmental Protection. 2008, Tallahassee: Bureau of Invasive Plant ManagementGoogle Scholar
- Short F, Carruthers T, Dennison W, Waycott M: Global seagrass distribution and diversity: A bioregional model. J Exp Mar Biol Ecol. 2007, 350: 3-20. 10.1016/j.jembe.2007.06.012.Google Scholar
- Mukai H: Biogeography of the tropical seagrasses in the western Pacific. Mar Freshw Res. 1993, 44: 1-17.Google Scholar
- Benzie JAH: Genetic structure of marine organisms and SE Asian biogeography. Biogeography and geological Evolution of SE Asia. Edited by: Hall R, Holloway JD. 1998, Netherlands: Backhuys Publishers, 197-209.Google Scholar
- McCoy ED, Heck KL: Biogeography of corals, seagrasses, and mangroves; an alternative to the center of origin concept. Syst Biol. 1976, 25: 201-210.Google Scholar
- Den HC: The sea-grasses of the world. 1970, Amsterdam: North Holland Publishing CoGoogle Scholar
- Kuo J, Kanamoto Z, Iizumi H, Aioi K, Mukai H: Seagrasses from the Nansei Islands, Southern Japanese Archipelago: species composition, distribution and biogeography. Mar Ecol. 2006, 27: 290-298. 10.1111/j.1439-0485.2006.00098.x.Google Scholar
- Waycott M, Freshwater DW, York RA, Calladine A, Kenworthy WJ: Evolutionary trends in the seagrass genus Halophila (thouars): Insights from molecular phylogeny. Bull Mar Sci. 2002, 71: 1299-1308.Google Scholar
- Jacobs SWL, Frank KA: Notes on Vallisneria (Hydrocharitaceae) in Australia, with descriptions of two new species. Telopea. 1997, 7: 111-118.Google Scholar
- Les DH, Jacobs SWL, Tippery NP, Chen L, Moody ML, Wilstermann-Hildebrand M: Systematics of Vallisneria (Hydrocharitaceae). Syst Bot. 2008, 33: 49-65. 10.1600/036364408783887483.Google Scholar
- Delph LF: Sex allocation: evolution to and from dioecy. Curr Biol. 2009, 19: R249-R251. 10.1016/j.cub.2009.01.048.PubMedGoogle Scholar
- Canovas FG, Mota CF, Serrao EA, Pearson GA: Driving south: a multi-gene phylogeny of the brown algal family Fucaceae reveals relationships and recent drivers of a marine radiation. BMC Evol Biol. 2011, 11: 371-10.1186/1471-2148-11-371.PubMedPubMed CentralGoogle Scholar
- Chen JM, Robert GW, Wang QF: Evolution of aquatic life-forms in Alismatidae: Phylogenetic estimation from chloroplast rbcL gene sequence data. Isr J Plant Sci. 2004, 52: 323-329. 10.1560/N1B6-7T7E-1BL8-V9XC.Google Scholar
- Zurawski G, Clegg M: Evolution of higher-plant chloroplast DNA-encoded genes: Implications for structure-function and phylogenetic studies. Annu Rev Plant Physiol. 1987, 38: 391-418. 10.1146/annurev.pp.38.060187.002135.Google Scholar
- Les DH, Schneider EL, Padgett DJ, Soltis PS, Soltis DE, Zanis M: Phylogeny, classification and floral evolution of water lilies (Nymphaeaceae; Nymphaeales): A synthesis of non-molecular, rbcL, matK, and 18S rDNA data. Syst Bot. 1999, 24: 28-46. 10.2307/2419384.Google Scholar
- Bult C, Kallersjo M, Suh Y: Amplification and sequencing of 16/18S rDNA from gel-purified total plant DNA. Plant Mol Biol Rep. 1992, 10: 273-284. 10.1007/BF02668360.Google Scholar
- Sass C, Little DP, Stevenson DW, Specht CD: DNA barcoding in the cycadales: Testing the potential of proposed barcoding markers for species identification of cycads. PLoS One. 2007, 2: e1154-10.1371/journal.pone.0001154.PubMedPubMed CentralGoogle Scholar
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25: 4876-4882. 10.1093/nar/25.24.4876.PubMedPubMed CentralGoogle Scholar
- Escobar JS, Scornavacca C, Cenci A, Guilhaumon C, Santoni S, Douzery EJP, Ranwez V, Glemin S, David J: Multigenic phylogeny and analysis of tree incongruences in Triticeae (Poaceae). BMC Evol Biol. 2011, 11: 181-10.1186/1471-2148-11-181.PubMedPubMed CentralGoogle Scholar
- Geisler JH, McGowen MR, Yang G, Gatesy J: A supermatrix analysis of genomic, morphological, and paleontological data from crown Cetacea. BMC Evol Biol. 2011, 11: 112-10.1186/1471-2148-11-112.PubMedPubMed CentralGoogle Scholar
- Stamatakis A: RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006, 22: 2688-2690. 10.1093/bioinformatics/btl446.PubMedGoogle Scholar
- Cronquist A: An integrated system of classification of flowering plants. 1981, New York: Columbia University PressGoogle Scholar
- Nylander JAA: MrModeltest v2. Program distributed by the author. 2004, Evolutionary Biology Centre, Uppsala University, SwedenGoogle Scholar
- Rambaut A, Drummond AJ: Tracer v1.4. 2007, Available from http://beast.bio.ed.ac.uk/TracerGoogle Scholar
- Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.PubMedPubMed CentralGoogle Scholar
- Drummond AJ, Ho SY, Phillips MJ, Rambaut A: Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006, 4: e88-10.1371/journal.pbio.0040088.PubMedPubMed CentralGoogle Scholar
- Smith SA, Beaulieu JM, Donoghue MJ: An uncorrelated relaxed-clock analysis suggests an earlier origin for flowering plants. Proc Natl Acad Sci USA. 2010, 107: 5897-5902. 10.1073/pnas.1001225107.PubMedPubMed CentralGoogle Scholar
- Adamson EAS, Hurwood DA, Mather PB: A reappraisal of the evolution of Asian snakehead fishes (Pisces, Channidae) using molecular data from multiple genes and fossil calibration. Mol Phylogenet Evol. 2010, 56: 707-717. 10.1016/j.ympev.2010.03.027.PubMedGoogle Scholar
- Popp M, Mirre V, Brochmann C: A single Mid-Pleistocene long-distance dispersal by a bird can explain the extreme bipolar disjunction in crowberries (Empetrum). Proc Natl Acad Sci USA. 2011, 108: 6520-6525. 10.1073/pnas.1012249108.PubMedPubMed CentralGoogle Scholar
- Eiseman NJ, Mcmillan C: A new species of seagrass, Halophila johnsonii from the Atlantic coast of Florida. Aquat Bot. 1980, 9: 15-19.Google Scholar
- Johnstone IM: Papua-New-Guinea seagrasses and aspects of the biology and growth of Enhalus-acoroides (Lf) Royle. Aquat Bot. 1979, 7: 197-208.Google Scholar
- Kenworthy WJ, Wyllie-Echeverria S, Coles R, Pergent G, Pergent-Martini C: Seagrass conservation biology: An interdisciplinary science for protection of the seagrass biome. Seagrasses: Biology, Ecology and Conservation. Edited by: Larkum AWD, Orth RJ, Duarte CM. 2006, Netherlands: Springer, 595-623.Google Scholar
- Cook CDK: Aquatic plant book. 1990, Hague: SPB Academic PublishingGoogle Scholar
- Wingfield RA, Murphy KJ, Hollingsworth P, Gaywood MJ: The Ecology of Najas Flexilis. 2004, Scottish Natural Heritage Commissioned Report No. 017Google Scholar
- Gambi MC, Barbieri F, Bianchi CN: New record of the alien seagrass Halophila stipulacea (Hydrocharitaceae) in the western Mediterranean: A further clue to changing Mediterranean Sea biogeography. JMBA2 Biodiv Rec. 2009, 2: e84-Google Scholar
- Wang QF, Haynes RR, Hellquist CB: Alismataceae and Butomaceae. Flora of China. Edited by: Wu ZY, Peter HR. 2010, Beijing and St. Louis: Science Press and Missouri Botanical Garden Press, 23: 84-90.Google Scholar
- Ronquist F: Dispersal-vicariance analysis: A new approach to the quantification of historical biogeography. Syst Biol. 1997, 46: 195-203. 10.1093/sysbio/46.1.195.Google Scholar
- Maddison WP, Maddison DR: Mesquite: a modular system for evolutionary analysis. Version 2.74. 2010, [http://mesquiteproject.org]Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.