Comparative Analysis of the Complete Chloroplast Genome Sequences of Three Closely Related East-Asian Wild Roses (Rosa sect. Synstylae; Rosaceae)

Comparative analysis the chloroplast genomes of Celastrus (Celastraceae) species: Provide insights into molecular evolution, species identification and phylogenetic relationships

BACKGROUND

The genus Celastrus is an important medicinal plant resource. The similarity of morphology and the lack of complete chloroplast genome analysis have significantly impeded the exploration of species identification, molecular evolution and phylogeny of Celastrus.

PURPOSE

In order to resolve the phylogenic controversy of Celastrus species, the chloroplast genome comparative analysis was performed to provide genetic evidence.

METHODS

In this study, we collected and sequenced ten chloroplast genomes of Celastrus species from China and downloaded three chloroplast genomes from the databases. The chloroplast genomes were compared and analyzed to explore their characteristics and evolution. Furthermore, the phylogenetic relationships of Celastrus species were inferred based on the whole chloroplast genomes and protein-coding genes.

RESULTS

All the 13 Celastrus species chloroplast genomes showed a typical quadripartite structure with genome sizes ranging from 155,113 to 157,366 bp. The intron loss of the rps16 gene occurred in all the 13 Celastrus species. The GC content, gene sequence, repeat types and codon bias pattern were highly conserved. Ten highly variation regions were identified, which can be used as potential DNA markers in molecular identification of Celastrus species. Eight genes, including accD, atp4, ndhB, rpoC1, rbcL, rpl2, rpl20 and ycf1, were detected to experience positive selection. Phylogenetic analysis showed that Celastrus was a monophyletic group and Tripterygium was the closest sister-group. Noteworthy, C. gemmatus Loes. and C. orbiculatus Thunb. can be discriminated using the chloroplast genome as a super barcode. The comparative and phylogenetic analysis results proposed that C. tonkinensis Pitard. was the synonym of C. hindsii Benth.

CONCLUSION

The comparative analysis of the Celastrus chloroplast genomes can provide comprehensive genetic evidence for molecular evolution, species identification and phylogenetic relationships.

Full-length transcriptome and RNA-Seq analyses reveal the resistance mechanism of sesame in response to Corynespora cassiicola

BACKGROUND

Corynespora leaf spot is a common leaf disease occurring in sesame, and the disease causes leaf yellowing and even shedding, which affects the growth quality of sesame. At present, the mechanism of sesame resistance to this disease is still unclear. Understanding the resistance mechanism of sesame to Corynespora leaf spot is highly important for the control of infection. In this study, the leaves of the sesame resistant variety (R) and the sesame susceptible variety (S) were collected at 0-48 hpi for transcriptome sequencing, and used a combined third-generation long-read and next-generation short-read technology approach to identify some key genes and main pathways related to resistance.

RESULTS

The gene expression levels of the two sesame varieties were significantly different at 0, 6, 12, 24, 36 and 48 hpi, indicating that the up-regulation of differentially expressed genes in the R might enhanced the resistance. Moreover, combined with the phenotypic observations of sesame leaves inoculated at different time points, we found that 12 hpi was the key time point leading to the resistance difference between the two sesame varieties at the molecular level. The WGCNA identified two modules significantly associated with disease resistance, and screened out 10 key genes that were highly expressed in R but low expressed in S, which belonged to transcription factors (WRKY, AP2/ERF-ERF, and NAC types) and protein kinases (RLK-Pelle_DLSV, RLK-Pelle_SD-2b, and RLK-Pelle_WAK types). These genes could be the key response factors in the response of sesame to infection by Corynespora cassiicola. GO and KEGG enrichment analysis showed that specific modules could be enriched, which manifested as enrichment in biologically important pathways, such as plant signalling hormone transduction, plant-pathogen interaction, carbon metabolism, phenylpropanoid biosynthesis, glutathione metabolism, MAPK and other stress-related pathways.

CONCLUSIONS

This study provides an important resource of genes contributing to disease resistance and will deepen our understanding of the regulation of disease resistance, paving the way for further molecular breeding of sesame.

The complete chloroplast genome sequence of Rosa 'Limoncello' (Rosales: Rosaceae)

Rosa 'Limoncello' finds applications in gardening and landscaping. In this study, we assembled and annotated the complete chloroplast genome of this variety for the first time. The length of its chloroplast genome was 156,493 bp, containing two short inverted repeat regions of 26,052 bp, each separated by a large single-copy region of 85,649 bp and a small single-copy region of 18,740 bp. The chloroplast DNA of R. 'Limoncello' consisted of 135 genes, including 90 protein-coding genes, eight ribosomal RNA genes, and 37 transfer RNA genes. On comparing the complete chloroplast sequence of R. 'Limoncello' with that of other Rosa species, R. 'Limoncello' was found to be closely related to Rosa cymosa. Thus, information on the chloroplast genome sequence of this rose variety can facilitate phylogenetic studies of the genus Rosa.

Comparative analysis of the chloroplast genomes of Rosa species and RNA editing analysis

BACKGROUND

The genus Rosa (Rosaceae) contains approximately 200 species, most of which have high ecological and economic values. Chloroplast genome sequences are important for studying species differentiation, phylogeny, and RNA editing.

RESULTS

In this study, the chloroplast genomes of three Rosa species, Rosa hybrida, Rosa acicularis, and Rosa rubiginosa, were assembled and compared with other reported Rosa chloroplast genomes. To investigate the RNA editing sites in R. hybrida (commercial rose cultivar), we mapped RNA-sequencing data to the chloroplast genome and analyzed their post-transcriptional features. Rosa chloroplast genomes presented a quadripartite structure and had highly conserved gene order and gene content. We identified four mutation hotspots (ycf3-trnS, trnT-trnL, psbE-petL, and ycf1) as candidate molecular markers for differentiation in the Rosa species. Additionally, 22 chloroplast genomic fragments with a total length of 6,192 bp and > 90% sequence similarity with their counterparts were identified in the mitochondrial genome, representing 3.96% of the chloroplast genome. Phylogenetic analysis including all sections and all subgenera revealed that the earliest divergence in the chloroplast phylogeny roughly distinguished species of sections Pimpinellifoliae and Rosa and subgenera Hulthemia. Moreover, DNA- and RNA-sequencing data revealed 19 RNA editing sites, including three synonymous and 16 nonsynonymous, in the chloroplast genome of R. hybrida that were distributed among 13 genes.

CONCLUSIONS

The genome structure and gene content of Rosa chloroplast genomes are similar across various species. Phylogenetic analysis based on the Rosa chloroplast genomes has high resolution. Additionally, a total of 19 RNA editing sites were validated by RNA-Seq mapping in R. hybrida. The results provide valuable information for RNA editing and evolutionary studies of Rosa and a basis for further studies on genomic breeding of Rosa species.

The first chloroplast sequence of Rosa davurica Pall. var. Davurica

Rosa davurica Pall. var. davurica is a member of the plant family Rosaceae. Although R. davurica has high application value, its chloroplast genome sequence has not been reported. This study aims to reveal the genetic characteristics of the chloroplast genome of Rosa roxburghii. The length of its total chloroplast DNA is 156,971 bp, with 37.22% G/C content. Its chloroplast genome has two inverted repeat (IRa and IRb) regions totaling 26,051 bp which are separated by a large single copy (LSC) region of 86,032 bp and a small single copy (SSC) region of 18,837 bp. The genome contains 131 independent genes (86 protein-coding, 37 tRNA, and 8 rRNA), and there are 18 repeated genes within the IR region. Among these genes, 17 genes contained one or two introns. The phylogenetic analysis showed that R. davurica was relatively close to other Rosa species, such as the Rosa hybrid.

Complete chloroplast genomes of five Cuscuta species and their evolutionary significance in the Cuscuta genus

BACKGROUND

Cuscuta, a parasitic plant species in the Convolvulaceae family, grows in many countries and regions. However, the relationship between some species is still unclear. Therefore, more studies are needed to assess the variation of the chloroplast (cp) genome in Cuscuta species and their relationship with subgenera or sections, thus, providing important information on the evolution of Cuscuta species.

