Complete plastid genome sequences of three Rosids (Castanea, Prunus, Theobroma): evidence for at least two independent transfers of rpl22 to the nucleus

Complete chloroplast genomes of the hemiparasitic genus Cymbaria: Insights into comparative analysis, development of molecular markers, and phylogenetic relationships

The hemiparasitic tribe Cymbarieae (Orobanchaceae) plays a crucial role in elucidating the initial stage of the transition from autotrophism to heterotrophism. However, the complete chloroplast genome of the type genus Cymbaria has yet to be reported. In addition, the traditional Mongolian medicine Cymbaria daurica is frequently subjected to adulteration or substitution because of the minor morphological differences with Cymbaria mongolica. In this study, the complete chloroplast genomes of the two Cymbaria species were assembled and annotated, and those of other published 52 Orobanchaceae species were retrieved for comparative analyses. We found that the Cymbaria chloroplast genomes are characterized by pseudogenization or loss of stress-relevant genes (ndh) and a unique rbcL-matK inversion. Unlike the high variability observed in holoparasites, Cymbaria and other hemiparasites exhibit high similarity to autotrophs in genome size, guanine-cytosine (GC) content, and intact genes. Notably, four pairs of specific DNA barcodes were developed and validated to distinguish the medicinal herb from its adulterants. Phylogenetic analyses revealed that the genus Cymbaria and the Schwalbea-Siphonostegia clade are grouped into the tribe Cymbarieae, which forms a sister clade to the remaining Orobanchaceae parasitic lineages. Moreover, the diversification of monophyletic Cymbaria occurred during the late Miocene (6.72 Mya) in the Mongol-Chinese steppe region. Our findings provide valuable genetic resources for studying the phylogeny of Orobanchaceae and plant parasitism, and genetic tools to validate the authenticity of the traditional Mongolian medicine "Xinba.".

Comparative chloroplast genomes study of five officinal Ardisia Species: Unraveling interspecific diversity and evolutionary insights in Ardisia

Ardisia S.W. (Primulaceae), naturally distributed in tropical and subtropical regions, has edible and medicinal values and is prevalent in clinical and daily use in China. More genetic information for distinct species delineation is needed to support the development and utilization of the genus Ardisia. We sequenced, annotated, and compared the chloroplast genomes of five Ardisia species: A. brunnescens, A. pusilla, A. squamulosa, A. crenata, and A. brevicaulis in this study. We found a typical quadripartite structure in all five chloroplast genomes, with lengths ranging from 155,045 to 156,943 bp. Except for A. pusilla, which lacked the ycf15 gene, the other four Ardisia species contained 114 unique genes, including 79 protein-coding genes, 30 tRNAs, and four rRNAs. In addition, the rps19 pseudogene gene was present only in A. brunnescens. Five highly variable DNA barcodes were identified for five Ardisia species, including trnT-GGU-psbD, trnT-UGU-trnL-UAA, rps4-trnT-UGU, rpl32-trnL-UAG, and rpoB-trnC-GAA. The RNA editiing sites of protein-coding genes in the five Ardisia plastome were characterized and compared, and 274 (A. crenata)-288 (A. brevicaulis) were found. The results of the phylogenetic analysis were consistent with the morphological classification. Sequence alignment and phylogenetic analysis showed that ycf15 genes were highly divergent in Primulaceae. Reconstructions of ancestral character states indicated that leaf margin morphology is critical for classifying the genus Ardisia, with a rodent-like character being the most primitive. These results provide valuable information on the taxonomy and evolution of Ardisia plants.

Chloroplast Genomes Evolution and Phylogenetic Relationships of Caragana species

Caragana sensu lato (s.l.) includes approximately 100 species that are mainly distributed in arid and semi-arid regions. Caragana species are ecologically valuable for their roles in windbreaking and sand fixation. However, the taxonomy and phylogenetic relationships of the genus Caragana are still unclear. In this study, we sequenced and assembled the chloroplast genomes of representative species of Caragana and reconstructed robust phylogenetic relationships at the section level. The Caragana chloroplast genome has lost the inverted repeat region and wascategorized in the inverted repeat loss clade (IRLC). The chloroplast genomes of the eight species ranged from 128,458 bp to 135,401 bp and contained 110 unique genes. All the Caragana chloroplast genomes have a highly conserved structure and gene order. The number of long repeats and simple sequence repeats (SSRs) showed significant variation among the eight species, indicating heterogeneous evolution in Caragana. Selective pressure analysis of the genes revealed that most of the protein-coding genes evolved under purifying selection. The phylogenetic analyses indicated that each section forms a clade, except the section Spinosae, which was divided into two clades. This study elucidated the evolution of the chloroplast genome within the widely distributed genus Caragana. The detailed information obtained from this study can serve as a valuable resource for understanding the molecular dynamics and phylogenetic relationships within Caragana.

The complete chloroplast genome sequence of Prunus salicina cultivar 'Zuili' (Rosaceae)

'Zuili' is a distinguished plum (Prunus salicina Lindley 1830) originating from Tongxiang City, Zhejiang Province, China, and a nationally recognized geographical indication product. In this study, we reported the complete chloroplast genome of P. salicina cultivar 'Zuili'. The genome has a circular structure of 157,935 bp containing a large single-copy (LSC) region of 86,133 bp, a small single-copy (SSC) region of 19,028 bp, and two inverted repeats (IRs) of 26,387 by each. It harbors 130 genes (111 unique genes), including 85 protein-coding genes (78 are unique), eight ribosomal RNA genes (four are unique), and 37 transfer genes (29 are unique). The phylogenetic analysis based on whole chloroplast genomes showed 'Zuili' was clustered with Prunus salicina cultivar 'Wuyuecui' (MW406461.1) and 'No. 2 Guofeng' (MW406472.1). This study provides valuable information that can contribute to the identification and further evolutionary analysis of Prunus salicina cultivar 'Zuili'.

Plastid phylogenomics of the tribe potentilleae (Rosaceae)

The tribe Potentilleae comprises approximately 1700 species in 13 genera, making it one of the largest of the 16 tribes in Rosaceae. Our understanding of the composition and relationships among members of Potentilleae has advanced dramatically with the application of molecular markers in the last two decades. Yet there is still much work remaining toward a robust phylogenetic framework for the entire Potentilleae and a comprehensive genus-level dating framework for the tribe. The goals of the present study were to establish a phylogenetic framework for Potentilleae, infer the origin and diversification of the tribe using a temporal framework, and explore the taxonomic implications in light of the updated phylogenetic framework. We used the plastome sequences from 158 accessions representing 139 taxa covering all 13 recognized genera of the tribe to reconstruct the Potentilleae phylogeny. High phylogenetic resolution was recovered along the Potentilleae backbone. Two major clades were recovered within Potentilleae, corresponding to the two subtribes Fragariinae and Potentillinae. Within Fragariinae, two subclades were recovered. In one subclade, Sibbaldia sensu stricto is sister to a clade containing Sibbaldianthe, Comarum, Farinopsis, and Alchemilla sensu lato. In the other subclade, Fragaria is sister to a clade comprising Chamaerhodos, Chamaecallis, Drymocallis, Dasiphora, and Potaninia. Within Potentillinae, Argentina is sister to Potentilla sensu stricto. Within Potentilla sensu stricto, clade Himalaya is sister to Alba, and the Himalaya-Alba clade together is sister to a clade comprising Reptans, Potentilla ancistrifolia Bunge, Fragarioides, Ivesioid, and Argentea. Divergence time estimates indicated that tribe Potentilleae originated during the middle Eocene, and subtribes Fragariinae and Potentillinae diverged around the Eocene-Oligocene transition, and divergence times dated for Potentilleae genera ranged from the early Miocene to the late Pleistocene.

