Application of Chloroplast Phylogenomics to Resolve Species Relationships Within the Plant Genus Amaranthus

Complete Chloroplast Genome of Alternanthera sessilis and Comparative Analysis with Its Congeneric Invasive Weed Alternanthera philoxeroides

Alternanthera sessilis is considered the closest relative to the invasive weed Alternanthera philoxeroides in China, making it an important native species for studying the invasive mechanisms and adaptations of A. philoxeroides. Chloroplasts play a crucial role in a plant's environmental adaptation, with their genomes being pivotal in the evolution and adaptation of both invasive and related species. However, the chloroplast genome of A. sessilis has remained unknown until now. In this study, we sequenced and assembled the complete chloroplast genome of A. sessilis using high-throughput sequencing. The A. sessilis chloroplast genome is 151,935 base pairs long, comprising two inverted repeat regions, a large single copy region, and a small single copy region. This chloroplast genome contains 128 genes, including 8 rRNA-coding genes, 37 tRNA-coding genes, 4 pseudogenes, and 83 protein-coding genes. When compared to the chloroplast genome of the invasive weed A. philoxeroides and other Amaranthaceae species, we observed significant variations in the ccsA, ycf1, and ycf2 regions in the A. sessilis chloroplast genome. Moreover, two genes, ccsA and accD, were found to be undergoing rapid evolution due to positive selection pressure. The phylogenetic trees were constructed for the Amaranthaceae family, estimating the time of independent species formation between A. philoxeroides and A. sessilis to be approximately 3.5186-8.8242 million years ago. These findings provide a foundation for understanding the population variation within invasive species among the Alternanthera genus.

The Complete Chloroplast Genomes of Blepharoglossum elegans and B. grossum and Comparative Analysis with Related Species (Orchidaceae, Malaxideae)

Blepharoglossum is a rare orchid genus of the Malaxidinae primarily distributed in tropical Pacific islands, with several species occurring in the Taiwan and Hainan Islands of China. Currently, the monophyletic status of Blepharoglossum has been challenged, and the phylogenetic relationships among its allied groups have remained unresolved with traditional DNA markers. In this study, we initially sequenced and annotated the chloroplast (cp) genomes of two Blepharoglossum species, Blepharoglossum elegans (Lindl.) L. Li and Blepharoglossum grossum (Rchb.f.) L. Li. These cp genomes of Blepharoglossum share the typical quadripartite and circular structure. Each of the genomes encodes a total of 133 functional genes, including 87 protein-coding genes (CDS), 38 tRNA genes and 8 rRNA genes. By comparing the sequence differences between these two cp genomes, it was found that they are relatively conserved in terms of overall gene content and gene arrangement. However, a total of 684 SNPs and 2664 indels were still identified, with ycf1, clpP, and trnK-UUU protein-coding genes having the highest number of SNPs and indels. In further comparative analyses among the six cp genomes in Malaxidinae, significant sequence divergences were identified in the intergenic regions, namely rps16-trnQ-UUG, trnS-GCU-trnG-GCC, rpoB-trnC-GCA, trnE-UUC-trnT-GGU, trnF-GAA-trnV-UAC, atpB-rbcL, petA-psbJ, psbE-petL, psbB-psbT, trnN-GUU-rpl32, trnV-GAC-rps7, and rps7-trnL-CAA, and five coding regions, including matK, and rpoC2, ycf1, and two ycf2 genes. Phylogenetic analysis indicated that Blepharoglossum and Oberonia form a highly supported sister group relationship. Our results are consistent with previous studies and present increased resolution among major clades.

Characterization of the Chloroplast Genome of Argyranthemum frutescens and a Comparison with Other Species in Anthemideae

