psbE-psbL and ndhA Intron, the Promising Plastid DNA Barcode of Fagopyrum

Chloroplast genome characterization of Uncaria guianensis and Uncaria tomentosa and evolutive dynamics of the Cinchonoideae subfamily

Uncaria species are used in traditional medicine and are considered of high therapeutic value and economic importance. This work describes the assembly and annotation of the chloroplast genomes of U. guianensis and U. tomentosa, as well as a comparative analysis. The genomes were sequenced on MiSeq Illumina, assembled with NovoPlasty, and annotated using CHLOROBOX GeSeq. Addictionaly, comparative analysis were performed with six species from NCBI databases and primers were designed in Primer3 for hypervariable regions based on the consensus sequence of 16 species of the Rubiaceae family and validated on an in-silico PCR in OpenPrimeR. The genome size of U. guianensis and U. tomentosa was 155,505 bp and 156,390 bp, respectively. Both Species have 131 genes and GC content of 37.50%. The regions rpl32-ccsA, ycf1, and ndhF-ccsA showed the three highest values of nucleotide diversity within the species of the Rubiaceae family and within the Uncaria genus, these regions were trnH-psbA, psbM-trnY, and rps16-psbK. Our results indicates that the primer of the region ndhA had an amplification success for all species tested and can be promising for usage in the Rubiaceae family. The phylogenetic analysis recovered a congruent topology to APG IV. The gene content and the chloroplast genome structure of the analyzed species are conserved and most of the genes are under negative selection. We provide the cpDNA of Neotropical Uncaria species, an important genomic resource for evolutionary studies of the group.

Comparative plastome genomics, taxonomic delimitation and evolutionary divergences of Tetraena hamiensis var. qatarensis and Tetraena simplex (Zygophyllaceae)

The Zygophyllum and Tetraena genera are intriguingly important ecologically and medicinally. Based on morphological characteristics, T. hamiensis var. qatarensis, and T. simplex were transferred from Zygophyllum to Tetraena with the least genomic datasets available. Hence, we sequenced the T. hamiensis and T. simplex and performed in-depth comparative genomics, phylogenetic analysis, and estimated time divergences. The complete plastomes ranged between 106,720 and 106,446 bp-typically smaller than angiosperms plastomes. The plastome circular genomes are divided into large single-copy regions (~ 80,964 bp), small single-copy regions (~ 17,416 bp), and two inverted repeats regions (~ 4170 bp) in both Tetraena species. An unusual shrinkage of IR regions 16-24 kb was identified. This resulted in the loss of 16 genes, including 11 ndh genes which encode the NADH dehydrogenase subunits, and a significant size reduction of Tetraena plastomes compared to other angiosperms. The inter-species variations and similarities were identified using genome-wide comparisons. Phylogenetic trees generated by analyzing the whole plastomes, protein-coding genes, matK, rbcL, and cssA genes exhibited identical topologies, indicating that both species are sisters to the genus Tetraena and may not belong to Zygophyllum. Similarly, based on the entire plastome and proteins coding genes datasets, the time divergence of Zygophyllum and Tetraena was 36.6 Ma and 34.4 Ma, respectively. Tetraena stem ages were 31.7 and 18.2 Ma based on full plastome and protein-coding genes. The current study presents the plastome as a distinguishing and identification feature among the closely related Tetraena and Zygophyllum species. It can be potentially used as a universal super-barcode for identifying plants.

Complete chloroplast genome structure of four Ulmus species and Hemiptelea davidii and comparative analysis within Ulmaceae species

