Plastome phylogenomics and new fossil evidence from Dominican amber shed light on the evolutionary history of the Neotropical fern genus Pecluma.

PREMISE: Molecular studies based on chloroplast markers have questioned the monophyly of the fern genus Pecluma (Polypodioideae, Polypodiaceae), which has several species of Polypodium nested within it. We explored the delimitation of Pecluma and its biogeographic pattern by evaluating the phylogenetic position of four Polypodium species not sequenced thus far and integrating the first fossil evidence of Pecluma. METHODS: Using herbarium material, we applied a genome-skimming approach to obtain a phylogenetic hypothesis of Polypodioideae; assessed the combination of character states observed in the fossil from Miocene Dominican amber using a previously published phylogeny of Polypodioideae based on four plastid markers as framework; calculated divergence times; and conducted an ancestral area estimation. RESULTS: Within Polypodioideae, Pecluma was recovered as sister to Phlebodium. Three of the newly sequenced species-Polypodium otites, P. pinnatissimum, and P. ursipes-were recovered with maximum support within the Pecluma clade, whereas P. christensenii remained within Polypodium. The closest combination of character states of the fossil was found within Pecluma. Our biogeographic analyses suggest an Eocene origin of the genus in South America, with several subsequent Oligocene and Miocene colonization events to Mexico-Central America and to the West Indies. CONCLUSIONS: Although the circumscription of Pecluma is still challenging, our results elucidate the origin and age of the genus. The newly described fossil, Pecluma hispaniolae sp. nov., supports the hypothesis that the epiphytic communities of the Greater Antilles exhibit a constant generic composition since the Miocene. We propose new combinations (Pecluma otites, Pecluma pinnatissima, and Pecluma ursipes) to accommodate three species previously classified in Polypodium.

Plastome evolution of Engelhardia facilitates phylogeny of Juglandaceae.

BACKGROUND: Engelhardia (Juglandaceae) is a genus of significant ecological and economic importance, prevalent in the tropics and subtropics of East Asia. Although previous efforts based on multiple molecular markers providing profound insights into species delimitation and phylogeography of Engelhardia, the maternal genome evolution and phylogeny of Engelhardia in Juglandaceae still need to be comprehensively evaluated. In this study, we sequenced plastomes from 14 samples of eight Engelhardia species and the outgroup Rhoiptelea chiliantha, and incorporated published data from 36 Juglandaceae and six outgroup species to test phylogenetic resolution. Moreover, comparative analyses of the plastomes were conducted to investigate the plastomes evolution of Engelhardia and the whole Juglandaceae family. RESULTS: The 13 Engelhardia plastomes were highly similar in genome size, gene content, and order. They exhibited a typical quadripartite structure, with lengths from 161,069 bp to 162,336 bp. Three mutation hotspot regions (TrnK-rps16, ndhF-rpl32, and ycf1) could be used as effective molecular markers for further phylogenetic analyses and species identification. Insertion and deletion (InDels) may be an important driving factor for the evolution of plastomes in Juglandoideae and Engelhardioideae. A total of ten codons were identified as the optimal codons in Juglandaceae. The mutation pressure mostly contributed to shaping codon usage. Seventy-eight protein-coding genes in Juglandaceae experienced relaxed purifying selection, only rpl22 and psaI genes showed positive selection (Ka/Ks > 1). Phylogenetic results fully supported Engelhardia as a monophyletic group including two sects and the division of Juglandaceae into three subfamilies. The Engelhardia originated in the Late Cretaceous and diversified in the Late Eocene, and Juglandaceae originated in the Early Cretaceous and differentiated in Middle Cretaceous. The phylogeny and divergence times didn't support rapid radiation occurred in the evolution history of Engelhardia. CONCLUSION: Our study fully supported the taxonomic treatment of at the section for Engelhardia species and three subfamilies for Juglandaceae and confirmed the power of phylogenetic resolution using plastome sequences. Moreover, our results also laid the foundation for further studying the course, tempo and mode of plastome evolution of Engelhardia and the whole Juglandaceae family.

Phylogenomic data resolved the deep relationships of Gymnogynoideae (Selaginellaceae).