RESULTS

In the present study, we identified the whole cp genomes of C. epithymum, C. europaea, C. gronovii, C. chinensis and C. japonica, and then constructed a phylogenetic tree of 23 Cuscuta species based on the complete genome sequences and protein-coding genes. The complete cp genome sequences of C. epithymum and C. europaea were 96,292 and 97,661 bp long, respectively, and lacked an inverted repeat region. Most cp genomes of Cuscuta spp. have tetragonal and circular structures except for C. epithymum, C. europaea, C. pedicellata and C. approximata. Based on the number of genes and the structure of cp genome and the patterns of gene reduction, we found that C. epithymum and C. europaea belonged to subgenus Cuscuta. Most of the cp genomes of the 23 Cuscuta species had single nucleotide repeats of A and T. The inverted repeat region boundaries among species were similar in the same subgenera. Several cp genes were lost. In addition, the numbers and types of the lost genes in the same subgenus were similar. Most of the lost genes were related to photosynthesis (ndh, rpo, psa, psb, pet, and rbcL), which could have gradually caused the plants to lose the ability to photosynthesize.

CONCLUSION

Our results enrich the data on cp. genomes of genus Cuscuta. This study provides new insights into understanding the phylogenetic relationships and variations in the cp genome of Cuscuta species.

Intraspecific Chloroplast Genome Variation and Domestication Origins of Major Cultivars of Styphnolobium japonicum

Styphnolobium japonicum is a significant resource of ornamental and medicinal plants. In this study, we employed high-throughput sequencing to assemble nine chloroplast genomes of S. japonicum. We compared and reconstructed the phylogenetic relationships of these genomes, along with three publicly available chloroplast genomes. Our results showed that the length of the 12 S. japonicum chloroplast genomes ranged from 158,613 bp to 158,837 bp, all containing 129 unique functional genes. The genetic diversity within S. japonicum chloroplast genomes was relatively low, with π = 0.00029, Theta-W = 0.00028, and an indel frequency of 0.62 indels/1 kb. Among the four regions, the SSC region exhibited the highest genetic diversity and indel frequency, while the IR region had the lowest. Non-coding regions displayed greater genetic variation compared to coding regions, with a few highly variable regions identified. The phylogenetic tree constructed revealed that the major cultivars of S. japonicum originated from two genetic 'sources. S. japonicum 'JinhuaiJ2' had an independent origin and showed close relatedness to S. japonicum var. violacea, S. japonicum var. japonicum, and S. japonicum f. oligophylla. On the other hand, other major cultivars shared a common genetic origin and were closely related to S. japonicum f. pendula. This study highlights the variability of chloroplast genomes within S. japonicum and provides insights into the genetic origins of major cultivars and their relationships with different varieties and forma.

Extinction Threat to a Previously Undescribed Species of Gall Wasp (Hymenoptera: Cynipidae) and Two Associated Parasitoid Species (Hymenoptera: Braconidae and Eulophidae) on a Threatened Rose

Diplolepis ogawai Abe and Ide sp. nov. (Hymenoptera: Cynipidae) induces galls on Rosa hirtula (Regel) Nakai (Rosales: Rosaceae), which is endemic to a restricted area of Honshu, the main island of Japan. The gall is induced mainly on the leaf of R. hirtula in spring and the mature gall falls to the ground in early summer. The gall-inducing wasp emerges from the gall on the ground in the following spring, suggesting that D. ogawai is univoltine. From spring to summer, the braconid Syntomernus flavus Samartsev and Ku and the eulophid Aprostocetus sp. are parasitic on the larva of D. ogawai in the gall, and the adult wasp of both parasitoid species emerges from the gall on the ground in summer. For S. flavus, this is the first distribution record in Japan and the first host record. Since R. hirtula is threatened with extinction by succession and deforestation, D. ogawai and its two parasitoid wasp species are considered to be at risk of coextinction with the threatened rose. In the event that the population size of this rose species is further reduced, D. ogawai and its parasitoids may -become extinct prior to the extinction of R. hirtula. To conserve these three wasp species associated with R. hirtula, protection of remnant vegetation where individuals of this threatened rose species grow is necessary.

Complete Chloroplast Genome Sequence of the Long Blooming Cultivar Camellia 'Xiari Qixin': Genome Features, Comparative and Phylogenetic Analysis

The camellia flower is a famous woody plant with a long-cultivated history and high ornamental value. It is extensively planted and utilized around the world and owns a massive germplasm resource. Camellia 'Xiari Qixin' belongs to one of the typical cultivars in the four seasons camellia hybrids series. Due to its long flowering period, this kind of cultivar is identified as a precious resource of camellia flowers. In this study, the complete chloroplast genome sequence of C. 'Xiari Qixin' was first reported. Its whole chloroplast genome is 157,039 bp in length with an overall GC content of 37.30%, composed of a large single copy region (LSC, 86,674 bp), a small single copy region (SSC, 18,281 bp), and a pair of inverted repeat regions (IRs, 26,042 bp each). A total of 134 genes were predicted in this genome, including 8 ribosomal RNA genes, 37 transfer RNA genes, and 89 protein-coding genes. In addition, 50 simple sequence repeats (SSRs) and 36 long repeat sequences were detected. By comparing C. 'Xiari Qixin' and seven Camellia species on the chloroplast genome, seven mutation hotspot regions were identified, including psbK, trnS (GCU)-trnG(GCC), trnG(GCC), petN-psbM, trnF(GAA)-ndhJ, trnP(UGG)-psaJ, and ycf1. Phylogenetic analysis of 30 chloroplast genomes showed that the genetic relationship between C. 'Xiari Qixin' and Camellia azalea is quite close in evolution. These results could not only provide a valuable database for determining the maternal origin of Camellia cultivars, but also contribute to the exploration of the phylogenetic relationship and utilization of germplasm resources for Camellia.

The complete chloroplast genome sequence of Berchemia racemosa Siebold & Zucc. (Rhamnaceae), a rare plant species in Korea

Berchemia racemosa Siebold & Zucc., 1845 is a rare species distributed in restricted areas in the western Korean peninsula. In this study, the complete chloroplast genome (plastome) of B. racemosa was sequenced and assembled by Illumina paired-end sequencing. The plastome of B. racemosa was 161,187 bp in length and was quadripartite in structure, including a large single-copy (LSC) region of 89,503 bp, a small single-copy (SSC) region of 18,214 bp, and two inverted repeats of 26,735 bp. The GC content was 37.2%. The plastome of B. racemosa contains 130 genes, including eight ribosomal RNA (rRNA) genes, 37 transfer RNA (tRNA) genes, and 85 protein-coding genes. Phylogenetic analysis using complete genome sequences showed that B. racemosa is most closely related to Berchemia flavescens.

The Complete Chloroplast Genome Sequence of Machilus chuanchienensis (Lauraceae): Genome Structure and Phylogenetic Analysis

Machilus chuanchienensis is an ecological tree distributed in southwestern China. It has a significant valuation with making Hawk tea using its leaves, an ethnic traditional tea-like beverage with a long history in Chinese tea culture. The whole chloroplast (cp) genome is an ideal model for the phylogenetic study of Lauraceae because of its simple structure and highly conserved features. There have been numerous reports of complete cp genome sequences in Lauraceae, but little is known about M. chuanchienensis. Here, the next-generation sequencing (NGS) was used to sequence the M. chuanchienensis cp genome. Then, a comprehensive comparative genome analysis was performed. The results revealed that the M. chuanchienensis's cp genome measured 152,748 base pairs (bp) with a GC content of 39.15% and coded 126 genes annotated, including comprising eight ribosomal RNA (rRNA), 36 transporter RNA (tRNA), and 82 protein-coding genes. In addition, the cp genome presented a typical quadripartite structure comprising a large single-copy (LSC; 93,811) region, a small single-copy (SSC; 18,803) region, and the inverted repeats (IRs; 20,067) region and contained 92 simple sequence repeat (SSR) locus in total. Phylogenetic relationships of 37 species indicated that M. chuanchienensis was a sister to M. balansae, M. melanophylla, and M. minutiflora. Further research on this crucial species may benefit significantly from these findings.

