The complete chloroplast genome of Ginkgo biloba reveals the mechanism of inverted repeat contraction

Accumulation of Large Lineage-Specific Repeats Coincides with Sequence Acceleration and Structural Rearrangement in Plantago Plastomes

Repeats can mediate rearrangements and recombination in plant mitochondrial genomes and plastid genomes. While repeat accumulations are linked to heightened evolutionary rates and complex structures in specific lineages, debates persist regarding the extent of their influence on sequence and structural evolution. In this study, 75 Plantago plastomes were analyzed to investigate the relationships between repeats, nucleotide substitution rates, and structural variations. Extensive repeat accumulations were associated with significant rearrangements and inversions in the large inverted repeats (IRs), suggesting that repeats contribute to rearrangement hotspots. Repeats caused infrequent recombination that potentially led to substoichiometric shifting, supported by long-read sequencing. Repeats were implicated in elevating evolutionary rates by facilitating localized hypermutation, likely through DNA damage and repair processes. This study also observed a decrease in nucleotide substitution rates for loci translocating into IRs, supporting the role of biased gene conversion in maintaining lower substitution rates. Combined with known parallel changes in mitogenomes, it is proposed that potential dysfunction in nuclear-encoded genes associated with DNA replication, recombination, and repair may drive the evolution of Plantago organellar genomes. These findings contribute to understanding how repeats impact organellar evolution and stability, particularly in rapidly evolving plant lineages.

Characterization of the complete chloroplast genome of Wolffia arrhiza and comparative genomic analysis with relative Wolffia species

Lemnoideae, commonly referred to as the duckweed, are aquatic plants found worldwide. Wolffia species are known for their extreme reduction in size and complexity, lacking both roots and leaves, and they hold the distinction of being the smallest plants among angiosperms. Interestingly, it belongs to the Araceae family, despite its apparent morphological differences from land plants in the same family. Traditional morphological methods have limitations in classifying these plants, making molecular-level information essential. The chloroplast genome of Wolffia arrhiza is revealed that a total length of 169,602 bp and a total GC content of 35.78%. It follows the typical quadripartite structure, which includes a large single copy (LSC, 92,172 bp) region, a small single copy (SSC, 13,686 bp) region, and a pair of inverted repeat (IR, 31,872 bp each) regions. There are 131 genes characterized, comprising 86 Protein-Coding Genes, 37 Transfer RNA (tRNA) genes, and 8 ribosomal RNA (rRNA) genes. Moreover, 48 simple sequence repeats and 32 long repeat sequences were detected. Comparative analysis between W. arrhiza and six other Lemnoideae species identified 12 hotspots of high nucleotide diversity. In addition, a phylogenetic analysis was performed using 14 species belonging to the Araceae family and one external species as an outgroup. This analysis unveiled W. arrhiza and Wolffia globosa as closely related sister species. Therefore, this research has revealed the complete chloroplast genome data of W. arrhiza, offering a more detailed understanding of its evolutionary position and phylogenetic categorization within the Lemnoideae subfamily.

The Development of Plant Genome Sequencing Technology and Its Conservation and Application in Endangered Gymnosperms

Genome sequencing is widely recognized as a fundamental pillar in genetic research and legal studies of biological phenomena, providing essential insights for genetic investigations and legal analyses of biological events. The field of genome sequencing has experienced significant progress due to rapid improvements in scientific and technological developments. These advancements encompass not only significant improvements in the speed and quality of sequencing but also provide an unparalleled opportunity to explore the subtle complexities of genomes, particularly in the context of rare species. Such a wide range of possibilities has successfully supported the validation of plant gene functions and the refinement of precision breeding methodologies. This expanded scope now includes a comprehensive exploration of the current state and conservation efforts of gymnosperm gene sequencing, offering invaluable insights into their genomic landscapes. This comprehensive review elucidates the trajectory of development and the diverse applications of genome sequencing. It encompasses various domains, including crop breeding, responses to abiotic stress, species evolutionary dynamics, biodiversity, and the unique challenges faced in the conservation and utilization of gymnosperms. It highlights both ongoing challenges and the unveiling of forthcoming developmental trajectories.

The complete plastid genome provides insight into maternal plastid inheritance mode of the living fossil plant Ginkgo biloba

In the current research, we focus on uniparental inheritance of chloroplast genome of the living fossil plant, Ginkgo biloba L., one of the gymnosperms, using genomic data.•Our results provide strong genomic evidence to support plastid maternal inheritance mode of G. biloba, which is different from most other gymnosperms.•The combination of manually genetic crosses and genomic data is proved to be an efficient way to investigate the inheritance mode of chloroplasts genome in land plants.•The current research also provides a case study for further research of plastid inheritance in gymnosperms using genomic techniques, which will contribute to a better understanding of cytologically uniparental inheritance mode and evolutionary mechanism of plastids in both gymnosperms and angiosperms.

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.

Molecular Characterization and Genetic Diversity of Ginkgo (Ginkgo biloba L.) Based on Insertions and Deletions (InDel) Markers

As a "living fossil", ginkgo (Ginkgo biloba L.) has significant ornamental, medicinal, and timber value. However, the breeding improvement of ginkgo was limited by the lack of enough excellent germplasms and suitable molecular markers. Here, we characterized numerous polymorphic insertion/deletion (InDel) markers using RAD-seq in 12 different ginkgo cultivars. The total of 279,534 InDels identified were unequally distributed across 12 chromosomes in the ginkgo genome. Of these, 52.56% (146,919) and 47.44% (132,615) were attributed to insertions and deletions, respectively. After random selection and validation, 26 pairs of polymorphic primers were used for molecular diversity analysis in 87 ginkgo cultivars and clones. The average values of observed heterozygosity and polymorphism information were 0.625 and 0.517, respectively. The results of population structure analyses were similar to those of neighbor-joining and principal component analyses, which divided all germplasms into two distinct groups. Moreover, 11 ginkgo core collections accounted for approximately 12.64% of the total ginkgo germplasms obtained, representing well the allelic diversity of all original germplasms. Therefore, these InDels can be used for germplasm management and genetic diversity analyses in ginkgo and the core collections will be used effectively for ginkgo genetic improvement.

Insights into phylogenetic relationships in Pinus inferred from a comparative analysis of complete chloroplast genomes

BACKGROUND

Pinus is the largest genus of Pinaceae and the most primitive group of modern genera. Pines have become the focus of many molecular evolution studies because of their wide use and ecological significance. However, due to the lack of complete chloroplast genome data, the evolutionary relationship and classification of pines are still controversial. With the development of new generation sequencing technology, sequence data of pines are becoming abundant. Here, we systematically analyzed and summarized the chloroplast genomes of 33 published pine species.

RESULTS

Generally, pines chloroplast genome structure showed strong conservation and high similarity. The chloroplast genome length ranged from 114,082 to 121,530 bp with similar positions and arrangements of all genes, while the GC content ranged from 38.45 to 39.00%. Reverse repeats showed a shrinking evolutionary trend, with IRa/IRb length ranging from 267 to 495 bp. A total of 3,205 microsatellite sequences and 5,436 repeats were detected in the studied species chloroplasts. Additionally, two hypervariable regions were assessed, providing potential molecular markers for future phylogenetic studies and population genetics. Through the phylogenetic analysis of complete chloroplast genomes, we offered novel opinions on the genus traditional evolutionary theory and classification.

CONCLUSION

We compared and analyzed the chloroplast genomes of 33 pine species, verified the traditional evolutionary theory and classification, and reclassified some controversial species classification. This study is helpful in analyzing the evolution, genetic structure, and the development of chloroplast DNA markers in Pinus.

Comparative and phylogenetic analysis of complete chloroplast genomes from five Artemisia species

