Insight into infrageneric circumscription through complete chloroplast genome sequences of two Trillium species

Seven Complete Chloroplast Genomes from Symplocos: Genome Organization and Comparative Analysis

In the present study, chloroplast genome sequences of four species of Symplocos (S. chinensis for. pilosa, S. prunifolia, S. coreana, and S. tanakana) from South Korea were obtained by Ion Torrent sequencing and compared with the sequences of three previously reported Symplocos chloroplast genomes from different species. The length of the Symplocos chloroplast genome ranged from 156,961 to 157,365 bp. Overall, 132 genes including 87 functional genes, 37 tRNA genes, and eight rRNA genes were identified in all Symplocos chloroplast genomes. The gene order and contents were highly similar across the seven species. The coding regions were more conserved than the non-coding regions, and the large single-copy and small single-copy regions were less conserved than the inverted repeat regions. We identified five new hotspot regions (rbcL, ycf4, psaJ, rpl22, and ycf1) that can be used as barcodes or species-specific Symplocos molecular markers. These four novel chloroplast genomes provide basic information on the plastid genome of Symplocos and enable better taxonomic characterization of this genus.

Comparative analysis and implications of the chloroplast genomes of three thistles (Carduus L., Asteraceae)

BACKGROUND

Carduus, commonly known as plumeless thistles, is a genus in the Asteraceae family that exhibits both medicinal value and invasive tendencies. However, the genomic data of Carduus (i.e., complete chloroplast genomes) have not been sequenced.

METHODS

We sequenced and assembled the chloroplast genome (cpDNA) sequences of three Carduus species using the Illumina Miseq sequencing system and Geneious Prime. Phylogenetic relationships between Carduus and related taxa were reconstructed using Maximum Likelihood and Bayesian Inference analyses. In addition, we used a single nucleotide polymorphism (SNP) in the protein coding region of the matK gene to develop molecular markers to distinguish C. crispus from C. acanthoides and C. tenuiflorus.

RESULTS

The cpDNA sequences of C. crispus, C. acanthoides, and C. tenuiflorus ranged from 152,342 bp to 152,617 bp in length. Comparative genomic analysis revealed high conservation in terms of gene content (including 80 protein-coding, 30 tRNA, and four rRNA genes) and gene order within the three focal species and members of subfamily Carduoideae. Despite their high similarity, the three species differed with respect to the number and content of repeats in the chloroplast genome. Additionally, eight hotspot regions, including psbI-trnS_GCU, trnE_UUC-rpoB, trnR_UCU-trnG_UCC, psbC-trnS_UGA, trnT_UGU-trnL_UAA, psbT-psbN, petD-rpoA, and rpl16-rps3, were identified in the study species. Phylogenetic analyses inferred from 78 protein-coding and non-coding regions indicated that Carduus is polyphyletic, suggesting the need for additional studies to reconstruct relationships between thistles and related taxa. Based on a SNP in matK, we successfully developed a molecular marker and protocol for distinguishing C. crispus from the other two focal species. Our study provides preliminary chloroplast genome data for further studies on plastid genome evolution, phylogeny, and development of species-level markers in Carduus.

Species Identification of Dracaena Using the Complete Chloroplast Genome as a Super-Barcode

The taxonomy and nomenclature of Dracaena plants are much disputed, particularly for several Dracaena species in Asia. However, neither morphological features nor common DNA regions are ideal for identification of Dracaena spp. Meanwhile, although multiple Dracaena spp. are sources of the rare traditional medicine dragon's blood, the Pharmacopoeia of the People's Republic of China has defined Dracaena cochinchinensis as the only source plant. The inaccurate identification of Dracaena spp. will inevitably affect the clinical efficacy of dragon's blood. It is therefore important to find a better method to distinguish these species. Here, we report the complete chloroplast (CP) genomes of six Dracaena spp., D. cochinchinensis, D. cambodiana, D. angustifolia, D. terniflora, D. hokouensis, and D. elliptica, obtained through high-throughput Illumina sequencing. These CP genomes exhibited typical circular tetramerous structure, and their sizes ranged from 155,055 (D. elliptica) to 155,449 bp (D. cochinchinensis). The GC content of each CP genome was 37.5%. Furthermore, each CP genome contained 130 genes, including 84 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. There were no potential coding or non-coding regions to distinguish these six species, but the maximum likelihood tree of the six Dracaena spp. and other related species revealed that the whole CP genome can be used as a super-barcode to identify these Dracaena spp. This study provides not only invaluable data for species identification and safe medical application of Dracaena but also an important reference and foundation for species identification and phylogeny of Liliaceae plants.

