Plastid Phylogenomic Analyses Resolve Tofieldiaceae as the Root of the Early Diverging Monocot Order Alismatales

The complete chloroplast genome sequence of Lemna turionifera (Araceae)

Lemna turionifera is native to North America and northern Asia, with significant potential for industrial wastewater remediation. The complete nucleotide sequence of the L. turionifera chloroplast genome (cpDNA) was determined. The cpDNA is a circular molecule of 166,606 bp and containing a pair of inverted repeats (IRs) measuting 31,663 bp each. These IRs are flanked by a small single-copy region of 13,542 bp and a large single-copy region of 89,738 bp. The chloroplast genome of L. turionifera consisted of 112 unique genes, including 78 protein-encoding genes, 30 tRNA genes, and four rRNA genes. The phylogenetic analysis utilizing cpDNA provided a well-supported resolution of the relationships among subfamilies within the Araceae family. Our findings indicated a close relationship between L. turionifera and a clade consisting of L. minor, L. japonica, and L. gibba. The availability of the complete chloroplast genome sequence of L. turionifera presents valuable data for future phylogenetic investigations within the Lemnaceae family.

An updated phylogeny and adaptive evolution within Amaranthaceae s.l. inferred from multiple phylogenomic datasets

Amaranthaceae s.l. is a widely distributed family consisting of over 170 genera and 2000 species. Previous molecular phylogenetic studies have shown that Amaranthaceae s.s. and traditional Chenopodiaceae form a monophyletic group (Amaranthaceae s.l.), however, the relationships within this evolutionary branch have yet to be fully resolved. In this study, we assembled the complete plastomes and full-length ITS of 21 Amaranthaceae s.l. individuals and compared them with 38 species of Amaranthaceae s.l. Through plastome structure and sequence alignment analysis, we identified a reverse complementary region approximately 5200 bp long in the genera Atriplex and Chenopodium. Adaptive evolution analysis revealed significant positive selection in eight genes, which likely played a driving role in the evolution of Amaranthaceae s.l., as demonstrated by partitioned evolutionary analysis. Furthermore, we found that about two-thirds of the examined species lack the ycf15 gene, potentially associated with natural selection pressures from their adapted habitats. The phylogenetic tree indicated that some genera (Chenopodium, Halogeton, and Subtr. Salsolinae) are paraphyletic lineages. Our results strongly support the clustering of Amaranthaceae s.l. with monophyletic traditional Chenopodiaceae (Clades I and II) and Amaranthaceae s.s. After a comprehensive analysis, we determined that cytonuclear conflict, gene selection by adapted habitats, and incomplete lineage sorting (ILS) events were the primary reasons for the inconsistent phylogeny of Amaranthaceae s.l. During the last glacial period, certain species within Amaranthaceae s.l. underwent adaptations to different environments and began to differentiate rapidly. Since then, these species may have experienced morphological and genetic changes distinct from those of other genera due to intense selection pressure.

The complete chloroplast genomes of three Alismataceae species, including the medicinally important Alisma orientale

Alismataceae is one of the early diverged families of monocotyledonous plants. We report the complete chloroplast genomes of three Alisma species, including Alisma orientale (Sam.) Juzep. 1934, A. subcordatum Raf. 1908, and A. triviale Pursh 1813, of which A. orientale is a traditional Chinese medical plant used widely to treat diuretics, diabetes, hepatitis, and inflammation. We sequenced the complete chloroplast genomes with the Illumina Nova-Seq 6000 platform using herbarium collections. The chloroplast genomes of A. orientale, A. subcordatum and A. triviale are 159,861 bp, 160,180 bp, and 159,727 bp in length, respectively. The three chloroplast genomes each contain 113 genes, including four rRNAs, 30 tRNAs genes, and 79 protein-coding genes, and the average GC content is 36.0%. Based on the whole chloroplast genomes of 19 species of Alismataceae and the close allies, the medicinally important A. orientale was found to be closely related to another medicinal plant Alisma plantago-aquatica L. 1753 in the phylogenetic analysis. The genus Alisma was supported to be monophyletic.

The complete chloroplast genome of Philodendron hederaceum (Jacq.) Schott 1829 (Alismatales: Araceae)

