Plastome Evolution and Phylogeny of Orchidaceae, With 24 New Sequences

Genome of the early spider-orchid Ophrys sphegodes provides insights into sexual deception and pollinator adaptation

Pollinator-driven evolution of floral traits is thought to be a major driver of angiosperm speciation and diversification. Ophrys orchids mimic female insects to lure male pollinators into pseudocopulation. This strategy, called sexual deception, is species-specific, thereby providing strong premating reproductive isolation. Identifying the genomic architecture underlying pollinator adaptation and speciation may shed light on the mechanisms of angiosperm diversification. Here, we report the 5.2 Gb chromosome-scale genome sequence of Ophrys sphegodes. We find evidence for transposable element expansion that preceded the radiation of the O. sphegodes group, and for gene duplication having contributed to the evolution of chemical mimicry. We report a highly differentiated genomic candidate region for pollinator-mediated evolution on chromosome 2. The Ophrys genome will prove useful for investigations into the repeated evolution of sexual deception, pollinator adaptation and the genomic architectures that facilitate evolutionary radiations.

Plastid genome evolution in leafless members of the orchid subfamily Orchidoideae, with a focus on Degranvillea dermaptera

PREMISE

Leafless, heterotrophic plants are prime examples of organismal modification, the genomic consequences of which have received considerable interest. In particular, plastid genomes (plastomes) are being sequenced at a high rate, allowing continual refinement of conceptual models of reductive evolution in heterotrophs. However, numerous sampling gaps exist, hindering the ability to conduct comprehensive phylogenomic analyses in these plants.

METHODS

Using floral tissue from an herbarium specimen, we sequenced and analyzed the plastome of Degranvillea dermaptera, a rarely collected, leafless orchid species from South America about which little is known, including its phylogenetic affinities.

RESULTS

The plastome is the most reduced of those sequenced among the orchid subfamily Orchidoideae. In Degranvillea, it has lost the majority of genes found in leafy autotrophic species, is structurally rearranged, and has similar gene content to the most reduced plastomes among the orchids. We found strong evidence for the placement of Degranvillea within the subtribe Spiranthinae using models that explicitly account for heterotachy, or lineage-specific evolutionary rate variation over time. We further found evidence of relaxed selection on several genes and of correlations among substitution rates and several other "traits" of the plastome among leafless members of orchid subfamily Orchidoideae.

CONCLUSIONS

Our findings advance knowledge on the phylogenetic relationships and paths of plastid genome evolution among the orchids, which have experienced more independent transitions to heterotrophy than any other plant family. This study demonstrates the importance of herbarium collections in comparative genomics of poorly known species of conservation concern.

Plastid phylogenomics reveals evolutionary relationships in the mycoheterotrophic orchid genus Dipodium and provides insights into plastid gene degeneration

The orchid genus Dipodium R.Br. (Epidendroideae) comprises leafy autotrophic and leafless mycoheterotrophic species, with the latter confined to sect. Dipodium. This study examined plastome degeneration in Dipodium in a phylogenomic and temporal context. Whole plastomes were reconstructed and annotated for 24 Dipodium samples representing 14 species and two putatively new species, encompassing over 80% of species diversity in sect. Dipodium. Phylogenomic analysis based on 68 plastid loci including a broad outgroup sampling across Orchidaceae found that sect. Leopardanthus is the sister lineage to sect. Dipodium. Dipodium ensifolium, the only leafy autotrophic species in sect. Dipodium, was found to be a sister to all leafless, mycoheterotrophic species, supporting a single evolutionary origin of mycoheterotrophy in the genus. Divergence-time estimations found that Dipodium arose ca. 33.3 Ma near the lower boundary of the Oligocene and that crown diversification commenced in the late Miocene, ca. 11.3 Ma. Mycoheterotrophy in the genus was estimated to have evolved in the late Miocene, ca. 7.3 Ma, in sect. Dipodium. The comparative assessment of plastome structure and gene degradation in Dipodium revealed that plastid ndh genes were pseudogenised or physically lost in all Dipodium species, including in leafy autotrophic species of both Dipodium sections. Levels of plastid ndh gene degradation were found to vary among species as well as within species, providing evidence of relaxed selection for retention of the NADH dehydrogenase complex within the genus. Dipodium exhibits an early stage of plastid genome degradation, as all species were found to have retained a full set of functional photosynthesis-related genes and housekeeping genes. This study provides important insights into plastid genome degradation along the transition from autotrophy to mycoheterotrophy in a phylogenomic and temporal context.

Phylogenetic position and plastid genome structure of Vietorchis, a mycoheterotrophic genus of Orchidaceae (subtribe Orchidinae) endemic to Vietnam

The orchid genus Vietorchis comprises three species, all discovered in the 21 century. Each of these species is achlorophyllous, mycoheterotrophic and is known to be endemic to Vietnam. The type species of the genus, V. aurea, occurs in a single location in northern Vietnam within a lowland limestone karstic area. Vietorchis furcata and V. proboscidea, in contrast, are confined to mountains of southern Vietnam, far away from any limestone formations. Taxonomic placement of Vietorchis remained uncertain for the reason of inconclusive morphological affinities. At the same time, the genus has never been included into molecular phylogenetic studies. We investigate the phylogenetic relationships of two species of Vietorchis (V. aurea and V. furcata) based on three DNA datasets: (1) a dataset comprising two nuclear regions, (2) a dataset comprising two plastid regions, and (3) a dataset employing data on the entire plastid genomes. Our phylogenetic reconstructions support the placement of Vietorchis into the subtribe Orchidinae (tribe Orchideae, subfamily Orchidoideae). This leads to a conclusion that the previously highlighted similarities in the rhizome morphology between Vietorchis and certain mycoheterotrophic genera of the subfamilies Epidendroideae and Vanilloideae are examples of a convergence. Vietorchis is deeply nested within Orchidinae, and therefore the subtribe Vietorchidinae is to be treated as a synonym of Orchidinae. In the obtained phylogenetic reconstructions, Vietorchis is sister to the photosynthetic genus Sirindhornia. Sirindhornia is restricted to limestone mountains, which allows to speculate that association with limestone karst is plesiomorphic for Vietorchis. Flower morphology is concordant with the molecular data in placing Vietorchis into Orchidinae and strongly supports the assignment of the genus to one of the two major clades within this subtribe. Within this clade, however, Vietorchis shows no close structural similarity with any of its genera; in particular, the proximity between Vietorchis and Sirindhornia has never been proposed. Finally, we assembled the plastid genome of V. furcata, which is 65969 bp long and contains 45 unique genes, being one of the most reduced plastomes in the subfamily Orchidoideae. The plastome of Vietorchis lacks any rearrangements in comparison with the closest studied autotrophic species, and possesses substantially contracted inverted repeats. No signs of positive selection acting on the protein-coding plastid sequences were detected.

Partial mycoheterotrophy in the leafless orchid Eulophia zollingeri specialized on wood-decaying fungi

Although the absence of normal leaves is often considered a sign of full heterotrophy, some plants remain at least partially autotrophic despite their leafless habit. Leafless orchids with green stems and capsules probably represent a late evolutionary stage toward full mycoheterotrophy and serve as valuable models for understanding the pathways leading to this nutritional strategy. In this study, based on molecular barcoding and isotopic analysis, we explored the physiological ecology of the leafless orchid Eulophia zollingeri, which displays green coloration, particularly during its fruiting phase. Although previous studies had shown that E. zollingeri, in its adult stage, is associated with Psathyrellaceae fungi and exhibits high 13C isotope signatures similar to fully mycoheterotrophic orchids, it remained uncertain whether this symbiotic relationship is consistent throughout the orchid's entire life cycle and whether the orchid relies exclusively on mycoheterotrophy for its nutrition during the fruiting season. Our study has demonstrated that E. zollingeri maintains a specialized symbiotic relationship with Psathyrellaceae fungi throughout all life stages. However, isotopic analysis and chlorophyll data have shown that the orchid also engages in photosynthesis to meet its carbon needs, particularly during the fruiting stage. This research constitutes the first discovery of partial mycoheterotrophy in leafless orchids associated with saprotrophic non-rhizoctonia fungi.