Structural Diversities and Phylogenetic Signals in Plastomes of the Early-Divergent Angiosperms: A Case Study in Saxifragales

As representative of the early-divergent groups of angiosperms, Saxifragales is extremely divergent in morphology, comprising 15 families. Within this order, our previous case studies observed significant structural diversities among the plastomes of several lineages, suggesting a possible role in elucidating their deep phylogenetic relationships. Here, we collected 208 available plastomes from 11 constituent families to explore the evolutionary patterns among Saxifragales. With thorough comparisons, the losses of two genes and three introns were found in several groups. Notably, 432 indel events have been observed from the introns of all 17 plastomic intron-containing genes, which could well play an important role in family barcoding. Moreover, numerous heterogeneities and strong intrafamilial phylogenetic implications were revealed in pttRNA (plastomic tRNA) structures, and the unique structural patterns were also determined for five families. Most importantly, based on the well-supported phylogenetic trees, evident phylogenetic signals were detected in combinations with the identified pttRNAs features and intron indels, demonstrating abundant lineage-specific characteristics for Saxifragales. Collectively, the results reported here could not only provide a deeper understanding into the evolutionary patterns of Saxifragales, but also provide a case study for exploring the plastome evolution at a high taxonomic level of angiosperms.

Complete chloroplast genomes of Sorbus sensu stricto (Rosaceae): comparative analyses and phylogenetic relationships

BACKGROUND

Sorbus sensu stricto (Sorbus s.s.) is a genus with important economical values because of its beautiful leaves, and flowers and especially the colorful fruits. It belongs to the tribe Maleae of the family Rosaceae, and comprises about 90 species mainly distributed in China. There is on-going dispute about its infrageneric classification and species delimitation as the species are morphologically similar. With the aim of shedding light on the circumscription of taxa within the genus, phylogenetic analyses were performed using 29 Sorbus s.s. chloroplast (cp) genomes (16 newly sequenced) representing two subgenera and eight sections.

RESULTS

The 16 cp genomes newly sequenced range between 159,646 bp and 160,178 bp in length. All the samples examined and 22 taxa re-annotated in Sorbus sensu lato (Sorbus s.l.) contain 113 unique genes with 19 of these duplicated in the inverted repeat (IR). Six hypervariable regions including trnR-atpA, petN-psbM, rpl32-trnL, trnH-psbA, trnT-trnL and ndhC-trnV were screened and 44-53 SSRs and 14-31 dispersed repeats were identified as potential molecular markers. Phylogenetic analyses under ML/BI indicated that Sorbus s.l. is polyphyletic, but Sorbus s.s. and the other five segregate genera, Aria, Chamaemespilus, Cormus, Micromeles and Torminalis are monophyletic. Two major clades and four sub-clades resolved with full-support within Sorbus s.s. are not consistent with the existing infrageneric classification. Two subgenera, subg. Sorbus and subg. Albocarmesinae are supported as monophyletic when S. tianschanica is transferred to subg. Albocarmesinae from subg. Sorbus and S. hupehensis var. paucijuga transferred to subg. Sorbus from subg. Albocarmesinae, respectively. The current classification at sectional level is not supported by analysis of cp genome phylogeny.

CONCLUSION

Phylogenomic analyses of the cp genomes are useful for inferring phylogenetic relationships in Sorbus s.s. Though genome structure is highly conserved in the genus, hypervariable regions and repeat sequences used are the most promising molecule makers for population genetics, species delimitation and phylogenetic studies.

Comprehensive analysis of complete chloroplast genome and phylogenetic aspects of ten Ficus species

BACKGROUND

The large genus Ficus comprises approximately 800 species, most of which possess high ornamental and ecological values. However, its evolutionary history remains largely unknown. Plastome (chloroplast genome) analysis had become an essential tool for species identification and for unveiling evolutionary relationships between species, genus and other rank groups. In this work we present the plastomes of ten Ficus species.

RESULTS

The complete chloroplast (CP) genomes of eleven Ficus specimens belonging to ten species were determined and analysed. The full length of the Ficus plastome was nearly 160 kbp with a similar overall GC content, ranging from 35.88 to 36.02%. A total of 114 unique genes, distributed in 80 protein-coding genes, 30 tRNAs, and 4 rRNAs, were annotated in each of the Ficus CP genome. In addition, these CP genomes showed variation in their inverted repeat regions (IR). Tandem repeats and mononucleotide simple sequence repeat (SSR) are widely distributed across the Ficus CP genome. Comparative genome analysis showed low sequence variability. In addition, eight variable regions to be used as potential molecular markers were proposed for future Ficus species identification. According to the phylogenetic analysis, these ten Ficus species were clustered together and further divided into three clades based on different subgenera. Simultaneously, it also showed the relatedness between Ficus and Morus.

CONCLUSION

The chloroplast genome structure of 10 Ficus species was similar to that of other angiosperms, with a typical four-part structure. Chloroplast genome sizes vary slightly due to expansion and contraction of the IR region. And the variation of noncoding regions of the chloroplast genome is larger than that of coding regions. Phylogenetic analysis showed that these eleven sampled CP genomes were divided into three clades, clustered with species from subgenus Urostigma, Sycomorus, and Ficus, respectively. These results support the Berg classification system, in which the subgenus Ficus was further decomposed into the subgenus Sycomorus. In general, the sequencing and analysis of Ficus plastomes, especially the ones of species with no or limited sequences available yet, contribute to the study of genetic diversity and species evolution of Ficus, while providing useful information for taxonomic and phylogenetic studies of Ficus.

Infrageneric Plastid Genomes of Cotoneaster (Rosaceae): Implications for the Plastome Evolution and Origin of C. wilsonii on Ulleung Island

Cotoneaster is a taxonomically and ornamentally important genus in the family Rosaceae; however, phylogenetic relationships among its species are complicated owing to insufficient morphological diagnostic characteristics and hybridization associated with polyploidy and apomixis. In this study, we sequenced the complete plastomes of seven Cotoneaster species (C. dielsianus, C. hebephyllus, C. integerrimus, C. mongolicus, C. multiflorus, C. submultiflorus, and C. tenuipes) and included the available complete plastomes in a phylogenetic analysis to determine the origin of C. wilsonii, which is endemic to Ulleung Island, Korea. Furthermore, based on 15 representative lineages within the genus, we carried out the first comparative analysis of Cotoneaster plastid genomes to gain an insight into their molecular evolution. The plastomes were highly conserved, with sizes ranging from 159,595 bp (C. tenuipes) to 160,016 bp (C. hebephyllus), and had a GC content of 36.6%. The frequency of codon usage showed similar patterns among the 15 Cotoneaster species, and 24 of the 35 protein-coding genes were predicted to undergo RNA editing. Eight of the 76 common protein-coding genes, including ccsA, matK, ndhD, ndhF, ndhK, petA, rbcL, and rpl16, were positively selected, implying their potential roles in adaptation and speciation. Of the 35 protein-coding genes, 24 genes (15 photosynthesis-related, seven self-replications, and three others) were found to harbor RNA editing sites. Furthermore, several mutation hotspots were identified, including trnG-UCC/trnR-UCU/atpA and trnT-UGU/trnL-UAA. Maximum likelihood analysis based on 57 representative plastomes of Cotoneaster and two Heteromeles plastomes as outgroups revealed two major lineages within the genus, which roughly correspond to two subgenera, Chaenopetalum and Cotoneaster. The Ulleung Island endemic, C. wilsonii, shared its most recent common ancestor with two species, C. schantungensis and C. zabelii, suggesting its potential origin from geographically close members of the subgenus Cotoneaster, section Integerrimi.

Comparative Chloroplast Genome Analysis of Wax Gourd (Benincasa hispida) with Three Benincaseae Species, Revealing Evolutionary Dynamic Patterns and Phylogenetic Implications

