Complete plastid genome sequences of Drimys, Liriodendron, and Piper: implications for the phylogenetic relationships of magnoliids

The adaptive evolution of Quercus section Ilex using the chloroplast genomes of two threatened species

Chloroplast (cp) genome sequences have been extensively used for phylogenetic and evolutionary analyses, as many have been sequenced in recent years. Identification of Quercus is challenging because many species overlap phenotypically owing to interspecific hybridization, introgression, and incomplete lineage sorting. Therefore, we wanted to gain a better understanding of this genus at the level of the maternally inherited chloroplast genome. Here, we sequenced, assembled, and annotated the cp genomes of the threatened Quercus marlipoensis (160,995 bp) and Q. kingiana (161,167 bp), and mined these genomes for repeat sequences and codon usage bias. Comparative genomic analyses, phylogenomics, and selection pressure analysis were also performed in these two threatened species along with other species of Quercus. We found that the guanine and cytosine content of the two cp genomes were similar. All 131 annotated genes, including 86 protein-coding genes, 37 transfer RNA genes, and 8 ribosomal RNA genes, had the same order in the two species. A strong A/T bias was detected in the base composition of simple sequence repeats. Among the 59 synonymous codons, the codon usage pattern of the cp genomes in these two species was more inclined toward the A/U ending. Comparative genomic analyses indicated that the cp genomes of Quercus section Ilex are highly conserved. We detected eight highly variable regions that could be used as molecular markers for species identification. The cp genome structure was consistent and different within and among the sections of Quercus. The phylogenetic analysis showed that section Ilex was not monophyletic and was divided into two groups, which were respectively nested with section Cerris and section Cyclobalanopsis. The two threatened species sequenced in this study were grouped into the section Cyclobalanopsis. In conclusion, the analyses of cp genomes of Q. marlipoensis and Q. kingiana promote further study of the taxonomy, phylogeny and evolution of these two threatened species and Quercus.

Plastome structure and phylogenetic relationships of genus Hydrocotyle (apiales): provide insights into the plastome evolution of Hydrocotyle

BACKGROUND

The genus Hydrocotyle Tourn. ex L. is a key group for further study on the evolution of Apiales, comprising around 170 species globally. Previous studies mainly focused on separate sections and provided much information about this genus, but its infrageneric relationships are still confusing. In addition, the genetic basis of its adaptive evolution remains poorly understood. To investigate the phylogeny and evolution of the genus, we selected ten representative species covering two of three diversity distribution centers and exhibiting rich morphology diversity. Comparative plastome analysis was conducted to clarify the structural character of Hydrocotyle plastomes. Positive selection analyses were implemented to assess the evolution of the genus. Phylogenetic inferences with protein-coding sequences (CDS) of Hydrocotyle and 17 related species were also performed.

RESULTS

Plastomes within Hydrocotyle were generally conservative in structure, gene order, and size. A total of 14 regions (rps16-trnK, trnQ-rps16, atpI-atpH, trnC-petN-psbM, ycf3-trnS, accD-psaI-ycf4, petA-psbJ, rps12-rpl20, rpl16 intron, rps3-rpl16 intron, rps9-rpl22, ndhF-rpl32, ndhA intron, and ycf1a) were recognized as hotspot regions within the genus, which suggested to be promising DNA barcodes for global phylogenetic analysis of Hydrocotyle. The ycf15 gene was suggested to be a protein-coding gene for Hydrocotyle species, and it could be used as a DNA barcode to identify Hydrocotyle. In phylogenetic analysis, three monophyletic clades (Clade I, II, III) were identified with evidence of rapid radiation speciation within Clade I. The selective pressure analysis detected that six CDS genes (ycf1b, matK, atpF, accD, rps14, and psbB) of Hydrocotyle species were under positive selection. Within the genus, the last four genes were conservative, suggesting a relation to the unique evolution of the genus in Apiales. Seven genes (atpE, matK, psbH, ycf1a, ycf1b, rpoA, and ycf2) were detected to be under some degree of positive selection in different taxa within the genus Hydrocotyle, indicating their role in the adaptive evolution of species.

CONCLUSIONS

Our study offers new insights into the phylogeny and adaptive evolution of Hydrocotyle. The plastome sequences could significantly enhance phylogenetic resolution and provide genomic resources and potential DNA markers useful for future studies of the genus.

Complete chloroplast genome of Asarum chungbuensis (C.S. Yook & J.G. Kim) B.U. Oh 2005 (Aristolochiaceae), a Korean endemic species

Asarum chungbuensis, a species endemic to Korea, has a limited distribution across the Korean Peninsula and is used in traditional medicine. Despite its importance, the genome structure, genetic composition, and phylogenetic relationships based on its chloroplast genome have not been documented. In this study, the complete chloroplast genome of A. chungbuensis was newly assembled. The chloroplast genome is 190,179 base pairs (bp) long, and the overall GC content (%) of the plastid was 36.8%. The chloroplast genome size of A. chungbuensis is longer than that of the normal chloroplast genome (160 kb) because of an inverted small single-copy (SSC) duplication that incorporates the SSC into an inverted repeat (IR) region. By extension, this duplication event causes this chloroplast genome to lack an SSC, unlike other formal structures. The chloroplast genome, with a tripartite structure, consisted of a single-copy region of 93,351 bp with a 34.6% GC content and two IR regions, each with a length of 48,414 bp and a 38.8% GC content. Additionally, it was found to have 113 genes, including 79 PCG genes, 30 tRNA genes, and four rRNA genes. Phylogenetic analysis revealed that A. chungbuensis was grouped with A. heterotropoides var. seoulense, which diverged from the clade comprising A. koreanum and A. patens. The newly sequenced A. chungbuensis chloroplast genome could provide valuable genomic information for determining unique genome structures, especially for the assessment of genetic diversity, phylogenetic relationships, species conservation, and biogeographic studies of the genus Asarum.

Phylogenomic analysis and dynamic evolution of chloroplast genomes of Clematis nannophylla

Clematis nannophylla is a perennial shrub of Clematis with ecological, ornamental, and medicinal value, distributed in the arid and semi-arid areas of northwest China. This study successfully determined the chloroplast (cp) genome of C. nannophylla, reconstructing a phylogenetic tree of Clematis. This cp genome is 159,801 bp in length and has a typical tetrad structure, including a large single-copy, a small single-copy, and a pair of reverse repeats (IRa and IRb). It contains 133 unique genes, including 89 protein-coding, 36 tRNA, and 8 rRNA genes. Additionally, 66 simple repeat sequences, 50 dispersed repeats, and 24 tandem repeats were found; many of the dispersed and tandem repeats were between 20-30 bp and 10-20 bp, respectively, and the abundant repeats were located in the large single copy region. The cp genome was relatively conserved, especially in the IR region, where no inversion or rearrangement was observed, further revealing that the coding regions were more conserved than the noncoding regions. Phylogenetic analysis showed that C. nannophylla is more closely related to C. fruticosa and C. songorica. Our analysis provides reference data for molecular marker development, phylogenetic analysis, population studies, and cp genome processes to better utilise C. nannophylla.

Characterization and Comparative Analysis of Complete Chloroplast Genomes of Four Bromus (Poaceae, Bromeae) Species

Bromus (Poaceae Bromeae) is a forage grass with high adaptability and ecological and economic value. Here, we sequenced Bromus ciliatus, Bromus benekenii, Bromus riparius, and Bromus rubens chloroplast genomes and compared them with four previously described species. The genome sizes of Bromus species ranged from 136,934 bp (Bromus vulgaris) to 137,189 bp (Bromus ciliates, Bromus biebersteinii), with a typical quadripartite structure. The studied species had 129 genes, consisting of 83 protein-coding, 38 tRNA-coding, and 8 rRNA-coding genes. The highest GC content was found in the inverted repeat (IR) region (43.85-44.15%), followed by the large single-copy (LSC) region (36.25-36.65%) and the small single-copy (SSC) region (32.21-32.46%). There were 33 high-frequency codons, with those ending in A/U accounting for 90.91%. A total of 350 simple sequence repeats (SSRs) were identified, with single-nucleotide repeats being the most common (61.43%). A total of 228 forward and 141 palindromic repeats were identified. No reverse or complementary repeats were detected. The sequence identities of all sequences were very similar, especially with respect to the protein-coding and inverted repeat regions. Seven highly variable regions were detected, which could be used for molecular marker development. The constructed phylogenetic tree indicates that Bromus is a monophyletic taxon closely related to Triticum. This comparative analysis of the chloroplast genome of Bromus provides a scientific basis for species identification and phylogenetic studies.