Organelle Genomes of Epipogium roseum Provide Insight into the Evolution of Mycoheterotrophic Orchids

Epipogium roseum, commonly known as one of the ghost orchids due to its rarity and almost transparent color, is a non-photosynthetic and fully mycoheterotrophic plant. Given its special nutritional strategies and evolutionary significance, the mitogenome was first characterized, and three plastomes sampled from Asia were assembled. The plastomes were found to be the smallest among Orchidaceae, with lengths ranging from 18,339 to 19,047 bp, and exhibited high sequence variety. For the mitogenome, a total of 414,552 bp in length, comprising 26 circular chromosomes, were identified. A total of 54 genes, including 38 protein-coding genes, 13 tRNA genes, and 3 rRNA genes, were annotated. Multiple repeat sequences spanning a length of 203,423 bp (45.47%) were discovered. Intriguingly, six plastid regions via intracellular gene transfer and four plastid regions via horizontal gene transfer to the mitogenome were observed. The phylogenomics, incorporating 90 plastomes and 56 mitogenomes, consistently revealed the sister relationship of Epipogium and Gastrodia, with a bootstrap percentage of 100%. These findings shed light on the organelle evolution of Orchidaceae and non-photosynthetic plants.

The complete chloroplast genome of Chrysoglossum ornatum (Epidendroideae, Orchidaceae) and its phylogenetic analysis

Chrysoglossum ornatum Blume, the type species of Chrysoglossum Blume, belongs to the tribe Collabieae of the subfamily Epidendroideae of Orchidaceae. In this study, we sequenced, assembled, and analyzed the complete chloroplast genome of C. ornatum. The result showed that the complete chloroplast genome of C. ornatum was 158,175 bp in size, consisting of a large single-copy (LSC) region of 87,235 bp, a small single-copy (SSC) region of 18,384 bp, and a pair of inverted repeats (IRs) of 26,278 bp. The chloroplast genome encoded 113 unique genes, comprising 80 protein-coding genes, 29 tRNA genes, and four rRNA genes. Phylogenetic analysis inferred from the complete chloroplast genome indicated that Chrysoglossum was closely related to Collabium Blume. This study provides genomic resources helpful for further phylogenetic and biodiversity research on Chrysoglossum.

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.

The complete chloroplast genome sequences of four Liparis species (Orchidaceae) and phylogenetic implications

Liparis Richard (Malaxideae, Epidendroideae) is a large and diverse genus of the family Orchidaceae, the taxonomy of which is complicated and controversial. In this study, we sequenced, assembled and analyzed four complete chloroplast genomes of Liparis species including L. kumokiri, L. makinoana, L. pauliana, and L. viridiflora, and evaluated their phylogenetic relationships with related species for the first time. These four chloroplast genomes (size range 153,095 to 158,239 bp) possess typical quadripartite structures that consist of a large single copy (LSC, 83,533-86,752 bp), a small single copy (SSC, 17,938-18,156 bp) and a pair of inverted repeats (IRs, 26,421-26,933 bp). The genomes contain 133 genes, including 87 protein coding genes, 38 tRNAs and 8 rRNA genes. The genome arrangements, gene contents, gene order, long repeats and simple sequence repeats were similar with small differences observed among these four chloroplast genomes. Five highly variable regions including ycf1, ndhA, ndhF, trnQ and trnK were identified from the comparative analysis with other nine related Liparis species, which had the potential to be used as DNA markers for species identification and phylogenetic studies of Liparis species. Phylogenetic analysis based on the complete chloroplast genome sequences strongly supported the polyphyly of Liparis and its further division into three branches. These results provided valuable information to illustrate the complicated taxonomy, phylogeny and evolution process of the Liparis genus.

Comparative and phylogenetic analyses of nine complete chloroplast genomes of Orchidaceae

The orchid family has 200,000 species and 700 genera, and it is found worldwide in the tropics and subtropics. In China, there are 1247 species and subspecies of orchids belonging to the Orchidaceae family. Orchidaceae is one of the most diverse plant families in the world, known for their lush look, remarkable ecological tolerance, and capability for reproduction. It has significant decorative and therapeutic value. In terms of evolution, the orchid family is one of the more complicated groups, but up until now, little has been known about its affinities. This study examined the properties of 19 chloroplast (cp) genomes, of which 11 had previously been published and nine had only recently been revealed. Following that, topics such as analysis of selection pressure, codon usage, amino acid frequencies, repeated sequences, and reverse repeat contraction and expansion are covered. The Orchidaceae share similar cp chromosomal characteristics, and we have conducted a preliminary analysis of their evolutionary connections. The cp genome of this family has a typical tepartite structure and a high degree of consistency across species. Platanthera urceolata with more tandem repeats of the cp genome. Similar cp chromosomal traits can be seen in the orchidaceae. Galearis roborowskyi, Neottianthe cucullata, Neottianthe monophylla, Platanthera urceolata and Ponerorchis compacta are the closest cousins, according to phylogenetic study.

Biotechnological approaches for in vitro propagation, conservation and secondary metabolites production in Bulbophyllum, an endangered orchid genus: a review

Bulbophyllum represents the largest genus in the family Orchidaceae. The orchid species of this genus are widely used in the traditional medicine systems in different Asian countries such as China, India, Indonesia and Thailand. Studies on the secondary metabolites of Bulbophyllum have revealed the presence of important phytochemicals such as phenols, flavonoids, alkaloids, tannins, triterpenoids, sesquiterpenoids, steroids and glycosides. Some species of Bulbophyllum are reported to be of horticultural importance for their unique flowers. Habitat destruction and unsustainable utilization of different species of Bulbophyllum have led to a decline in the natural populations. The present review provides insights into the phytochemistry and ethnomedicinal uses of different species of Bulbophyllum, and highlights the biotechnological approaches developed for its conservation and sustainable utilization. Overall, the details provided in the present review can potentially be used for genome editing and biotechnological advances to develop plants with improved traits, which will be essential for the judicious utilization of the Bulbophyllum species so as to conserve and save the populations in the wild.

Characteristics and Comparative Analysis of Seven Complete Plastomes of Trichoglottis s.l. (Aeridinae, Orchidaceae)

Trichoglottis exhibits a range of rich variations in colors and shapes of flower and is a valuable ornamental orchid genus. The genus Trichoglottis has been expanded by the inclusion of Staurochilus, but this Trichoglottis sensu lato (s.l.) was recovered as a non-monophyletic genus based on molecular sequences from one or a few DNA regions. Here, we present phylogenomic data sets, incorporating complete plastome sequences from seven species (including five species sequenced in this study) of Trichoglottis s.l. (including two species formerly treated as Staurochilus), to compare plastome structure and to reconstruct the phylogenetic relationships of this genus. The seven plastomes possessed the typical quadripartite structure of angiosperms and ranged from 149,402 bp to 149,841 bp with a GC content of 36.6-36.7%. These plastomes contain 120 genes, which comprise 74 protein-coding genes, 38 tRNA genes, and 8 rRNA genes, all ndh genes were pseudogenized or lost. A total of 98 (T. philippinensis) to 134 (T. ionosma) SSRs and 33 (T. subviolacea) to 46 (T. ionosma) long repeats were detected. The consistent and robust phylogenetic relationships of Trichoglottis were established using a total of 25 plastid genomes from the Aeridinae subtribe. The genus Trichoglottis s.l. was strongly supported as a monophyletic group, and two species formerly treated as Staurochilus were revealed as successively basal lineages. In addition, five mutational hotspots (trnNGUU-rpl32, trnLUAA, trnSGCU-trnGUCC, rbcL-accD, and trnTGGU-psbD) were identified based on the ranking of PI values. Our research indicates that plastome data is a valuable source for molecular identification and evolutionary studies of Trichoglottis and its related genera.

Characterization and Comparative Analysis of the Complete Plastomes of Five Epidendrum (Epidendreae, Orchidaceae) Species