The unresolved phylogenetic framework within the Selaginellaceae subfamily Gymnogynoideae (ca. 130 species) has hindered our comprehension of the diversification and evolution of Selaginellaceae, one of the most important lineages in land plant evolution. Here, based on plastid and nuclear data extracted from genomic sequencing of more than 90% species of all genera except two in Gymnogynoideae, a phylogenomic study focusing on the contentious relationships among the genera in Gymnogynoideae was conducted. Our major results included the following: (1) Only single-copy region (named NR) and only one ribosomal operon was firstly found in Afroselaginella among vascular plants, the plastome structure of Gymnogynoideae is diverse among the six genera, and the direct repeats (DR) type is inferred as the ancestral state in the subfamily; (2) The first strong evidence was found to support Afroselaginella as a sister to Megaloselaginella. Alternative placements of Ericetorum and Gymnogynum were detected, and their relationships were investigated by analyzing the variation of phylogenetic signals; and (3) The most likely genus-level relationships in Gymnogynoideae might be: ((Bryodesma, Lepidoselaginella), (((Megaloselaginella, Afroselaginella), Ericetorum), Gymnogynum)), which was supported by maximum likelihood phylogeny based on plastid datasets, maximum likelihood, and Bayesian inference based on SCG dataset and concatenated nuclear and plastid datasets and the highest proportion of phylogenetic signals of plastid genes.

Comparative chloroplast genomes: insights into the identification and phylogeny of rapid radiation genus Rhodiola.

Rhodiola L. is a genus exhibiting rapid radiation and represents a typical case for studying plastid gene adaptation in species that spread from high altitudes to low altitudes. In this study, 23 samples of 18 Rhodiola species were collected from the Qinghai-Tibetan Plateau and five scattered alpine areas, and the plastid genomes (plastomes) of these species were sequenced, annotated, and compared between high-altitude and widely distributed groups. The plastomes of Rhodiola were found to be highly conserved in terms of gene size, content, and order but highly variable in several lineage-specific features, such as codon usage bias, IR boundary shifting, and distinct repeat sequence structures binding to SSRs. Codon usage in the genes of photosystem II exhibited an obvious preference, reflecting significant environmental adaptation pressures. In this study, three repeat regions compounded with trinucleotide and mononucleotide repeats were found for the first time in R. forrestii, R. himalensis, and R. yunnanensis. High-variability regions such as ndhF, ycf1, trnH-psbA, and rpoC1-rpoB were screened, laying the foundation for the precise identification of these species. The phylogenetic analysis revealed the occurrence of cyto-nuclear discordance, likely originating from the frequent interspecific hybridization events observed within Rhodiola species during rapid radiation. Dioecious and hermaphrodite species can be broadly categorized into two subclades, probably they have different environmental adaptation strategies in response to climate change. In addition, the phylogenetic tree supported the monophyly of R. forrestii and R. yunnanensis, which compose R. Sect. Pseudorhodiola. In conclusion, plastome data enrich the genetic information available for the Rhodiola genus and may provide insight into species migration events during climate change.

Comparative analysis of plastid genomes reveals rearrangements, repetitive sequence features, and phylogeny in the Annonaceae.

The Annonaceae stands as the most species rich family in the Magnoliales, a basal group of angiosperms. Widely distributed in tropical and subtropical regions, it holds significant ecological and economic value. The plastid genome (plastome) is often employed in studies related to plant phylogenetics, comparative genomics, evolutionary biology, and genetic engineering. Nonetheless, research progress on plastid genomics in the Annonaceae has been relatively slow. In this study, we analyzed the structure and repetitive sequence features of plastomes from 28 Annonaceae species. Among them, Mitrephora tomentosa and Desmos chinensis were newly sequenced, with sizes of 160,157 bp and 192,167 bp, and GC contents of 38.3% and 38.4%, respectively. The plastome size in the Annonaceae ranged from 158,837 bp to 202,703 bp, with inverted repeat (IR) region sizes ranging from 64,621 bp to 25,861 bp. Species exhibiting expansion in the IR region showed an increase in plastome size and gene number, frequent boundary changes, different expansion modes (bidirectional or unidirectional), and an increase in repetitive sequences. Specifically, a large number of dispersed repetitive sequences lead to an increase in the size of the LSC region in Goniothalamus tamirensis. Phylogenetic analysis revealed Annonoideae and Malmeoideae as monophyletic groups and sister clades, with Cananga odorata outside of them, followed by Anaxagorea javanica. This research uncovers the structural variation characteristics of plastomes in the Annonaceae, providing valuable information for understanding the phylogeny and plastome evolution of Annonaceae.

