The Complete Plastome Sequences of Four Orchid Species: Insights into the Evolution of the Orchidaceae and the Utility of Plastomic Mutational Hotspots

Whole chloroplast genome sequence and phylogenetic analysis of Calanthe discolor (Orchidaceae)

The orchid Calanthe discolor, which has high ornamental and medicinal value, is mainly distributed in Zhejiang, Jiangsu, and southeast Hubei Provinces of China, as well as in Japan and the southern Korean peninsula. In this study, the whole chloroplast genome sequence of C. discolor was first assembled using high-throughput Illumina paired-end technology, providing data to evaluate the evolution of this species. The C. discolor chloroplast genome was158,286 bp long, including a large single-copy region of 87,095 bp, a small single-copy region of 18,407 bp, and two copies of a repeat region (26,392-bp each). The overall G + C content was 41.2%. A total of 133 genes were predicted from the genome, including 87 protein-coding genes, eight ribosomal RNAs, 38 transfer RNAs. Phylogenetic analysis indicated a close relationship between C. discolor and C. bicolor.

Plastome Evolution, Phylogenomics, and DNA Barcoding Investigation of Gastrochilus (Aeridinae, Orchidaceae), with a Focus on the Systematic Position of Haraella retrocalla

Gastrochilus is an orchid genus containing about 70 species in tropical and subtropical Asia with high morphological diversity. The phylogenetic relationships among this genus have not been fully resolved, and the plastome evolution has not been investigated either. In this study, five plastomes of Gastrochilus were newly reported, and sixteen plastomes of Gastrochilus were used to conduct comparative and phylogenetic analyses. Our results showed that the Gastrochilus plastomes ranged from 146,183 to 148,666 bp, with a GC content of 36.7-36.9%. There were 120 genes annotated, consisting of 74 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. No contraction and expansion of IR borders, gene rearrangements, or inversions were detected. Additionally, the repeat sequences and codon usage bias of Gastrochilus plastomes were highly conserved. Twenty hypervariable regions were selected as potential DNA barcodes. The phylogenetic relationships within Gastrochilus were well resolved based on the whole plastome, especially among main clades. Furthermore, both molecular and morphological data strongly supported Haraella retrocalla as a member of Gastrochilus (G. retrocallus).

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.

Comparative Analysis of Chloroplast Pan-Genomes and Transcriptomics Reveals Cold Adaptation in Medicago sativa

Alfalfa (Medicago sativa) is a perennial forage legume that is widely distributed all over the world; therefore, it has an extremely complex genetic background. Though population structure and phylogenetic studies have been conducted on a large group of alfalfa nuclear genomes, information about the chloroplast genomes is still lacking. Chloroplast genomes are generally considered to be conservative and play an important role in population diversity analysis and species adaptation in plants. Here, 231 complete alfalfa chloroplast genomes were successfully assembled from 359 alfalfa resequencing data, on the basis of which the alfalfa chloroplast pan-genome was constructed. We investigated the genetic variations of the alfalfa chloroplast genome through comparative genomic, genetic diversity, phylogenetic, population genetic structure, and haplotype analysis. Meanwhile, the expression of alfalfa chloroplast genes under cold stress was explored through transcriptome analysis. As a result, chloroplast genomes of 231 alfalfa lack an IR region, and the size of the chloroplast genome ranges from 125,192 bp to 126,105 bp. Using population structure, haplotypes, and construction of a phylogenetic tree, it was found that alfalfa populations could be divided into four groups, and multiple highly variable regions were found in the alfalfa chloroplast genome. Transcriptome analysis showed that tRNA genes were significantly up-regulated in the cold-sensitive varieties, while rps7, rpl32, and ndhB were down-regulated, and the editing efficiency of ycf1, ycf2, and ndhF was decreased in the cold-tolerant varieties, which may be due to the fact that chloroplasts store nutrients through photosynthesis to resist cold. The huge number of genetic variants in this study provide powerful resources for molecular markers.

Comparative analysis and characterization of the chloroplast genome of Krascheninnikovia ceratoides (Amarathaceae): a xerophytic semi-shrub exhibiting drought resistance and high-quality traits

BACKGROUND

Krascheninnikovia ceratoides, a perennial halophytic semi-shrub belonging to the genus Krascheninnikovia (Amarathaceae), possesses noteworthy ecological, nutritional, and economic relevance. This species is primarily distributed across arid, semi-arid, and saline-alkaline regions of the Eurasian continent, encompassing Inner Mongolia, Xinjiang, Qinghai, Gansu, Ningxia, and Tibet.

RESULTS

We reported the comprehensive chloroplast (cp) genome of K. ceratoides, characterized by a circular conformation spanning 151,968 bp with a GC content of 36.60%. The cp genome encompassed a large single copy (LSC, 84,029 bp), a small single copy (SSC, 19,043 bp), and a pair of inverted repeats (IRs) regions (24,448 bp each). This genome harbored 128 genes and encompassed 150 simple sequence repeats (SSRs). Through comparative analyses involving cp genomes from other Cyclolobeae (Amarathaceae) taxa, we observed that the K. ceratoides cp genome exhibited high conservation, with minor divergence events in protein-coding genes (PCGs) accD, matK, ndhF, ndhK, ycf1, and ycf2. Phylogenetic reconstructions delineated K. ceratoides as the sister taxon to Atriplex, Chenopodium, Dysphania, and Suaeda, thus constituting a robust clade. Intriguingly, nucleotide substitution ratios (Ka/Ks) between K. ceratoides and Dysphania species for ycf1 and ycf2 genes surpassed 1.0, indicating the presence of positive selection pressure on these loci.

CONCLUSIONS

The findings of this study augment the genomic repository for the Amarathaceae family and furnish crucial molecular instruments for subsequent investigations into the ecological adaptation mechanisms of K. ceratoides within desert ecosystems.

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.

Phylogenetics Study to Compare Chloroplast Genomes in Four Magnoliaceae Species

Magnoliaceae, a family of perennial woody plants, contains several endangered species whose taxonomic status remains ambiguous. The study of chloroplast genome information can help in the protection of Magnoliaceae plants and confirmation of their phylogenetic relationships. In this study, the chloroplast genomes were sequenced, assembled, and annotated in Woonyoungia septentrionalis and three Michelia species (Michelia champaca, Michelia figo, and Michelia macclurei). Comparative analyses of genomic characteristics, repetitive sequences, and sequence differences were performed among the four Magnoliaceae plants, and phylogenetic relationships were constructed with twenty different magnolia species. The length of the chloroplast genomes varied among the four studied species ranging from 159,838 bp (Woonyoungia septentrionalis) to 160,127 bp (Michelia macclurei). Four distinct hotspot regions were identified based on nucleotide polymorphism analysis. They were petA-psbJ, psbJ-psbE, ndhD-ndhE, and rps15-ycf1. These gene fragments may be developed and utilized as new molecular marker primers. By using Liriodendron tulipifera and Liriodendron chinense as outgroups reference, a phylogenetic tree of the four Magnoliaceae species and eighteen other Magnoliaceae species was constructed with the method of Shared Coding Sequences (CDS). Results showed that the endangered species, W. septentrionalis, is relatively genetically distinct from the other three species, indicating the different phylogenetic processes among Magnoliaceae plants. Therefore, further genetic information is required to determine the relationships within Magnoliaceae. Overall, complete chloroplast genome sequences for four Magnoliaceae species reported in this paper have shed more light on phylogenetic relationships within the botanical group.

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.

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

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