The complete chloroplast genome of Oxytropis ochrocephala Bunge 1874 (Fabaceae) and its phylogenetic analysis

This study presents the first-ever complete chloroplast (cp) genome sequence of Oxytropis ochrocephala Bunge 1874, a member of the Fabaceae family. The cp genome spans 126,996 base pairs and includes 109 genes, comprising 76 protein-coding genes, 29 tRNA genes, and four rRNA genes. Notably, the genome lacks an inverted repeat (IR) region. Additionally, we constructed phylogenetic trees for 34 species within Trib. Galegeae, employing both maximum likelihood (ML) and Bayesian inference (BI) methods. These analyses robustly support the monophyly of the Oxytropis species, evidenced by high bootstrap values (BP = 100) and posterior probabilities (PP = 1). Within this clade, O. ochrocephala exhibits a sister relationship with other Oxytropis species. Our findings provide valuable insights into the genetic makeup and evolutionary relationships of O. ochrocephala within the Galegeae tribe.

A comparative analysis of plastome evolution in autotrophic Piperales

PREMISE

Many plastomes of autotrophic Piperales have been reported to date, describing a variety of differences. Most studies focused only on a few species or a single genus, and extensive, comparative analyses have not been done. Here, we reviewed publicly available plastome reconstructions for autotrophic Piperales, reanalyzed publicly available raw data, and provided new sequence data for all previously missing genera. Comparative plastome genomics of >100 autotrophic Piperales were performed.

METHODS

We performed de novo assemblies to reconstruct the plastomes of newly generated sequence data. We used Sanger sequencing and read mapping to verify the assemblies and to bridge assembly gaps. Furthermore, we reconstructed the phylogenetic relationships as a foundation for comparative plastome genomics.

RESULTS

We identified a plethora of assembly and annotation issues in published plastome data, which, if unattended, will lead to an artificial increase of diversity. We were able to detect patterns of missing and incorrect feature annotation and determined that the inverted repeat (IR) boundaries were the major source for erroneous assembly. Accounting for the aforementioned issues, we discovered relatively stable junctions of the IRs and the small single-copy region (SSC), whereas the majority of plastome variations among Piperales stems from fluctuations of the boundaries of the IR and the large single-copy (LSC) region.

CONCLUSIONS

This study of all available plastomes of autotrophic Piperales, expanded by new data for previously missing genera, highlights the IR-LSC junctions as a potential marker for discrimination of various taxonomic levels. Our data indicates a pseudogene-like status for cemA and ycf15 in various Piperales. Based on a review of published data, we conclude that incorrect IR-SSC boundary identification is the major source for erroneous plastome assembly. We propose a gold standard for assembly and annotation of high-quality plastomes based on de novo assembly methods and appropriate references for gene annotation.

Complete Chloroplast Genomes of Four Oaks from the Section Cyclobalanopsis Improve the Phylogenetic Analysis and Understanding of Evolutionary Processes in the Genus Quercus

Quercus is a valuable genus ecologically, economically, and culturally. They are keystone species in many ecosystems. Species delimitation and phylogenetic studies of this genus are difficult owing to frequent hybridization. With an increasing number of genetic resources, we will gain a deeper understanding of this genus. In the present study, we collected four Quercus section Cyclobalanopsis species (Q. poilanei, Q. helferiana, Q. camusiae, and Q. semiserrata) distributed in Southeast Asia and sequenced their complete genomes. Following analysis, we compared the results with those of other species in the genus Quercus. These four chloroplast genomes ranged from 160,784 bp (Q. poilanei) to 161,632 bp (Q. camusiae) in length, with an overall guanine and cytosine (GC) content of 36.9%. Their chloroplast genomic organization and order, as well as their GC content, were similar to those of other Quercus species. We identified seven regions with relatively high variability (rps16, ndhk, accD, ycf1, psbZ-trnG-GCC, rbcL-accD, and rpl32-trnL-UAG) which could potentially serve as plastid markers for further taxonomic and phylogenetic studies within Quercus. Our phylogenetic tree supported the idea that the genus Quercus forms two well-differentiated lineages (corresponding to the subgenera Quercus and Cerris). Of the three sections in the subgenus Cerris, the section Ilex was split into two clusters, each nested in the other two sections. Moreover, Q. camusiae and Q. semiserrata detected in this study diverged first in the section Cyclobalanopsis and mixed with Q. engleriana in the section Ilex. In particular, 11 protein coding genes (atpF, ndhA, ndhD, ndhF, ndhK, petB, petD, rbcL, rpl22, ycf1, and ycf3) were subjected to positive selection pressure. Overall, this study enriches the chloroplast genome resources of Quercus, which will facilitate further analyses of phylogenetic relationships in this ecologically important tree genus.

The complete chloroplast genome of Hedychium flavum Roxb., an ornamental, edible and medicinal plant

Hedychium flavum Roxb. 1820 is a perennial herb mainly distributed in China, India, Myanmar and Thailand with ornamental, edible and medicinal value. It is extensively cultivated as a source of aromatic essential oils, ornamental plant, food flavorings and vegetables, and folk medicine. In this study, we sequence the complete chloroplast genome of H. flavum by de novo assembly. The assembled genome has a typical quadripartite circular structure with 163,909 bp in length, containing a large single-copy region (LSC, 88,589 bp), a small single-copy region (SSC, 15,762 bp), and two inverted repeat regions (IRs, 29,779 bp). The cp genome contains 133 genes, including 87 protein-coding genes, 38 tRNA genes and 8 rRNA genes. Phylogenetic analysis based on the complete cp genome shows a close affinity of H. flavum and H. neocarneum with 100% bootstrap support. This study will provide useful genetic resource for further phylogenetic analysis of the genus Hedychium and Zingiberaceae.

The complete chloroplast genome of Chloranthus nervosus Collett ex Hemsl. 1890 (Chloranthaceae)

In this study, we sequenced and assembled the complete chloroplast genome of Chloranthus nervosus Collett ex Hemsl. 1890. The total length of the complete chloroplast sequence was found to be 158,002 bp. It consisted of a large single-copy (LSC) region of 87,127 bp, a small single-copy (SSC) region of 18,541 bp, and a pair of inverted repeat (IR) regions, each with a length of 26,167 bp. The overall GC content of the complete chloroplast genome was 38.9%, with the LSC region, SSC region, and IR regions exhibiting GC contents of 37.4%, 34.1%, and 43.1%, respectively. The annotation of the chloroplast genome revealed a total of 131 genes, comprising 86 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Phylogenetic analysis revealed that the seven sampled species of Chloranthus were divided into two clades. Within the clade characterized by long filamentous anther connectives, C. nervosus showed the closest relation to C. japonicus. These findings validated the previous preliminary results on the phylogenetic relationships of the seven species of Chloranthus with strong support.

Comparison of chloroplast genomes and phylogenetic analysis of four species in Quercus section Cyclobalanopsis

The identification in Quercus L. species was considered to be difficult all the time. The fundamental phylogenies of Quercus have already been discussed by morphological and molecular means. However, the morphological characteristics of some Quercus groups may not be consistent with the molecular results (such as the group Helferiana), which may lead to blurring of species relationships and prevent further evolutionary researches. To understand the interspecific relationships and phylogenetic positions, we sequenced and assembled the CPGs (160,715 bp-160842 bp) of four Quercus section Cyclobalanopsis species by Illumina pair-end sequencing. The genomic structure, GC content, and IR/SC boundaries exhibited significant conservatism. Six highly variable hotspots were detected in comparison analysis, among which rpoC1, clpP and ycf1 could be used as molecular markers. Besides, two genes (petA, ycf2) were detected to be under positive selection pressure. The phylogenetic analysis showed: Trigonobalanus genus and Fagus genus located at the base of the phylogeny tree; The Quercus genus species were distincted to two clades, including five sections. All Compound Trichome Base species clustered into a single branch, which was in accordance with the results of the morphological studies. But neither of group Gilva nor group Helferiana had formed a monophyly. Six Compound Trichome Base species gathered together in pairs to form three branch respectively (Quercus kerrii and Quercus chungii; Quercus austrocochinchinensis with Quercus gilva; Quercus helferiana and Quercus rex). Due to a low support rate (0.338) in the phylogeny tree, the interspecies relationship between the two branches differentiated by this node remained unclear. We believe that Q. helferiana and Q. kerrii can exist as independent species due to their distance in the phylogeny tree. Our study provided genetic information in Quercus genus, which could be applied to further studies in taxonomy and phylogenetics.