Philodendron hederaceum (Jacq.) Schott 1829, a species of the Araceae family, is a foliage plant of ornamental value. The complete chloroplast genome sequence of Philodendron hederaceum was obtained by the de novo assembly of NovaSeq 6000 (Illumina Co., San Diego, CA) paired-end short reads and Oxford Nanopore long reads. The complete chloroplast genome of P. hederaceum was 168,609 bp in length, with a large single-copy (LSC) region of 94,393 bp, a small single-copy (SSC) region of 25,774 bp, and a pair of identical inverted repeat regions (IRs) of 24,221 bp. The genome contained a total of 129 genes, including 85 protein-coding genes, 36 transfer RNA (tRNA) genes, and eight ribosomal RNA (rRNA) genes. The phylogenetic analysis of P. hederaceum with 19 related species and two outgroup species revealed the closest taxonomical relationship with Philodendron lanceolatum in the Araceae family.

Comparative Plastid Genome and Phylogenomic Analyses of Potamogeton Species

Potamogetonaceae are aquatic plants divided into six genera. The largest genus in the family is Potamogeton, which is morphologically diverse with many hybrids and polyploids. Potamogetonaceae plastomes were conserved in genome size (155,863 bp-156,669 bp), gene contents (113 genes in total, comprising 79 protein-coding genes and 30 tRNA and 4 rRNA genes), and GC content (36.5%). However, we detected a duplication of the trnH gene in the IR region of the Potamogeton crispus and P. maakianus plastomes. A comparative analysis of Alismatales indicated that the plastomes of Potamogetonaceae, Cymodaceae, and Ruppiaceae have experienced a 6-kb inversion of the rbcL-trnV region and the ndh complex has been lost in the Najas flexilis plastome. Five divergent hotspots (rps16-trnQ, atpF intron, rpoB-trnC, trnC-psbM, and ndhF-rpl32) were identified among the Potamogeton plastomes, which will be useful for species identification. Phylogenetic analyses showed that the family Potamogetonaceae is a well-defined with 100% bootstrap support and divided into two different clades, Potamogeton and Stuckenia. Compared to the nucleotide substitution rates among Alismatales, we found neutral selection in all plastid genes of Potamogeton species. Our results reveal the complete plastome sequences of Potamogeton species, and will be helpful for taxonomic identification, the elucidation of phylogenetic relationships, and the plastome structural analysis of aquatic plants.

Ruppia mongolica (Ruppiaceae), a new species from Inner Mongolia (China), based on morphological and genetic data

Ruppia mongolica Y. Zou & X.W. Xu, a new species from Inner Mongolia, China, is described and illustrated. The phylogenetic position of the new species within the genus was analyzed based on eight chloroplast DNA fragments and an ingroup sampling of all Eurasian species of Ruppia. The results showed that R. mongolica formed a separate branch between R. sinensis and the clade of R. maritima, R. brevipedunculata, R. drepanensis, and R. cirrhosa. Based on molecular and geographical evidence, our study reveals that R. mongolica is closely related to R. sinensis and R. brevipedunculata but differs from the former in the length and shape of the peduncle and seed size, and from the latter in the length of the peduncle, number of carpels per flower, and seed size. In addition, the karyotype analysis revealed that R. mongolica is octoploid, which is first reported within Ruppia, further supporting R. mongolica as a new species.

The complete chloroplast genome sequence of Zannichellia palustris (Potamogetonaceae)

Zannichellia palustris Linnaeus 1753 is a cosmopolitan submerged species capable of rapidly responding to environmental changes, with potential applications in the ecological treatment of heavy metal pollution in water bodies. This study aimed to characterize the complete chloroplast genome of Z. palustris, which has not been reported previously. The chloroplast genome of Z. palustris displays a quadripartite structure with a length of 155,262 base pairs (bp), comprising a large single copy (LSC) region of 85,397 bp, a small single copy (SSC) region of 18,057 bp, and a pair of inverted repeat (IR) regions of 25,904 bp. The GC content of the genome is 35.8%, with corresponding values of 33.4% for the LSC, 28.2% for the SSC, and 42.5% for the IR regions. The genome contained 130 genes, including 85 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Phylogenetic analysis within the order Alismatales revealed that Z. palustris clusters with the clade of Potamogeton perfoliatus, P. crispus and Stuckenia pectinata.

Complete Chloroplast Genomes Provide Insights Into Evolution and Phylogeny of Campylotropis (Fabaceae)