The climate changes promoted the chloroplast genomic evolution of Dendrobium orchids among multiple photosynthetic pathways

Dendrobium orchids have multiple photosynthetic pathways, which can be used as a model system for studying the evolution of crassulacean acid metabolism (CAM). In this study, based on the results of the net photosynthetic rates (P_n), we classified Dendrobium species into three photosynthetic pathways, then employed and compared their chloroplast genomes. The Dendrobium chloroplast genomes have typical quartile structures, ranging from 150,841-153,038 bp. The apparent differences in GC content, sequence variability, and IR junctions of SSC/IR_B junctions (J_SBs) were measured within chloroplast genomes among different photosynthetic pathways. The phylogenetic analysis has revealed multiple independent CAM origins among the selected Dendrobium species. After counting insertions and deletions (InDels), we found that the occurrence rates and distribution densities among different photosynthetic pathways were inconsistent. Moreover, the evolution patterns of chloroplast genes in Dendrobium among three photosynthetic pathways were also diversified. Considering the diversified genome structure variations and the evolution patterns of protein-coding genes among Dendrobium species, we proposed that the evolution of the chloroplast genomes was disproportional among different photosynthetic pathways. Furthermore, climatic correlation revealed that temperature and precipitation have influenced the distribution among different photosynthetic pathways and promoted the foundation of CAM pathway in Dendrobium orchids. Based on our study, we provided not only new insights into the CAM evolution of Dendrobium but also provided beneficial genetic data resources for the further systematical study of Dendrobium.

Comparative chloroplast genomics of 34 species in subtribe Swertiinae (Gentianaceae) with implications for its phylogeny

BACKGROUND

Subtribe Swertiinae, a medicinally significant and highly speciose Subtribe of family Gentianaceae. Despite previous extensive studies based on both morphology and molecular data, intergeneric and infrageneric relationships within subtribe Swertiinae remain controversial.

METHODS

Here, we employed four newly generated Swertia chloroplast genomes with thirty other published genomes to elucidate their genomic characteristics.

RESULTS

The 34 chloroplast genomes were small and ranged in size from 149,036 to 154,365 bp, each comprising two inverted repeat regions (size range 25,069-26,126 bp) that separated large single-copy (80,432-84,153 bp) and small single-copy (17,887-18,47 bp) regions, and all the chloroplast genomes showed similar gene orders, contents, and structures. These chloroplast genomes contained 129-134 genes each, including 84-89 protein-coding genes, 37 tRNAs, and 8 rRNAs. The chloroplast genomes of subtribe Swertiinae appeared to have lost some genes, such as rpl33, rpl2 and ycf15 genes. Comparative analyses revealed that two mutation hotspot regions (accD-psaI and ycf1) could serve as effective molecular markers for further phylogenetic analyses and species identification in subtribe Swertiinae. Positive selection analyses showed that two genes (ccsA and psbB) had high Ka/Ks ratios, indicating that chloroplast genes may have undergone positive selection in their evolutionary history. Phylogenetic analysis showed that the 34 subtribe Swertiinae species formed a monophyletic clade, with Veratrilla, Gentianopsis and Pterygocalyx located at the base of the phylogenetic tree. Some genera of this subtribe, however, were not monophyletic, including Swertia, Gentianopsis, Lomatogonium, Halenia, Veratrilla and Gentianopsis. In addition, our molecular phylogeny was consistent with taxonomic classification of subtribe Swertiinae in the Roate group and Tubular group. The results of molecular dating showed that the divergence between subtrib Gentianinae and subtrib Swertiinae was estimated to occur in 33.68 Ma. Roate group and Tubular group in subtribe Swertiinae approximately diverged in 25.17 Ma.

CONCLUSION

Overall, our study highlighted the taxonomic utility of chloroplast genomes in subtribe Swertiinae, and the genetic markers identified here will facilitate future studies on the evolution, conservation, population genetics, and phylogeography of subtribe Swertiinae species.

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

Chloroplast phylogenomics and the taxonomy of Saxifraga section Ciliatae (Saxifragaceae)

Comprising ca. 200 species, Saxifraga sect. Ciliatae is the most species-rich section of Saxifraga s.str., whose center of diversity is in the Tibeto-Himalayan region. The infra-sectional classification of sect. Ciliatae is still in debate due to the high level of species richness, as well as remarkable variations of habitat, morphology, physiology and life cycles. Subdivisions of this section proposed in various taxonomic systems have not been adequately tested in previous phylogenetic studies, partly due to low taxonomic sampling density, but also to the use of few DNA markers. In order to achieve a more robust infra-sectional classification of sect. Ciliatae, complete chloroplast genomes of 94 taxa from this section were analyzed, of which 93 were newly sequenced, assembled and annotated. The length of the 94 plastomes of sect. Ciliatae taxa range from 143,479 to 159,938 bp, encoding 75 to 79 unique protein-coding genes (PCGs). Analyses of the 94 plastomes revealed high conservation in structural organization, gene arrangement, and gene content. Gene loss and changes of IR boundaries were detected but in extremely low frequency. The molecular phylogenetic tree from concatenated PCGs and complete chloroplast genome sequences exhibited high resolution and support values and confirms that sect. Ciliatae is monophyletic. Three well-supported clades were revealed within the section that agree relatively well with the subsectional taxonomy of Gornall (1987), but some minor modifications should be made. Firstly, the monotypic subsection Cinerascentes should be abandoned and its constituent species, S. cinerascens, assigned to subsect. Gemmiparae. Secondly, subsections Rosulares and Serpyllifoliae should be merged and become subsect. Rosulares. Section Ciliatae thus comprises: subsect. Hirculoideae Engl. & Irmsch.; subsect. Rosulares Gornall; subsect. Gemmiparae Engl. & Irmsch.; subsect. Flagellares (C. B. Clarke) Engl. & Irmsch. and subsect. Hemisphaericae (Engl. & Irmsch.) Gornall.

In-depth analysis of genomes and functional genomics of orchid using cutting-edge high-throughput sequencing

