Plastomes of Bletilla (Orchidaceae) and Phylogenetic Implications

Phylogenomics and plastomics offer new evolutionary perspectives on Kalanchoideae (Crassulaceae)

BACKGROUND AND AIMS

Kalanchoideae is one of three subfamilies within Crassulaceae and contains four genera. Despite previous efforts, the phylogeny of Kalanchoideae remains inadequately resolved with persistent issues including low support, unstructured topologies and polytomies. This study aimed to address two central objectives: (1) resolving the pending phylogenetic questions within Kalanchoideae by using organelle-scale 'barcodes' (plastomes) and nuclear data; and (2) investigating interspecific diversity patterns among Kalanchoideae plastomes.

METHODS

To explore the plastome evolution in Kalanchoideae, we newly sequenced 38 plastomes representing all four constituent genera (Adromischus, Cotyledon, Kalanchoe and Tylecodon). We performed comparative analyses of plastomic features, including GC and gene contents, gene distributions at the IR (inverted repeat) boundaries, nucleotide divergence, plastomic tRNA (pttRNA) structures and codon aversions. Additionally, phylogenetic inferences were inferred using both the plastomic dataset (79 genes) and nuclear dataset (1054 genes).

KEY RESULTS

Significant heterogeneities were observed in plastome lengths among Kalanchoideae, strongly correlated with LSC (large single copy) lengths. Informative diversities existed in the gene content at SSC/IRa (small single copy/inverted repeat a), with unique patterns individually identified in Adromischus leucophyllus and one major Kalanchoe clade. The ycf1 gene was assessed as a shared hypervariable region among all four genera, containing nine lineage-specific indels. Three pttRNAs exhibited unique structures specific to Kalanchoideae and the genera Adromischus and Kalanchoe. Moreover, 24 coding sequences revealed a total of 41 lineage-specific unused codons across all four constituent genera. The phyloplastomic inferences clearly depicted internal branching patterns in Kalanchoideae. Most notably, by both plastid- and nuclear-based phylogenies, our research offers the first evidence that Kalanchoe section Eukalanchoe is not monophyletic.

CONCLUSIONS

This study conducted comprehensive analyses on 38 newly reported Kalanchoideae plastomes. Importantly, our results not only reconstructed well-resolved phylogenies within Kalanchoideae, but also identified highly informative unique markers at the subfamily, genus and species levels. These findings significantly enhance our understanding of the evolutionary history of Kalanchoideae.

Comparative plastome analysis of the sister genera Ceratocephala and Myosurus (Ranunculaceae) reveals signals of adaptive evolution to arid and aquatic environments

BACKGROUND

Expansion and contraction of inverted repeats can cause considerable variation of plastid genomes (plastomes) in angiosperms. However, little is known about whether structural variations of plastomes are associated with adaptation to or occupancy of new environments. Moreover, adaptive evolution of angiosperm plastid genes remains poorly understood. Here, we sequenced the complete plastomes for four species of xerophytic Ceratocephala and hydrophytic Myosurus, as well as Ficaria verna. By an integration of phylogenomic, comparative genomic, and selection pressure analyses, we investigated evolutionary patterns of plastomes in Ranunculeae and their relationships with adaptation to dry and aquatic habitats.

RESULTS

Owing to the significant contraction of the boundary of IRA/LSC towards the IRA, plastome sizes and IR lengths of Myosurus and Ceratocephala are smaller within Ranunculeae. Compared to other Ranunculeae, the Myosurus plastome lost clpP and rps16, one copy of rpl2 and rpl23, and one intron of rpoC1 and rpl16, and the Ceratocephala plastome added an infA gene and lost one copy of rpl2 and two introns of clpP. A total of 11 plastid genes (14%) showed positive selection, two genes common to Myosurus and Ceratocephala, seven in Ceratocephala only, and two in Myosurus only. Four genes showed strong signals of episodic positive selection. The rps7 gene of Ceratocephala and the rpl32 and ycf4 genes of Myosurus showed an increase in the rate of variation close to 3.3 Ma.

CONCLUSIONS

The plastomic structure variations as well as the positive selection of two plastid genes might be related to the colonization of new environments by the common ancestor of Ceratocephala and Myosurus. The seven and two genes under positive selection might be related to the adaptation to dry and aquatic habitats in Ceratocephala and Myosurus, respectively. Moreover, intensified aridity and frequent sea-level fluctuations, as well as global cooling, might have favored an increased rate of change in some genes at about 3.3 Ma, associated with adaptation to dry and aquatic environments, respectively. These findings suggest that changing environments might have influenced structural variations of plastomes and fixed new mutations arising on some plastid genes owing to adaptation to specific habitats.