Artemisia Linn. is a large genus within the family Asteraceae that includes several important medicinal plants. Because of their similar morphology and chemical composition, traditional identification methods often fail to distinguish them. Therefore, developing an effective identification method for Artemisia species is an urgent requirement. In this study, we analyzed 15 chloroplast (cp) genomes, including 12 newly sequenced genomes, from 5 Artemisia species. The cp genomes from the five Artemisia species had a typical quadripartite structure and were highly conserved across species. They had varying lengths of 151,132-151,178 bp, and their gene content and codon preferences were similar. Mutation hotspot analysis identified four highly variable regions, which can potentially be used as molecular markers to identify Artemisia species. Phylogenetic analysis showed that the five Artemisia species investigated in this study were sister branches to each other, and individuals of each species formed a monophyletic clade. This study shows that the cp genome can provide distinguishing features to help identify closely related Artemisia species and has the potential to serve as a universal super barcode for plant identification.

Does IR-loss promote plastome structural variation and sequence evolution?

Plastids are one of the main distinguishing characteristics of the plant cell. The plastid genome (plastome) of most autotrophic seed plants possesses a highly conserved quadripartite structure containing a large single-copy (LSC) and a small single-copy (SSC) region separated by two copies of the inverted repeat (termed as IR_A and IR_B). The IRs have been inferred to stabilize the plastid genome via homologous recombination-induced repair mechanisms. IR loss has been documented in seven autotrophic flowering plant lineages and two autotrophic gymnosperm lineages, and the plastomes of these species (with a few exceptions) are rearranged to a great extent. However, some plastomes containing normal IRs also show high structural variation. Therefore, the role of IRs in maintaining plastome stability is still controversial. In this study, we first integrated and compared genome structure and sequence evolution of representative plastomes of all nine reported IR-lacking lineages and those of their closest relative(s) with canonical inverted repeats (CRCIRs for short) to explore the role of the IR in maintaining plastome structural stability and sequence evolution. We found the plastomes of most IR-lacking lineages have experienced significant structural rearrangement, gene loss and duplication, accumulation of novel small repeats, and acceleration of synonymous substitution compared with those of their CRCIRs. However, the IR-lacking plastomes show similar structural variation and sequence evolution rate, and even less rearrangement distance, dispersed repeat number, tandem repeat number, indels frequency and GC3 content than those of IR-present plastomes with variation in Geraniaceae. We argue that IR loss is not a driver of these changes but is instead itself a consequence of other processes that more broadly shape both structural and sequence-level plastome evolution.

Comparative Analysis of the Complete Chloroplast Genomes of Eight Ficus Species and Insights into the Phylogenetic Relationships of Ficus

The genus Ficus is an evergreen plant, the most numerous species in the family Moraceae, and is often used as a food and pharmacy source. The phylogenetic relationships of the genus Ficus have been debated for many years due to the overlapping phenotypic characters and morphological similarities between the genera. In this study, the eight Ficus species (Ficus altissima, Ficus auriculata, Ficus benjamina, Ficus curtipes, Ficus heteromorpha, Ficus lyrata, Ficus microcarpa, and Ficus virens) complete chloroplast (cp) genomes were successfully sequenced and phylogenetic analyses were made with other Ficus species. The result showed that the eight Ficus cp genomes ranged from 160,333 bp (F. heteromorpha) to 160,772 bp (F. curtipes), with a typical quadripartite structure. It was found that the eight Ficus cp genomes had similar genome structures, containing 127 unique genes. The cp genomes of the eight Ficus species contained 89−104 SSR loci, which were dominated by mono-nucleotides repeats. Moreover, we identified eight hypervariable regions (trnS-GCU_trnG-UCC, trnT-GGU_psbD, trnV-UAC_trnM-CAU, clpP_psbB, ndhF_trnL-UAG, trnL-UAG_ccsA, ndhD_psaC, and ycf1). Phylogenetic analyses have shown that the subgenus Ficus and subgenus Synoecia exhibit close affinities and based on the results, we prefer to merge the subgenus Synoecia into the subgenus Ficus. At the same time, new insights into the subgeneric classification of the Ficus macrophylla were provided. Overall, these results provide useful data for further studies on the molecular identification, phylogeny, species identification and population genetics of speciation in the Ficus genus.

The complete chloroplast genome of Persicaria perfoliata and comparative analysis with Four Medicinal Plants of Polygonaceae

Polygonaceae is a large family of medicinal herbs that includes many species used as traditional Chinese medicine, such as Persicaria perfoliate. Here, we sequenced the complete cp genome of P. perfoliata using Illumina sequencing technology with the purpose to provide a method to facilitate accurate identification. After being annotated, the cp genome of P. perfoliata was compared with Fagopyrum tataricum, Persicaria chinensis, Fagopyrum dibotrys and Fallopia multiflora. The complete cp genome of P. perfoliata is 160,730 bp in length, containing a small single copy (SSC) region of 12,927 bp, a large single copy (LSC) region of 85,433 bp and a pair of inverted repeats (IR) regions of 62,370 bp. A total of 131 genes were annotated, including eight rRNA genes, 34 tRNA genes and 84 protein-coding genes. Forty-two simple sequence repeats and fifty-five repeat sequences were identified. Mutational hot spots analyses indicated that five genes (matK, ndhF, ccsA, cemA, rpl20) could be selected as candidates for molecular markers. Moreover, phylogenetic analysis showed that all the Polygonaceae species formed a monophyletic clade, and P. perfoliata showed the closest relationship with P. chinense. The study provides valuable molecular information to accurately identify P. perfoliata and assist in its development and application.

Fifteen complete chloroplast genomes of Trapa species (Trapaceae): insight into genome structure, comparative analysis and phylogenetic relationships

BACKGROUND

Trapa L. is a floating-leaved aquatic plant with important economic and ecological values. However, the species identification and phylogenetic relationship within Trapa are still controversial, which necessitates the need for plastid genome information of Trapa. In this study, complete chloroplast genomes of 13 Trapa species/taxa were sequenced and annotated. Combined with released sequences, comparative analyses of chloroplast genomes were performed on the 15 Trapa species/taxa for the first time.

RESULTS

The Trapa chloroplast genomes exhibited typical quadripartite structures with lengths from 155,453 to 155,559 bp. The gene orders and contents within Trapa were conservative, but several changes were found in the microstructure. The intron loss of rpl2, also detected in Lythraceae, was found in all Trapa species/taxa, suggesting close genetic relationship between Lythraceae and Trapaceae. Notably, two small-seed species (T. incisa and T. maximowiczii) showed the smallest genome size with 155,453 and 155,477 bp, respectively. Each cp genome contained the same 130 genes consisting of 85 protein-coding genes, 37 tRNA genes and 8 rRNA genes. Trapa species/taxa showed 37 (T. incisa and T. maximowiczii) to 41 (T. sibirica) long repeats, including forward, palindromic, reversed and complementary repeats. There were 110 (T. quadrispinosa) to 123 (T. incisa and T. maximowiczii) SSR (simple sequence repeat) loci in Trapa chloroplast genomes. Comparative analyses revealed that two hotspot regions (atpA-atpF and rps2-rpoC2) in Trapa chloroplast genomes could be served as potential molecular markers. Three phylogenetic analyses (ML, MP and BI) consistently showed that there were two clusters within Trapa, including large- and small-seed species/taxa, respectively; for the large-seed Trapa, they clustered according to their geographical origin and tubercle morphology on the surface of seeds.

CONCLUSION

In summary, we have acquired the sequences of 13 Trapa chloroplast genomes, and performed the comparative analyses within Trapa for the first time. The results have helped us better identify the Trapa species/taxa and deepen the understanding of genetic basis and phylogenetic relationship of Trapa, which will facilitate the effective management and utilization of the important genetic resources in the future.

Comparative Analyses of 3,654 Plastid Genomes Unravel Insights Into Evolutionary Dynamics and Phylogenetic Discordance of Green Plants