New insights into the plastome evolution of Lauraceae using herbariomics

BACKGROUND

The family Lauraceae possesses ca. 50 genera and 2,500-3,000 species that are distributed in the pantropics. Only half of the genera of the family were represented in previously published plastome phylogenies because of the difficulty of obtaining research materials. Plastomes of Hypodaphnideae and the Mezilaurus group, two lineages with unusual phylogenetic positions, have not been previously reported and thus limit our full understanding on the plastome evolution of the family. Herbariomics, promoted by next generation sequencing technology, can make full use of herbarium specimens, and provides opportunities to fill the sampling gap.

RESULTS

In this study, we sequenced five new plastomes (including four genera which are reported for the first time, viz. Chlorocardium, Hypodaphnis, Licaria and Sextonia) from herbarium specimens using genome skimming to conduct a comprehensive analysis of plastome evolution of Lauraceae as a means of sampling representatives of all major clades of the family. We identified and recognized six types of plastomes and revealed that at least two independent loss events at the IR-LSC boundary and an independent expansion of SSC occurred in the plastome evolution of the family. Hypodaphnis possesses the ancestral type of Lauraceae with trnI-CAU, rpl23 and rpl2 duplicated in the IR regions (Type-I). The Mezilaurus group shares the same plastome structure with the core Lauraceae group in the loss of trnI-CAU, rpl23 and rpl2 in the IRa region (Type-III). Two new types were identified in the Ocotea group: (1) the insertion of trnI-CAU between trnL-UAG and ccsA in the SSC region of Licaria capitata and Ocotea bracteosa (Type-IV), and (2) trnI-CAU and pseudogenizated rpl23 inserted in the same region of Nectandra angustifolia (Type-V). Our phylogeny suggests that Lauraceae are divided into nine major clades largely in accordance with the plastome types. The Hypodaphnideae are the earliest diverged lineage supported by both robust phylogeny and the ancestral plastome type. The monophyletic Mezilaurus group is sister to the core Lauraceae.

CONCLUSIONS

By using herbariomics, we built a more complete picture of plastome evolution and phylogeny of the family, thus providing a convincing case for further use of herbariomics in phylogenetic studies of the Lauraceae.

Phylogenetic and comparative analyses of Hydnora abyssinica plastomes provide evidence for hidden diversity within Hydnoraceae

BACKGROUND

To date, plastid genomes have been published for all but two holoparasitic angiosperm families. However, only a single or a few plastomes represent most of these families. Of the approximately 40 genera of holoparasitic angiosperms, a complete plastid genome sequence is available for only about half. In addition, less than 15 species are currently represented with more than one published plastid genome, most of which belong to the Orobanchaceae. Therefore, a significant portion of the holoparasitic plant plastome diversity remains unexplored. This limited information could hinder potential evolutionary pattern recognition as well as the exploration of inter- and intra-species plastid genome diversity in the most extreme holoparasitic angiosperms.

RESULTS

Here, we report the first plastomes of Kenyan Hydnora abyssinica accessions. The plastomes have a typical quadripartite structure and encode 24 unique genes. Phylogenetic tree reconstruction recovers the Kenyan accessions as monophyletic and together in a clade with the Namibian H. abyssinica accession and the recently published H. arabica from Oman. Hydnora abyssinica as a whole however is recovered as non-monophyletic, with H. arabica nested within. This result is supported by distinct structural plastome synapomorphies as well as pairwise distance estimates that reveal hidden diversity within the Hydnora species in Africa.

CONCLUSION

We propose to increase efforts to sample widespread holoparasitic species for their plastid genomes, as is the case with H. abyssinica, which is widely distributed in Africa. Morphological reinvestigation and further molecular data are needed to fully investigate the diversity of H. abyssinica along the entire range of distribution, as well as the diversity of currently synonymized taxa.

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.

Insights into phylogenetic relationships in Pinus inferred from a comparative analysis of complete chloroplast genomes

BACKGROUND

Pinus is the largest genus of Pinaceae and the most primitive group of modern genera. Pines have become the focus of many molecular evolution studies because of their wide use and ecological significance. However, due to the lack of complete chloroplast genome data, the evolutionary relationship and classification of pines are still controversial. With the development of new generation sequencing technology, sequence data of pines are becoming abundant. Here, we systematically analyzed and summarized the chloroplast genomes of 33 published pine species.

RESULTS

Generally, pines chloroplast genome structure showed strong conservation and high similarity. The chloroplast genome length ranged from 114,082 to 121,530 bp with similar positions and arrangements of all genes, while the GC content ranged from 38.45 to 39.00%. Reverse repeats showed a shrinking evolutionary trend, with IRa/IRb length ranging from 267 to 495 bp. A total of 3,205 microsatellite sequences and 5,436 repeats were detected in the studied species chloroplasts. Additionally, two hypervariable regions were assessed, providing potential molecular markers for future phylogenetic studies and population genetics. Through the phylogenetic analysis of complete chloroplast genomes, we offered novel opinions on the genus traditional evolutionary theory and classification.

CONCLUSION

We compared and analyzed the chloroplast genomes of 33 pine species, verified the traditional evolutionary theory and classification, and reclassified some controversial species classification. This study is helpful in analyzing the evolution, genetic structure, and the development of chloroplast DNA markers in Pinus.

Comparitive Analysis of the Chloroplast Genomes of Three Houpoea Plants

The genus Houpoea belongs to the family Magnoliaceae, and the species in this genus have important medicinal values. However, the investigation of the correlation between the evolution of the genus and its phylogeny has been severely hampered by the unknown range of species within the genus and the paucity of research on its chloroplast genome. Thus, we selected three species of Houpoea: Houpoea officinalis var officinalis (OO), Houpoea officinalis var. biloba (OB), and Houpoea rostrata (R). With lengths of 160,153 bp (OO), 160,011 bp (OB), and 160,070 bp (R), respectively, the whole chloroplast genomes (CPGs) of these three Houpoea plants were acquired via Illumina sequencing technology, and the findings were annotated and evaluated. These three chloroplast genomes were revealed by the annotation findings to be typical tetrads. A total of 131, 132, and 120 different genes were annotated. The CPGs of the three species had 52, 47, and 56 repeat sequences, which were primarily found in the ycf2 gene. A useful tool for identifying species is the approximately 170 simple sequence repeats (SSRs) that have been found. The border area of the reverse repetition region (IR) was studied, and it was shown that across the three Houpoea plants, it is highly conservative, with only changes between H. rostrata and the other two plants observed. Numerous highly variable areas (rps3-rps19, rpl32-trnL, ycf1, ccsA, etc.) have the potential to serve as the barcode label for Houpoea, according to an examination of mVISTA and nucleotide diversity (Pi). Phylogenetic relation indicates that Houpoea is a monophyletic taxon, and its genus range and systematic position are consistent with the Magnoliaceae system of Sima Yongkang-Lu Shugang, including five species and varieties of H. officinalis var. officinalis, H. rostrata, H. officinalis var. biloba, Houpoea obovate, and Houpoea tripetala, which evolved and differentiated from the ancestors of Houpoea to the present Houpoea in the above order. This study provides valuable information on the genus Houpoea, enriches the CPG information on Houpoea genus, and provides genetic resources for the further classification of and phylogenetic research on Houpoea.

Phylogenetic studies of magnoliids: Advances and perspectives

Magnoliids are the largest flowering plant clades outside of the eudicots and monocots, which are distributed worldwide and have high economic, ornamental and ecological values. Eudicots, monocots and magnoliids are the three major clades of Mesangiospermae, and their phylogenetic relationship is one of the most interesting issues. In recent years, with the continuous accumulation of genomic information, the evolutionary status of magnoliids has become a hot spot in plant phylogenetic research. Although great efforts have been made to study the evolution of magnoliids using molecular data from several representative species such as nuclear genome, plastid genome, mitochondrial genome, and transcriptome, the results of current studies on the phylogenetic status of magnoliids are inconsistent. Here, we systematically describe the current understanding of the molecular research on magnoliid phylogeny and review the differences in the evolutionary state of magnoliids. Understanding the research approaches and limitations of magnoliid phylogeny can guide research strategies to further improve the study of the phylogenetic evolution of magnoliids.