The implication of plastid transcriptome analysis in petaloid monocotyledons: A case study of Lilium lancifolium (Liliaceae, Liliales)

Transcriptome data provide useful information for studying the evolutionary history of angiosperms. Previously, different genomic events (i.e., duplication, deletion, and pseudogenization) were discovered in the plastid genome of Liliales; however, the effects of these events have not addressed because of the lack of transcriptome data. In this study, we completed the plastid genome (cpDNA) and generated transcriptome data of Lilium lancifolium. Consequently, the cpDNA of L. lancifolium is 152,479 bp in length, which consists of one large single copy (81,888 bp), one small single copy (17,607 bp), and two inverted repeat regions (26,544 bp). The comparative genomic analysis of newly sequenced cpDNA and transcriptome data revealed 90 RNA editing sites of which two positions are located in the rRNA coding region of L. lancifolium. A further check on the secondary structure of rRNA showed that RNA editing causes notable structural changes. Most of the RNA editing contents are C-to-U conversions, which result in nonsynonymous substitutions. Among coding regions, ndh genes have the highest number of RNA editing sites. Our study provided the first profiling of plastid transcriptome analyses in Liliales and fundamental information for further studies on post-transcription in this order as well as other petaloid monocotyledonous species.

The newly developed single nucleotide polymorphism (SNP) markers for a potentially medicinal plant, Crepidiastrum denticulatum (Asteraceae), inferred from complete chloroplast genome data

Medicinal effects of Crepidiastrum denticulatum have been previously reported. However, the genomic resources of this species and its applications have not been studied. In this study, based on the next generation sequencing method (Miseq sequencing system), we characterize the chloroplast genome of C. denticulatum which contains a large single copy (84,112 bp) and a small single copy (18,519 bp), separated by two inverted repeat regions (25,074 bp). This genome consists of 80 protein-coding gene, 30 tRNAs, and four rRNAs. Notably, the trnT_GGU is pseudogenized because of a small insertion within the coding region. Comparative genomic analysis reveals a high similarity among Asteraceae taxa. However, the junctions between LSC, SSC, and IRs locate in different positions within rps19 and ycf1 among examined species. Also, we describe a newly developed single nucleotide polymorphism (SNP) marker for C. denticulatum based on amplification-refractory mutation system (ARMS) technique. The markers, inferred from SNP in rbcL and matK genes, show effectiveness to recognize C. denticulatum from other related taxa through simple PCR protocol. The chloroplast genome-based molecular markers are effective to distinguish a potentially medicinal species, C. denticulatum, from other related taxa. Additionally, the complete chloroplast genome of C. denticulatum provides initial genomic data for further studies on phylogenomics, population genetics, and evolutionary history of Crepidiastrum as well as other taxa in Asteraceae.

Comparative chloroplast genomes of Paris Sect. Marmorata: insights into repeat regions and evolutionary implications

Background

Species of Paris Sect. Marmorata are valuable medicinal plants to synthesize steroidal saponins with effective pharmacological therapy. However, the wild resources of the species are threatened by plundering exploitation before the molecular genetics studies uncover the genomes and evolutionary significance. Thus, the availability of complete chloroplast genome sequences of Sect. Marmorata is necessary and crucial to the understanding the plastome evolution of this section and facilitating future population genetics studies. Here, we determined chloroplast genomes of Sect. Marmorata, and conducted the whole chloroplast genome comparison.

Results

This study presented detailed sequences and structural variations of chloroplast genomes of Sect. Marmorata. Over 40 large repeats and approximately 130 simple sequence repeats as well as a group of genomic hotspots were detected. Inverted repeat contraction of this section was inferred via comparing the chloroplast genomes with the one of P. verticillata. Additionally, almost all the plastid protein coding genes were found to prefer ending with A/U. Mutation bias and selection pressure predominately shaped the codon bias of most genes. And most of the genes underwent purifying selection, whereas photosynthetic genes experienced a relatively relaxed purifying selection.

Conclusions

Repeat sequences and hotspot regions can be scanned to detect the intraspecific and interspecific variability, and selected to infer the phylogenetic relationships of Sect. Marmorata and other species in subgenus Daiswa. Mutation and natural selection were the main forces to drive the codon bias pattern of most plastid protein coding genes. Therefore, this study enhances the understanding about evolution of Sect. Marmorata from the chloroplast genome, and provide genomic insights into genetic analyses of Sect. Marmorata.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5281-x) contains supplementary material, which is available to authorized users.

Reconfiguration of the plastid genome in Lamprocapnos spectabilis: IR boundary shifting, inversion, and intraspecific variation