The plastid organelle is essential for many vital cellular processes and the growth and development of plants. The availability of a large number of complete plastid genomes could be effectively utilized to understand the evolution of the plastid genomes and phylogenetic relationships among plants. We comprehensively analyzed the plastid genomes of Viridiplantae comprising 3,654 taxa from 298 families and 111 orders and compared the genomic organizations in their plastid genomic DNA among major clades, which include gene gain/loss, gene copy number, GC content, and gene blocks. We discovered that some important genes that exhibit similar functions likely formed gene blocks, such as the psb family presumably showing co-occurrence and forming gene blocks in Viridiplantae. The inverted repeats (IRs) in plastid genomes have doubled in size across land plants, and their GC content is substantially higher than non-IR genes. By employing three different data sets [all nucleotide positions (nt123), only the first and second codon positions (nt12), and amino acids (AA)], our phylogenomic analyses revealed Chlorokybales + Mesostigmatales as the earliest-branching lineage of streptophytes. Hornworts, mosses, and liverworts forming a monophylum were identified as the sister lineage of tracheophytes. Based on nt12 and AA data sets, monocots, Chloranthales and magnoliids are successive sister lineages to the eudicots + Ceratophyllales clade. The comprehensive taxon sampling and analysis of different data sets from plastid genomes recovered well-supported relationships of green plants, thereby contributing to resolving some long-standing uncertainties in the plant phylogeny.

Loss of the IR region in conifer plastomes: Changes in the selection pressure and substitution rate of protein-coding genes

Plastid genomes (plastomes) have a quadripartite structure, but some species have drastically reduced or lost inverted repeat (IR) regions. IR regions are important for genome stability and the evolution rate. In the evolutionary process of gymnosperms, the typical IRs of conifers were lost, possibly affecting the evolutionary rate and selection pressure of genomic protein-coding genes. In this study, we selected 78 gymnosperm species (51 genera, 13 families) for evolutionary analysis. The selection pressure analysis results showed that negative selection effects were detected in all 50 common genes. Among them, six genes in conifers had higher ω values than non-conifers, and 12 genes had lower ω values. The evolutionary rate analysis results showed that 9 of 50 common genes differed between conifers and non-conifers. It is more obvious that in non-conifers, the rates of psbA (trst, trsv, ratio, dN, dS, and ω) were 2.6- to 3.1-fold of conifers. In conifers, trsv, ratio, dN, dS, and ω of ycf2 were 1.2- to 3.6-fold of non-conifers. In addition, the evolution rate of ycf2 in the IR was significantly reduced. psbA is undergoing dynamic change, with an abnormally high evolution rate as a small portion of it enters the IR region. Although conifers have lost the typical IR regions, we detected no change in the substitution rate or selection pressure of most protein-coding genes due to gene function, plant habitat, or newly acquired IRs.

The dynamic history of gymnosperm plastomes: Insights from structural characterization, comparative analysis, phylogenomics, and time divergence

Gymnosperms are among the most endangered groups of plant species; they include ginkgo, pines (Conifers I), cupressophytes (Conifers II), cycads, and gnetophytes. The relationships among the five extant gymnosperm groups remain equivocal. We analyzed 167 available gymnosperm plastomes and investigated their diversity and phylogeny. We found that plastome size, structure, and gene order were highly variable in the five gymnosperm groups, of which Parasitaxus usta (Vieill.) de Laub. and Macrozamia mountperriensis F.M.Bailey had the smallest and largest plastomes, respectively. The inverted repeats (IRs) of the five groups were shown to have evolved through distinctive evolutionary scenarios. The IRs have been lost in all conifers but retained in cycads and gnetophytes. A positive association between simple sequence repeat (SSR) abundance and plastome size was observed, and the SSRs with the most variation were found in Pinaceae. Furthermore, the number of repeats was negatively correlated with IR length; thus, the highest number of repeats was detected in Conifers I and II, in which the IRs had been lost. We constructed a phylogeny based on 29 shared genes from 167 plastomes. With the plastome tree and 13 calibrations, we estimated the tree height between present-day angiosperms and gymnosperms to be ∼380 million years ago (mya). The placement of Gnetales in the tree agreed with the Gnetales-other gymnosperms hypothesis. The divergence between Ginkgo and cycads was estimated as ∼284 mya; the crown age of the cycads was 251 mya. Our time-calibrated plastid-based phylogenomic tree provides a framework for comparative studies of gymnosperm evolution.

Novel Structural Variation and Evolutionary Characteristics of Chloroplast tRNA in Gossypium Plants

Cotton is one of the most important fiber and oil crops in the world. Chloroplast genomes harbor their own genetic materials and are considered to be highly conserved. Transfer RNAs (tRNAs) act as "bridges" in protein synthesis by carrying amino acids. Currently, the variation and evolutionary characteristics of tRNAs in the cotton chloroplast genome are poorly understood. Here, we analyzed the structural variation and evolution of chloroplast tRNA (cp tRNA) based on eight diploid and two allotetraploid cotton species. We also investigated the nucleotide evolution of chloroplast genomes in cotton species. We found that cp tRNAs in cotton encoded 36 or 37 tRNAs, and 28 or 29 anti-codon types with lengths ranging from 60 to 93 nucleotides. Cotton chloroplast tRNA sequences possessed specific conservation and, in particular, the Ψ-loop contained the conserved U-U-C-X3-U. The cp tRNAs of Gossypium L. contained introns, and cp tRNA^Ile contained the anti-codon (C-A-U), which was generally the anti-codon of tRNA^Met. The transition and transversion analyses showed that cp tRNAs in cotton species were iso-acceptor specific and had undergone unequal rates of evolution. The intergenic region was more variable than coding regions, and non-synonymous mutations have been fixed in cotton cp genomes. On the other hand, phylogeny analyses indicated that cp tRNAs of cotton were derived from several inferred ancestors with greater gene duplications. This study provides new insights into the structural variation and evolution of chloroplast tRNAs in cotton plants. Our findings could contribute to understanding the detailed characteristics and evolutionary variation of the tRNA family.

Comparative transcriptome analysis infers bulb derived in vitro cultures as a promising source for sipeimine biosynthesis in Fritillaria cirrhosa D. Don (Liliaceae, syn. Fritillaria roylei Hook.) - High value Himalayan medicinal herb

Fritillaria cirrhosa D. Don (Liliaceae, syn. Fritillaria roylei Hook.) is a critically endangered medicinal herb of immense importance due to its pharmaceutical bioactive compound, especially sipeimine, used for the treatment of chronic respiratory disorders. However, the industrial demand for sipeimine solely depends on its endangered natural habitat. Therefore; there is an utmost need for its biodiversity conservation as well as for the sustainable utilization of phytochemicals. Plant cell culture and transcriptomics-based molecular bioprospection of key regulatory genes involved in sipeimine biosynthesis as such will play a crucial role in exploring the unexplored traits, that are in supply crisis or nearly in extinction stage. De novo comparative transcriptome sequencing of the bulb (in vivo), callus, and regenerated plantlets (in vitro) resulted in more than 150 million high-quality paired-end clean reads that assembled into final 31,428 transcripts. Functional annotation and unigenes classification with multiple public databases such as KEGG, Refseq, Uniprot, TAIR, GO, and COG, etc. along with chemical structures and functional biocatalytic activity analysis of different steroidal alkaloids facilitated the identification of 30 unigenes specific to sipeimine biosynthesis. Additionally, ABC transporters and TFs like bHLH, MYC, MYB, and WRKY suggests their possible role in metabolite translocation and regulation in vivo as well as in vitro tissues. Differential gene expression and quantitative analysis revealed that the MVA pathway probably the predominant route for 5C intermediate (IPP & DMAPP) biosynthesis. Further, the genes involved in the downstream biosynthesis pathway viz. SQLE, CAS1, SMT1, SMO1, SMO2, SC5DL, DHCR7, DHCR24, CYP710A, 3β-HSD, CYP90D2, and CYP374A6 shown similar expression pattern with RNA-Seq and qRT-PCR findings. The positive correlation between higher expression of proposed biosynthetic pathway genes and relatively higher accumulation of sipeimine in differentiated naturally grown bulb tissues (in vivo), undifferentiated cells (callus), and de-differentiated tissues i.e. regenerated plantlets (in vitro) has been evident from the present study. Comprehensive genomic resources created in F. cirrhosa will provide strong evidence of bulb derived in vitro culture as an alternative promising source for steroidal alkaloids biosynthesis and metabolite upscaling through genetic and metabolic engineering.