The chloroplast genome of black pepper (Piper nigrum L.) and its comparative analysis with related Piper species

Piper nigrum, also known as black pepper, is an economically and ecologically important crop of the genus Piper. It has been titled as the king of spices due to its wide consumption throughout the world. In the present investigation, the chloroplast genome of P. nigrum has been assembled from a whole genome sequence by integrating the short and long reads generated through Illumina and PacBio platforms, respectively. The chloroplast genome was observed to be 161,522 bp in size, having a quadripartite structure with a large single copy (LSC) region of 89,153 bp and a small single copy (SSC) region of 18,255 bp separated by a copy of inverted repeats (IRs), each 27,057 bp in length. Taking into consideration all the duplicated genes, a total of 131 genes were observed, which included 81 protein-coding genes, 37 tRNAs, 4 rRNAs, and 1 pseudogene. Individually, the LSC region consisted of 83 genes, the SSC region had 13 genes, and 18 genes were present in each IR region. Additionally, 216 SSRs were detected and 11 of these were validated through amplification in 12 species of Piper. The features of the chloroplast genome have been compared with those of the genus Piper. Our results provide useful insights into evolutionary and molecular studies of black pepper which will contribute to its further genetic improvement and breeding.

The complete chloroplast genome of Elsholtzia fruticosa (D. Don) Rehd. (Labiatae), an ornamental plant with high medicinal value

Elsholtzia fruticosa is an ornamental plant with high medicinal value. In this study, we sequenced and analyzed the complete chloroplast (cp) genome of the species. The complete cp sequence is 151,550 bp, including the large single-copy (LSC) region of 82,778 bp, the small single-copy (SSC) region of 17,492 bp, and a pair of invert repeats (IRs) regions of 25,640 bp. It encodes 132 unique genes in total, including 87 protein-coding genes, 37 transfer RNA genes (tRNAs), and eight ribosomal RNA genes (rRNAs). The comparative analysis of complete cp genomes showed that the genomic structure and gene order of E. fruticosa cps were conserved. The sequences of rps15, rps19, ycf1, ycf3, ycf15, psbL, psaI, trnG-UCC, trnS-GCU, trnR-UCU, trnL-UAG, trnP-UG, and trnL-UAA serve as hotspots for developing the DNA barcoding of Elsholtzia species. There are 49 SSR loci in the cp genome of E. fruticosa, among which the repeat numbers of mononucleotide, dinucleotide, trinucleotide, tetranucleotide, and pentanucleotides SSR are 37, 9, 3, 0, and 0, respectively. A total of 50 repeats were detected, including 15 forward repeats, seven reverse repeats, 26 palindromic repeats, and two complementary repeats. Phylogenetic analysis based on the complete cp genome and protein-coding DNA sequences of 26 plants indicates that E. fruticosa has a dose relationship with E. splendens and E. byeonsanensis.

Comparison of the chloroplast genomes and phylogenomic analysis of Elaeocarpaceae

Background

Elaeocarpaceae is a vital family in tropical and subtropical forests. Compared with the important position of Elaeocarpaceae species in forest ecosystem and the concern of medicinal value, the most research on Elaeocarpaceae are classification and taxonomy. Molecular systematics has corrected the morphological misjudgment, and it belongs to Oxalidales. Phylogenetic and divergence time estimates of Elaeocarpaceae is mostly constructed by using chloroplast gene fragments. At present, although there are reports on the chloroplast structure of Elaeocarpaceae, a comprehensive analysis of the chloroplast structure of Elaeocarpaceae is lacking.

Methods

To understand the variation in chloroplast sequence size and structure in Elaeocarpaceae, the chloroplast genomes of nine species were sequenced using the Illumina HiSeq 2500 platform and further assembled and annotated with Elaeocarpus japonicus and Sloanea sinensis (family Elaeocarpaceae) as references. A phylogenomic tree was constructed based on the complete chloroplast genomes of the 11 species representing five genera of Elaeocarpaceae. Chloroplast genome characteristics were examined by using Circoletto and IRscope software.

Results

The results revealed the following: (a) The 11 sequenced chloroplast genomes ranged in size from 157,546 to 159,400 bp. (b) The chloroplast genomes of Elaeocarpus, Sloanea, Crinodendron and Vallea lacked the rpl32 gene in the small single-copy (SSC) region. The large single-copy (LSC) region of the chloroplast genomes lacked the ndhK gene in Elaeocarpus, Vallea stipularis, and Aristotelia fruticosa. The LSC region of the chloroplast genomes lacked the infA gene in genus Elaeocarpus and Crinodendron patagua. (c) Through inverted repeat (IR) expansion and contraction analysis, a significant difference was found between the LSC/IRB and IRA/LSC boundaries among these species. Rps3 was detected in the neighboring regions of the LSC and IRb regions in Elaeocarpus. (d) Phylogenomic analysis revealed that the genus Elaeocarpus is closely related to Crinodendron patagua on an independent branch and Aristotelia fruticosa is closely related to Vallea stipularis, forming a clade with the genus Sloanea. Structural comparisons showed that Elaeocarpaceae diverged at 60 Mya, the genus Elaeocarpus diverged 53 Mya and that the genus Sloanea diverged 0.44 Mya. These results provide new insight into the evolution of the Elaeocarpaceae.

The complete plastid genome sequence of Chloranthus fortunei (A. Gray) Solms-Laub. in Chloranthaceae

Chloranthus fortunei (A. Gray) Solms-Laub. is a perennial herb in a basal angiosperm family Chloranthaceae. Here, we reported the complete plastid genome of C. fortunei using Illumina short-read data. The total genome size was 157,063 bp in length, containing 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. The gene content and order were consistent with previously reported Chloranthus plastid genomes. The overall GC content of the C. fortunei plastid genome was 39.0%. In the phylogenetic result, genus Chloranthus was monophyletic and divided into two subclades: C. japonicus+C. angustifolius+C. fortunei, and C. henryi+C. spicatus+C. erectus. Our phylogenetic result was consistent with previous phylogenetic studies, and was supported by a previously proposed infrageneric classification of the genus Chloranthus.

Comparative analysis of the complete chloroplast genomes of six threatened subgenus Gynopodium (Magnolia) species

BACKGROUND

The subgenus Gynopodium belonging to genus Magnolia have high ornamental, economic, and ecological value. Subgenus Gynopodium contains eight species, but six of these species are threatened. No studies to date have characterized the characteristics of the chloroplast genomes (CPGs) within subgenus Gynopodium species. In this study, we compared the structure of CPGs, identified the mutational hotspots and resolved the phylogenetic relationship of subgenus Gynopodium.

RESULTS

The CPGs of six subgenus Gynopodium species ranged in size from 160,027 bp to 160,114 bp. A total of 131 genes were identified, including 86 protein-coding genes, eight ribosomal RNA genes, and 37 transfer RNA genes. We detected neither major expansions or contractions in the inverted repeat region, nor rearrangements or insertions in the CPGs of six subgenus Gynopodium species. A total of 300 large repeat sequences (forward, reverse, and palindrome repeats), 847 simple sequence repeats, and five highly variable regions were identified. One gene (ycf1) and four intergenic regions (psbA-trnH-GUG, petA-psbJ, rpl32-trnL-UAG, and ccsA-ndhD) were identified as mutational hotspots by their high nucleotide diversity (Pi) values (≥ 0.004), which were useful for species discrimination. Maximum likelihood and Bayesian inference trees were concordant and indicated that Magnoliaceae consisted of two genera Liriodendron and Magnolia. Six species of subgenus Gynopodium clustered as a monophyletic clade, forming a sister clade with subgenus Yulania (BS = 100%, PP = 1.00). Due to the non-monophyly of subgenus Magnolia, subgenus Gynopodium should be treated as a section of Magnolia. Within section Gynopodium, M. sinica diverged first (posterior probability = 1, bootstrap = 100), followed by M. nitida, M. kachirachirai and M. lotungensis. M. omeiensis was sister to M. yunnanensis (posterior probability = 0.97, bootstrap = 50).

