The complete chloroplast genome sequence of watercress (Nasturtium officinale R. Br.): Genome organization, adaptive evolution and phylogenetic relationships in Cardamineae

Comprehensive Analysis of the Complete Chloroplast Genome of Cinnamomum daphnoides (Lauraceae), An Endangered Island Endemic Plant

Cinnamomum daphnoides (Siebold & Zucc 1846) is a rare and endangered island species with a unique Sino-Japanese distribution pattern. However, inormation regarding the species' chloroplast (cp) genome, structural features, and the phylogenetic relationship is still lacking. We utilized high-throughput sequencing technology to assemble and annotate the first cp genome of C. daphnoides (GenBank OR654104), followed by genomic characterization and phylogenetic analysis to fill the gaps in this species' cp genome. Our analysis showed that the cp genome has a quadripartite structure spanning 152,765 bp with a GC content of 39.15%. The genome encodes 126 genes, which include 36 tRNA genes, 8 rRNA genes, and 82 mRNA genes. Specifically, 44 genes are related to photosynthesis, 59 are associated with self-replication, six are other genes, and four have unknown functionality. The Codon usage bias in the genome exhibits a preference for A/U bases. We identified 29 interspaced repeat sequences that belonging to three types of repeat sequences. A total of 217 cpSSR loci were detected with single nucleotide repeats (59.91%) being the most frequent loci, mainly composed of A/T repeats. Our selection pressure analysis revealed that the ycf2 gene experienced strong positive selection (Ka/Ks = 1.81, P > 0.844). Further, we identified three highly variable fragments (psbM, psbT, and ycf1) that can be utilized as specific DNA barcoding markers for species definition and population genetic studies. We conducted boundary analysis, which showed that the structure and gene sequence of the two species were highly conserved. Finally, our phylogenetic analysis supports that C. daphnoides is close to C. cassia in the Cinnamomum genes, indicating that the two species share a common ancestry. Overall, providing genomic information on C. daphnoides will be beneficial for the conservation and utilization of endangered plant genetic resources. It will also serve as a reference for the identification of species and the phylogenetic analysis of Cinnamomum. This information will be useful in future research.

Complete chloroplast genome sequence of Artemisia littoricola (Asteraceae) from Dokdo Island Korea: genome structure, phylogenetic analysis, and biogeography study

The Asteraceae family, particularly the Artemisia genus, presents taxonomic challenges due to limited morphological characteristics and frequent natural hybridization. Molecular tools, such as chloroplast genome analysis, offer solutions for accurate species identification. In this study, we sequenced and annotated the chloroplast genome of Artemisia littoricola sourced from Dokdo Island, employing comparative analyses across six diverse Artemisia species. Our findings reveal conserved genome structures with variations in repeat sequences and junction boundaries. Notably, the chloroplast genome of A. littoricola spans 150,985 bp, consistent with other Artemisia species, and comprises 131 genes, including 86 protein-coding, 37 tRNA, and 8 rRNA genes. Among these genes, 16 possess a single intron, while clpP and ycf3 exhibit two introns each. Furthermore, 18 genes display duplicated copies within the IR regions. Moreover, the genome possesses 42 Simple Sequence Repeats (SSRs), predominantly abundant in A/T content and located within intergenic spacer regions. The analysis of codon usage revealed that the codons for leucine were the most frequent, with a preference for ending with A/U. While the chloroplast genome exhibited conservation overall, non-coding regions showed lower conservation compared to coding regions, with the Inverted Repeat (IR) region displaying higher conservation than single-copy regions. Phylogenetic analyses position A. littoricola within subgenus Dracunculus, indicating a close relationship with A. scoparia and A. desertorum. Additionally, biogeographic reconstructions suggest ancestral origins in East Asia, emphasizing Mongolia, China (North East and North Central and South Central China), and Korea. This study underscores the importance of chloroplast genomics in understanding Artemisia diversity and evolution, offering valuable insights into taxonomy, evolutionary patterns, and biogeographic history. These findings not only enhance our understanding of Artemisia's intricate biology but also contribute to conservation efforts and facilitate the development of molecular markers for further research and applications in medicine and agriculture.

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.

The chloroplast genome of Chrozophora sabulosa Kar. & Kir. and its exploration in the evolutionary position uncertainty of genus Chrozophora

Chrozophora sabulosa Kar. & Kir. is a biennial herbaceous plant that belongs to the Euphorbiaceae family and has medicinal properties. This research aimed to identify the genetic characteristics and phylogenetic position of the Chrozophora genus within the Euphorbiaceae family. The evolutionary position of the Chrozophora genus was previously unknown due to insufficient research. Therefore, to determine the evolutionary link between C. sabulosa and other related species, we conducted a study using the NGS Illumina platform to sequence the C. sabulosa chloroplast (cp.) genome. The study results showed that the genome was 156,488 bp in length. It had a quadripartite structure consisting of two inverted repeats (IRb and IRa) of 24,649-bp, separated by an 87,696-bp LSC region and a 19,494-bp SSC region. The CP genome contained 113 unique genes, including four rRNA genes, 30 tRNA genes, and 79 CDS genes. In the second copy of the inverted repeat, there were 18 duplicated genes. The C. sabulosa lacks the petD, petB, rpl2, and rps16 intron. The analysis of simple sequence repeats (SSRs) revealed 93 SSR loci of 22 types and 78 oligonucleotide repeats of four kinds. The phylogenetic investigation showed that the Chrozophora genus evolved paraphyletically from other members of the Euphorbiaceae family. To support the phylogenetic findings, we selected species from the Euphorbiaceae and Phyllanthaceae families to compare with C. sabulosa for Ks and Ka substitution rates, InDels investigation, IR contraction and expansion, and SNPs analysis. The results of these comparative studies align with the phylogenetic findings. We identified six highly polymorphic regions shared by both families, which could be used as molecular identifiers for the Chrozophora genus (rpl33-rps18, rps18-rpl20, rps15-ycf1, ndhG-ndhI, psaI-ycf4, petA-psbJ). The cp. genome sequence of C. sabulosa reveals the evolution of plastid sequences in Chrozophora species. This is the first time the cp. genome of a Chrozophora genus has been sequenced, serving as a foundation for future sequencing of other species within the Chrozophoreae tribe and facilitating in-depth taxonomic research. The results of this research will also aid in identifying new Chrozophora species.

Genomic exploration of Sesuvium sesuvioides: comparative study and phylogenetic analysis within the order Caryophyllales from Cholistan desert, Pakistan

BACKGROUND

The Aizoaceae family's Sesuvium sesuvioides (Fenzl) Verdc is a medicinal species of the Cholistan desert, Pakistan. The purpose of this study was to determine the genomic features and phylogenetic position of the Sesuvium genus in the Aizoaceae family. We used the Illumina HiSeq2500 and paired-end sequencing to publish the complete chloroplast sequence of S. sesuvioides.

RESULTS

The 155,849 bp length cp genome sequence of S. sesuvioides has a 36.8% GC content. The Leucine codon has the greatest codon use (10.6%), 81 simple sequence repetitions of 19 kinds, and 79 oligonucleotide repeats. We investigated the phylogeny of the order Caryophyllales' 27 species from 23 families and 25 distinct genera. The maximum likelihood tree indicated Sesuvium as a monophyletic genus, and sister to Tetragonia. A comparison of S. sesuvioides, with Sesuvium portulacastrum, Mesembryanthemum crystallinum, Mesembryanthemum cordifolium, and Tetragonia tetragonoides was performed using the NCBI platform. In the comparative investigation of genomes, all five genera revealed comparable cp genome structure, gene number and composition. All five species lacked the rps15 gene and the rpl2 intron. In most comparisons with S. sesuvioides, transition substitutions (Ts) were more frequent than transversion substitutions (Tv), producing Ts/Tv ratios larger than one, and the Ka/Ks ratio was lower than one. We determined ten highly polymorphic regions, comprising rpl22, rpl32-trnL-UAG, trnD-GUC-trnY-GUA, trnE-UUC-trnT-GGU, trnK-UUU-rps16, trnM-CAU-atpE, trnH-GUG-psbA, psaJ-rpl33, rps4-trnT-UGU, and trnF-GAA-ndhJ.

CONCLUSION

The whole S. sesuvioides chloroplast will be examined as a resource for in-depth taxonomic research of the genus when more Sesuvium and Aizoaceae species are sequenced in the future. The chloroplast genomes of the Aizoaceae family are well preserved, with little alterations, indicating the family's monophyletic origin. This study's highly polymorphic regions could be utilized to build realistic and low-cost molecular markers for resolving taxonomic discrepancies, new species identification, and finding evolutionary links among Aizoaceae species. To properly comprehend the evolution of the Aizoaceae family, further species need to be sequenced.

Comparative and phylogenetic analyses of nine complete chloroplast genomes of Orchidaceae