Epidendrum, one of the three largest genera of Orchidaceae, exhibits significant horticultural and ornamental value and serves as an important research model in conservation, ecology, and evolutionary biology. Given the ambiguous identification of germplasm and complex evolutionary relationships within the genus, the complete plastome of this genus (including five species) were firstly sequenced and assembled to explore their characterizations. The plastomes exhibited a typical quadripartite structure. The lengths of the plastomes ranged from 147,902 bp to 150,986 bp, with a GC content of 37.16% to 37.33%. Gene annotation revealed the presence of 78-82 protein-coding genes, 38 tRNAs, and 8 rRNAs. A total of 25-38 long repeats and 130-149 SSRs were detected. Analysis of relative synonymous codon usage (RSCU) indicated that leucine (Leu) was the most and cysteine (Cys) was the least. The consistent and robust phylogenetic relationships of Epidendrum and its closely related taxa were established using a total of 43 plastid genomes from the tribe Epidendreae. The genus Epidendrum was supported as a monophyletic group and as a sister to Cattleya. Meanwhile, four mutational hotspots (trnCGCA-petN, trnDGUC-trnYGUA, trnSGCU-trnGUCC, and rpl32-trnLUAG) were identified for further phylogenetic studies. Our analysis demonstrates the promising utility of plastomes in inferring the phylogenetic relationships of Epidendrum.

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.

Ancient Horizontal Gene Transfers from Plastome to Mitogenome of a Nonphotosynthetic Orchid, Gastrodia pubilabiata (Epidendroideae, Orchidaceae)

Gastrodia pubilabiata is a nonphotosynthetic and mycoheterotrophic orchid belonging to subfamily Epidendroideae. Compared to other typical angiosperm species, the plastome of G. pubilabiata is dramatically reduced in size to only 30,698 base pairs (bp). This reduction has led to the loss of most photosynthesis-related genes and some housekeeping genes in the plastome, which now only contains 19 protein coding genes, three tRNAs, and three rRNAs. In contrast, the typical orchid species contains 79 protein coding genes, 30 tRNAs, and four rRNAs. This study decoded the entire mitogenome of G. pubilabiata, which consisted of 44 contigs with a total length of 867,349 bp. Its mitogenome contained 38 protein coding genes, nine tRNAs, and three rRNAs. The gene content of G. pubilabiata mitogenome is similar to the typical plant mitogenomes even though the mitogenome size is twice as large as the typical ones. To determine possible gene transfer events between the plastome and the mitogenome individual BLASTN searches were conducted, using all available orchid plastome sequences and flowering plant mitogenome sequences. Plastid rRNA fragments were found at a high frequency in the mitogenome. Seven plastid protein coding gene fragments (ndhC, ndhJ, ndhK, psaA, psbF, rpoB, and rps4) were also identified in the mitogenome of G. pubilabiata. Phylogenetic trees using these seven plastid protein coding gene fragments suggested that horizontal gene transfer (HGT) from plastome to mitogenome occurred before losses of photosynthesis related genes, leading to the lineage of G. pubilabiata. Compared to species phylogeny of the lineage of orchid, it was estimated that HGT might have occurred approximately 30 million years ago.

Extreme plastomes in holoparasitic Balanophoraceae are not the norm

BACKGROUND

Balanophoraceae plastomes are known for their highly condensed and re-arranged nature alongside the most extreme nucleotide compositional bias known to date, culminating in two independent reconfigurations of their genetic code. Currently, a large portion of the Balanophoraceae diversity remains unexplored, hindering, among others, evolutionary pattern recognition. Here, we explored newly sequenced plastomes of Sarcophyte sanguinea and Thonningia sanguinea. The reconstructed plastomes were analyzed using various methods of comparative genomics based on a representative taxon sampling.

RESULTS

Sarcophyte, recovered sister to the other sampled Balanophoraceae s. str., has plastomes up to 50% larger than those currently published. Its gene set contains five genes lost in any other species, including matK. Five cis-spliced introns are maintained. In contrast, the Thonningia plastome is similarly reduced to published Balanophoraceae and retains only a single cis-spliced intron. Its protein-coding genes show a more biased codon usage compared to Sarcophyte, with an accumulation of in-frame TAG stop codons. Structural plastome comparison revealed multiple, previously unknown, structural rearrangements within Balanophoraceae.

CONCLUSIONS

For the "minimal plastomes" of Thonningia, we propose a genetic code change identical to sister genus Balanophora. Sarcophyte however differs drastically from our current understanding on Balanophoraceae plastomes. With a less-extreme nucleotide composition, there is no evidence for an altered genetic code. Using comparative genomics, we identified a hotspot for plastome reconfiguration in Balanophoraceae. Based on previously published and newly identified structural reconfigurations, we propose an updated model of evolutionary plastome trajectories for Balanophoraceae, illustrating a much greater plastome diversity than previously known.

Habitat-related plastome evolution in the mycoheterotrophic Neottia listeroides complex (Orchidaceae, Neottieae)

BACKGROUND

Mycoheterotrophs, acquiring organic carbon and other nutrients from mycorrhizal fungi, have evolved repeatedly with substantial plastid genome (plastome) variations. To date, the fine-scale evolution of mycoheterotrophic plastomes at the intraspecific level is not well-characterized. A few studies have revealed unexpected plastome divergence among species complex members, possibly driven by various biotic/abiotic factors. To illustrate evolutionary mechanisms underlying such divergence, we analyzed plastome features and molecular evolution of 15 plastomes of Neottia listeroides complex from different forest habitats.

RESULTS

These 15 samples of Neottia listeroides complex split into three clades according to their habitats approximately 6 million years ago: Pine Clade, including ten samples from pine-broadleaf mixed forests, Fir Clade, including four samples from alpine fir forests and Fir-willow Clade with one sample. Compared with those of Pine Clade members, plastomes of Fir Clade members show smaller size and higher substitution rates. Plastome size, substitution rates, loss and retention of plastid-encoded genes are clade-specific. We propose to recognized six species in N. listeroides complex and slightly modify the path of plastome degradation.

CONCLUSIONS

Our results provide insight into the evolutionary dynamics and discrepancy of closely related mycoheterotrophic orchid lineages at a high phylogenetic resolution.

Comprehensive phylogenetic analyses of Orchidaceae using nuclear genes and evolutionary insights into epiphytism

Orchidaceae (with >28,000 orchid species) are one of the two largest plant families, with economically and ecologically important species, and occupy global and diverse niches with primary distribution in rainforests. Among orchids, 70% grow on other plants as epiphytes; epiphytes contribute up to ~50% of the plant diversity in rainforests and provide food and shelter for diverse animals and microbes, thereby contributing to the health of these ecosystems. Orchids account for over two-thirds of vascular epiphytes and provide an excellent model for studying evolution of epiphytism. Extensive phylogenetic studies of Orchidaceae and subgroups have ;been crucial for understanding relationships among many orchid lineages, although some uncertainties remain. For example, in the largest subfamily Epidendroideae with nearly all epiphytic orchids, relationships among some tribes and many subtribes are still controversial, hampering evolutionary analyses of epiphytism. Here we obtained 1,450 low-copy nuclear genes from 610 orchid species, including 431 with newly generated transcriptomes, and used them for the reconstruction of robust Orchidaceae phylogenetic trees with highly supported placements of tribes and subtribes. We also provide generally well-supported phylogenetic placements of 131 genera and 437 species that were not sampled by previous plastid and nuclear phylogenomic studies. Molecular clock analyses estimated the Orchidaceae origin at ~132 million years ago (Ma) and divergences of most subtribes from 52 to 29 Ma. Character reconstruction supports at least 14 parallel origins of epiphytism; one such origin was placed at the most recent common ancestor of ~95% of epiphytic orchids and linked to modern rainforests. Ten occurrences of rapid increase in the diversification rate were detected within Epidendroideae near and after the K-Pg boundary, contributing to ~80% of the Orchidaceae diversity. This study provides a robust and the largest family-wide Orchidaceae nuclear phylogenetic tree thus far and new insights into the evolution of epiphytism in vascular plants.

Evolutionary differences in gene loss and pseudogenization among mycoheterotrophic orchids in the tribe Vanilleae (subfamily Vanilloideae)

INTRODUCTION

Galeola lindleyana is a mycoheterotrophic orchid belonging to the tribe Vanilleae within the subfamily Vanilloideae.