Comparative genomics on chloroplasts of Rubus (Rosaceae).

Rubus, the largest genus in Rosaceae, contains over 1400 species that distributed in multiple habitats across the world, with high species diversity in the temperate regions of Northern Hemisphere. Multiple Rubus species are cultivated for their valuable fruits. However, the intrageneric classification and phylogenetic relationships are still poorly understood. In this study, we sequenced, assembled, and characterized 17 plastomes of Rubus, and conducted comparative genomics integrating with 47 previously issued plastomes of this genus. The 64 plastomes of Rubus exhibited typical quadripartite structure with sizes ranging from 155,144 to 156,700 bp, and contained 132 genes including 87 protein-coding genes, 37 tRNA genes and eight rRNA genes. All plastomes are conservative in the gene order, the frequency of different types of long repeats and simple sequence repeats (SSRs), the codon usage, and the selection pressure of protein-coding genes. However, there are also some differences in the Rubus plastomes, including slight contraction and expansion of the IRs, a variation in the numbers of SSRs and long repeats, and some genes in certain clades undergoing intensified or relaxed purifying selection. Phylogenetic analysis based on whole plastomes showed that the monophyly of Rubus was strongly supported and resolved it into six clades corresponding to six subgenera. Moreover, we identified 12 highly variable regions that could be potential molecular markers for phylogenetic, population genetic, and barcoding studies. Overall, our study provided insight into plastomic structure and sequence diversification of Rubus, which could be beneficial for future studies on identification, evolution, and phylogeny in this genus.

Hybridization, polyploidization, and morphological convergence make dozens of taxa into one chaotic genetic pool: a phylogenomic case of the Ficus erecta species complex (Moraceae).

The Ficus erecta complex, characterized by its morphological diversity and frequent interspecific overlap, shares pollinating fig wasps among several species. This attribute, coupled with its intricate phylogenetic relationships, establishes it as an exemplary model for studying speciation and evolutionary patterns. Extensive researches involving RADseq (Restriction-site associated DNA sequencing), complete chloroplast genome data, and flow cytometry methods were conducted, focusing on phylogenomic analysis, genetic structure, and ploidy detection within the complex. Significantly, the findings exposed a pronounced nuclear-cytoplasmic conflict. This evidence, together with genetic structure analysis, confirmed that hybridization within the complex is a frequent occurrence. The ploidy detection revealed widespread polyploidy, with certain species exhibiting multiple ploidy levels, including 2×, 3×, and 4×. Of particular note, only five species (F. abelii, F. erecta, F. formosana, F. tannoensis and F. vaccinioides) in the complex were proved to be monophyletic. Species such as F. gasparriniana, F. pandurata, and F. stenophylla were found to encompass multiple phylogenetically distinct lineages. This discovery, along with morphological comparisons, suggests a significant underestimation of species diversity within the complex. This study also identified F. tannoensis as an allopolyploid species originating from F. vaccinioide and F. erecta. Considering the integration of morphological, molecular systematics, and cytological evidences, it is proposed that the scope of the F. erecta complex should be expanded to the entire subsect. Frutescentiae. This would redefine the complex as a continuously evolving group comprising at least 33 taxa, characterized by blurred species boundaries, frequent hybridization and polyploidization, and ambiguous genetic differentiation.

Comparative plastomes of eight subgenus Chamaesyce plants and system authentication of Euphorbiae Humifusae Herba.

