Plastome evolution and phylogeny of subtribe Aeridinae (Vandeae, Orchidaceae)

Characterization of Angraecum (Angraecinae, Orchidaceae) Plastomes and Utility of Sequence Variability Hotspots

Angraecum, commonly known as Darwin's orchid, is the largest genus of Angraecinae (Orchidaceae). This genus exhibits a high morphological diversity, making it as a good candidate for macroevolutionary studies. In this study, four complete plastomes of Angraecum were firstly reported and the potential variability hotspots were explored. The plastomes possessed the typical quadripartite structure and ranged from 150,743 to 151,818 base pair (bp), with a guanine-cytosine (GC) content of 36.6-36.9%. The plastomes all contained 120 genes, consisting of 74 protein-coding genes (CDS), 38 transfer RNA (tRNA) genes and 8 ribosomal RNA (rRNA) genes; all ndh genes were pseudogenized or lost. A total of 30 to 46 long repeats and 55 to 63 SSRs were identified. Relative synonymous codon usage (RSCU) analysis indicated a high degree of conservation in codon usage bias. The Ka/Ks ratios of most genes were lower than 1, indicating that they have undergone purifying selection. Based on the ranking of Pi (nucleotide diversity) values, five regions (trnSGCU-trnGGCC, ycf1-trnNGGU, trnNGUU-rpl32, psaC-ndhE and trnSGCU-trnGGCC) and five protein-coding genes (rpl32, rps16, psbK, rps8, and ycf1) were identified. The consistent and robust phylogenetic relationships of Angraecum were established based on a total of 40 plastomes from the Epidendroideae subfamily. The genus Angraecum was strongly supported as a monophyletic group and sister to Aeridinae. Our study provides an ideal system for investigating molecular identification, plastome evolution and DNA barcoding for Angraecum.

Comparative Analysis of Luisia (Aeridinae, Orchidaceae) Plastomes Shed Light on Plastomes Evolution and Barcodes Investigation

Luisia, a genus of the subtribe Aeridinae of Orchidaceae, comprises ca. 40 species. Members of Luisia exhibit unique morphological characteristics and represent a valuable ornamental orchid genus. However, due to the scarcity of distinct morphological characters, species identification within this genus is ambiguous and controversial. In the present study, next-generation sequencing (NGS) methods were used to assemble the plastomes of five Luisia species and compare them with one publicly available Luisia plastid genome data. The plastomes of Luisia possessed a quadripartite structure, with sizes ranging from 146,243 bp to 147,430 bp. The plastomes of six Luisia species contained a total of 120 genes, comprising 74 protein-coding genes, 38 tRNA genes and eight rRNA genes. Notably, all ndh genes were pseudogenized or lost. An analysis of codon usage bias showed that leucine (Leu) exhibited the highest frequency, while cysteine (Cys) exhibited the lowest frequency. A total of 57 to 64 SSRs and 42 to 49 long repeats were identified. Five regions and five coding sequences were identified for DNA barcodes, based on the nucleotide diversity (Pi) analysis. The species of Luisia constituted a monophyletic group and were sister to Paraphalaenopsis with strong support. Our study deepens the understanding of species identification, plastome evolution and the phylogenetic positions of Luisia.

Comparative and phylogenetic analysis of Chiloschista (Orchidaceae) species and DNA barcoding investigation based on plastid genomes

BACKGROUND

Chiloschista (Orchidaceae, Aeridinae) is an epiphytic leafless orchid that is mainly distributed in tropical or subtropical forest canopies. This rare and threatened orchid lacks molecular resources for phylogenetic and barcoding analysis. Therefore, we sequenced and assembled seven complete plastomes of Chiloschista to analyse the plastome characteristics and phylogenetic relationships and conduct a barcoding investigation.

RESULTS

We are the first to publish seven Chiloschista plastomes, which possessed the typical quadripartite structure and ranged from 143,233 bp to 145,463 bp in size. The plastomes all contained 120 genes, consisting of 74 protein-coding genes, 38 tRNA genes and eight rRNA genes. The ndh genes were pseudogenes or lost in the genus, and the genes petG and psbF were under positive selection. The seven Chiloschista plastomes displayed stable plastome structures with no large inversions or rearrangements. A total of 14 small inversions (SIs) were identified in the seven Chiloschista plastomes but were all similar within the genus. Six noncoding mutational hotspots (trnNGUU-rpl32 > rpoB-trnCGCA > psbK-psbI > psaC-rps15 > trnEUUC-trnTGGU > accD-psaI) and five coding sequences (ycf1 > rps15 > matK > psbK > ccsA) were selected as potential barcodes based on nucleotide diversity and species discrimination analysis, which suggested that the potential barcode ycf1 was most suitable for species discrimination. A total of 47-56 SSRs and 11-14 long repeats (> 20 bp) were identified in Chiloschista plastomes, and they were mostly located in the large single copy intergenic region. Phylogenetic analysis indicated that Chiloschista was monophyletic. It was clustered with Phalaenopsis and formed the basic clade of the subtribe Aeridinae with a moderate support value. The results also showed that seven Chiloschista species were divided into three major clades with full support.