METHODS

In this study, the G. lindleyana plastome was assembled and annotated, and compared with other Vanilleae orchids, revealing the evolutionary variations between the photoautotrophic and mycoheterotrophic plastomes.

RESULTS

The G. lindleyana plastome was found to include 32 protein-coding genes, 16 tRNA genes and four ribosomal RNA genes, including 11 pseudogenes. Almost all of the genes encoding photosynthesis have been lost physically or functionally, with the exception of six genes encoding ATP synthase and psaJ in photosystem I. The length of the G. lindleyana plastome has decreased to 100,749 bp, while still retaining its typical quadripartite structure. Compared with the photoautotrophic Vanilloideae plastomes, the inverted repeat (IR) regions and the large single copy (LSC) region of the mycoheterotrophic orchid's plastome have contracted, while the small single copy (SSC) region has expanded significantly. Moreover, the difference in length between the two ndhB genes was found to be 682 bp, with one of them spanning the IRb/SSC boundary. The Vanilloideae plastomes were varied in their structural organization, gene arrangement, and gene content. Even the Cyrtosia septentrionalis plastome which was found to be closest in length to the G. lindleyana plastome, differed in terms of its gene arrangement and gene content. In the LSC region, the psbA, psbK, atpA and psaB retained in the G. lindleyana plastome were missing in the C. septentrionalis plastome, while, the matK, rps16, and atpF were incomplete in the C. septentrionalis plastome, yet still complete in that of the G. lindleyana. Lastly, compared with the G. lindleyana plastome, a 15 kb region located in the SSC area between ndhB-rrn16S was found to be inverted in the C. septentrionalis plastome. These changes in gene content, gene arrangment and gene structure shed light on the polyphyletic evolution of photoautotrophic orchid plastomes to mycoheterotrophic orchid plastomes.

DISCUSSION

Thus, this study's decoding of the mycoheterotrophic G. lindleyana plastome provides valuable resource data for future research and conservation of endangered orchids.

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

Identification of Laportea bulbifera using the complete chloroplast genome as a potentially effective super-barcode

Laportea bulbifera, a Miao medicine grown in karst areas, has exerted a unique curative effect on skin itching in the elderly, with an annual sales of > 100 million Yuan. Owing to the shortage of resources and large morphological variations in L. bulbifera, it is difficult to identify the species correctly using only traditional methods, which seriously affects the safety of drug usage for patients. This study obtained the complete high-quality L. bulbifera chloroplast (cp) genome, using second- and third-generation high-throughput sequencing. The cp genome was 149,911 bp in length, with a typical quadripartite structure. A total of 127 genes were annotated, including 83 protein-coding genes, 36 tRNA genes, and 8 rRNA genes. There was an inverted small single copy (SSC) structure in the L. bulbifera cp genome, one large-scale rearrangement of ~ 39 kb excised in the SSC and IR regions. The complete cp genome sequence is used as a potentially effective super-barcode and the highly variable regions (ycf1, matK, and ndhD) can be used as potentially specific barcodes to accurately distinguish L. bulbifera from counterfeits and closely related species. This study is important for the identification of L. bulbifera and lays a theoretical foundation for elucidating the phylogenetic relationship of the species.

The ecological adaptation of the unparalleled plastome character evolution in slipper orchids

Plastomes may have undergone adaptive evolution in the process of plant adaptation to diverse environments, whereby species may differ in plastome characters. Cypripedioideae successfully colonized distinct environments and could be an ideal group for studying the interspecific variation and adaptive evolution of plastomes. Comparative study of plastomes, ancestral state reconstruction, phylogenetic-based analysis, ecological niche modelling, and selective pressure analysis were conducted to reveal the evolutionary patterns of plastomes in Cypripedioideae and their relationship with environmental factors. The plastomes of the three evolved genera had reduced plastome size, increased GC content, and compacted gene content compared to the basal group. Variations in plastome size and GC content are proved to have clear relationships with climate regions. Furthermore, ecological niche modelling revealed that temperature and water factors are important climatic factors contributing to the distributional difference which is directly correlated with the climate regions. The temperature-sensitive genes ndh genes, infA, and rpl20 were found to be either lost/pseudogenized or under positive selection in the evolved groups. Unparalleled plastome character variations were discovered in slipper orchids. Our study indicates that variations in plastome characters have adaptive consequences and that temperature and water factors are important climatic factors that affect plastome evolution. This research highlights the expectation that plants can facilitate adaptation to different environmental conditions with the changes in plastome and has added critical insight for understanding the process of plastome evolution in plants.

The complete chloroplast genome of Dendrobium moschatum (Buch.-Ham.) Sw. 1805 (Orchidaceae)

The morphological characteristic of Dendrobium moschatum (Buch.-Ham.) Sw. 1805 is very distinctive among Dendrobium Sw. 1799, and it has high medicinal and ornamental values. Here, we reported the first complete chloroplast genome of D. moschatum. The complete genome of D. moschatum was 159,701 bp in length with 130 genes, including 38 tRNA, 8 rRNA, and 84 protein-coding genes. Phylogenetic analysis showed that D. moschatum was strongly allied with D. denneanum Kerr. 1933.

Phylogenomics and plastome evolution of a Brazilian mycoheterotrophic orchid, Pogoniopsis schenckii

PREMISE

Pogoniopsis likely represents an independent photosynthesis loss in orchids. We use phylogenomic data to better identify the phylogenetic placement of this fully mycoheterotrophic taxon, and investigate its molecular evolution.

METHODS

We performed likelihood analysis of plastid and mitochondrial phylogenomic data to localize the position of Pogoniopsis schenckii in orchid phylogeny, and investigated the evolution of its plastid genome.

RESULTS

All analyses place Pogoniopsis in subfamily Epidendroideae, with strongest support from mitochondrial data, which also place it near tribe Sobralieae with moderately strong support. Extreme rate elevation in Pogoniopsis plastid genes broadly depresses branch support; in contrast, mitochondrial genes are only mildly rate elevated and display very modest and localized reductions in bootstrap support. Despite considerable genome reduction, including loss of photosynthesis genes and multiple translation apparatus genes, gene order in Pogoniopsis plastomes is identical to related autotrophs, apart from moderately shifted inverted repeat (IR) boundaries. All cis-spliced introns have been lost in retained genes. Two plastid genes (accD, rpl2) show significant strengthening of purifying selection. A retained plastid tRNA gene (trnE-UUC) of Pogoniopsis lacks an anticodon; we predict that it no longer functions in translation but retains a secondary role in heme biosynthesis.

CONCLUSIONS

Slowly evolving mitochondrial genes clarify the placement of Pogoniopsis in orchid phylogeny, a strong contrast with analysis of rate-elevated plastome data. We documented the effects of the novel loss of photosynthesis: for example, despite massive gene loss, its plastome is fully colinear with other orchids, and it displays only moderate shifts in selective pressure in retained genes.

Comparative and phylogenetic analysis of complete chloroplast genomes from five Artemisia species

Artemisia Linn. is a large genus within the family Asteraceae that includes several important medicinal plants. Because of their similar morphology and chemical composition, traditional identification methods often fail to distinguish them. Therefore, developing an effective identification method for Artemisia species is an urgent requirement. In this study, we analyzed 15 chloroplast (cp) genomes, including 12 newly sequenced genomes, from 5 Artemisia species. The cp genomes from the five Artemisia species had a typical quadripartite structure and were highly conserved across species. They had varying lengths of 151,132-151,178 bp, and their gene content and codon preferences were similar. Mutation hotspot analysis identified four highly variable regions, which can potentially be used as molecular markers to identify Artemisia species. Phylogenetic analysis showed that the five Artemisia species investigated in this study were sister branches to each other, and individuals of each species formed a monophyletic clade. This study shows that the cp genome can provide distinguishing features to help identify closely related Artemisia species and has the potential to serve as a universal super barcode for plant identification.

Evolution patterns of NBS genes in the genus Dendrobium and NBS-LRR gene expression in D. officinale by salicylic acid treatment

BACKGROUND