Argyranthemum frutescens, which belongs to the Anthemideae (Asteraceae), is widely cultivated as an ornamental plant. In this study, the complete chloroplast genome of A. frutescens was obtained based on the sequences generated by Illumina HiSeq. The chloroplast genome of A. frutescens was 149,626 base pairs (bp) in length, containing a pair of inverted repeats (IR, 24,510 bp) regions separated by a small single-copy (SSC, 18,352 bp) sequence and a large single-copy (LSC, 82,254 bp) sequence. The genome contained 132 genes, consisting of 85 coding DNA sequences, 37 tRNA genes, and 8 rRNA genes, with nineteen genes duplicated in the IR region. A comparison chloroplast genome analysis among ten species from the tribe of Anthemideae revealed that the chloroplast genome size varied, but the genome structure, gene content, and oligonucleotide repeats were highly conserved. Highly divergent regions, e.g., ycf1, trnK-psbK, petN-psbM intronic, were detected. Phylogenetic analysis supported Argyranthemum as a separate genus. The findings of this study will be helpful in the exploration of the phylogenetic relationships of the tribe of Anthemideae and contribute to the breeding improvement of A. frutescens.

Recovery of chloroplast genomes from medieval millet grains excavated from the Areni-1 cave in southern Armenia

Panicum miliaceum L. was domesticated in northern China at least 7000 years ago and was subsequentially adopted in many areas throughout Eurasia. One such locale is Areni-1 an archaeological cave site in Southern Armenia, where vast quantities archaeobotanical material were well preserved via desiccation. The rich botanical material found at Areni-1 includes P. miliaceum grains that were identified morphologically and¹⁴C dated to the medieval period (873 ± 36 CE and 1118 ± 35 CE). To investigate the demographic and evolutionary history of the Areni-1 millet, we used ancient DNA extraction, hybridization capture enrichment, and high throughput sequencing to assemble three chloroplast genomes from the medieval grains and then compared these sequences to 50 modern P. miliaceum chloroplast genomes. Overall, the chloroplast genomes contained a low amount of diversity with domesticated accessions separated by a maximum of 5 SNPs and little inference on demography could be made. However, in phylogenies the chloroplast genomes separated into two clades, similar to what has been reported for nuclear DNA from P. miliaceum. The chloroplast genomes of two wild (undomesticated) accessions of P. miliaceum contained a relatively large number of variants, 11 SNPs, not found in the domesticated accessions. These results demonstrate that P. miliaceum grains from archaeological sites can preserve DNA for at least 1000 years and serve as a genetic resource to study the domestication of this cereal crop.

Plastome structure of 8 Calanthe s.l. species (Orchidaceae): comparative genomics, phylogenetic analysis

Background

Calanthe (Epidendroideae, Orchidaceae) is a pantropical genus distributed in Asia and Africa. Its species are of great importance in terms of economic, ornamental and medicinal values. However, due to limited and confusing delimitation characters, the taxonomy of the Calanthe alliance (Calanthe, Cephalantheropsis, and Phaius) has not been sufficiently resolved. Additionally, the limited genomic information has shown incongruences in its systematics and phylogeny. In this study, we used illumina platform sequencing, performed a de novo assembly, and did a comparative analysis of 8 Calanthe group species' plastomes: 6 Calanthe and 2 Phaius species. Phylogenetic analyses were used to reconstruct the relationships of the species as well as with other species of the family Orchidaceae.

Results

The complete plastomes of the Calanthe group species have a quadripartite structure with varied sizes ranging between 150,105bp-158,714bp, including a large single-copy region (LSC; 83,364bp- 87,450bp), a small single-copy region (SSC; 16,297bp -18,586bp), and a pair of inverted repeat regions (IRs; 25,222bp - 26,430bp). The overall GC content of these plastomes ranged between 36.6-36.9%. These plastomes encoded 131-134 differential genes, which included 85-88 protein-coding genes, 37-38 tRNA genes, and 8 rRNA genes. Comparative analysis showed no significant variations in terms of their sequences, gene content, gene order, sequence repeats and the GC content hence highly conserved. However, some genes were lost in C. delavayi (P. delavayi), including ndhC, ndhF, and ndhK genes. Compared to the coding regions, the non-coding regions had more sequence repeats hence important for species DNA barcoding. Phylogenetic analysis revealed a paraphyletic relationship in the Calanthe group, and confirmed the position of Phaius delavayi in the genus Calanthe as opposed to its previous placement in Phaius.