CONCLUSION

This study was the first to analyse the plastome characteristics of the genus Chiloschista in Orchidaceae, and the results showed that Chiloschista plastomes have conserved plastome structures. Based on the plastome hotspots of nucleotide diversity, several genes and noncoding regions are suitable for phylogenetic and population studies. Chiloschista may provide an ideal system to investigate the dynamics of plastome evolution and DNA barcoding investigation for orchid studies.

Comparative Phylogenetic Analysis for Aerides (Aeridinae, Orchidaceae) Based on Six Complete Plastid Genomes

Aerides Lour. (Orchidaceae, Aeridinae) is a group of epiphytic orchids with high ornamental value, mainly distributed in tropical and subtropical forests, that comprises approximately 20 species. The species are of great value in floriculture and garden designing because of their beautiful flower shapes and colors. Although the morphological boundaries of Aerides are clearly defined, the relationship between Aerides and other closely related genera is still ambiguous in terms of phylogeny. To better understand their phylogenetic relationships, this study used next-generation sequencing technology to investigate the phylogeny and DNA barcoding of this taxonomic unit using genetic information from six Aerides plastid genomes. The quadripartite-structure plastomes ranged from 147,244 bp to 148,391 bp and included 120 genes. Among them, 74 were protein coding genes, 38 were tRNA genes and 8 were rRNA genes, while the ndh genes were pseudogenized or lost. Four non-coding mutational hotspots (rpl20-rpl33, psbM, petB, rpoB-trnCGCA, Pi > 0.06) were identified. A total of 71-77 SSRs and 19-46 long repeats (>30 bp) were recognized in Aerides plastomes, which were mostly located in the large single-copy region. Phylogenetic analysis indicated that Aerides was monophylic and sister to Renanthera. Moreover, our results confirmed that six Aerides species can be divided into three major clades. These findings provide assistance for species identification and DNA barcoding investigation in Aerides, as well as contributes to future research on the phylogenomics of Orchidaceae.

Complete Chloroplast Genomes and Comparative Analyses of Three Paraphalaenopsis (Aeridinae, Orchidaceae) Species

Paraphalaenopsis, a genus of perennial herbs from the family Orchidaceae, contains a number of ornamental species. However, there is no information on the chloroplast genomes of Paraphalaenopsis, which limits our studies of this genus. In this study, we reported the chloroplast genomes of three species of Paraphalaenopsis (P. labukensis, P. denevel, and P. laycockii 'Semi-alba') and performed comprehensive comparative analysis. These three chloroplast genomes showed a typical quadripartile structure. Their lengths ranged from 147,311 bp to 149,240 bp. Each genome contained 120 unique genes, including 74 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. Comparative analysis revealed major differences in sequence divergence in the three chloroplast genomes. In addition, six hypervariable regions were identified (psbM-trnD^GUC, psbB, ccsA, trnK^UUU, trnS^GCU-trnG^UCC, rps16-trnQ^UUG) that can be used as DNA molecular markers. Phylogenetic relationships were determined using the chloroplast genomes of 28 species from 12 genera of Aeridinae. Results suggested that Paraphalaenopsis was a clade of Aeridinae that was sister to the Holcoglossum-Vanda clade, with 100% bootstrap support within Aeridinae. The findings of this study provided the foundation for future studies on the phylogenetic analysis of Aeridinae.

Complete chloroplast genome structural characterization of two Phalaenopsis (Orchidaceae) species and comparative analysis with their alliance

BACKGROUND

The taxonomy and infrageneric delimitation of Phalaenopsis Blume has been significantly disputed due to some overlapping morphological features between species related, which needed further evidence for clarification. The structural characterization of complete chloroplast genomes of P. storbatiana and P. wilsonii were analyzed and compared with those of related taxa to provide a better understanding of their genomic information on taxonomy and phylogeny.

RESULTS