Comparative plastomic analysis and insights into the phylogeny of Salvia (Lamiaceae)

Salvia is the largest genus of Lamiaceae, with almost 1000 species, and has been divided into 11 subgenera. Salvia subg. Glutinaria, native to East Asia, is particularly important because of its potential medicinal value. However, the interspecific relationships of this subgenus have not been resolved and the plastomes of Salvia have rarely been studied. In the current study, we compared plastid genome structure and organization of 19 species of Salvia (14 newly sequenced and 5 previously published). Our comparative analysis showed that all Salvia plastomes examined have a quadripartite structure typical of most angiosperms and contain an identical set of 114 unique genes (80 protein-coding genes, 4 rRNA genes, and 30 tRNA genes). The plastome structure of all Salvia species is highly conserved like other Lamiaceae plastomes. Gene content, gene order, and GC content were highly similar in these plastomes. The inverted repeats/single copy region (IR/SC) boundaries of Salvia are highly conserved, and IR contraction only occurred in two species (Salvia mekongensis and S. rosmarinus). In Salvia, sequence divergence was higher in non-coding regions than in coding regions. We found that using large single copy (LSC) and small single copy regions (SSC) with exclusion of the rapidly evolving sites produced the highest resolution in phylogenetic analysis of Salvia, suggesting that using suitable informative sites to build trees is more conducive in phylogenetic research. This study assembled a powerful matrix data set for studying the phylogeny of Salvia, resolving the interspecific relationship of Salvia subg. Glutinaria. The newly sequenced plastid genomes will also enrich the plastome database of Salvia, providing the scientific basis for the development and utilization of germplasm resources of this large and important genus.

Uncovering dynamic evolution in the plastid genome of seven Ligusticum species provides insights into species discrimination and phylogenetic implications

Ligusticum L., one of the largest members in Apiaceae, encompasses medicinally important plants, the taxonomic statuses of which have been proved to be difficult to resolve. In the current study, the complete chloroplast genomes of seven crucial plants of the best-known herbs in Ligusticum were presented. The seven genomes ranged from 148,275 to 148,564 bp in length with a highly conserved gene content, gene order and genomic arrangement. A shared dramatic decrease in genome size resulted from a lineage-specific inverted repeat (IR) contraction, which could potentially be a promising diagnostic character for taxonomic investigation of Ligusticum, was discovered, without affecting the synonymous rate. Although a higher variability was uncovered in hotspot divergence regions that were unevenly distributed across the chloroplast genome, a concatenated strategy for rapid species identification was proposed because separate fragments inadequately provided variation for fine resolution. Phylogenetic inference using plastid genome-scale data produced a concordant topology receiving a robust support value, which revealed that L. chuanxiong had a closer relationship with L. jeholense than L. sinense, and L. sinense cv. Fuxiong had a closer relationship to L. sinense than L. chuanxiong, for the first time. Our results not only furnish concrete evidence for clarifying Ligusticum taxonomy but also provide a solid foundation for further pharmaphylogenetic investigation.

Chloroplast Genome Sequences and Comparative Analyses of Combretaceae Mangroves with Related Species

In the Combretaceae family, only two species of Lumnitzera and one species of Laguncularia belong to mangroves. Among them, Lumnitzera littorea (Jack) Voigt. is an endangered mangrove plant in China for the limited occurrence and seed abortion. In contrast, Lumnitzera racemosa Willd. is known as the most widespread mangrove plant in China. Laguncularia racemosa C. F. Gaertn., an exotic mangrove in China, has the fast growth and high adaptation ability. To better understand the phylogenetic positions of these mangroves in Combretaceae and in Myrtales and to provide information for studies on evolutionary adaptation for intertidal habitat, the complete chloroplast (cp) genomes of Lu. racemosa and La. racemosa were sequenced. Furthermore, we present here the results from the assembly and annotation of the two cp genomes, which were further subjected to the comparative analysis with Lu. littorea cp genomes we published before and other eleven closely related species within Myrtales. The chloroplast genomes of the three Combretaceae mangrove species: Lu. littorea, Lu. racemosa, and La. racemosa are 159,687 bp, 159,473 bp, and 158,311 bp in size. All three cp genomes host 130 genes including 85 protein-coding genes, 37 tRNAs, and 4 rRNAs. A comparative analysis of those three genomes revealed the high similarity of genes in coding-regions and conserved gene order in the IR and LSC/SSC regions. The differences between Lumnitzera and Laguncularia cp genomes are the locations of rps19 and rpl2 genes in the IR/SC boundary regions. Investigating the effects of selection events on shared protein-coding genes showed a relaxed selection had acted on the ycf2, ycf1, and matK genes of Combretaceae mangroves compared to the nonmangrove species Eucalyptus aromaphloia. The phylogenetic analysis based on the whole chloroplast genome sequence with one outgroup species strongly supported three Combretaceae mangroves together with other two Combretaceae species formed a cluster in Combretaceae. This study is the first report on the comparative analysis of three Combretaceae mangrove chloroplast genomes, which will provide the significant information for understanding photosynthesis and evolution in Combretaceae mangrove plants.

Complete Chloroplast Genomes and Comparative Analyses of L. chinensis, L. anhuiensis, and L. aurea (Amaryllidaceae)

The genus Lycoris (about 20 species) includes important medicinal and ornamental plants. Due to the similar morphological features and insufficient genomic resources, germplasm identification and molecular phylogeny analysis are very limited. Here, we sequenced the complete chloroplast genomes of L. chinensis, L. anhuiensis, and L. aurea; they have very similar morphological traits that make it difficult to identify. The full length of their cp genomes was nearly 158k bp with the same guanine-cytosine content of 37.8%. A total of 137 genes were annotated, including 87 protein-coding genes, 42 tRNAs, and eight rRNAs. A comparative analysis revealed the conservation in sequence size, GC content, and gene content. Some variations were observed in repeat structures, gene expansion on the IR-SC (Inverted Repeat-Single-Copy) boundary regions. Together with the cpSSR (chloroplast simple sequence repeats), these genetic variations are useful to develop molecular markers for germplasm identification. Phylogenetic analysis showed that seven Lycoris species were clustered into a monophyletic group, and closed to Narcissus in Amaryllidaceae. L. chinensis, L. anhuiensis, and L. longituba were clustered together, suggesting that they were very likely to be derived from one species, and had the same ancestor with L. squamigera. Our results provided information on the study of genetic diversity, origins or relatedness of native species, and the identification of cultivars.

The characteristics of the chloroplast genome of the Michelia martini

The chloroplast genome sequence of Michelia martini was sequenced using high-throughput sequencing technology and analyzed phylogenetically in the present study. The complete chloroplast genome was 159,819 bp in length, including a large single copy region (LSC) of 88,078 bp and small single copy region (SSC) of 18,801 bp, and a pair of inverted repeat regions (IR) of 26,470 bp. The contents of CG in the chloroplast genome were 39.29%. The sequence contained 128 unique genes, including 81 protein-coding genes, 38 tRNA genes and 8 rRNA genes. The phylogenetic analysis revealed that M. martini is closely related to Magnolia maudiae.

The complete chloroplast genome sequence of Michelia compressa

Michelia compressa is an evergreen ornamental tree species. The high-throughput sequencing technology was used to sequence and assemble the chloroplast genome of Michelia compressa. Results showed that the chloroplast genome is 160,061 bp in length, of which the inverted repeats sequence (IRs) is 26,581 bp, the large single-copy region (LSC) and the small single copy region (SSC) are 88,097 bp and 18,802 bp, respectively. The GC content of the plastome was 39.2%, with 43.2%, 37.9% and 34.2% in IRs, LSC and SSC, respectively. A total of 132 genes are annotated, 86 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. This study enriched the Michelia compressa genomic information which provides the basis for rational exploitation and utilization of germplasm resources.

Plastome Structural Conservation and Evolution in the Clusioid Clade of Malpighiales