High-throughput sequencing technology has been facilitated the development of new methodologies and approaches for studying the origin and evolution of plant genomes and subgenomes, population domestication, and functional genomics. Orchids have tens of thousands of members in nature. Many of them have promising application potential in the extension and conservation of the ecological chain, the horticultural use of ornamental blossoms, and the utilization of botanical medicines. However, a large-scale gene knockout mutant library and a sophisticated genetic transformation system are still lacking in the improvement of orchid germplasm resources. New gene editing tools, such as the favored CRISPR-Cas9 or some base editors, have not yet been widely applied in orchids. In addition to a large variety of orchid cultivars, the high-precision, high-throughput genome sequencing technology is also required for the mining of trait-related functional genes. Nowadays, the focus of orchid genomics research has been directed to the origin and classification of species, genome evolution and deletion, gene duplication and chromosomal polyploidy, and flower morphogenesis-related regulation. Here, the progressing achieved in orchid molecular biology and genomics over the past few decades have been discussed, including the evolution of genome size and polyploidization. The frequent incorporation of LTR retrotransposons play important role in the expansion and structural variation of the orchid genome. The large-scale gene duplication event of the nuclear genome generated plenty of recently tandem duplicated genes, which drove the evolution and functional divergency of new genes. The evolution and loss of the plastid genome, which mostly affected genes related to photosynthesis and autotrophy, demonstrated that orchids have experienced more separate transitions to heterotrophy than any other terrestrial plant. Moreover, large-scale resequencing provide useful SNP markers for constructing genetic maps, which will facilitate the breeding of novel orchid varieties. The significance of high-throughput sequencing and gene editing technologies in the identification and molecular breeding of the trait-related genes in orchids provides us with a representative trait-improving gene as well as some mechanisms worthy of further investigation. In addition, gene editing has promise for the improvement of orchid genetic transformation and the investigation of gene function. This knowledge may provide a scientific reference and theoretical basis for orchid genome studies.

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.

The complete chloroplast genome and phylogenetic relationship of Apocynum pictum (Apocynaceae), a Central Asian shrub and second-class national protected species of western China

Apocynum pictum of the dogbane family, Apocynaceae, is a perennial semi-shrub species of ecological, medicinal, and economic value. It is mainly distributed in semi-arid, saline-alkaline, and desert regions of Xinjiang, Qinghai, and Gansu of western China and adjacent regions from Kazakhstan and Mongolia. Here, we reported the complete chloroplast (cp) genome of A. pictum for the first time, and we found that it had a circular structure with an estimated length of 150,749 bp and a GC content of 38.3%. The cp genome was composed of a large single copy (LSC), a single small single copy (SSC), and two inverted repeat (IR) regions, which were 81,888 bp, 17,251 bp and 25,805 bp long, respectively. The cp genome of A. pictum encoded 134 genes and contained 66 simple sequence repeats (SSRs). A comparative analysis with other cp genomes from Apocynaceae indicated that the cp genome of A. pictum was very conserved, except for subtle differences occurring in the protein-coding genes accD, ndhF, rpl22, rpl32, rpoC2, ycf1 and ycf2. A phylogenetic reconstruction showed that A. pictum and A. venetum were sister species, forming a strongly supported clade with Trachelospermum. Interestingly, nucleotide substitution ratios (Ka/Ks) between A. pictum and A. venetum on accD and ndhF were >1.0, suggesting positive selective pressure on these genes. Our result enriches the genomic resources for the diverse dogbane family and provides critical molecular resources to develop future studies on ecological adaptation to desert habitats in Apocynum.

Comparative Analyses of Chloroplast Genomes Provide Comprehensive Insights into the Adaptive Evolution of Paphiopedilum (Orchidaceae)

An elucidation of how the selection pressures caused by habitat environments affect plant plastid genomes and lead to the adaptive evolution of plants, is a very intense area of research in evolutionary biology. The genus Paphiopedilum is a predominant group of orchids that includes over 66 species with high horticultural and ornamental value. However, owing to the destructive exploitation and habitat deterioration of wild germplasm resources of Paphiopedilum, it needs more molecular genetic resources and studies on this genus. The chloroplast is cytoplasmically inherited and often used in evolutionary studies. Thus, for this study, we newly sequenced, assembled and annotated five chloroplast genomes of the Paphiopedilum species. The size of these genomes ranged from 155,886 bp (P. henryanum) to 160,503 bp (P. ‘GZSLKY’ Youyou) and they contained 121–122 genes, which consisted of 76 protein coding genes, eight ribosomal RNAs, and 37–38 transfer RNAs. Combined with the other 14 Paphiopedilum species, the characteristics of the repeat sequences, divergent hotspot regions, and the condo usage bias were evaluated and identified, respectively. The gene transfer analysis showed that some fragments of the ndh and ycf gene families were shared by both the chloroplast and nucleus. Although the genomic structure and gene content was conserved, there was a significant boundary shift caused by the inverted repeat (IR) expansion and small single copy (SSC) contraction. The lower GC content and loss of ndh genes could be the result of adaptive evolutionary responses to its unique habitats. The genes under positive selection, including accD, matK, psbM, rpl20, rps12, ycf1, and ycf2 might be regarded as potential candidate genes for further study, which significantly contribute to the adaptive evolution of Paphiopedilum.

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.

Impact of culture vessels on micro-morphological features of in vitro Dendrobium Sabin Blue orchid

Different designs of the plant tissue culture vessel, such as size, material, and shape, may alter its microenvironment atmosphere. The present study was conducted on protocorm-like bodies (PLBs) of Dendrobium Sabin Blue orchid to determine the development of PLBs on plastic and glass culture vessels of different sizes. PLBs were cultured in half-strength Murashige and Skoog (MS) medium with the same initial weight of 0.5 g in 10 replicates. The growth index of the PLBs was calculated after 11 weeks to study their growth in every vessel; additionally, biochemical analysis was performed to determine carbohydrate content, proline concentration, and photosynthesis pigments in the PLBs. Scanning electron microscopy (SEM) was performed to study stomata development on PLBs in each vessel, and histological analyses were conducted to study the cell structure. Overall, the PLBs cultured in a large 470 ml plastic vessel showed successful growth with a high growth index, high carbohydrate content, low-stress condition, and high chlorophyll content. SEM confirmed that the presence of trichome and rhizoid in PLBs cultured in the 470 ml plastic vessel. Histological analysis showed the formation of the shoot on the PLBs and the presence of starch granules. Thus, the use of plastic as a culture vessel provides a good impact for culturing PLBs and has low cost.

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.

Advances in Breeding, Bioprospecting, and In Vitro Culture of Laelia Orchid Species

Orchids (Orchidaceae) are plants that are highly appreciated by their beautiful flowers worldwide. Moreover, they represent a source of metabolites with applications in medicine and biotechnology. Within the Orchidaceae family, the Laelia genus is a group of orchid species from the Neotropics and is probably one of the most representative genera of America. Laelia orchids are cultivated by their splendid flowers and are widely used in orchid breeding. Here, we revise the use of the Laelia genus in orchid breeding and metabolite bioprospecting. We also analyze the use of plant tissue culture (PTC) as an alternative to conventional propagation and as a strategy for the recovery of those Laelia species threatened with extinction. We summarize and discuss the recent advances in the application of different PTC techniques for mass multiplication based on asymbiotic germination, organogenesis, protocorm-like bodies development, and somatic embryogenesis, and the advances of in vitro conservation by cryoconservation and the use of slow-growth promoting hormones. Finally, we suggest future directions and venues in research for Laelia species.

Geographical Distribution and Relationship with Environmental Factors of Paphiopedilum Subgenus Brachypetalum Hallier (Orchidaceae) Taxa in Southwest China