The genus Campylotropis Bunge (Desmodieae, Papilionoideae) comprises about 37 species distributed in temperate and tropical Asia. Despite the great potential in soil conservation, horticulture, and medicine usage, little is known about the evolutionary history and phylogenetic relationships of Campylotropis due to insufficient genetic resources. Here, we sequenced and assembled 21 complete chloroplast genomes of Campylotropis species. In combination with the previously published chloroplast genomes of C. macrocarpa and closely related species, we conducted comparative genomics and phylogenomic analysis on these data. Comparative analysis of the genome size, structure, expansion and contraction of inverted repeat (IR) boundaries, number of genes, GC content, and pattern of simple sequence repeats (SSRs) revealed high similarities among the Campylotropis chloroplast genomes. The activities of long sequence repeats contributed to the variation in genome size and gene content in Campylotropis chloroplast genomes. The Campylotropis chloroplast genomes showed moderate sequence variation, and 13 highly variable regions were identified for species identification and further phylogenetic studies. We also reported one more case of matK pseudogene in the legume family. The phylogenetic analysis confirmed the monophyly of Campylotropis and the sister relationship between Lespedeza and Kummerowia, the latter two genera were then sister to Campylotropis. The intrageneric relationships of Campylotropis based on genomic scale data were firstly reported in this study. The two positively selected genes (atpF and rps19) and eight fast-evolving genes identified in this study may help us to understand the adaptation of Campylotropis species. Overall, this study enhances our understanding of the chloroplast genome evolution and phylogenetic relationships of Campylotropis.

Refined Interpretation of the Pistillate Flower in Ceratophyllum Sheds Fresh Light on Gynoecium Evolution in Angiosperms

Molecular phylogenetic analyses have revealed a superclade of mesangiosperms with five extant lineages: monocots, eudicots, magnoliids, Ceratophyllum and Chloranthaceae. Both Ceratophyllum and Chloranthaceae are ancient lineages with a long fossil record; their precise placement within mesangiosperms is uncertain. Morphological studies have suggested that they form a clade together with some Cretaceous fossils, including Canrightia, Montsechia and Pseudoasterophyllites. Apart from Canrightia, members of this clade share unilocular gynoecia commonly interpreted as monomerous with ascidiate carpels. Alternatively, the gynoecium of Ceratophyllum has also been interpreted as syncarpous with a single fertile carpel (pseudomonomerous). We investigate patterns of morphological, anatomical and developmental variation in gynoecia of three Ceratophyllum species to explore the controversial interpretation of its gynoecium as either monomerous or pseudomonomerous. We use an angiosperm-wide morphological data set and contrasting tree topologies to estimate the ancestral gynoecium type in both Ceratophyllum and mesangiosperms. Gynoecia of all three Ceratophyllum species possess a small (sometimes vestigial) glandular appendage on the abaxial side and an occasionally bifurcating apex. The ovary is usually unilocular with two procambium strands, but sometimes bilocular and/or with three strands in C. demersum. None of the possible phylogenetic placements strongly suggest apocarpy in the stem lineage of Ceratophyllum. Rescoring Ceratophyllum as having two united carpels affects broader-scale reconstructions of the ancestral gynoecium in mesangiosperms. Our interpretation of the glandular appendage as a tepal or staminode homologue makes the Ceratophyllum ovary inferior, thus resembling (semi)inferior ovaries of most Chloranthaceae and potentially related fossils Canrightia and Zlatkocarpus. The entire structure of the flower of Ceratophyllum suggests strong reduction following a long and complex evolutionary history. The widely accepted notion that apocarpy is ancestral in mesangiosperms (and angiosperms) lacks robust support, regardless of which modes of carpel fusion are considered. Our study highlights the crucial importance of incorporating fossils into large-scale analyses to understand character evolution.

Phylogenomic Analyses of Alismatales Shed Light into Adaptations to Aquatic Environments

Land plants first evolved from freshwater algae, and flowering plants returned to water as early as the Cretaceous and multiple times subsequently. Alismatales is the largest clade of aquatic angiosperms including all marine angiosperms, as well as terrestrial plants. We used Alismatales to explore plant adaptations to aquatic environments by analyzing a data set that included 95 samples (89 Alismatales species) covering four genomes and 91 transcriptomes (59 generated in this study). To provide a basis for investigating adaptations, we assessed phylogenetic conflict and whole-genome duplication (WGD) events in Alismatales. We recovered a relationship for the three main clades in Alismatales as (Tofieldiaceae, Araceae) + core Alismatids. We also found phylogenetic conflict among the three main clades that was best explained by incomplete lineage sorting and introgression. Overall, we identified 18 putative WGD events across Alismatales. One of them occurred at the most recent common ancestor of core Alismatids, and three occurred at seagrass lineages. We also found that lineage and life-form were both important for different evolutionary patterns for the genes related to freshwater and marine adaptation. For example, several light- or ethylene-related genes were lost in the seagrass Zosteraceae, but are present in other seagrasses and freshwater species. Stomata-related genes were lost in both submersed freshwater species and seagrasses. Nicotianamine synthase genes, which are important in iron intake, expanded in both submersed freshwater species and seagrasses. Our results advance the understanding of the adaptation to aquatic environments and WGDs using phylogenomics.