CONCLUSION

The CPGs and characteristics information provided by our study could be useful in species identification, conservation genetics and resolving phylogenetic relationships of Magnoliaceae species.

Comparative analysis of the chloroplast genomes of eight Piper species and insights into the utilization of structural variation in phylogenetic analysis

With more than 2000 species, Piper is regarded as having high medicinal, cosmetic, and edible value. There also remain some taxonomic and evolutionary uncertainties about the genus. This study performed chloroplast genome sequencing of eight poorly studied Piper species and a comparative analysis with black pepper (Piper nigrum). All examined species were highly similar in gene content, with 79 protein-coding genes, 24 tRNAs, and four rRNAs. They also harbored significant structural differences: The number of SSRs ranged from 63 to 87, over 10,000 SNPs were detected, and over 1,000 indels were found. The spatial distribution of structural differences was uneven, with the IR and LSC being relatively more conserved and the SSC region highly variable. Such structural variations of the chloroplast genome can help in evaluating the phylogenetic relationships between species, deciding some hard-to-distinguish evolutionary relationships, or eliminating improper markers. The SSC region may be evolving at high speed, and some species showed a high degree of sequence variation in the SSC region, which seriously affected marker sequence detection. Conversely, CDS sequences tended to lack variation, and some CDSs can serve as ideal markers for phylogenetic reconstruction. All told, this study provides an effective strategy for selecting chloroplast markers, analyzing difficult-to-distinguish phylogenetic relationships and avoiding the taxonomic errors caused by high degree of sequence variations.

Chloroplast Genomes of Genus Tilia: Comparative Genomics and Molecular Evolution

Tilia is a complex genus in the family Malvaceae that has high ecological and economical values. Owing to the lack of sufficient distinguishable morphological and molecular characteristics, interspecific relationships in this genus are not clear. Chloroplast (cp) genomes are small, meanwhile most angiosperms usually undergo matrilineal inheritance. Consequently, they can be used in molecular systematics and phylogenetic analyses. Here, we sequenced and assembled cp genomes from T. endochrysea, T. tomentosa, T. miqueliana, T. americana and T. cordata, and compared them with those of seven previously reported Tilia species. Similar gene contents, gene orders and GC contents existed among the 12 cp genomes, which ranged from 162,564 to 162,855 bp and encoded 113 unique genes. Abundant simple sequence repeats (119–127) and dispersed repeats (97–135) were detected in Tilia cp genomes. In total, 11 hypervariable regions were identified that could be suitable for species identification and phylogenetic studies. A phylogenetic analysis of Malvaceae based on 5 hypervariable genes (matK + ndhF + rpoB + rpoC2+ycf1) revealed that all eight subfamilies were monophyletic groups. Additionally, the genus Tilia was divided into three groups on the basis of all 521 molecular variation loci. The current study provides valuable insights into the genomic evolution of the genus Tilia.

Structural plastome evolution in holoparasitic Hydnoraceae with special focus on inverted and direct repeats

Plastome condensation during adaptation to a heterotrophic lifestyle is generally well understood and lineage-independent models have been derived. However, understanding the evolutionary trajectories of comparatively old heterotrophic lineages, that are on the cusp of a minimal plastomes, is essential to complement and expand current knowledge. We study Hydnoraceae, one of the oldest and least investigated parasitic angiosperm lineages. Plastome comparative genomics, using seven out of eight known species of the genus Hydnora and three species of Prosopanche, reveal a high degree of structural similarity and shared gene content; contrasted by striking dissimilarities with respect to repeat content (inverted and direct repeats). We identified varying IR content and positions, likely resulting from multiple, independent evolutionary events and a direct repeat gain in Prosopanche. Considering different evolutionary trajectories and based on a fully resolved and supported species-level phylogenetic hypothesis, we describe three possible, distinct models to explain the Hydnoraceae plastome states. For comparative purposes we also report the first plastid genomes for the closely related autotrophic genera Lactoris (Lactoridaceae) and Thottea (Aristolochiaceae).

Fifteen complete chloroplast genomes of Trapa species (Trapaceae): insight into genome structure, comparative analysis and phylogenetic relationships

BACKGROUND

Trapa L. is a floating-leaved aquatic plant with important economic and ecological values. However, the species identification and phylogenetic relationship within Trapa are still controversial, which necessitates the need for plastid genome information of Trapa. In this study, complete chloroplast genomes of 13 Trapa species/taxa were sequenced and annotated. Combined with released sequences, comparative analyses of chloroplast genomes were performed on the 15 Trapa species/taxa for the first time.

RESULTS

The Trapa chloroplast genomes exhibited typical quadripartite structures with lengths from 155,453 to 155,559 bp. The gene orders and contents within Trapa were conservative, but several changes were found in the microstructure. The intron loss of rpl2, also detected in Lythraceae, was found in all Trapa species/taxa, suggesting close genetic relationship between Lythraceae and Trapaceae. Notably, two small-seed species (T. incisa and T. maximowiczii) showed the smallest genome size with 155,453 and 155,477 bp, respectively. Each cp genome contained the same 130 genes consisting of 85 protein-coding genes, 37 tRNA genes and 8 rRNA genes. Trapa species/taxa showed 37 (T. incisa and T. maximowiczii) to 41 (T. sibirica) long repeats, including forward, palindromic, reversed and complementary repeats. There were 110 (T. quadrispinosa) to 123 (T. incisa and T. maximowiczii) SSR (simple sequence repeat) loci in Trapa chloroplast genomes. Comparative analyses revealed that two hotspot regions (atpA-atpF and rps2-rpoC2) in Trapa chloroplast genomes could be served as potential molecular markers. Three phylogenetic analyses (ML, MP and BI) consistently showed that there were two clusters within Trapa, including large- and small-seed species/taxa, respectively; for the large-seed Trapa, they clustered according to their geographical origin and tubercle morphology on the surface of seeds.

CONCLUSION

In summary, we have acquired the sequences of 13 Trapa chloroplast genomes, and performed the comparative analyses within Trapa for the first time. The results have helped us better identify the Trapa species/taxa and deepen the understanding of genetic basis and phylogenetic relationship of Trapa, which will facilitate the effective management and utilization of the important genetic resources in the future.

Plastid Phylogenomic Analysis of Tordylieae Tribe (Apiaceae, Apioideae)

Based on the nrDNA ITS sequence data, the Tordylieae tribe is recognized as monophyletic with three major lineages: the subtribe Tordyliinae, the Cymbocarpum clade, and the Lefebvrea clade. Recent phylogenomic investigations showed incongruence between the nuclear and plastid genome evolution in the tribe. To assess phylogenetic relations and structure evolution of plastomes in Tordylieae, we generated eleven complete plastome sequences using the genome skimming approach and compared them with the available data from this tribe and close relatives. Newly assembled plastomes had lengths ranging from 141,148 to 150,103 base pairs and contained 122-127 genes, including 79-82 protein-coding genes, 35-37 tRNAs, and 8 rRNAs. We observed substantial differences in the inverted repeat length and gene content, accompanied by a complex picture of multiple JLA and JLB shifts. In concatenated phylogenetic analyses, Tordylieae plastomes formed at least three not closely related lineages with plastomes of the Lefebvrea clade as a sister group to plastomes from the Selineae tribe. The newly obtained data have increased our knowledge on the range of plastome variability in Apiaceae.

Phylogenetic relationships in the Sorghum genus based on sequencing of the chloroplast and nuclear genes

Sorghum [Sorghum bicolor (L.) Moench] is an important food crop with a diverse gene pool residing in its wild relatives. A total of 15 sorghum accessions from the unexploited wild gene pool of the Sorghum genus, representing the five subgenera, were sequenced, and the complete chloroplast genomes and 99 common single-copy concatenated nuclear genes were assembled. Annotation of the chloroplast genomes identified a total of 81 protein-coding genes, 38 tRNA, and four rRNA genes. The gene content and gene order among the species was identical. A total of 153 nonsynonymous amino acid changes in 40 genes were identified across the species. Phylogenetic analysis of both the whole chloroplast genome and nuclear genes revealed a similar topology with two distinct clades within the genus. The species within the subgenera Eusorghum, Chaetosorghum, and Heterosorghum clustered in one clade, whereas the species within the subgenera Parasorghum and Stiposorghum clustered in a second clade. However, the subgenera Parasorghum and Stiposorghum were not monophyletic, suggesting the need for further research to resolve the relationships within this group. The close relationship between the two monotypic subgenera Chaetosorghum and Heterosorghum suggests that species within these subgenera could be considered as one group. This analysis provides an improved understanding of the genetic relationships within the Sorghum genus and defines diversity in wild sorghum species that may be useful for crop improvement.