Euphorbiae Humifusae Herba (EHH) was provided with medicinal and edible uses, but frequently was adulterated with its closely related species. Hence, this study sought to identify EHH via an integrated approach comprising data from its morphological evaluation, HPLC analysis, comparative plastomes analysis and allele-specific PCR identification. First, the morphological characteristics of 8 subgenus Chamaesyce plants were summarized. Then, HPLC analysis showed that 18 batches of EHH were adulterated or unqualified. Furthermore, the plastomes of the 8 subg. Chamaesyce species were analyzed. Phylogenetic analysis revealed a sister relationship among the 8 subg. Chamaesyce species. The allele-specific PCR authentication was developed by the nucleotide polymorphisms (SNPs) and insertions or deletions (InDels) analysis. The results of allele-specific PCR showed that 27 batches of EHH were adulterated, indicating that the superior sensitivity of molecular authentication over the other methods used. This study provided a reference for rational use and phylogenetic research of EHH.

Extreme Reconfiguration of Plastid Genomes in Papaveraceae: Rearrangements, Gene Loss, Pseudogenization, IR Expansion, and Repeats.

The plastid genomes (plastomes) of angiosperms are typically highly conserved, with extreme reconfiguration being uncommon, although reports of such events have emerged in some lineages. In this study, we conducted a comprehensive comparison of the complete plastomes from twenty-two species, covering seventeen genera from three subfamilies (Fumarioideae, Hypecooideae, and Papaveroideae) of Papaveraceae. Our results revealed a high level of variability in the plastid genome size of Papaveraceae, ranging from 151,864 bp to 219,144 bp in length, which might be triggered by the expansion of the IR region and a large number of repeat sequences. Moreover, we detected numerous large-scale rearrangements, primarily occurring in the plastomes of Fumarioideae and Hypecooideae. Frequent gene loss or pseudogenization were also observed for ndhs, accD, clpP, infA, rpl2, rpl20, rpl32, rps16, and several tRNA genes, particularly in Fumarioideae and Hypecooideae, which might be associated with the structural variation in their plastomes. Furthermore, we found that the plastomes of Fumarioideae exhibited a higher GC content and more repeat sequences than those of Papaveroideae. Our results showed that Papaveroideae generally displayed a relatively conserved plastome, with the exception of Eomecon chionantha, while Fumarioideae and Hypecooideae typically harbored highly reconfigurable plastomes, showing high variability in the genome size, gene content, and gene order. This study provides insights into the plastome evolution of Papaveraceae and may contribute to the development of effective molecular markers.

On Pattern-Cladistic Analyses Based on Complete Plastid Genome Sequences.

The fundamental Hennigian principle, grouping solely on synapomorphy, is seldom used in modern phylogenetics. In the submitted paper, we apply this principle in reanalyzing five datasets comprising 197 complete plastid genomes (plastomes). We focused on the latter because plastome-based DNA sequence data gained dramatic popularity in molecular systematics during the last decade. We show that pattern-cladistic analyses based on complete plastid genome sequences can successfully resolve affinities between plant taxa, simultaneously simplifying both the genomic and analytical frameworks of phylogenetic studies. We developed "Matrix to Newick" (M2N), a program to represent the standard molecular alignment of plastid genomes in the form of trees or relationships directly. Thus, massive plastome-based DNA sequence data can be successfully represented in a relational form rather than as a standard molecular alignment. Application of methods of median supertree construction (the Average Consensus method has been used as an example in this study) or Maximum Parsimony analysis to relational representations of plastome sequence data may help systematist to avoid the complicated assumption-based frameworks of Maximum Likelihood or Bayesian phylogenetics that are most used today in massive plastid sequence data analyses. We also found that significant amounts of pure genomic information that typically accommodate the majority of current plastid phylogenomic studies can be effectively dropped by systematists if they focus on the pattern-cladistics or relational analyses of plastome-based molecular data. The proposed pattern-cladistic approach is a powerful and straightforward heuristic alternative to modern plastome-based phylogenetics.

Comparative Plastomes of Curcuma alismatifolia (Zingiberaceae) Reveal Diversified Patterns among 56 Different Cut-Flower Cultivars.