Dendrobium officinale Kimura et Migo, which contains rich polysaccharides, flavonoids and alkaloids, is a Traditional Chinese Medicine (TCM) with important economic benefits, while various pathogens have brought huge losses to its industrialization. NBS gene family is the largest class of plant disease resistance (R) genes, proteins of which are widely distributed in the upstream and downstream of the plant immune systems and are responsible for receiving infection signals and regulating gene expression respectively. It is of great significance for the subsequent disease resistance breeding of D. officinale to identify NBS genes by using the newly published high-quality chromosome-level D. officinale genome.

RESULTS

In this study, a total of 655 NBS genes were uncovered from the genomes of D. officinale, D. nobile, D. chrysotoxum, V. planifolia, A. shenzhenica, P. equestris and A. thaliana. The phylogenetic results of CNL-type protein sequences showed that orchid NBS-LRR genes have significantly degenerated on branches a and b. The Dendrobium NBS gene homology analysis showed that the Dendrobium NBS genes have two obvious characteristics: type changing and NB-ARC domain degeneration. Because the NBS-LRR genes have both NB-ARC and LRR domains, 22 D. officinale NBS-LRR genes were used for subsequent analyses, such as gene structures, conserved motifs, cis-elements and functional annotation analyses. All these results suggested that D. officinale NBS-LRR genes take part in the ETI system, plant hormone signal transduction pathway and Ras signaling pathway. Finally, there were 1,677 DEGs identified from the salicylic acid (SA) treatment transcriptome data of D. officinale. Among them, six NBS-LRR genes (Dof013264, Dof020566, Dof019188, Dof019191, Dof020138 and Dof020707) were significantly up-regulated. However, only Dof020138 was closely related to other pathways from the results of WGCNA, such as pathogen identification pathways, MAPK signaling pathways, plant hormone signal transduction pathways, biosynthetic pathways and energy metabolism pathways.

CONCLUSION

Our results revealed that the NBS gene degenerations are common in the genus Dendrobium, which is the main reason for the diversity of NBS genes, and the NBS-LRR genes generally take part in D. officinale ETI system and signal transduction pathways. In addition, the D. officinale NBS-LRR gene Dof020138, which may have an important breeding value, is indirectly activated by SA in the ETI system.

Phylogenomics, plastome degradation and mycoheterotrophy evolution of Neottieae (Orchidaceae), with emphasis on the systematic position and Loess Plateau-Changbai Mountains disjunction of Diplandrorchis

BACKGROUND

Mycoheterotrophy is a unique survival strategy adapted to dense forests and has attracted biologists' attention for centuries. However, its evolutionary origin and related plastome degradation are poorly understood. The tribe Neottieae contains various nutrition types, i.e., autotrophy, mixotrophy, and mycoheterotrophy. Here, we present a comprehensive phylogenetic analysis of the tribe based on plastome and nuclear ITS data. We inferred the evolutionary shift of nutrition types, constructed the patterns of plastome degradation, and estimated divergence times and ancestral ranges. We also used an integration of molecular dating and ecological niche modeling methods to investigate the disjunction between the Loess Plateau and Changbai Mountains in Diplandrorchis, a mycoheterotrophic genus endemic to China that was included in a molecular phylogenetic study for the first time.

RESULTS

Diplandrorchis was imbedded within Neottia and formed a clade with four mycoheterotrophic species. Autotrophy is the ancestral state in Neottieae, mixotrophy independently originated at least five times, and three shifts from mixotrophy to mycoheterotrophy independently occurred. The five mixotrophic lineages possess all plastid genes or lost partial/all ndh genes, whereas each of the three mycoheterotroph lineages has a highly reduced plastome: one lost part of its ndh genes and a few photosynthesis-related genes, and the other two lost almost all ndh, photosynthesis-related, rpo, and atp genes. These three mycoheterotrophic lineages originated at about 26.40 Ma, 25.84 Ma, and 9.22 Ma, respectively. Diplandrorchis had presumably a wide range in the Pliocene and migrated southward in the Pleistocene.

CONCLUSIONS

The Pleistocene climatic fluctuations and the resultant migration resulted in the Loess Plateau-Changbai Mountains disjunction of Diplandrorchis. In the evolution of mycoheterotrophic lineages, the loss of plastid-encoded genes and plastome degradation are staged and irreversible, constraining mycoheterotrophs to inhabit understories with low light levels. Accordingly, the rise of local forests might have promoted the origin of conditions in which mycoheterotrophy is advantageous.

Danxiaorchismangdangshanensis (Orchidaceae, Epidendroideae), a new species from central Fujian Province based on morphological and genomic data

Danxiaorchismangdangshanensis, a new mycoheterotrophic species from Fujian Province, China, is described and illustrated. The new species is morphologically similar to D.singchiana, but its callus of labellum is a less distinctive Y-shape with three auricles on the apex, four pollinia that are narrowly elliptic in shape and equal in size, and it lacks fine roots. The plastome of D.mangdangshanensis is highly degraded. Phylogenetic analyses distinguished D.mangdangshanensis from its congeners, D.singchiana and D.yangii, with strong support based on nrITS + matK and plastomes, respectively.

Characterization of the plastome of Physaliscordata and comparative analysis of eight species of Physalis sensu stricto

In this study, we sequenced, assembled, and annotated the plastome of Physaliscordata Mill. and compared it with seven species of the genus Physalis sensu stricto. Sequencing, annotating, and comparing plastomes allow us to understand the evolutionary mechanisms associated with physiological functions, select possible molecular markers, and identify the types of selection that have acted in different regions of the genome. The plastome of P.cordata is 157,000 bp long and presents the typical quadripartite structure with a large single-copy (LSC) region of 87,267 bp and a small single-copy (SSC) region of 18,501 bp, which are separated by two inverted repeat (IRs) regions of 25,616 bp each. These values are similar to those found in the other species, except for P.angulata L. and P.pruinosa L., which presented an expansion of the LSC region and a contraction of the IR regions. The plastome in all Physalis species studied shows variation in the boundary of the regions with three distinct types, the percentage of the sequence identity between coding and non-coding regions, and the number of repetitive regions and microsatellites. Four genes and 10 intergenic regions show promise as molecular markers and eight genes were under positive selection. The maximum likelihood analysis showed that the plastome is a good source of information for phylogenetic inference in the genus, given the high support values and absence of polytomies. In the Physalis plastomes analyzed here, the differences found, the positive selection of genes, and the phylogenetic relationships do not show trends that correspond to the biological or ecological characteristics of the species studied.

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.

Plastomes of Bletilla (Orchidaceae) and Phylogenetic Implications

The genus Bletilla is a small genus of only five species distributed across Asia, including B. chartacea, B. foliosa, B. formosana, B. ochracea and B. striata, which is of great medicinal importance. Furthermore, this genus is a member of the key tribe Arethuseae (Orchidaceae), harboring an extremely complicated taxonomic history. Recently, the monophyletic status of Bletilla has been challenged, and the phylogenetic relationships within this genus are still unclear. The plastome, which is rich in both sequence and structural variation, has emerged as a powerful tool for understanding plant evolution. Along with four new plastomes, this work is committed to exploring plastomic markers to elucidate the phylogeny of Bletilla. Our results reveal considerable plastomic differences between B. sinensis and the other three taxa in many aspects. Most importantly, the specific features of the IR junction patterns, novel pttRNA structures and codon aversion motifs can serve as useful molecular markers for Bletilla phylogeny. Moreover, based on maximum likelihood and Bayesian inference methods, our phylogenetic analyses based on two datasets of Arethuseae strongly imply that Bletilla is non-monophyletic. Accordingly, our findings from this study provide novel potential markers for species identification, and shed light on the evolution of Bletilla and Arethuseae.

Plastome structure of 8 Calanthe s.l. species (Orchidaceae): comparative genomics, phylogenetic analysis

Background

Calanthe (Epidendroideae, Orchidaceae) is a pantropical genus distributed in Asia and Africa. Its species are of great importance in terms of economic, ornamental and medicinal values. However, due to limited and confusing delimitation characters, the taxonomy of the Calanthe alliance (Calanthe, Cephalantheropsis, and Phaius) has not been sufficiently resolved. Additionally, the limited genomic information has shown incongruences in its systematics and phylogeny. In this study, we used illumina platform sequencing, performed a de novo assembly, and did a comparative analysis of 8 Calanthe group species' plastomes: 6 Calanthe and 2 Phaius species. Phylogenetic analyses were used to reconstruct the relationships of the species as well as with other species of the family Orchidaceae.