The clusioid clade of Malpighiales is comprised of five families: Bonnetiaceae, Calophyllaceae, Clusiaceae, Hypericaceae and Podostemaceae. Recent studies have found the plastome structure of Garcinia mangostana L. from Clusiaceae was conserved, while plastomes of five riverweed species from Podostemaceae showed significant structural variations. The diversification pattern of plastome structure of the clusioid clade worth a thorough investigation. Here we determined five complete plastomes representing four families of the clusioid clade. Our results found that the plastomes of the early diverged three families (Clusiaceae, Bonnetiaceae and Calophyllaceae) in the clusioid clade are relatively conserved, while the plastomes of the other two families show significant variations. The Inverted Repeat (IR) regions of Tristicha trifaria and Marathrum foeniculaceum (Podostemaceae) are greatly reduced following the loss of the ycf1 and ycf2 genes. An inversion over 50 kb spanning from trnK-UUU to rbcL in the LSC region is shared by Cratoxylum cochinchinense (Hypericaceae), T. trifaria and Ma. foeniculaceum (Podostemaceae). The large inversed colinear block in Hypericaceae and Podostemaceae contains all the genes in the 50-kb inversed colinear block in a clade of Papilionoideae, with two extra genes (trnK-UUU and matK) at one end. Another endpoint of both inversions in the two clusioids families and Papilionoideae is located between rbcL and accD. This study greatly helped to clarify the plastome evolution in the clusioid clade.

Expanded inverted repeat region with large scale inversion in the first complete plastid genome sequence of Plantago ovata

Plantago ovata (Plantaginaceae) is an economically and medicinally important species, however, least is known about its genomics and evolution. Here, we report the first complete plastome genome of P. ovata and comparison with previously published genomes of related species from Plantaginaceae. The results revealed that P. ovata plastome size was 162,116 bp and that it had typical quadripartite structure containing a large single copy region of 82,084 bp and small single copy region of 5,272 bp. The genome has a markedly higher inverted repeat (IR) size of 37.4 kb, suggesting large-scale inversion of 13.8 kb within the expanded IR regions. In addition, the P. ovata plastome contains 149 different genes, including 43 tRNA, 8 rRNA, and 98 protein-coding genes. The analysis revealed 139 microsatellites, of which 71 were in the non-coding regions. Approximately 32 forward, 34 tandem, and 17 palindromic repeats were detected. The complete genome sequences, 72 shared genes, matK gene, and rbcL gene from related species generated the same phylogenetic signals, and phylogenetic analysis revealed that P. ovata formed a single clade with P. maritima and P. media. The divergence time estimation as employed in BEAST revealed that P. ovata diverged from P. maritima and P. media about 11.0 million years ago (Mya; 95% highest posterior density, 10.06-12.25 Mya). In conclusion, P. ovata had significant variation in the IR region, suggesting a more stable P. ovata plastome genome than that of other Plantaginaceae species.

Comparative Genomic Analysis Reveals the Mechanism Driving the Diversification of Plastomic Structure in Taxaceae Species

Inverted repeat (IR) regions in the plastomes from land plants induce homologous recombination, generating isomeric plastomes. While the plastomes of Taxaceae species often lose one of the IR regions, considerable isomeric plastomes were created in Taxaceae species with a hitherto unclarified mechanism. To investigate the detailed mechanism underpinning the IR-independent genesis of plastomic diversity, we sequenced four Taxaceae plastomes, including Taxus cuspidata Siebold & Zuccarini, Taxus fauna Nan Li & R. R. Mill, and two individuals of Taxus wallichiana Zuccarini. Then we compared these structures with those of previously reported Taxaceae plastomes. Our analysis identified four distinct plastome forms that originated from the rearrangements of two IR-flanking inverted fragments. The presence of isomeric plastomes was then verified in T. cuspidata individuals. Both rearrangement analyses and phylogenetic results indicated that Taxaceae were separated into two clades, one including Taxus and Pseudotaxus and another formed by Amentotaxus and Torreya. Our reconstructed scenario suggests that the minimum number of inversion events required for the transformation of the plastome of Cephalotaxus oliveri Masters into the diversified Taxaceae plastomes ranged from three to six. To sum up, our study reveals a distinct pattern and the mechanism driving the structural diversification of Taxaceae plastomes, which will advance our understanding of the maintenance of plastomic diversity and complexity in conifers.

The complete chloroplast genome sequence of Manglietia yuyuanensis (Magnoliaceae)

Manglietia yuyuanensis is an important afforestation and excellent broad-leaved tree species in southern China. In this study, we assembled the complete chloroplast genome of M. yuyuanensis based on the Illumina sequences, sequence analysis showed the genome was 160,078 bp in length presenting a typical quadripartite structure and contains an inverted repeat region (IR, 26,467 b), a small single-copy (SSC) region, and a large single-copy (LSC) region (18,785 and 88,359 bp, respectively). The overall GC content was 39.27%. The sequence contained 128 unique genes, including 81 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. The maximum-likelihood (ML) phylogenetic analysis revealed that M. maudiae was closely related to Manglietia insignis. A phylogenetic analysis revealed that M. yuyuanensis is closely related to Manglietia glaucifolia, with the genus Manglietia.

Resequencing 545 ginkgo genomes across the world reveals the evolutionary history of the living fossil

As Charles Darwin anticipated, living fossils provide excellent opportunities to study evolutionary questions related to extinction, competition, and adaptation. Ginkgo (Ginkgo biloba L.) is one of the oldest living plants and a fascinating example of how people have saved a species from extinction and assisted its resurgence. By resequencing 545 genomes of ginkgo trees sampled from 51 populations across the world, we identify three refugia in China and detect multiple cycles of population expansion and reduction along with glacial admixture between relict populations in the southwestern and southern refugia. We demonstrate multiple anthropogenic introductions of ginkgo from eastern China into different continents. Further analyses reveal bioclimatic variables that have affected the geographic distribution of ginkgo and the role of natural selection in ginkgo’s adaptation and resilience. These investigations provide insights into the evolutionary history of ginkgo trees and valuable genomic resources for further addressing various questions involving living fossil species.

Ginkgo is one of the living fossils from the plant kingdom. Here, authors conduct population genomics analyses to reveal its refugia and demographic history, and provide evidence of multiple anthropogenic introductions of ginkgo from eastern China into different continents.

Intergeneric Relationships within the Family Salicaceae s.l. based on Plastid Phylogenomics

Many Salicaceae s.l. plants are recognized for their important role in the production of products such as wood, oils, and medicines, and as a model organism in life studies. However, the difference in plastid sequence, phylogenetic relationships, and lineage diversification of the family Salicaceae s.l. remain poorly understood. In this study, we compare 24 species representing 18 genera of the family. Simple sequence repeats (SSRs) are considered effective molecular markers for plant species identification and population genetics. Among them, a total of 1798 SSRs were identified, among which mononucleotide repeat was the most common with 1455 accounts representing 80.92% of the total. Most of the SSRs are located in the non-coding region. We also identified five other types of repeats, including 1750 tandems, 434 forward, 407 palindromic, 86 reverse, and 30 complementary repeats. The species in Salicaceae s.l. have a conserved plastid genome. Each plastome presented a typical quadripartite structure and varied in size due to the expansion and contraction of the inverted repeat (IR) boundary, lacking major structural variations, but we identified six divergence hotspot regions. We obtained phylogenetic relationships of 18 genera in Salicaceae s.l. and the 24 species formed a highly supported lineage. Casearia was identified as the basal clade. The divergence time between Salicaceae s.l. and the outgroup was estimated as ~93 Mya; Salix, and Populus diverged around 34 Mya, consistent with the previously reported time. Our research will contribute to a better understanding of the phylogenetic relationships among the members of the Salicaceae s.l.

Complete Chloroplast Genome of Pinus densiflora Siebold & Zucc. and Comparative Analysis with Five Pine Trees