Comparative genome sequence and phylogenetic analysis of chloroplast for evolutionary relationship among Pinus species

Genus Pinus is a widely dispersed genus of conifer plants in the Northern Hemisphere. However, the inadequate accessibility of genomic knowledge limits our understanding of molecular phylogeny and evolution of Pinus species. In this study, the evolutionary features of complete plastid genome and the phylogeny of the Pinus genus were studied. A total of thirteen divergent hotspot regions (trnk-UUU, matK, trnQ-UUG, atpF, atpH, rpoC1, rpoC2, rpoB, ycf2, ycf1, trnD-GUC, trnY-GUA, and trnH-GUG) were identified that would be utilized as possible genetic markers for determination of phylogeny and population genetics analysis of Pinus species. Furthermore, seven genes (petD, psaI, psaM, matK, rps18, ycf1, and ycf2) with positive selection site in Pinus species were identified. Based on the whole genome this phylogenetic study showed that twenty-four Pinus species form a significant genealogical clade. Divergence time showed that the Pinus species originated about 100 million years ago (MYA) (95% HPD, 101.76.35–109.79 MYA), in lateral stages of Cretaceous. Moreover, two of the subgenera are consequently originated in 85.05 MYA (95% HPD, 81.04–88.02 MYA). This study provides a phylogenetic relationship and a chronological framework for the future study of the molecular evolution of the Pinus species.

Comparative plastid genomics of Mazaceae: focusing on a new recognized genus, Puchiumazus

Six Mazaceae plastomes were assembled in this study, including the newly recognized genus, Puchiumazus. Comparative plastid genomic analysis provided new insights into Mazaceae. The phylogenetic categorization of Mazus lanceifolius (Mazaceae) has long been uncertain. In 2021, the scholars Bo Li, D. G. Zhang, and C. L. Xiang republished M. lanceifolius as a new species Puchiumazus lanceifolius, within a new genus Puchiumazus. However, there is little plastome information on Mazaceae. Following the publishing of the new genus Puchiumazus, it is now necessary to study the Mazaceae plastome features to comprehensively understand this young family. The Mazaceae plastomes all have a typical quartile structure. The plastomes have a size ranging from 152,388 to 154,252 bp, and each plastome contains 112 unique genes, including 78 protein-coding genes, 4 rRNA genes, and 30 tRNA genes. A comparative analysis showed that these plastome sequences are highly conserved. Furthermore, we identified four relatively hypervariable regions (trnQ-UUC-psbK, trnS-GCU- trnS-CGA, trnT-UGU-trnL-UAA and ycf1) that can be used as potential DNA barcodes for the identification of this clade. Phylogenetic relationships based on the whole plastome sequences of 25 samples of 14 genera of Lamiales placed M. lanceifolius in the basal clade of the family Mazaceae, with 100% bootstrap support. In summary, the M. lanceifolius results indicate that a new monotype genus (Puchiumazus) should be established at the whole-plastome level. This study provides plastid genomic resources for exploring the phylogeny of Mazaceae.

Research Progress in Plant Molecular Systematics of Lauraceae

Simple Summary

Lauraceae, as an angiosperm group with high ecological and economic value, has been widely studied. With the development of science and technology, the research of Lauraceae has changed from morphology to molecular systematics. Our paper reviewed the molecular phylogeny of Lauraceae in recent years. From the aspects of gene fragments, chloroplast genome and DNA barcodes, we mainly discussed the establishment of Cinnamomeae in the ‘Core Lauraceae’ and the tribal controversial genera (Actinodaphne, Sassafras, Caryodaphnopsis, Neocinnamomum and Cassytha). We think that the whole genome and inflorescence characteristics are the breakthrough to solve the tribal problem of Lauraceae. Using reliable molecular and morphological evidence to reconstruct the phylogenetic relationship of Lauraceae will provide an important theoretical basis for the rational utilization of Lauraceae resources, the development of potential resources and the protection of rare plants.

Abstract

Lauraceae is a large family of woody plants with high ecological and economic value. The tribal and generic division and phylogenetic relationship of Lauraceae have long been controversial. Based on morphological and molecular evidence, phylogenetic relationships within the Cinnamomeae, Laureae and Perseeae tribes, also called ‘the Core Lauraceae’, have arisen particular attention. In this review, we comprehensively collated the literatures on the phylogeny of Lauraceae published in recent years and summarized progress made in molecular systematic researches employing gene fragments, chloroplast genomes and DNA barcodings analyses. We clarified the phylogenetic relationships and main controversies of ‘the Core Lauraceae’, the systemic position of fuzzy genera (Neocinnamomum, Caryodaphnopsis and Cassytha) and the development of chloroplast genome and DNA barcodes. We further suggested and proposed the whole genome analysis and different inflorescence types would be possible to provide more information for further research on phylogenetic relationships and taxonomy of Lauraceae.