Curcuma alismatifolia (Zingiberaceae) is an ornamental species with high economic value due to its recent rise in popularity among floriculturists. Cultivars within this species have mixed genetic backgrounds from multiple hybridization events and can be difficult to distinguish via morphological and histological methods alone. Given the need to improve identification resources, we carried out the first systematic study using plastomic data wherein genomic evolution and phylogenetic relationships from 56 accessions of C. alismatifolia were analyzed. The newly assembled plastomes were highly conserved and ranged from 162,139 bp to 164,111 bp, including 79 genes that code for proteins, 30 tRNA genes, and 4 rRNA genes. The A/T motif was the most common of SSRs in the assembled genomes. The Ka/Ks values of most genes were less than 1, and only two genes had Ka/Ks values above 1, which were rps15 (1.15), and ndhl (1.13) with petA equal to 1. The sequence divergence between different varieties of C. alismatifolia was large, and the percentage of variation in coding regions was lower than that in the non-coding regions. Such data will improve cultivar identification, marker assisted breeding, and preservation of germplasm resources.

More than a spiny morphology: plastome variation in the prickly pear cacti (Opuntieae).

BACKGROUND: Plastid genomes (plastomes) have long been recognized as highly conserved in their overall structure, size, gene arrangement and content among land plants. However, recent studies have shown that some lineages present unusual variations in some of these features. Members of the cactus family are one of these lineages, with distinct plastome structures reported across disparate lineages, including gene losses, inversions, boundary movements or loss of the canonical inverted repeat (IR) region. However, only a small fraction of cactus diversity has been analysed so far. METHODS: Here, we investigated plastome features of the tribe Opuntieae, the remarkable prickly pear cacti, which represent one of the most diverse and important lineages of Cactaceae. We assembled de novo the plastome of 43 species, representing a comprehensive sampling of the tribe, including all seven genera, and analysed their evolution in a phylogenetic comparative framework. Phylogenomic analyses with different datasets (full plastome sequences and genes only) were performed, followed by congruence analyses to assess signals underlying contentious nodes. KEY RESULTS: Plastomes varied considerably in length, from 121 to 162 kbp, with striking differences in the content and size of the IR region (contraction and expansion events), including a lack of the canonical IR in some lineages and the pseudogenization or loss of some genes. Overall, nine different types of plastomes were reported, deviating in the presence of the IR region or the genes contained in the IR. Overall, plastome sequences resolved phylogenetic relationships within major clades of Opuntieae with high bootstrap values but presented some contentious nodes depending on the dataset analysed (e.g. whole plastome vs. genes only). Congruence analyses revealed that most plastidial regions lack phylogenetic resolution, while few markers are supporting the most likely topology. Likewise, alternative topologies are driven by a handful of plastome markers, suggesting recalcitrant nodes in the phylogeny. CONCLUSIONS: Our study reveals a dynamic nature of plastome evolution across closely related lineages, shedding light on peculiar features of plastomes. Variation of plastome types across Opuntieae is remarkable in size, structure and content and can be important for the recognition of species in some major clades. Unravelling connections between the causes of plastome variation and the consequences for species biology, physiology, ecology, diversification and adaptation is a promising and ambitious endeavour in cactus research. Although plastome data resolved major phylogenetic relationships, the generation of nuclear genomic data is necessary to confront these hypotheses and assess the recalcitrant nodes further.

Comparative plastomes of Pueraria montana var. lobata (Leguminosae: Phaseoleae) and closely related taxa: insights into phylogenomic implications and evolutionary divergence.

BACKGROUND: Pueraria montana var. lobata (kudzu) is an important food and medicinal crop in Asia. However, the phylogenetic relationships between Pueraria montana var. lobata and the other two varieties (P. montana var. thomsonii and P. montana var. montana) remain debated. Although there is increasing evidence showing that P. montana var. lobata adapts to various environments and is an invasive species in America, few studies have systematically investigated the role of the phylogenetic relationships and evolutionary patterns of plastomes between P. montana var. lobata and its closely related taxa. RESULTS: 26 newly sequenced chloroplast genomes of Pueraria accessions resulted in assembled plastomes with sizes ranging from 153,360 bp to 153,551 bp. Each chloroplast genome contained 130 genes, including eight rRNA genes, 37 tRNA genes, and 85 protein-coding genes. For 24 newly sequenced accessions of these three varieties of P. montana, we detected three genes and ten noncoding regions with higher nucleotide diversity (π). After incorporated publically available chloroplast genomes of Pueraria and other legumes, 47 chloroplast genomes were used to construct phylogenetic trees, including seven P. montana var. lobata, 14 P. montana var. thomsonii and six P. montana var. montana. Phylogenetic analysis revealed that P. montana var. lobata and P. montana var. thomsonii formed a clade, while all sampled P. montana var. montana formed another cluster based on cp genomes, LSC, SSC and protein-coding genes. Twenty-six amino acid residues were identified under positive selection with the site model. We also detected six genes (accD, ndhB, ndhC, rpl2, rpoC2, and rps2) that account for among-site variation in selective constraint under the clade model between accessions of the Pueraria montana var. lobata clade and the Pueraria montana var. montana clade. CONCLUSION: Our data provide novel comparative plastid genomic insights into conservative gene content and structure of cp genomes pertaining to P. montana var. lobata and the other two varieties, and reveal an important phylogenetic clue and plastid divergence among related taxa of P. montana come from loci that own moderate variation and underwent modest selection.