Benincasa hispida (wax gourd) is an important Cucurbitaceae crop, with enormous economic and medicinal importance. Here, we report the de novo assembly and annotation of the complete chloroplast genome of wax gourd with 156,758 bp in total. The quadripartite structure of the chloroplast genome comprises a large single-copy (LSC) region with 86,538 bp and a small single-copy (SSC) region with 18,060 bp, separated by a pair of inverted repeats (IRa and IRb) with 26,080 bp each. Comparison analyses among B. hispida and three other species from Benincaseae presented a significant conversion regarding nucleotide content, genome structure, codon usage, synonymous and non-synonymous substitutions, putative RNA editing sites, microsatellites, and oligonucleotide repeats. The LSC and SSC regions were found to be much more varied than the IR regions through a divergent analysis of the species within Benincaseae. Notable IR contractions and expansions were observed, suggesting a difference in genome size, gene duplication and deletion, and the presence of pseudogenes. Intronic gene sequences, such as trnR-UCU–atpA and atpH–atpI, were observed as highly divergent regions. Two types of phylogenetic analysis based on the complete cp genome and 72 genes suggested sister relationships between B. hispida with the Citrullus, Lagenaria, and Cucumis. Variations and consistency with previous studies regarding phylogenetic relationships are discussed. The cp genome of B. hispida provides valuable genetic information for the detection of molecular markers, research on taxonomic discrepancies, and the inference of the phylogenetic relationships of Cucurbitaceae.

Complete chloroplast genomes of Rubus species (Rosaceae) and comparative analysis within the genus

BACKGROUND

Rubus is the largest genus of the family Rosaceae and is valued as medicinal, edible, and ornamental plants. Here, we sequenced and assembled eight chloroplast (cp) genomes of Rubus from the Dabie Mountains in Central China. Fifty-one Rubus species were comparatively analyzed for the cp genomes including the eight newly discovered genomes and forty-three previously reported in GenBank database (NCBI).

RESULTS

The eight newly obtained cp genomes had the same quadripartite structure as the other cp genomes in Rubus. The length of the eight plastomes ranged from 155,546 bp to 156,321 bp with similar GC content (37.0 to 37.3%). The results indicated 133-134 genes were annotated for the Rubus plastomes, which contained 88 or 89 protein coding genes (PCGs), 37 transfer RNA genes (tRNAs), and eight ribosomal RNA genes (rRNAs). Among them, 16 (or 18) of the genes were duplicated in the IR region. Structural comparative analysis results showed that the gene content and order were relatively preserved. Nucleotide variability analysis identified nine hotspot regions for genomic divergence and multiple simple sequences repeats (SSRs), which may be used as markers for genetic diversity and phylogenetic analysis. Phylogenetic relationships were highly supported within the family Rosaceae, as evidenced by sub-clade taxa cp genome sequences.

CONCLUSION

Thus, the whole plastome may be used as a super-marker in phylogenetic studies of this genus.

The investigation of intraspecific characteristics and comparative analyses of the complete mitochondrial genome of Stegobium paniceum (Linnaeus, 1758) (Coleoptera: Ptinidae) assembled from public NGS raw reads of the black truffle, Tuber melanosporum

Due to the rapid development of NGS technologies, a huge amount of NGS raw reads have been accumulated in public repositories, such as the Short Read Archive of NCBI. We successfully rescued the complete mitochondrial genome of Stegobium paniceum, a drug store beetle, from public NGS raw reads of truffle generated from the whole genome project. The circular mitogenome of S. paniceum is 15,474 bp long including 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNAs, and a single large non-coding region of 803 bp. Intraspecific transfer RNAs structure and sequence variations were investigated and simple sequence repeats identified from three S. paniceum mitochondrial genomes were compared showing their diversities as fundamental data to utilize them in various aspects including developing efficient molecular markers in the family, Ptinidae. Phylogenetic analysis of 23 Bostrichoidea mitochondrial genomes presented better species identification based on phylogenetic analyses and the optimal options for constructing phylogenetic trees based on Bostrichoidea mitochondrial genomes. Our results present not only utilization of public NGS raw read sequences but also intraspecific features of S. paniceum mitochondrial genomes and comparative analysis of Bostrichoidea mitochondrial genomes in various aspects.

Chloroplast Phylogenomic Analyses Reveal a Maternal Hybridization Event Leading to the Formation of Cultivated Peanuts

Peanuts (Arachis hypogaea L.) offer numerous healthy benefits, and the production of peanuts has a prominent role in global food security. As a result, it is in the interest of society to improve the productivity and quality of peanuts with transgenic means. However, the lack of a robust phylogeny of cultivated and wild peanut species has limited the utilization of genetic resources in peanut molecular breeding. In this study, a total of 33 complete peanut plastomes were sequenced, analyzed and used for phylogenetic analyses. Our results suggest that sect. Arachis can be subdivided into two lineages. All the cultivated species are contained in Lineage I with AABB and AA are the two predominant genome types present, while species in Lineage II possess diverse genome types, including BB, KK, GG, etc. Phylogenetic studies also indicate that all allotetraploid cultivated peanut species have been derived from a possible maternal hybridization event with one of the diploid Arachis duranensis accessions being a potential AA sub-genome ancestor. In addition, Arachis monticola, a tetraploid wild species, is placed in the same group with all the cultivated peanuts, and it may represent a transitional species, which has been through the recent hybridization event. This research could facilitate a better understanding of the taxonomic status of various Arachis species/accessions and the evolutionary relationship among them, and assists in the correct and efficient use of germplasm resources in breeding efforts to improve peanuts for the benefit of human beings.

Variation and Evolution of the Whole Chloroplast Genomes of Fragaria spp. (Rosaceae)

Species identification is vital for protecting species diversity and selecting high-quality germplasm resources. Wild Fragaria spp. comprise rich and excellent germplasm resources; however, the variation and evolution of the whole chloroplast (cp) genomes in the genus Fragaria have been ignored. In the present study, 27 complete chloroplast genomes of 11 wild Fragaria species were sequenced using the Illumina platform. Then, the variation among complete cp genomes of Fragaria was analyzed, and phylogenetic relationships were reconstructed from those genome sequences. There was an overall high similarity of sequences, with some divergence. According to analysis with mVISTA, non-coding regions were more variable than coding regions. Inverted repeats (IRs) were observed to contract or expand to different degrees, which resulted in different sizes of cp genomes. Additionally, five variable loci, trnS-trnG, trnR-atpA, trnC-petN, rbcL-accD, and psbE-petL, were identified that could be used to develop DNA barcoding for identification of Fragaria species. Phylogenetic analyses based on the whole cp genomes supported clustering all species into two groups (A and B). Group A species were mainly distributed in western China, while group B contained several species from Europe and Americas. These results support allopolyploid origins of the octoploid species F. chiloensis and F. virginiana and the tetraploid species F. moupinensis and F. tibetica. The complete cp genomes of these Fragaria spp. provide valuable information for selecting high-quality Fragaria germplasm resources in the future.

The complete chloroplast genome sequence of Rubus hirsutus Thunb. and a comparative analysis within Rubus species

Rubus hirsutus is a type of tonifying kidney-essence herb that belongs to the Rosaceae family, and has been commonly used to treat multiple diseases, such as polyuria, impotence, and infertility. In this study, we determined the complete chloroplast sequence of R. hirsutus and conduced a comparative analysis within the genus Rubus. The assembled chloroplast (cp.) genome is 156,380 bp in length with a GC content of 37.0% and shares a conserved quadripartite structure within the other cp. genomes in this genus. A total of 132 unique genes were annotated in the cp. genome of R. hirsutus, which contained 87 protein-coding genes, 37 tRNAs, and eight rRNAs. Seventeen duplicated genes were identified in the inverted repeats region. Furthermore, 70 simple sequence repeats and 35 long repeats were detected in total in the R. hirsutus chloroplast genome. Eight mutational hotspots were identified in the cp. genome of this species with higher nucleotide variations in non-coding regions than those of coding regions. Furthermore, the gene order, codon usage, and repeat sequence distribution were highly consistent in Rubus according to the results of a comparative analysis. A phylogenetic analysis indicated that there was a sister relationship between R. hirsutus and R. chingii. Overall, the complete chloroplast genome of R. hirsutus and the comparative analysis will help to further the evolutionary study, conservation, phylogenetic reconstruction, and development of molecular barcodes for the genus Rubus.