The complete plastome sequences of five Aponogeton species (Aponogetonaceae): Insights into the structural organization and mutational hotspots

Members of the aquatic plant genus Aponogeton are widely used commercially in aquariums because of their variable leaf shape and unique inflorescences. However, due to extensive similarity between species in this genus, morphological characters are generally inadequate for taxonomic classification. Currently, molecular makers available for taxonomic and phylogenetic studies of Aponogeton are limited. One approach to clarifying relationships between species in these complex groups is to use divergence hotspot regions within the genome. Here, we sequenced and analyzed the plastomes of five Aponogeton species collected from China, Zambia, and Kenya, and subsequently screened these plastomes for divergent DNA hotspots. The five plastomes are circular structures with sizes ranging from 154,167 bp to 154,860 bp. The Large and the Small Single Copies are separated by two Inverted Repeats. One hundred and thirteen unique genes were identified including 79 protein-coding, 30 tRNA, and four rRNA genes. We found that the most abundant repeats in all but one species were mononucleotide repeats (A/T) and that there were 23 potential RNA ending sites. Interestingly, a ~3 kb inversion, which includes the accD gene, was detected within the Asian species of Aponogeton. The inversion may be related to more frequent exchanges between this region and the nuclear genome. Furthermore, we detected mutational hotspot sites among the five Aponogeton species. Three of these hotspots are intergenic spacer regions (accD-psaI, rbcL-accD and trnH-GUG-psbA) that might be suitable for use as barcodes to resolve intra-generic relationships. We also identified four highly variable protein-coding genes (ccsA, rpl22, rps16 and ycf1) may be used as barcodes to resolve the higher-level phylogenies. Our study will provide valuable molecular resources for the taxonomic and phylogenomic study of the complex genus Aponogeton.

Complete Chloroplast Genome Characterization of Oxalis Corniculata and Its Comparison with Related Species from Family Oxalidaceae

Oxalis corniculata L. (family Oxalidaceae) is a small creeper wood sorrel plant that grows well in moist climates. Despite being medicinally important, little is known about the genomics of this species. Here, we determined the complete chloroplast genome sequence of O. corniculata for the first time and compared it with other members of family Oxalidaceae. The genome was 152,189 bp in size and comprised of a pair of 25,387 bp inverted repeats (IR) that separated a large 83,427 bp single copy region (LSC) and a small 16,990 bp single copy region (SSC). The chloroplast genome of O. corniculata contains 131 genes with 83 protein coding genes, 40 tRNA genes, and 8 rRNA genes. The analysis revealed 46 microsatellites, of which 6 were present in coding sequences (CDS) regions, 34 in the LSC, 8 in the SSC, and 2 in the single IR region. Twelve palindromic repeats, 30 forward repeats, and 32 tandem repeats were also detected. Chloroplast genome comparisons revealed an overall high degree of sequence similarity between O. corniculata and O. drummondii and some divergence in the intergenic spacers of related species in Oxalidaceae. Furthermore, the seven most divergent genes (ccsA, clpP, rps8, rps15, rpl22, matK, and ycf1) among genomes were observed. Phylogenomic characterization on the basis of 60 shared genes revealed that O. corniculata is closely related to O. drummondii. The complete O. corniculata genome sequenced in the present study is a valuable resource for investigating the population and evolutionary genetics of family Oxalidaceae and can be used to identify related species.

Complete chloroplast genomes of four Physalis species (Solanaceae): lights into genome structure, comparative analysis, and phylogenetic relationships

BACKGROUND

Physalis L. is a genus of herbaceous plants of the family Solanaceae, which has important medicinal, edible, and ornamental values. The morphological characteristics of Physalis species are similar, and it is difficult to rapidly and accurately distinguish them based only on morphological characteristics. At present, the species classification and phylogeny of Physalis are still controversial. In this study, the complete chloroplast (cp) genomes of four Physalis species (Physalis angulata, P. alkekengi var. franchetii, P. minima and P. pubescens) were sequenced, and the first comprehensive cp genome analysis of Physalis was performed, which included the previously published cp genome sequence of Physalis peruviana.

RESULTS