Backbone phylogeny and adaptive evolution of Pleurospermum s. l.: New insights from phylogenomic analyses of complete plastome data.

Pleurospermum is a taxonomically challenging taxon of Apiaceae, as its circumscription and composition remain controversial for morphological similarities with several related genera, leading to a dispute between Pleurospermum in the broad sense and strict sense. While evidence from previous molecular studies recognized plural branching lineages within the Pleurospermum s. l., it did not support the latest delimitation of Pleurospermum s. str. by only two closely related northern species. So far, no proper delimitation for Pleurospermum has come up, and many of the plural taxa in Pleurospermum s. l. remain unresolved, which may be due to poor phylogenetic resolution yielded barely from ITS sequences. Herein, we newly assembled 40 complete plastomes from 36 species of Pleurospermum s. l. and related genera, 34 of which were first reported and generated a well-resolved backbone phylogeny in a framework of the subfamily Apioideae. From the phylogeny with greatly improved resolution, a total of six well-supported monophyletic lineages within Pleurospermum s. l. were recognized falling in different major clades of Apioideae. Combining morphological characteristics with phylogenetic inference, we suggested to re-delimit the Pleurospermum s. str. by introducing nine species mainly from the Himalayan regions and proposed its boundary features; the remaining species were suggested to be excluded from Pleurospermum to incorporate into their more related taxa being revealed. On this basis, the plastome comparison revealed not only the high conservatism but also the mild differences among lineages in plastome structure and gene evolution. Overall, our study provided a backbone phylogeny essential for further studies of the taxonomically difficult taxa within Pleurospermum s. l.

Plastome evolution in the East Asian lobelias (Lobelioideae) using phylogenomic and comparative analyses.

Lobelia species, as rich source of the alkaloid lobeline which has been shown to have important biological activity, have been used in folk medicine throughout East Asia to treat various diseases. However, Lobelia is a complex and varied genus in East Asia and is thus difficult to identify. Genomic resources would aid identification, however the availability of such information is poor, preventing a clear understanding of their evolutionary history from being established. To close this gap in the available genomic data, in this study, 17 plastomes of East Asian lobelias were newly sequenced and assembled. Although the plastomes of Lobelia sect. Hypsela, L. sect. Speirema, and L. sect. Rhynchopetalum shared the gene structure, the inverted repeat (IR)/large single copy (LSC) boundaries, genome size, and the number of repeats were variable, indicating the non-conservative nature of plastome evolution within these sections. However, the genomes of the Lobelia sect. Delostemon and L. sect. Stenotium showed rearrangements, revealing that these two sections might have undergone different evolutionary histories. We assessed nine hotspot genes and 27-51 simple sequence repeat motifs, which will also serve as valuable DNA barcode regions in future population genetics studies and for the delineation of plant species. Our phylogenetic analysis resolved the evolutionary positions of the five sections in agreement with previous evolutionary trees based on morphological features. Although phylogenetic reconstruction of Lobelioideae based on the rpoC2 gene has rarely been performed, our results indicated that it contains a considerable amount of phylogenetic information and offers great promise for further phylogenetic analysis of Lobelioideae. Our site-specific model identified 173 sites under highly positive selections. The branch-site model exhibited 11 positive selection sites involving four genes in the East Asian branches. These four genes may play critical roles in the adaptation of East Asian taxa to diverse environments. Our study is the first to detect plastome organization, phylogenetic utility, and signatures of positive selection in the plastomes of East Asian lobelias, which will help to further advance taxonomic and evolutionary studies and the utilization of medicinal plant resources.