Results

The complete plastomes of the Calanthe group species have a quadripartite structure with varied sizes ranging between 150,105bp-158,714bp, including a large single-copy region (LSC; 83,364bp- 87,450bp), a small single-copy region (SSC; 16,297bp -18,586bp), and a pair of inverted repeat regions (IRs; 25,222bp - 26,430bp). The overall GC content of these plastomes ranged between 36.6-36.9%. These plastomes encoded 131-134 differential genes, which included 85-88 protein-coding genes, 37-38 tRNA genes, and 8 rRNA genes. Comparative analysis showed no significant variations in terms of their sequences, gene content, gene order, sequence repeats and the GC content hence highly conserved. However, some genes were lost in C. delavayi (P. delavayi), including ndhC, ndhF, and ndhK genes. Compared to the coding regions, the non-coding regions had more sequence repeats hence important for species DNA barcoding. Phylogenetic analysis revealed a paraphyletic relationship in the Calanthe group, and confirmed the position of Phaius delavayi in the genus Calanthe as opposed to its previous placement in Phaius.

Conclusion

This study provides a report on the complete plastomes of 6 Calanthe and 2 Phaius species and elucidates the structural characteristics of the plastomes. It also highlights the power of plastome data to resolve phylogenetic relationships and clarifies taxonomic disputes among closely related species to improve our understanding of their systematics and evolution. Furthermore, it also provides valuable genetic resources and a basis for studying evolutionary relationships and population genetics among orchid species.

Plastome structure of 8 Calanthe s.l. species (Orchidaceae): comparative genomics, phylogenetic analysis

BACKGROUND

Calanthe (Epidendroideae, Orchidaceae) is a pantropical genus distributed in Asia and Africa. Its species are of great importance in terms of economic, ornamental and medicinal values. However, due to limited and confusing delimitation characters, the taxonomy of the Calanthe alliance (Calanthe, Cephalantheropsis, and Phaius) has not been sufficiently resolved. Additionally, the limited genomic information has shown incongruences in its systematics and phylogeny. In this study, we used illumina platform sequencing, performed a de novo assembly, and did a comparative analysis of 8 Calanthe group species' plastomes: 6 Calanthe and 2 Phaius species. Phylogenetic analyses were used to reconstruct the relationships of the species as well as with other species of the family Orchidaceae.

RESULTS

The complete plastomes of the Calanthe group species have a quadripartite structure with varied sizes ranging between 150,105bp-158,714bp, including a large single-copy region (LSC; 83,364bp- 87,450bp), a small single-copy region (SSC; 16,297bp -18,586bp), and a pair of inverted repeat regions (IRs; 25,222bp - 26,430bp). The overall GC content of these plastomes ranged between 36.6-36.9%. These plastomes encoded 131-134 differential genes, which included 85-88 protein-coding genes, 37-38 tRNA genes, and 8 rRNA genes. Comparative analysis showed no significant variations in terms of their sequences, gene content, gene order, sequence repeats and the GC content hence highly conserved. However, some genes were lost in C. delavayi (P. delavayi), including ndhC, ndhF, and ndhK genes. Compared to the coding regions, the non-coding regions had more sequence repeats hence important for species DNA barcoding. Phylogenetic analysis revealed a paraphyletic relationship in the Calanthe group, and confirmed the position of Phaius delavayi in the genus Calanthe as opposed to its previous placement in Phaius.

CONCLUSION

This study provides a report on the complete plastomes of 6 Calanthe and 2 Phaius species and elucidates the structural characteristics of the plastomes. It also highlights the power of plastome data to resolve phylogenetic relationships and clarifies taxonomic disputes among closely related species to improve our understanding of their systematics and evolution. Furthermore, it also provides valuable genetic resources and a basis for studying evolutionary relationships and population genetics among orchid species.

Uncovering haplotype diversity in cultivated Mexican vanilla species

PREMISE

Although vanilla is one of the best-known spices, there is only a limited understanding of its biology and genetics within Mexico, where its cultivation originated and where phenotypic variability is high. This study aims to augment our understanding of vanilla's genetic resources by assessing species delimitation and genetic, geographic, and climatic variability within Mexican cultivated vanilla.

METHODS

Using nuclear and plastid DNA sequence data from 58 Mexican samples collected from three regions and 133 ex situ accessions, we assessed species monophyly using phylogenetic analyses and genetic distances. Intraspecific genetic variation was summarized through the identification of haplotypes. Within the primarily cultivated species, Vanilla planifolia, haplotype relationships were further verified using plastome and rRNA gene sequences. Climatic niche and haplotype composition were assessed across the landscape.

RESULTS

Three species (Vanilla planifolia, V. pompona, and V. insignis) and 13 haplotypes were identified among Mexican vanilla. Within V. planifolia haplotypes, hard phylogenetic incongruences between plastid and nuclear sequences suggest past hybridization events. Eight haplotypes consisted exclusively of Mexican samples. The dominant V. planifolia haplotype occurred throughout all three regions as well as outside of its country of origin. Haplotype richness was found to be highest in regions around Papantla and Chinantla.

CONCLUSIONS

Long histories of regional cultivation support the consideration of endemic haplotypes as landraces shaped by adaptation to local conditions and/or hybridization. Results may aid further genomic investigations of vanilla's genetic resources and ultimately support the preservation of genetic diversity within the economically important crop.

Orchid Phylotranscriptomics: The Prospects of Repurposing Multi-Tissue Transcriptomes for Phylogenetic Analysis and Beyond

The Orchidaceae is rivaled only by the Asteraceae as the largest plant family, with the estimated number of species exceeding 25,000 and encompassing more than 700 genera. To gain insights into the mechanisms driving species diversity across both global and local scales, well-supported phylogenies targeting different taxonomic groups and/or geographical regions will be crucial. High-throughput sequencing technologies have revolutionized the field of molecular phylogenetics by simplifying the process of obtaining genome-scale sequence data. Consequently, there has been an explosive growth of such data in public repositories. Here we took advantage of this unprecedented access to transcriptome data from predominantly non-phylogenetic studies to assess if it can be repurposed to gain rapid and accurate phylogenetic insights across the orchids. Exhaustive searches revealed transcriptomic data for more than 100 orchid species spanning 5 subfamilies, 13 tribes, 21 subtribes, and 50 genera that were amendable for exploratory phylotranscriptomic analysis. Next, we performed re-assembly of the transcriptomes before strategic selection of the final samples based on a gene completeness evaluation. Drawing on these data, we report phylogenetic analyses at both deep and shallow evolutionary scales via maximum likelihood and shortcut coalescent species tree methods. In this perspective, we discuss some key outcomes of this study and conclude by highlighting other complementary, albeit rarely explored, insights beyond phylogenetic analysis that repurposed multi-tissue transcriptome can offer.

Plastid RNA editing reduction accompanied with genetic variations in Cymbidium, a genus with diverse lifestyle modes

Recent sequencing efforts have broadly uncovered the evolutionary trajectory of plastid genomes (plastomes) of flowering plants in diverse habitats, yet our knowledge of the evolution of plastid posttranscriptional modifications is limited. In this study, we generated 11 complete plastomes and performed ultra-deep transcriptome sequencing to investigate the co-evolution of plastid RNA editing and genetic variation in Cymbidium, a genus with diverse trophic lifestyles. Genome size and gene content is reduced in terrestrial and green mycoheterotrophic orchids relative to their epiphytic relatives. This could be partly due to extensive losses and pseudogenization of ndh genes for the plastid NADH dehydrogenase-like complex, but independent pseudogenization of ndh genes has also occurred in the epiphyte C. mannii, which was reported to use strong crassulacean acid metabolism photosynthesis. RNA editing sites are abundant but variable in number among Cymbidium plastomes. The nearly twofold variation in editing abundance is mainly due to extensive reduction of ancestral editing sites in ndh transcripts of terrestrial, mycoheterotrophic, and C. mannii plastomes. The co-occurrence of editing reduction and pseudogenization in ndh genes suggests functional constraints on editing machinery may be relaxed, leading to nonrandom loss of ancestral edited sites via reduced editing efficiency. This study represents the first systematic examination of RNA editing evolution linked to plastid genome variation in a single genus. We also propose an explanation for how genomic and posttranscriptional variations might be affected by lifestyle-associated ecological adaptation strategies in Cymbidium.