The determination of the geographic distributions of orchid species and their relationships with environmental factors are considered fundamental to their conservation. Paphiopedilum subgenus Brachypetalum is one of the most primitive, ornamental, and threatened groups of Orchidaceae. However, little is known about the distribution of Brachypetalum orchids and how they are influenced by environmental factors. In this study, we developed a database on the geographical distribution of Brachypetalum orchids based on a large-scale field investigation in the Guangxi, Guizhou, and Yunnan provinces of southwest China (2019–2020). Using this database, we first adopted the nonparametric Mann–Whitney U test to analyze the differences in the geographical distributions and growth environments of Brachypetalum orchids. In addition, we also used the method of principal component analysis (PCA) to explore distribution patterns of Brachypetalum orchids in relation to environmental factors (topography, climate, anthropogenic disturbance, productivity, and soil) in southwest China. Our results indicated that Brachypetalum orchid species were mainly distributed in the karst limestone habitats of southwest China. In general, there were 194 existing localities with the occurrence of seven target orchids in the investigated area. Of the discovered species in our study, 176 locations (~90.7%) were distributed primarily in the karst habitat. Among them, the range of 780–1267 m was the most concentrated elevation of Brachypetalum orchids. In addition, the findings also suggested that the distribution of Brachypetalum orchids in southwest China was relatively scattered in geographical space. However, the density of the distribution of Brachypetalum orchids was high, between 104° and 108° E and between 25° and 26° N. The results of the Mann–Whitney U test revealed that there are obviously different geographical distributions and growth environments of Brachypetalum in southwest China. More specifically, we found some extremely significant differences (p < 0.001) in elevation, mean diurnal range, precipitation of coldest quarter, solar radiation, and exchangeable Ca2+ between the provinces of southwest China. The PCA analysis revealed that elevation, solar radiation, temperature (mean diurnal range, annual temperature range) and precipitation (precipitation seasonality, precipitation of the warmest quarter) were found to be the most significant factors in determining Brachypetalum orchids’ distribution. These findings have implications in assessing conservation effectiveness and determining niche breadth to better protect the populations of these Brachypetalum orchid species in the future.

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.

Comparative Chloroplast Genomics and Phylogenetic Analysis of Zygophyllum (Zygophyllaceae) of China

The genus Zygophyllum comprises over 150 species within the plant family Zygophyllaceae. These species predominantly grow in arid and semiarid areas, and about 20 occur in northwestern China. In this study, we sampled 24 individuals of Zygophyllum representing 15 species and sequenced their complete chloroplast (cp) genomes. For comparison, we also sequenced cp genomes of two species of Peganum from China representing the closely allied family, Nitrariaceae. The 24 cp genomes of Zygophyllum were smaller and ranged in size from 104,221 to 106,286 bp, each containing a large single-copy (LSC) region (79,245-80,439 bp), a small single-copy (SSC) region (16,285-17,146 bp), and a pair of inverted repeat (IR) regions (3,792-4,466 bp). These cp genomes contained 111-112 genes each, including 74-75 protein-coding genes (PCGs), four ribosomal RNA genes, and 33 transfer RNA genes, and all cp genomes showed similar gene order, content, and structure. The cp genomes of Zygophyllum appeared to lose some genes such as ndh genes and rRNA genes, of which four rRNA genes were in the SSC region, not in the IR regions. However, the SC and IR regions had greater similarity within Zygophyllum than between the genus and Peganum. We detected nine highly variable intergenic spacers: matK-trnQ, psaC-rps15, psbZ-trnG, rps7-trnL, rps15-trnN, trnE-trnT, trnL-rpl32, trnQ-psbK, and trnS-trnG. Additionally, we identified 156 simple sequence repeat (cpSSR) markers shared among the genomes of the 24 Zygophyllum samples and seven cpSSRs that were unique to the species of Zygophyllum. These markers may be useful in future studies on genetic diversity and relationships of Zygophyllum and closely related taxa. Using the sequenced cp genomes, we reconstructed a phylogeny that strongly supported the division of Chinese Zygophyllum into herbaceous and shrubby clades. We utilized our phylogenetic results along with prior morphological studies to address several remaining taxonomic questions within Zygophyllum. Specifically, we found that Zygophyllum kaschgaricum is included within Zygophyllum xanthoxylon supporting the present treatment of the former genus Sarcozygium as a subgenus within Zygophyllum. Our results provide a foundation for future research on the genetic resources of Zygophyllum.

Phylogenetic analysis and development of molecular markers for five medicinal Alpinia species based on complete plastome sequences

BACKGROUND

Alpinia species are widely used as medicinal herbs. To understand the taxonomic classification and plastome evolution of the medicinal Alpinia species and correctly identify medicinal products derived from Alpinia species, we systematically analyzed the plastome sequences from five Alpinia species. Four of the Alpinia species: Alpinia galanga (L.) Willd., Alpinia hainanensis K.Schum., Alpinia officinarum Hance, and Alpinia oxyphylla Miq., are listed in the Chinese pharmacopeia. The other one, Alpinia nigra (Gaertn.) Burtt, is well known for its medicinal values.

RESULTS

The four Alpinia species: A. galanga, A. nigra, A. officinarum, and A. oxyphylla, were sequenced using the Next-generation sequencing technology. The plastomes were assembled using Novoplasty and annotated using CPGAVAS2. The sizes of the four plastomes range from 160,590 bp for A. galanga to 164,294 bp for A. nigra, and display a conserved quadripartite structure. Each of the plastomes encodes a total of 111 unique genes, including 79 protein-coding, 28 tRNA, and four rRNA genes. In addition, 293-296 SSRs were detected in the four plastomes, of which the majority are mononucleotides Adenine/Thymine and are found in the noncoding regions. The long repeat analysis shows all types of repeats are contained in the plastomes, of which palindromic repeats occur most frequently. The comparative genomic analyses revealed that the pair of the inverted repeats were less divergent than the single-copy region. Analysis of sequence divergence on protein-coding genes showed that two genes (accD and ycf1) had undergone positive selection. Phylogenetic analysis based on coding sequence of 77 shared plastome genes resolves the molecular phylogeny of 20 species from Zingiberaceae. In particular, molecular phylogeny of four sequenced Alpinia species (A. galanga, A. nigra, A. officinarum, and A. oxyphylla) based on the plastome and nuclear sequences showed congruency. Furthermore, a comparison of the four newly sequenced Alpinia plastomes and one previously reported Alpinia plastomes (accession number: NC_048461) reveals 59 highly divergent intergenic spacer regions. We developed and validated two molecular markers Alpp and Alpr, based on two regions: petN-psbM and psaJ-rpl33, respectively. The discrimination success rate was 100 % in validation experiments.

CONCLUSIONS

The results from this study will be invaluable for ensuring the effective and safe uses of Alpinia medicinal products and for the exploration of novel Alpinia species to improve human health.

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.

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.

Reference-based assembly of chloroplast genome from leaf transcriptome data of Pterocarpus santalinus

Chloroplast genome sequencing is an essential tool to understand genome evolution and phylogenetic relationship. The available methods for constructing chloroplast genome include chloroplast enrichment followed by long overlapping PCR or extraction and assembly of chloroplast-specific reads from whole-genome datasets. In the present study, we propose an alternate strategy of extraction and assembly of chloroplast-specific reads from leaf transcriptome data of Pterocarpus santalinus using bowtie2 aligner program. The assembled genome was compared with the published chloroplast genome of P. santalinus for genome size, number of predicted genes, microsatellite repeat motifs, and nucleotide repeats. A near-complete chloroplast genome was assembled from the transcriptome reads. The proposed method requires less computational time and know-how, limited virtual memory, and is cost-effective when compared to whole-genome sequencing. Assembly of Cp genome from transcriptome data will enhance the resolution of phylogenetic studies through comparative plastome analysis, facilitate accurate species/genotype discrimination and accelerate the development of transplastomic plants with enhanced biotic and abiotic tolerance.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02943-0.

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.