Pinus densiflora (Korean red pine) is widely distributed in East Asia and considered one of the most important species in Korea. In this study, the complete chloroplast genome of P. densiflora was sequenced by combining the advantages of Oxford Nanopore MinION and Illumina MiSeq. The sequenced genome was then compared with that of a previously published conifer plastome. The chloroplast genome was found to be circular and comprised of a quadripartite structure, including 113 genes encoding 73 proteins, 36 tRNAs and 4 rRNAs. It had short inverted repeat regions and lacked ndh gene family genes, which is consistent with other Pinaceae species. The gene content of P. densiflora was found to be most similar to that of P. sylvestris. The newly attempted sequencing method could be considered an alternative method for obtaining accurate genetic information, and the chloroplast genome sequence of P. densiflora revealed in this study can be used in the phylogenetic analysis of Pinus species.

Characterization of the complete chloroplast genome of Michelia maudiae (Magnoliaceae)

Michelia maudiae is an evergreen ornamental shrub with strong aromatic flowers with great important for horticulture. In this study, the complete chloroplast genome of M. maudiae was assembled based on the Illumina sequences. The genome was 160,154 bp in length presenting a typical quadripartite structure and contains an inverted repeat region (IR, 26,575 bp), a large single copy region (LSC, 88,200 bp) and a small single copy region (SSC, 18,804 bp). The overall GC content was 39.24%. A total of 130 genes were identified, including 83 protein-coding genes, 38 transfer RNA genes, and eight rRNA genes. Twelve gene species contained one or two introns. The maximum-likelihood phylogenetic analysis revealed that M. maudiae was closely related with Manglietia insignis.

The complete chloroplast genome sequence of Cynanchum forrestii Schltr. (Asclepiadaceae) and its phylogenetic analysis

Cynanchum forrestii is a folk medicinal plant in southwest China. In this study, we sequenced complete chloroplast (cp) genome sequence of C. forrestii to investigate its phylogenetic relationship. The whole cp genome of C. forrestii was 159,917 bp in length with 43.5% overall GC content, including a large single-copy (LSC) region of 91,189 bp and a small single-copy (SSC) region of 19,972 bp, which were separated by a pair of inverted repeats (IRs) of 24,378 bp. The cp genome contained 112 genes, including 78 protein-coding genes, 30 tRNA genes, and 4 rRNA genes. The phylogenetic analysis based on cp genome sequences showed that Cynanchum was closely related with Asclepias and Calotropis.

The newly developed single nucleotide polymorphism (SNP) markers for a potentially medicinal plant, Crepidiastrum denticulatum (Asteraceae), inferred from complete chloroplast genome data

Medicinal effects of Crepidiastrum denticulatum have been previously reported. However, the genomic resources of this species and its applications have not been studied. In this study, based on the next generation sequencing method (Miseq sequencing system), we characterize the chloroplast genome of C. denticulatum which contains a large single copy (84,112 bp) and a small single copy (18,519 bp), separated by two inverted repeat regions (25,074 bp). This genome consists of 80 protein-coding gene, 30 tRNAs, and four rRNAs. Notably, the trnT_GGU is pseudogenized because of a small insertion within the coding region. Comparative genomic analysis reveals a high similarity among Asteraceae taxa. However, the junctions between LSC, SSC, and IRs locate in different positions within rps19 and ycf1 among examined species. Also, we describe a newly developed single nucleotide polymorphism (SNP) marker for C. denticulatum based on amplification-refractory mutation system (ARMS) technique. The markers, inferred from SNP in rbcL and matK genes, show effectiveness to recognize C. denticulatum from other related taxa through simple PCR protocol. The chloroplast genome-based molecular markers are effective to distinguish a potentially medicinal species, C. denticulatum, from other related taxa. Additionally, the complete chloroplast genome of C. denticulatum provides initial genomic data for further studies on phylogenomics, population genetics, and evolutionary history of Crepidiastrum as well as other taxa in Asteraceae.

Identification and Phylogenetic Analysis of the Complete Chloroplast Genomes of Three Ephedra Herbs Containing Ephedrine

Ephedrae Herba and Ephedrae Radix et Rhizoma (Mahuang) have been used as Chinese herbal medicines. Ephedra plants mainly live in deserts and have good governance of desertification. Despite their important medicinal and environmental protection value, dietary supplements containing ephedrine from Ephedra species may threaten the health of people. Morphological resemblance amongst species causes difficulty in identifying the original species of Ephedra herbs. Chloroplast (CP) genome shows good prospects in identification and phylogenetic analysis. This study introduced the structures of the CP genomes of three Ephedra species and analysed their phylogenetic relationships. Three complete CP genomes of Ephedra showed four-part annular structures, namely, two single-copy regions and two inverted repeat regions. The entire CP genomes of three Ephedra species in terms of size were 109,550 bp (E. sinica), 109,667 bp (E. intermedia), and 109,558 bp (E. equisetina). Each CP genome of the three Ephedra species encoded 118 genes, including 73 protein-coding genes, 37 tRNA genes and 8 ribosomal RNA genes. Eleven high-variation regions were screened through mVISTA to be potential specific DNA barcodes for identifying Ephedra species. Maximum likelihood and maximum parsimony trees showed that CP genomes could be used to identify Ephedra species. The Ephedra species had a close phylogenetic relationship with Gnetum species and Welwitschia mirabilis. This research provided valuable information for the identification and phylogenetic analysis of gymnosperms and drug safety of Ephedra.

Mobile Elements Shape Plastome Evolution in Ferns

Plastid genomes display remarkable organizational stability over evolutionary time. From green algae to angiosperms, most plastid genomes are largely collinear, with only a few cases of inversion, gene loss, or, in extremely rare cases, gene addition. These plastome insertions are mostly clade-specific and are typically of nuclear or mitochondrial origin. Here, we expand on these findings and present the first family-level survey of plastome evolution in ferns, revealing a novel suite of dynamic mobile elements. Comparative plastome analyses of the Pteridaceae expose several mobile open reading frames that vary in sequence length, insertion site, and configuration among sampled taxa. Even between close relatives, the presence and location of these elements is widely variable when viewed in a phylogenetic context. We characterize these elements and refer to them collectively as Mobile Open Reading Frames in Fern Organelles (MORFFO). We further note that the presence of MORFFO is not restricted to Pteridaceae, but is found across ferns and other plant clades. MORFFO elements are regularly associated with inversions, intergenic expansions, and changes to the inverted repeats. They likewise appear to be present in mitochondrial and nuclear genomes of ferns, indicating that they can move between genomic compartments with relative ease. The origins and functions of these mobile elements are unknown, but MORFFO appears to be a major driver of structural genome evolution in the plastomes of ferns, and possibly other groups of plants.

Reconfiguration of the plastid genome in Lamprocapnos spectabilis: IR boundary shifting, inversion, and intraspecific variation

We generated a complete plastid genome (plastome) sequence for Lamprocapnos spectabilis, providing the first complete plastome from the subfamily Fumarioideae (Papaveraceae). The Lamprocapnos plastome shows large differences in size, structure, gene content, and substitution rates compared with two sequenced Papaveraceae plastomes. We propose a model that explains the major rearrangements observed, involving at least six inverted repeat (IR) boundary shifts and five inversions, generating a number of gene duplications and relocations, as well as a two-fold expansion of the IR and miniaturized small single-copy region. A reduction in the substitution rates for genes transferred from the single-copy regions to the IR was observed. Accelerated substitution rates of plastid accD and clpP were detected in the Lamprocapnos plastome. The accelerated substitution rate for the accD gene was correlated with a large insertion of amino acid repeat (AAR) motifs in the middle region, but the forces driving the higher substitution rate of the clpP gene are unclear. We found a variable number of AARs in Lamprocapnos accD and ycf1 genes within individuals, and the repeats were associated with coiled-coil regions. In addition, comparative analysis of three Papaveraceae plastomes revealed loss of rps15 in Papaver, and functional replacement to the nucleus was identified.

Plastome characteristics of Cannabaceae