Conclusion

This study provides a report on the complete plastomes of 6 Calanthe and 2 Phaius species and elucidates the structural characteristics of the plastomes. It also highlights the power of plastome data to resolve phylogenetic relationships and clarifies taxonomic disputes among closely related species to improve our understanding of their systematics and evolution. Furthermore, it also provides valuable genetic resources and a basis for studying evolutionary relationships and population genetics among orchid species.

Plastome structure of 8 Calanthe s.l. species (Orchidaceae): comparative genomics, phylogenetic analysis

BACKGROUND

Calanthe (Epidendroideae, Orchidaceae) is a pantropical genus distributed in Asia and Africa. Its species are of great importance in terms of economic, ornamental and medicinal values. However, due to limited and confusing delimitation characters, the taxonomy of the Calanthe alliance (Calanthe, Cephalantheropsis, and Phaius) has not been sufficiently resolved. Additionally, the limited genomic information has shown incongruences in its systematics and phylogeny. In this study, we used illumina platform sequencing, performed a de novo assembly, and did a comparative analysis of 8 Calanthe group species' plastomes: 6 Calanthe and 2 Phaius species. Phylogenetic analyses were used to reconstruct the relationships of the species as well as with other species of the family Orchidaceae.

RESULTS

The complete plastomes of the Calanthe group species have a quadripartite structure with varied sizes ranging between 150,105bp-158,714bp, including a large single-copy region (LSC; 83,364bp- 87,450bp), a small single-copy region (SSC; 16,297bp -18,586bp), and a pair of inverted repeat regions (IRs; 25,222bp - 26,430bp). The overall GC content of these plastomes ranged between 36.6-36.9%. These plastomes encoded 131-134 differential genes, which included 85-88 protein-coding genes, 37-38 tRNA genes, and 8 rRNA genes. Comparative analysis showed no significant variations in terms of their sequences, gene content, gene order, sequence repeats and the GC content hence highly conserved. However, some genes were lost in C. delavayi (P. delavayi), including ndhC, ndhF, and ndhK genes. Compared to the coding regions, the non-coding regions had more sequence repeats hence important for species DNA barcoding. Phylogenetic analysis revealed a paraphyletic relationship in the Calanthe group, and confirmed the position of Phaius delavayi in the genus Calanthe as opposed to its previous placement in Phaius.

CONCLUSION

This study provides a report on the complete plastomes of 6 Calanthe and 2 Phaius species and elucidates the structural characteristics of the plastomes. It also highlights the power of plastome data to resolve phylogenetic relationships and clarifies taxonomic disputes among closely related species to improve our understanding of their systematics and evolution. Furthermore, it also provides valuable genetic resources and a basis for studying evolutionary relationships and population genetics among orchid species.

A rapid, simple, and reliable assay to authenticate Peruvian kiwicha (A. caudatus) for food applications

Amaranth has acquired great economic impact as functional food, with species originating from Mexico dominating global trade. In contrast, the Peruvian A. caudatus (kiwicha) has been vastly neglected, although it is endowed with very promising nutritive traits. Morphological plasticity and taxonomic ambiguities render authentication of Amaranth difficult, such that the identity of commercial samples is often unclear. To safeguard the authenticity of kiwicha and, thus, consumer safety, we characterised a germplasm collection of 84 Amaranth accessions on both, the morphological and the genetic level. We show that kiwicha can be delineated phenotypically from other species by its late flowering, taller posture, and lower grain yields. Instead, flower and seed color, often used as proxy for identity, do not qualify as taxonomic markers. Using the plastidic barcoding marker psbA-trnH igs we were able to identify a specific Single Nucleotide Polymorphism (SNP) that separated kiwicha from all other species of Amaranth. This allowed us to develop a sequencing-free authentication assay using an Amplified Refractory Mutation System (ARMS) strategy. As a result kiwicha in commercial samples can be authenticated by a single duplex-PCR yielding a diagnostic side band reporting A. caudatus against all other species of Amaranthus. This fingerprinting assay will help to develop the nutritive potential of kiwicha and to safeguard seed material for A. caudatus against adulteration by the far more prevalent species from Mexico.