The chloroplast genome evolution of Venus slipper (Paphiopedilum): IR expansion, SSC contraction, and highly rearranged SSC regions

BACKGROUND

Paphiopedilum is the largest genus of slipper orchids. Previous studies showed that the phylogenetic relationships of this genus are not well resolved, and sparse taxon sampling documented inverted repeat (IR) expansion and small single copy (SSC) contraction of the chloroplast genomes of Paphiopedilum.

RESULTS

Here, we sequenced, assembled, and annotated 77 plastomes of Paphiopedilum species (size range of 152,130 - 164,092 bp). The phylogeny based on the plastome resolved the relationships of the genus except for the phylogenetic position of two unstable species. We used phylogenetic and comparative genomic approaches to elucidate the plastome evolution of Paphiopedilum. The plastomes of Paphiopedilum have a conserved genome structure and gene content except in the SSC region. The large single copy/inverted repeat (LSC/IR) boundaries are relatively stable, while the boundaries of the inverted repeat and small single copy region (IR/SSC) varied among species. Corresponding to the IR/SSC boundary shifts, the chloroplast genomes of the genus experienced IR expansion and SSC contraction. The IR region incorporated one to six genes of the SSC region. Unexpectedly, great variation in the size, gene order, and gene content of the SSC regions was found, especially in the subg. Parvisepalum. Furthermore, Paphiopedilum provides evidence for the ongoing degradation of the ndh genes in the photoautotrophic plants. The estimated substitution rates of the protein coding genes show accelerated rates of evolution in clpP, psbH, and psbZ. Genes transferred to the IR region due to the boundary shift also have higher substitution rates.

CONCLUSIONS

We found IR expansion and SSC contraction in the chloroplast genomes of Paphiopedilum with dense sampling, and the genus shows variation in the size, gene order, and gene content of the SSC region. This genus provides an ideal system to investigate the dynamics of plastome evolution.

The chloroplast genome evolution of Venus slipper (Paphiopedilum): IR expansion, SSC contraction, and highly rearranged SSC regions

BACKGROUND

Paphiopedilum is the largest genus of slipper orchids. Previous studies showed that the phylogenetic relationships of this genus are not well resolved, and sparse taxon sampling documented inverted repeat (IR) expansion and small single copy (SSC) contraction of the chloroplast genomes of Paphiopedilum.

RESULTS

Here, we sequenced, assembled, and annotated 77 plastomes of Paphiopedilum species (size range of 152,130 - 164,092 bp). The phylogeny based on the plastome resolved the relationships of the genus except for the phylogenetic position of two unstable species. We used phylogenetic and comparative genomic approaches to elucidate the plastome evolution of Paphiopedilum. The plastomes of Paphiopedilum have a conserved genome structure and gene content except in the SSC region. The large single copy/inverted repeat (LSC/IR) boundaries are relatively stable, while the boundaries of the inverted repeat and small single copy region (IR/SSC) varied among species. Corresponding to the IR/SSC boundary shifts, the chloroplast genomes of the genus experienced IR expansion and SSC contraction. The IR region incorporated one to six genes of the SSC region. Unexpectedly, great variation in the size, gene order, and gene content of the SSC regions was found, especially in the subg. Parvisepalum. Furthermore, Paphiopedilum provides evidence for the ongoing degradation of the ndh genes in the photoautotrophic plants. The estimated substitution rates of the protein coding genes show accelerated rates of evolution in clpP, psbH, and psbZ. Genes transferred to the IR region due to the boundary shift also have higher substitution rates.

CONCLUSIONS

We found IR expansion and SSC contraction in the chloroplast genomes of Paphiopedilum with dense sampling, and the genus shows variation in the size, gene order, and gene content of the SSC region. This genus provides an ideal system to investigate the dynamics of plastome evolution.

Plastid phylogenomics resolves ambiguous relationships within the orchid family and provides a solid timeframe for biogeography and macroevolution

Recent phylogenomic analyses based on the maternally inherited plastid organelle have enlightened evolutionary relationships between the subfamilies of Orchidaceae and most of the tribes. However, uncertainty remains within several subtribes and genera for which phylogenetic relationships have not ever been tested in a phylogenomic context. To address these knowledge-gaps, we here provide the most extensively sampled analysis of the orchid family to date, based on 78 plastid coding genes representing 264 species, 117 genera, 18 tribes and 28 subtribes. Divergence times are also provided as inferred from strict and relaxed molecular clocks and birth-death tree models. Our taxon sampling includes 51 newly sequenced plastid genomes produced by a genome skimming approach. We focus our sampling efforts on previously unplaced clades within tribes Cymbidieae and Epidendreae. Our results confirmed phylogenetic relationships in Orchidaceae as recovered in previous studies, most of which were recovered with maximum support (209 of the 262 tree branches). We provide for the first time a clear phylogenetic placement for Codonorchideae within subfamily Orchidoideae, and Podochilieae and Collabieae within subfamily Epidendroideae. We also identify relationships that have been persistently problematic across multiple studies, regardless of the different details of sampling and genomic datasets used for phylogenetic reconstructions. Our study provides an expanded, robust temporal phylogenomic framework of the Orchidaceae that paves the way for biogeographical and macroevolutionary studies.

Three Novel Biphenanthrene Derivatives and a New Phenylpropanoid Ester from Aerides multiflora and Their α-Glucosidase Inhibitory Activity

A phytochemical investigation on the whole plants of Aerides multiflora revealed the presence of three new biphenanthrene derivatives named aerimultins A-C (1-3) and a new natural phenylpropanoid ester dihydrosinapyl dihydroferulate (4), together with six known compounds (5-10). The structures of the new compounds were elucidated by analysis of their spectroscopic data. All of the isolates were evaluated for their α-glucosidase inhibitory activity. Aerimultin C (3) showed the most potent activity. The other compounds, except for compound 4, also exhibited stronger activity than the positive control acarbose. Compound 3 showed non-competitive inhibition of the enzyme as determined from a Lineweaver-Burk plot. This study is the first phytochemical and biological investigation of A. multiflora.

Plastome structure and adaptive evolution of Calanthe s.l. species