The orchid family has 200,000 species and 700 genera, and it is found worldwide in the tropics and subtropics. In China, there are 1247 species and subspecies of orchids belonging to the Orchidaceae family. Orchidaceae is one of the most diverse plant families in the world, known for their lush look, remarkable ecological tolerance, and capability for reproduction. It has significant decorative and therapeutic value. In terms of evolution, the orchid family is one of the more complicated groups, but up until now, little has been known about its affinities. This study examined the properties of 19 chloroplast (cp) genomes, of which 11 had previously been published and nine had only recently been revealed. Following that, topics such as analysis of selection pressure, codon usage, amino acid frequencies, repeated sequences, and reverse repeat contraction and expansion are covered. The Orchidaceae share similar cp chromosomal characteristics, and we have conducted a preliminary analysis of their evolutionary connections. The cp genome of this family has a typical tepartite structure and a high degree of consistency across species. Platanthera urceolata with more tandem repeats of the cp genome. Similar cp chromosomal traits can be seen in the orchidaceae. Galearis roborowskyi, Neottianthe cucullata, Neottianthe monophylla, Platanthera urceolata and Ponerorchis compacta are the closest cousins, according to phylogenetic study.

Comparative analysis of chloroplast genome and new insights into phylogenetic relationships of Ajuga and common adulterants

Introduction

The potential contamination of herbal medicinal products poses a significant concern for consumer health. Given the limited availability of genetic information concerning Ajuga species, it becomes imperative to incorporate supplementary molecular markers to enhance and ensure accurate species identification.

Methods

In this study, the chloroplast (cp) genomes of seven species of the genus Ajuag were sequenced, de novo assembled and characterized.

Results

exhibiting lengths ranging from 150,342 bp to 150,472 bp, encompassing 86 - 88 protein-coding genes (PCGs), 35 - 37 transfer RNA, and eight ribosomal RNA. The repetitive sequences, codon uses, and cp genomes of seven species were highly conserved, and PCGs were the reliable molecular markers for investigating the phylogenetic relationship within the Ajuga genus. Moreover, four mutation hotspot regions (accD-psaI, atpH-atpI, ndhC-trnV(UAC), and ndhF-rpl23) were identified within cp genomes of Ajuga, which could help distinguish A. bracteosa and its contaminants. Based on cp genomes and PCGs, the phylogenetic tree preliminary confirmed the position of Ajuga within the Lamiaceae family. It strongly supported a sister relationship between Subsect. Genevense and Subsect. Biflorae, suggesting the merger of Subsect. Biflorae and Subsect. Genevenses into one group rather than maintaining separate categorizations. Additionally, molecular clock analysis estimated the divergence time of Ajuga to be around 7.78 million years ago.

Discussion

The species authentication, phylogeny, and evolution analyses of the Ajuga species may benefit from the above findings.

A Comprehensive Analysis of Chloroplast Genome Provides New Insights into the Evolution of the Genus Chrysosplenium

Chrysosplenium, a perennial herb in the family Saxifragaceae, prefers to grow in low light and moist environments and is divided into two sections of Alternifolia and Oppositifolia based on phyllotaxy. Although there has been some progress in the phylogeny of Chrysosplenium over the years, the phylogenetic position of some species is still controversial. In this study, we assembled chloroplast genomes (cp genomes) of 34 Chrysosplenium species and performed comparative genomic and phylogenetic analyses in combination with other cp genomes of previously known Chrysosplenium species, for a total of 44 Chrysosplenium species. The comparative analyses revealed that cp genomes of Chrysosplenium species were more conserved in terms of genome structure, gene content and arrangement, SSRs, and codon preference, but differ in genome size and SC/IR boundaries. Phylogenetic analysis showed that cp genomes effectively improved the phylogenetic support and resolution of Chrysosplenium species and strongly supported Chrysosplenium species as a monophyletic taxon and divided into three branches. The results also showed that the sections of Alternifolia and Oppositifolia were not monophyletic with each other, and that C. microspermum was not clustered with other Chrysosplenium species with alternate leaves, but with C. sedakowii into separate branches. In addition, we identified 10 mutational hotspot regions that could serve as potential DNA barcodes for Chrysosplenium species identification. In contrast to Peltoboykinia, the clpP and ycf2 genes of Chrysosplenium were subjected to positive selection and had multiple significant positive selection sites. We further detected a significant positive selection site on the petG gene between the two sections of Chrysosplenium. These evolutionary characteristics may be related to the growth environment of Chrysosplenium species. This study enriches the cp genomes of Chrysosplenium species and provides a reference for future studies on its evolution and origin.

Comparative analysis of Diospyros (Ebenaceae) plastomes: Insights into genomic features, mutational hotspots, and adaptive evolution

Diospyros (Ebenaceae) is a widely distributed genus of trees and shrubs from pantropical to temperate regions, with numerous species valued for their fruits (persimmons), timber, and medicinal values. However, information regarding their plastomes and chloroplast evolution is scarce. The present study performed comparative genomic and evolutionary analyses on plastomes of 45 accepted Diospyros species, including three newly sequenced ones. Our study showed a highly conserved genomic structure across the Diospyros species, with 135-136 encoding genes, including 89 protein-coding genes, 1-2 pseudogenes (Ψycf1 for all, Ψrps19 for a few), 37 tRNA genes and 8 rRNA genes. Comparative analysis of Diospyros identified three intergenic regions (ccsA-ndhD, rps16-psbK and petA-psbJ) and five genes (rpl33, rpl22, petL, psaC and rps15) as the mutational hotspots in these species. Phylogenomic analysis identified the phylogenetic position of three newly sequenced ones and well supported a monophylogenetic (sub)temperate taxa and four clades in the pantropical taxa. The analysis codon usage identified 30 codons with relative synonymous codon usage (RSCU) values >1 and 29 codons ending with A and U bases. A total of three codons (UUA, GCU, and AGA) with highest RSCU values were identified as the optimal codons. Effective number of codons (ENC)-plot indicated the significant role of mutational pressure in shaping codon usage, while most protein-coding genes in Diospyros experienced relaxed purifying selection (d N/d S < 1). Additionally, the psbH gene showed positive selection (d N/d S > 1) in the (sub)temperate species. Thus, the results provide a meaningful foundation for further elaborating Diospyros's genetic architecture and taxonomy, enriching genetic diversity and conserving genetic resources.

The complete chloroplast genome sequence of Eugenia klotzschiana O. Berg unveils the evolutionary dynamics in plastomes of Myrteae DC. tribe (Myrtaceae)

Myrteae is the most diversified tribe in the Myrtaceae family and has great ecological and economic importance. Here, we performed the assembly and annotation of the chloroplast genome of Eugenia klotzschiana O. Berg and used this in a comparative analysis with other 13 species from the Myrteae tribe. The E. klotzschiana plastome exhibited a length of 158,977 bp and a very conserved structure and gene composition when compared with other Myrteae genomes. We identified 34 large repetitive sequences and 94 SSR repeats in E. klotzschiana plastome. The trnT-trnL, rpl32-trnL, ndhF-rpl32, psbE-petL, and ycf1 regions were identified as mutational hotspots. A negative selection signal was detected in 74 protein-coding genes while neutral selection was detected in two genes (rps12 and psaI). Furthermore, 222 RNA editing sites were identified in the E. klotzschiana plastome. We also obtained a plastome-based Myrtales phylogenetic tree, including E. klotzschiana for the first time in a molecular phylogeny, recovering its sister relationship for all other Eugenia species. Our results illuminate how evolution shaped the chloroplast genome structure and composition in the Myrteae tribe, especially in the E. klotzschiana plastome.

Genetic relationship of Pleione based on the chloroplast genome

Pleione (Orchidaceae) is a very famous horticultural plant with a high international reputation for its unique flower shape and abundant color. The small difference in morphological characteristics among Pleione species caused by weak reproductive isolation and easy hybridization makes the taxonomic status of individual species very confusing. Chloroplast (cp) genome information is of great significance for the study of plant phylogeny and taxonomy. In this study, the cp genomes of Pleione were sequenced and the complete cp structure and sequence characteristics of 19 species were compared and analyzed. The cp genome of Pleione species exhibited a conserved tetrad structure and each species encoded 135 protein-coding genes, 38 tRNA and 8 RNA genes. The cp genome sizes of 19 Pleione were 157964-159269 bp and the length of LSC, SSC and IR were 85953-87172 bp, 18499-18712 bp, 26459-26756 bp, respectively. Palindromic and forward repeats accounted for a high proportion and the SSRs were mainly mononucleotide repeats in Pleione. Analysis of chloroplast sequence differences indicated that the differences in coding regions were smaller than those in non-coding regions, and the variation in LSC and SSC regions was greater than that in IR regions. Phylogenetic analysis showed that all Pleione species inferred from the cp genome were clustered together and received high support. However, the genetic relationship of Pleione plants is different from the current update system of this genus. Therefore, the demarcation of Pleione interspecific relationships still needs further investigation due to the lack of sufficient evidence. The cp genome serves as valuable information for the identification of Pleione species and the study of phylogenetic relationships.

The complete chloroplast genome of Hibiscus syriacus using long-read sequencing: Comparative analysis to examine the evolution of the tribe Hibisceae