Phylogenetic Analysis of Wild Species and the Maternal Origin of Cultivars in the Genus Lilium Using 114 Plastid Genomes

Lilies are one of the most important ornamental flowers worldwide with approximately 100 wild species and numerous cultivars, but the phylogenetic relationships among wild species and their contributions to these cultivars are poorly resolved. We collected the major Lilium species and cultivars and assembled their plastome sequences. Our phylogenetic reconstruction using 114 plastid genomes, including 70 wild species representing all sections and 42 cultivars representing six hybrid divisions and two outgroups, uncovered well-supported genetic relationships within Lilium. The wild species were separated into two distinct groups (groups A and B) associated with geographical distribution, which further diversified into eight different clades that were phylogenetically well supported. Additional support was provided by the distributions of indels and single-nucleotide variants, which were consistent with the topology. The species of sections Archelirion, Sinomartagon III, and Leucolirion 6a and 6b were the maternal donors for Oriental hybrids, Asiatic hybrids, Trumpet hybrids, and Longiflorum hybrids, respectively. The maternal donors of the OT hybrids originated from the two sections Archelirion and Leucolirion 6a, and LA hybrids were derived from the two sections Leucolirion 6b and Sinomartagon. Our study provides an important basis for clarifying the infrageneric classification and the maternal origin of cultivars in Lilium.

Prospects and potentials of underutilized leafy Amaranths as vegetable use for health-promotion

Climate change causes environmental variation worldwide, which is one of the most serious threats to global food security. In addition, more than 2 billion people in the world are reported to suffer from serious malnutrition, referred to as 'hidden hunger.' Dependence on only a few crops could lead to the loss of genetic diversity and high fragility of crop breeding in systems adapting to global scale climate change. The exploitation of underutilized species and genetic resources, referred to as orphan crops, could be a useful approach for resolving the issue of adaptability to environmental alteration, biodiversity preservation, and improvement of nutrient quality and quantity to ensure food security. Moreover, the use of these alternative crops will help to increase the human health benefits and the income of farmers in developing countries. In this review, we highlight the potential of orphan crops, especially amaranths, for use as vegetables and health-promoting nutritional components. This review highlights promising diversified sources of amaranth germplasms, their tolerance to abiotic stresses, and their nutritional, phytochemical, and antioxidant values for vegetable purposes. Betalains (betacyanins and betaxanthins), unique antioxidant components in amaranth vegetables, are also highlighted regarding their chemodiversity across amaranth germplasms and their stability and degradation. In addition, we discuss the physiological functions, antioxidant, antilipidemic, anticancer, and antimicrobial activities, as well as the biosynthesis pathway, molecular, biochemical, genetics, and genomic mechanisms of betalains in detail.

Genetic variation and structure of complete chloroplast genome in alien monoecious and dioecious Amaranthus weeds