Calanthe s.l. is the most diverse group in the tribe Collabieae (Orchidaceae), which are pantropical in distribution. Illumina sequencing followed by de novo assembly was used in this study, and the plastid genetic information of Calanthe s.l. was used to investigate the adaptive evolution of this taxon. Herein, the complete plastome of five Calanthe s.l. species (Calanthe davidii, Styloglossum lyroglossa, Preptanthe rubens, Cephalantheropsis obcordata, and Phaius tankervilliae) were determined, and the two other published plastome sequences of Calanthe s.l. were added for comparative analyses to examine the evolutionary pattern of the plastome in the alliance. The seven plastomes ranged from 150,181 bp (C. delavayi) to 159,014 bp (C. davidii) in length and were all mapped as circular structures. Except for the three ndh genes (ndhC, ndhF, and ndhK) lost in C. delavayi, the remaining six species contain identical gene orders and numbers (115 gene). Nucleotide diversity was detected across the plastomes, and we screened 14 mutational hotspot regions, including 12 non-coding regions and two gene regions. For the adaptive evolution investigation, three species showed positive selected genes compared with others, C. obcordata (cemA), S. lyroglossa (infA, ycf1 and ycf2) and C. delavayi (nad6 and ndhB). Six genes were under site-specific positive selection in Calanthe s.l., namely, accD, ndhB, ndhD, rpoC2, ycf1, and ycf2, most of which are involved in photosynthesis. These results, including the new plastomes, provide resources for the comparative plastome, breeding, and plastid genetic engineering of orchids and flowering plants.

The large single-copy (LSC) region functions as a highly effective and efficient molecular marker for accurate authentication of medicinal Dendrobium species

Having great medicinal values, Dendrobium species of "Fengdou" (DSFs) are a taxonomically complex group in Dendrobium genus including many closely related and recently diverged species. Traditionally used DNA markers have been proved to be insufficient in authenticating many species of this group. Here, we investigated 101 complete plastomes from 23 DSFs, comprising 72 newly sequenced and 29 documented, which all exhibited well-conserved genomic organization and gene order. Plastome-wide comparison showed the co-occurrence of single nucleotide polymorphisms (SNPs) and insertions/deletions (indels), which can be explained by both the repeat-associated and indel-associated mutation hypotheses. Moreover, guanine-cytosine (GC) content was found to be negatively correlated with the three divergence variables (SNPs, indels and repeats), indicating that GC content may reflect the level of the local sequence divergence. Our species authentication analyses revealed that the relaxed filtering strategies of sequence alignment had no negative impact on species identification. By assessing the maximum likelihood (ML) trees inferred from different datasets, we found that the complete plastome and large single-copy (LSC) datasets both successfully identified all 23 DSFs with the maximum bootstrap values. However, owing to the high efficiency of LSC in species identification, we recommend using LSC for accurate authentication of DSFs.

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

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

The Complete Chloroplast Genome of Arabidopsis thaliana Isolated in Korea (Brassicaceae): An Investigation of Intraspecific Variations of the Chloroplast Genome of Korean A. thaliana

Arabidopsis thaliana (L.) Heynh. is a model organism of plant molecular biology. More than 1,700 whole genome sequences have been sequenced, but no Korean isolate genomes have been sequenced thus far despite the fact that many A. thaliana isolated in Japan and China have been sequenced. To understand the genetic background of Korean natural A. thaliana (named as 180404IB4), we presented its complete chloroplast genome, which is 154,464 bp long and has four subregions: 85,164 bp of large single copy (LSC) and 17,781 bp of small single copy (SSC) regions are separated by 26,257 bp of inverted repeat (IRs) regions including 130 genes (85 protein-coding genes, eight rRNAs, and 37 tRNAs). Fifty single nucleotide polymorphisms (SNPs) and 14 insertion and deletions (INDELs) are identified between 180404IB4 and Col0. In addition, 101 SSRs and 42 extendedSSRs were identified on the Korean A. thaliana chloroplast genome, indicating a similar number of SSRs on the rest five chloroplast genomes with a preference of sequence variations toward the SSR region. A nucleotide diversity analysis revealed two highly variable regions on A. thaliana chloroplast genomes. Phylogenetic trees with three more chloroplast genomes of East Asian natural isolates show that Korean and Chinese natural isolates are clustered together, whereas two Japanese isolates are not clustered, suggesting the need for additional investigations of the chloroplast genomes of East Asian isolates.

The Complete Chloroplast Genome Sequences of 14 Curcuma Species: Insights Into Genome Evolution and Phylogenetic Relationships Within Zingiberales

Zingiberaceae is taxonomically complex family where species are perennial herb. However, lack of chloroplast genomic information severely hinders our understanding of Zingiberaceae species in the research of evolution and phylogenetic relationships. In this study, the complete chloroplast (cp) genomes of fourteen Curcuma species were assembled and characterized using next-generation sequencing. We compared the genome features, repeat sequences, sequence divergence, and constructed the phylogenetic relationships of the 25 Zingiberaceae species. In each Zingiberaceae species, the 25 complete chloroplast genomes ranging from 155,890 bp (Zingiber spectabile) to 164,101 bp (Lanxangia tsaoko) contained 111 genes consisting of 77 protein coding genes, 4 ribosomal RNAs and 30 transfer RNAs. These chloroplast genomes are similar to most angiosperm that consisted of a four-part circular DNA molecules. Moreover, the characteristics of the long repeats sequences and simple sequence repeats (SSRs) were found. Six divergent hotspots regions (matK-trnk, Rps16-trnQ, petN-psbM, rpl32, ndhA, and ycf1) were identified in the 25 Zingiberaceae chloroplast genomes, which could be potential molecular markers. In addition to Wurfbainia longiligularis, the ψycf1 was discovered among the 25 Zingiberaceae species. The shared protein coding genes from 52 Zingiberales plants and four other family species as out groups were used to construct phylogenetic trees distinguished by maximum parsimony (MP), maximum likelihood (ML) and Bayesian inference (BI) and showed that Musaceae was the basal group in Zingiberales, and Curcuma had a close relationship with Stahlianthu. Besides this, Curcuma flaviflora was clustered together with Zingiber. Its distribution area (Southeast Asia) overlaps with the latter. Maybe hybridization occur in related groups within the same region. This may explain why Zingiberaceae species have a complex phylogeny, and more samples and genetic data were necessary to confirm their relationship. This study provide the reliable information and high-quality chloroplast genomes and genome resources for future Zingiberaceae research.

The plastid NAD(P)H dehydrogenase-like complex: structure, function and evolutionary dynamics

The thylakoid NAD(P)H dehydrogenase-like (NDH) complex is a large protein complex that reduces plastoquinone and pumps protons into the lumen generating protonmotive force. In plants, the complex consists of both nuclear and chloroplast-encoded subunits. Despite its perceived importance for stress tolerance and ATP generation, chloroplast-encoded NDH subunits have been lost numerous times during evolution in species occupying seemingly unrelated environmental niches. We have generated a phylogenetic tree that reveals independent losses in multiple phylogenetic lineages, and we use this tree as a reference to discuss possible evolutionary contexts that may have relaxed selective pressure for retention of ndh genes. While we are still yet unable to pinpoint a singular specific lifestyle that negates the need for NDH, we are able to rule out several long-standing explanations. In light of this, we discuss the biochemical changes that would be required for the chloroplast to dispense with NDH functionality with regards to known and proposed NDH-related reactions.

Comparative Plastid Genomics of Neotropical Bulbophyllum (Orchidaceae; Epidendroideae)