Complete Chloroplast Genome of Hypericum perforatum and Dynamic Evolution in Hypericum (Hypericaceae)

Hypericum perforatum (St. John's Wort) is a medicinal plant from the Hypericaceae family. Here, we sequenced the whole chloroplast genome of H. perforatum and compared the genome variation among five Hypericum species to discover dynamic changes and elucidate the mechanisms that lead to genome rearrangements in the Hypericum chloroplast genomes. The H. perforatum chloroplast genome is 139,725 bp, exhibiting a circular quadripartite structure with two copies of inverted repeats (IRs) separating a large single-copy region and a small single-copy region. The H. perforatum chloroplast genome encodes 106 unique genes, including 73 protein-coding genes, 29 tRNAs, and 4 rRNAs. Hypericum chloroplast genomes exhibit genome rearrangement and significant variations among species. The genome size variation among the five Hypericum species was remarkably associated with the expansion or contraction of IR regions and gene losses. Three genes-trnK-UUU, infA, and rps16-were lost, and three genes-rps7, rpl23, and rpl32-were pseudogenized in Hypericum. All the Hypericum chloroplast genomes lost the two introns in clpP, the intron in rps12, and the second intron in ycf3. Hypericum chloroplast genomes contain many long repeat sequences, suggesting a role in facilitating rearrangements. Most genes, according to molecular evolution assessments, are under purifying selection.

Sequencing and Phylogenetic Analysis of the Chloroplast Genome of Three Apricot Species

The production and quality of apricots in China is currently limited by the availability of germplasm resource characterizations, including identification at the species and cultivar level. To help address this issue, the complete chloroplast genomes of Prunus armeniaca L., P. sibirica L. and kernel consumption apricot were sequenced, characterized, and phylogenetically analyzed. The three chloroplast (cp) genomes ranged from 157,951 to 158,224 bp, and 131 genes were identified, including 86 protein-coding genes, 37 rRNAs, and 8 tRNAs. The GC content ranged from 36.70% to 36.75%. Of the 170 repetitive sequences detected, 42 were shared by all three species, and 53-57 simple sequence repeats were detected with AT base preferences. Comparative genomic analysis revealed high similarity in overall structure and gene content as well as seven variation hotspot regions, including psbA-trnK-UUU, rpoC1-rpoB, rpl32-trnL-UAG, trnK-rps16, ndhG-ndhI, ccsA-ndhD, and ndhF-trnL. Phylogenetic analysis showed that the three apricot species clustered into one group, and the genetic relationship between P. armeniaca and kernel consumption apricot was the closest. The results of this study provide a theoretical basis for further research on the genetic diversity of apricots and the development and utilization of molecular markers for the genetic engineering and breeding of apricots.

Comparative and phylogenetic analysis of complete chloroplast genomes from seven Neocinnamomum taxa (Lauraceae)

The genus Neocinnamomum is considered to be one of the most enigmatic groups in Lauraceae, mainly distributed in tropical and subtropical regions of Southeast Asia. The genus contains valuable oilseed and medicinal tree species. However, there are few studies on the genus Neocinnamomum at present, and its interspecific relationship is still unclear. In order to explore the genetic structure and evolutionary characteristics of the Neocinnamomum chloroplast genome and to resolve the species relationships within the genus, comparative genomic and phylogenetic analyses were performed on the whole chloroplast genome sequences of 51 samples representing seven Neocinnamomum taxa. The whole Neocinnamomum chloroplast genome size ranged from 150,753-150,956 bp, with a GC content of 38.8%-38.9%. A total of 128 genes were annotated within the Neocinnamomum chloroplast genome, including 84 protein coding genes, 8 rRNA genes, and 36 tRNA genes. Between 71-82 SSRs were detected, among which A/T base repeats were the most common. The chloroplast genome contained a total of 31 preferred codons. Three highly variable regions, trnN-GUU-ndhF, petA-psbJ, and ccsA-ndhD, were identified with Pi values > 0.004. Based on the whole chloroplast genome phylogenetic tree, the phylogenetic relationships among the seven Neocinnamomum taxa were determined. N. delavayi and N. fargesii were the most closely related species, and N. lecomtei was identified as the most basal taxon. In this study, the characteristics and sequence variation of the chloroplast genomes of seven Neocinnamomum taxa were revealed, and the genetic relationship among the species was clarified. The results of this study will provide a reference for subsequent molecular marker development and phylogenetic research of Neocinnamomum.

Insights into cryptic speciation of quillworts in China

Cryptic species are commonly misidentified because of high morphological similarities to other species. One group of plants that may harbor large numbers of cryptic species is the quillworts (Isoëtes spp.), an ancient aquatic plant lineage. Although over 350 species of Isoëtes have been reported globally, only ten species have been recorded in China. The aim of this study is to better understand Isoëtes species diversity in China. For this purpose, we systematically explored the phylogeny and evolution of Isoëtes using complete chloroplast genome (plastome) data, spore morphology, chromosome number, genetic structure, and haplotypes of almost all Chinese Isoëtes populations. We identified three ploidy levels of Isoëtes in China-diploid (2n = 22), tetraploid (2n = 44), and hexaploid (2n = 66). We also found four megaspore and microspore ornamentation types in diploids, six in tetraploids, and three in hexaploids. Phylogenetic analyses confirmed that I. hypsophila as the ancestral group of the genus and revealed that Isoëtes diploids, tetraploids, and hexaploids do not form monophyletic clades. Most individual species possess a single genetic structure; however, several samples have conflicting positions on the phylogenetic tree based on SNPs and the tree based on plastome data. All 36 samples shared 22 haplotypes. Divergence time analysis showed that I. hypsophila diverged in the early Eocene (∼48.05 Ma), and most other Isoëtes species diverged 3-20 Ma. Additionally, different species of Isoëtes were found to inhabit different water systems and environments along the Yangtze River. These findings provide new insights into the relationships among Isoëtes species in China, where highly similar morphologic populations may harbor many cryptic species.

The complete chloroplast genome sequence of Castanopsis fordii Hance (Fagaceae)

Castanopsis fordii Hance 1884 is a typical evergreen broad-leaved forest plant in the south subtropical and middle subtropical regions of China. It has high utilization value in wood production and soil erosion protection. Here, we first reported and characterized the complete chloroplast (cp) genome sequence of C. fordii based on Illumina paired-end sequencing data. The complete cp genome sequence of C. fordii was 160,853 base pairs (bp) in length which contained two inverted repeats (IRs) of 25,699 bp separated by a large single-copy (LSC) and a small single copy (SSC) of 90,474 bp and 18,981 bp, respectively. The cpDNA contained 129 genes, comprising 85 protein-coding genes, 36 tRNA genes, 8 rRNA genes. The overall GC content of the plastome was 36.8%. Phylogenetic analysis base on 14 chloroplast genomes indicated that C. fordii was closely related to the species C. tibetana and C. concinna in Fagaceae.

Complete chloroplast genome of Isoetes orientalis (Isoetaceae), an endangered quillwort from China

Isoetes orientalis is an endangered hexaploidy species of Isoetaceae in China and the complete chloroplast genome of this species has not been reported. In the present study, a complete chloroplast genome of Isoetes orientalis (Isoetaceae) was assembled and annotated. This chloroplast genome has a circular structure of 145,504 bp in length, comprising a pair of inverted repeat (IR) regions of 13,207 bp each, a large single-copy (LSC) region of 91,864 bp, and a small single-copy (SSC) region of 27,226 bp. The chloroplast genome contains 136 genes, including 84 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Phylogenetic analysis showed that I. orientalis was closely related to I. sinensis. These results provide additional resources for future studies on Isoetes from China and across the globe.

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.

"Omics" insights into plastid behavior toward improved carotenoid accumulation

Plastids are a group of diverse organelles with conserved carotenoids synthesizing and sequestering functions in plants. They optimize the carotenoid composition and content in response to developmental transitions and environmental stimuli. In this review, we describe the turbulence and reforming of transcripts, proteins, and metabolic pathways for carotenoid metabolism and storage in various plastid types upon organogenesis and external influences, which have been studied using approaches including genomics, transcriptomics, proteomics, and metabonomics. Meanwhile, the coordination of plastid signaling and carotenoid metabolism including the effects of disturbed carotenoid biosynthesis on plastid morphology and function are also discussed. The "omics" insight extends our understanding of the interaction between plastids and carotenoids and provides significant implications for designing strategies for carotenoid-biofortified crops.

Characterization, comparative phylogenetic, and gene transfer analyses of organelle genomes of Rhododendron × pulchrum

Rhododendron × pulchrum, an important horticultural species, is widely distributed in Europe, Asia, and North America. To analyze the phylogenetic and organelle genome information of R. × pulchrum and its related species, the organelle genome of R. × pulchrum was sequenced and assembled. The complete mitochondrial genome showed lineage DNA molecules, which were 816,410 bp long and contained 64 genes, namely 24 transfer RNA (tRNA) genes, 3 ribosomal RNA (rRNA) genes, and 37 protein-coding genes. The chloroplast genome of R. × pulchrum was reassembled and re-annotated; the results were different from those of previous studies. There were 42 and 46 simple sequence repeats (SSR) identified from the mitochondrial and chloroplast genomes of R. × pulchrum, respectively. Five genes (nad1, nad2, nad4, nad7, and rps3) were potentially useful molecular markers. The R. × pulchrum mitochondrial genome collinear alignment among five species of the Ericaceae showed that the mitochondrial genomes of these related species have a high degree of homology with R. × pulchrum in this gene region, and the most conservative genes were trnC-GCA, trnD-GUC, trnM-CAU, trnN-GUU, trnY-GUA, atp4, nad4, nad2, nad5, ccmC, and rrn26. The phylogenetic trees of mitochondrial genome showed that R. simsii was a sister to R. × pulchrum. The results verified that there was gene rearrangement between R. × pulchrum and R. simsii mitochondrial genomes. The codon usage bias of 10 Ericaceae mitochondrial genes and 7 Rhododendron chloroplast genes were influenced by mutation, while other genes codon usages had undergone selection. The study identified 13 homologous fragments containing gene sequences between the chloroplast and mitochondrial genomes of R. × pulchrum. Overall, our results illustrate the organelle genome information could explain the phylogenetics of plants and could be used to develop molecular markers and genetic evolution. Our study will facilitate the study of population genetics and evolution in Rhododendron and other genera in Ericaceae.

The analysis of genetic structure and characteristics of the chloroplast genome in different Japanese apricot germplasm populations

BACKGROUND

Chloroplast (cp) genomes are generally considered to be conservative and play an important role in population diversity analysis in plants, but the characteristics and diversity of the different germplasm populations in Japanese apricot are still not clear.

RESULTS

A total of 146 cp genomes from three groups of wild, domesticated, and bred accessions of Japanese apricot were sequenced in this study. The comparative genome analysis revealed that the 146 cp genomes were divided into 41 types, and ranged in size from 157,886 to 158,167 bp with a similar structure and composition to those of the genus Prunus. However, there were still minor differences in the cp genome that were mainly caused by the contraction and expansion of the IR region, and six types of SSR in which mono-nucleotide repeats were the most dominant type of repeats in the cp genome. The genes rpl33 and psbI, and intergenic regions of start-psbA, rps3-rpl22, and ccsA-ndhD, showed the highest nucleotide polymorphism in the whole cp genome. A total of 325 SNPs were detected in the 146 cp genomes, and more than 70% of the SNPs were in region of large single-copy (LSC). The SNPs and haplotypes in the cp genome indicated that the wild group had higher genetic diversity than the domesticated and bred groups. In addition, among wild populations, Southwest China, including Yunnan, Tibet, and Bijie of Guizhou, had the highest genetic diversity. The genetic relationship of Japanese apricot germplasm resources in different regions showed a degree of correlation with their geographical distribution.

CONCLUSION

Comparative analysis of chloroplast genomes of 146 Japanese apricot resources was performed to analyze the used to explore the genetic relationship and genetic diversity among Japanese apricot resources with different geographical distributions, providing some reference for the origin and evolution of Japanese apricot.

Chloroplast Genomes of Genus Tilia: Comparative Genomics and Molecular Evolution

Tilia is a complex genus in the family Malvaceae that has high ecological and economical values. Owing to the lack of sufficient distinguishable morphological and molecular characteristics, interspecific relationships in this genus are not clear. Chloroplast (cp) genomes are small, meanwhile most angiosperms usually undergo matrilineal inheritance. Consequently, they can be used in molecular systematics and phylogenetic analyses. Here, we sequenced and assembled cp genomes from T. endochrysea, T. tomentosa, T. miqueliana, T. americana and T. cordata, and compared them with those of seven previously reported Tilia species. Similar gene contents, gene orders and GC contents existed among the 12 cp genomes, which ranged from 162,564 to 162,855 bp and encoded 113 unique genes. Abundant simple sequence repeats (119–127) and dispersed repeats (97–135) were detected in Tilia cp genomes. In total, 11 hypervariable regions were identified that could be suitable for species identification and phylogenetic studies. A phylogenetic analysis of Malvaceae based on 5 hypervariable genes (matK + ndhF + rpoB + rpoC2+ycf1) revealed that all eight subfamilies were monophyletic groups. Additionally, the genus Tilia was divided into three groups on the basis of all 521 molecular variation loci. The current study provides valuable insights into the genomic evolution of the genus Tilia.

The Manchurian Walnut Genome: Insights into Juglone and Lipid Biosynthesis

BACKGROUND

Manchurian walnut (Juglans mandshurica Maxim.) is a tree with multiple industrial uses and medicinal properties in the Juglandaceae family (walnuts and hickories). J. mandshurica produces juglone, which is a toxic allelopathic agent and has potential utilization value. Furthermore, the seed of J. mandshurica is rich in various unsaturated fatty acids and has high nutritive value.

FINDINGS

Here, we present a high-quality chromosome-scale reference genome assembly and annotation for J. mandshurica (n = 16) with a contig N50 of 21.4 Mb by combining PacBio high-fidelity reads with high-throughput chromosome conformation capture data. The assembled genome has an estimated sequence size of 548.7 Mb and consists of 657 contigs, 623 scaffolds, and 40,453 protein-coding genes. In total, 60.99% of the assembled genome consists of repetitive sequences. Sixteen super-scaffolds corresponding to the 16 chromosomes were assembled, with a scaffold N50 length of 33.7 Mb and a BUSCO complete gene percentage of 98.3%. J. mandshurica displays a close sequence relationship with Juglans cathayensis, with a divergence time of 13.8 million years ago. Combining the high-quality genome, transcriptome, and metabolomics data, we constructed a gene-to-metabolite network and identified 566 core and conserved differentially expressed genes, which may be involved in juglone biosynthesis. Five CYP450 genes were found that may contribute to juglone accumulation. NAC, bZip, NF-YA, and NF-YC are positively correlated with the juglone content. Some candidate regulators (e.g., FUS3, ABI3, LEC2, and WRI1 transcription factors) involved in the regulation of lipid biosynthesis were also identified.

CONCLUSIONS

Our genomic data provide new insights into the evolution of the walnut genome and create a new platform for accelerating molecular breeding and improving the comprehensive utilization of these economically important tree species.

A Comparative Analysis of the Chloroplast Genomes of Four Polygonum Medicinal Plants

Polygonum is a generalized genus of the Polygonaceae family that includes various herbaceous plants. In order to provide aid in understanding the evolutionary and phylogenetic relationship in Polygonum at the chloroplast (cp) genome-scale level, we sequenced and annotated the complete chloroplast genomes of four Polygonum species using next-generation sequencing technology and CpGAVAS. Then, repeat sequences, IR contractions, and expansion and transformation sites of chloroplast genomes of four Polygonum species were studied, and a phylogenetic tree was built using the chloroplast genomes of Polygonum. The results indicated that the chloroplast genome construction of Polygonum also displayed characteristic four types of results, comparable to the published chloroplast genome of recorded angiosperms. The chloroplast genomes of the four Polygonum plants are highly consistent in genome size (159,015 bp-163,461 bp), number of genes (112 genes, including 78 protein-coding genes, 30 tRNA genes, and 4 rRNA genes), gene types, gene order, codon usage, and repeat sequence distribution, which identifies the high preservation among the Polygonum chloroplast genomes. The Polygonum phylogenetic tree was recreated by a full sequence of the chloroplast genome, which illustrates that the P. bistorta, P. orientale, and P. perfoliatum are divided into the same branch, and P. aviculare belongs to Fallopia. The precise system site of lots base parts requires further verification, but the study would provide a basis for developing the available genetic resources and evolutionary relationships of Polygonum.

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.

Gene loss, genome rearrangement, and accelerated substitution rates in plastid genome of Hypericum ascyron (Hypericaceae)

BACKGROUND

Comparative genomic analysis exhibits dynamic evolution of plastid genome (plastome) in the clusioid clade of Malpighiales, which comprise five families, including multiple inversions and gene losses. Little is known about the plastome evolution in Hypericaceae, a large family in the clade. Only the plastome of one species, Cratoxylum cochinchinense, has been published.

RESULTS

We generated a complete plastome sequence for Hypericum ascyron, providing the first complete plastome from the tribe Hypericeae (Hypericaceae). The H. ascyron plastome exhibits dynamic changes in gene and intron content, structure, and sequence divergence compared to the C. cochinchinense plastome from the tribe Cratoxyleae (Hypericaceae). Transcriptome data determined the evolutionary fate of the missing plastid genes infA, rps7, rps16, rpl23, and rpl32 in H. ascyron. Putative functional transfers of infA, rps7, and rpl32 were detected to the nucleus, whereas rps16 and rpl23 were substituted by nuclear-encoded homologs. The plastid rpl32 was integrated into the nuclear-encoded SODcp gene. Our findings suggested that the transferred rpl32 had undergone subfunctionalization by duplication rather than alternative splicing. The H. ascyron plastome rearrangements involved seven inversions, at least three inverted repeat (IR) boundary shifts, which generated gene relocations and duplications. Accelerated substitution rates of plastid genes were observed in the H. ascyron plastome compared with that of C. cochinchinense plastid genes. The higher substitution rates in the accD and clpP were correlated with structural change, including a large insertion of amino acids and losses of two introns, respectively. In addition, we found evidence of positive selection of the clpP, matK, and rps3 genes in the three branches related to H. ascyron. In particular, the matK gene was repeatedly under selection within the family Hypericaceae. Selective pressure in the H. ascyron matK gene was associated with the loss of trnK-UUU and relocation into the IR region.

CONCLUSIONS

The Hypericum ascyron plastome sequence provides valuable information for improving the understanding of plastome evolution among the clusioid of the Malpighiales. Evidence for intracellular gene transfer from the plastid to the nucleus was detected in the nuclear transcriptome, providing insight into the evolutionary fate of plastid genes in Hypericaceae.

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.

Plastid phylogenomics provides novel insights into the infrafamilial relationship of Polypodiaceae

The polygrammoids (Polypodiaceae) are the most species-rich and diversified epiphytic fern lineages, and hold an important role to understand the deep diverging events and rapid adaptation to changing environments in the plant tree of life. Despite progress in the phylogeny of this group of ferns in previous multilocus phylogenetic studies, uncertainty remains especially in backbone relationships among closely related clades, and the phylogenetic placement of recalcitrant species or lineages. Here, we investigated the deep phylogenetic relationships within Polypodiaceae by sampling all major lineages and using 81 plastid genomes (plastomes), of which 70 plastomes were newly sequenced with high-throughput sequencing technology. Based on parsimony, maximum-likelihood, Bayesian and multispecies coalescent analyses of genome skimming data, we achieved a better resolution of the backbone phylogeny of Polypodiaceae. Using simulated data matrices, we detected that potential phylogenetic artefacts, such as long-branch attraction and insufficient taxonomic sampling, may have a confounding impact on the incongruence of phylogenetic inferences. Furthermore, our phylogenetic analyses offer greater resolution than previous multilocus studies, providing a robust framework for future phylogenetic implications on the subfamilial taxonomy of Polypodiaceae. Our phylogenomic study not only demonstrates the advantage of a character-rich plastome dataset for resolving the recalcitrant lineages that have undergone rapid radiation, but also sheds new light on integrative explorations understanding the evolutionary history of large fern groups in the genomic era.

Evolution of Rosaceae Plastomes Highlights Unique Cerasus Diversification and Independent Origins of Fruiting Cherry

Rosaceae comprises numerous types of economically important fruits, ornamentals, and timber. The lack of plastome characteristics has blocked our understanding of the evolution of plastome and plastid genes of Rosaceae crops. Using comparative genomics and phylogenomics, we analyzed 121 Rosaceae plastomes of 54 taxa from 13 genera, predominantly including Cerasus (true cherry) and its relatives. To our knowledge, we generated the first comprehensive map of genomic variation across Rosaceae plastomes. Contraction/expansion of inverted repeat regions and sequence losses of the two single-copy regions underlie large genomic variations in size among Rosaceae plastomes. Plastid protein-coding genes were characterized with a high proportion (over 50%) of synonymous variants and insertion-deletions with multiple triplets. Five photosynthesis-related genes were specially selected in perennial woody trees. Comparative genomic analyses implied divergent evolutionary patterns between pomaceous and drupaceous trees. Across all examined plastomes, unique and divergent evolution was detected in Cerasus plastomes. Phylogenomic analyses and molecular dating highlighted the relatively distant phylogenetic relationship between Cerasus and relatives (Microcerasus, Amygdalus, Prunus, and Armeniaca), which strongly supported treating the monophyletic true cherry group as a separate genus excluding dwarf cherry. High genetic differentiation and distinct phylogenetic relationships implied independent origins and domestication between fruiting cherries, particularly between Prunus pseudocerasus (Cerasus pseudocerasus) and P. avium (C. avium). Well-resolved maternal phylogeny suggested that cultivated P. pseudocerasus originated from Longmenshan Fault zone, the eastern edge of Himalaya-Hengduan Mountains, where it was subjected to frequent genomic introgression between its presumed wild ancestors and relatives.

Phylogenetic relationships in the Sorghum genus based on sequencing of the chloroplast and nuclear genes

Sorghum [Sorghum bicolor (L.) Moench] is an important food crop with a diverse gene pool residing in its wild relatives. A total of 15 sorghum accessions from the unexploited wild gene pool of the Sorghum genus, representing the five subgenera, were sequenced, and the complete chloroplast genomes and 99 common single-copy concatenated nuclear genes were assembled. Annotation of the chloroplast genomes identified a total of 81 protein-coding genes, 38 tRNA, and four rRNA genes. The gene content and gene order among the species was identical. A total of 153 nonsynonymous amino acid changes in 40 genes were identified across the species. Phylogenetic analysis of both the whole chloroplast genome and nuclear genes revealed a similar topology with two distinct clades within the genus. The species within the subgenera Eusorghum, Chaetosorghum, and Heterosorghum clustered in one clade, whereas the species within the subgenera Parasorghum and Stiposorghum clustered in a second clade. However, the subgenera Parasorghum and Stiposorghum were not monophyletic, suggesting the need for further research to resolve the relationships within this group. The close relationship between the two monotypic subgenera Chaetosorghum and Heterosorghum suggests that species within these subgenera could be considered as one group. This analysis provides an improved understanding of the genetic relationships within the Sorghum genus and defines diversity in wild sorghum species that may be useful for crop improvement.

Whole genome based insights into the phylogeny and evolution of the Juglandaceae

BACKGROUND

The walnut family (Juglandaceae) contains commercially important woody trees commonly called walnut, wingnut, pecan and hickory. Phylogenetic relationships and diversification within the Juglandaceae are classic and hot scientific topics that have been elucidated by recent fossil, morphological, molecular, and (paleo) environmental data. Further resolution of relationships among and within genera is still needed and can be achieved by analysis of the variation of chloroplast, mtDNA, and nuclear genomes.

RESULTS

We reconstructed the backbone phylogenetic relationships of Juglandaceae using organelle and nuclear genome data from 27 species. The divergence time of Juglandaceae was estimated to be 78.7 Mya. The major lineages diversified in warm and dry habitats during the mid-Paleocene and early Eocene. The plastid, mitochondrial, and nuclear phylogenetic analyses all revealed three subfamilies, i.e., Juglandoideae, Engelhardioideae, Rhoipteleoideae. Five genera of Juglandoideae were strongly supported. Juglandaceae were estimated to have originated during the late Cretaceous, while Juglandoideae were estimated to have originated during the Paleocene, with evidence for rapid diversification events during several glacial and geological periods. The phylogenetic analyses of organelle sequences and nuclear genome yielded highly supported incongruence positions for J. cinerea, J. hopeiensis, and Platycarya strobilacea. Winged fruit were the ancestral condition in the Juglandoideae, but adaptation to novel dispersal and regeneration regimes after the Cretaceous-Paleogene boundary led to the independent evolution of zoochory among several genera of the Juglandaceae.

CONCLUSIONS

A fully resolved, strongly supported, time-calibrated phylogenetic tree of Juglandaceae can provide an important framework for studying classification, diversification, biogeography, and comparative genomics of plant lineages. Our addition of new, annotated whole chloroplast genomic sequences and identification of their variability informs the study of their evolution in walnuts (Juglandaceae).

Whole genome based insights into the phylogeny and evolution of the Juglandaceae

BACKGROUND

The walnut family (Juglandaceae) contains commercially important woody trees commonly called walnut, wingnut, pecan and hickory. Phylogenetic relationships and diversification within the Juglandaceae are classic and hot scientific topics that have been elucidated by recent fossil, morphological, molecular, and (paleo) environmental data. Further resolution of relationships among and within genera is still needed and can be achieved by analysis of the variation of chloroplast, mtDNA, and nuclear genomes.

RESULTS

We reconstructed the backbone phylogenetic relationships of Juglandaceae using organelle and nuclear genome data from 27 species. The divergence time of Juglandaceae was estimated to be 78.7 Mya. The major lineages diversified in warm and dry habitats during the mid-Paleocene and early Eocene. The plastid, mitochondrial, and nuclear phylogenetic analyses all revealed three subfamilies, i.e., Juglandoideae, Engelhardioideae, Rhoipteleoideae. Five genera of Juglandoideae were strongly supported. Juglandaceae were estimated to have originated during the late Cretaceous, while Juglandoideae were estimated to have originated during the Paleocene, with evidence for rapid diversification events during several glacial and geological periods. The phylogenetic analyses of organelle sequences and nuclear genome yielded highly supported incongruence positions for J. cinerea, J. hopeiensis, and Platycarya strobilacea. Winged fruit were the ancestral condition in the Juglandoideae, but adaptation to novel dispersal and regeneration regimes after the Cretaceous-Paleogene boundary led to the independent evolution of zoochory among several genera of the Juglandaceae.

CONCLUSIONS

A fully resolved, strongly supported, time-calibrated phylogenetic tree of Juglandaceae can provide an important framework for studying classification, diversification, biogeography, and comparative genomics of plant lineages. Our addition of new, annotated whole chloroplast genomic sequences and identification of their variability informs the study of their evolution in walnuts (Juglandaceae).

Evolutionary patterns of nucleotide substitution rates in plastid genomes of Quercus

Molecular evolution, including nucleotide substitutions, plays an important role in understanding the dynamics and mechanisms of species evolution. Here, we sequenced whole plastid genomes (plastomes) of Quercus fabri, Quercus semecarpifolia, Quercus engleriana, and Quercus phellos and compared them with 14 other Quercus plastomes to explore their evolutionary relationships using 67 shared protein-coding sequences. While many previously identified evolutionary relationships were found, our findings do not support previous research which retrieve Quercus subg. Cerris sect. Ilex as a monophyletic group, with sect. Ilex found to be polyphyletic and composed of three strongly supported lineages inserted between sections Cerris and Cyclobalanposis. Compared with gymnosperms, Quercus plastomes showed higher evolutionary rates (Dn/Ds = 0.3793). Most protein-coding genes experienced relaxed purifying selection, and the high Dn value (0.1927) indicated that gene functions adjusted to environmental changes effectively. Our findings suggest that gene interval regions play an important role in Quercus evolution. We detected greater variation in the intergenic regions (trnH-psbA, trnK_UUU-rps16, trnfM_CAU-rps14, trnS_GCU-trnG_GCC, and atpF-atpH), intron losses (petB and petD), and pseudogene loss and degradation (ycf15). Additionally, the loss of some genes suggested the existence of gene exchanges between plastid and nuclear genomes, which affects the evolutionary rate of the former. However, the connective mechanism between these two genomes is still unclear.

A Repertory of Rearrangements and the Loss of an Inverted Repeat Region in Passiflora Chloroplast Genomes

Chloroplast genomes (cpDNA) in angiosperms are usually highly conserved. Although rearrangements have been observed in some lineages, such as Passiflora, the mechanisms that lead to rearrangements are still poorly elucidated. In the present study, we obtained 20 new chloroplast genomes (18 species from the genus Passiflora, and Dilkea retusa and Mitostemma brevifilis from the family Passifloraceae) in order to investigate cpDNA evolutionary history in this group. Passiflora cpDNAs vary in size considerably, with ∼50 kb between shortest and longest. Large inverted repeat (IR) expansions were identified, and at the extreme opposite, the loss of an IR was detected for the first time in Passiflora, a rare event in angiosperms. The loss of an IR region was detected in Passiflora capsularis and Passiflora costaricensis, a species in which occasional biparental chloroplast inheritance has previously been reported. A repertory of rearrangements such as inversions and gene losses were detected, making Passiflora one of the few groups with complex chloroplast genome evolution. We also performed a phylogenomic study based on all the available cp genomes and our analysis implies that there is a need to reconsider the taxonomic classifications of some species in the group.

Characterization and Dynamics of Intracellular Gene Transfer in Plastid Genomes of Viola (Violaceae) and Order Malpighiales

Functional gene transfer from organelles to the nucleus, known as intracellular gene transfer (IGT), is an ongoing process in flowering plants. The complete plastid genomes (plastomes) of two Ulleung island endemic violets, Viola ulleungdoensis and V. woosanensis, were characterized, revealing a lack of the plastid-encoded infA, rpl32, and rps16 genes. In addition, functional replacement of the three plastid-encoded genes in the nucleus was confirmed within the genus Viola and the order Malpighiales. Three strategies for the acquisition of a novel transit peptide for successful IGT were identified in the genus Viola. Nuclear INFA acquired a novel transit peptide with very low identity between these proteins, whereas the nuclear RPL32 gene acquired an existing transit peptide via fusion with the nuclear-encoded plastid-targeted SOD gene (Cu-Zn superoxide dismutase superfamily) as one exon, and translated both proteins in the cytosol using alternative mRNA splicing. Nuclear RPS16 contains an internal transit peptide without an N-terminal extension. Gene loss or pseudogenization of the plastid-borne rpl32 and rps16 loci was inferred to occur in the common ancestor of the genus Viola based on an infrageneric phylogenetic framework in Korea. Although infA was lost in the common ancestor of the order Malpighiales, the rpl32 and rps16 genes were lost multiple times independently within the order. Our current study sheds additional light on plastid genome composition and IGT mechanisms in the violet genus and in the order Malpighiales.

Reference-based assembly of chloroplast genome from leaf transcriptome data of Pterocarpus santalinus

Chloroplast genome sequencing is an essential tool to understand genome evolution and phylogenetic relationship. The available methods for constructing chloroplast genome include chloroplast enrichment followed by long overlapping PCR or extraction and assembly of chloroplast-specific reads from whole-genome datasets. In the present study, we propose an alternate strategy of extraction and assembly of chloroplast-specific reads from leaf transcriptome data of Pterocarpus santalinus using bowtie2 aligner program. The assembled genome was compared with the published chloroplast genome of P. santalinus for genome size, number of predicted genes, microsatellite repeat motifs, and nucleotide repeats. A near-complete chloroplast genome was assembled from the transcriptome reads. The proposed method requires less computational time and know-how, limited virtual memory, and is cost-effective when compared to whole-genome sequencing. Assembly of Cp genome from transcriptome data will enhance the resolution of phylogenetic studies through comparative plastome analysis, facilitate accurate species/genotype discrimination and accelerate the development of transplastomic plants with enhanced biotic and abiotic tolerance.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02943-0.

Rampant Nuclear Transfer and Substitutions of Plastid Genes in Passiflora

Gene losses in plastid genomes (plastomes) are often accompanied by functional transfer to the nucleus or substitution of an alternative nuclear-encoded gene. Despite the highly conserved gene content in plastomes of photosynthetic land plants, recent gene loss events have been documented in several disparate angiosperm clades. Among these lineages, Passiflora lacks several essential ribosomal genes, rps7, rps16, rpl20, rpl22, and rpl32, the two largest plastid genes, ycf1 and ycf2, and has a highly divergent rpoA. Comparative transcriptome analyses were performed to determine the fate of the missing genes in Passiflora. Putative functional transfers of rps7, rpl22, and rpl32 to nucleus were detected, with the nuclear transfer of rps7, representing a novel event in angiosperms. Plastid-encoded rps7 was transferred into the intron of a nuclear-encoded plastid-targeted thioredoxin m-type gene, acquiring its plastid transit peptide (TP). Plastid rpl20 likely experienced a novel substitution by a duplicated, nuclear-encoded mitochondrial-targeted rpl20 that has a similar gene structure. Additionally, among rosids, evidence for a third independent transfer of rpl22 in Passiflora was detected that gained a TP from a nuclear gene containing an organelle RNA recognition motif. Nuclear transcripts representing rpoA, ycf1, and ycf2 were not detected. Further analyses suggest that the divergent rpoA remains functional and that the gene is under positive or purifying selection in different clades. Comparative analyses indicate that alternative translocon and motor protein complexes may have substituted for the loss of ycf1 and ycf2 in Passiflora.

Comparative Analysis of the Complete Chloroplast Genomes of Four Chestnut Species (Castanea)

Chloroplast (cp) DNA genomes are traditional workhorses for studying the evolution of species and reconstructing phylogenetic relationships in plants. Species of the genus Castanea (chestnuts and chinquapins) are valued as a source of nuts and timber wherever they grow, and chestnut species hybrids are common. We compared the cp genomes of C. mollissima, C. seguinii, C. henryi, and C. pumila. These cp genomes ranged from 160,805 bp to 161,010 bp in length, comprising a pair of inverted repeat (IR) regions (25,685 to 25,701 bp) separated by a large single-copy (LSC) region (90,440 to 90,560 bp) and a small single-copy (SSC) region (18,970 to 19,049 bp). Each cp genome encoded the same 113 genes; 82–83 protein-coding genes, 30 transfer RNA genes, and four ribosomal RNA genes. There were 18 duplicated genes in the IRs. Comparative analysis of cp genomes revealed that rpl22 was absent in all analyzed species, and the gene ycf1 has been pseudo-genized in all Chinese chestnuts except C. pumlia. We analyzed the repeats and nucleotide substitutions in these plastomes and detected several highly variable regions. The phylogenetic analyses based on plastomes confirmed the monophyly of Castanea species.

The chloroplast genome of Prunus zhengheensis: Genome comparative and phylogenetic relationships analysis

Prunus zhengheensis is a novel species originated in Fujian province, China. However, there is no further information available on its classification and molecular biology study. In this study, we first report the complete chloroplast (cp) genome sequence of P. zhengheensis. The cp genome of P. zhengheensis is 158,106 bp and GC content is 36.73%, is a circular structure composed of LSC (large single copy), SSC (small single copy), and IR (inverted repeat) regions, with the size of the three regions being 86,321 bp, 18,999 bp and 26,393 bp, respectively. The cp genome of P. zhengheensis contains 130 genes, and 242 SSRs are identified in the cp genome. The comparative analysis of cp genomes in eight Prunus plants demonstrates the subtle divergences occur in the protein-coding gene rps18, rps12, psbF, rpl33, matK, and rbcL, and that the K_A/K_S nucleotide substitution ratio of the ndhF of P. zhengheensis and P. armeniaca is 1.79636. The phylogenetic results indicate that the P. zhengheensis is closely related to P. mume, compared to other species of Prunus. Our research results provide the important genomic information for molecular phylogeny of P. zhengheensis.

Phylogenetic relationships and molecular evolution of woody forest tree family Aceraceae based on plastid phylogenomics and nuclear gene variations

The forest tree family Aceraceae is widespread in the northern hemisphere and it has ecological and economic importance. However, the phylogenetic relationships and classifications within the family are still controversial due to transitional intraspecific morphological characteristics and introgression hybridization among species. In this study, we determined the evolutionary relationships and molecular evolution of Aceraceae based on plastid phylogenomics and two nuclear gene variations. Phylogenetic analysis based on the plastid genomes suggested that Aceraceae species can be divided into two larger sub-clades corresponding to the two genera Acer and Dipteronia. Conjoint analysis of the plastid and nuclear gene sequences supported the classification with two genera in the family. Molecular dating showed that the two genera diverged 60.2 million years ago, which is generally consistently with previously reported results. Divergence hotspots and positively selected genes identified in the plastid genomes could be useful genetic resources in Aceraceae.

The evolutionary fate of rpl32 and rps16 losses in the Euphorbia schimperi (Euphorbiaceae) plastome

Gene transfers from mitochondria and plastids to the nucleus are an important process in the evolution of the eukaryotic cell. Plastid (pt) gene losses have been documented in multiple angiosperm lineages and are often associated with functional transfers to the nucleus or substitutions by duplicated nuclear genes targeted to both the plastid and mitochondrion. The plastid genome sequence of Euphorbia schimperi was assembled and three major genomic changes were detected, the complete loss of rpl32 and pseudogenization of rps16 and infA. The nuclear transcriptome of E. schimperi was sequenced to investigate the transfer/substitution of the rpl32 and rps16 genes to the nucleus. Transfer of plastid-encoded rpl32 to the nucleus was identified previously in three families of Malpighiales, Rhizophoraceae, Salicaceae and Passifloraceae. An E. schimperi transcript of pt SOD-1-RPL32 confirmed that the transfer in Euphorbiaceae is similar to other Malpighiales indicating that it occurred early in the divergence of the order. Ribosomal protein S16 (rps16) is encoded in the plastome in most angiosperms but not in Salicaceae and Passifloraceae. Substitution of the E. schimperi pt rps16 was likely due to a duplication of nuclear-encoded mitochondrial-targeted rps16 resulting in copies dually targeted to the mitochondrion and plastid. Sequences of RPS16-1 and RPS16-2 in the three families of Malpighiales (Salicaceae, Passifloraceae and Euphorbiaceae) have high sequence identity suggesting that the substitution event dates to the early divergence within Malpighiales.

Plastome analysis unveils Inverted Repeat (IR) expansion and positive selection in Sea Lavenders (Limonium, Plumbaginaceae, Limonioideae, Limonieae)

The genus Limonium, commonly known as Sea Lavenders, is one of the most species-rich genera of the family Plumbaginaceae. In this study, two new plastomes for the genus Limonium, viz. L. tetragonum and L. bicolor, were sequenced and compared to available Limonium plastomes, viz. L. aureum and L. tenellum, to understand the gene content and structural variations within the family. The loss of the rpl16 intron and pseudogenisation of rpl23 was observed. This study reports, for the first time, expansion of the IRs to include the ycf1 gene in Limonium plastomes, incongruent with previous studies. Two positively selected genes, viz. ndhF and ycf2, were identified. Furthermore, putative barcodes are proposed for the genus, based on the nucleotide diversity of four Limonium plastomes.

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

Background

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

Results

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

Conclusions

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

Supplementary Information

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

Phylogenetic relationships in Chinese oaks (Fagaceae, Quercus): Evidence from plastid genome using low-coverage whole genome sequencing

China is a second center of oak diversity but with less intensively systematic studies. Here, with 49 species representing all four sections in China, we firstly gave insight into the comprehensive phylogenetic relationships of Chinese oaks based on 54 complete plastid genomes. Our results recovered a robust phylogenetic framework and provided strong support for most nodes. The phylogenetic tree supported Quercus section Ilex as not monophyletic, in which Quercus section Cyclobalanopsis and Quercus section Cerris were nested. Most likely, incomplete lineage sorting and/or introgression among ancestral lineages in these three sections resulted in this complex pattern. The current distribution, diversification and molecular differentiation of Q. sect. Ilex in China are likely consequences of local adaptation to the geographic and paleoclimatic changes, which were driven by the uplift of Tibetan Plateau, the Hengduan Mountains and the Himalayas.

Complete chloroplast genome of a rare and endangered plant species Osteomeles subrotunda: genomic features and phylogenetic relationships with other Rosaceae plants

Osteomeles subrotunda is a rare and endangered plant species with extremely small populations. In our study, we sequenced the complete chloroplast (CP) genome of O. subrotunda and described its structural organization, and performed comparative genomic analyses with other Rosaceae CP genomes. The plastome of O. subrotunda was 159,902 bp in length with 36.6% GC content and contained a pair of inverted repeats of 26,367 bp which separated a large single-copy region of 87,933 bp and a small single-copy region of 19,235 bp. The CP genome included 130 genes, of which 85 were protein-coding genes, 37 were transfer RNAs, and eight were ribosomal RNAs. Two genes, rps19 and ycf1, which are located at the borders of IRB/SSC and IRB/LSC, were presumed to be pseudogenes. A total of 61 SSRs were detected, of which, 59 loci were mono-nucleotide repeats, and two were di-nucleotide repeats. The phylogenic analysis indicated that the 14 Rosaceae species were divided into three groups, among which O. subrotunda grouped with P. rupicola, E. japonica, P. pashia, C. japonica, S. torminalis, and M. florentina, and it was found to be a sister clade to C. japonica. Our newly sequenced CP genome of O. subrotunda will provide essential data for further studies on population genetics and biodiversity.

Green giant-a tiny chloroplast genome with mighty power to produce high-value proteins: history and phylogeny

Free-living cyanobacteria were entrapped by eukaryotic cells ~2 billion years ago, ultimately giving rise to chloroplasts. After a century of debate, the presence of chloroplast DNA was demonstrated in the 1960s. The first chloroplast genomes were sequenced in the 1980s, followed by ~100 vegetable, fruit, cereal, beverage, oil and starch/sugar crop chloroplast genomes in the past three decades. Foreign genes were expressed in isolated chloroplasts or intact plant cells in the late 1980s and stably integrated into chloroplast genomes, with typically maternal inheritance shown in the 1990s. Since then, chloroplast genomes conferred the highest reported levels of tolerance or resistance to biotic or abiotic stress. Although launching products with agronomic traits in important crops using this concept has been elusive, commercial products developed include enzymes used in everyday life from processing fruit juice, to enhancing water absorption of cotton fibre or removal of stains as laundry detergents and in dye removal in the textile industry. Plastid genome sequences have revealed the framework of green plant phylogeny as well as the intricate history of plastid genome transfer events to other eukaryotes. Discordant historical signals among plastid genes suggest possible variable constraints across the plastome and further understanding and mitigation of these constraints may yield new opportunities for bioengineering. In this review, we trace the evolutionary history of chloroplasts, status of autonomy and recent advances in products developed for everyday use or those advanced to the clinic, including treatment of COVID-19 patients and SARS-CoV-2 vaccine.

Comparative chloroplast genomes: insights into the evolution of the chloroplast genome of Camellia sinensis and the phylogeny of Camellia

BACKGROUND

Chloroplast genome resources can provide useful information for the evolution of plant species. Tea plant (Camellia sinensis) is among the most economically valuable member of Camellia. Here, we determined the chloroplast genome of the first natural triploid Chinary type tea ('Wuyi narcissus' cultivar of Camellia sinensis var. sinensis, CWN) and conducted the genome comparison with the diploid Chinary type tea (Camellia sinensis var. sinensis, CSS) and two types of diploid Assamica type teas (Camellia sinensis var. assamica: Chinese Assamica type tea, CSA and Indian Assamica type tea, CIA). Further, the evolutionary mechanism of the chloroplast genome of Camellia sinensis and the relationships of Camellia species based on chloroplast genome were discussed.

RESULTS

Comparative analysis showed the evolutionary dynamics of chloroplast genome of Camellia sinensis were the repeats and insertion-deletions (indels), and distribution of the repeats, indels and substitutions were significantly correlated. Chinese tea and Indian tea had significant differences in the structural characteristic and the codon usage of the chloroplast genome. Analysis of sequence characterized amplified region (SCAR) using sequences of the intergenic spacers (trnE/trnT) showed none of 292 different Camellia sinensis cultivars had similar sequence characteristic to triploid CWN, but the other four Camellia species did. Estimations of the divergence time showed that CIA diverged from the common ancestor of two Assamica type teas about 6.2 Mya (CI: 4.4-8.1 Mya). CSS and CSA diverged to each other about 0.8 Mya (CI: 0.4-1.5 Mya). Moreover, phylogenetic clustering was not exactly consistent with the current taxonomy of Camellia.

CONCLUSIONS

The repeat-induced and indel-induced mutations were two important dynamics contributed to the diversification of the chloroplast genome in Camellia sinensis, which were not mutually exclusive. Chinese tea and Indian tea might have undergone different selection pressures. Chloroplast transfer occurred during the polyploid evolution in Camellia sinensis. In addition, our results supported the three different domestication origins of Chinary type tea, Chinese Assamica type tea and Indian Assamica type tea. And, the current classification of some Camellia species might need to be further discussed.

Plastid Genomes of Flowering Plants: Essential Principles

The plastid genome (plastome ) has proved a valuable source of data for evaluating evolutionary relationships among angiosperms. Through basic and applied approaches, plastid transformation technology offers the potential to understand and improve plant productivity, providing food, fiber, energy, and medicines to meet the needs of a burgeoning global population. The growing genomic resources available to both phylogenetic and biotechnological investigations is allowing novel insights and expanding the scope of plastome research to encompass new species. In this chapter, we present an overview of some of the seminal and contemporary research that has contributed to our current understanding of plastome evolution and attempt to highlight the relationship between evolutionary mechanisms and the tools of plastid genetic engineering.

Phylogeny of Euclidieae (Brassicaceae) based on plastome and nuclear ribosomal DNA data

Euclidieae, a morphologically diverse tribe in the family Brassicaceae (Cruciferae), consists of 29 genera and more than 150 species distributed mainly in Asia. Prior phylogenetic analyses on Euclidieae are inadequate. In this study, sequence data from the plastid genome and nuclear ribosomal DNA of 72 species in 27 genera of Euclidieae were used to infer the inter- and intra-generic relationships within. The well-resolved and strongly supported plastome phylogenies revealed that Euclidieae could be divided into five clades. Both Cymatocarpus and Neotorularia are polyphyletic in nuclear and plastome phylogenies. Besides, the conflicts of systematic positions of three species of Braya and two species of Solms-laubachia s.l. indicated that hybridization and or introgression might have happened during the evolutionary history of the tribe. Results from divergence-time analyses suggested an early Miocene origin of Euclidieae, and it probably originated from the Central Asia, Pamir Plateau and West Himalaya. In addition, multiple ndh genes loss and pseudogenization were detected in eight species based on comparative genomic study.

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.

Evolutionary directions of single nucleotide substitutions and structural mutations in the chloroplast genomes of the family Calycanthaceae

BACKGROUND

Chloroplast genome sequence data is very useful in studying/addressing the phylogeny of plants at various taxonomic ranks. However, there are no empirical observations on the patterns, directions, and mutation rates, which are the key topics in chloroplast genome evolution. In this study, we used Calycanthaceae as a model to investigate the evolutionary patterns, directions and rates of both nucleotide substitutions and structural mutations at different taxonomic ranks.

RESULTS

There were 2861 polymorphic nucleotide sites on the five chloroplast genomes, and 98% of polymorphic sites were biallelic. There was a single-nucleotide substitution bias in chloroplast genomes. A → T or T → A (2.84%) and G → C or C → G (3.65%) were found to occur significantly less frequently than the other four transversion mutation types. Synonymous mutations kept balanced pace with nonsynonymous mutations, whereas biased directions appeared between transition and transversion mutations and among transversion mutations. Of the structural mutations, indels and repeats had obvious directions, but microsatellites and inversions were non-directional. Structural mutations increased the single nucleotide mutations rates. The mutation rates per site per year were estimated to be 0.14-0.34 × 10- 9 for nucleotide substitution at different taxonomic ranks, 0.64 × 10- 11 for indels and 1.0 × 10- 11 for repeats.

CONCLUSIONS

Our direct counts of chloroplast genome evolution events provide raw data for correctly modeling the evolution of sequence data for phylogenetic inferences.

Complete Chloroplast Genome Characterization of Oxalis Corniculata and Its Comparison with Related Species from Family Oxalidaceae

Oxalis corniculata L. (family Oxalidaceae) is a small creeper wood sorrel plant that grows well in moist climates. Despite being medicinally important, little is known about the genomics of this species. Here, we determined the complete chloroplast genome sequence of O. corniculata for the first time and compared it with other members of family Oxalidaceae. The genome was 152,189 bp in size and comprised of a pair of 25,387 bp inverted repeats (IR) that separated a large 83,427 bp single copy region (LSC) and a small 16,990 bp single copy region (SSC). The chloroplast genome of O. corniculata contains 131 genes with 83 protein coding genes, 40 tRNA genes, and 8 rRNA genes. The analysis revealed 46 microsatellites, of which 6 were present in coding sequences (CDS) regions, 34 in the LSC, 8 in the SSC, and 2 in the single IR region. Twelve palindromic repeats, 30 forward repeats, and 32 tandem repeats were also detected. Chloroplast genome comparisons revealed an overall high degree of sequence similarity between O. corniculata and O. drummondii and some divergence in the intergenic spacers of related species in Oxalidaceae. Furthermore, the seven most divergent genes (ccsA, clpP, rps8, rps15, rpl22, matK, and ycf1) among genomes were observed. Phylogenomic characterization on the basis of 60 shared genes revealed that O. corniculata is closely related to O. drummondii. The complete O. corniculata genome sequenced in the present study is a valuable resource for investigating the population and evolutionary genetics of family Oxalidaceae and can be used to identify related species.

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.