Amaranthus is a complex taxon with economic importance as well as harmful weeds. We studied the genetic variation and structure of the chloroplast genomes of 22 samples from 17 species of three subgenera. It was found that the length of the chloroplast genome of Amaranthus varied from 149,949 bp of A. polygonoides to 150,757 bp of A. albus. The frequencies of SNPs and InDels in chloroplast genomes were 1.79% and 2.86%, and the variation mainly occurred in the non-coding regions. The longest InDel was 387 bp, which occurred on ycf2, followed by 384 bp InDel on psbM-trnD. Two InDels in ndhE-I on the SSC make the three subgenera clearly distinguished. In LSC, SSC and IRs regions, there were four 30 bp forward and reverse repeats, and the repeats in SSC and LSC were in nearly opposite positions in circular genome structure, and almost divided the circular genome into symmetrical structures. In the topological tree constructed by chloroplast genome, species in subgen. Amaranthus and subgen. Acnida form monophyletic branches separately and cluster together. A. albus, A. blitoides and A. polygonoides were separated from subgen. Albersia, and the rest of subgen. Albersia were clustered into a monophyletic branch. The rpoC2, ycf1, ndhF-rpl32 were good at distinguishing most amaranths. The trnk-UUU-atpF, trnT-UGU-atpB, psbE-clpP, rpl14-rps19, and ndhF-D can distinguish several similar species. In general, the chloroplast genome is of certain value for the identification of the similar species of Amaranthus, which provides more evidence for clarifying the phylogenetic relationships within the genus.

Comparative plastome analysis of Musaceae and new insights into phylogenetic relationships

Background

Musaceae is an economically important family consisting of 70-80 species. Elucidation of the interspecific relationships of this family is essential for a more efficient conservation and utilization of genetic resources for banana improvement. However, the scarcity of herbarium specimens and quality molecular markers have limited our understanding of the phylogenetic relationships in wild species of Musaceae. Aiming at improving the phylogenetic resolution of Musaceae, we analyzed a comprehensive set of 49 plastomes for 48 species/subspecies representing all three genera of this family.

Results

Musaceae plastomes have a relatively well-conserved genomic size and gene content, with a full length ranging from 166,782 bp to 172,514 bp. Variations in the IR borders were found to show phylogenetic signals to a certain extent in Musa. Codon usage bias analysis showed different preferences for the same codon between species and three genera and a common preference for A/T-ending codons. Among the two genes detected under positive selection (dN/dS > 1), ycf2 was indicated under an intensive positive selection. The divergent hotspot analysis allowed the identification of four regions (ndhF-trnL, ndhF, matK-rps16, and accD) as specific DNA barcodes for Musaceae species.

Bayesian and maximum likelihood phylogenetic analyses using full plastome resulted in nearly identical tree topologies with highly supported relationships between species. The monospecies genus Musella is sister to Ensete, and the genus Musa was divided into two large clades, which corresponded well to the basic number of n = x = 11 and n = x =10/9/7, respectively. Four subclades were divided within the genus Musa. A dating analysis covering the whole Zingiberales indicated that the divergence of Musaceae family originated in the Palaeocene (59.19 Ma), and the genus Musa diverged into two clades in the Eocene (50.70 Ma) and then started to diversify from the late Oligocene (29.92 Ma) to the late Miocene. Two lineages (Rhodochlamys and Australimusa) radiated recently in the Pliocene /Pleistocene periods.

Conclusions

The plastome sequences performed well in resolving the phylogenetic relationships of Musaceae and generated new insights into its evolution. Plastome sequences provided valuable resources for population genetics and phylogenetics at lower taxon.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08454-3.

Characterization of the first complete chloroplast genome of Amaranthus hybridus (Caryophyllales: Amaranthaceae) with phylogenetic implications

In the present study, the complete chloroplast genome of Amaranthus hybridus was sequenced and assembled. The complete chloroplast genome of Amaranthus hybridus is 150,709 in size, with the GC content of 36.56%. The chloroplast genome of Amaranthus hybridus contained 86 protein-coding genes (PCGs), eight ribosomal RNA (rRNA) genes, and 37 transfer RNA (tRNA) genes. Phylogenetic analysis based on combined chloroplast gene dataset indicated that the Amaranthus hybridus exhibited a close relationship with A. hypochondriacus and A. caudatus.

Identification of Amaranthus Species Using Visible-Near-Infrared (Vis-NIR) Spectroscopy and Machine Learning Methods