Comparative Analysis of Genomic and Transcriptome Sequences Reveals Divergent Patterns of Codon Bias in Wheat and Its Ancestor Species

The synonymous codons usage shows a characteristic pattern of preference in each organism. This codon usage bias is thought to have evolved for efficient protein synthesis. Synonymous codon usage was studied in genes of the hexaploid wheat Triticum aestivum (AABBDD) and its progenitor species, Triticum urartu (AA), Aegilops tauschii (DD), and Triticum turgidum (AABB). Triticum aestivum exhibited stronger usage bias for G/C-ending codons than did the three progenitor species, and this bias was especially higher compared to T. turgidum and Ae. tauschii. High GC content is a primary factor influencing codon usage in T. aestivum. Neutrality analysis showed a significant positive correlation (p<0.001) between GC12 and GC3 in the four species with regression line slopes near zero (0.16–0.20), suggesting that the effect of mutation on codon usage was only 16–20%. The GC3s values of genes were associated with gene length and distribution density within chromosomes. tRNA abundance data indicated that codon preference corresponded to the relative abundance of isoaccepting tRNAs in the four species. Both mutation and selection have affected synonymous codon usage in hexaploid wheat and its progenitor species. GO enrichment showed that GC biased genes were commonly enriched in physiological processes such as photosynthesis and response to acid chemical. In some certain gene families with important functions, the codon usage of small parts of genes has changed during the evolution process of T. aestivum.

Complete chloroplast genome sequences of the medicinal plant Piper hancei

Chloroplast genome sequences have been used in phylogenetic and population genetics studies. Here, we assembled the chloroplast genome of Piper hancei Maxim. that is a traditional Chinese medicine. The genome length was 161,476 bp and included a pair of inverted repeats of 27,058 bp, a large single-copy region of 89,144 bp and a small single-copy region of 18,216 bp. It contained 113 different genes, including 79 protein-coding genes, 30 transfer RNA (tRNA), and four ribosomal RNA genes. Moreover, we also identified 82 SSRs. The phylogenetic inference based on the whole chloroplast genome of 20 taxa showed P. hancei was sister to P. kadsura.

Unprecedented Intraindividual Structural Heteroplasmy in Eleocharis (Cyperaceae, Poales) Plastomes

Plastid genomes (plastomes) of land plants have a conserved quadripartite structure in a gene-dense unit genome consisting of a large inverted repeat that separates two single copy regions. Recently, alternative plastome structures were suggested in Geraniaceae and in some conifers and Medicago the coexistence of inversion isomers has been noted. In this study, plastome sequences of two Cyperaceae, Eleocharis dulcis (water chestnut) and Eleocharis cellulosa (gulf coast spikerush), were completed. Unlike the conserved plastomes in basal groups of Poales, these Eleocharis plastomes have remarkably divergent features, including large plastome sizes, high rates of sequence rearrangements, low GC content and gene density, gene duplications and losses, and increased repetitive DNA sequences. A novel finding among these features was the unprecedented level of heteroplasmy with the presence of multiple plastome structural types within a single individual. Illumina paired-end assemblies combined with PacBio single-molecule real-time sequencing, long-range polymerase chain reaction, and Sanger sequencing data identified at least four different plastome structural types in both Eleocharis species. PacBio long read data suggested that one of the four E. dulcis plastome types predominates.

Phylogenetics and comparative plastome genomics of two of the largest genera of angiosperms, Piper and Peperomia (Piperaceae)

Biological radiations provide unique opportunities to understand the evolution of biodiversity. One such radiation is the pepper plant family Piperaceae, an early-diverging and mega-diverse lineage that could serve as a model to study the diversification of angiosperms. However, traditional genetic markers lack sufficient variation for such studies, and testing hypotheses on poorly resolved phylogenetic frameworks becomes challenging. Limited genomic data is available for Piperaceae, which contains two of the largest genera of angiosperms, Piper (>2100 species) and Peperomia (>1300 species). To address this gap, we used genome skimming to assemble and annotate whole plastomes (152-161kbp) and >5kbp nuclear ribosomal DNA region from representatives of Piper and Peperomia. We conducted phylogenetic and comparative genomic analyses to study plastome evolution and investigate the role of hybridization in this group. Plastome phylogenetic trees were well resolved and highly supported, with a hard incongruence observed between plastome and nuclear phylogenetic trees suggesting hybridization in Piper. While all plastomes of Piper and Peperomia had the same gene content and order, there were informative structural differences between them. First, ycf1 was more variable and longer in Piper than Peperomia, extending well into the small single copy region by thousands of base pairs. We also discovered previously unknown structural variation in 14 out of 25 Piper taxa, tandem duplication of the trnH-GUG gene resulting in an expanded large single copy region. Other early-diverging angiosperms have a duplicated trnH-GUG, but the specific rearrangement we found is unique to Piper and serves to refine knowledge of relationships among early-diverging angiosperms. Our study demonstrates that genome skimming is an efficient approach to produce plastome assemblies for comparative genomics and robust phylogenies of species-rich plant genera.

The complete plastome sequences of Pseudowintera colorata and Tasmannia lanceolata (Winteraceae - Canellales)

In this study, we report the complete plastome sequences of two Winteraceae taxa, Pseudowintera colorata (MT555077) and Tasmannia lanceolata (MT555078). Both plastomes show typical quadripartite structure. The plastome size of P. colorata is 161,675 bp, which consists of 89,583 bp large single-copy (LSC), 18,606 bp small single-copy (SSC), and 26,743 bp inverted repeat (IR) regions. The plastome size of T. lanceolata is 160,424 bp, which consist of 88,589 bp LSC, 18,351 bp SSC, and 26,742 bp IR regions. Both plastomes contain 113 genes, including 79 protein-coding, 30 tRNA, and four rRNA genes. Sixteen genes contain one intron and two genes (clpP and ycf3) have two introns. Ninety-three and 89 simple sequence repeat (SSR) loci are scattered in the P. colorata and T. lanceolata plastomes, respectively. Our phylogenetic tree shows the relationship of (T. lanceolate (P. colorata,Drimys granadensis)) in the Winteraceae. The Canellales (incl. Winteraceae) are the sister group of Piperales.

Comparative and Phylogenetic Analysis of Complete Chloroplast Genomes in Eragrostideae (Chloridoideae, Poaceae)

Eragrostideae Stapf, the second-largest tribe in Chloridoideae (Poaceae), is a taxonomically complex tribe. In this study, chloroplast genomes of 13 Eragrostideae species were newly sequenced and used to resolve the phylogenetic relationships within Eragrostideae. Including seven reported chloroplast genomes from Eragrostideae, the genome structure, number and type of genes, codon usage, and repeat sequences of 20 Eragrostideae species were analyzed. The length of these chloroplast genomes varied from 130,773 bp to 135,322 bp. These chloroplast genomes showed a typical quadripartite structure, including a large single-copy region (77,993–80,643 bp), a small single-copy region (12,410–12,668 bp), and a pair of inverted repeats region (19,394–21,074 bp). There were, in total, 129–133 genes annotated in the genome, including 83–87 protein-coding genes, eight rRNA genes, and 38 tRNA genes. Forward and palindromic repeats were the most common repeat types. In total, 10 hypervariable regions (rpl22, rpoA, ndhF, matK, trnG–UCC-trnT–GGU, ndhF–rpl32, ycf4–cemA, rpl32–trnL–UAG, trnG–GCC–trnfM–CAU, and ccsA–ndhD) were found, which can be used as candidate molecular markers for Eragrostideae. Phylogenomic studies concluded that Enneapogon diverged first, and Eragrostis including Harpachne is the sister to Uniola. Furthermore, Harpachne harpachnoides is considered as a species of Eragrostis based on morphological and molecular evidence. In addition, the interspecies relationships within Eragrostis are resolved based on complete chloroplast genomes. This study provides useful chloroplast genomic information for further phylogenetic analysis of Eragrostideae.