Hibiscus syriacus, a member of the tribe Hibisceae, is considered an important ornamental and medicinal plant in east Asian countries. Here, we sequenced and assembled the complete chloroplast genome of H. syriacus var. Baekdansim using the PacBio long-read sequencing platform. A quadripartite structure with 161,026 base pairs was obtained, consisting of a pair of inverted repeats (IRA and IRB) with 25,745 base pairs, separated by a large single-copy region of 89,705 base pairs and a short single-copy region of 19,831 base pairs. This chloroplast genome had 79 protein-coding genes, 30 transfer RNA genes, 4 ribosomal RNA genes, and 109 simple sequence repeat regions. Among them, ndhD and rpoC1, containing traces of RNA-editing events associated with adaptive evolution, were identified by analysis of putative RNA-editing sites. Codon usage analysis revealed a preference for A/U-terminated codons. Furthermore, the codon usage pattern had a clustering tendency similar to that of the phylogenetic analysis of the tribe Hibisceae. This study provides clues for understanding the relationships and refining the taxonomy of the tribe Hibisceae.

Plastid Phylogenomic Analyses Reveal the Taxonomic Position of Peucedanum franchetii

Peucedanum franchetii is a famous folk medicinal plant in China. However, the taxonomy of the P. franchetii has not been sufficiently resolved. Due to similar morphological features between P. franchetii and Ligusticopsis members, the World Flora Online (WFO) Plant List suggested that this species transformed into the genus Ligusticopsis and merged with Ligusticopsis likiangensis. However, both species are obviously diverse in leaf shape, bracts, and bracteoles. To check the taxonomic position of P. franchetii, we newly sequenced and assembled the plastome of P. franchetii and compared it with nine other plastomes of the genus Ligusticopsis. Ten plastomes were highly conserved and similar in gene order, codon bias, RNA editing sites, IR borders, and SSRs. Nevertheless, 10 mutation hotspot regions (infA, rps8, matK, ndhF, rps15, psbA-trnH, rps2-rpoC2, psbA-trnK, ycf2-trnL, and ccsA-ndhD) were still detected. In addition, both phylogenetic analyses based on plastome data and ITS sequences robustly supported that P. franchetii was not clustered with members of Peucedanum but nested in Ligusticopsis. P. franchetii was sister to L. likiangensis in the ITS topology but clustered with L. capillacea in the plastome tree. These findings implied that P. franchetii should be transferred to genus Ligusticopsis and not merged with L. likiangensis, but as an independent species, which was further verified by morphological evidences. Therefore, transferring P. franchetii under the genus Ligusticopsis as an independent species was reasonable, and a new combination was presented.

Characterization of the Complete Chloroplast Genome of Four Species in Callerya

BACKGROUND

Callerya reticulata (Bentham) Schot, Callerya dielsiana (Harms) P.K. Loc ex Z. Wei & Pedley, Callerya nitida var. hirsutissima (Z. Wei) X.Y. Zhu, and Callerya nitida (Bentham) R. Geesink, which belongs to the Leguminosae family, are important medicinal plants in China. The genus Callerya includes 26 species, 18 species are distributed in China, and the vine stems of some species are used as traditional medicinal herbs because they have important pharmacological activity. Due to the high similarity of appearance, it is difficult to identify them in the market by appearance alone. Therefore, circulating of Callerya-related materia medica on the market is confusing, sometimes even leading to drug safety problems. It is urgent to develop molecular methods for their identification.

OBJECTIVE

To sequence and analyze the complete chloroplast (cp) genomes of C. reticulata, C. dielsiana, C. nitida var. hirsutissima, and C. nitida and to analyze their cp genome differences as a basis for seeking easier DNA barcoding for their identification.

METHOD

After using Illumina high-throughput sequencing and nanopore sequencing to obtain the genome data, some bioinformatics software was used to assembly and analyze the molecular structure of cp genomes.

RESULTS

The complete cp genomes of the four species were circular molecules, which ranged from 130 435 to 132 546 bp, and GC contents ranged from 33.89% to 34.89%. Each of them includes a large single-copy region, a small single-copy region, and without large inverted repeat regions.

CONCLUSIONS

These results suggested that highly variable regions of the four cp genomes would provide useful plastid markers, which could be used as a potential genomic resource to resolve phylogenetic questions and provide a reference for mining specific DNA barcodes of these species.

HIGHLIGHTS

Our study provided highly effective molecular markers for subsequent phylogenetic analysis, species identification, and biogeographic analysis of Callerya.

The entire chloroplast genome sequence of Asparagus setaceus (Kunth) Jessop: Genome structure, gene composition, and phylogenetic analysis in Asparagaceae

Asparagus setaceus (Kunth) Jessop is a horticultural plant of the genus Asparagus. Herein, the whole chloroplast (cp) genome of A. setaceus was sequenced with PacBio and Illumina sequencing systems. The cp genome shows a characteristic quadripartite structure with 158,076 bp. In total, 135 genes were annotated, containing 89 protein-coding, 38 tRNA, and 8 rRNA genes. Contrast with the previous cp genome of A. setaceus registered in NCBI, we identified 7 single-nucleotide polymorphisms and 15 indels, mostly situated in noncoding areas. Meanwhile, 36 repeat structures and 260 simple sequence repeats were marked out. A bias for A/T-ending codons was shown in this cp genome. Furthermore, we predicted 78 RNA-editing sites in 29 genes, which were all for C-to-U transitions. And it was also proven that positive selection was exerted on the rpoC1 gene of A. setaceus with the K a/K s data. Meanwhile, a conservative gene order and highly similar sequences of protein-coding genes were revealed within Asparagus species. Phylogenetic tree analysis indicated that A. setaceus was a sister to Asparagus cochinchinensis. Taken together, our released genome provided valuable information for the gene composition, genetics comparison, and the phylogeny studies of A. setaceus.

Characterization of the Evolutionary Pressure on Anisodus tanguticus Maxim. with Complete Chloroplast Genome Sequence

Anisodus tanguticus Maxim. (Solanaceae), a traditional endangered Tibetan herb, is endemic to the Qinghai-Tibet Plateau. Here, we report the de novo assembled chloroplast (cp) genome sequences of A. tanguticus (155,765 bp). The cp contains a pair of inverted repeated (IRa and IRb) regions of 25,881 bp that are separated by a large single copy (LSC) region (86,516 bp) and a small single copy SSC (17,487 bp) region. A total of 132 functional genes were annotated in the cp genome, including 87 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Moreover, 199 simple sequence repeats (SSR) and 65 repeat structures were detected. Comparative plastome analyses revealed a conserved gene order and high similarity of protein-coding sequences. The A. tanguticus cp genome exhibits contraction and expansion, which differs from Przewalskia tangutica and other related Solanaceae species. We identified 30 highly polymorphic regions, mostly belonging to intergenic spacer regions (IGS), which may be suitable for the development of robust and cost-effective markers for inferring the phylogeny of the genus Anisodus and family Solanaceae. Analysis of the Ka/Ks ratios of the Hyoscyameae tribe revealed significant positive selection exerted on the cemA, rpoC2, and clpP genes, which suggests that protein metabolism may be an important strategy for A. tanguticus and other species in Hyoscyameae in adapting to the adverse environment on the Qinghai-Tibetan Plateau. Phylogenetic analysis revealed that A. tanguticus clustered closer with Hyoscyamus niger than P. tangutica. Our results provide reliable genetic information for future exploration of the taxonomy and phylogenetic evolution of the Hyoscyameae tribe and related species.

Comparative Analysis of Complete Chloroplast Genomes of Nine Species of Litsea (Lauraceae): Hypervariable Regions, Positive Selection, and Phylogenetic Relationships

Litsea is a group of evergreen trees or shrubs in the laurel family, Lauraceae. Species of the genus are widely used for a wide range of medicinal and industrial aspects. At present, most studies related to the gene resources of Litsea are restricted to morphological analyses or features of individual genomes, and currently available studies of select molecular markers are insufficient. In this study, we assembled and annotated the complete chloroplast genomes of nine species in Litsea, carried out a series of comparative analyses, and reconstructed phylogenetic relationships within the genus. The genome length ranged from 152,051 to 152,747 bp and a total of 128 genes were identified. High consistency patterns of codon bias, repeats, divergent analysis, single nucleotide polymorphisms (SNP) and insertions and deletions (InDels) were discovered across the genus. Variations in gene length and the presence of the pseudogene ycf1^Ψ, resulting from IR contraction and expansion, are reported. The hyper-variable gene rpl16 was identified for its exceptionally high Ka/Ks and Pi values, implying that those frequent mutations occurred as a result of positive selection. Phylogenetic relationships were recovered for the genus based on analyses of full chloroplast genomes and protein-coding genes. Overall, both genome sequences and potential molecular markers provided in this study enrich the available genomic resources for species of Litsea. Valuable genomic resources and divergent analysis are also provided for further research of the evolutionary patterns, molecular markers, and deeper phylogenetic relationships of Litsea.

The chloroplast genome of Salix floderusii and characterization of chloroplast regulatory elements