The feasibility of rapid and non-destructive classification of six different Amaranthus species was investigated using visible-near-infrared (Vis-NIR) spectra coupled with chemometric approaches. The focus of this research would be to use a handheld spectrometer in the field to classify six Amaranthus sp. in different geographical regions of South Korea. Spectra were obtained from the adaxial side of the leaves at 1.5 nm intervals in the Vis-NIR spectral range between 400 and 1075 nm. The obtained spectra were assessed with four different preprocessing methods in order to detect the optimum preprocessing method with high classification accuracy. Preprocessed spectra of six Amaranthus sp. were used as input for the machine learning-based chemometric analysis. All the classification results were validated using cross-validation to produce robust estimates of classification accuracies. The different combinations of preprocessing and modeling were shown to have a classification accuracy of between 71% and 99.7% after the cross-validation. The combination of Savitzky-Golay preprocessing and Support vector machine showed a maximum mean classification accuracy of 99.7% for the discrimination of Amaranthus sp. Considering the high number of spectra involved in this study, the growth stage of the plants, varying measurement locations, and the scanning position of leaves on the plant are all important. We conclude that Vis-NIR spectroscopy, in combination with appropriate preprocessing and machine learning methods, may be used in the field to effectively classify Amaranthus sp. for the effective management of the weedy species and/or for monitoring their food applications.

Genome-wide microsatellites in amaranth: development, characterization, and cross-species transferability

Amaranth (Amaranthus spp.) belonging to Amaranthaceae, is known as "the crop of the future" because of its incredible nutritional quality. Amaranthus spp. (> 70) have a huge diversity in terms of their plant morphology, production and nutritional quality; however, these species are not well characterized at molecular level due to unavailability of robust and reproducible molecular markers, which is essential for crop improvement programs. In the present study, 13,051 genome-wide microsatellite motifs were identified and subsequently utilized for marker development using A. hypochondriacus (L.) genome (JPXE01.1). Out of those, 1538 motifs were found with flanking sequences suitable for primer designing. Among designed primers, 225 were utilized for validation of which 119 (52.89%) primers were amplified. Cross-species transferability and evolutionary relatedness among ten species of Amaranthus (A. hypochondriacus, A. caudatus, A. retroflexus, A. cruentus, A. tricolor, A. lividus, A. hybridus, A. viridis, A. edulis, and A. dubius) were also studied using 45 microsatellite motifs. The maximum (86.67%) and minimum (28.89%) cross-species transferability were observed in A. caudatus and A. dubius, respectively, that indicated high variability present across the Amaranthus spp. Total 97 alleles were detected among 10 species of Amaranthus. The averages of major allele frequency, gene diversity, heterozygosity and PIC were 0.733, 0.347, 0.06, and 0.291, respectively. Nei's genetic dissimilarity coefficients ranged from 0.0625 (between A. tricolor and A. hybridus) to 0.7918 (between A. viridis and A. lividus). The phylogenetic tree grouped ten species into three major clusters. Genome-wide development of microsatellite markers and their transferability revealed relationships among amaranth species which ultimately can be useful for species identification, DNA fingerprinting, and QTLs/gene(s) identification.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02930-5.

Genetics and breeding for climate change in Orphan crops

Climate change is rapidly changing how we live, what we eat and produce, the crops we breed and the target traits. Previously underutilized orphan crops that are climate resilient are receiving much attention from the crops research community, as they are often the only crops left in the field after periods of extreme weather conditions. There are several orphan crops with incredible resilience to biotic and abiotic stresses. Some are nutritious, while others provide good sources of biofuel, medicine and other industrial raw materials. Despite these benefits, orphan crops are still lacking in important genetic and genomic resources that could be used to fast track their improvement and make their production profitable. Progress has been made in generating draft genomes of at least 28 orphan crops over the last decade, thanks to the reducing cost of sequencing. The implementation of a structured breeding program that takes advantage of additional modern crop improvement tools such as genomic selection, speed breeding, genome editing, high throughput phenotyping and breeding digitization would make rapid improvement of these orphan crops possible, but would require coordinated research investment. Other production challenges such as lack of adequate germplasm conservation, poor/non-existent seed systems and agricultural extension services, as well as poor marketing channels will also need to be improved if orphan crops were to be profitable. We review the importance of breeding orphan crops under the increasing effects of climate change, highlight existing gaps that need to be addressed and share some lessons to be learned from major crops.