Plastome Characterization and Phylogenomic Analysis Yield New Insights into the Evolutionary Relationships among the Species of the Subgenus Bryocles (Hosta; Asparagaceae) in East Asia

The genus Hosta, which has a native distribution in temperate East Asia and a number of species ranging from 23 to 40, represents a taxonomically important and ornamentally popular plant. Despite its taxonomic and horticultural importance, the genus Hosta has remained taxonomically challenging owing to insufficient diagnostic features, continuous morphological variation, and the process of hybridization and introgression, making species circumscription and phylogenetic inference difficult. In this study, we sequenced 11 accessions of Hosta plastomes, including members of three geographically defined subgenera, Hosta, Bryocles, and Giboshi, determined the characteristics of plastomes, and inferred their phylogenetic relationships. We found highly conserved plastomes among the three subgenera, identified several mutation hotspots that can be used as barcodes, and revealed the patterns of codon usage bias and RNA editing sites. Five positively selected plastome genes (rbcL, rpoB, rpoC2, rpl16, and rpl20) were identified. Phylogenetic analysis suggested (1) the earliest divergence of subg. Hosta, (2) non-monophyly of subg. Bryocles and its two sections (Lamellatae and Stoloniferae), (3) a sister relationship between H. sieboldiana (subg. Giboshi) and H. ventricosa (subg. Bryocles), and (4) reciprocally monophyletic and divergent lineages of H. capitata in Korea and Japan, requiring further studies of their taxonomic distinction.

A Comparative Analyses of the Complete Mitochondrial Genomes of Fungal Endosymbionts in Sogatella furcifera, White-Backed Planthoppers

Sogatella furcifera Horvath, commonly known as the white-backed planthoppers (WBPH), is an important pest in East Asian rice fields. Fungal endosymbiosis is widespread among planthoppers in the infraorder Fulgoromorpha and suborder Auchenorrhyncha. We successfully obtained complete mitogenome of five WBPH fungal endosymbionts, belonging to the Ophiocordycipitaceae family, from next-generation sequencing (NGS) reads obtained from S. furcifera samples. These five mitogenomes range in length from 55,390 bp to 55,406 bp, which is shorter than the mitogenome of the fungal endosymbiont found in Ricania speculum, black planthoppers. Twenty-eight protein-coding genes (PCGs), 12 tRNAs, and 2 rRNAs were found in the mitogenomes. Two single-nucleotide polymorphisms, two insertions, and three deletions were identified among the five mitogenomes, which were fewer in number than those of four species of Ophiocordycipitaceae, Ophiocordyceps sinensis, Hirsutella thompsonii, Hirsutella rhossiliensis, and Tolypocladium inflatum. Noticeably short lengths (up to 18 bp) of simple sequence repeats were identified in the five WBPH fungal endosymbiont mitogenomes. Phylogenetic analysis based on conserved PCGs across 25 Ophiocordycipitaceae mitogenomes revealed that the five mitogenomes were clustered with that of R. speculum, forming an independent clade. In addition to providing the full mitogenome sequences, obtaining complete mitogenomes of WBPH endosymbionts can provide insights into their phylogenetic positions without needing to isolate the mtDNA from the host. This advantage is of value to future studies involving fungal endosymbiont mitogenomes.

Characterization and comparative analysis among plastome sequences of eight endemic Rubus (Rosaceae) species in Taiwan

Genus Rubus represents the second largest genus of the family Rosaceae in Taiwan, with 41 currently recognized species across three subgenera (Chamaebatus, Idaoeobatus, and Malochobatus). Despite previous morphological and cytological studies, little is known regarding the overall phylogenetic relationships among the Rubus species in Taiwan, and their relationships to congeneric species in continental China. We characterized eight complete plastomes of Taiwan endemic Rubus species: subg. Idaeobatus (R. glandulosopunctatus, R. incanus, R. parviaraliifolius, R rubroangustifolius, R. taitoensis, and R. taiwanicolus) and subg. Malachobatus (R. kawakamii and R. laciniastostipulatus) to determine their phylogenetic relationships. The plastomes were highly conserved and the size of the complete plastome sequences ranged from 155,566 to 156,236 bp. The overall GC content ranged from 37.0 to 37.3%. The frequency of codon usage showed similar patterns among species, and 29 of the 73 common protein-coding genes were positively selected. The comparative phylogenomic analysis identified four highly variable intergenic regions (rps16/trnQ, petA/psbJ, rpl32/trnL-UAG, and trnT-UGU/trnL-UAA). Phylogenetic analysis of 31 representative complete plastomes within the family Rosaceae revealed three major lineages within Rubus in Taiwan. However, overall phylogenetic relationships among endemic species require broader taxon sampling to gain new insights into infrageneric relationships and their plastome evolution.

Plastome Characterization and Phylogenomics of East Asian Beeches with a Special Emphasis on Fagus multinervis on Ulleung Island, Korea

Beech trees of the genus Fagus (Fagaceae) are monoecious and distributed in the Northern Hemisphere. They represent an important component of mixed broad-leaved evergreen-deciduous forests and are an economically important source of timber. Despite their ecological and economical importance, however, little is known regarding the overall plastome evolution among Fagus species in East Asia. In particular, the taxonomic position and status of F. multinervis, a beech species endemic to Ulleung Island of Korea, remains unclear even today. Therefore, in this study, we characterized four newly completed plastomes of East Asian Fagus species (one accession each of F. crenata and F. multinervis and two accessions of F. japonica). Moreover, we performed phylogenomic analyses comparing these four plastomes with F. sylvatica (European beech) plastome. The four plastomes were highly conserved, and their size ranged from 158,163 to 158,348 base pair (bp). The overall GC content was 37.1%, and the sequence similarity ranged from 99.8% to 99.99%. Codon usage patterns were similar among species, and 7 of 77 common protein-coding genes were under positive selection. Furthermore, we identified five highly variable hotspot regions of the Fagus plastomes (ccsA/ndhD, ndhD/psaC, ndhF/rpl32, trnS-GCU/trnG-UCC, and ycf1). Phylogenetic analysis revealed the monophyly of Fagus as well as early divergence of the subgenus Fagus and monophyletic Engleriana. Finally, phylogenetic results supported the taxonomic distinction of F. multinervis from its close relatives F. engleriana and F. japonica. However, the sister species and geographic origin of F. multinervis on Ulleung Island could not be determined.

The Complete Chloroplast Genome of Arabidopsis thaliana Isolated in Korea (Brassicaceae): An Investigation of Intraspecific Variations of the Chloroplast Genome of Korean A. thaliana

Arabidopsis thaliana (L.) Heynh. is a model organism of plant molecular biology. More than 1,700 whole genome sequences have been sequenced, but no Korean isolate genomes have been sequenced thus far despite the fact that many A. thaliana isolated in Japan and China have been sequenced. To understand the genetic background of Korean natural A. thaliana (named as 180404IB4), we presented its complete chloroplast genome, which is 154,464 bp long and has four subregions: 85,164 bp of large single copy (LSC) and 17,781 bp of small single copy (SSC) regions are separated by 26,257 bp of inverted repeat (IRs) regions including 130 genes (85 protein-coding genes, eight rRNAs, and 37 tRNAs). Fifty single nucleotide polymorphisms (SNPs) and 14 insertion and deletions (INDELs) are identified between 180404IB4 and Col0. In addition, 101 SSRs and 42 extendedSSRs were identified on the Korean A. thaliana chloroplast genome, indicating a similar number of SSRs on the rest five chloroplast genomes with a preference of sequence variations toward the SSR region. A nucleotide diversity analysis revealed two highly variable regions on A. thaliana chloroplast genomes. Phylogenetic trees with three more chloroplast genomes of East Asian natural isolates show that Korean and Chinese natural isolates are clustered together, whereas two Japanese isolates are not clustered, suggesting the need for additional investigations of the chloroplast genomes of East Asian isolates.

