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.

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

Three complete plastome sequences from the families of Lamiaceae, Mazaceae, and Phrymaceae (Lamiales).

In this study, we report the first complete plastome sequence of Vitex rotundifolia (Lamiaceae) (MT937186). In addition, the plastome sequences of Phryma leptostachya subsp. asiatica (Phrymaceae) (153,324 bp; MT948145) and Mazus pumilus (Mazaceae) (152,847 bp; MT937187) are also included. The gene orders and structures of the three plastomes are collinear with those of the typical plastome of angiosperm. The plastome size of V. rotundifolia is 154,370 bp in length and consists of a large single-copy region of 85,079 bp and a small single-copy region of 17,917 bp, which are separated by a pair of 25,687 bp-long inverted repeat regions. In addition, the plastome sizes of P. leptostachya subsp. asiatica and M. pumilus are 153,324 bp and 152,847 bp, respectively. The three plastomes contain 113 genes, including 79 protein-coding, 30 tRNA, and four rRNA genes. Sixteen genes contain one intron and two genes have two introns. A total of 41 simple sequence repeat loci was identified in the V. rotundifolia plastome. Phylogenetic analysis shows that Viticoideae is a sister group of the last of Lamiaceae except Nepetoideae. The Mazaceae are a sister group of Lamiaceae, while Phrymaceae form a sister group to the Paulowniaceae-Orobanchaceae clade.

Corrigendum: Plastome Evolution and Phylogeny of Orchidaceae, With 24 New Sequences.

[This corrects the article DOI: 10.3389/fpls.2020.00022.].

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.

Extensive Losses of Photosynthesis Genes in the Plastome of a Mycoheterotrophic Orchid, Cyrtosia septentrionalis (Vanilloideae: Orchidaceae).

Cyrtosia septentrionalis is an achlorophyllous mycoheterotrophic orchid in the subfamily Vanilloideae (Orchidaceae). This article reports C. septentrionalis's complete plastome sequence and compare it with other orchid plastomes with a same mycoheterotrophic nutritional mode. The C. septentrionalis plastome has decreased to 96,859 bp in length, but it still maintains a quadripartite structure. The C. septentrionalis plastome contains 38 protein-coding genes, 25 tRNA genes, and four ribosomal RNA genes. Most genes related to photosynthesis have been lost, whereas the majority of housekeeping genes remain; this pattern corresponds to the end of stage 3 gene degradation. The inverted repeat regions of the C. septentrionalis plastome have decreased to 10,414 bp and mainly contain the gene ycf2. A block consisting of four rrn genes and rps7 and rps12 has shifted to a small single-copy region. As a result, the small single-copy region was found to be expanded, despite the loss of all ndh genes in the region. Three inversion mutations are required to explain the C. septentrionalis plastome's current gene order. The species is endangered, and these results have implications for its conservation.

The complete plastome sequence of the endangered orchid Kuhlhasseltia nakaiana (Orchidaceae).

In this study, we report the complete plastome sequence of Kuhlhasseltia nakaiana (F.Maek.) Ormerod (Orchidaceae) (NCBI acc. no. KY354041), an endangered plant species protected by the national law of Korea. The gene order and number in the K. nakaiana plastome were similar to a typical orchid plastome. The complete plastome was 147,614 bp in length and consisted of a large single copy region of 81,617 bp and a small single copy region of 13,673 bp, separated by two inverted repeats of 26,162 bp. The plastome contained 103 genes, of which 69 were protein-coding genes, 30 were tRNA genes, and four were rRNA genes. Fourteen genes contained one intron and two genes (clpP and ycf3) had two introns. The AT content of the plastome was 60.5%. A total of 74 simple sequence repeat regions were identified from the plastome. Phylogenetic analysis determined that K. nakaiana was a member of the tribe Cranichideae and revealed the sister group relationship between K. nakaiana and Ludisia discolor within the tribe Cranichideae.

The complete plastome sequence of the endangered orchid Habenaria radiata (Orchidaceae).

In this study, we determined the complete chloroplast sequence of Habenaria radiata (Thunb.) Spreng. (Orchidaceae) (NCBI acc. no. KX871237), an endangered plant species protected by the national law of Korea. The gene order and gene content of the H. radiata plastome are similar to those of typical angiosperm plastomes. The 11 ndh genes, which are usually lost in orchid plastomes, are intact in the H. radiata plastome. The complete plastome is 155,353 bp in length and consists of a large single copy of 84,833 bp and a small single copy of 17,718 bp, separated by two inverted repeats of 26,401 bp. The plastome contains 113 genes, of which 79 are protein-coding genes, 30 are tRNA genes, and four are rRNA genes. Sixteen genes contain one intron and two genes (clpP, ycf3) have two introns. A total of 76 simple sequence repeat (SSR) loci, which consist of 58 mono-SSR, 17 di-SSR, and 1 tri-SSR, are scattered along the H. radiata plastome. Some of these plastome SSR and high sequence divergent regions may be useful for development of genetic markers for the populations of H. radiata and other congeneric species. Phylogenetic analysis identified the sister relationship between H. radiata and H. pantlingiana within the tribe Orchideae.