A new perspective on codon usage, selective pressure, and phylogenetic implications of the plastomes in the Telephium clade (Crassulaceae)

The Telephium clade of the Crassulaceae family contains many medicinal, ornamental, and ecologically restorative plants. However, the phylogenetic relationships within the clade remain debated, and comprehensive analyses of codon usage and selection pressure in Telephium plastomes are limited. In this study, we assembled and annotated four plastomes and performed extensive analyses. The plastomes exhibited a typical quadripartite structure and high conservation. The lengths ranged from 151,357 bp to 151,641 bp with 134 genes identified. The GC content was the highest within IR, followed by LSC, and lowest in the SSC region. Meanwhile, a unique inversion was observed within the LSC region of Meterostachys sikokianus. Polymorphisms analysis revealed minimum nucleotide diversity in the IR regions, with over ten highly polymorphic regions identified. Phylogenetically, two subclades formed within the monophyletic Telephium clade, with Umbilicus as the sister group to the remaining Hylotelephium subclade members. Notably, no significant positive selection was found among the 79 plastid genes, which showed varying evolutionary patterns. However, 19 genes contained codons under positive selection. The specific functions of these sites require further investigation. Synonymous codon usage was biased and conserved across the tested plastomes, shaped by natural selection, mutations and other factors of varying influence. We also identified 34 taxon-specific codon aversion motifs from 49 plastid genes. Our plastomic analyses elucidate phylogenetic relationships and evolutionary patterns in this medicinal clade, providing a foundation for further research on these ecologically and pharmaceutically important plants.

Comparative chloroplast genomics of 24 species shed light on the genome evolution and phylogeny of subtribe Coelogyninae (Orchidaceae)

BACKGROUND

The orchids of the subtribe Coelogyninae are among the most morphologically diverse and economically important groups within the subfamily Epidendroideae. Previous molecular studies have revealed that Coelogyninae is an unambiguously monophyletic group. However, intergeneric and infrageneric relationships within Coelogyninae are largely unresolved. There has been long controversy over the classification among the genera within the subtribe.

RESULTS

The complete chloroplast (cp.) genomes of 15 species in the subtribe Coelogyninae were newly sequenced and assembled. Together with nine available cp. genomes in GenBank from representative clades of the subtribe, we compared and elucidated the characteristics of 24 Coelogyninae cp. genomes. The results showed that all cp. genomes shared highly conserved structure and contained 135 genes arranged in the same order, including 89 protein-coding genes, 38 tRNAs, and eight rRNAs. Nevertheless, structural variations in relation to particular genes at the IR/SC boundary regions were identified. The diversification pattern of the cp. genomes showed high consistency with the phylogenetic placement of Coelogyninae. The number of different types of SSRs and long repeats exhibited significant differences in the 24 Coelogyninae cp. genomes, wherein mononucleotide repeats (A/T), and palindromic repeats were the most abundant. Four mutation hotspot regions (ycf1a, ndhF-rp132, psaC-ndhE, and rp132-trnL) were determined, which could serve as effective molecular markers. Selection pressure analysis revealed that three genes (ycf1a, rpoC2 and ycf2 genes) might have experienced apparent positive selection during the evolution. Using the alignments of whole cp. genomes and protein-coding sequences, this study presents a well-resolved phylogenetic framework of Coelogyninae.

CONCLUSION

The inclusion of 55 plastid genome data from a nearly complete generic-level sampling provide a comprehensive view of the phylogenetic relationships among genera and species in subtribe Coelogyninae and illustrate the diverse genetic variation patterns of plastid genomes in this species-rich plant group. The inferred relationships and informally recognized major clades within the subtribe are presented. The genetic markers identified here will facilitate future studies on the genetics and phylogeny of subtribe Coelogyninae.

The Complete Chloroplast Genomes of Blepharoglossum elegans and B. grossum and Comparative Analysis with Related Species (Orchidaceae, Malaxideae)