Infrageneric Plastid Genomes of Cotoneaster (Rosaceae): Implications for the Plastome Evolution and Origin of C. wilsonii on Ulleung Island

Cotoneaster is a taxonomically and ornamentally important genus in the family Rosaceae; however, phylogenetic relationships among its species are complicated owing to insufficient morphological diagnostic characteristics and hybridization associated with polyploidy and apomixis. In this study, we sequenced the complete plastomes of seven Cotoneaster species (C. dielsianus, C. hebephyllus, C. integerrimus, C. mongolicus, C. multiflorus, C. submultiflorus, and C. tenuipes) and included the available complete plastomes in a phylogenetic analysis to determine the origin of C. wilsonii, which is endemic to Ulleung Island, Korea. Furthermore, based on 15 representative lineages within the genus, we carried out the first comparative analysis of Cotoneaster plastid genomes to gain an insight into their molecular evolution. The plastomes were highly conserved, with sizes ranging from 159,595 bp (C. tenuipes) to 160,016 bp (C. hebephyllus), and had a GC content of 36.6%. The frequency of codon usage showed similar patterns among the 15 Cotoneaster species, and 24 of the 35 protein-coding genes were predicted to undergo RNA editing. Eight of the 76 common protein-coding genes, including ccsA, matK, ndhD, ndhF, ndhK, petA, rbcL, and rpl16, were positively selected, implying their potential roles in adaptation and speciation. Of the 35 protein-coding genes, 24 genes (15 photosynthesis-related, seven self-replications, and three others) were found to harbor RNA editing sites. Furthermore, several mutation hotspots were identified, including trnG-UCC/trnR-UCU/atpA and trnT-UGU/trnL-UAA. Maximum likelihood analysis based on 57 representative plastomes of Cotoneaster and two Heteromeles plastomes as outgroups revealed two major lineages within the genus, which roughly correspond to two subgenera, Chaenopetalum and Cotoneaster. The Ulleung Island endemic, C. wilsonii, shared its most recent common ancestor with two species, C. schantungensis and C. zabelii, suggesting its potential origin from geographically close members of the subgenus Cotoneaster, section Integerrimi.

Effects of geographic isolation on the Bulbophyllum chloroplast genomes

BACKGROUND

Because chloroplast (cp) genome has more conserved structures than nuclear genome and mitochondrial genome, it is a useful tool in estimating the phylogenetic relationships of plants. With a series of researches for cp genomes, there have been comprehensive understandings about the cp genome features. The genus Bulbophyllum widely distributed in Asia, South America, Australia and other places. Therefore, it is an excellent type genus for studying the effects of geographic isolation.

RESULTS

In this study, the cp genomes of nine Bulbophyllum orchids were newly sequenced and assembled using the next-generation sequencing technology. Based on 19 Asian (AN) and eight South American (SA) Bulbophyllum orchids, the cp genome features of AN clade and SA clade were compared. Comparative analysis showed that there were considerable differences in overall cp genome features between two clades in three aspects, including basic cp genome features, SSC/IR_B junctions (J_SBs) and mutational hotspots. The phylogenetic analysis and divergence time estimation results showed that the AN clade has diverged from the SA clade in the late Oligocene (21.50-30.12 mya). After estimating the occurrence rates of the insertions and deletions (InDels), we found that the change trends of cp genome structures between two clades were different under geographic isolation. Finally, we compared selective pressures on cp genes and found that long-term geographic isolation made AN and SA Bulbophyllum cp genes evolved variably.

CONCLUSION

The results revealed that the overall structural characteristics of Bulbophyllum cp genomes diverged during the long-term geographic isolation, and the crassulacean acid metabolism (CAM) pathway may play an important role in the Bulbophyllum species evolution.

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.

Ethnobotany, Phytochemistry, Biological Activities, and Health-Promoting Effects of the Genus Bulbophyllum

The genus Bulbophyllum is of scientific interest due to the phytochemical components and diverse biological activities found across species of the genus. Most Bulbophyllum species are epiphytic and located in habitats that range from subtropical dry forests to wet montane cloud forests. In many cultures, the genus Bulbophyllum has a religious, protective, ornamenting, cosmetic, and medicinal role. Detailed investigations into the molecular pharmacological mechanisms and numerous biological effects of Bulbophyllum spp. remain ambiguous. The review focuses on an in-depth discussion of studies containing data on phytochemistry and preclinical pharmacology. Thus, the purpose of this review was to summarize the therapeutic potential of Bulbophyllum spp. biocompounds. Data were collected from several scientific databases such as PubMed and ScienceDirect, other professional websites, and traditional medicine books to obtain the necessary information. Evidence from pharmacological studies has shown that various phytoconstituents in some Bulbophyllum species have different biological health-promoting activities such as antimicrobial, antifungal, antioxidant, anti-inflammatory, anticancer, and neuroprotective. No toxicological effects have been reported to date. Future clinical trials are needed for the clinical confirmation of biological activities proven in preclinical studies. Although orchid species are cultivated for ornamental purposes and have a wide traditional use, the novelty of this review is a summary of biological actions from preclinical studies, thus supporting ethnopharmacological data.

Complete pan-plastome sequences enable high resolution phylogenetic classification of sugar beet and closely related crop wild relatives

Background

As the major source of sugar in moderate climates, sugar-producing beets (Beta vulgaris subsp. vulgaris) have a high economic value. However, the low genetic diversity within cultivated beets requires introduction of new traits, for example to increase their tolerance and resistance attributes – traits that often reside in the crop wild relatives. For this, genetic information of wild beet relatives and their phylogenetic placements to each other are crucial. To answer this need, we sequenced and assembled the complete plastome sequences from a broad species spectrum across the beet genera Beta and Patellifolia, both embedded in the Betoideae (order Caryophyllales). This pan-plastome dataset was then used to determine the wild beet phylogeny in high-resolution.

Results

We sequenced the plastomes of 18 closely related accessions representing 11 species of the Betoideae subfamily and provided high-quality plastome assemblies which represent an important resource for further studies of beet wild relatives and the diverse plant order Caryophyllales. Their assembly sizes range from 149,723 bp (Beta vulgaris subsp. vulgaris) to 152,816 bp (Beta nana), with most variability in the intergenic sequences. Combining plastome-derived phylogenies with read-based treatments based on mitochondrial information, we were able to suggest a unified and highly confident phylogenetic placement of the investigated Betoideae species.

Our results show that the genus Beta can be divided into the two clearly separated sections Beta and Corollinae. Our analysis confirms the affiliation of B. nana with the other Corollinae species, and we argue against a separate placement in the Nanae section. Within the Patellifolia genus, the two diploid species Patellifolia procumbens and Patellifolia webbiana are, regarding the plastome sequences, genetically more similar to each other than to the tetraploid Patellifolia patellaris. Nevertheless, all three Patellifolia species are clearly separated.

Conclusion

In conclusion, our wild beet plastome assemblies represent a new resource to understand the molecular base of the beet germplasm. Despite large differences on the phenotypic level, our pan-plastome dataset is highly conserved. For the first time in beets, our whole plastome sequences overcome the low sequence variation in individual genes and provide the molecular backbone for highly resolved beet phylogenomics. Hence, our plastome sequencing strategy can also guide genomic approaches to unravel other closely related taxa.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08336-8.

Comparative Analysis of the Complete Chloroplast Genomes of Nine Paphiopedilum Species