We generated a complete plastid genome (plastome) sequence for Lamprocapnos spectabilis, providing the first complete plastome from the subfamily Fumarioideae (Papaveraceae). The Lamprocapnos plastome shows large differences in size, structure, gene content, and substitution rates compared with two sequenced Papaveraceae plastomes. We propose a model that explains the major rearrangements observed, involving at least six inverted repeat (IR) boundary shifts and five inversions, generating a number of gene duplications and relocations, as well as a two-fold expansion of the IR and miniaturized small single-copy region. A reduction in the substitution rates for genes transferred from the single-copy regions to the IR was observed. Accelerated substitution rates of plastid accD and clpP were detected in the Lamprocapnos plastome. The accelerated substitution rate for the accD gene was correlated with a large insertion of amino acid repeat (AAR) motifs in the middle region, but the forces driving the higher substitution rate of the clpP gene are unclear. We found a variable number of AARs in Lamprocapnos accD and ycf1 genes within individuals, and the repeats were associated with coiled-coil regions. In addition, comparative analysis of three Papaveraceae plastomes revealed loss of rps15 in Papaver, and functional replacement to the nucleus was identified.

Molecular Phylogeny and Dating of Forsythieae (Oleaceae) Provide Insight into the Miocene History of Eurasian Temperate Shrubs

Tribe Forsythieae (Oleaceae), containing two genera (Abeliophyllum and Forsythia) and 13 species, is economically important plants used as ornamentals and in traditional medicine. This tribe species occur primarily in mountainous regions of Eurasia with the highest species diversity in East Asia. Here, we examine 11 complete chloroplast genome and nuclear cycloidea2 (cyc2) DNA sequences of 10 Forsythia species and Abeliophyllum distichum using Illumina platform to provide the phylogeny and biogeographic history of the tribe. The chloroplast genomes of the 11 Forsythieae species are highly conserved, except for a deletion of about 400 bp in the accD-psaI region detected only in Abeliophyllum. Within Forsythieae species, analysis of repetitive sequences revealed a total of 51 repeats comprising 26 forward repeats, 22 palindromic repeats, and 3 reverse repeats. Of those, 19 repeats were common and 32 were unique to one or more Forsythieae species. Our phylogenetic analyses supported the monophyly of Forsythia and its sister group is Abeliophyllum using the concatenated dataset of 78 chloroplast genes. Within Forsythia, Forsythia likiangensis and F. giraldiana were basal lineages followed by F. europaea; the three species are characterized by minutely serrate or entire leaf margins. The remaining species, which are distributed in East Asia, formed two major clades. One clade included F. ovata, F. velutina, and F. japonica; they are morphologically supported by broadly ovate leaves. Another clade of F. suspensa, F. saxatilis, F. viridissima, and F. koreana characterized by lanceolate leaves (except F. suspensa which have broad ovate leaves). Although cyc2 phylogeny is largely congruent to chloroplast genome phylogeny, we find the discordance between two phylogenies in the position of F. ovata suggesting that introgression of the chloroplast genome from one species into the nuclear background of another by interspecific hybridization in East Asian Forsythia species. Molecular dating and biogeographic reconstructions suggest an origin of the Forsythieae species in East China in the Miocene. Distribution patterns in Forsythia indicated that the species were radially differentiated from East China, and the speciation of the European F. europaea was the result of both vicariance and dispersal in the late Miocene to Pliocene.

Complete structure and variation of the chloroplast genome of Agropyron cristatum (L.) Gaertn

Agropyron cristatum (L.) Gaertner, a perennial grass in the tribe Triticeae (Poaceae), is a wild relative of cereal crops that is suitable for genetic improvement. In this study, we first sequenced the complete chloroplast (cp) genome of Ag. cristatum using Hiseq4000 PE150. The Ag. cristatum chloroplast genome is 135,554bp in length, has a typical quadripartite structure and contains 76 protein-coding genes, 29 tRNA genes and four rRNA genes. The cp genome of Ag. cristatum was used for comparison with other seven Triticeae species. One large variable region (800bp), which primarily contained the rpl23 (non-reciprocally translocated from IRs) and accD genes, was detected between rbcL gene and psaI gene within LSC region. The deletion of the accD and translocated rpl23 genes in Ag. cristatum indicated an independent gene-loss events or an additional divergence in Triticeae. Analyses of the dn/ds ratio and K2-P's genetic distance for 76 protein-coding genes showed that genes with evolutionary divergence might suffer from the effect of sequence regional constraints or gene functional constraints in Triticeae species. Our research will generally contribute to the knowledge of plastid genome evolution in Triticeae.