The Physalis cp genomes exhibited typical quadripartite and circular structures, and were relatively conserved in their structure and gene synteny. However, the Physalis cp genomes showed obvious variations at four regional boundaries, especially those of the inverted repeat and the large single-copy regions. The cp genomes' lengths ranged from 156,578 bp to 157,007 bp. A total of 114 different genes, 80 protein-coding genes, 30 tRNA genes, and 4 rRNA genes, were observed in four new sequenced Physalis cp genomes. Differences in repeat sequences and simple sequence repeats were detected among the Physalis cp genomes. Phylogenetic relationships among 36 species of 11 genera of Solanaceae based on their cp genomes placed Physalis in the middle and upper part of the phylogenetic tree, with a monophyletic evolution having a 100% bootstrap value.

CONCLUSION

Our results enrich the data on the cp genomes of the genus Physalis. The availability of these cp genomes will provide abundant information for further species identification, increase the taxonomic and phylogenetic resolution of Physalis, and assist in the investigation and utilization of Physalis plants.

New Insights Into the Plastome Evolution of the Millettioid/Phaseoloid Clade (Papilionoideae, Leguminosae)

The Millettioid/Phaseoloid (MP) clade from the subfamily Papilionoideae (Leguminosae) consists of six tribes and ca. 3,000 species. Previous studies have revealed some plastome structural variations (PSVs) within this clade. However, many deep evolutionary relationships within the clade remain unresolved. Due to limited taxon sampling and few genetic markers in previous studies, our understanding of the evolutionary history of this clade is limited. To address this issue, we sampled 43 plastomes (35 newly sequenced) representing all the six tribes of the MP clade to examine genomic structural variations and phylogenetic relationships. Plastomes of the species from the MP clade were typically quadripartite (size ranged from 140,029 to 160,040 bp) and contained 109-111 unique genes. We revealed four independent gene losses (ndhF, psbI, rps16, and trnS-GCU), multiple IR-SC boundary shifts, and six inversions in the tribes Desmodieae, Millettieae, and Phaseoleae. Plastomes of the species from the MP clade have experienced significant variations which provide valuable information on the evolution of the clade. Plastid phylogenomic analyses using Maximum Likelihood and Bayesian methods yielded a well-resolved phylogeny at the tribal and generic levels within the MP clade. This result indicates that plastome data is useful and reliable data for resolving the evolutionary relationships of the MP clade. This study provides new insights into the phylogenetic relationships and PSVs within this clade.

The complete chloroplast genomes of two species in threatened monocot genus Caldesia in China

Caldesia is a genus in the family Alismataceae mainly found in the tropical and temperate regions of the Northern hemisphere. In China, two species, Caldesia parnassifolia, and Caldesia grandis are recorded as critically endangered in sporadic regions. Available protection of the genetic resource of these threatened species has been impeded due to limited genomic information. Here, we sequence the whole chloroplast (cp) genome of the two Caldesia species using high throughput sequencing technology. The whole cp genomes of C. parnassifolia and C. grandis were 167,647 bp and 168,500 bp, respectively with a typical quadripartite structure. There were 115 unique genes with 81 protein-coding genes, 31 tRNA genes, and four rRNA genes. Both species showed a GC content of 37.1%. A duplication of two tRNA genes and a ~ 6 kb inversion region in the LSC was noted in both species. Mononucleotide simple sequence repeats (SSRs) A/T were most abundant for both Caldesia species. High nucleotide variability was recorded in ycf1 gene and trnK-UUU/rps16 intergenic spacer region. All RNA editing conversions were C-U in 23 and 24 protein-coding genes for C. parnassifolia and C. grandis, respectively. Phylogenetic analysis placed both Caldesia species as sister to Sagittaria lichuanensis. This study will be useful for further evolutionary, systematic researches and conservation of the genus Caldesia.

Complete Chloroplast Genomes and Comparative Analysis of Sequences Evolution among Seven Aristolochia (Aristolochiaceae) Medicinal Species

Aristolochiaceae, comprising about 600 species, is a unique plant family containing aristolochic acids (AAs). In this study, we sequenced seven species of Aristolochia, and retrieved eleven chloroplast (cp) genomes published for comparative genomics analysis and phylogenetic constructions. The results show that the cp genomes had a typical quadripartite structure with conserved genome arrangement and moderate divergence. The cp genomes range from 159,308 bp to 160,520 bp in length and have a similar GC content of 38.5%–38.9%. A total number of 113 genes were identified, including 79 protein-coding genes, 30 tRNAs and four rRNAs. Although genomic structure and size were highly conserved, the IR-SC boundary regions were variable between these seven cp genomes. The trnH-GUG genes, are one of major differences between the plastomes of the two subgenera Siphisia and Aristolochia. We analyzed the features of nucleotide substitutions, distribution of repeat sequences and simple sequences repeats (SSRs), positive selections in the cp genomes, and identified 16 hotspot regions for genomes divergence that could be utilized as potential markers for phylogeny reconstruction. Phylogenetic relationships of the family Aristolochiaceae inferred from the 18 cp genome sequences were consistent and robust, using maximum parsimony (MP), maximum likelihood (ML), and Bayesian analysis (BI) methods.