In this study, the chloroplast (cp) genomes of Hemiptelea davidii, Ulmus parvifolia, Ulmus lamellosa, Ulmus castaneifolia, and Ulmus pumila ‘zhonghuajinye’ were spliced, assembled and annotated using the Illumina HiSeq PE150 sequencing platform, and then compared to the cp genomes of other Ulmus and Ulmaceae species. The results indicated that the cp genomes of the five sequenced species showed a typical tetrad structure with full lengths ranging from 159,113 to 160,388 bp. The large single copy (LSC), inverted repeat (IR), and small single copy (SSC) lengths were in the range of 87,736–88,466 bp, 26,317–26,622 bp and 18,485–19,024 bp, respectively. A total of 130–131 genes were annotated, including 85–86 protein-coding genes, 37 tRNA genes and eight rRNA genes. The GC contents of the five species were similar, ranging from 35.30 to 35.62%. Besides, the GC content was different in different region and the GC content in IR region was the highest. A total of 64-133 single sequence repeat (SSR) loci were identified among all 21 Ulmaceae species. The (A)_n and (T)_n types of mononucleotide were highest in number, and the lengths were primarily distributed in 10–12 bp, with a clear AT preference. A branch-site model and a Bayes Empirical Bayes analysis indicated that the rps15 and rbcL had the positive selection sites. Besides, the analysis of mVISTA and sliding windows got a lot of hotspots such as trnH/psbA, rps16/trnQ, trnS/trnG, trnG/trnR and rpl32/trnL, which could be utilized as potential markers for the species identification and phylogeny reconstruction within Ulmus in the further studies. Moreover, the evolutionary tree of Ulmaceae species based on common protein genes, whole cp genome sequences and common genes in IR region of the 23 Ulmaceae species were constructed using the ML method. The results showed that these Ulmaceae species were divided into two branches, one that included Ulmus, Zelkova and Hemiptelea, among which Hemiptelea was the first to differentiate and one that included Celtis, Trema, Pteroceltis, Gironniera and Aphananthe. Besides, these variations found in this study could be used for the classification, identification and phylogenetic study of Ulmus species. Our study provided important genetic information to support further investigations into the phylogenetic development and adaptive evolution of Ulmus and Ulmaceae species.

Plastome comparison and phylogenomics of Fagopyrum (Polygonaceae): insights into sequence differences between Fagopyrum and its related taxa

BACKGROUND

Fagopyrum (Polygonaceae) is a small plant lineage comprised of more than fifteen economically and medicinally important species. However, the phylogenetic relationships of the genus are not well explored, and the characteristics of Fagopyrum chloroplast genomes (plastomes) remain poorly understood so far. It restricts the comprehension of species diversity in Fagopyrum. Therefore, a comparative plastome analysis and comprehensive phylogenomic analyses are required to reveal the taxonomic relationship among species of Fagopyrum.

RESULTS

In the current study, 12 plastomes were sequenced and assembled from eight species and two varieties of Fagopyrum. In the comparative analysis and phylogenetic analysis, eight previously published plastomes of Fagopyrum were also included. A total of 49 plastomes of other genera in Polygonaceae were retrieved from GenBank and used for comparative analysis with Fagopyrum. The variation of the Fagopyrum plastomes is mainly reflected in the size and boundaries of inverted repeat/single copy (IR/SC) regions. Fagopyrum is a relatively basal taxon in the phylogenomic framework of Polygonaceae comprising a relatively smaller plastome size (158,768-159,985 bp) than another genus of Polygonaceae (158,851-170,232 bp). A few genera of Polygonaceae have nested distribution of the IR/SC boundary variations. Although most species of Fagopyrum show the same IRb/SC boundary with species of Polygonaceae, only a few species show different IRa/SC boundaries. The phylogenomic analyses of Fagopyrum supported the cymosum and urophyllum groups and resolved the systematic position of subclades within the urophyllum group. Moreover, the repeat sequence types and numbers were found different between groups of Fagopyrum. The plastome sequence identity showed significant differences between intra-group and inter-group.

CONCLUSIONS

The deletions of intergenic regions cause a short length of Fagopyrum plastomes, which may be the main reason for plastome size diversity in Polygonaceae species. The phylogenomic reconstruction combined with the characteristics comparison of plastomes supports grouping within Fagopyrum. The outcome of these genome resources may facilitate the taxonomy, germplasm resources identification as well as plant breeding of Fagopyrum.

DNA Barcoding and Phylogenomic Analysis of the Genus Fritillaria in China Based on Complete Chloroplast Genomes

The Fritillaria is an extremely complicated genus in taxonomy and phylogeny, which contains numerous medicinal species in China. Both traditional characteristic-based taxonomy and universal DNA barcodes (ITS, trnH-psbA, and rbcL) are difficult to effectively identify the species. Here, we generated a large dataset of chloroplast genomes from multiple accessions per species of Fritillaria to evaluate their effectiveness in species discrimination. Moreover, phylogeny of species in China was explored based on the complete chloroplast genomes, and then divergence times of each node were estimated. The results showed that all 21 species in Fritillaria here (including two suspicious species) could be correctly discriminated using cpDNA genomes except F. cirrhosa, which suggested that DNA super-barcode could greatly enhance species discriminatory resolution for complicated genera. Furthermore, four regions (ycf1, matK-trnG-GCC, rpoC1, and matK) gained remarkably higher resolution than that of other plastid regions, but only matK might be suitable to identify Fritillaria species in consideration of its lengths. Phylogenomic analysis showed that the subgenus Fritillaria in China was divided into four major clades with obvious geographic structure. Among them, Clade I, mainly distributed in southwest China, was a young and complicated group. Moreover, according to the analysis, taxonomic treatments of the two suspicious species, namely "F. omeiensis" and "F. hupehensis" in Flora of China (2000) are questionable and might need further revision. Molecular dating revealed that both origin and divergence of subgenus Fritillaria, as well as its four major clades, were significantly associated with geological and climatic fluctuations during the Middle to Late Miocene. This study would enrich case studies of DNA super-barcode and provide new insights on speciation, lineage diversification, and biogeography of the Fritillaria in China.

