The complete plastome sequence of Rubus takesimensis endemic to Ulleung Island, Korea: Insights into molecular evolution of anagenetically derived species in Rubus (Rosaceae)

Characteristics and phylogenetic analysis of the complete chloroplast genome of Rubus chingii Hu 1925 from the family Rosaceae

Rubus chingii Hu 1925 is an important medicinal vine shrub in the Rosaceae family, widely distributed in China and Japan. In this study, the complete chloroplast genome of R. chingii was sequenced and identified. The chloroplast genome was 155,563 bp in size with a total GC content of 37.06%. Two 25,749-bp inverted repeat (IRA and IRB) regions divided the genome as four sections, with the remainder forming a large single-copy (LSC, 85,322 bp) and a small single-copy (SSC, 18,743 bp) regions. This genome contained a total of 131 genes, of which 86 were protein-coding genes, 37 tRNA genes, and eight rRNA genes. The phylogenetic analysis showed that R. chingii, along with several other R. longisepalus, R. tsangii, R. hirsutus, R. taiwanicola, R. rubroangustifolius, and R. glandulosopunctatus, formed the monophylic group. Interestingly, the chloroplast genome structure we reported was different from the previously reported structure and provided richer phylogenetic analysis information in the Rubus genus compared to previous studies. The genome information reported in this paper will provide some useful information for further investigation on the evolution of the family Rosaceae.

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.

Characterization of the complete chloroplast genome sequence of Potentilla gageodoensis (rosaceae), endemic to the continental islands of Korea

This study reports the complete chloroplast genome sequence of a continental island endemic, Potentilla gageodoensis. The total plastome size was 156,397 bp, comprising one large single-copy (LSC; 85,768 bp), one small single-copy (SSC; 18,589 bp), and two inverted repeat (IR) regions (IR_a and IR_b, each with 26,020 bp). The overall GC content was 36.92%, and the plastome contained 131 genes, comprising 84 protein-coding genes with two pseudogenes (infA and ycf1), 37 transfer RNA genes, and eight ribosomal RNA genes. Phylogenetic analysis performed using 27 representative Rosoideae plastomes suggests that the genus Potentilla is not monophyletic and that P. gageodoensis is sister to the clade containing four taxa of Potentilla (P. freyniana, P. freyniana var. chejuensis, P. stolonifera, and P. stolonifera var. quelpaertensis). The present study reveals the taxonomic distinction of P. gageodoensis from its congeneric species in Korea and the plastome sequence obtained from this study can be used to study phylogenetic relationships and taxonomic status.

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.

Plastome Diversity and Phylogenomic Relationships in Asteraceae

Plastid genomes are in general highly conserved given their slow evolutionary rate, and thus large changes in their structure are unusual. However, when specific rearrangements are present, they are often phylogenetically informative. Asteraceae is a highly diverse family whose evolution is long driven by polyploidy (up to 48x) and hybridization, both processes usually complicating systematic inferences. In this study, we generated one of the most comprehensive plastome-based phylogenies of family Asteraceae, providing information about the structure, genetic diversity and repeat composition of these sequences. By comparing the whole-plastome sequences obtained, we confirmed the double inversion located in the long single-copy region, for most of the species analyzed (with the exception of basal tribes), a well-known feature for Asteraceae plastomes. We also showed that genome size, gene order and gene content are highly conserved along the family. However, species representative of the basal subfamily Barnadesioideae-as well as in the sister family Calyceraceae-lack the pseudogene rps19 located in one inverted repeat. The phylogenomic analysis conducted here, based on 63 protein-coding genes, 30 transfer RNA genes and 21 ribosomal RNA genes from 36 species of Asteraceae, were overall consistent with the general consensus for the family's phylogeny while resolving the position of tribe Senecioneae and revealing some incongruences at tribe level between reconstructions based on nuclear and plastid DNA data.

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.

The complete chloroplast genome of Castanopsis sieboldii (Makino) Hatus (Fagaceae)

Castanopsis sieboldii (Makino) Hatus is an evergreen tree that distributes in Eastern Asia including Islands of Korea and Japan. The chloroplast genome of C. sieboldii was successfully sequenced. Its length is 160,705 bp long (GC ratio is 36.8%) and has four subregions: 90,821 bp of large single copy (34.6%) and 19,014 bp of small single copy (30.8%) regions are separated by 25,075 bp of inverted repeat (42.8%) regions including 134 genes (89 protein-coding genes, eight rRNAs, and 37 tRNAs). Interspecific variations of Castanopsis are at a moderate level in comparison to those of the other genera. Phylogenetic trees show that C. sieboldii chloroplast genome was clustered with the other two Castanopsis species.