Chloroplast Genomic Resource of Paris for Species Discrimination

Paris is famous in China for its medicinal value and has been included in the Chinese Pharmacopoeia. Inaccurate identification of these species could confound their effective exploration, conservation, and domestication. Due to the plasticity of the morphological characteristics, correct identification among Paris species remains problematic. In this regard, we report the complete chloroplast genome of P. thibetica and P. rugosa to develop highly variable molecular markers. Comparing three chloroplast genomes, we sought out the most variable regions to develop the best cpDNA barcodes for Paris. The size of Paris chloroplast genome ranged from 162,708 to 163,200 bp. A total of 134 genes comprising 81 protein coding genes, 45 tRNA genes and 8 rRNA genes were observed in all three chloroplast genomes. Eight rapidly evolving regions were detected, as well as the difference of simple sequence repeats (SSR) and repeat sequence. Two regions of the coding gene ycf1, ycf1a and ycf1b, evolved the quickest and were proposed as core barcodes for Paris. The complete chloroplast genome sequences provide more integrated and adequate information for better understanding the phylogenetic pattern and improving efficient discrimination during species identification.

A Dynamic Tandem Repeat in Monocotyledons Inferred from a Comparative Analysis of Chloroplast Genomes in Melanthiaceae

Chloroplast genomes (cpDNA) are highly valuable resources for evolutionary studies of angiosperms, since they are highly conserved, are small in size, and play critical roles in plants. Slipped-strand mispairing (SSM) was assumed to be a mechanism for generating repeat units in cpDNA. However, research on the employment of different small repeated sequences through SSM events, which may induce the accumulation of distinct types of repeats within the same region in cpDNA, has not been documented. Here, we sequenced two chloroplast genomes from the endemic species Heloniopsis tubiflora (Korea) and Xerophyllum tenax (USA) to cover the gap between molecular data and explore "hot spots" for genomic events in Melanthiaceae. Comparative analysis of 23 complete cpDNA sequences revealed that there were different stages of deletion in the rps16 region across the Melanthiaceae. Based on the partial or complete loss of rps16 gene in cpDNA, we have firstly reported potential molecular markers for recognizing two sections (Veratrum and Fuscoveratrum) of Veratrum. Melathiaceae exhibits a significant change in the junction between large single copy and inverted repeat regions, ranging from trnH_GUG to a part of rps3. Our results show an accumulation of tandem repeats in the rpl23-ycf2 regions of cpDNAs. Small conserved sequences exist and flank tandem repeats in further observation of this region across most of the examined taxa of Liliales. Therefore, we propose three scenarios in which different small repeated sequences were used during SSM events to generate newly distinct types of repeats. Occasionally, prior to the SSM process, point mutation event and double strand break repair occurred and induced the formation of initial repeat units which are indispensable in the SSM process. SSM may have likely occurred more frequently for short repeats than for long repeat sequences in tribe Parideae (Melanthiaceae, Liliales). Collectively, these findings add new evidence of dynamic results from SSM in chloroplast genomes which can be useful for further evolutionary studies in angiosperms. Additionally, genomics events in cpDNA are potential resources for mining molecular markers in Liliales.

Analysis of Complete Chloroplast Genome Sequences Improves Phylogenetic Resolution in Paris (Melanthiaceae)

The genus Paris in the broad concept is an economically important group in the monocotyledonous family Melanthiaceae (tribe Parideae). The phylogeny of Paris was controversial in previous morphology-based classification and molecular phylogeny. Here, the complete cp genomes of eleven Paris taxa were sequenced, to better understand the evolutionary relationships among these plants and the mutation patterns in their chloroplast (cp) genomes. Comparative analyses indicated that the overall cp genome structure among the Paris taxa is quite similar. The triplication of trnI-CAU was found only in the cp genomes of P. quadrifolia and P. verticillata. Phylogenetic analyses based on the complete cp genomes did not resolve Paris as a monophyletic group, instead providing evidence supporting division of the twelve taxa into two segregate genera: Paris sensu strict and Daiswa. The sister relationship between Daiswa and Trillium was well supported. We recovered two fully supported lineages with divergent distribution in Daiswa; however, none of the previously recognized sections in Daiswa was resolved as monophyletic using plastome data, suggesting that the infrageneric relationships and biogeography of Daiswa species require further investigation. Ten highly divergent DNA regions, suitable for species identification, were detected among the 12 cp genomes. This study is the first successful attempt to provide well-supported evolutionary relationships in Paris based on phylogenomic analyses. The findings highlight the potential of the whole cp genomes for improving resolution in phylogeny as well as species identification in phylogenetically and taxonomically difficult plant genera.