Comprehensive Analyses of the Complete Mitochondrial Genome of Figulus binodulus (Coleoptera: Lucanidae)

Figulus binodulus Waterhouse is a small stag beetle distributed in East Asia. We determined the first mitochondrial genome of F. binodulus of which is 16,261-bp long including 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNAs, and a single large noncoding region of 1,717 bp. Gene order of F. binodulus is identical to the ancestral insect mitochondrial gene order as in most other stag beetle species. All of 22 tRNAs could be shaped into typical clover-leaf structure except trnSer1. Comparative analyses of 21 Lucanidae mitochondrial genomes was conducted in aspect of their length and AT-GC ratio. Nucleotide diversities analyses provide that cox1 and cox2 in Lucanidae are less diverse than those of Scarabaeoidea. Fifty simple sequence repeats (SSRs) were identified on F. binodulus mitochondrial genome. Comparative analysis of SSRs among five mitochondrial genomes displayed similar trend along with SSR types. Figulus binodulus was sister to all other available family Lucanidae species in the phylogenetic tree.

Molecular Evidence for Hybrid Origin and Phenotypic Variation of Rosa Section Chinenses

Rosa sect. Chinenses (Rosaceae) is an important parent of modern rose that is widely distributed throughout China and plays an important role in breeding and molecular biological research. R. sect. Chinenses has variable morphological traits and mixed germplasm. However, the taxonomic status and genetic background of sect. Chinenses varieties remain unclear. In this study, we collected germplasm resources from sect. Chinenses varieties with different morphological traits. Simple sequence repeat (SSR) markers, chloroplast markers, and single copy nuclear markers were used to explore the genetic background of these germplasm resources. We described the origin of hybridization of rose germplasm resources by combining different molecular markers. The results showed that the flower and hip traits of different species in R. sect. Chinenses were significantly different. The SSR analysis showed that the two wild type varieties have different genetic backgrounds. The double petal varieties of R. sect. Chinenses could be hybrids of two wild type varieties. A phylogenetic analysis showed that the maternal inheritance of sect. Chinenses varieties had two different origins. To some extent, variation in the morphological traits of double petal species of R. sect. Chinenses reflects the influence of cultivation process. This study emphasizes that different genetic markers vary in their characteristics. Therefore, analyzing different genetic markers in could provide an insight into highly heterozygous species.

ycf1-ndhF genes, the most promising plastid genomic barcode, sheds light on phylogeny at low taxonomic levels in Prunus persica

Background

Chloroplast genome sequencing is becoming a valuable process for developing several DNA barcodes. At present, plastid DNA barcode for systematics and evolution in flowering plant rely heavily on the use of non-coding genes. The present study was performed to verify the novelty and suitability of the two hotspot barcode plastid coding gene ycf1 and ndhF, to estimate the rate of molecular evolution in the Prunus genus at low taxonomic levels.

Results

Here, 25 chloroplast genomes of Prunus genus were selected for sequences annotation to search for the highly variable coding DNA barcode regions. Among them, 5 genera were of our own data, including the ornamental, cultivated, and wild haplotype, while 20 genera have been downloaded from the GenBank database. The results indicated that the two hotspot plastid gene ycf1 and ndhF were the most variable regions within the coding genes in Prunus with an average of 3268 to 3416 bp in length, which have been predicted to have the highest nucleotide diversity, with the overall transition/transversion bias (R = 1.06). The ycf1-ndhF structural domains showed a positive trend evident in structure variation among the 25 specimens tested, due to the variant overlap^’s gene annotation and insertion or deletion with a broad trend of the full form of IGS sequence. As a result, the principal component analysis (PCA) and the ML tree data drew an accurate monophyletic annotations cluster in Prunus species, offering unambiguous identification without overlapping groups between peach, almond, and cherry.

Conclusion

To this end, we put forward the domain of the two-locus ycf1-ndhF genes as the most promising coding plastid DNA barcode in P. persica at low taxonomic levels. We believe that the discovering of further variable loci with high evolutionary rates is extremely useful and potential uses as a DNA barcode in P. persica for further phylogeny study and species identification.

Molecular markers from the chloroplast genome of rose provide a complementary tool for variety discrimination and profiling

The rose is one of the most important ornamental woody plants because of its extensive use and high economic value. Herein, we sequenced a complete chloroplast genome of the miniature rose variety Rosa 'Margo Koster' and performed comparative analyses with sequences previously published for other species in the Rosaceae family. The chloroplast genome of Rosa 'Margo Koster', with a size of 157,395 bp, has a circular quadripartite structure typical of angiosperm chloroplast genomes and contains a total of 81 protein-coding genes, 30 tRNA genes and 4 rRNA genes. Conjunction regions in the chloroplast genome of Rosa 'Margo Koster' were verified and manually corrected by Sanger sequencing. Comparative genome analysis showed that the IR contraction and expansion events resulted in rps19 and ycf1 pseudogenes. The phylogenetic analysis within the Rosa genus showed that Rosa 'Margo Koster' is closer to Rosa odorata than to other Rosa species. Additionally, we identified and screened highly divergent sequences and cpSSRs and compared their power to discriminate rose varieties by Sanger sequencing and capillary electrophoresis. The results showed that 15 cpSSRs are polymorphic, but their discriminating power is only moderate among a set of rose varieties. However, more than 150 single nucleotide variations (SNVs) were discovered in the flanking region of cpSSRs, and the results indicated that these SNVs have a higher divergence and stronger power for profiling rose varieties. These findings suggest that nucleotide mutations in the chloroplast genome may be an effective and powerful tool for rose variety discrimination and DNA profiling. These molecular markers in the chloroplast genome sequence of Rosa spp. will facilitate population and phylogenetic studies and other related studies of this species.

Characterization and Comparison of Two Complete Plastomes of Rosaceae Species (Potentilla dickinsii var. glabrata and Spiraea insularis) Endemic to Ulleung Island, Korea

Potentilla dickinsii var. glabrata and Spiraea insularis in the family Rosaceae are species endemic to Ulleung Island, Korea, the latter of which is listed as endangered. In this study, we characterized the complete plastomes of these two species and compared these with previously reported plastomes of other Ulleung Island endemic species of Rosaceae (Cotoneaster wilsonii, Prunus takesimensis, Rubus takesimensis, and Sorbus ulleungensis). The highly conserved complete plastomes of P. dickinsii var. glabrata and S. insularis are 158,637 and 155,524 base pairs with GC contents of 37% and 36.9%, respectively. Comparative phylogenomic analysis identified three highly variable intergenic regions (trnT-UGU/trnL-UAA, rpl32/trnL-UAG, and ndhF/rpl32) and one variable genic region (ycf1). Only 14 of the 75 protein-coding genes have been subject to strong purifying selection. Phylogenetic analysis of 23 representative plastomes within the Rosaceae supported the monophyly of Potentilla and the sister relationship between Potentilla and Fragaria and indicated that S. insularis is sister to a clade containing Cotoneaster, Malus, Pyrus, and Sorbus. The plastome resources generated in this study will contribute to elucidating the plastome evolution of insular endemic Rosaceae on Ulleung Island and also in assessing the genetic consequences of anagenetic speciation for various endemic lineages on the island.

Comparative analyses of 32 complete plastomes of Tef (Eragrostis tef ) accessions from Ethiopia: phylogenetic relationships and mutational hotspots

