Orchidaceae is the 3rd largest family of angiosperms, an evolved young branch of monocotyledons. This family contains a number of economically-important horticulture and flowering plants. However, the limited availability of genomic information largely hindered the study of molecular evolution and phylogeny of Orchidaceae. In this study, we determined the evolutionary characteristics of whole chloroplast (cp) genomes and the phylogenetic relationships of the family Orchidaceae. We firstly characterized the cp genomes of four orchid species: Cremastra appendiculata, Calanthe davidii, Epipactis mairei, and Platanthera japonica. The size of the chloroplast genome ranged from 153,629 bp (C. davidi) to 160,427 bp (E. mairei). The gene order, GC content, and gene compositions are similar to those of other previously-reported angiosperms. We identified that the genes of ndhC, ndhI, and ndhK were lost in C. appendiculata, in that the ndh I gene was lost in P. japonica and E. mairei. In addition, the four types of repeats (forward, palindromic, reverse, and complement repeats) were examined in orchid species. E. mairei had the highest number of repeats (81), while C. davidii had the lowest number (57). The total number of Simple Sequence Repeats is at least 50 in C. davidii, and, at most, 78 in P. japonica. Interestingly, we identified 16 genes with positive selection sites (the psbH, petD, petL, rpl22, rpl32, rpoC1, rpoC2, rps12, rps15, rps16, accD, ccsA, rbcL, ycf1, ycf2, and ycf4 genes), which might play an important role in the orchid species’ adaptation to diverse environments. Additionally, 11 mutational hotspot regions were determined, including five non-coding regions (ndhB intron, ccsA-ndhD, rpl33-rps18, ndhE-ndhG, and ndhF-rpl32) and six coding regions (rps16, ndhC, rpl32, ndhI, ndhK, and ndhF). The phylogenetic analysis based on whole cp genomes showed that C. appendiculata was closely related to C. striata var. vreelandii, while C. davidii and C. triplicate formed a small monophyletic evolutionary clade with a high bootstrap support. In addition, five subfamilies of Orchidaceae, Apostasioideae, Cypripedioideae, Epidendroideae, Orchidoideae, and Vanilloideae, formed a nested evolutionary relationship in the phylogenetic tree. These results provide important insights into the adaptive evolution and phylogeny of Orchidaceae.

Pinus L. is the largest genus of conifers and provides a classical model for studying species divergence and phylogenetic evolution by gymnosperms. However, our poor understanding of sequence divergence in the whole plastid genomes of Pinus species severely hinders studies of their evolution and phylogeny. Thus, we analyzed the sequences of 97 Pinus plastid genomes, including four newly sequenced genomes and 93 previously published plastomes, to explore the evolution and phylogenetic relationships in the genus Pinus. The complete chloroplast genomes of Pinus species ranged in size from 109 640 bp (P. cembra L.) to 121 976 bp (P. glabra Walter), and these genomes comprised circular DNA molecules in a similar manner to those of most gymnosperms. We identified 9108 repeats where most of the repeats comprised the dispersed type with 3983 (44%), followed by tandem repeats with 2999 (33%), and then palindromic repeats with 2126 (23%). Sixteen divergence hotspot regions were identified in Pinus plastid genomes, which could be useful molecular markers for future population genetics studies. Phylogenetic analysis showed that Pinus species could be divided into two diverged clades comprising the subgenera Strobus (single needle section) and Pinus (double needles section). Molecular dating suggested that the genus Pinus originated approximately 130.38 Mya during the late Cretaceous. The two subgenera subsequently split 85.86 Mya, which was largely consistent with the other molecular results based on partial DNA markers. These findings provide important insights into the sequence variations and phylogenetic evolution of Pinus plastid genomes.

The complete nucleotide sequence of the Taiwan red pine Pinus taiwanensis Hayata chloroplast genome (cpDNA) is determined in this study. The genome is composed of 119,741 bp in length, containing a pair of very short inverted repeat (IRa and IRb) regions of 495 bp, which was divided by a large single-copy (LSC) region of 65,670 bp and a small single-copy (SSC) region of 53,080 bp in length. The cpDNA contained 115 genes, including 74 protein-coding genes (73 PCG species), 4 ribosomal RNA genes (four rRNA species) and 37 tRNA genes (22 tRNA species). Out of these genes, 12 harbored a single intron, and one (rps12) contained a couple of introns. The overall AT content of the Taiwan red pine cpDNA is 61.5%, while the corresponding values of the LSC, SSC and IR regions are 62.2%, 60.6% and 63.6%, respectively. A maximum parsimony phylogenetic analysis suggested that the genus Pinus, Picea, Abies and Larix were strongly supported as monophyletic, and the cpDNA of P. taiwanensis is closely related to that of P. thunbergii.