Comparative analysis of chloroplast genome structure and molecular dating in Myrtales

BACKGROUND

Myrtales is a species rich branch of Rosidae, with many species having important economic, medicinal, and ornamental value. At present, although there are reports on the chloroplast structure of Myrtales, a comprehensive analysis of the chloroplast structure of Myrtales is lacking. Phylogenetic and divergence time estimates of Myrtales are mostly constructed by using chloroplast gene fragments, and the support for relationships is low. A more reliable method to reconstruct the species divergence time and phylogenetic relationships is by using whole chloroplast genomes. In this study, we comprehensively analyzed the structural characteristics of Myrtales chloroplasts, compared variation hotspots, and reconstructed the species differentiation time of Myrtales with four fossils and one secondary calibration point.

RESULTS

A total of 92 chloroplast sequences of Myrtales, representing six families, 16 subfamilies and 78 genera, were obtained including nine newly sequenced chloroplasts by whole genome sequencing. Structural analyses showed that the chloroplasts range in size between 152,214-171,315 bp and exhibit a typical four part structure. The IR region is between 23,901-36,747 bp, with the large single copy region spanning 83,691-91,249 bp and the small single copy region spanning 11,150-19,703 bp. In total, 123-133 genes are present in the chloroplasts including 77-81 protein coding genes, four rRNA genes and 30-31 tRNA genes. The GC content was 36.9-38.9%, with the average GC content being 37%. The GC content in the LSC, SSC and IR regions was 34.7-37.3%, 30.6-36.8% and 39.7-43.5%, respectively. By analyzing nucleotide polymorphism of the chloroplast, we propose 21 hypervariable regions as potential DNA barcode regions for Myrtales. Phylogenetic analyses showed that Myrtales and its corresponding families are monophyletic, with Combretaceae and the clade of Onagraceae + Lythraceae (BS = 100%, PP = 1) being sister groups. The results of molecular dating showed that the crown of Myrtales was most likely to be 104.90 Ma (95% HPD = 87.88-114.18 Ma), and differentiated from the Geraniales around 111.59 Ma (95% HPD = 95.50-118.62 Ma).

CONCLUSIONS

The chloroplast genome structure of Myrtales is similar to other angiosperms and has a typical four part structure. Due to the expansion and contraction of the IR region, the chloroplast genome sizes in this group are slightly different. The variation of noncoding regions of the chloroplast genome is larger than those of coding regions. Phylogenetic analysis showed that Combretaceae and Onagraceae + Lythraceae were well supported as sister groups. Molecular dating indicates that the Myrtales crown most likely originated during the Albian age of the Lower Cretaceous. These chloroplast genomes contribute to the study of genetic diversity and species evolution of Myrtales, while providing useful information for taxonomic and phylogenetic studies of Myrtales.

Structural Mutations in the Organellar Genomes of Valeriana sambucifolia f. dageletiana (Nakai. ex Maekawa) Hara Show Dynamic Gene Transfer

Valeriana sambucifolia f. dageletiana (Nakai. ex Maekawa) Hara is a broad-leaved valerian endemic to Ulleung Island, a noted hot spot of endemism in Korea. However, despite its widespread pharmacological use, this plant remains comparatively understudied. Plant cells generally contain two types of organellar genomes (the plastome and the mitogenome) that have undergone independent evolution, which accordingly can provide valuable information for elucidating the phylogenetic relationships and evolutionary histories of terrestrial plants. Moreover, the extensive mega-data available for plant genomes, particularly those of plastomes, can enable researchers to gain an in-depth understanding of the transfer of genes between different types of genomes. In this study, we analyzed two organellar genomes (the 155,179 bp plastome and the 1,187,459 bp mitogenome) of V. sambucifolia f. dageletiana and detected extensive changes throughout the plastome sequence, including rapid structural mutations associated with inverted repeat (IR) contraction and genetic variation. We also described features characterizing the first reported mitogenome sequence obtained for a plant in the order Dipsacales and confirmed frequent gene transfer in this mitogenome. We identified eight non-plastome-originated regions (NPRs) distributed within the plastome of this endemic plant, for six of which there were no corresponding sequences in the current nucleotide sequence databases. Indeed, one of these unidentified NPRs unexpectedly showed certain similarities to sequences from bony fish. Although this is ostensibly difficult to explain, we suggest that this surprising association may conceivably reflect the occurrence of gene transfer from a bony fish to the plastome of an ancestor of V. sambucifolia f. dageletiana mediated by either fungi or bacteria.

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.

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.