Blepharoglossum is a rare orchid genus of the Malaxidinae primarily distributed in tropical Pacific islands, with several species occurring in the Taiwan and Hainan Islands of China. Currently, the monophyletic status of Blepharoglossum has been challenged, and the phylogenetic relationships among its allied groups have remained unresolved with traditional DNA markers. In this study, we initially sequenced and annotated the chloroplast (cp) genomes of two Blepharoglossum species, Blepharoglossum elegans (Lindl.) L. Li and Blepharoglossum grossum (Rchb.f.) L. Li. These cp genomes of Blepharoglossum share the typical quadripartite and circular structure. Each of the genomes encodes a total of 133 functional genes, including 87 protein-coding genes (CDS), 38 tRNA genes and 8 rRNA genes. By comparing the sequence differences between these two cp genomes, it was found that they are relatively conserved in terms of overall gene content and gene arrangement. However, a total of 684 SNPs and 2664 indels were still identified, with ycf1, clpP, and trnK-UUU protein-coding genes having the highest number of SNPs and indels. In further comparative analyses among the six cp genomes in Malaxidinae, significant sequence divergences were identified in the intergenic regions, namely rps16-trnQ-UUG, trnS-GCU-trnG-GCC, rpoB-trnC-GCA, trnE-UUC-trnT-GGU, trnF-GAA-trnV-UAC, atpB-rbcL, petA-psbJ, psbE-petL, psbB-psbT, trnN-GUU-rpl32, trnV-GAC-rps7, and rps7-trnL-CAA, and five coding regions, including matK, and rpoC2, ycf1, and two ycf2 genes. Phylogenetic analysis indicated that Blepharoglossum and Oberonia form a highly supported sister group relationship. Our results are consistent with previous studies and present increased resolution among major clades.

The characteristics of the complete chloroplast genome of Staurogyne concinnula (Hance) O. Kuntze (Acanthaceae)

Staurogyne concinnula (Hance) O. Kuntze (Acanthaceae) is an important ornamental herb mainly distributed in the southern region of China, including Fujian, Guangdong, Hainan, and Taiwan provinces. However, the complete chloroplast genome of S. concinnula, which could serve as a genetic resource for studies on its taxonomy and evolution, is poorly studied at present. In this study, we reported the complete chloroplast genome of S. concinnula that was assembled using high-throughput sequencing data. The chloroplast genome was 153,783 bp long, with a typical quadripartite structure containing a small single-copy region (SSC; 17,855 bp), a large single-copy region (LSC; 84,636 bp) and a pair of inverted repeats (IRs; each 25,646 bp). The overall GC content of the chloroplast genome was 38.04%. A total of 86 protein-coding genes (PCGs), 8 rRNA genes, and 37 tRNA genes were predicted. Phylogenetic analysis based on the combined sequences of 86 PCGs with the other 16 closely related species of Acanthaceae indicated that S. concinnula is closely related to Avicennia marina. The genomic data and finding from the phylogenetic studies of S. concinnula could provide useful information and give light to in-depth studies on the evolution pattern of the understudied species, as well as Staurogyne.

Structural Diversities and Phylogenetic Signals in Plastomes of the Early-Divergent Angiosperms: A Case Study in Saxifragales

As representative of the early-divergent groups of angiosperms, Saxifragales is extremely divergent in morphology, comprising 15 families. Within this order, our previous case studies observed significant structural diversities among the plastomes of several lineages, suggesting a possible role in elucidating their deep phylogenetic relationships. Here, we collected 208 available plastomes from 11 constituent families to explore the evolutionary patterns among Saxifragales. With thorough comparisons, the losses of two genes and three introns were found in several groups. Notably, 432 indel events have been observed from the introns of all 17 plastomic intron-containing genes, which could well play an important role in family barcoding. Moreover, numerous heterogeneities and strong intrafamilial phylogenetic implications were revealed in pttRNA (plastomic tRNA) structures, and the unique structural patterns were also determined for five families. Most importantly, based on the well-supported phylogenetic trees, evident phylogenetic signals were detected in combinations with the identified pttRNAs features and intron indels, demonstrating abundant lineage-specific characteristics for Saxifragales. Collectively, the results reported here could not only provide a deeper understanding into the evolutionary patterns of Saxifragales, but also provide a case study for exploring the plastome evolution at a high taxonomic level of angiosperms.

Ten Plastomes of Crassula (Crassulaceae) and Phylogenetic Implications

The genus Crassula is the second-largest genus in the family Crassulaceae, with about 200 species. As an acknowledged super-barcode, plastomes have been extensively utilized for plant evolutionary studies. Here, we first report 10 new plastomes of Crassula. We further focused on the structural characterizations, codon usage, aversion patterns, and evolutionary rates of plastomes. The IR junction patterns-IRb had 110 bp expansion to rps19-were conservative among Crassula species. Interestingly, we found the codon usage patterns of matK gene in Crassula species are unique among Crassulaceae species with elevated ENC values. Furthermore, subgenus Crassula species have specific GC-biases in the matK gene. In addition, the codon aversion motifs from matK, pafI, and rpl22 contained phylogenetic implications within Crassula. The evolutionary rates analyses indicated all plastid genes of Crassulaceae were under the purifying selection. Among plastid genes, ycf1 and ycf2 were the most rapidly evolving genes, whereas psaC was the most conserved gene. Additionally, our phylogenetic analyses strongly supported that Crassula is sister to all other Crassulaceae species. Our findings will be useful for further evolutionary studies within the Crassula and Crassulaceae.