Salix floderusii is a rare alpine tree species in the Salix genus. Unfortunately, no extensive germplasm identification, molecular phylogeny, and chloroplast genomics of this plant have been conducted. We sequenced the chloroplast (cp) genome of S. floderusii for the first time using second-generation sequencing technology. The cp genome was 155,540 bp long, including a large single-copy region (LSC, 84,401 bp), a small single-copy region (SSC, 16,221 bp), and inverted repeat regions (IR, 54,918 bp). A total of 131 genes were identified, including 86 protein genes, 37 tRNA genes, and 8 rRNA genes. The S. floderusii cp genome contains 1 complement repeat, 24 forward repeats, 17 palindromic repeats, and 7 reverse repeats. Analysis of the IR borders showed that the IRa and IRb regions of S. floderusii and Salix caprea were shorter than those of Salix cinerea, which may affect plastome evolution. Furthermore, four highly variable regions were found, including the rpl22 coding region, psbM/trnD-GUC non-coding region, petA/psbJ non-coding region, and ycf1 coding region. These high variable regions can be used as candidate molecular markers and as a reference for identifying future Salix species. In addition, phylogenetic analysis indicated that the cp genome of S. floderusii is sister to Salix cupularis and belongs to the Subgenus Vetrix. Genes (Sf-trnI, Sf-PpsbA, aadA, Sf-TpsbA, Sf-trnA) obtained via cloning were inserted into the pBluescript II SK (+) to yield the cp expression vectors, which harbored the selectable marker gene aadA. The results of a spectinomycin resistance test indicated that the cp expression vector had been successfully constructed. Moreover, the aadA gene was efficiently expressed under the regulation of predicted regulatory elements. The present study provides a solid foundation for establishing subsequent S. floderusii cp transformation systems and developing strategies for the genetic improvement of S. floderusii.

Complete Chloroplast Genome Features of Dendrocalamusfarinosus and Its Comparison and Evolutionary Analysis with Other Bambusoideae Species

Dendrocalamus farinosus is one of the essential bamboo species mainly used for food and timber in the southwestern region of China. In this study, the complete chloroplast (cp) genome of D. farinosus is sequenced, assembled, and the phylogenetic relationship analyzed. The cp genome has a circular and quadripartite structure, has a total length of 139,499 bp and contains 132 genes: 89 protein-coding genes, eight rRNAs and 35 tRNAs. The repeat analyses showed that three types of repeats (palindromic, forward and reverse) are present in the genome. A total of 51 simple sequence repeats are identified in the cp genome. The comparative analysis between different species belonging to Dendrocalamus revealed that although the cp genomes are conserved, many differences exist between the genomes. The analysis shows that the non-coding regions were more divergent than the coding regions, and the inverted repeat regions are more conserved than the single-copy regions. Moreover, these results also indicate that rpoC2 may be used to distinguish between different bamboo species. Phylogenetic analysis results supported that D. farinosus was closely related to D. latiflorus. Furthermore, these bamboo species’ geographical distribution and rhizome types indicate two evolutionary pathways: one is from the tropics to the alpine zone, and the other is from the tropics to the warm temperate zone. Our study will be helpful in the determination of the cp genome sequences of D. farinosus, and provides new molecular data to understand the Bambusoideae evolution.

The chloroplast genome of Farsetia hamiltonii Royle, phylogenetic analysis, and comparative study with other members of Clade C of Brassicaceae

BACKGROUND

Farsetia hamiltonii Royle is a medicinally important annual plant from the Cholistan desert that belongs to the tribe Anastaticeae and clade C of the Brassicaceae family. We provide the entire chloroplast sequence of F.hamiltonii, obtained using the Illumina HiSeq2500 and paired-end sequencing. We compared F. hamiltonii to nine other clade C species, including Farsetia occidentalis, Lobularia libyca, Notoceras bicorne, Parolinia ornata, Morettia canescens, Cochlearia borzaeana, Megacarpaea polyandra, Biscutella laevigata, and Iberis amara. We conducted phylogenetic research on the 22 Brassicaceae species, which included members from 17 tribes and six clades.

RESULTS

The chloroplast genome sequence of F.hamiltonii of 154,802 bp sizes with 36.30% GC content and have a typical structure comprised of a Large Single Copy (LSC) of 83,906 bp, a Small Single Copy (SSC) of 17,988 bp, and two copies of Inverted Repeats (IRs) of 26,454 bp. The genomes of F. hamiltonii and F. occidentalis show shared amino acid frequencies and codon use, RNA editing sites, simple sequence repeats, and oligonucleotide repeats. The maximum likelihood tree revealed Farsetia as a monophyletic genus, closely linked to Morettia, with a bootstrap score of 100. The rate of transversion substitutions (Tv) was higher than the rate of transition substitutions (Ts), resulting in Ts/Tv less than one in all comparisons with F. hamiltonii, indicating that the species are closely related. The rate of synonymous substitutions (Ks) was greater than non-synonymous substitutions (Ka) in all comparisons with F. hamiltonii, with a Ka/Ks ratio smaller than one, indicating that genes underwent purifying selection. Low nucleotide diversity values range from 0.00085 to 0.08516, and IR regions comprise comparable genes on junctions with minimal change, supporting the conserved status of the selected chloroplast genomes of the clade C of the Brassicaceae family. We identified ten polymorphic regions, including rps8-rpl14, rps15-ycf1, ndhG-ndhI, psbK-psbI, ccsA-ndhD, rpl36-rps8, petA-psbJ, ndhF-rpl32, psaJ-rpl3, and ycf1 that might be exploited to construct genuine and inexpensive to solve taxonomic discrepancy and understand phylogenetic relationship amongst Brassicaceae species.

CONCLUSION

The entire chloroplast sequencing of F. hamiltonii sheds light on the divergence of genic chloroplast sequences among members of the clade C. When other Farsetia species are sequenced in the future, the full F. hamiltonii chloroplast will be used as a source for comprehensive taxonomical investigations of the genus. The comparison of F. hamiltonii and other clade C species adds new information to the phylogenetic data and evolutionary processes of the clade. The results of this study will also provide further molecular uses of clade C chloroplasts for possible plant genetic modifications and will help recognise more Brassicaceae family species.

Complete chloroplast genome of Ilex dabieshanensis: Genome structure, comparative analyses with three traditional Ilex tea species, and its phylogenetic relationships within the family Aquifoliaceae

Ilex dabieshanensis K. Yao & M. B. Deng is not only a highly valued tree species for landscaping, it is also a good material for making kuding tea due to its anti-inflammatory and lipid-lowering medicinal properties. Utilizing next-generation and long-read sequencing technologies, we assembled the whole chloroplast genome of I. dabieshanensis. The genome was 157,218 bp in length, exhibiting a typical quadripartite structure with a large single copy (LSC: 86,607 bp), a small single copy (SSC: 18,427 bp) and a pair of inverted repeat regions (IRA and IRB: each of 26,092 bp). A total of 121 predicted genes were encoded, including 113 distinctive (79 protein-coding genes, 30 tRNAs, and 4 rRNAs) and 8 duplicated (8 protein-coding genes) located in the IR regions. Overall, 132 SSRs and 43 long repeats were detected and could be used as potential molecular markers. Comparative analyses of four traditional Ilex tea species (I. dabieshanensis, I. paraguariensis, I. latifolia and I. cornuta) revealed seven divergent regions: matK-rps16, trnS-psbZ, trnT-trnL, atpB-rbcL, petB-petD, rpl14-rpl16, and rpl32-trnL. These variations might be applicable for distinguishing different species within the genus Ilex. Phylogenetic reconstruction strongly suggested that I. dabieshanensis formed a sister clade to I. cornuta and also showed a close relationship to I. latifolia. The generated chloroplast genome information in our study is significant for Ilex tea germplasm identification, phylogeny and genetic improvement.

The complete chloroplast genome sequence of the medicinal plant Abrus pulchellus subsp. cantoniensis: genome structure, comparative and phylogenetic relationship analysis

Abrus pulchellus subsp. cantoniensis, an endemic medicinal plant in southern China, is clinically used to treat jaundice hepatitis, cholecystitis, stomachache and breast carbuncle. Here, we assembled and analyzed the first complete chloroplast (cp) genome of A. pulchellus subsp. cantoniensis. The A. pulchellus subsp. cantoniensis cp genome size is 156,497 bp with 36.5% GC content. The cp genome encodes 130 genes, including 77 protein-coding genes, 30 tRNA genes and four rRNA genes, of which 19 genes are duplicated in the inverted repeats (IR) regions. A total of 30 codons exhibited codon usage bias with A/U-ending. Moreover, 53 putative RNA editing sites were predicted in 20 genes, all of which were cytidine to thymine transitions. Repeat sequence analysis identified 45 repeat structures and 125 simple-sequence repeats (SSRs) in A. pulchellus subsp. cantoniensis cp genome. In addition, 19 mononucleotides (located in atpB, trnV-UAC, ycf3, atpF, rps16, rps18, clpP, rpl16, trnG-UCC and ndhA) and three compound SSRs (located in ndhA, atpB and rpl16) showed species specificity between A. pulchellus subsp. cantoniensis and Abrus precatorius, which might be informative sources for developing molecular markers for species identification. Furthermore, phylogenetic analysis inferred that A. pulchellus subsp. cantoniensis was closely related to A. precatorius, and the genus Abrus formed a subclade with Canavalia in the Millettioid/Phaseoloid clade. These data provide a valuable resource to facilitate the evolutionary relationship and species identification of this species.