Evolutionary and phylogenetic aspects of the chloroplast genome of Chaenomeles species

Chaenomeles (family Rosaceae) is a genus of five diploid species of deciduous spiny shrubs that are native to Central Asia and Japan. It is an important horticultural crop (commonly known as flowering quinces) in Europe and Asia for its high yield in fruits that are rich in juice, aroma, and dietary fiber. Therefore, the development of effective genetic markers of Chaenomeles species is advantageous for crop improvement through breeding and selection. In this study, we successfully assembled and analyzed the chloroplast genome of five Chaenomeles species. The chloroplast genomes of the five Chaenomeles species were very similar with no structural or content rearrangements among them. The chloroplast genomes ranged from 159,436 to 160,040 bp in length and contained a total of 112 unique genes, including 78 protein-coding genes, 30 tRNAs, and 4 rRNAs. Three highly variable regions, including trnR-atpA, trnL-F, and rpl32-ccsA, were identified. Phylogenetic analysis based on the complete chloroplast genome showed that Chaenomeles forms a monophyletic clade and had a close relationship with the genera Docynia and Malus. Analyses for phylogenetic relationships and the development of available genetic markers in future could provide valuable information regarding genetics and breeding mechanisms of the Chaenomeles 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.

Comparison of Whole Plastome Sequences between Thermogenic Skunk Cabbage Symplocarpus renifolius and Nonthermogenic S. nipponicus (Orontioideae; Araceae) in East Asia

Symplocarpus, a skunk cabbage genus, includes two sister groups, which are drastically different in life history traits and thermogenesis, as follows: The nonthermogenic summer flowering S. nipponicus and thermogenic early spring flowering S. renifolius. Although the molecular basis of thermogenesis and complete chloroplast genome (plastome) of thermogenic S. renifolius have been well characterized, very little is known for that of S. nipponicus. We sequenced the complete plastomes of S. nipponicus sampled from Japan and Korea and compared them with that of S. renifolius sampled from Korea. The nonthermogenic S. nipponicus plastomes from Japan and Korea had 158,322 and 158,508 base pairs, respectively, which were slightly shorter than the thermogenic plastome of S. renifolius. No structural or content rearrangements between the species pairs were found. Six highly variable noncoding regions (psbC/trnS, petA/psbJ, trnS/trnG, trnC/petN, ycf4/cemA, and rpl3/rpl22) were identified between S. nipponicus and S. renifolius and 14 hot-spot regions were also identified at the subfamily level. We found a similar total number of SSR (simple sequence repeat) motifs in two accessions of S. nipponicus sampled from Japan and Korea. Phylogenetic analysis supported the basal position of subfamily Orontioideae and the monophyly of genus Symplocarpus, and also revealed an unexpected evolutionary relationship between S. nipponicus and S. renifolius.

Rubus ulmifolius Schott fruits: A detailed study of its nutritional, chemical and bioactive properties

There is a growing interest in wild edible species that represent a source of several health-promoting compounds, providing a potential strategy to diversify and enrich the daily diet. The aim of the present work was to characterize the nutritional and chemical composition of Rubus ulmifolius Schott fruits. Furthermore, their antimicrobial activity, non-anthocyanin and anthocyanin phenolic profile were also determined. According to the obtained results, R. ulmifolius fruits exhibited a high concentration in carbohydrates and a low fat content, in comparison with the other nutrients and non-nutrients detected in this sample. The colour parameters demonstrated differences in a* and b* parameters after lyophilisation process. Glucose and fructose were the most abundant free sugars detected and quinic acid showed the highest content compared to the other five organic acids identified. The fatty acids profile revealed 25 compounds, being mostly represented by polyunsaturated fatty acids and evidencing linolenic and α-linolenic acid as the most abundant. All tocopherol isoforms were detected, revealing γ-tocopherol with highest amount. Cyanidin-3-O-glucoside, ellagic acid pentoside, ellagic acid glucuronide and sanguiin H-10 were the main phenolic compounds present, which could be related to the antimicrobial activity (MIC values ranging between 5 and 20 mg/mL) revealed by R. ulmifolius fruits. These results showed that this fruit is a good source of nutrients as also non-nutrient compounds, with human health benefits.