Paphiopedilum is known as “lady’s or Venus” slipper orchids due to its prominent shoe-shaped labellum, with high ornamental value. Phylogenetic relationships among some species in Paphiopedilum genus cannot be effectively determined by morphological features alone or through the analysis of nuclear or chloroplast DNA fragments. In order to provide aid in understanding the evolutionary and phylogenetic relationship in Paphiopedilum at chloroplast (cp) genome-scale level, the complete cp genomes of six Paphiopedilum species were newly sequenced in this study, and three other published cp genome sequences of Paphiopedilum were included in the comparative analyses. The cp genomes of the six Paphiopedilum species ranged from 154,908 bp (P. hirsutissimum) to 161,300 bp (P. victoria-mariae) in size, all constituting four-part annular structures. Analyses of the nucleotide substitutions, insertions/deletions, and simple sequence repeats in the cp genomes were conducted. Ten highly variable regions that could serve as potential DNA barcodes or phylogenetic markers for this diverse genus were identified. Sequence variations in the non-coding regions were greater than that in the conserved protein-coding regions, as well as in the large single copy (LSC) and small single copy (SSC) regions than in the inverted repeat (IR) regions. Phylogenetic analysis revealed that all Paphiopedilum species clustered in one monophyletic clade in the Cypripedioideae subfamily and then subdivided into seven smaller branches corresponding to different subgenus or sections of the genus, with high bootstrap supports, indicate that cp genome sequencing can be an effective means in resolving the complex relationship in Paphiopedilum.

Ecological and phylogenetic constraints determine the stage of anthetic ovule development in orchids

PREMISE

Unlike most flowering plants, orchid flowers have under-developed ovules that complete development only after pollination. Classical studies reported variation in the stage in which ovule development is arrested, but the extent of this variation and its evolutionary and ecological significance are unclear.

METHODS

Here, we used light microscopy to observe ovule development at anthesis for 39 species not previously studied and surveyed the literature gaining information on 94 orchid species. Tropical and temperate members of all five orchid subfamilies as well as species with contrasting pollination strategies (rewarding versus deceptive) and life forms (epiphytic versus terrestrial) were represented. We analyzed the data using statistical comparisons and a phylogenetic generalized least square (PGLS) analysis.

RESULTS

Apostasioideae, the sister to the rest of the orchids, have mature ovules similar to other Asparagales, while under-differentiated ovules are present in the other subfamilies. Ovule developmental stages showed high variation even among closely related groups. Ovules were more developed in terrestrial than in epiphytic, in temperate than in tropical, and in rewarding than in deceptive pollination orchid species. This latter comparison was also significant in the PGLS analysis.

CONCLUSIONS

These results suggest that ovule developmental stage in orchids can be shaped by ecological factors, such as seasonality and pollination strategy, and can be selected for optimizing female reproductive investment.

Comparative plastomics of Amaryllidaceae: inverted repeat expansion and the degradation of the ndh genes in Strumaria truncata Jacq

Amaryllidaceae is a widespread and distinctive plant family contributing both food and ornamental plants. Here we present an initial survey of plastomes across the family and report on both structural rearrangements and gene losses. Most plastomes in the family are of similar gene arrangement and content however some taxa have shown gains in plastome length while in several taxa there is evidence of gene loss. Strumaria truncata shows a substantial loss of ndh family genes while three other taxa show loss of cemA, which has been reported only rarely. Our sparse sampling of the family has detected sufficient variation to suggest further sampling across the family could be a rich source of new information on plastome variation and evolution.

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.

Morphological and genomic evidence for a new species of Corallorhiza (Orchidaceae: Epidendroideae) from SW China

Corallorhiza sinensis, a new species of mycoheterotrophic orchid from western Sichuan, China, is described and illustrated based on molecular and morphological evidence. It is morphologically similar to Corallorhiza trifida, but can be distinguished by bigger flowers, both sepals and petals with 3 veins, and longer lateral lobes of lip. To distinguish the new Corallorhiza species and explore its phylogenetic position within subtribe Calypsoinae, this study employed sequences of the nuclear ribosomal DNA (nrDNA) and whole plastome assembled from the genome skimming approach. The plastome is 148,124 bp in length, including a pair of inverted repeats (IRs) of 26,165 bp, a large single-copy region (LSC) of 82,207 bp, and a small single-copy region (SSC) of 13,587 bp. Further, phylogenetic analyses were performed using nrDNA sequence and 79 coding sequences (CDSs) from 26 complete plastomes of subtribe Calypsoinae. The phylogenetic tree based on nrDNA sequence suggested that Corallorhiza is a monophyletic group, and strongly support C. sinensis as sister to the rest species of Corallorhiza. The plastid tree showed that 10 Corallorhiza species grouped into two clades and C. sinensis is most closely related to the North American C. striata and C. bentleyi instead of Oreorchis foliosa and O . angustata in the same clade. The plastid tree and nrDNA tree indicate Oreorchis is a paraphyletic. Although the topological conflicts are displayed between plastome and nrDNA phylogenies of C. sinensis, it is still the most closely related to Corallorhiza. Comparative analysis showed that C. sinensis populations are characteristic of the intermediate morphological traits between Corallorhiza and Oreorchis. The finding of this new species from China shed new light on the phylogeny of Oreorchis and Corallorhiza.

Plastid phylogenomics of Pleurothallidinae (Orchidaceae): Conservative plastomes, new variable markers, and comparative analyses of plastid, nuclear, and mitochondrial data

We present the first comparative plastome study of Pleurothallidinae with analyses of structural and molecular characteristics and identification of the ten most-variable regions to be incorporated in future phylogenetic studies. We sequenced complete plastomes of eight species in the subtribe and compared phylogenetic results of these to parallel analyses of their nuclear ribosomal DNA operon (26S, 18S, and 5.8S plus associated spacers) and partial mitochondrial genome sequences (29–38 genes and partial introns). These plastomes have the typical quadripartite structure for which gene content is similar to those of other orchids, with variation only in the composition of the ndh genes. The independent loss of ndh genes had an impact on which genes border the inverted repeats and thus the size of the small single-copy region, leading to variation in overall plastome length. Analyses of 68 coding sequences indicated the same pattern of codon usage as in other orchids, and 13 protein-coding genes under positive selection were detected. Also, we identified 62 polymorphic microsatellite loci and ten highly variable regions, for which we designed primers. Phylogenomic analyses showed that the top ten mutational hotspots represent well the phylogenetic relationships found with whole plastome sequences. However, strongly supported incongruence was observed among plastid, nuclear ribosomal DNA operon, and mitochondrial DNA trees, indicating possible occurrence of incomplete lineage sorting and/or introgressive hybridization. Despite the incongruence, the mtDNA tree retrieved some clades found in other analyses. These results, together with performance in recent studies, support a future role for mitochondrial markers in Pleurothallidinae phylogenetics.

Plastid phylogenomics improves resolution of phylogenetic relationship in the Cheirostylis and Goodyera clades of Goodyerinae (Orchidoideae, Orchidaceae)

Goodyerinae are one of phylogenetically unresolved groups of Orchidaceae. The lack of resolution achieved through the analyses of previous molecular sequences from one or a few markers has long confounded phylogenetic estimation and generic delimitation. Here, we present large-scale phylogenomic data to compare the plastome structure of the two main clades (Goodyera and Cheirostylis) in this subtribe and further adopt two strategies, combining plastid coding sequences and the whole plastome, to investigate phylogenetic relationships. A total of 46 species in 16 genera were sampled, including 39 species in 15 genera sequenced in this study. The plastomes of heterotrophic species are not drastically reduced in overall size, but display a pattern congruent with a loss of photosynthetic function. The plastomes of autotrophic species ranged from 147 to 165 kb and encoded from 132 to 137 genes. Three unusual structural features were detected: a 1.0-kb inversion in the large single-copy region of Goodyera schlechtendaliana; the loss and/or pseudogenization of ndh genes only in two species, Cheirostylis chinensis and C. montana; and the expansion of inverted repeat regions and contraction of small single-copy region in Hetaeria oblongifolia. Phylogenomic analyses provided improved resolution for phylogenetic relationships. All genera were recovered as monophyletic, except for Goodyera and Hetaeria, which were each recovered as non-monophyletic. Nomenclatural changes are needed until the broader sampling and biparental inherited markers. This study provides a phylogenetic framework of Goodyerinae and insight into plastome evolution of Orchidaceae.