Complete chloroplast genome sequences of Dioscorea: Characterization, genomic resources, and phylogenetic analyses

Dioscorea L., the largest genus of the family Dioscoreaceae with over 600 species, is not only an important food but also a medicinal plant. The identification and classification of Dioscorea L. is a rather difficult task. In this study, we sequenced five Dioscorea chloroplast genomes, and analyzed with four other chloroplast genomes of Dioscorea species from GenBank. The Dioscorea chloroplast genomes displayed the typical quadripartite structure of angiosperms, which consisted of a pair of inverted repeats separated by a large single-copy region, and a small single-copy region. The location and distribution of repeat sequences and microsatellites were determined, and the rapidly evolving chloroplast genome regions (trnK-trnQ, trnS-trnG, trnC-petN, trnE-trnT, petG-trnW-trnP, ndhF, trnL-rpl32, and ycf1) were detected. Phylogenetic relationships of Dioscorea inferred from chloroplast genomes obtained high support even in shortest internodes. Thus, chloroplast genome sequences provide potential molecular markers and genomic resources for phylogeny and species identification.

Intergeneric Relationships within the Early-Diverging Angiosperm Family Nymphaeaceae Based on Chloroplast Phylogenomics

The order Nymphaeales, consisting of three families with a record of eight genera, has gained significant interest from botanists, probably due to its position as a basal angiosperm. The phylogenetic relationships within the order have been well studied; however, a few controversial nodes still remain in the Nymphaeaceae. The position of the Nuphar genus and the monophyly of the Nymphaeaceae family remain uncertain. This study adds to the increasing number of the completely sequenced plastid genomes of the Nymphaeales and applies a large chloroplast gene data set in reconstructing the intergeneric relationships within the Nymphaeaceae. Five complete chloroplast genomes were newly generated, including a first for the monotypic Euryale genus. Using a set of 66 protein-coding genes from the chloroplast genomes of 17 taxa, the phylogenetic position of Nuphar was determined and a monophyletic Nymphaeaceae family was obtained with convincing statistical support from both partitioned and unpartitioned data schemes. Although genomic comparative analyses revealed a high degree of synteny among the chloroplast genomes of the ancient angiosperms, key minor variations were evident, particularly in the contraction/expansion of the inverted-repeat regions and in RNA-editing events. Genome structure, and gene content and arrangement were highly conserved among the chloroplast genomes. The intergeneric relationships defined in this study are congruent with those inferred using morphological data.

Evolutionary Comparisons of the Chloroplast Genome in Lauraceae and Insights into Loss Events in the Magnoliids

Available plastomes of the Lauraceae show similar structure and varied size, but there has been no systematic comparison across the family. In order to understand the variation in plastome size and structure in the Lauraceae and related families of magnoliids, we here compare 47 plastomes, 15 newly sequenced, from 27 representative genera. We reveal that the two shortest plastomes are in the parasitic Lauraceae genus Cassytha, with lengths of 114,623 (C. filiformis) and 114,963 bp (C. capillaris), and that they have lost NADH dehydrogenase (ndh) genes in the large single-copy region and one entire copy of the inverted repeat (IR) region. The plastomes of the core Lauraceae group, with lengths from 150,749 bp (Nectandra angustifolia) to 152,739 bp (Actinodaphne trichocarpa), have lost trnI-CAU, rpl23, rpl2, a fragment of ycf2, and their intergenic regions in IRb region, whereas the plastomes of the basal Lauraceae group, with lengths from 157,577 bp (Eusideroxylon zwageri) to 158,530 bp (Beilschmiedia tungfangensis), have lost rpl2 in IRa region. The plastomes of Calycanthus (Calycanthaceae, Laurales) have lost rpl2 in IRb region, but the plastome of Caryodaphnopsis henryi (Lauraceae) remain intact, as do those of the nonLaurales magnoliid genera Piper, Liriodendron, and Magnolia. On the basis of our phylogenetic analysis and structural comparisons, different loss events occurred in different lineages of the Laurales, and fragment loss events in the IR regions have largely driven the contraction of the plastome in the Lauraceae. These results provide new insights into the evolution of the Lauraceae as well as the magnoliids as a whole.

Comparative Plastid Genomes of Primula Species: Sequence Divergence and Phylogenetic Relationships