Cannabaceae is an economically important family that includes ten genera and ca. 117 accepted species. To explore the structure and size variation of their plastomes, we sequenced ten plastomes representing all ten genera of Cannabaceae. Each plastome possessed the typical angiosperm quadripartite structure and contained a total of 128 genes. The Inverted Repeat (IR) regions in five plastomes had experienced small expansions (330-983 bp) into the Large Single-Copy (LSC) region. The plastome of Chaetachme aristata has experienced a 942-bp IR contraction and lost rpl22 and rps19 in its IRs. The substitution rates of rps19 and rpl22 decreased after they shifted from the LSC to IR. A 270-bp inversion was detected in the Parasponia rugosa plastome, which might have been mediated by 18-bp inverted repeats. Repeat sequences, simple sequence repeats, and nucleotide substitution rates varied among these plastomes. Molecular markers with more than 13% variable sites and 5% parsimony-informative sites were identified, which may be useful for further phylogenetic analysis and species identification. Our results show strong support for a sister relationship between Gironniera and Lozanell (BS = 100). Celtis, Cannabis-Humulus, Chaetachme-Pteroceltis, and Trema-Parasponia formed a strongly supported clade, and their relationships were well resolved with strong support (BS = 100). The availability of these ten plastomes provides valuable genetic information for accurately identifying species, clarifying taxonomy and reconstructing the intergeneric phylogeny of Cannabaceae.

Ginkgo biloba's footprint of dynamic Pleistocene history dates back only 390,000 years ago

BACKGROUND

At the end of the Pliocene and the beginning of Pleistocene glaciation and deglaciation cycles Ginkgo biloba went extinct all over the world, and only few populations remained in China in relict areas serving as sanctuary for Tertiary relict trees. Yet the status of these regions as refuge areas with naturally existing populations has been proven not earlier than one decade ago. Herein we elaborated the hypothesis that during the Pleistocene cooling periods G. biloba expanded its distribution range in China repeatedly. Whole plastid genomes were sequenced, assembled and annotated, and sequence data was analyzed in a phylogenetic framework of the entire gymnosperms to establish a robust spatio-temporal framework for gymnosperms and in particular for G. biloba Pleistocene evolutionary history.

RESULTS

Using a phylogenetic approach, we identified that Ginkgoatae stem group age is about 325 million years, whereas crown group radiation of extant Ginkgo started not earlier than 390,000 years ago. During repeated warming phases, Gingko populations were separated and isolated by contraction of distribution range and retreated into mountainous regions serving as refuge for warm-temperate deciduous forests. Diversification and phylogenetic splits correlate with the onset of cooling phases when Ginkgo expanded its distribution range and gene pools merged.

CONCLUSIONS

Analysis of whole plastid genome sequence data representing the entire spatio-temporal genetic variation of wild extant Ginkgo populations revealed the deepest temporal footprint dating back to approximately 390,000 years ago. Present-day directional West-East admixture of genetic diversity is shown to be the result of pronounced effects of the last cooling period. Our evolutionary framework will serve as a conceptual roadmap for forthcoming genomic sequence data, which can then provide deep insights into the demographic history of Ginkgo.

Complete Chloroplast Genome of Pinus massoniana (Pinaceae): Gene Rearrangements, Loss of ndh Genes, and Short Inverted Repeats Contraction, Expansion

The chloroplast genome (CPG) of Pinus massoniana belonging to the genus Pinus (Pinaceae), which is a primary source of turpentine, was sequenced and analyzed in terms of gene rearrangements, ndh genes loss, and the contraction and expansion of short inverted repeats (IRs). P. massoniana CPG has a typical quadripartite structure that includes large single copy (LSC) (65,563 bp), small single copy (SSC) (53,230 bp) and two IRs (IRa and IRb, 485 bp). The 108 unique genes were identified, including 73 protein-coding genes, 31 tRNAs, and 4 rRNAs. Most of the 81 simple sequence repeats (SSRs) identified in CPG were mononucleotides motifs of A/T types and located in non-coding regions. Comparisons with related species revealed an inversion (21,556 bp) in the LSC region; P. massoniana CPG lacks all 11 intact ndh genes (four ndh genes lost completely; the five remained truncated as pseudogenes; and the other two ndh genes remain as pseudogenes because of short insertions or deletions). A pair of short IRs was found instead of large IRs, and size variations among pine species were observed, which resulted from short insertions or deletions and non-synchronized variations between "IRa" and "IRb". The results of phylogenetic analyses based on whole CPG sequences of 16 conifers indicated that the whole CPG sequences could be used as a powerful tool in phylogenetic analyses.

accD nuclear transfer of Platycodon grandiflorum and the plastid of early Campanulaceae

BACKGROUND

Campanulaceae species are known to have highly rearranged plastid genomes lacking the acetyl-CoA carboxylase (ACC) subunit D gene (accD), and instead have a nuclear (nr)-accD. Plastid genome information has been thought to depend on studies concerning Trachelium caeruleum and genome announcements for Adenophora remotiflora, Campanula takesimana, and Hanabusaya asiatica. RNA editing information for plastid genes is currently unavailable for Campanulaceae. To understand plastid genome evolution in Campanulaceae, we have sequenced and characterized the chloroplast (cp) genome and nr-accD of Platycodon grandiflorum, a basal member of Campanulaceae.

RESULTS

We sequenced the 171,818 bp cp genome containing a 79,061 bp large single-copy (LSC) region, a 42,433 bp inverted repeat (IR) and a 7840 bp small single-copy (SSC) region, which represents the cp genome with the largest IR among species of Campanulaceae. The genome contains 110 genes and 18 introns, comprising 77 protein-coding genes, four RNA genes, 29 tRNA genes, 17 group II introns, and one group I intron. RNA editing of genes was detected in 18 sites of 14 protein-coding genes. Platycodon has an IR containing a 3' rps12 operon, which occurs in the middle of the LSC region in four other species of Campanulaceae (T. caeruleum, A. remotiflora, C. takesimana, and H. asiatica), but lacks accD, clpP, infA, and rpl23, as has been found in these four species. Platycodon nr-accD contains about 3.2 kb intron between nr-accD.e1 and nr-accD.e2 at the same insertion point as in other Campanulaceae. The phylogenies of the plastid genomes and accD show that Platycodon is basal in the Campanulaceae clade, indicating that IR disruption in Campanulaceae occurred after the loss of accD, clpP, infA, and rpl23 in the cp genome, which occurred during plastid evolution in Campanulaceae.

CONCLUSIONS

The plastid genome of P. grandiflorum lacks the rearrangement of the IR found in T. caeruleum, A. remotiflora, C. takesimana, and H. asiatica. The absence of accD, clpP, infA, and rpl23 in the plastid genome is a synapomorphic characteristic of Campanulaceae. The chloroplast genome phylogeny supports the hypothesis that chloroplast genomic arrangement occurred after accD nuclear transfer and loss of the four genes in the plastid of early Campanulaceae as a lineage of asterids.

The complete chloroplast genome sequence of the medicinal and economic plant woad Isatis indigotica (Brassicaceae)

Chloroplast (cp) genome sequences provide a valuable source for phylogenetic analysis. It becomes a popular useful tool for population and phylogeny in recently report. Here, the complete chloroplast genome of the Isatis indigotica has been reconstructed from the whole-genome Illumina sequencing data. The circular genome is 156,670 bp in size and comprises a pair of inverted repeat (IR) regions of 26,995 bp each, a large single-copy (LSC) region of 84,907 bp and a small single-copy (SSC) region of 17,773 bp. The total GC content is 36.5%, while the corresponding values of the LSC, SSC and IR region are 34.2%, 29.7% and 42.3%, respectively. The chloroplast genome contains 140 genes, including 94 protein-coding genes. The phylogenetic analysis showed a strong sister relationship with Raphanus sativus in Brassicaceae. These findings provide a foundation for further investigation of cp genome evolution in I. indigotica and other higher plants.

Expansion of inverted repeat does not decrease substitution rates in Pelargonium plastid genomes

For species with minor inverted repeat (IR) boundary changes in the plastid genome (plastome), nucleotide substitution rates were previously shown to be lower in the IR than the single copy regions (SC). However, the impact of large-scale IR expansion/contraction on plastid nucleotide substitution rates among closely related species remains unclear. We included plastomes from 22 Pelargonium species, including eight newly sequenced genomes, and used both pairwise and model-based comparisons to investigate the impact of the IR on sequence evolution in plastids. Ten types of plastome organization with different inversions or IR boundary changes were identified in Pelargonium. Inclusion in the IR was not sufficient to explain the variation of nucleotide substitution rates. Instead, the rate heterogeneity in Pelargonium plastomes was a mixture of locus-specific, lineage-specific and IR-dependent effects. Our study of Pelargonium plastomes that vary in IR length and gene content demonstrates that the evolutionary consequences of retaining these repeats are more complicated than previously suggested.