Chloroplast genome variation and phylogenetic relationships of Atractylodes species

Background

Atractylodes DC is the basic original plant of the widely used herbal medicines “Baizhu” and “Cangzhu” and an endemic genus in East Asia. Species within the genus have minor morphological differences, and the universal DNA barcodes cannot clearly distinguish the systemic relationship or identify the species of the genus. In order to solve these question, we sequenced the chloroplast genomes of all species of Atractylodes using high-throughput sequencing.

Results

The results indicate that the chloroplast genome of Atractylodes has a typical quadripartite structure and ranges from 152,294 bp (A. carlinoides) to 153,261 bp (A. macrocephala) in size. The genome of all species contains 113 genes, including 79 protein-coding genes, 30 transfer RNA genes and four ribosomal RNA genes. Four hotspots, rpl22-rps19-rpl2, psbM-trnD, trnR-trnT^(GGU), and trnT^(UGU)-trnL, and a total of 42–47 simple sequence repeats (SSR) were identified as the most promising potentially variable makers for species delimitation and population genetic studies. Phylogenetic analyses of the whole chloroplast genomes indicate that Atractylodes is a clade within the tribe Cynareae; Atractylodes species form a monophyly that clearly reflects the relationship within the genus.

Conclusions

Our study included investigations of the sequences and structural genomic variations, phylogenetics and mutation dynamics of Atractylodes chloroplast genomes and will facilitate future studies in population genetics, taxonomy and species identification.

Supplementary Information

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

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.

Sequencing and Analysis of Chrysanthemum carinatum Schousb and Kalimeris indica. The Complete Chloroplast Genomes Reveal Two Inversions and rbcL as Barcoding of the Vegetable

Chrysanthemum carinatum Schousb and Kalimeris indica are widely distributed edible vegetables and the sources of the Chinese medicine Asteraceae. The complete chloroplast (cp) genome of Asteraceae usually occurs in the inversions of two regions. Hence, the cp genome sequences and structures of Asteraceae species are crucial for the cp genome genetic diversity and evolutionary studies. Hence, in this paper, we have sequenced and analyzed for the first time the cp genome size of C. carinatum Schousb and K. indica, which are 149,752 bp and 152,885 bp, with a pair of inverted repeats (IRs) (24,523 bp and 25,003) separated by a large single copy (LSC) region (82,290 bp and 84,610) and a small single copy (SSC) region (18,416 bp and 18,269), respectively. In total, 79 protein-coding genes, 30 distinct transfer RNA (tRNA) genes, four distinct rRNA genes and two pseudogenes were found not only in C. carinatum Schousb but also in the K. indica cp genome. Fifty-two (52) and fifty-nine (59) repeats, and seventy (70) and ninety (90) simple sequence repeats (SSRs) were found in the C. carinatum Schousb and K. indica cp genomes, respectively. Codon usage analysis showed that leucine, isoleucine, and serine are the most frequent amino acids and that the UAA stop codon was the significantly favorite stop codon in both cp genomes. The two inversions, the LSC region ranging from trnC-GCA to trnG-UCC and the whole SSC region were found in both of them. The complete cp genome comparison with other Asteraceae species showed that the coding area is more conservative than the non-coding area. The phylogenetic analysis revealed that the rbcL gene is a good barcoding marker for identifying different vegetables. These results give an insight into the identification, the barcoding, and the understanding of the evolutionary model of the Asteraceae cp genome.