It was shown that chloroplast genomes of Phalaenopsis storbatiana and P. wilsonii had a typical quadripartite structure with conserved genome arrangements and moderate divergence. The chloroplast genomes of P. storbatiana and P. wilsonii were 145,885 bp and 145,445 bp in length, respectively, and shared a similar GC content of 36.8%. Gene annotations of two species revealed 109 single-copy genes consistently. In addition, 20 genes duplicated in the inverted regions, 16 genes each possessed one or more introns, and five ndh (NA (D)H dehydrogenase) genes were observed in both. Comparative analysis of the total cp genomes of P. storbatiana and P. wilsonii with those of other six related Phalaenopsis species confirmed the stable sequence identity for coding and non-coding regions and higher sequence variation in SC regions than IR regions. Most of their protein-coding genes had a high degree of codon preference. Moreover, 45 genes were discovered with significantly positive selection. However, different amplifications in IR regions were observed in these eight species. Phylogenetic analysis based on CDS from 60 species representing main clades in Orchidaceae indicated that Phalaenopsis species including P. stobartiana and P. wilsonii formed a monophyletic clade with high bootstrap nested in tribe Vandeae of Epidendroideae, which was consistent with those from previous studies.

CONCLUSIONS

The results could provide insight into understanding the plastome evolution and phylogenetic relationships of Phalaenopsis.

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).

Comparative phylogenetic analysis of complete plastid genomes of Renanthera (Orchidaceae)

Owing to its attractive flower shape and color, Renanthera (Orchidaceae), comprising about 19 species, has significant ornamental value as a houseplant, in floral design and in landscape gardens. Two species of Renanthera are categorized as endangered and critically endangered in China’s Red List and international trade in these orchids is currently strictly monitored by the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES). This paper reports on the de novo assembled and annotated plastome of four species of Renanthera; R. citrina, R. coccinea, R. imschootiana, and R. philippinensis. The length of the plastome sequences ranged from 144,673 bp (R. imschootiana) to 149,007 bp (R. coccinea) with GC content of 36.6–36.7%. The plastomes showed a typical quadripartite structure, including a large single-copy (84,241–86,404 bp), a small single-copy (11,468–12,167 bp), and a pair of inverted repeats (24,482–25,715 bp) regions. Of the 120 genes detected, 74 were protein coding, 38 were tRNA, and eight were rRNA genes. The plastome of Renanthera is rather conserved, but nucleotide variations that could distinguish them apart are noticeable—the total number of tandem repeats ranged from 62 (in R. imschootiana) to 74 (in R. citrina); while the number of long repeats ranged from 21 (in R. imschootiana and R. philippinensis) to 43 (in R. citrina). Three hypervariable regions (psbI-trnS-GCU, trnG-GCC, rpl32) were identified. Phylogenetic analyses based on the CDS using maximum likelihood (ML) and Bayesian inference (BI) revealed that Renanthera is closely related to Holcoglossum, Neofinetia, Pendulorchis, and Vanda. The relationship between the four species of Renanthera was fully resolved; a monophyletic clade was formed and R. coccinea was recorded as the first to diverge from the rest. The genetic data obtained from this study could serve as a useful resource for species identification in Renanthera as well as contribute to future research on the phylogenomics of Orchidaceae.

Comparative and phylogenetic analyses of six Kenya Polystachya (Orchidaceae) species based on the complete chloroplast genome sequences

BACKGROUND

Polystachya Hook. is a large pantropical orchid genus (c. 240 species) distributed in Africa, southern Asia and the Americas, with the center of diversity in Africa. Previous studies on species of this genus have not obtained the complete chloroplast genomes, structures and variations. Additionally, the phylogenetic position of the genus in the Orchidaceae is still controversial and uncertain. Therefore, in this study, we sequenced the complete plastomes of six Kenya Polystachya species based on genome skimming, subjected them to comparative genomic analysis, and reconstructed the phylogenetic relationships with other Orchidaceae species.

RESULTS

The results exhibited that the chloroplast genomes had a typical quadripartite structure with conserved genome arrangement and moderate divergence. The plastomes of the six Polystachya species ranged from 145,484 bp to 149,274 bp in length and had an almost similar GC content of 36.9-37.0%. Gene annotation revealed 106-109 single-copy genes. In addition, 19 genes are duplicated in the inverted regions, and 16 genes each possessd one or more introns. Although no large structural variations were observed among the Polystachya plastomes, about 1 kb inversion was found in Polystachya modesta and all 11 ndh genes in the Polystachya plastomes were lost or pseudogenized. Comparative analysis of the overall sequence identity among six complete chloroplast genomes confirmed that for both coding and non-coding regions in Polystachya, SC regions exhibit higher sequence variation than IRs. Furthermore, there were various amplifications in the IR regions among the six Polystachya species. Most of the protein-coding genes of these species had a high degree of codon preference. We screened out SSRs and found seven relatively highly variable loci. Moreover, 13 genes were discovered with significant positive selection. Phylogenetic analysis showed that the six Polystachya species formed a monophyletic clade and were more closely related to the tribe Vandeae. Phylogenetic relationships of the family Orchidaceae inferred from the 85 chloroplast genome sequences were generally consistent with previous studies and robust.