Abundant RNA editing sites of chloroplast protein-coding genes in Ginkgo biloba and an evolutionary pattern analysis

BACKGROUND

RNA editing is a posttranscriptional modification process that alters the RNA sequence so that it deviates from the genomic DNA sequence. RNA editing mainly occurs in chloroplasts and mitochondrial genomes, and the number of editing sites varies in terrestrial plants. Why and how RNA editing systems evolved remains a mystery. Ginkgo biloba is one of the oldest seed plants and has an important evolutionary position. Determining the patterns and distribution of RNA editing in the ancient plant provides insights into the evolutionary trend of RNA editing, and helping us to further understand their biological significance.

RESULTS

In this paper, we investigated 82 protein-coding genes in the chloroplast genome of G. biloba and identified 255 editing sites, which is the highest number of RNA editing events reported in a gymnosperm. All of the editing sites were C-to-U conversions, which mainly occurred in the second codon position, biased towards to the U_A context, and caused an increase in hydrophobic amino acids. RNA editing could change the secondary structures of 82 proteins, and create or eliminate a transmembrane region in five proteins as determined in silico. Finally, the evolutionary tendencies of RNA editing in different gene groups were estimated using the nonsynonymous-synonymous substitution rate selection mode.

CONCLUSIONS

The G. biloba chloroplast genome possesses the highest number of RNA editing events reported so far in a seed plant. Most of the RNA editing sites can restore amino acid conservation, increase hydrophobicity, and even influence protein structures. Similar purifying selections constitute the dominant evolutionary force at the editing sites of essential genes, such as the psa, some psb and pet groups, and a positive selection occurred in the editing sites of nonessential genes, such as most ndh and a few psb genes.

Genes Translocated into the Plastid Inverted Repeat Show Decelerated Substitution Rates and Elevated GC Content

Plant chloroplast genomes (plastomes) are characterized by an inverted repeat (IR) region and two larger single copy (SC) regions. Patterns of molecular evolution in the IR and SC regions differ, most notably by a reduced rate of nucleotide substitution in the IR compared to the SC region. In addition, the organization and structure of plastomes is fluid, and rearrangements through time have repeatedly shuffled genes into and out of the IR, providing recurrent natural experiments on how chloroplast genome structure can impact rates and patterns of molecular evolution. Here we examine four loci (psbA, ycf2, rps7, and rps12 exon 2-3) that were translocated from the SC into the IR during fern evolution. We use a model-based method, within a phylogenetic context, to test for substitution rate shifts. All four loci show a significant, 2- to 3-fold deceleration in their substitution rate following translocation into the IR, a phenomenon not observed in any other, nontranslocated plastid genes. Also, we show that after translocation, the GC content of the third codon position and of the noncoding regions is significantly increased, implying that gene conversion within the IR is GC-biased. Taken together, our results suggest that the IR region not only reduces substitution rates, but also impacts nucleotide composition. This finding highlights a potential vulnerability of correlating substitution rate heterogeneity with organismal life history traits without knowledge of the underlying genome structure.

Comparison of intraspecific, interspecific and intergeneric chloroplast diversity in Cycads

Cycads are among the most threatened plant species. Increasing the availability of genomic information by adding whole chloroplast data is a fundamental step in supporting phylogenetic studies and conservation efforts. Here, we assemble a dataset encompassing three taxonomic levels in cycads, including ten genera, three species in the genus Cycas and two individuals of C. debaoensis. Repeated sequences, SSRs and variations of the chloroplast were analyzed at the intraspecific, interspecific and intergeneric scale, and using our sequence data, we reconstruct a phylogenomic tree for cycads. The chloroplast was 162,094 bp in length, with 133 genes annotated, including 87 protein-coding, 37 tRNA and 8 rRNA genes. We found 7 repeated sequences and 39 SSRs. Seven loci showed promising levels of variations for application in DNA-barcoding. The chloroplast phylogeny confirmed the division of Cycadales in two suborders, each of them being monophyletic, revealing a contradiction with the current family circumscription and its evolution. Finally, 10 intraspecific SNPs were found. Our results showed that despite the extremely restricted distribution range of C. debaoensis, using complete chloroplast data is useful not only in intraspecific studies, but also to improve our understanding of cycad evolution and in defining conservation strategies for this emblematic group.

Five Complete Chloroplast Genome Sequences from Diospyros: Genome Organization and Comparative Analysis

Diospyros is the largest genus in Ebenaceae, comprising more than 500 species with remarkable economic value, especially Diospyros kaki Thunb., which has traditionally been an important food resource in China, Korea, and Japan. Complete chloroplast (cp) genomes from D. kaki, D. lotus L., D. oleifera Cheng., D. glaucifolia Metc., and Diospyros 'Jinzaoshi' were sequenced using Illumina sequencing technology. This is the first cp genome reported in Ebenaceae. The cp genome sequences of Diospyros ranged from 157,300 to 157,784 bp in length, presenting a typical quadripartite structure with two inverted repeats each separated by one large and one small single-copy region. For each cp genome, 134 genes were annotated, including 80 protein-coding, 31 tRNA, and 4 rRNA unique genes. In all, 179 repeats and 283 single sequence repeats were identified. Four hypervariable regions, namely, intergenic region of trnQ_rps16, trnV_ndhC, and psbD_trnT, and intron of ndhA, were identified in the Diospyros genomes. Phylogenetic analyses based on the whole cp genome, protein-coding, and intergenic and intron sequences indicated that D. oleifera is closely related to D. kaki and could be used as a model plant for future research on D. kaki; to our knowledge, this is proposed for the first time. Further, these analyses together with two large deletions (301 and 140 bp) in the cp genome of D. 'Jinzaoshi', support its placement as a new species in Diospyros. Both maximum parsimony and likelihood analyses for 19 taxa indicated the basal position of Ericales in asterids and suggested that Ebenaceae is monophyletic in Ericales.

The complete plastid genome sequence of Picea jezoensis (Pinaceae: Piceoideae)

The nucleotide sequence of the complete chloroplast genome of P. jezoensis was completed. The total genome size was 124 146 bp, containing a pair of very short inverted repeats (IRa and IRb) of 422 bp, which were separated by large single copy (LSC) and small single copy (SSC) with 66 956 bp and 56 346 bp, respectively. The overall GC contents of the plastid genome were determined as 38.8%. One hundred fifteen genes including 68 peptide-encoding genes, 35 tRNA genes, four rRNA genes, six open-reading frames, and two pseudogenes were annotated. In these genes, 15 genes contained only one or two introns. Phylogenetic analyses using maximum likelihood (ML) methods were performed from fully sequenced Gymnosperms and other species of dataset composed of 69 protein-coding genes.

YCF1: A Green TIC: Response to the de Vries et al. Commentary

This response to a recent Commentary article by de Vries et al. highlights critical errors in the annotation and identification of Ycf1 homologs in the sequenced chloroplast genomes. Contrary to what is reported by de Vries et al., the majority of chloroplast genomes sequenced to date appear to have retained a typical Ycf1 sequence (i.e., including the N-terminal 6TM domain and a variable hydrophilic C-terminal domain) as my group previously reported. Our evidence continues to support the model that Ycf1 forms an essential component of a "green TIC" that is largely conserved among the Chlorophyta and land plants. Since the establishment of this green TIC with Tic20 as the core component, some cases of loss of Ycf1 during the evolution of the green lineages might be regarded as modifications or alterations of the complex. Here, I discuss our working model that the presence of an alternative "nonphotosynthetic-type" or "ancestral-type" TIC might explain other (or specific) cases of the lack of Ycf1, not only in early lineages, including Glaucophyta and Rhodophyta, but also in the grasses.