Complete chloroplast genome features of the model heavy metal hyperaccumulator Arabis paniculata Franch and its phylogenetic relationships with other Brassicaceae species

Arabis paniculata Franch (Brassicaceae) has been widely used for the phytoremediation of heavy mental, owing to its hyper tolerance of extreme Pb, Zn, and Cd concentrations. However, studies on its genome or plastid genome are scarce. In the present study, we obtained the complete chloroplast (cp) genome of A. paniculata via de novo assembly through the integration of Illumina reads and PacBio subreads. The cp genome presents a typical quadripartite cycle with a length of 153,541 bp, and contains 111 unigenes, with 79 protein-coding genes, 28 tRNAs and 4 rRNAs. Codon usage analysis showed that the codons for leucine were the most frequent codons and preferentially ended with A/U. Synonymous (Ks) and non-synonymous (Ka) substitution rate analysis indicated that the unigenes, ndhF and rpoC2, related to "NADH-dehydrogenase" and "RNA polymerase" respectively, underwent the lowest purifying selection pressure. Phylogenetic analysis demonstrated that Arabis flagellosa and A. hirsuta are more similar to each other than to A. paniculata, and Arabis is the closest relative of Draba among all Brassicaceae genera. These findings provide valuable information for the optimal exploitation of this model species as a heavy-metal hyperaccumulator.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01151-1.

Characterization of the chloroplast genome of Lonicera ruprechtiana Regel and comparison with other selected species of Caprifoliaceae

Lonicera ruprechtiana Regel is widely used as a greening tree in China and also displays excellent pharmacological activities. The phylogenetic relationship between L. ruprechtiana and other members of Caprifoliaceae remains unclear. In this study, the complete cp genome of L. ruprechtiana was identified using high-throughput Illumina pair-end sequencing data. The circular cp genome was 154,611 bp long and has a large single-copy region of 88,182 bp and a small single-copy region of 18,713 bp, with the two parts separated by two inverted repeat (IR) regions (23,858 bp each). A total of 131 genes were annotated, including 8 ribosomal RNAs, 39 transfer RNAs, and 84 protein-coding genes (PCGs). In addition, 49 repeat sequences and 55 simple sequence repeat loci of 18 types were also detected. Codon usage analysis demonstrated that the Leu codon is preferential for the A/U ending. Maximum-likelihood phylogenetic analysis using 22 Caprifoliaceae species revealed that L. ruprechtiana was closely related to Lonicera insularis. Comparison of IR regions revealed that the cp genome of L. ruprechtiana was largely conserved with that of congeneric species. Moreover, synonymous (Ks) and non-synonymous (Ka) substitution rate analysis showed that most genes were under purifying selection pressure; ycf3, and some genes associated with subunits of NADH dehydrogenase, subunits of the cytochrome b/f complex, and subunits of the photosystem had been subjected to strong purifying selection pressure (Ka/Ks < 0.1). This study provides useful genetic information for future study of L. ruprechtiana evolution.

Complete chloroplast genome of Hordeum brevisubulatum: Genome organization, synonymous codon usage, phylogenetic relationships, and comparative structure analysis

BACKGROUND

Hordeum brevisubulatum, known as fine perennial forage, is used for soil salinity improvement in northern China. Chloroplast (cp) genome is an ideal model for assessing its genome evolution and the phylogenetic relationships. We de novo sequenced and analyzed the cp genome of H. brevisubulatum, providing a fundamental reference for further studies in genetics and molecular breeding.

RESULTS

The cp genome of H. brevisubulatum was 137,155 bp in length with a typical quadripartite structure. A total of 130 functional genes were annotated and the gene of accD was lost in the process of evolution. Among all the annotated genes, 16 different genes harbored introns and the genes of ycf3 and rps12 contained two introns. Parity rule 2 (PR2) plot analysis showed that majority of genes had a bias toward T over A in the coding strand in all five Hordeum species, and a slight G over C in the other four Hordeum species except for H. bogdanil. Additionally, 52 dispersed repeat sequences and 182 simple sequence repeats were identified. Moreover, some unique SSRs of each species could be used as molecular markers for further study. Compared to the other four Hordeum species, H. brevisubulatum was most closely related to H. bogdanii and its cp genome was relatively conserved. Moreover, inverted repeat regions (IRa and IRb) were less divergent than other parts and coding regions were relatively conserved compared to non-coding regions. Main divergence was presented at the SSC/IR border.

CONCLUSIONS

This research comprehensively describes the architecture of the H. brevisubulatum cp genome and improves our understanding of its cp biology and genetic diversity, which will facilitate biological discoveries and cp genome engineering.

Comparative and phylogenetic analyses of eleven complete chloroplast genomes of Dipterocarpoideae

BACKGROUND

In South-east Asia, Dipterocarpoideae is predominant in most mature forest communities, comprising around 20% of all trees. As large quantity and high quality wood are produced in many species, Dipterocarpoideae plants are the most important and valuable source in the timber market. The d-borneol is one of the essential oil components from Dipterocarpoideae (for example, Dryobalanops aromatica or Dipterocarpus turbinatus) and it is also an important traditional Chinese medicine (TCM) formulation known as "Bingpian" in Chinese, with antibacterial, analgesic and anti-inflammatory effects and can enhance anticancer efficiency.

METHODS

In this study, we analyzed 20 chloroplast (cp) genomes characteristics of Dipterocarpoideae, including eleven newly reported genomes and nine cp genomes previously published elsewhere, then we explored the chloroplast genomic features, inverted repeats contraction and expansion, codon usage, amino acid frequency, the repeat sequences and selective pressure analyses. At last, we constructed phylogenetic relationships of Dipterocarpoideae and found the potential barcoding loci.

RESULTS

The cp genome of this subfamily has a typical quadripartite structure and maintains a high degree of consistency among species. There were slightly more tandem repeats in cp genomes of Dipterocarpus and Vatica, and the psbH gene was subjected to positive selection in the common ancestor of all the 20 species of Dipterocarpoideae compared with three outgroups. Phylogenetic tree showed that genus Shorea was not a monophyletic group, some Shorea species and genus Parashorea are placed in one clade. In addition, the rpoC2 gene can be used as a potential marker to achieve accurate and rapid species identification in subfamily Dipterocarpoideae.

CONCLUSIONS

Dipterocarpoideae had similar cp genomic features and psbM, rbcL, psbH may function in the growth of Dipterocarpoideae. Phylogenetic analysis suggested new taxon treatment is needed for this subfamily indentification. In addition, rpoC2 is potential to be a barcoding gene to TCM distinguish.

Gene Loss, Pseudogenization in Plastomes of Genus Allium (Amaryllidaceae), and Putative Selection for Adaptation to Environmental Conditions

Amaryllidaceae is a large family with more than 1,600 species, belonging to 75 genera. The largest genus—Allium—is vast, comprising about a thousand species. Allium species (as well as other members of the Amaryllidaceae) are widespread and diversified, they are adapted to a wide range of habitats from shady forests to open habitats like meadows, steppes, and deserts. The genes present in chloroplast genomes (plastomes) play fundamental roles for the photosynthetic plants. Plastome traits could thus be associated with geophysical abiotic characteristics of habitats. Most chloroplast genes are highly conserved and are used as phylogenetic markers for many families of vascular plants. Nevertheless, some studies revealed signatures of positive selection in chloroplast genes of many plant families including Amaryllidaceae. We have sequenced plastomes of the following nine Allium (tribe Allieae of Allioideae) species: A. zebdanense, A. moly, A. victorialis, A. macleanii, A. nutans, A. obliquum, A. schoenoprasum, A. pskemense, A. platyspathum, A. fistulosum, A. semenovii, and Nothoscordum bivalve (tribe Leucocoryneae of Allioideae). We compared our data with previously published plastomes and provided our interpretation of Allium plastome genes’ annotations because we found some noteworthy inconsistencies with annotations previously reported. For Allium species we estimated the integral evolutionary rate, counted SNPs and indels per nucleotide position as well as compared pseudogenization events in species of three main phylogenetic lines of genus Allium to estimate whether they are potentially important for plant physiology or just follow the phylogenetic pattern. During examination of the 38 species of Allium and the 11 of other Amaryllidaceae species we found that rps16, rps2, infA, ccsA genes have lost their functionality multiple times in different species (regularly evolutionary events), while the pseudogenization of other genes was stochastic events. We found that the “normal” or “pseudo” state of rps16, rps2, infA, ccsA genes correlates well with the evolutionary line of genus the species belongs to. The positive selection in various NADH dehydrogenase (ndh) genes as well as in matK, accD, and some others were found. Taking into account known mechanisms of coping with excessive light by cyclic electron transport, we can hypothesize that adaptive evolution in genes, coding subunits of NADH-plastoquinone oxidoreductase could be driven by abiotic factors of alpine habitats, especially by intensive light and UV radiation.