CONCLUSIONS

Our study is the initial report of the complete chloroplast genomes of the six Polystachya species, elucidates the structural characteristics of the chloroplast genome of Polystachya, and filters out highly variable sequences that can contribute to the development of DNA markers for use in the study of genetic variability and evolutionary studies in Polystachya. In addition, the phylogenetic results strongly support that the genus of Polystachya is a part of the tribe Vandeae.

Chloroplast genomic diversity in Bulbophyllum section Macrocaulia (Orchidaceae, Epidendroideae, Malaxideae): Insights into species divergence and adaptive evolution

Bulbophyllum is the largest genus in Orchidaceae with a pantropical distribution. Due to highly significant diversifications, it is considered to be one of the most taxonomically and phylogenetically complex taxa. The diversification pattern and evolutionary adaptation of chloroplast genomes are poorly understood in this species-rich genus, and suitable molecular markers are necessary for species determination and phylogenetic analysis. A natural Asian section Macrocaulia was selected to estimate the interspecific divergence of chloroplast genomes in this study. Here, we sequenced the complete chloroplast genome of four Bulbophyllum species, including three species from section Macrocaulia. The four chloroplast genomes had a typical quadripartite structure with a genome size ranged from 156,182 to 158,524 bp. The chloroplast genomes included 113 unique genes encoding 79 proteins, 30 tRNAs and 4 rRNAs. Comparison of the four chloroplast genomes showed that the three species from section Macrocaulia had similar structure and gene contents, and shared a number of indels, which mainly contribute to its monophyly. In addition, interspecific divergence level was also great. Several exclusive indels and polymorphism SSR loci might be used for taxonomical identification and determining interspecific polymorphisms. A total of 20 intergenic regions and three coding genes of the most variable hotspot regions were proposed as candidate effective molecular markers for future phylogenetic relationships at different taxonomical levels and species divergence in Bulbophyllum. All of chloroplast genes in four Bulbophyllum species were under purifying selection, while 13 sites within six genes exhibited site-specific selection. A whole chloroplast genome phylogenetic analysis based on Maximum Likelihood, Bayesian and Parsimony methods all supported the monophyly of section Macrocaulia and the genus of Bulbophyllum. Our findings provide valuable molecular markers to use in accurately identifying species, clarifying taxonomy, and resolving the phylogeny and evolution of the genus Bulbophyllum. The molecular markers developed in this study will also contribute to further research of conservation of Bulbophyllum species.

The complete chloroplast genome sequence of Oxystophyllum changjiangense (Orchidaceae)

Oxystophyllum changjiangense has high economic value due to its wide applications in horticultural and medicinal fields. Here, the first chloroplast genome of O. changjiangense was sequenced and reported. The chloroplast genome displayed the typical quadripartite structure containing a pair of inverted repeats (IR), a long single-copy region (LSC), and a short single-copy region (SSC). Total 110 genes were found including 76 protein-coding genes, 30 tRNA genes, and four rRNA genes. Phylogeny analysis showed O. changjiangense has a close relationship with Phaius species.

Plastome evolution in the Caesalpinia group (Leguminosae) and its application in phylogenomics and populations genetics

The chloroplast genomes of Caesalpinia group species are structurally conserved, but sequence level variation is useful for both phylogenomic and population genetic analyses. Variation in chloroplast genomes (plastomes) has been an important source of information in plant biology. The Caesalpinia group has been used as a model in studies correlating ecological and genomic variables, yet its intergeneric and infrageneric relationships are not fully solved, despite densely sampled phylogenies including nuclear and plastid loci by Sanger sequencing. Here, we present the de novo assembly and characterization of plastomes from 13 species from the Caesalpinia group belonging to eight genera. A comparative analysis was carried out with 13 other plastomes previously available, totalizing 26 plastomes and representing 15 of the 26 known Caesalpinia group genera. All plastomes showed a conserved quadripartite structure and gene repertoire, except for the loss of four ndh genes in Erythrostemon gilliesii. Thirty polymorphic regions were identified for inter- or intrageneric analyses. The 26 aligned plastomes were used for phylogenetic reconstruction, revealing a well-resolved topology, and dividing the Caesalpinia group into two fully supported clades. Sixteen microsatellite (cpSSR) loci were selected from Cenostigma microphyllum for primer development and at least two were cross-amplified in different Leguminosae subfamilies by in vitro or in silico approaches. Four loci were used to assess the genetic diversity of C. microphyllum in the Brazilian Caatinga. Our results demonstrate the structural conservation of plastomes in the Caesalpinia group, offering insights into its systematics and evolution, and provides new genomic tools for future phylogenetic, population genetics, and phylogeographic studies.