The dynamic history of plastome structure across aquatic subclass Alismatidae.

BACKGROUND: The rapidly increasing availability of complete plastomes has revealed more structural complexity in this genome under different taxonomic levels than expected, and this complexity provides important evidence for understanding the evolutionary history of angiosperms. To explore the dynamic history of plastome structure across the subclass Alismatidae, we sampled and compared 38 complete plastomes, including 17 newly assembled, representing all 12 recognized families of Alismatidae. RESULT: We found that plastomes size, structure, repeat elements, and gene content were highly variable across the studied species. Phylogenomic relationships among families were reconstructed and six main patterns of variation in plastome structure were revealed. Among these, the inversion from rbcL to trnV-UAC (Type I) characterized a monophyletic lineage of six families, but independently occurred also in Caldesia grandis. Three independent ndh gene loss events were uncovered across the Alismatidae. In addition, we detected a positive correlation between the number of repeat elements and the size of plastomes and IR in Alismatidae. CONCLUSION: In our study, ndh complex loss and repeat elements likely contributed to the size of plastomes in Alismatidae. Also, the ndh loss was more likely related to IR boundary changes than the adaptation of aquatic habits. Based on existing divergence time estimation, the Type I inversion may have occurred during the Cretaceous-Paleogene in response to the extreme paleoclimate changes. Overall, our findings will not only allow exploring the evolutionary history of Alismatidae plastome, but also provide an opportunity to test if similar environmental adaptations result in convergent restructuring in plastomes.

Evolutionary Patterns of the Chloroplast Genome in Vanilloid Orchids (Vanilloideae, Orchidaceae).

The Vanilloideae (vanilloids) is one of five subfamilies of Orchidaceae and is composed of fourteen genera and approximately 245 species. In this study, the six new chloroplast genomes (plastomes) of vanilloids (two Lecanorchis, two Pogonia, and two Vanilla species) were decoded, and then the evolutionary patterns of plastomes were compared to all available vanilloid plastomes. Pogonia japonica has the longest plastome, with 158,200 bp in genome size. In contrast, Lecanorchis japonica has the shortest plastome with 70,498 bp in genome size. The vanilloid plastomes have regular quadripartite structures, but the small single copy (SSC) region was drastically reduced. Two different tribes of Vanilloideae (Pogonieae and Vanilleae) showed different levels of SSC reductions. In addition, various gene losses were observed among the vanilloid plastomes. The photosynthetic vanilloids (Pogonia and Vanilla) showed signs of stage 1 degradation and had lost most of their ndh genes. The other three species (one Cyrotsia and two Lecanorchis), however, had stage 3 or stage 4 degradation and had lost almost all the genes in their plastomes, except for some housekeeping genes. The Vanilloideae were located between the Apostasioideae and Cypripedioideae in the maximum likelihood tree. A total of ten rearrangements were found among ten Vanilloideae plastomes when compared to the basal Apostasioideae plastomes. The four sub-regions of the single copy (SC) region shifted into an inverted repeat (IR) region, and the other four sub-regions of the IR region shifted into the SC regions. Both the synonymous (dS) and nonsynonymous (dN) substitution rates of IR in-cooperated SC sub-regions were decelerated, while the substitution rates of SC in-cooperated IR sub-regions were accelerated. A total of 20 protein-coding genes remained in mycoheterotrophic vanilloids. Almost all these protein genes show accelerated base substitution rates compared to the photosynthetic vanilloids. Two of the twenty genes in the mycoheterotrophic species faced strong "relaxed selection" pressure (p-value < 0.05).

Sequence Analysis of the Plastomes of Two Tibetan Medicinal Plants of the Family Papaveraceae.