Molecular evolution of chloroplast genomes in Monsteroideae (Araceae)

This study provides broad insight into the chloroplast genomes of the subfamily Monsteroideae. The identified polymorphic regions may be suitable for designing unique and robust molecular markers for phylogenetic inference. Monsteroideae is the third largest subfamily (comprises 369 species) and one of the early diverging lineages of the monocot plant family Araceae. The phylogeny of this important subfamily is not well resolved at the species level due to scarcity of genomic resources and suitable molecular markers. Here, we report annotated chloroplast genome sequences of four Monsteroideae species: Spathiphyllum patulinervum, Stenospermation multiovulatum, Monstera adansonii, and Rhaphidophora amplissima. The quadripartite chloroplast genomes (size range 163,335-164,751 bp) consist of a pair of inverted repeats (25,270-25,931 bp), separating a small single copy region (21,448-22,346 bp) from a large single copy region (89,714-91,841 bp). The genomes contain 114 unique genes, including four rRNA genes, 80 protein-coding genes, and 30 tRNA genes. Gene features, amino acid frequencies, codon usage, GC contents, oligonucleotide repeats, and inverted repeats dynamics exhibit similarities among the four genomes. Higher rate of synonymous substitutions was observed as compared to non-synonymous substitutions in 76 protein-coding genes. Positive selection was observed in seven protein-coding genes, including psbK, ndhK, ndhD, rbcL, accD, rps8, and ycf2. Our included species of Araceae showed the monophyly in Monsteroideae and other subfamilies. We report 30 suitable polymorphic regions. The polymorphic regions identified here might be suitable for designing unique and robust markers for inferring the phylogeny and phylogeography among closely related species within the genus Spathiphyllum and among distantly related species within the subfamily Monsteroideae. The chloroplast genomes presented here are a valuable contribution towards understanding the molecular evolutionary dynamics in the family Araceae.

The complete chloroplast genome of greater duckweed (Spirodela polyrhiza 7498) using PacBio long reads: insights into the chloroplast evolution and transcription regulation

Background

Duckweeds (Lemnaceae) are aquatic plants distributed all over the world. The chloroplast genome, as an efficient solar-powered reactor, is an invaluable resource to study biodiversity and to carry foreign genes. The chloroplast genome sequencing has become routine and less expensive with the delivery of high-throughput sequencing technologies, allowing us to deeply investigate genomics and transcriptomics of duckweed organelles.

Results

Here, the complete chloroplast genome of Spirodela polyrhiza 7498 (SpV2) is assembled by PacBio sequencing. The length of 168,956 bp circular genome is composed of a pair of inverted repeats of 31,844 bp, a large single copy of 91,210 bp and a small single copy of 14,058 bp. Compared to the previous version (SpV1) assembled from short reads, the integrity and quality of SpV2 are improved, especially with the retrieval of two repeated fragments in ycf2 gene. There are a number of 107 unique genes, including 78 protein-coding genes, 25 tRNA genes and 4 rRNA genes. With the evidence of full-length cDNAs generated from PacBio isoform sequencing, seven genes (ycf3, clpP, atpF, rpoC1, rpl2, rps12 and ndhA) are detected to contain type-II introns. The ndhA intron has 50% more sequence divergence than the species-barcoding marker of atpF-atpH, showing the potential power to discriminate close species. A number of 37 RNA editing sites are recognized to have cytosine (C) to uracil (U) substitutions, eight of which are newly defined including six from the intergenic regions and two from the coding sequences of rpoC2 and ndhA genes. In addition, nine operon classes are identified using transcriptomic data. It is found that the operons contain multiple subunit genes encoding the same functional complexes comprising of ATP synthase, photosynthesis system, ribosomal proteins, et.al., which could be simultaneously transcribed and coordinately translated in response to the cell stimuli.

Conclusions

The understanding of the chloroplast genomics and the transcriptomics of S.polyrhiza would greatly facilitate the study of phylogenetic evolution and the application of genetically engineering duckweeds.