Chloroplast genome features of Moricandia arvensis (Brassicaceae), a C3-C4 intermediate photosynthetic species

Moricandia arvensis, a plant species originating from the Mediterranean, has been classified as a rare C3-C4 intermediate species, and it is a possible bridge during the evolutionary process from C3 to C4 plant photosynthesis in the family Brassicaceae. Understanding the genomic structure, gene order, and gene content of chloroplasts (cp) of such species can provide a glimpse into the evolution of photosynthesis. In the present study, we obtained a well-annotated cp genome of M. arvensis using long PacBio and short Illumina reads with a de novo assembly strategy. The M. arvensis cp genome was a quadripartite circular molecule with the length of 153,312 bp, including two inverted repeats (IR) regions of 26,196 bp, divided by a small single copy (SSC) region of 17,786 bp and a large single copy (LSC) region of 83,134 bp. We detected 112 unigenes in this genome, comprising 79 protein-coding genes, 29 tRNAs, and four rRNAs. Forty-nine long repeat sequences and 51 simple sequence repeat (SSR) loci of 15 repeat types were identified. The analysis of Ks (synonymous) and Ka (non-synonymous) substitution rates indicated that the genes associated with “subunits of ATP synthase” (atpB), “subunits of NADH-dehydrogenase” (ndhG and ndhE), and “self-replication” (rps12 and rpl16) showed relatively higher Ka/Ks values than those of the other genes. The gene content, gene order, and LSC/IR/SSC boundaries and adjacent genes of the M. arvensis cp genome were highly conserved compared to those in related C3 species. Our phylogenetic analysis demonstrated that M. arvensis was clustered into a subclade with cultivated Brassica species and Raphanus sativus, indicating that M. arvensis was not involved in an independent evolutionary origin event. These results will open the way for further studies on the evolutionary process from C3 to C4 photosynthesis and hopefully provide guidance for utilizing M. arvensis as a resource for improvinng photosynthesis efficiency in cultivated Brassica species.

Characterization of the Complete Chloroplast Genome of Buddleja Lindleyana

BACKGROUND

Buddleja lindleyana Fort., which belongs to the Loganiaceae with a distribution throughout the tropics, is widely used as an ornamental plant in China. Buddleja contains several morphologically similar species, which need to be identified by molecular identification. But there is little molecular research on the genus Buddleja.

OBJECTIVE

Using molecular biology techniques to sequence and analyze the complete chloroplast (cp) genome of B. lindleyana.

METHODS

According to next-generation sequencing to sequence the genome data, a series of bioinformatics software were used to assembly and analysis the molecular structure of cp genome of B. lindleyana.

RESULTS

The complete cp genome of B. lindleyana is a circular 154,487-bp-long molecule with a GC content of 38.1%. It has a familiar quadripartite structure, including a large single-copy region (LSC; 85,489 bp), a small single-copy region (SSC; 17,898bp) and a pair of inverted repeats (IRs; 25,550 bp). A total of 133 genes were identified in the genome, including 86 protein-coding genes, 37 tRNA genes, 8 rRNA genes and 2 pseudogenes.

CONCLUSIONS

These results suggested that B. lindelyana cp genome could be used as a potential genomic resource to resolve the phylogenetic positions and relationships of Loganiaceae, and will offer valuable information for future research in the identification of Buddleja species and will conduce to genomic investigations of these species.

Chloroplast genomes of two Mediterranean Bupleurum species and the phylogenetic relationship inferred from combined analysis with East Asian species

The chloroplast genomes of Mediterranean Bupleurum species are reported for the first time. Phylogenetic analysis supports the species as a basal clade of Bupleurum with divergence time at 35.40 Ma. Bupleurum is one of the most species-rich genus with high medicinal value in Apiaceae. Although infrageneric classifications of Bupleurum have been the subject of numerous studies, it still remains controversial. Chloroplast genome information will prove essential in advancing our understanding on phylogenetic study. Here we report cp genomes of two woody Bupleurum species (Bupleurum gibraltaricum and B. fruticosum) endemic to Mediterranean. The complete cp genomes of the two species were 157,303 and 157,391 bp in size, respectively. They encoded 114 unique genes including 30 tRNA genes, 4 rRNA genes and 80 protein coding genes. Genome structure, distributions of SDRs and SSRs, gene content exhibited similarities among Bupleurum species. High variable hotspots were detected in eight intergenic spacers and four genes. Most of genes were under purifying selection with two exceptions: atpF and clpP. The phylogenetic analysis based on 80 coding genes revealed that the genus was divided into 2 distinct clades corresponding to the 2 subgenera (subg. Penninervia, subg. Bupleurum) with divergence time at the end of collision of India with Eurasia. Most species diversified mainly during the later period of uplift of Qinghai-Tibetan Plateau. The cp genomes of the two Bupleurum species can be significant complementary to insights into the cp genome characteristics of this genus. The comparative chloroplast genomes and phylogenetic analysis advances our understanding of the evolution of cp genomes and phylogeny in Bupleurum.

Chloroplast genome sequence of Chongming lima bean (Phaseolus lunatus L.) and comparative analyses with other legume chloroplast genomes

Background

Lima bean (Phaseolus lunatus L.) is a member of subfamily Phaseolinae belonging to the family Leguminosae and an important source of plant proteins for the human diet. As we all know, lima beans have important economic value and great diversity. However, our knowledge of the chloroplast genome level of lima beans is limited.

Results

The chloroplast genome of lima bean was obtained by Illumina sequencing technology for the first time. The Cp genome with a length of 150,902 bp, including a pair of inverted repeats (IRA and IRB 26543 bp each), a large single-copy (LSC 80218 bp) and a small single-copy region (SSC 17598 bp). In total, 124 unique genes including 82 protein-coding genes, 34 tRNA genes, and 8 rRNA genes were identified in the P. lunatus Cp genome. A total of 61 long repeats and 290 SSRs were detected in the lima bean Cp genome. It has a typical 50 kb inversion of the Leguminosae family and an 70 kb inversion to subtribe Phaseolinae. rpl16, accD, petB, rsp16, clpP, ndhA, ndhF and ycf1 genes in coding regions was found significant variation, the intergenic regions of trnk-rbcL, rbcL-atpB, ndhJ-rps4, psbD-rpoB, atpI-atpA, atpA-accD, accD-psbJ, psbE-psbB, rsp11-rsp19, ndhF-ccsA was found in a high degree of divergence. A phylogenetic analysis showed that P. lunatus appears to be more closely related to P. vulgaris, V.unguiculata and V. radiata.

Conclusions

The characteristics of the lima bean Cp genome was identified for the first time, these results will provide useful insights for species identification, evolutionary studies and molecular biology research.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07467-8.

Complete chloroplast genome features and phylogenetic analysis of Eruca sativa (Brassicaceae)

Eruca sativa Mill. (Brassicaceae) is an important edible vegetable and a potential medicinal plant due to the antibacterial activity of its seed oil. Here, the complete chloroplast (cp) genome of E. sativa was de novo assembled with a combination of long PacBio reads and short Illumina reads. The E. sativa cp genome had a quadripartite structure that was 153,522 bp in size, consisting of one large single-copy region of 83,320 bp and one small single-copy region of 17,786 bp which were separated by two inverted repeat (IRa and IRb) regions of 26,208 bp. This complete cp genome harbored 113 unique genes: 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. Forty-nine long repetitive sequences and 69 simple sequence repeats were identified in the E. sativa cp genome. A codon usage analysis of the E. sativa cp genome showed a bias toward codons ending in A/T. The E. sativa cp genome was similar in size, gene composition, and linearity of the structural region when compared with other Brassicaceae cp genomes. Moreover, the analysis of the synonymous (Ks) and non-synonymous (Ka) substitution rates demonstrated that protein-coding genes generally underwent purifying selection pressure, expect ycf1, ycf2, and rps12. A phylogenetic analysis determined that E. sativa is evolutionarily close to important Brassica species, indicating that it may be possible to transfer favorable E. sativa alleles into other Brassica species. Our results will be helpful to advance genetic improvement and breeding of E. sativa, and will provide valuable information for utilizing E. sativa as an important resource to improve other Brassica species.

Chloroplast Genome of Rambutan and Comparative Analyses in Sapindaceae