Compared to traditional DNA markers, genome-scale datasets can provide mass information to effectively address historically difficult phylogenies. Primula is the largest genus in the family Primulaceae, with members distributed mainly throughout temperate and arctic areas of the Northern Hemisphere. The phylogenetic relationships among Primula taxa still maintain unresolved, mainly due to intra- and interspecific morphological variation, which was caused by frequent hybridization and introgression. In this study, we sequenced and assembled four complete plastid genomes (Primula handeliana, Primula woodwardii, Primula knuthiana, and Androsace laxa) by Illumina paired-end sequencing. A total of 10 Primula species (including 7 published plastid genomes) were analyzed to investigate the plastid genome sequence divergence and their inferences for the phylogeny of Primula. The 10 Primula plastid genomes were similar in terms of their gene content and order, GC content, and codon usage, but slightly different in the number of the repeat. Moderate sequence divergence was observed among Primula plastid genomes. Phylogenetic analysis strongly supported that Primula was monophyletic and more closely related to Androsace in the Primulaceae family. The phylogenetic relationships among the 10 Primula species showed that the placement of P. knuthiana–P. veris clade was uncertain in the phylogenetic tree. This study indicated that plastid genome data were highly effective to investigate the phylogeny.

The complete chloroplast genome sequence of Epipremnum aureum and its comparative analysis among eight Araceae species

Epipremnum aureum is an important foliage plant in the Araceae family. In this study, we have sequenced the complete chloroplast genome of E. aureum by using Illumina Hiseq sequencing platforms. This genome is a double-stranded circular DNA sequence of 164,831 bp that contains 35.8% GC. The two inverted repeats (IRa and IRb; 26,606 bp) are spaced by a small single-copy region (22,868 bp) and a large single-copy region (88,751 bp). The chloroplast genome has 131 (113 unique) functional genes, including 86 (79 unique) protein-coding genes, 37 (30 unique) tRNA genes, and eight (four unique) rRNA genes. Tandem repeats comprise the majority of the 43 long repetitive sequences. In addition, 111 simple sequence repeats are present, with mononucleotides being the most common type and di- and tetranucleotides being infrequent events. Positive selection pressure on rps12 in the E. aureum chloroplast has been demonstrated via synonymous and nonsynonymous substitution rates and selection pressure sites analyses. Ycf15 and infA are pseudogenes in this species. We constructed a Maximum Likelihood phylogenetic tree based on the complete chloroplast genomes of 38 species from 13 families. Those results strongly indicated that E. aureum is positioned as the sister of Colocasia esculenta within the Araceae family. This work may provide information for further study of the molecular phylogenetic relationships within Araceae, as well as molecular markers and breeding novel varieties by chloroplast genetic-transformation of E. aureum in particular.

Phylogenic study of Lemnoideae (duckweeds) through complete chloroplast genomes for eight accessions

Background

Phylogenetic relationship within different genera of Lemnoideae, a kind of small aquatic monocotyledonous plants, was not well resolved, using either morphological characters or traditional markers. Given that rich genetic information in chloroplast genome makes them particularly useful for phylogenetic studies, we used chloroplast genomes to clarify the phylogeny within Lemnoideae.

Methods

DNAs were sequenced with next-generation sequencing. The duckweeds chloroplast genomes were indirectly filtered from the total DNA data, or directly obtained from chloroplast DNA data. To test the reliability of assembling the chloroplast genome based on the filtration of the total DNA, two methods were used to assemble the chloroplast genome of Landoltia punctata strain ZH0202. A phylogenetic tree was built on the basis of the whole chloroplast genome sequences using MrBayes v.3.2.6 and PhyML 3.0.

Results

Eight complete duckweeds chloroplast genomes were assembled, with lengths ranging from 165,775 bp to 171,152 bp, and each contains 80 protein-coding sequences, four rRNAs, 30 tRNAs and two pseudogenes. The identity of L. punctata strain ZH0202 chloroplast genomes assembled through two methods was 100%, and their sequences and lengths were completely identical. The chloroplast genome comparison demonstrated that the differences in chloroplast genome sizes among the Lemnoideae primarily resulted from variation in non-coding regions, especially from repeat sequence variation. The phylogenetic analysis demonstrated that the different genera of Lemnoideae are derived from each other in the following order: Spirodela, Landoltia, Lemna, Wolffiella, and Wolffia.

Discussion

This study demonstrates potential of whole chloroplast genome DNA as an effective option for phylogenetic studies of Lemnoideae. It also showed the possibility of using chloroplast DNA data to elucidate those phylogenies which were not yet solved well by traditional methods even in plants other than duckweeds.