Eragrostis tef is an important cereal crop in Ethiopia with excellent storage properties, high–quality food, and the unique ability to thrive in extreme environmental conditions. However, the application of advanced molecular tools for breeding and conservation of these species is extremely limited. Therefore, developing chloroplast genome resources and high-resolution molecular markers are valuable to E. tef population and biogeographic studies. In the current study, we assembled and compared the complete plastomes of 32 E. tef accessions. The size of the plastomes ranged from 134,349 to 134,437 bp with similar GC content (∼38.3%). Genomes annotations revealed 112 individual genes, including 77 protein-coding, 31 tRNA, and 4 rRNA genes. Comparison of E. tef plastomes revealed a low degree of intraspecific sequence variations and no structural differentiations. Furthermore, we found 34 polymorphic sites (13 cpSSRs, 12 InDels, and 9 SNPs) that can be used as valuable DNA barcodes. Among them, the majority (88%) of the polymorphic sites were identified in the noncoding genomic regions. Nonsynonymous (ka) and synonymous (ks) substitution analysis showed that all PCGs were under purifying selection (ka/ks <1). The phylogenetic analyses of the whole plastomes and polymorphic region sequences were able to distinguish the accession from the southern population, indicating its potential to be used as a super-barcode. In conclusion, the newly generated plastomes and polymorphic markers developed here could be a useful genomic resource in molecular breeding, population genetics and the biogeographical study of E. tef.

The complete chloroplast genome and characteristics analysis of Callistemon rigidus R.Br

Callistemon rigidus R.Br. one of the traditional Chinese medicinal plants, is acrid-flavored and mild-natured, with the prominent effects reducing swelling, resolving phlegm, and dispelling rheumatism. Clinically, it has been commonly used to treat cold, cough and asthma, pain and swelling from impact injuries, eczema, rheumatic arthralgia. The chloroplast genome study on Callistemon rigidus R.Br. is a few seen. This study demonstrates the data collected from the assembly and annotation of the chloroplast (cp) genome of Callistemon rigidus R.Br., followed by furthers comparative analysis with the cp genomes of closely related species. C. rigidus R.Br. showed a cp genome in the size of 158, 961 bp long with 36.78% GC content, among which a pair of inverted repeats (IRs) of 26, 671 bp separated a large single-copy (LSC) region of 87, 162 bp and a small single-copy (SSC) region of 18, 457 bp. Altogether 131 genes were hosted, including 37 transfer RNAs, 8 ribosomal RNAs, and 86 protein-coding genes. 284 simple sequence repeats (SSRs) were also marked out. A comparative analysis of the genome structure and the sequence data of closely related species unveiled the conserved gene order in the IR and LSC/SSC regions, a quite constructive finding for future phylogenetic research. Overall, this study providing C. rigidus R.Br. genomic resources could positively contribute to the evolutionary study and the phylogenetic reconstruction of Myrtaceae.

The complete chloroplast genome of Agrimonia pilosa Ledeb. isolated in Korea (Rosaceae): investigation of intraspecific variations on its chloroplast genomes

The complete chloroplast genome of Agrimonia pilosa Ledeb. isolated in Korea is 155,125 bp long (GC ratio is 36.9%) and has four subregions: 84,458 bp of large single copy (34.9%) and 18,737 bp of small single copy (30.4%) regions are separated by 25,965 bp of inverted repeat (42.6%) regions including 129 genes (84 protein-coding genes, eight rRNAs, and 37 tRNAs). 258 SNPs and 542 INDELs were identified as intraspecific variations against the partial genome (KY419942). Phylogenetic trees show that our chloroplast genome was clustered with the previous A. pilosa chloroplast genome.

Complete Chloroplast Genomes of 14 Mangroves: Phylogenetic and Comparative Genomic Analyses

Mangroves are a group of plant species that occupy the coastal intertidal zone and are major components of this ecologically important ecosystem. Mangroves belong to about twenty diverse families. Here, we sequenced and assembled chloroplast genomes of 14 mangrove species from eight families spanning five rosid orders and one asterid order: Fabales (Pongamia pinnata), Lamiales (Avicennia marina), Malpighiales (Excoecaria agallocha, Bruguiera sexangula, Kandelia obovata, Rhizophora stylosa, and Ceriops tagal), Malvales (Hibiscus tiliaceus, Heritiera littoralis, and Thespesia populnea), Myrtales (Laguncularia racemosa, Sonneratia ovata, and Pemphis acidula), and Sapindales (Xylocarpus moluccensis). These chloroplast genomes range from 149 kb to 168 kb in length. A conserved structure of two inverted repeats (IRa and IRb, ~25.8 kb), one large single-copy region (LSC, ~89.0 kb), and one short single-copy region (SSC, ~18.9 kb) as well as ~130 genes (85 protein-coding, 37 tRNAs, and 8 rRNAs) was observed. We found the lowest divergence in the IR regions among the four regions. We also identified simple sequence repeats (SSRs), which were found to be variable in numbers. Most chloroplast genes are highly conserved, with only four genes under positive selection or relaxed pressure. Combined with publicly available chloroplast genomes, we carried out phylogenetic analysis and confirmed the previously reported phylogeny within rosids, including the positioning of obscure families in Malpighiales. Our study reports 14 mangrove chloroplast genomes and illustrates their genome features and evolution.

Long-reads reveal that Rhododendron delavayi plastid genome contains extensive repeat sequences, and recombination exists among plastid genomes of photosynthetic Ericaceae

BACKGROUND

Rhododendron delavayi Franch. var. delavayi is a wild ornamental plant species in Guizhou Province, China. The lack of its plastid genome information seriously hinders the further application and conservation of the valuable resource.

METHODS

The complete plastid genome of R. delavayi was assembled from long sequence reads. The genome was then characterized, and compared with those of other photosynthetic Ericaceae species.

RESULTS

The plastid genome of R. delavayi has a typical quadripartite structure, and a length of 202,169 bp. It contains a large number of repeat sequences and shows preference for codon usage. The comparative analysis revealed the irregular recombination of gene sets, including rearrangement and inversion, in the large single copy region. The extreme expansion of the inverted repeat region shortened the small single copy, and expanded the full length of the genome. In addition, consistent with traditional taxonomy, R. delavayi with nine other species of the same family were clustered into Ericaceae based on the homologous protein-coding sequences of the plastid genomes. Thus, the long-read assembly of the plastid genome of R. delavayi would provide basic information for the further study of the evolution, genetic diversity, and conservation of R. delavayi and its relatives.

The chloroplasts genomic analyses of Rosa laevigata, R. rugosa and R. canina

BACKGROUND

Many species of the genus Rosa have been used as ornamental plants and traditional medicines. However, industrial development of roses is hampered due to highly divergent characteristics.

METHODS

We analyzed the chloroplast (cp) genomes of Rosa laevigata, R. rugosa and R. canina, including the repeat sequences, inverted-repeat (IR) contractions and expansions, and mutation sites.

RESULTS

The size of the cp genome of R. laevigata, R. rugosa and R. canina was between 156 333 bp and 156 533 bp, and contained 113 genes (30 tRNA genes, 4 rRNA genes and 79 protein-coding genes). The regions with a higher degree of variation were screened out (trnH-GUU, trnS-GCU, trnG-GCC, psbA-trnH, trnC-GCA,petN, trnT-GGU, psbD, petA, psbJ, ndhF, rpl32,psaC and ndhE). Such higher-resolution loci lay the foundation of barcode-based identification of cp genomes in Rosa genus. A phylogenetic tree of the genus Rosa was reconstructed using the full sequences of the cp genome. These results were largely in accordance with the current taxonomic status of Rosa.

CONCLUSIONS

Our data: (i) reveal that cp genomes can be used for the identification and classification of Rosa species; (ii) can aid studies on molecular identification, genetic transformation, expression of secondary metabolic pathways and resistant proteins; (iii) can lay a theoretical foundation for the discovery of disease-resistance genes and cultivation of Rosa species.

Comparison of the complete plastomes and the phylogenetic analysis of Paulownia species