Rambutan (Nephelium lappaceum L.) is an important fruit tree that belongs to the family Sapindaceae and is widely cultivated in Southeast Asia. We sequenced its chloroplast genome for the first time and assembled 161,321 bp circular DNA. It is characterized by a typical quadripartite structure composed of a large (86,068 bp) and small (18,153 bp) single-copy region interspersed by two identical inverted repeats (IRs) (28,550 bp). We identified 132 genes including 78 protein-coding genes, 29 tRNA and 4 rRNA genes, with 21 genes duplicated in the IRs. Sixty-three simple sequence repeats (SSRs) and 98 repetitive sequences were detected. Twenty-nine codons showed biased usage and 49 potential RNA editing sites were predicted across 18 protein-coding genes in the rambutan chloroplast genome. In addition, coding gene sequence divergence analysis suggested that ccsA, clpP, rpoA, rps12, psbJ and rps19 were under positive selection, which might reflect specific adaptations of N. lappaceum to its particular living environment. Comparative chloroplast genome analyses from nine species in Sapindaceae revealed that a higher similarity was conserved in the IR regions than in the large single-copy (LSC) and small single-copy (SSC) regions. The phylogenetic analysis showed that N. lappaceum chloroplast genome has the closest relationship with that of Pometia tomentosa. The understanding of the chloroplast genomics of rambutan and comparative analysis of Sapindaceae species would provide insight into future research on the breeding of rambutan and Sapindaceae evolutionary studies.

Chloroplast Genome Evolution in Four Montane Zingiberaceae Taxa in China

Chloroplasts are critical to plant survival and adaptive evolution. The comparison of chloroplast genomes could provide insight into the adaptive evolution of closely related species. To identify potential adaptive evolution in the chloroplast genomes of four montane Zingiberaceae taxa (Cautleya, Roscoea, Rhynchanthus, and Pommereschea) that inhabit distinct habitats in the mountains of Yunnan, China, the nucleotide sequences of 13 complete chloroplast genomes, including five newly sequenced species, were characterized and compared. The five newly sequenced chloroplast genomes (162,878-163,831 bp) possessed typical quadripartite structures, which included a large single copy (LSC) region, a small single copy (SSC) region, and a pair of inverted repeat regions (IRa and IRb), and even though the structure was highly conserved among the 13 taxa, one of the rps19 genes was absent in Cautleya, possibly due to expansion of the LSC region. Positive selection of rpoA and ycf2 suggests that these montane species have experienced adaptive evolution to habitats with different sunlight intensities and that adaptation related to the chloroplast genome has played an important role in the evolution of Zingiberaceae taxa.

Chloroplast Genomes of Two Species of Cypripedium: Expanded Genome Size and Proliferation of AT-Biased Repeat Sequences

The size of the chloroplast genome (plastome) of autotrophic angiosperms is generally conserved. However, the chloroplast genomes of some lineages are greatly expanded, which may render assembling these genomes from short read sequencing data more challenging. Here, we present the sequencing, assembly, and annotation of the chloroplast genomes of Cypripedium tibeticum and Cypripedium subtropicum. We de novo assembled the chloroplast genomes of the two species with a combination of short-read Illumina data and long-read PacBio data. The plastomes of the two species are characterized by expanded genome size, proliferated AT-rich repeat sequences, low GC content and gene density, as well as low substitution rates of the coding genes. The plastomes of C. tibeticum (197,815 bp) and C. subtropicum (212,668 bp) are substantially larger than those of the three species sequenced in previous studies. The plastome of C. subtropicum is the longest one of Orchidaceae to date. Despite the increase in genome size, the gene order and gene number of the plastomes are conserved, with the exception of an ∼75 kb large inversion in the large single copy (LSC) region shared by the two species. The most striking is the record-setting low GC content in C. subtropicum (28.2%). Moreover, the plastome expansion of the two species is strongly correlated with the proliferation of AT-biased non-coding regions: the non-coding content of C. subtropicum is in excess of 57%. The genus provides a typical example of plastome expansion induced by the expansion of non-coding regions. Considering the pros and cons of different sequencing technologies, we recommend hybrid assembly based on long and short reads applied to the sequencing of plastomes with AT-biased base composition.

Sequencing and comparative analysis of the chloroplast genome of Ribes odoratum provide insights for marker development and phylogenetics in Ribes

Ribes odoratum, commonly known as clove currant, is a perennial deciduous shrub noted for its long-lasting fragrant flowers and edible fruits. Owing to its ornamental values, this species has been widely used in city gardening and urban landscaping. Here, the complete cp genome of R. odoratum was de novo assembled for the first time. The plastome is 157,152 bp in length, with a GC content of 38.2%. The cp genome featured a typical quadripartite structure, consisting of a pair of inverted repeat regions of 25,961 bp, separated by a large single copy region of 86,896 bp, and a small single copy region of 18,333 bp. A total of 131 genes were annotated in the plastome, including 86 protein coding genes, 37 tRNA genes, and 8 rRNA genes. 56 SSRs were identified, among which, 82.35% were located in the intergenic regions. A strong A/T bias in base composition was observed in these cpSSRs. In addition, 49 repeats of different sizes and types were also found in the plastome. Through comparison, seven divergence hotspots were identified between the cp genomes of R. odoratum and R. fasciculatum var. chinense. Sequences of these divergent regions could be developed as potential markers for species delimitation in further studies. We re-investigated the relationship aomong 32 Saxifragales species through plastome-based phylogenywhich revealed that R. odoratum as a sister of R. fasciculatum var. chinense. Thus, our study provides genomic resources and valuable reference for marker development and phylogenomics in Ribes.

SUPPLEMENTARY INFORMATION

The online version of this article (10.1007/s12298-021-00932-4) contains supplementary material, which is available to authorized users.

Complete chloroplast genomes of three important species, Abelmoschus moschatus, A. manihot and A. sagittifolius: Genome structures, mutational hotspots, comparative and phylogenetic analysis in Malvaceae

Abelmoschus is an economically and phylogenetically valuable genus in the family Malvaceae. Owing to coexistence of wild and cultivated form and interspecific hybridization, this genus is controversial in systematics and taxonomy and requires detailed investigation. Here, we present whole chloroplast genome sequences and annotation of three important species: A. moschatus, A. manihot and A. sagittifolius, and compared with A. esculentus published previously. These chloroplast genome sequences ranged from 163121 bp to 163453 bp in length and contained 132 genes with 87 protein-coding genes, 37 transfer RNA and 8 ribosomal RNA genes. Comparative analyses revealed that amino acid frequency and codon usage had similarity among four species, while the number of repeat sequences in A. esculentus were much lower than other three species. Six categories of simple sequence repeats (SSRs) were detected, but A. moschatus and A. manihot did not contain hexanucleotide SSRs. Single nucleotide polymorphisms (SNPs) of A/T, T/A and C/T were the largest number type, and the ratio of transition to transversion was from 0.37 to 0.55. Abelmoschus species showed relatively independent inverted-repeats (IR) boundary traits with different boundary genes compared with the other related Malvaceae species. The intergenic spacer regions had more polymorphic than protein-coding regions and intronic regions, and thirty mutational hotpots (≥200 bp) were identified in Abelmoschus, such as start-psbA, atpB-rbcL, petD-exon2-rpoA, clpP-intron1 and clpP-exon2.These mutational hotpots could be used as polymorphic markers to resolve taxonomic discrepancies and biogeographical origin in genus Abelmoschus. Moreover, phylogenetic analysis of 33 Malvaceae species indicated that they were well divided into six subfamilies, and genus Abelmoschus was a well-supported clade within genus Hibiscus.

The Complete Plastid Genome of Artocarpus camansi: A High Degree of Conservation of the Plastome Structure in the Family Moraceae

Understanding the plastid genome is extremely important for the interpretation of the genetic mechanisms associated with essential physiological and metabolic functions, the identification of possible marker regions for phylogenetic or phylogeographic analyses, and the elucidation of the modes through which natural selection operates in different regions of this genome. In the present study, we assembled the plastid genome of Artocarpus camansi, compared its repetitive structures with Artocarpus heterophyllus, and searched for evidence of synteny within the family Moraceae. We also constructed a phylogeny based on 56 chloroplast genes to assess the relationships among three families of the order Rosales, that is, the Moraceae, Rhamnaceae, and Cannabaceae. The plastid genome of A. camansi has 160,096 bp, and presents the typical circular quadripartite structure of the Angiosperms, comprising a large single copy (LSC) of 88,745 bp and a small single copy (SSC) of 19,883 bp, separated by a pair of inverted repeat (IR) regions each with a length of 25,734 bp. The total GC content was 36.0%, which is very similar to Artocarpus heterophyllus (36.1%) and other moraceous species. A total of 23,068 codons and 80 SSRs were identified in the A. camansi plastid genome, with the majority of the SSRs being mononucleotide (70.0%). A total of 50 repeat structures were observed in the A. camansi plastid genome, in contrast with 61 repeats in A. heterophyllus. A purifying selection signal was found in 70 of the 79 protein-coding genes, indicating that they have all been highly conserved throughout the evolutionary history of the genus. The comparative analysis of the structural characteristics of the chloroplast among different moraceous species found a high degree of similarity in the sequences, which indicates a highly conserved evolutionary model in these plastid genomes. The phylogenetic analysis also recovered a high degree of similarity between the chloroplast genes of A. camansi and A. heterophyllus, and reconfirmed the hypothesis of the intense conservation of the plastome in the family Moraceae.

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 Chloroplast Genome of Carya illinoinensis: Genome Structure, Adaptive Evolution, and Phylogenetic Analysis