Pantropical Bulbophyllum, with ∼2,200 species, is one of the largest genera in Orchidaceae. Although phylogenetics and taxonomy of the ∼60 American species in the genus are generally well understood, some species complexes need more study to clearly delimit their component species and provide information about their evolutionary history. Previous research has suggested that the plastid genome includes phylogenetic markers capable of providing resolution at low taxonomic levels, and thus it could be an effective tool if these divergent regions can be identified. In this study, we sequenced the complete plastid genome of eight Bulbophyllum species, representing five of six Neotropical taxonomic sections. All plastomes conserve the typical quadripartite structure, and, although the general structure of plastid genomes is conserved, differences in ndh-gene composition and total length were detected. Total length was determined by contraction and expansion of the small single-copy region, a result of an independent loss of the seven ndh genes. Selection analyses indicated that protein-coding genes were generally well conserved, but in four genes, we identified 95 putative sites under positive selection. Furthermore, a total of 54 polymorphic simple sequence repeats were identified, for which we developed amplification primers. In addition, we propose 10 regions with potential to improve phylogenetic analyses of Neotropical Bulbophyllum species.

The complete chloroplast genome and characteristics analysis of Callistemon rigidus R.Br

Callistemon rigidus R.Br. one of the traditional Chinese medicinal plants, is acrid-flavored and mild-natured, with the prominent effects reducing swelling, resolving phlegm, and dispelling rheumatism. Clinically, it has been commonly used to treat cold, cough and asthma, pain and swelling from impact injuries, eczema, rheumatic arthralgia. The chloroplast genome study on Callistemon rigidus R.Br. is a few seen. This study demonstrates the data collected from the assembly and annotation of the chloroplast (cp) genome of Callistemon rigidus R.Br., followed by furthers comparative analysis with the cp genomes of closely related species. C. rigidus R.Br. showed a cp genome in the size of 158, 961 bp long with 36.78% GC content, among which a pair of inverted repeats (IRs) of 26, 671 bp separated a large single-copy (LSC) region of 87, 162 bp and a small single-copy (SSC) region of 18, 457 bp. Altogether 131 genes were hosted, including 37 transfer RNAs, 8 ribosomal RNAs, and 86 protein-coding genes. 284 simple sequence repeats (SSRs) were also marked out. A comparative analysis of the genome structure and the sequence data of closely related species unveiled the conserved gene order in the IR and LSC/SSC regions, a quite constructive finding for future phylogenetic research. Overall, this study providing C. rigidus R.Br. genomic resources could positively contribute to the evolutionary study and the phylogenetic reconstruction of Myrtaceae.

A comparative plastomics approach reveals available molecular markers for the phylogeographic study of Dendrobium huoshanense, an endangered orchid with extremely small populations

Comparative plastomics approaches have been used to identify available molecular markers for different taxonomic level studies of orchid species. However, the adoption of such methods has been largely limited in phylogeographic studies. Therefore, in this study, Dendrobium huoshanense, an endangered species with extremely small populations, was used as a model system to test whether the comparative plastomic approaches could screen available molecular markers for the phylogeographic study. We sequenced two more plastomes of D. huoshanense and compared them with our previously published one. A total of 27 mutational hotspot regions and six polymorphic cpSSRs have been screened for the phylogeographic studies of D. huoshanense. The cpDNA haplotype data revealed that the existence of haplotype distribution center was located in Dabieshan Mts. (Huoshan). The genetic diversity and phylogenetic analyses showed that the populations of D. huoshanense have been isolated and evolved independently for long period. On the contrary, based on cpSSR data, the genetic structure analysis revealed a mixed structure among the populations in Anhui and Jiangxi province, which suggested that the hybridization or introgression events have occurred among the populations of D. huoshanense. These results indicated that human activities have played key roles in shaping the genetic diversity and distributional patterns of D. huoshanense. According to our results, both two markers showed a high resolution for the phylogeographic studies of D. huoshanense. Therefore, we put forth that comparative plastomic approaches could revealed available molecular markers for phylogeographic study, especially for the species with extremely small populations.

Plastome Evolution and Phylogeny of Orchidaceae, With 24 New Sequences

In order to understand the evolution of the orchid plastome, we annotated and compared 124 complete plastomes of Orchidaceae representing all the major lineages in their structures, gene contents, gene rearrangements, and IR contractions/expansions. Forty-two of these plastomes were generated from the corresponding author's laboratory, and 24 plastomes-including nine genera (Amitostigma, Bulbophyllum, Dactylorhiza, Dipodium, Galearis, Gymnadenia, Hetaeria, Oreorchis, and Sedirea)-are new in this study. All orchid plastomes, except Aphyllorchis montana, Epipogium aphyllum, and Gastrodia elata, have a quadripartite structure consisting of a large single copy (LSC), two inverted repeats (IRs), and a small single copy (SSC) region. The IR region was completely lost in the A. montana and G. elata plastomes. The SSC is lost in the E. aphyllum plastome. The smallest plastome size was 19,047 bp, in E. roseum, and the largest plastome size was 178,131 bp, in Cypripedium formosanum. The small plastome sizes are primarily the result of gene losses associated with mycoheterotrophic habitats, while the large plastome sizes are due to the expansion of noncoding regions. The minimal number of common genes among orchid plastomes to maintain minimal plastome activity was 15, including the three subunits of rpl (14, 16, and 36), seven subunits of rps (2, 3, 4, 7, 8, 11, and 14), three subunits of rrn (5, 16, and 23), trnC-GCA, and clpP genes. Three stages of gene loss were observed among the orchid plastomes. The first was ndh gene loss, which is widespread in Apostasioideae, Vanilloideae, Cypripedioideae, and Epidendroideae, but rare in the Orchidoideae. The second stage was the loss of photosynthetic genes (atp, pet, psa, and psb) and rpo gene subunits, which are restricted to Aphyllorchis, Hetaeria, Hexalectris, and some species of Corallorhiza and Neottia. The third stage was gene loss related to prokaryotic gene expression (rpl, rps, trn, and others), which was observed in Epipogium, Gastrodia, Lecanorchis, and Rhizanthella. In addition, an intermediate stage between the second and third stage was observed in Cyrtosia (Vanilloideae). The majority of intron losses are associated with the loss of their corresponding genes. In some orchid taxa, however, introns have been lost in rpl16, rps16, and clpP(2) without their corresponding gene being lost. A total of 104 gene rearrangements were counted when comparing 116 orchid plastomes. Among them, many were concentrated near the IRa/b-SSC junction area. The plastome phylogeny of 124 orchid species confirmed the relationship of {Apostasioideae [Vanilloideae (Cypripedioideae (Orchidoideae, Epidendroideae))]} at the subfamily level and the phylogenetic relationships of 17 tribes were also established. Molecular clock analysis based on the whole plastome sequences suggested that Orchidaceae diverged from its sister family 99.2 mya, and the estimated divergence times of five subfamilies are as follows: Apostasioideae (79.91 mya), Vanilloideae (69.84 mya), Cypripedioideae (64.97 mya), Orchidoideae (59.16 mya), and Epidendroideae (59.16 mya). We also released the first nuclear ribosomal (nr) DNA unit (18S-ITS1-5.8S-ITS2-28S-NTS-ETS) sequences for the 42 species of Orchidaceae. Finally, the phylogenetic tree based on the nrDNA unit sequences is compared to the tree based on the 42 identical plastome sequences, and the differences between the two datasets are discussed in this paper.

The complete chloroplast genome of Magnolia polytepala: Comparative analyses offer implication for genetics and phylogeny of Yulania