Plastid Genomes of Flowering Plants: Essential Principles

The plastid genome (plastome ) has proved a valuable source of data for evaluating evolutionary relationships among angiosperms. Through basic and applied approaches, plastid transformation technology offers the potential to understand and improve plant productivity, providing food, fiber, energy, and medicines to meet the needs of a burgeoning global population. The growing genomic resources available to both phylogenetic and biotechnological investigations is allowing novel insights and expanding the scope of plastome research to encompass new species. In this chapter, we present an overview of some of the seminal and contemporary research that has contributed to our current understanding of plastome evolution and attempt to highlight the relationship between evolutionary mechanisms and the tools of plastid genetic engineering.

Variation in Plastome Sizes Accompanied by Evolutionary History in Monogenomic Triticeae (Poaceae: Triticeae)

To investigate the pattern of chloroplast genome variation in Triticeae, we comprehensively analyzed the indels in protein-coding genes and intergenic sequence, gene loss/pseudonization, intron variation, expansion/contraction in inverted repeat regions, and the relationship between sequence characteristics and chloroplast genome size in 34 monogenomic Triticeae plants. Ancestral genome reconstruction suggests that major length variations occurred in four-stem branches of monogenomic Triticeae followed by independent changes in each genus. It was shown that the chloroplast genome sizes of monogenomic Triticeae were highly variable. The chloroplast genome of Pseudoroegneria, Dasypyrum, Lophopyrum, Thinopyrum, Eremopyrum, Agropyron, Australopyrum, and Henradia in Triticeae had evolved toward size reduction largely because of pseudogenes elimination events and length deletion fragments in intergenic. The Aegilops/Triticum complex, Taeniatherum, Secale, Crithopsis, Herteranthelium, and Hordeum in Triticeae had a larger chloroplast genome size. The large size variation in major lineages and their subclades are most likely consequences of adaptive processes since these variations were significantly correlated with divergence time and historical climatic changes. We also found that several intergenic regions, such as petN-trnC and psbE-petL containing unique genetic information, which can be used as important tools to identify the maternal relationship among Triticeae species. Our results contribute to the novel knowledge of plastid genome evolution in Triticeae.

Chloroplast genome features of an important medicinal and edible plant: Houttuynia cordata (Saururaceae)

Houttuynia cordata (Saururaceae), an ancient and relic species, has been used as an important medicinal and edible plant in most parts of Asia. However, because of the lack of genome information and reliable molecular markers, studies on its population structure, or phylogenetic relationships with other related species are still rare. Here, we de novo assembled the complete chloroplast (cp) genome of H. cordata using the integration of the long PacBio and short Illumina reads. The cp genome of H. cordata showed a typical quadripartite cycle of 160,226 bp. This included a pair of inverted repeats (IRa and IRb) of 26,853 bp, separated by a large single-copy (LSC) region of 88,180 bp and a small single-copy (SSC) region of 18,340 bp. A total of 112 unique genes, including 79 protein-coding genes, 29 tRNA genes, and four rRNA genes, were identified in this cp genome. Eighty-one genes were located on the LSC region, 13 genes were located on the SSC region, and 17 two-copy genes were located on the IR region. Additionally, 48 repeat sequences and 86 SSR loci, which can be used as genomic markers for population structure analysis, were also detected. Phylogenetic analysis using 21 cp genomes of the Piperales family demonstrated that H. cordata had a close relationship with the species within the Aristolochia genus. Moreover, the results of mVISTA analysis and comparisons of IR regions demonstrated that the cp genome of H. cordata was conserved with that of the Aristolochia species. Our results provide valuable information for analyzing the genetic diversity and population structure of H. cordata, which can contribute to further its genetic improvement and breeding.

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.

Primers for complete chloroplast genome sequencing in Magnolia

PREMISE

A new set of primers was developed for sequencing of whole chloroplast genomes of Magnolia species and gap-filling of unfinished genomes.

METHODS AND RESULTS

Two hundred and fifty primers were newly designed based on two previously reported chloroplast genomes from two different genera in Magnoliaceae. A total of 134 primer pairs, including the ones developed in this study and 18 previously reported ones, were enough to cover the entire chloroplast genome sequences in Magnoliaceae. Four species from different sections of Magnolia (M. dealbata, M. fraseri var. pyramidata, M. liliiflora, and M. odora) were used to show the general application of these primers to chloroplast genome sequencing in Magnolia.

CONCLUSIONS

Using the developed primers, four Magnolia chloroplast genomes were successfully assembled. These results show the utility of these primers across Magnolia and their potential use for phylogenetic studies, DNA barcoding, and population genetics in this group.

The complete chloroplast genome of tulip tree, Liriodendron tulifipera L. (Magnoliaceae): investigation of intra-species chloroplast variations

Liriodendron tulifipera L. belongs to Magnoliaceae which is one of the basal angiosperm families. To understand intra-species variations on chloroplast genome in Liriodendron genus, we presented complete chloroplast genome of L. tulifipera, which is 156,387 bp long and has four subregions: 85,606 bp of large single copy (LSC) and 18,778 bp of small single copy (SSC) regions are separated by 26,002 bp of inverted repeat (IR) regions including 129 genes (84 coding genes, 8 rRNAs, and 37 tRNAs). The overall GC content of the chloroplast genome is 37.0% and those in the LSC, SSC, and IR regions are 34.9%, 30.5%, and 42.8%, respectively. Twelve single nucleotide polymorphisms (SNPs) located in one region and one insertion and deletion are found between two L. tulifipera genomes. INDEL Phylogenetic trees show that two L. tulifipera chloroplasts are clustered together and are sister to Magnolia species.

Complete Chloroplast Genomes and Comparative Analysis of Sequences Evolution among Seven Aristolochia (Aristolochiaceae) Medicinal Species

Aristolochiaceae, comprising about 600 species, is a unique plant family containing aristolochic acids (AAs). In this study, we sequenced seven species of Aristolochia, and retrieved eleven chloroplast (cp) genomes published for comparative genomics analysis and phylogenetic constructions. The results show that the cp genomes had a typical quadripartite structure with conserved genome arrangement and moderate divergence. The cp genomes range from 159,308 bp to 160,520 bp in length and have a similar GC content of 38.5%–38.9%. A total number of 113 genes were identified, including 79 protein-coding genes, 30 tRNAs and four rRNAs. Although genomic structure and size were highly conserved, the IR-SC boundary regions were variable between these seven cp genomes. The trnH-GUG genes, are one of major differences between the plastomes of the two subgenera Siphisia and Aristolochia. We analyzed the features of nucleotide substitutions, distribution of repeat sequences and simple sequences repeats (SSRs), positive selections in the cp genomes, and identified 16 hotspot regions for genomes divergence that could be utilized as potential markers for phylogeny reconstruction. Phylogenetic relationships of the family Aristolochiaceae inferred from the 18 cp genome sequences were consistent and robust, using maximum parsimony (MP), maximum likelihood (ML), and Bayesian analysis (BI) methods.

Novel genetic code and record-setting AT-richness in the highly reduced plastid genome of the holoparasitic plant Balanophora

Plastid genomes (plastomes) vary enormously in size and gene content among the many lineages of nonphotosynthetic plants, but key lineages remain unexplored. We therefore investigated plastome sequence and expression in the holoparasitic and morphologically bizarre Balanophoraceae. The two Balanophora plastomes examined are remarkable, exhibiting features rarely if ever seen before in plastomes or in any other genomes. At 15.5 kb in size and with only 19 genes, they are among the most reduced plastomes known. They have no tRNA genes for protein synthesis, a trait found in only three other plastid lineages, and thus Balanophora plastids must import all tRNAs needed for translation. Balanophora plastomes are exceptionally compact, with numerous overlapping genes, highly reduced spacers, loss of all cis-spliced introns, and shrunken protein genes. With A+T contents of 87.8% and 88.4%, the Balanophora genomes are the most AT-rich genomes known save for a single mitochondrial genome that is merely bloated with AT-rich spacer DNA. Most plastid protein genes in Balanophora consist of ≥90% AT, with several between 95% and 98% AT, resulting in the most biased codon usage in any genome described to date. A potential consequence of its radical compositional evolution is the novel genetic code used by Balanophora plastids, in which TAG has been reassigned from stop to tryptophan. Despite its many exceptional properties, the Balanophora plastome must be functional because all examined genes are transcribed, its only intron is correctly trans-spliced, and its protein genes, although highly divergent, are evolving under various degrees of selective constraint.