Research Highlights: For the first time, the complete chloroplast (cp) genome of Carya illinoinensis cv. ‘Pawnee’ was de novo assembled. Comprehensive analysis the cp genome of C. illinoinensis revealed potential cpDNA markers for intraspecies identification, genes involved in adaptation, and its phylogenetic position. Background and Objectives: C. illinoinensis is an economically important nut tree in the family Juglandaceae. Cp-derived markers are helpful for genetic research, but they still need to be developed in C. illinoinensis. Additionally, the adaptation and phylogenetic relationships of C. illinoinensis have not been revealed based on the complete cp genome. Materials and Methods: Chloroplast genomic DNA of C. illinoinensis cv. ‘Pawnee’ was extracted and subjected to Illumina sequencing. Results: The cp genome is 160,819 bp in size, exhibiting a typical quadripartite structure with a large single copy (LSC) of 90,022 bp, a small single copy (SSC) of 18,791 bp, and a pair of inverted repeats (IRA and IRB) regions of 26,003 bp each. The genome was predicted to encode 112 unique genes, including 79 protein-coding genes, 29 tRNAs, and four rRNAs, with 19 duplicates in the IR regions. In total, 213 SSRs and 44 long repeats were identified in the cp genome. A comparison of two different C. illinoinensis genotypes, ‘Pawnee’ and 87MX3-2.11, obtained 143 SNPs and 74 indels. The highly variable regions such as atpF, clpP, and ndhA genes, and matK-rps16, trnS-trnG, and trnT-psbD intergenic spacers might be helpful for future intraspecific identification. Positive selection was acting on the ccsA and rps12 cp genes based on the Ka/Ks ratios. Phylogenetic analysis indicated that C. illinoinensis forms a sister clade to Asian Carya species, represented by C. kweichowensis and Annamocarya sinensis. Conclusions: The genome information in our study will have significance for further research on the intraspecies identification and genetic improvement of C. illinoinensis.

Comparative analysis of chloroplast genomes for five Dicliptera species (Acanthaceae): molecular structure, phylogenetic relationships, and adaptive evolution

The genus Dicliptera (Justicieae, Acanthaceae) consists of approximately 150 species distributed throughout the tropical and subtropical regions of the world. Newly obtained chloroplast genomes (cp genomes) are reported for five species of Dilciptera (D. acuminata, D. peruviana, D. montana, D. ruiziana and D. mucronata) in this study. These cp genomes have circular structures of 150,689-150,811 bp and exhibit quadripartite organizations made up of a large single copy region (LSC, 82,796-82,919 bp), a small single copy region (SSC, 17,084-17,092 bp), and a pair of inverted repeat regions (IRs, 25,401-25,408 bp). Guanine-Cytosine (GC) content makes up 37.9%-38.0% of the total content. The complete cp genomes contain 114 unique genes, including 80 protein-coding genes, 30 transfer RNA (tRNA) genes, and four ribosomal RNA (rRNA) genes. Comparative analyses of nucleotide variability (Pi) reveal the five most variable regions (trnY-GUA-trnE-UUC, trnG-GCC, psbZ-trnG-GCC, petN-psbM, and rps4-trnL-UUA), which may be used as molecular markers in future taxonomic identification and phylogenetic analyses of Dicliptera. A total of 55-58 simple sequence repeats (SSRs) and 229 long repeats were identified in the cp genomes of the five Dicliptera species. Phylogenetic analysis identified a close relationship between D. ruiziana and D. montana, followed by D. acuminata, D. peruviana, and D. mucronata. Evolutionary analysis of orthologous protein-coding genes within the family Acanthaceae revealed only one gene, ycf15, to be under positive selection, which may contribute to future studies of its adaptive evolution. The completed genomes are useful for future research on species identification, phylogenetic relationships, and the adaptive evolution of the Dicliptera species.

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.

The complete chloroplast genome sequence of Ulmus elongata (Ulmaceae)

Ulmus elongata is a species of Sect. Chaetoptelea (Liemb.) S chneid in Ulmaceae, and it is an endangered wild plant listed in the second class of the Protected Plants in China. The complete chloroplast genome (cp) of U. elongata was reported in this study. The result showed that the cp genome was 159,230 bp in length including a large single-copy (LSC) 87,718 bp and a small single-copy (SSC) 18,690 bp, which were separated by two inverted repeats (IRs) of 26,411 bp with the typical quadripartite structure, respectively. The genome encoded 132 genes, including 87 protein-coding genes, 37 tRNA genes, and eight rRNA genes. The GC content was 35.57%. Chloroplast sequences were used for constructing phylogenetic tree to determine the evolutionary status of U. elongata. The maximum-likelihood phylogenetic analysis showed that U. elongata was clustered with five other Ulmus species, and the relationship between Ulmus and Zelkova was closest. The success of cp genome assembly of U. elongata has laid a foundation for the study of chloroplast molecular biology and can effectively promote the study of genetic breeding and molecular evolution of U. elongata.

Comparative Analysis of the Complete Chloroplast Genomes in Allium Subgenus Cyathophora (Amaryllidaceae): Phylogenetic Relationship and Adaptive Evolution

Recent advances in molecular phylogenetics provide us with information of Allium L. taxonomy and evolution, such as the subgenus Cyathophora, which is monophyletic and contains five species. However, previous studies detected distinct incongruence between the nrDNA and cpDNA phylogenies, and the interspecies relationships of this subgenus need to be furtherly resolved. In our study, we newly assembled the whole chloroplast genome of four species in subgenus Cyathophora and two allied Allium species. The complete cp genomes were found to possess a quadripartite structure, and the genome size ranged from 152,913 to 154,174 bp. Among these cp genomes, there were subtle differences in the gene order, gene content, and GC content. Seven hotspot regions (infA, rps16, rps15, ndhF, trnG-UCC, trnC-GCA, and trnK-UUU) with nucleotide diversity greater than 0.02 were discovered. The selection analysis showed that some genes have elevated Ka/Ks ratios. Phylogenetic analysis depended on the complete chloroplast genome (CCG), and the intergenic spacer regions (IGS) and coding DNA sequences (CDS) showed same topologies with high support, which revealed that subgenus Cyathophora was a monophyletic group, containing four species, and A. cyathophorum var. farreri was sister to A. spicatum with 100% bootstrap value. Our study revealed selective pressure may exert effect on several genes of the six Allium species, which may be useful for them to adapt to their specific living environment. We have well resolved the phylogenetic relationship of species in the subgenus Cyathophora, which will contribute to future evolutionary studies or phylogeographic analysis of Allium.

The complete chloroplast genome sequence of yellow mustard (Sinapis alba L.) and its phylogenetic relationship to other Brassicaceae species

As a member of the large Brassicaceae family, yellow mustard (Sinapis alba L.) has been used as an important gene pool for the genetic improvement of cash crops in Brassicaceae. Understanding the phylogenetic relationship between Sinapis alba (S. alba) and other Brassicaceae crops can provide guidance on the introgression of its favorable alleles into related species. The chloroplast (cp) genome is an ideal model for assessing genome evolution and the phylogenetic relationships of complex angiosperm families. Herein, we de novo assembled the complete cp genome of S. alba by integrating the PacBio and Illumina sequencing platforms. A 153,760 bp quadripartite cycle without any gap was obtained, including a pair of inverted repeats (IRa and IRb) of 26,221 bp, separated by a large single copy (LSC) region of 83,506 bp and a small single copy (SSC) region of 17,821 bp. A total of 78 protein-coding genes, 30 tRNA genes, and four rRNA genes were identified in this cp genome, as were 89 simple sequence repeat (SSR) loci of 18 types. The codon usage analysis revealed a preferential use of the Leu codon with the A/U ending. The phylogenetic analysis using 82 Brassicaceae species demonstrated that S. alba had a close relationship with important Brassica and Raphanus species; moreover, it likely originated from a separate evolutionary pathway compared with the congeneric Sinapis arvensis. The synonymous (Ks) and non-synonymous (Ks) substitution rate analysis showed that genes encoding "Subunits of cytochrome b/f complex" were under the lowest purifying selection pressure, whereas those associated with "Maturase", "Subunit of acetyl-CoA", and "Subunits of NADH-dehydrogenase" underwent relatively higher purifying selection pressures. Our results provide valuable information for fully utilizing the S. alba cp genome as a potential genetic resource for the genetic improvement of Brassica and Raphanus species.

The complete chloroplast genome of Rhododendron delavayi (Ericaceae)

Rhododendron delavayi, as a member of Ericaceae family, has been widely used as an important garden flower. The cp genome of R. delavayi exhibited a typical quadripartite cycle with 193,798 bp, comprising of a pair of inverted repeats (IRa and IRb) of 15,494 bp intersected by a large single copy (LSC) region of 160,234 bp and a quite small single copy region of 2576 bp. Totally, 123 unique genes were assembled in this cp genome, including 80 protein genes, 35 tRNAs and 8 rRNAs. Out of these assembled genes, 88 genes (71.54%) were single copy. Phylogenetic analysis based on 14 cp genome of related species showed that the R. delavayi was closely related to Vaccinium oldhamii. This study provides important information for future evolution, genetic and molecular biology studies of Rhododendron.