Magnoliaceae is a primitive taxon in the angiosperms, comprising approximately 240 species in 2-17 genera. Many of them have been widely cultivated due to their horticultural and medicinal value. However, there are uncertainties and controversies about the delimitation of the genera except Liriodendron L. in this family. The Yulania taxa is also the focus of dispute at the genus and section levels. In this study, we compared ten Yulania plastomes, including the newly sequenced M. polytepala. The plastome-wide comparative analysis demonstrated that 1) Yulania cp genomes were highly conserved, and the majority differences existed in IR regions with the loss/retention of trnV-GAC or ycf15 gene, 2) mutational hotspots with high levels of nucleotide diversity (Pi > 0.02) existed in both coding (rpoA, and ycf1) and no-coding (ccsA-ndhD, ndhE-ndhG, ndhF-rpl32, petA-psbJ, rpl32-trnL, rps3-rps19, and trnH-psbA) regions among the genus Yulania. Combined with other data from Magnoliaceae plastomes, our reconstructed molecular phylogenetic tree revealed that Yulania is monophyletic, separated from the genus Magnolia L. (=Magnolia subg. Magnolia L.), but seems a sister of Michelia L. Moreover, M. polytepala which belongs to the genus Yulania is most closely related to M. liliiflora. All these results indicated that plastome data may contribute to investigating taxonomy, population genetics and phylogeny of Yulania.

Target Nuclear and Off-Target Plastid Hybrid Enrichment Data Inform a Range of Evolutionary Depths in the Orchid Genus Epidendrum

Universal angiosperm enrichment probe sets designed to enrich hundreds of putatively orthologous nuclear single-copy loci are increasingly being applied to infer phylogenetic relationships of different lineages of angiosperms at a range of evolutionary depths. Studies applying such probe sets have focused on testing the universality and performance of the target nuclear loci, but they have not taken advantage of off-target data from other genome compartments generated alongside the nuclear loci. Here we do so to infer phylogenetic relationships in the orchid genus Epidendrum and closely related genera of subtribe Laeliinae. Our aims are to: 1) test the technical viability of applying the plant anchored hybrid enrichment (AHE) method (Angiosperm v.1 probe kit) to our focal group, 2) mine plastid protein coding genes from off-target reads; and 3) evaluate the performance of the target nuclear and off-target plastid loci in resolving and supporting phylogenetic relationships along a range of taxonomical depths. Phylogenetic relationships were inferred from the nuclear data set through coalescent summary and site-based methods, whereas plastid loci were analyzed in a concatenated partitioned matrix under maximum likelihood. The usefulness of target and flanking non-target nuclear regions and plastid loci was assessed through the estimation of their phylogenetic informativeness. Our study successfully applied the plant AHE probe kit to Epidendrum, supporting the universality of this kit in angiosperms. Moreover, it demonstrated the feasibility of mining plastome loci from off-target reads generated with the Angiosperm v.1 probe kit to obtain additional, uniparentally inherited sequence data at no extra sequencing cost. Our analyses detected some strongly supported incongruences between nuclear and plastid data sets at shallow divergences, an indication of potential lineage sorting, hybridization, or introgression events in the group. Lastly, we found that the per site phylogenetic informativeness of the ycf1 plastid gene surpasses that of all other plastid genes and several nuclear loci, making it an excellent candidate for assessing phylogenetic relationships at medium to low taxonomic levels in orchids.

Characterization of sequence variability hotspots in Cranichideae plastomes (Orchidaceae, Orchidoideae)

This study reports complete plastome sequences for six species of Neotropical Cranichideae and focuses on identification of the most variable regions (hotspots) in this group of orchids. These structure of these six plastomes is relatively conserved, exhibiting lengths ranging between 142,599 to 154,562 bp with 36.7% GC on average and exhibiting typical quadripartite arrangement (LSC, SSC and two IRs). Variation detected in the LSC/IR and SSC/IR junctions is explained by the loss of ndhF and ycf1 length variation. For the two genera of epiphytic clade in Spiranthinae, almost whole sets of the ndh-gene family were missing. Eight mutation hotspots were identified based on nucleotide diversity, sequence variability and parsimony-informative sites. Three of them (rps16-trnQ, trnT-trnL, rpl32-trnL) seem to be universal hotspots in the family, and the other five (trnG-trnR, trnR-atpA, trnP-psaJ, rpl32-infA, and rps15-ycf1) are described for the first time as orchid molecular hotspots. These regions have much more variation than all those used previously in phylogenetics of the group and offer useful plastid markers for phylogenetic, barcoding and population genetic studies. The use of whole plastomes or exclusive no-gap matrices also positioned with high support the holomycotrophic Rhizanthella among Orchidoideae plastomes in model-based analyses, showing the utility of plastomes for phylogenetic placement of this unusual genus.

Complete Chloroplast Genome of Paphiopedilum delenatii and Phylogenetic Relationships among Orchidaceae

Paphiopedilum delenatii is a native orchid of Vietnam with highly attractive floral traits. Unfortunately, it is now listed as a critically endangered species with a few hundred individuals remaining in nature. In this study, we performed next-generation sequencing of P. delenatii and assembled its complete chloroplast genome. The whole chloroplast genome of P. delenatii was 160,955 bp in size, 35.6% of which was GC content, and exhibited typical quadripartite structure of plastid genomes with four distinct regions, including the large and small single-copy regions and a pair of inverted repeat regions. There were, in total, 130 genes annotated in the genome: 77 coding genes, 39 tRNA genes, 8 rRNA genes, and 6 pseudogenes. The loss of ndh genes and variation in inverted repeat (IR) boundaries as well as data of simple sequence repeats (SSRs) and divergent hotspots provided useful information for identification applications and phylogenetic studies of Paphiopedilum species. Whole chloroplast genomes could be used as an effective super barcode for species identification or for developing other identification markers, which subsequently serves the conservation of Paphiopedilum species.

Comparison and Phylogenetic Analysis of Chloroplast Genomes of Three Medicinal and Edible Amomum Species

Amomum villosum is an important medicinal and edible plant with several pharmacologically active volatile oils. However, identifying A. villosum from A. villosum var. xanthioides and A. longiligulare which exhibit similar morphological characteristics to A. villosum, is difficult. The main goal of this study, therefore, is to mine genetic resources and improve molecular methods that could be used to distinguish these species. A total of eight complete chloroplasts (cp) genomes of these Amomum species which were collected from the main producing areas in China were determined to be 163,608–164,069 bp in size. All genomes displayed a typical quadripartite structure with a pair of inverted repeat (IR) regions (29,820–29,959 bp) that separated a large single copy (LSC) region (88,680–88,857 bp) from a small single copy (SSC) region (15,288–15,369 bp). Each genome encodes 113 different genes with 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. More than 150 SSRs were identified in the entire cp genomes of these three species. The Sanger sequencing results based on 32 Amomum samples indicated that five highly divergent regions screened from cp genomes could not be used to distinguish Amomum species. Phylogenetic analysis showed that the cp genomes could not only accurately identify Amomum species, but also provide a solid foundation for the establishment of phylogenetic relationships of Amomum species. The availability of cp genome resources and the comparative analysis is beneficial for species authentication and phylogenetic analysis in Amomum.