BACKGROUND: With the rapid development of next-generation sequencing technology, more plants plastomes have been sequenced, further advancing species identification and phylogenetic studies. However, there are a few studies on the genetic and phylogenetic analysis of the plastomes of Dicranostigma lactucoides Hook. f. et Thoms. and Hypecoum leptocarpum Hook. f. et Thoms. METHODS: In this study, we sequenced and analyzed the plastomes of Dicranostigma lactucoides Hook. f. et Thoms. and Hypecoum leptocarpum Hook. f. et Thoms., and conducted a phylogenetic analysis using 13 related species. RESULTS: The results showed that the plastomes of both D. lactucoides and H. leptocarpum had a typical tetrad structure, with sizes of 166,819 bp and 163,282 bp, respectively. We annotated 133 genes for D. lactucoides and 120 genes for H. leptocarpum. A total of 72 and 43 simple repetitive sequences were detected in D. lactucoides and H. leptocarpum, respectively. Codon preference analysis showed that the relative usage frequency of codons and the relative abundance of synonymous codons used were the same for both plastomes. Nucleotide polymorphism analysis identified seven variant loci with high nucleotide diversity (Pi) values, all located in the large single copy (LSC) region. Inverted repeat (IR) boundary analysis revealed differences in gene types and locations on both sides of the boundary, except for the small single copy/inverted repeat a (SSC/IRa) boundary. The phylogenetic analysis showed the species clustered into two major groups, one with five genera (Hypecoum, Corydalis, Papaver, Meconopsis, and Dicranostigma) and the other with two genera (Coreanomecon; and Hylomecon). CONCLUSIONS: Comparative analysis of the plastome genomic characteristics and phylogeny of D. lactucoides and H. leptocarpum laid the foundation for identifying the above two species and the phylogenetic study and comprehensive exploitation of the Papaveraceae.

Comparative Analyses of Plastomes of Four Anubias (Araceae) Taxa, Tropical Aquatic Plants Endemic to Africa.

Anubias Schott (Araceae) have high ornamental properties as aquarium plants. However, the genus has difficulties in species identification, and the mechanism of its adaptation to the aquatic environment is unknown. To better identify species and understand the evolutionary history of Anubias, the plastomes of Anubias barteri Schott, A. barteri var. nana (Engl.) Crusio, and A. hastifolia Engl., were sequenced. The sizes of the plastomes of Anubias ranged from 169,841 bp to 170,037 bp. These plastomes were composed of conserved quadripartite circular structures and comprised 112 unique genes, including 78 protein-coding genes, 30 transfer RNA genes, and 4 ribosomal RNA genes. The comparative analysis of genome structure, repeat sequences, codon usage and RNA editing sites revealed high similarities among the Anubias plastomes, indicating the conservation of plastomes of Anubias. Three spacer regions with relatively high nucleotide diversity, trnL-CAA-ndhB, ycf1-ndhF, and rps15-ycf1, were found within the plastomes of Anubias. Phylogenetic analysis, based on 75 protein-coding genes, showed that Anubias was sister to Montrichardia arborescens (L.) Schott (BS = 99). In addition, four genes (ccsA, matK, ndhF, and ycf4) that contain sites undergoing positive selection were identified within the Anubias plastomes. These genes may play an important role in the adaptation of Anubias to the aquatic environment. The present study provides a valuable resource for further studies on species identification and the evolutionary history of Anubias.

Plastome phylogenomics and biogeography of the subfam. Polygonoideae (Polygonaceae).

Polygonaceae has a complex taxonomic history, although a few studies using plastid or nuclear DNA fragments have explored relationships within this family, intrafamilial relationships remain controversial. Here, we newly sequenced and annotated 17 plastomes representing 12 genera within Polygonaceae. Combined with previously published data, a total of 49 plastomes representing 22/46 Polygonaceae genera and 16/20 Polygonoideae genera were collected to infer the phylogeny of Polygonaceae, with an emphasis on Polygonoideae. Plastome comparisons revealed high conservation within Polygonoideae in structure and gene order. Phylogenetic analyses using both Maximum Likelihood and Bayesian methods revealed two major clades and seven tribes within Polygonoideae. BEAST and S-DIVA analyses suggested a Paleocene origin of Polygonoideae in Asia. While most genera of Polygonoideae originated and further diversified in Asia, a few genera experienced multiple long-distance dispersal events from Eurasia to North America after the Miocene, with a few dispersal events to the Southern Hemisphere also being detected. Both ancient vicariance and long-distance events have played important roles in shaping the current distribution pattern of Polygonoideae.