Green Biotechnology: A Brief Update on Plastid Genome Engineering

Plant genetic engineering has become an inevitable tool in the molecular breeding of crops. Significant progress has been made in the generation of novel plastid transformation vectors and optimized transformation protocols. There are several advantages of plastid genome engineering over conventional nuclear transformation. Some of the advantages include multigene engineering by expression of biosynthetic pathway genes as operons, extremely high-level expression of protein accumulation, lack of transgene silencing, etc. Transgene containment owing to maternal inheritance is another important advantage of plastid genome engineering. Chloroplast genome modification usually results in alteration of several thousand plastid genome copies in a cell. Several therapeutic proteins, edible vaccines, antimicrobial peptides, and industrially important enzymes have been successfully expressed in chloroplasts so far. Here, we critically recapitulate the latest developments in plastid genome engineering. Latest advancements in plastid genome sequencing are briefed. In addition, advancement of extending the toolbox for plastid engineering for selected applications in the area of molecular farming and production of industrially important enzyme is briefed.

Evolutionary Comparisons of the Chloroplast Genome in Lauraceae and Insights into Loss Events in the Magnoliids

Available plastomes of the Lauraceae show similar structure and varied size, but there has been no systematic comparison across the family. In order to understand the variation in plastome size and structure in the Lauraceae and related families of magnoliids, we here compare 47 plastomes, 15 newly sequenced, from 27 representative genera. We reveal that the two shortest plastomes are in the parasitic Lauraceae genus Cassytha, with lengths of 114,623 (C. filiformis) and 114,963 bp (C. capillaris), and that they have lost NADH dehydrogenase (ndh) genes in the large single-copy region and one entire copy of the inverted repeat (IR) region. The plastomes of the core Lauraceae group, with lengths from 150,749 bp (Nectandra angustifolia) to 152,739 bp (Actinodaphne trichocarpa), have lost trnI-CAU, rpl23, rpl2, a fragment of ycf2, and their intergenic regions in IRb region, whereas the plastomes of the basal Lauraceae group, with lengths from 157,577 bp (Eusideroxylon zwageri) to 158,530 bp (Beilschmiedia tungfangensis), have lost rpl2 in IRa region. The plastomes of Calycanthus (Calycanthaceae, Laurales) have lost rpl2 in IRb region, but the plastome of Caryodaphnopsis henryi (Lauraceae) remain intact, as do those of the nonLaurales magnoliid genera Piper, Liriodendron, and Magnolia. On the basis of our phylogenetic analysis and structural comparisons, different loss events occurred in different lineages of the Laurales, and fragment loss events in the IR regions have largely driven the contraction of the plastome in the Lauraceae. These results provide new insights into the evolution of the Lauraceae as well as the magnoliids as a whole.

Comparative chloroplast genomics between the invasive weed Mikania micrantha and its indigenous congener Mikania cordata: Structure variation, identification of highly divergent regions, divergence time estimation, and phylogenetic analysis

Mikania micrantha and Mikania cordata are the only two species in genus Mikania (Asteraceae) in China. They share very similar morphological and life-history characteristics but occupy quite different habitats. Most importantly, they generate totally different ecological consequences. While M. micrantha has become an exotic invasive weed, M. cordata exists as an indigenous species with no harmful effects on native plants or habitats. As a continuous study of our previously reported M. micrantha chloroplast (cp) genome, in this study we have further sequenced the M. cordata cp genome to (1) conduct a comparative genome analysis to gain insights into the mechanism of invasiveness; (2) develop cp markers to examine the population genetic adaptation of M. micrantha; and (3) screen variable genome regions of phylogenetic utility. The M. cordata chloroplast genome is 151,984 bp in length and displays a typical quadripartite structure. The number and distribution of protein coding genes, tRNA genes, and rRNA genes of M. cordata are identical to those of M. micrantha. The main difference lays in that the pseudogenization of ndhF and a 118-bp palindromic repeat only arises in M. cordata. Fourteen highly divergent regions, 235 base substitutions, and 58 indels were identified between the two cp genomes. Phylogenetic inferences revealed a sister relationship between M. micrantha and M. cordata whose divergence was estimated to occur around 1.78 million years ago (MYA). Twelve cpSSR loci were detected to be polymorphic and adopted to survey the genetic adaptation of M. micrantha populations. No cpSSR loci were found to undergo selection. Our results build a foundation to examine the invasive mechanism of Mikania weed.

Phylogenomics and barcoding of Panax: toward the identification of ginseng species

BACKGROUND

The economic value of ginseng in the global medicinal plant trade is estimated to be in excess of US$2.1 billion. At the same time, the evolutionary placement of ginseng (Panax ginseng) and the complex evolutionary history of the genus is poorly understood despite several molecular phylogenetic studies. In this study, we use a full plastome phylogenomic framework to resolve relationships in Panax and to identify molecular markers for species discrimination.

RESULTS

We used high-throughput sequencing of MBD2-Fc fractionated Panax DNA to supplement publicly available plastid genomes to create a phylogeny based on fully assembled and annotated plastid genomes from 60 accessions of 8 species. The plastome phylogeny based on a 163 kbp matrix resolves the sister relationship of Panax ginseng with P. quinquefolius. The closely related species P. vietnamensis is supported as sister of P. japonicus. The plastome matrix also shows that the markers trnC-rps16, trnS-trnG, and trnE-trnM could be used for unambiguous molecular identification of all the represented species in the genus.

CONCLUSIONS

MBD2 depletion reduces the cost of plastome sequencing, which makes it a cost-effective alternative to Sanger sequencing based DNA barcoding for molecular identification. The plastome phylogeny provides a robust framework that can be used to study the evolution of morphological characters and biosynthesis pathways of ginsengosides for phylogenetic bioprospecting. Molecular identification of ginseng species is essential for authenticating ginseng in international trade and it provides an incentive for manufacturers to create authentic products with verified ingredients.

The Complete Chloroplast Genome Sequence of Tree of Heaven (Ailanthus altissima (Mill.) (Sapindales: Simaroubaceae), an Important Pantropical Tree

Ailanthus altissima (Mill.) Swingle (Simaroubaceae) is a deciduous tree widely distributed throughout temperate regions in China, hence suitable for genetic diversity and evolutionary studies. Previous studies in A. altissima have mainly focused on its biological activities, genetic diversity and genetic structure. However, until now there is no published report regarding genome of this plant species or Simaroubaceae family. Therefore, in this paper, we first characterized A. altissima complete chloroplast genome sequence. The tree of heaven chloroplast genome was found to be a circular molecule 160,815 base pairs (bp) in size and possess a quadripartite structure. The A. altissima chloroplast genome contains 113 unique genes of which 79 and 30 are protein coding and transfer RNA (tRNA) genes respectively and also 4 ribosomal RNA genes (rRNA) with overall GC content of 37.6%. Microsatellite marker detection identified A/T mononucleotides as majority SSRs in all the seven analyzed genomes. Repeat analyses of seven Sapindales revealed a total of 49 repeats in A. altissima, Rhus chinensis, Dodonaea viscosa, Leitneria floridana, while Azadirachta indica, Boswellia sacra, and Citrus aurantiifolia had a total of 48 repeats. The phylogenetic analysis using protein coding genes revealed that A. altissima is a sister to Leitneria floridana and also suggested that Simaroubaceae is a sister to Rutaceae family. The genome information reported here could be further applied for evolution and invasion, population genetics, and molecular studies in this plant species and family.

The chloroplast genome sequence of Magnolia kobus DC. (Magnoliaceae)

As an endangered species, Magnolia kobus is distributed in Jeju island in Korea with only about 500–1000 individuals. In this study, we presented a complete chloroplast genome of M. kobus which is 159,443 bp and has four sub-regions: 87,484 bp of large single copy and 18,783 bp of small single copy regions are separated by 26,588 bp of inverted repeat regions including 113 genes (79 unique genes, four rRNAs and 30 tRNAs). Phylogenetic analysis using chloroplast genomes showed that M. kobus is a sister of M. insignis and M. laevifolia clade.