Paulownia species are important ecological, economic and ornamental species, but their phylogenetic relationship remains unclear, which seriously affects the development and utilization of these important resources. The complete chloroplast genomes of six Paulownia species were assembled by next-generation sequencing data. By adding two known Paulownia chloroplast genomes to these six assembled genomes, we performed the comparative analysis and phylogenetic tree reconstruction of Paulownia. The results indicated that the chloroplast genomes of Paulownia species ranged in size from 154,107 to 154,694 bp. These chloroplast genomes contained 117 unique functional genes, including 80 protein-coding genes, four rRNA genes, and 33 tRNA genes. Twelve hotspot regions, five protein-coding genes and seven noncoding regions, were identified in the chloroplast genomes that showed high levels of sequence variation. Additionally, positive selection was observed in three genes, rps2, rbcL and ndhG. The maximum likelihood (ML) and Bayesian (BI) analysis strongly supported the monophyletic origin of Paulownia species, which clustered into two major clades: One clade included P. coreana, P. tomentosa and P. kawakamii, while the other clade comprised the 5 other species including P. fargesii and P. australis. This study provides useful genetic information for phylogenetic reconstruction, taxonomic discrepancies, and studying species evolution and phylogeography in Paulownia.

Complete Chloroplast Genomes of Vachellia nilotica and Senegalia senegal: Comparative Genomics and Phylogenomic Placement in a New Generic System

Vachellia and Senegalia are the most important genera in the subfamily Mimosoideae (Fabaceae). Recently, species from both genera were separated from the long-characterized Acacia due to their macro-morphological characteristics. However, this morpho-taxonomic differentiation struggles to discriminate some species, for example, Vachellia nilotica and Senegalia senegal. Therefore, sequencing the chloroplast (cp) genomes of these species and determining their phylogenetic placement via conserved genes may help to validate the taxonomy. Hence, we sequenced the cp genomes of V. nilotica and S. senegal, and the results showed that the sizes of the genomes are 165.3 and 162.7 kb, respectively. The cp genomes of both species comprised large single-copy regions (93,849~91,791 bp) and pairs of inverted repeats (IR; 26,093~26,008 bp). The total numbers of genes found in the V. nilotica and S. senegal cp genomes were 135 and 132, respectively. Approximately 123:130 repeats and 290:281 simple sequence repeats were found in the S. senegal and V. nilotica cp genomes, respectively. Genomic characterization was undertaken by comparing these genomes with those of 17 species belonging to related genera in Fabaceae. A phylogenetic analysis of the whole genome dataset and 56 shared genes was undertaken by generating cladograms with the same topologies and placing both species in a new generic system. These results support the likelihood of identifying segregate genera from Acacia with phylogenomic disposition of both V. nilotica and S. senegal in the subfamily Mimosoideae. The current study is the first to obtain complete genomic information on both species and may help to elucidate the genome architecture of these species and evaluate the genetic diversity among species.

Comparison of the Whole-Plastome Sequence between the Bonin Islands Endemic Rubus boninensis and Its Close Relative, Rubus trifidus (Rosaceae), in the Southern Korean Peninsula

Rubus boninensis is a rare endemic species found on the Bonin Islands with a very restricted distribution. It is morphologically most closely related to Rubus trifidus, occurring widely in the southern Korean peninsula and Japan. This species pair provides a good example of anagenetic speciation on an oceanic island in the northwestern Pacific Ocean-R. trifidus as a continental progenitor and R. boninensis as an insular derivative species. In this study, we firstly characterized the complete plastome of R. boninensis and R. trifidus and compared this species pair to another anagenetically derived species pair (R. takesimensis-R. crataegifolius). The complete plastome of R. trifidus was 155,823 base pairs (bp) long, slightly longer (16 bp) than that of R. boninensis (155,807 bp). No structural or content rearrangements were found between the species pair. Eleven hotspot regions, including trnH/psbA, were identified between R. trifidus and R. boninensis. Phylogenetic analysis of 19 representative plastomes within the family Rosaceae suggested sister relationships between R. trifidus and R. boninensis, and between R. crataegifolius and R. takesimensis. The plastome resources generated by the present study will help elucidate plastome evolution and resolve phylogenetic relationships within highly complex and reticulated lineages of the genus Rubus.

Evolutionary analysis of six chloroplast genomes from three Persea americana ecological races: Insights into sequence divergences and phylogenetic relationships

Chloroplasts significantly influence species phylogenies because of their maternal inheritance and the moderate evolutionary rate of their genomes. Avocado, which is a member of the family Lauraceae, has received considerable attention from botanists, likely because of its position as a basal angiosperm. However, there is relatively little avocado genomic information currently available. In this study, six complete avocado chloroplast genomes from three ecological races were assembled to examine the sequence diversity among the three avocado ecological races. A comparative genomic analysis revealed that 515 simple sequence repeat loci and 176 repeats belonging to four other types were polymorphic across the six chloroplast genomes. Three highly variable regions (trnC-GCA-petN, petN-psbM, and petA-psbJ) were identified as highly informative markers. A phylogenetic analysis based on 79 common protein-coding genes indicated that the six examined avocado accessions from three ecological races form a monophyletic clade. The other three genera belonging to the Persea group clustered to form a sister clade with a high bootstrap value. These chloroplast genomes provide important genetic information for future attempts at identifying avocado races and for the related biological research.

Chloroplast genome sequences of Artemisia maritima and Artemisia absinthium: Comparative analyses, mutational hotspots in genus Artemisia and phylogeny in family Asteraceae

Artemisia L. is a complex genus of medicinal importance. Publicly available chloroplast genomes of few Artemisia species are insufficient to resolve taxonomic discrepancies at species level. We report chloroplast genome sequences of two further Artemisia species: A. maritima (151,061 bp) and A. absinthium (151,193 bp). Both genomes possess typical quadripartite structure comprising of a large single copy, a small single copy and a pair of long inverted repeats. The two genomes exhibited high similarities in genome sizes, gene synteny, GC content, synonymous and non-synonymous substitutions, codon usage, amino acids frequencies, RNA editing sites, microsatellites, and oligonucleotide repeats. Transition to transversion ratio was <1. Maximum likelihood tree showed Artemisia a monophyletic genus, sister to genus Chrysanthemum. We also identified 20 highly polymorphic regions including rpoC2-rps2, trnR-UCU-trnG-UCC, rps18-rpl20, and trnL-UAG-rpl32 that could be used to develop authentic and cost-effective markers to resolve taxonomic discrepancies and infer phylogenetic relationships among Artemisia species.

Comparison of Four Complete Chloroplast Genomes of Medicinal and Ornamental Meconopsis Species: Genome Organization and Species Discrimination

High-throughput sequencing of chloroplast genomes has been used to gain insight into the evolutionary relationships of plant species. In this study, we sequenced the complete chloroplast genomes of four species in the Meconopsis genus: M. racemosa, M. integrifolia (Maxim.) Franch, M. horridula and M. punicea. These plants grow in the wild and are recognized as having important medicinal and ornamental applications. The sequencing results showed that the size of the Meconopsis chloroplast genome ranges from 151864 to 153816 bp. A total of 127 genes comprising 90 protein-coding genes, 37 tRNA genes and 8 rRNA genes were observed in all four chloroplast genomes. Comparative analysis of the four chloroplast genomes revealed five hotspot regions (matK, rpoC2, petA, ndhF, and ycf1), which could potentially be used as unique molecular markers for species identification. In addition, the ycf1 gene may also be used as an effective molecular marker to distinguish Papaveraceae and determine the evolutionary relationships among plant species in the Papaveraceae family. Futhermore, these four genomes can provide valuable genetic information for other related studies.

The complete chloroplast genome of candidate new species from Rosa rugosa in Korea (Rosaceae)

Complete chloroplast genome of candidate new species from Rosa rugosa, named as Rosa angusta, is 156,989 bp long and has four subregions: 86,227 bp of large single copy (LSC) and 18,816 bp of small single copy (SSC) regions are separated by 25,793 bp of inverted repeat (IR) regions including 130 genes (85 protein-coding genes, eight rRNAs, and 37 tRNAs). The overall GC content of this chloroplast genome is 37.2% and in the LSC, SSC, and IR regions are 35.2%, 31.1%, and 42.8%, respectively. Phylogenetic trees show that R. angusta is close to R. rugosa with enough number of sequence variations.