Unparalleled GC content in the plastid DNA of Selaginella

Phylogenomic data resolved the deep relationships of Gymnogynoideae (Selaginellaceae)

The unresolved phylogenetic framework within the Selaginellaceae subfamily Gymnogynoideae (ca. 130 species) has hindered our comprehension of the diversification and evolution of Selaginellaceae, one of the most important lineages in land plant evolution. Here, based on plastid and nuclear data extracted from genomic sequencing of more than 90% species of all genera except two in Gymnogynoideae, a phylogenomic study focusing on the contentious relationships among the genera in Gymnogynoideae was conducted. Our major results included the following: (1) Only single-copy region (named NR) and only one ribosomal operon was firstly found in Afroselaginella among vascular plants, the plastome structure of Gymnogynoideae is diverse among the six genera, and the direct repeats (DR) type is inferred as the ancestral state in the subfamily; (2) The first strong evidence was found to support Afroselaginella as a sister to Megaloselaginella. Alternative placements of Ericetorum and Gymnogynum were detected, and their relationships were investigated by analyzing the variation of phylogenetic signals; and (3) The most likely genus-level relationships in Gymnogynoideae might be: ((Bryodesma, Lepidoselaginella), (((Megaloselaginella, Afroselaginella), Ericetorum), Gymnogynum)), which was supported by maximum likelihood phylogeny based on plastid datasets, maximum likelihood, and Bayesian inference based on SCG dataset and concatenated nuclear and plastid datasets and the highest proportion of phylogenetic signals of plastid genes.

Chloroplast genome structure analysis of Equisetum unveils phylogenetic relationships to ferns and mutational hotspot region

Equisetum is one of the oldest extant group vascular plants and is considered to be the key to understanding vascular plant evolution. Equisetum is distributed almost all over the world and has a high degree of adaptability to different environments. Despite the fossil record of horsetails (Equisetum, Equisetaceae) dating back to the Carboniferous, the phylogenetic relationship of this genus is not well, and the chloroplast evolution in Equisetum remains poorly understood. In order to fill this gap, we sequenced, assembled, and annotated the chloroplast genomes of 12 species of Equisetum, and compared them to 13 previously published vascular plants chloroplast genomes to deeply examine the plastome evolutionary dynamics of Equisetum. The chloroplast genomes have a highly conserved quadripartite structure across the genus, but these chloroplast genomes have a lower GC content than other ferns. The size of Equisetum plastomes ranges from 130,773 bp to 133,684 bp and they encode 130 genes. Contraction/expansion of IR regions and the number of simple sequences repeat regions underlie large genomic variations in size among them. Comparative analysis revealed we also identified 13 divergence hotspot regions. Additionally, the genes accD and ycf1 can be used as potential DNA barcodes for the identification and phylogeny of the genus Equisetum. Twelve photosynthesis-related genes were specifically selected in Equisetum. Comparative genomic analyses implied divergent evolutionary patterns between Equisetum and other ferns. Phylogenomic analyses and molecular dating revealed a relatively distant phylogenetic relationship between Equisetum and other ferns, supporting the division of pteridophyte into Lycophytes, Equisetaceae and ferns. The results show that the chloroplast genome can be used to solve phylogenetic problems within or between Equisetum species, and also provide genomic resources for the study of Equisetum systematics and evolution.

The complete chloroplast genome of an Arctic lycophyte, Huperzia arctica (Tolm.) Sipliv. 1973 (Lycopodiaceae)

Huperzia arctica (Tolm.) Sipliv. 1973 is a lycophyte species belonging to the Lycopodiaceae family, which is widely distributed in the Arctic region of Svalbard, Norway. To determine its taxonomic position, we sequenced the chloroplast genome of H. arctica. The complete chloroplast genome of H. arctica is 153,956 bp in length with 122 annotated genes, including 87 protein-coding genes, 31 tRNA genes, and 4 rRNA genes. To evaluate its evolutionary position, we performed phylogenetic analysis using 36 conserved orthologous protein-coding gene sequences from the chloroplast genomes of H. arctica and publicly available data from other Lycopodiaceae. H. arctica formed a monophyletic group with four other Huperzia species: H. lucidula, H. serrata, H. crispata, and H. javanica. However, it appeared as a separate species with a highly supported bootstrap value.

Assembly, annotation and analysis of the chloroplast genome of the Algarrobo tree Neltuma pallida (subfamily: Caesalpinioideae)

BACKGROUND

Neltuma pallida is a tree that grows in arid soils in northwestern Peru. As a predominant species of the Equatorial Dry Forest ecoregion, it holds significant economic and ecological value for both people and environment. Despite this, the species is severely threatened and there is a lack of genetic and genomic research, hindering the proposal of evidence-based conservation strategies.

RESULTS

In this work, we conducted the assembly, annotation, analysis and comparison of the chloroplast genome of a N. pallida specimen with those of related species. The assembled chloroplast genome has a length of 162,381 bp with a typical quadripartite structure (LSC-IRA-SSC-IRB). The calculated GC content was 35.97%. However, this is variable between regions, with a higher GC content observed in the IRs. A total of 132 genes were annotated, of which 19 were duplicates and 22 contained at least one intron in their sequence. A substantial number of repetitive sequences of different types were identified in the assembled genome, predominantly tandem repeats (> 300). In particular, 142 microsatellites (SSR) markers were identified. The phylogenetic reconstruction showed that N. pallida grouped with the other Neltuma species and with Prosopis cineraria. The analysis of sequence divergence between the chloroplast genome sequences of N. pallida, N. juliflora, P. farcta and Strombocarpa tamarugo revealed a high degree of similarity.

CONCLUSIONS

The N. pallida chloroplast genome was found to be similar to those of closely related species. With a size of 162,831 bp, it had the classical chloroplast quadripartite structure and GC content of 35.97%. Most of the 132 identified genes were protein-coding genes. Additionally, over 800 repetitive sequences were identified, including 142 SSR markers. In the phylogenetic analysis, N. pallida grouped with other Neltuma spp. and P. cineraria. Furthermore, N. pallida chloroplast was highly conserved when compared with genomes of closely related species. These findings can be of great potential for further diversity studies and genetic improvement of N. pallida.

Phylogeny, character evolution, and classification of Selaginellaceae (lycophytes)

Selaginella is the largest and most taxonomically complex genus in lycophytes. The fact that over 750 species are currently treated in a single genus makes Selaginellales/Selaginellaceae unique in pteridophytes. Here we assembled a dataset of six existing and newly sampled plastid and nuclear loci with a total of 684 accessions (74% increase of the earlier largest sampling) representing ca. 300 species to infer a new phylogeny. The evolution of 10 morphological characters is studied in the new phylogenetic context. Our major results include: (1) the nuclear and plastid phylogenies are congruent with each other and combined analysis well resolved and strongly supported the relationships of all but two major clades; (2) the Sinensis group is resolved as sister to S. subg. Pulviniella with strong support in two of the three analyses; (3) most morphological characters are highly homoplasious but some characters alone or combinations of characters well define the major clades in the family; and (4) an infrafamilial classification of Selaginellaceae is proposed and the currently defined Selaginella s.l. is split into seven subfamilies (corresponding to the current six subgenera + the Sinensis group) and 19 genera (the major diagnosable clades) with nine new species-poor genera. We support the conservation of Selaginella with a new type, S. flabellata, to minimize nomenclatural instability. We provide a key to subfamilies and genera, images illustrating their morphology, their morphological and geographical synopses, a list of constituent species, and necessary new combinations. This new classification will hopefully facilitate communication, promote further studies, and help conservation.

IRplus: An Augmented Tool to Detect Inverted Repeats in Plastid Genomes

High-throughput sequencing methods have increased the accessibility of plastid genomes, which are crucial for clarifying phylogenetic relationships. Current large sequencing efforts require software tools for routine display of their distinctive quadripartite structure, which is denoted by four junction sites. By concentrating on these junctions and their close vicinity, IRscope has emerged as the standard tool for detection of this structure and creating simplified comparative graphical maps of plastid genomes. Here, we provide an augmented version (IRplus) that encompasses a novel set of functions such as integrated error detection, flexible color schemes, and an upgraded method to detect inverted repeats in genomic sequences. Spanning across the plant tree of life, IRplus allows the quick visualization of various sets of plastid genomes and features, next to smooth interoperability with other widely used annotation file formats and platforms. The IRplus can be accessed at https://irscope.shinyapps.io/IRplus/, and source codes are freely available at https://github.com/AmiryousefiLab/IRplus.

Mitochondrial RNA editing sites affect the phylogenetic reconstruction of gymnosperms

•RNA editing sites may contain homoplasious signals that cause artifactual inferences in phylogenetic analyses.•Excluding RNA editing sites from gymnosperm mitochondrial genes restored the sister relationship of gnetophytes and Pinaceae.•Phylogenetic analysis based on mitochondrial genomic data should carefully evaluate the impact of RNA editing sites.

Evolution of 101 Apocynaceae plastomes and phylogenetic implications

Apocynaceae are one of the ten species-richest angiosperm families. However, the backbone phylogeny of the family is yet less well supported, and the evolution of plastome structure has not been thoroughly studied for the whole family. Herein, a total of 101 complete plastomes including 35 newly sequenced, 24 reassembled from public raw data and the rest from the NCBI GenBank database, representing 26 of 27 tribes of Apocynaceae, were used for comparative plastome analysis. Phylogenetic analyses were conducted using a combined plastid data matrix of 77 protein-coding genes from 162 taxa, encompassing all tribes and 41 of 49 subtribes of Apocynaceae. Plastome lengths ranged from 150,897 bp in Apocynum venetum to 178,616 bp in Hoya exilis. Six types of boundaries between the inverted repeat (IR) regions and single copy (SC) regions were identified. Different sizes of IR expansion were found in three lineages, including Alyxieae, Ceropegieae and Marsdenieae, suggesting multiple expansion events of the IRs over the SC regions in Apocynaceae. The IR regions of Marsdenieae evolved in two ways: expansion towards the large single copy (LSC) region in Lygisma + Stephanotis + Ruehssia + Gymnema (Cosmopolitan clade), and expansion towards both LSC and small single copy (SSC) region in Dischidia-Hoya alliance and Marsdenia (Asia-Pacific clade). Six coding genes and five non-coding regions were identified as highly variable, including accD, ccsA-ndhD, clpP, matK, ndhF, ndhG-ndhI, trnG(GCC)-trnfM(CAU), trnH(GUG)-psbA, trnY(GUA)-trnE(UUC), ycf1, and ycf2. Maximum likelihood and Bayesian phylogenetic analyses resulted in nearly identical tree topologies and produced a well-resolved backbone comprising 15 consecutive dichotomies that subdivided Apocynaceae into 15 clades. The subfamily Periplocoideae were embedded in the Apocynoid grade and were sister to the Echiteae-Odontadenieae-Mesechiteae clade with high support values. Three tribes (Melodineae, Vinceae, and Willughbeieae), the subtribe Amphineuriinae, and four genera (Beaumontia, Ceropegia, Hoya, and Stephanotis) were not resolved as monophyletic. Our work sheds light on the backbone phylogenetic relationships in the family Apocynaceae and offers insights into the evolution of Apocynaceae plastomes using the most densely sampled plastome dataset to date.

Mitogenome-based phylogenomics provides insights into the positions of the enigmatic sinensis group and the sanguinolenta group in Selaginellaceae (Lycophyte)

Spikemoss (Selaginellaceae) is one of the basal lineages of vascular plants. This family has a single genus Selaginella which consists of about 750 extant species. The phylogeny of Selaginellaceae has been extensively studied mainly based on plastid DNA and a few nuclear sequences. However, the placement of the enigmatic sinensis group is a long-term controversy because of the long branch in the plastid DNA phylogeny. The sanguinolenta group is also a phylogenetically problematic clade owing to two alternative positions resulted from different datasets. Here, we newly sequenced 34 mitochondrial genomes (mitogenomes) of individuals representing all seven subgenera and major clades in Selaginellaceae. We assembled the draft mitogenomes and annotated the genes and performed phylogenetic analyses based on the shared 17 mitochondrial genes. Our major results include: (1) all the assembled mitogenomes have complicated structures, unparalleled high GC content and a small gene content set, and the positive correlations among GC content, substitution rates and the number of RNA editing sites hold; (2) the sinensis group was well supported as a member of subg. Stachygynandrum; (3) the sanguinolenta group was strongly resolved as sister to all other Selaginella species except for subg. Selaginella. This study demonstrates the potential of mitogenome data in providing novel insights into phylogenetically recalcitrant problems.

Emergence of Novel RNA-Editing Sites by Changes in the Binding Affinity of a Conserved PPR Protein

RNA editing converts cytidines to uridines in plant organellar transcripts. Editing typically restores codons for conserved amino acids. During evolution, specific C-to-U editing sites can be lost from some plant lineages by genomic C-to-T mutations. By contrast, the emergence of novel editing sites is less well documented. Editing sites are recognized by pentatricopeptide repeat (PPR) proteins with high specificity. RNA recognition by PPR proteins is partially predictable, but prediction is often inadequate for PPRs involved in RNA editing. Here we have characterized evolution and recognition of a recently gained editing site. We demonstrate that changes in the RNA recognition motifs that are not explainable with the current PPR code allow an ancient PPR protein, QED1, to uniquely target the ndhB-291 site in Brassicaceae. When expressed in tobacco, the Arabidopsis QED1 edits 33 high-confident off-target sites in chloroplasts and mitochondria causing a spectrum of mutant phenotypes. By manipulating the relative expression levels of QED1 and ndhB-291, we show that the target specificity of the PPR protein depends on the RNA:protein ratio. Finally, our data suggest that the low expression levels of PPR proteins are necessary to ensure the specificity of editing site selection and prevent deleterious off-target editing.

A Comprehensive Evolutionary Study of Chloroplast RNA Editing in Gymnosperms: A Novel Type of G-to-A RNA Editing Is Common in Gymnosperms

Although more than 9100 plant plastomes have been sequenced, RNA editing sites of the whole plastome have been experimentally verified in only approximately 21 species, which seriously hampers the comprehensive evolutionary study of chloroplast RNA editing. We investigated the evolutionary pattern of chloroplast RNA editing sites in 19 species from all 13 families of gymnosperms based on a combination of genomic and transcriptomic data. We found that the chloroplast C-to-U RNA editing sites of gymnosperms shared many common characteristics with those of other land plants, but also exhibited many unique characteristics. In contrast to that noted in angiosperms, the density of RNA editing sites in ndh genes was not the highest in the sampled gymnosperms, and both loss and gain events at editing sites occurred frequently during the evolution of gymnosperms. In addition, GC content and plastomic size were positively correlated with the number of chloroplast RNA editing sites in gymnosperms, suggesting that the increase in GC content could provide more materials for RNA editing and facilitate the evolution of RNA editing in land plants or vice versa. Interestingly, novel G-to-A RNA editing events were commonly found in all sampled gymnosperm species, and G-to-A RNA editing exhibits many different characteristics from C-to-U RNA editing in gymnosperms. This study revealed a comprehensive evolutionary scenario for chloroplast RNA editing sites in gymnosperms, and reported that a novel type of G-to-A RNA editing is prevalent in gymnosperms.

Complete chloroplast genomes and comparative analyses of Hippeastrum ‘milady’, Hippeastrum albertii and Hippeastrum reticulatum (Amaryllidaceae)

Hippeastrum is a genus of ornamental plants with large, brightly colored flowers. Due to the very high seed-setting rate of the hybridization of Hippeastrum, the large population of hybrid progeny and the existence of superparent inheritance, it is difficult to trace the origin of the varieties collected from the market during breeding. In this study, we analyzed the chloroplast genomes of Hippeastrum ‘Milady’, H. alberti, and H. reticulatum using the Illumina NovaSeq sequencing platform and generated full-length sequences of 158,067, 158,067, and 158,522 bp, respectively. All three genomes had the typical tetrad structure. The large single copy, small single copy, and inverted repeat regions of H. reticulatum were observed to be respectively 277, 138, and 20 bp longer than the corresponding regions of H. ‘Milady’ and H. alberti. The results of comparative analysis of simple sequence repeats (SSRs), Ka/Ks ratios, codon preferences, and complete sequences of chloroplasts of these three taxa and 14 other plant species were as follows. First, the chloroplast genomes of H. ‘Milady’, H. alberti, and H. reticulatum contain 209, 209, and 211 SSR sites, respectively, most of which (123, 123, and 122, respectively) are single nucleotide repeats. Second, leucine, arginine, and serine are the most frequently used amino acids in the three chloroplast genomes. Third, H. ‘Milady’, H. alberti, and H. reticulatum are more closely related to Lycoris and Narcissus than to Allium and Agapanthus. Our results will provide information on the study of origins or relatedness of native species, and the identification of cultivars.

The evolution of extremely diverged plastomes in Selaginellaceae (lycophyte) is driven by repeat patterns and the underlying DNA maintenance machinery

Two factors are proposed to account for the unusual features of organellar genomes: the disruptions of organelle-targeted DNA replication, repair, and recombination (DNA-RRR) systems in the nuclear genome and repetitive elements in organellar genomes. Little is known about how these factors affect organellar genome evolution. The deep-branching vascular plant family Selaginellaceae is known to have a deficient DNA-RRR system and convergently evolved organellar genomes. However, we found that the plastid genome (plastome) of Selaginella sinensis has extremely accelerated substitution rates, a low GC content, pervasive repeat elements, a dynamic network structure, and it lacks direct or inverted repeats. Unexpectedly, its organelle DNA-RRR system is short of a plastid-targeted Recombinase A1 (RecA1) and a mitochondrion-targeted RecA3, in line with other explored Selaginella species. The plastome contains a large collection of short- and medium-sized repeats. Given the absence of RecA1 surveillance, we propose that these repeats trigger illegitimate recombination, accelerated mutation rates, and structural instability. The correlations between repeat quantity and architectural complexity in the Selaginella plastomes support these conclusions. We, therefore, hypothesize that the interplay of the deficient DNA-RRR system and the high repeat content has led to the extraordinary divergence of the S. sinensis plastome. Our study not only sheds new light on the mechanism of plastome divergence by emphasizing the power of cytonuclear integration, but it also reconciles the longstanding contradiction on the effects of DNA-RRR system disruption on genome structure evolution.

Plastid Phylogenomics and Plastomic Diversity of the Extant Lycophytes

Although extant lycophytes represent the most ancient surviving lineage of early vascular plants, their plastomic diversity has long been neglected. The ancient evolutionary history and distinct genetic diversity patterns of the three lycophyte families, each with its own characteristics, provide an ideal opportunity to investigate the interfamilial relationships of lycophytes and their associated patterns of evolution. To compensate for the lack of data on Lycopodiaceae, we sequenced and assembled 14 new plastid genomes (plastomes). Combined with other lycophyte plastomes available online, we reconstructed the phylogenetic relationships of the extant lycophytes based on 93 plastomes. We analyzed, traced, and compared the plastomic diversity and divergence of the three lycophyte families (Isoëtaceae, Lycopodiaceae, and Selaginellaceae) in terms of plastomic diversity by comparing their plastome sizes, GC contents, substitution rates, structural rearrangements, divergence times, ancestral states, RNA editings, and gene losses. Comparative analysis of plastid phylogenomics and plastomic diversity of three lycophyte families will set a foundation for further studies in biology and evolution in lycophytes and therefore in vascular plants.

Plastid phylogenomic analyses of the Selaginella sanguinolenta group (Selaginellaceae) reveal conflict signatures resulting from sequence types, outlier genes, and pervasive RNA editing

Different from the generally conserved plastomes (plastid genomes) of most land plants, the Selaginellaceae plastomes exhibit dynamic structure, high GC content and high substitution rates. Previous plastome analyses identified strong conflict on several clades in Selaginella, however the factors causing the conflictions and the impact on the phylogenetic inference have not been sufficiently investigated. Here, we dissect the distribution of phylogenetic signals and conflicts in Selaginella sanguinolenta group, the plastome of which is DR (direct repeats) structure and with genome-wide RNA editing. We analyzed the data sets including 22 plastomes representing all species of the S. sanguinolenta group, covering the entire geographical distribution from the Himalayas to Siberia and the Russian Far East regions. We recovered four different topologies by applying multispecies coalescent (ASTRAL) and concatenation methods (IQ-TREE and RAxML) on four data sets of PC (protein-coding genes), NC (non-coding sequences), PCN (the concatenated PC and NC), and RC (predicted RNA editing sites "C" were corrected by "T"), respectively. Six monophyletic clades, S. nummularifolia clade, S. rossii clade, S. sajanensis clade, S. sanguinolenta I clade, S. sanguinolenta II clade, and S. sanguinolenta III clade, were consistently resolved and supported by the characteristics of GC content, RNA editing frequency, and gene content. However, the relationships among these clades varied across the four topologies. To explore the underlying causes of the uncertainty, we compared the phylogenetic signals of the four topologies. We identified that the sequence types (coding versus non-coding), outlier genes (genes with extremely high |ΔGLS| values), and C-to-U RNA editing frequency in the protein-coding genes were responsible for the unstable phylogenomic relationship. We further revealed a significant positive correlation between the |ΔGLS| values and the variation coefficient of the RNA editing number. Our results demonstrated that the coalescent method performed better than the concatenation method in overcoming the problems caused by outlier genes and extreme RNA editing events. Our study particularly focused on the importance of exploring the plastid phylogenomic conflicts and suggested conducting concatenated analyses cautiously when adopting organelle genome data.

Anatomical and Biochemical Traits Related to Blue Leaf Coloration of Selaginella uncinata

Selaginella uncinata shows particularly rare blue leaves. Previous research has shown that structural interference by the cell wall of adaxial epidermal cells imparts blue coloration in leaves of S. uncinata; the objective of this study was to see whether anthocyanins might additionally contribute to this color, as changes in pH, and conjugation with metals and other flavonoids is also known to result in blue coloration in plants. We compared anatomical and biochemical traits of shade-grown (blue) S. uncinata leaves to high light (red) leaves of the same species and also to a non-blue (green) leaves of a congeneric S. kraussiana. By examining the anatomical structure, we found that the shape of adaxial epidermis of S. uncinata leaves was convex or lens-shaped on the lateral view and irregular circles with smooth embossment on the top view. These features were different from those of the abaxial and adaxial epidermis of S. kraussiana. We suspect that these structures increase the proportion of incident light entering the cell, deepening the leaf color, and therefore may be related to blue leaf color in S. uncinata. By examining biochemical traits, we found little difference in leaf pH value among the leaf types; all leaves contained several metal ions such as Mg, Fe, Mn, and copigments such as flavones. However, because there was no anthocyanin in blue S. uncinata leaves, we concluded that blue coloration in S. uncinata leaves is not caused by the three hypotheses of blue coloration: alkalization of the vacuole pH, metal chelation, or copigmentation with anthocyanins, but it may be related to the shape of the leaf adaxial epidermis.

Plastome structure, evolution, and phylogeny of Selaginella

As one of the earliest land plant lineages, Selaginella is important for studying land plant evolution. It is the largest genus of lycophytes containing 700-800 species. Some unique characters of Selaginella plastomes have been reported, but based only on 20 species. There have been no plastome phylogenies of Selaginella based on a relatively large sampling, and no efforts have been made to resolve the phylogeny of the enigmatic Sinensis group whose relationships have been unclear based on small datasets. Here we investigated the structures of 59 plastomes representing 51 species covering all six subgenera and 18 sections of Selaginella except two sections and including the intriguing Sinensis group for the first time. Our major results include: (1) the plastome size of Selaginella ranges tremendously from 78,492 bp to 187,632 bp; (2) there are numerous gene losses in Selaginella comparing with other lycophytes, Isoëtaceae and Lycopodiaceae; (3) the gene contents and plastome structures in Selaginella vary lineage-specifically and all infrageneric taxa are well supported in the plastome phylogeny; (4) the ndh gene family tends to lose or pseudogenize in those species with DR structure and without other short or medium repeats; (5) the short and medium repeat regions in SC mediate many conformations causing diverse and complex plastome structures, and six new conformations are discovered; (6) forty-eight species sampled have high GC content (>50%) but three species in the Sinensis group have ∼30% GC content in plastomes, similar to most vascular plants; (7) the Sinensis group is monophyletic, includes at least two subgroups, and has the smallest plastomes in land plants except some parasitic plants, and their plastomes do not contain any tRNAs; (8) the younger lineages in Selaginella tend to have higher GC content, whereas the older lineages tend to have lower GC content; and (9) because of incomplete genomic data and abnormal structures or some unknown reasons, even the concatenated plastomes could not well resolve the phylogenetic relationships in Selaginella with confidence, highlighting the difficulty in resolving the phylogeny and evolution of this particularly important land plant lineage.

Plastid phylogenomics of the Gynoxoid group (Senecioneae, Asteraceae) highlights the importance of motif-based sequence alignment amid low genetic distances

PREMISE

The genus Gynoxys and relatives form a species-rich lineage of Andean shrubs and trees with low genetic distances within the sunflower subtribe Tussilaginineae. Previous molecular phylogenetic investigations of the Tussilaginineae have included few, if any, representatives of this Gynoxoid group or reconstructed ambiguous patterns of relationships for it.

METHODS

We sequenced complete plastid genomes of 21 species of the Gynoxoid group and related Tussilaginineae and conducted detailed comparisons of the phylogenetic relationships supported by the gene, intron, and intergenic spacer partitions of these genomes. We also evaluated the impact of manual, motif-based adjustments of automatic DNA sequence alignments on phylogenetic tree inference.

RESULTS

Our results indicate that the inclusion of all plastid genome partitions is needed to infer well-supported phylogenetic trees of the Gynoxoid group. Whole plastome-based tree inference suggests that the genera Gynoxys and Nordenstamia are polyphyletic and form the core clade of the Gynoxoid group. This clade is sister to a clade of Aequatorium and Paragynoxys and also includes some but not all representatives of Paracalia.

CONCLUSIONS

The concatenation and combined analysis of all plastid genome partitions and the construction of manually-curated, motif-based DNA sequence alignments are found to be instrumental in the recovery of well-supported relationships of the Gynoxoid group. We demonstrate that the correct assessment of homology in genome-level plastid sequence data sets is crucial for subsequent phylogeny reconstruction and that the manual post-processing of multiple sequence alignments improves the reliability of such reconstructions amid low genetic distances between taxa.

Structural variation and evolution of chloroplast tRNAs in green algae

As one of the important groups of the core Chlorophyta (Green algae), Chlorophyceae plays an important role in the evolution of plants. As a carrier of amino acids, tRNA plays an indispensable role in life activities. However, the structural variation of chloroplast tRNA and its evolutionary characteristics in Chlorophyta species have not been well studied. In this study, we analyzed the chloroplast genome tRNAs of 14 species in five categories in the green algae. We found that the number of chloroplasts tRNAs of Chlorophyceae is maintained between 28-32, and the length of the gene sequence ranges from 71 nt to 91 nt. There are 23-27 anticodon types of tRNAs, and some tRNAs have missing anticodons that are compensated for by other types of anticodons of that tRNA. In addition, three tRNAs were found to contain introns in the anti-codon loop of the tRNA, but the analysis scored poorly and it is presumed that these introns are not functional. After multiple sequence alignment, the Ψ-loop is the most conserved structural unit in the tRNA secondary structure, containing mostly U-U-C-x-A-x-U conserved sequences. The number of transitions in tRNA is higher than the number of transversions. In the replication loss analysis, it was found that green algal chloroplast tRNAs may have undergone substantial gene loss during the course of evolution. Based on the constructed phylogenetic tree, mutations were found to accompany the evolution of the Green algae chloroplast tRNA. Moreover, chloroplast tRNAs of Chlorophyceae are consistent with those of monocotyledons and gymnosperms in terms of evolutionary patterns, sharing a common multi-phylogenetic pattern and rooted in a rich common ancestor. Sequence alignment and systematic analysis of tRNA in chloroplast genome of Chlorophyceae, clarified the characteristics and rules of tRNA changes, which will promote the evolutionary relationship of tRNA and the origin and evolution of chloroplast.

Integrated analysis of three newly sequenced fern chloroplast genomes: Genome structure and comparative analysis

BACKGROUND

Some ferns have medicinal properties and are used in therapeutic interventions. However, the classification and phylogenetic relationships of ferns remain incompletely reported. Considering that chloroplast genomes provide ideal information for species identification and evolution, in this study, three unpublished and one published ferns were sequenced and compared with other ferns to obtain comprehensive information on their classification and evolution.

MATERIALS AND METHODS

The complete chloroplast genomes of Dryopteris goeringiana (Kunze) Koidz, D. crassirhizoma Nakai, Athyrium brevifrons Nakai ex Kitagawa, and Polystichum tripteron (Kunze) Presl were sequenced using the Illumina HiSeq 4,000 platform. Simple sequence repeats (SSRs), nucleotide diversity analysis, and RNA editing were investigated in all four species. Genome comparison and inverted repeats (IR) boundary expansion and contraction analyses were also performed. The relationships among the ferns were studied by phylogenetic analysis based on the whole chloroplast genomes.

RESULTS

The whole chloroplast genomes ranged from 148,539 to 151,341 bp in size and exhibited typical quadripartite structures. Ten highly variable loci with parsimony informative (Pi) values of > 0.02 were identified. A total of 75-108 SSRs were identified, and only six SSRs were present in all four ferns. The SSRs contained a higher number of A + T than G + C bases. C-to-U conversion was the most common type of RNA editing event. Genome comparison analysis revealed that single-copy regions were more highly conserved than IR regions. IR boundary expansion and contraction varied among the four ferns. Phylogenetic analysis showed that species in the same genus tended to cluster together with and had relatively close relationships.

CONCLUSION

The results provide valuable information on fern chloroplast genomes that will be useful to identify and classify ferns, and study their phylogenetic relationships and evolution.

Complete chloroplast genomes shed light on phylogenetic relationships, divergence time, and biogeography of Allioideae (Amaryllidaceae)

Allioideae includes economically important bulb crops such as garlic, onion, leeks, and some ornamental plants in Amaryllidaceae. Here, we reported the complete chloroplast genome (cpDNA) sequences of 17 species of Allioideae, five of Amaryllidoideae, and one of Agapanthoideae. These cpDNA sequences represent 80 protein-coding, 30 tRNA, and four rRNA genes, and range from 151,808 to 159,998 bp in length. Loss and pseudogenization of multiple genes (i.e., rps2, infA, and rpl22) appear to have occurred multiple times during the evolution of Alloideae. Additionally, eight mutation hotspots, including rps15-ycf1, rps16-trnQ-UUG, petG-trnW-CCA, psbA upstream, rpl32-trnL-UAG, ycf1, rpl22, matK, and ndhF, were identified in the studied Allium species. Additionally, we present the first phylogenomic analysis among the four tribes of Allioideae based on 74 cpDNA coding regions of 21 species of Allioideae, five species of Amaryllidoideae, one species of Agapanthoideae, and five species representing selected members of Asparagales. Our molecular phylogenomic results strongly support the monophyly of Allioideae, which is sister to Amaryllioideae. Within Allioideae, Tulbaghieae was sister to Gilliesieae-Leucocoryneae whereas Allieae was sister to the clade of Tulbaghieae- Gilliesieae-Leucocoryneae. Molecular dating analyses revealed the crown age of Allioideae in the Eocene (40.1 mya) followed by differentiation of Allieae in the early Miocene (21.3 mya). The split of Gilliesieae from Leucocoryneae was estimated at 16.5 mya. Biogeographic reconstruction suggests an African origin for Allioideae and subsequent spread to Eurasia during the middle Eocene. Cool and arid conditions during the late Eocene led to isolation between African and Eurasian species. African Allioideae may have diverged to South American taxa in the late Oligocene. Rather than vicariance, long-distance dispersal is the most likely explanation for intercontinental distribution of African and South American Allioideae species.

Correlation Analysis Reveals an Important Role of GC Content in Accumulation of Deletion Mutations in the Coding Region of Angiosperm Plastomes

Variation in GC content is assumed to correlate with various processes, including mutation biases, recombination, and environmental parameters. To date, most genomic studies exploring the evolution of GC content have focused on nuclear genomes, but relatively few have concentrated on organelle genomes. We explored the mechanisms maintaining the GC content in angiosperm plastomes, with a particular focus on the hypothesis of phylogenetic dependence and the correlation with deletion mutations. We measured three genetic traits, namely, GC content, A/T tracts, and G/C tracts, in the coding region of plastid genomes for 1382 angiosperm species representing 350 families and 64 orders, and tested the phylogenetic signal. Then, we performed correlation analyses and revealed the variation in evolutionary rate of selected traits using RRphylo. The plastid GC content in the coding region varied from 28.10% to 43.20% across angiosperms, with a few non-photosynthetic species showing highly reduced values, highlighting the significance of functional constraints. We found strong phylogenetic signal in A/T tracts, but weak ones in GC content and G/C tracts, indicating adaptive potential. GC content was positively and negatively correlated with G/C and A/T tracts, respectively, suggesting a trade-off between these two deletion events. GC content evolved at various rates across the phylogeny, with significant increases in monocots and Lamiids, and a decrease in Fabids, implying the effects of some other factors. We hypothesize that variation in plastid GC content might be a mixed strategy of species to optimize fitness in fluctuating climates, partly through influencing the trade-off between AT → GC and GC → AT mutations.

Plastid Genomes of the Early Vascular Plant Genus Selaginella Have Unusual Direct Repeat Structures and Drastically Reduced Gene Numbers

The early vascular plants in the genus Selaginella, which is the sole genus of the Selaginellaceae family, have an important place in evolutionary history, along with ferns, as such plants are valuable resources for deciphering plant evolution. In this study, we sequenced and assembled the plastid genome (plastome) sequences of two Selaginella tamariscina individuals, as well as Selaginella stauntoniana and Selaginella involvens. Unlike the inverted repeat (IR) structures typically found in plant plastomes, Selaginella species had direct repeat (DR) structures, which were confirmed by Oxford Nanopore long-read sequence assembly. Comparative analyses of 19 lycophytes, including two Huperzia and one Isoetes species, revealed unique phylogenetic relationships between Selaginella species and related lycophytes, reflected by structural rearrangements involving two rounds of large inversions that resulted in dynamic changes between IR and DR blocks in the plastome sequence. Furthermore, we present other uncommon characteristics, including a small genome size, drastic reductions in gene and intron numbers, a high GC content, and extensive RNA editing. Although the 16 Selaginella species examined may not fully represent the genus, our findings suggest that Selaginella plastomes have undergone unique evolutionary events yielding genomic features unparalleled in other lycophytes, ferns, or seed plants.

Distinctive evolutionary pattern of organelle genomes linked to the nuclear genome in Selaginellaceae

Plastids and mitochondria are endosymbiotic organelles that store genetic information. The genomes of these organelles generally exhibit contrasting patterns regarding genome architecture and genetic content. However, they have similar genetic features in Selaginellaceae, and little is known about what causes parallel evolution. Here, we document the multipartite plastid genomes (plastomes) and the highly divergent mitochondrial genomes (mitogenomes) from spikemoss obtained by combining short- and long-reads. The 188-kb multipartite plastome has three ribosomal operon copies in the master genomic conformation, creating the alternative subgenomic conformation composed of 110- and 78-kb subgenomes. The long-read data indicated that the two different genomic conformations were present in almost equal proportions in the plastomes of Selaginella nipponica. The mitogenome of S. nipponica was assembled into 27 contigs with a total size of 110 kb. All contigs contained directly arranged repeats at both ends, which introduced multiple conformations. Our results showed that plastomes and mitogenomes share high tRNA losses, GC-biased nucleotides, elevated substitution rates and complicated organization. The exploration of nuclear-encoded organelle DNA replication, recombination and repair proteins indicated that, several single-targeted proteins, particularly plastid-targeted recombinase A1, have been lost in Selaginellaceae; conversely, the dual-targeted proteins remain intact. According to the reported function of recombinase A1, we propose that the plastomes of spikemoss often fail to pair homologous sequences during recombination, and the dual-targeted proteins play a key role in the convergent genetic features of plastomes and mitogenomes. Our results provide a distinctive evolutionary pattern of the organelle genomes in Selaginellaceae and evidence of their convergent evolution.

Episodic and guanine-cytosine-biased bursts of intragenomic and interspecific synonymous divergence in Ajugoideae (Lamiaceae) mitogenomes

Synonymous substitution rates in plant mitochondrial genomes vary by orders of magnitude among species, whereas synonymous rates among genes within a genome are generally consistent. Exceptionally, genes within the Ajuga reptans (Lamiaceae) mitochondrial genome exhibit unprecedented intragenomic heterogeneity in synonymous sequence divergence, but the biological mechanisms underlying this rate variation remain unclear. We tracked the origin and evolutionary trajectory of mitochondrial rate variations by dense sampling in Ajugoideae and found differences in the timing and magnitude of rate acceleration for particular genes. The most divergent genes accelerated earlier, retained a high rate across Ajugoideae, and are generally devoid of RNA editing, whereas moderately diverged genes accelerated later and retained relatively higher RNA editing frequency. The acceleration of mutation rates correlates with increased guanine-cytosine (GC) content, suggesting a key role for GC-biased gene conversion and/or repair after the breakage of ancestral gene clusters.

Extreme plastid RNA editing may confound phylogenetic reconstruction: A case study of Selaginella (lycophytes)

Cytidine-to-uridine (C-to-U) RNA editing is common in coding regions of organellar genomes throughout land plants. In most cases RNA editing alters translated amino acids or creates new start codons, potentially confounds phylogenetic reconstructions. In this study, we used the spike moss genus Selaginella (lycophytes), which has the highest frequency of RNA editing, as a model to test the effects of extreme RNA editing on phylogenetic reconstruction. We predicted the C-to-U RNA editing sites in coding regions of 18 Selaginella plastomes, and reconstructed the phylogenetic relationships within Selaginella based on three data set pairs consisted of plastome or RNA-edited coding sequences, first and second codon positions, and translated amino acid sequences, respectively. We predicted between 400 and 3100 RNA editing sites of 18 Selaginella plastomes. The numbers of RNA editing sites in plastomes were highly correlated with the GC content of first and second codon positions, but not correlated with the GC content of plastomes as a whole. Contrast phylogenetic analyses showed that there were substantial differences (e.g., the placement of clade B in Selaginella) between the phylogenies generated by the plastome and RNA-edited data sets. This empirical study provides evidence that extreme C-to-U RNA editing in the coding regions of organellar genomes alters the sequences used for phylogenetic reconstruction, and might even confound phylogenetic reconstruction. Therefore, RNA editing sites should be corrected when plastid or mitochondrial genes are used for phylogenetic studies, particularly in those lineages with abundant organellar RNA editing sites, such as hornworts, quillworts, spike mosses, and some seed plants.

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.

The dynamic evolution of mobile open reading frames in plastomes of Hymenophyllum Sm. and new insight on Hymenophyllum coreanum Nakai

In this study, four plastomes of Hymenophyllum, distributed in the Korean peninsula, were newly sequenced and phylogenomic analysis was conducted to reveal (1) the evolutionary history of plastomes of early-diverging fern species at the species level, (2) the importance of mobile open reading frames in the genus, and (3) plastome sequence divergence providing support for H. coreanum to be recognized as an independent species distinct from H. polyanthos. In addition, 1C-values of H. polyanthos and H. coreanum were measured to compare the genome size of both species and to confirm the diversification between them. The rrn16-trnV intergenic regions in the genus varied in length caused by Mobile Open Reading Frames in Fern Organelles (MORFFO). We investigated enlarged noncoding regions containing MORFFO throughout the fern plastomes and found that they were strongly associated with tRNA genes or palindromic elements. Sequence identity between plastomes of H. polyanthos and H. coreanum is quite low at 93.35% in the whole sequence and 98.13% even if the variation in trnV-rrn16 intergenic spacer was ignored. In addition, different genome sizes were found for these species based on the 1C-value. Consequently, there is no reason to consider them as a conspecies.

The Amount of RNA Editing Sites in Liverwort Organellar Genes Is Correlated with GC Content and Nuclear PPR Protein Diversity

RNA editing occurs in the organellar mRNAs of all land plants but the marchantioid liverworts, making liverworts a perfect group for studying the evolution of RNA editing. Here, we profiled the RNA editing of 42 exemplars spanning the ordinal phylogenetic diversity of liverworts, and screened for the nuclear-encoded pentatricopeptide repeat (PPR) proteins in the transcriptome assemblies of these taxa. We identified 7,428 RNA editing sites in 128 organellar genes from 31 non-marchantioid liverwort species, and characterized 25,059 PPR protein sequences. The abundance of organellar RNA editing sites varies greatly among liverwort lineages, genes, and codon positions, and shows strong positive correlations with the GC content of protein-coding genes, and the diversity of the PLS class of nuclear PPR proteins.

Comprehensive genomic analyses with 115 plastomes from algae to seed plants: structure, gene contents, GC contents, and introns

BACKGROUND

Chloroplasts are a common character in plants. The chloroplasts in each plant lineage have shaped their own genomes, plastomes, by structural changes and transferring many genes to nuclear genomes during plant evolution. Some plastid genes have introns that are mostly group II introns.

OBJECTIVE

This study aimed to get genomic and evolutionary insights on the plastomes from green algae to flowering plants.

METHODS

Plastomes of 115 species from green algae, bryophytes, pteridophytes (spore bearing vascular plants), gymnosperms, and angiosperms were mined from NCBI organelle genome database. Plastome structure, gene contents and GC contents were analyzed by the in-house developed Phyton code. Intronic features including presence/absence, length, intron phases were analyzed by manually in the annotated information in NCBI.

RESULTS

The canonical quadripartite structures were retained in most plastomes except of a few plastomes that had lost an invert repeat (IR). Expansion or reduction or deletion of IRs resulted in the length variation of the plastomes. The number of protein coding genes ranged from 40 to 92 with an average 79.43 ± 5.84 per plastome and gene losses were apparent in specific lineages. The number of trn genes ranged from 13 to 33 with an average 21.19 ± 2.42 per plastome. Ribosomal RNA genes, rrn, were located in the IRs so that they were present in a duplicate except of the species that had lost one of the IR. GC contents were variable from 24.9 to 51.0% with an average 38.21 ± 3.27%, indicating bias to high AT contents. Plastid introns were present in 18 protein coding genes, six trn genes, and one rrn gene. Intron losses occurred among the orthologous genes in different plant lineages. The plastid introns were long compared with the nuclear introns, which might be related with the spliceosome nuclear introns and self-splicing group II plastid introns. The trnK-UUU intron contained the maturase encoding matK gene except in the chlorophyte algae and monilophyte ferns in which the trnK-UUU was lost, but matK retained. There were many annotation artefacts in the intron positions in the NCBI database. In the analysis of intron phases, phase 0 introns were more frequent than those of phase 2 and 3 introns. Phase polymorphism was observed in the introns of clpP which was derived from nucleotide insertion. Plastid trn introns were long compared to the archaeal or eukaryotic nuclear tRNA introns. Of the six plastid trn introns, one was at the D loop and other five were at the anticodon loop. The insertion sites were conserved among the trn genes in archaea, eukaryotic nuclear and plastid tRNA genes.

CONCLUSIONS

Current study refurbrished the previous findings of structural variations, gene contents, and GC contents of the chloroplast genomes from green algae to flowering plants. The study also included some noble findings and discussions on the plastome introns including their length variations and phase variation. We also presented and corrected some false annotations on the introns in protein coding and tRNA genes in the genome database, which might be confirmed by the chloroplast transcriptome analysis in the future.

Plastome-based phylogenomics resolves the placement of the sanguinolenta group in the spikemoss of lycophyte (Selaginellaceae)

Selaginellaceae have been shown to be monophyletic in previous studies, and include only the single genus Selaginella. However, the two most recent classifications of the genus disagree in terms of the number of subgenera recognized, and the position of problematic clades such as the "sanguinolenta" group, which has been resolved in quite different positions in different studies. Here, we performed a plastid-genome based phylogenomic analysis of Selaginellaceae to address this problem. The sanguinolenta group, represented here by three species, was resolved as sister to the remaining members of subg. Stachygynandrum. Additionally, subg. Exaltatae, subg. Ericetorum, and subg. Gymnogynum in clade A clustered into a well supported monophyletic clade but with conflicting topology between subgenera inside, which is possibly attributed to the early divergence among them. We uncovered substantial variation in both synonymous (dS) and nonsynonymous (dN) substitution rate, and GC content in plastomes of Selaginellaceae. The values of dS, dN, and GC content were significantly higher than those of other lycophytes (Isoetaceae and Lycopodiaceae). We observed a significant positive correlation between the high GC content, and the elevated dS and dN rates. In addition, the dS and dN values inferred among branches of Selaginellaceae were extremely variable. Our data indicate that this unevenly distributed substitution rate likely reflected relaxed or intensified selection among different lineages, which is possibly related to the inconsistency of the subgeneric phylogenetic topologies of Selaginellaceae.

Robustaflavone-4'-dimethyl ether from Selaginella uncinata attenuated lipopolysaccharide-induced acute lung injury via inhibiting FLT3-mediated neutrophil activation

Neutrophils act as both messenger and effector which contributed to the pathogenesis of acute lung injury (ALI). Targeting neutrophils could be a novel strategy for prevention and treatment of ALI. Selaginella uncinata is widely used as an antitussive, antipyretic and anti-inflammatory herb to treat various pulmonary diseases, including lung cancer, asthma, pulmonary fibrosis and pneumonia. However, its effective constituents remain unknown. In the present study, the protective effects of flavonoids from S. uncinata (SUF) and its major compound robustaflavone-4'-dimethyl ether (RDE) against lipopolysaccharide (LPS)-induced ALI were investigated in mice and in neutrophils.The results showed that both SUF and RDE had the same inhibition on LPS-induced lung edema and neutrophil infiltration as well as the increased levels of IL-6, TNF-α, P-selectin and ICAM-1 in serum of LPS-challenged mice. Furthermore, RDE significantly inhibited inducible neutrophil activation in a concentration-dependent manner, and also reduced the levels of intracellular calcium as well as the expressions of CCR2. Rescue experiment showed that RDE suppressed FLT3 and its downstream p-p38 and p-AKT, which could be abolished by FLT3 agonist FLT3L but partly by MAPK agonist PDBu or AKT agonist SC79. Therefore, these results indicated that RDE as the main bioactive compound in SUF alleviated LPS-induced acute lung injury and inhibited neutrophil activation via inhibition of FLT3-mediatied AKT and MAPK pathways.

Dinoflagellates with relic endosymbiont nuclei as models for elucidating organellogenesis

Nucleomorphs are relic endosymbiont nuclei so far found only in two algal groups, cryptophytes and chlorarachniophytes, which have been studied to model the evolutionary process of integrating an endosymbiont alga into a host-governed plastid (organellogenesis). However, past studies suggest that DNA transfer from the endosymbiont to host nuclei had already ceased in both cryptophytes and chlorarachniophytes, implying that the organellogenesis at the genetic level has been completed in the two systems. Moreover, we have yet to pinpoint the closest free-living relative of the endosymbiotic alga engulfed by the ancestral chlorarachniophyte or cryptophyte, making it difficult to infer how organellogenesis altered the endosymbiont genome. To counter the above issues, we need novel nucleomorph-bearing algae, in which endosymbiont-to-host DNA transfer is on-going and for which endosymbiont/plastid origins can be inferred at a fine taxonomic scale. Here, we report two previously undescribed dinoflagellates, strains MGD and TGD, with green algal endosymbionts enclosing plastids as well as relic nuclei (nucleomorphs). We provide evidence for the presence of DNA in the two nucleomorphs and the transfer of endosymbiont genes to the host (dinoflagellate) genomes. Furthermore, DNA transfer between the host and endosymbiont nuclei was found to be in progress in both the MGD and TGD systems. Phylogenetic analyses successfully resolved the origins of the endosymbionts at the genus level. With the combined evidence, we conclude that the host-endosymbiont integration in MGD/TGD is less advanced than that in cryptophytes/chrorarachniophytes, and propose the two dinoflagellates as models for elucidating organellogenesis.

Evolutionary Model of Plastidial RNA Editing in Angiosperms Presumed from Genome-Wide Analysis of Amborella trichopoda

Amborella trichopoda is placed close to the base of the angiosperm lineage (basal angiosperm). By genome-wide RNA sequencing, we identified 184C-to-U RNA editing sites in the plastid genome of Amborella. This number is much higher than that observed in other angiosperms including maize (44 sites), rice (39 sites) and grape (115 sites). Despite the high frequency of RNA editing, the biased distribution of RNA editing sites in the genome, target codon preference and nucleotide preference adjacent to the edited cytidine are similar to that in other angiosperms, suggesting a common editing machinery. Consistent with this idea, the Amborella nuclear genome encodes 2-3 times more of the E- and DYW-subclass members of pentatricopeptide repeat proteins responsible for RNA editing site recognition in plant organelles. Among 165 editing sites in plastid protein coding sequences in Amborella, 100 sites were conserved at least in one out of 38 species selected to represent key branching points of the angiosperm phylogenetic tree. We assume these 100 sites represent at least a subset of the sites in the plastid editotype of ancestral angiosperms. We then mapped the loss and gain of editing sites on the phylogenetic tree of angiosperms. Our results support the idea that the evolution of angiosperms has led to the loss of RNA editing sites in plastids.

Directed Repeats Co-occur with Few Short-Dispersed Repeats in Plastid Genome of a Spikemoss, Selaginella vardei (Selaginellaceae, Lycopodiopsida)

BACKGROUND

It is hypothesized that the highly conserved inverted repeats (IR) structure of land plant plastid genomes (plastomes) is beneficial for stabilizing plastome organization, whereas the mechanism of the occurrence and stability maintenance of the recently reported direct repeats (DR) structure is yet awaiting further exploration. Here we describe the DR structure of the Selaginella vardei (Selaginellaceae) plastome, to elucidate the mechanism of DR occurrence and stability maintenance.

RESULTS

The plastome of S. vardei is 121,254 bp in length and encodes 76 genes, of which 62 encode proteins, 10 encode tRNAs, and four encode rRNAs. Unexpectedly, the two identical rRNA gene regions (13,893 bp) are arranged in a direct orientation (DR), rather than inverted. Comparing to the IR organization in Isoetes flaccida (Isoetaceae, Lycopodiopsida) plastome, a ca. 50-kb trnN-trnF inversion that spans one DR copy was found in the plastome of S. vardei, which might cause the orientation change. In addition, we find extremely rare short dispersed repeats (SDRs) in the plastomes of S. vardei and its closely related species S. indica.

CONCLUSIONS

We suggest that the ca. 50-kb inversion resulted in the DR structure, and the reduction in SDRs plays a key role in maintaining the stability of plastomes with DR structure by avoiding potential secondary recombination. We further confirmed the presence of homologous recombination between DR regions, which are able to generate subgenomes and form diverse multimers. Our study deepens the understanding of Selaginella plastomes and provides new insights into the diverse plastome structures in land plants.

Phylogeny of Chinese Allium Species in Section Daghestanica and Adaptive Evolution of Allium (Amaryllidaceae, Allioideae) Species Revealed by the Chloroplast Complete Genome

The genus Allium (Amaryllidaceae, Allioideae) is one of the largest monocotyledonous genera and it includes many economically important crops that are cultivated for consumption or medicinal uses. Recent advances in molecular phylogenetics have revolutionized our understanding of Allium taxonomy and evolution. However, the phylogenetic relationships in some Allium sections (such as the Allium section Daghestanica) and the genetic bases of adaptative evolution, remain poorly understood. Here, we newly assembled six chloroplast genomes from Chinese endemic species in Allium section Daghestanica and by combining these genomes with another 35 allied species, we performed a series of analyses including genome structure, GC content, species pairwise Ka/Ks ratios, and the SSR component, nucleotide diversity and codon usage. Positively selected genes (PSGs) were detected in the Allium lineage using the branch-site model. Comparison analysis of Bayesian and ML phylogeny on CCG (complete chloroplast genome), SCG (single copy genes) and CDS (coding DNA sequences) produced a well-resolved phylogeny of Allioideae plastid lineages, which illustrated several novel relationships with the section Daghestanica. In addition, six species in section Daghestanica showed highly conserved structures. The GC content and the GC3s content in Allioideae species exhibited lower values than studied non-Allioideae species, along with elevated pairwise Ka/Ks ratios. The rps2 gene was lost in all examined Allioideae species, and 10 genes with significant posterior probabilities for codon sites were identified in the positive selection analysis, seven of them are associated with photosynthesis. Our study uncovered a new species relationship in section Daghestanica and suggested that the selective pressure has played an important role in Allium adaptation and evolution, these results will facilitate our further understanding of evolution and adaptation of species in the genus Allium.

Extreme variation in rates of evolution in the plastid Clp protease complex

Eukaryotic cells represent an intricate collaboration between multiple genomes, even down to the level of multi-subunit complexes in mitochondria and plastids. One such complex in plants is the caseinolytic protease (Clp), which plays an essential role in plastid protein turnover. The proteolytic core of Clp comprises subunits from one plastid-encoded gene (clpP1) and multiple nuclear genes. TheclpP1 gene is highly conserved across most green plants, but it is by far the fastest evolving plastid-encoded gene in some angiosperms. To better understand these extreme and mysterious patterns of divergence, we investigated the history ofclpP1 molecular evolution across green plants by extracting sequences from 988 published plastid genomes. We find thatclpP1 has undergone remarkably frequent bouts of accelerated sequence evolution and architectural changes (e.g. a loss of introns andRNA-editing sites) within seed plants. AlthoughclpP1 is often assumed to be a pseudogene in such cases, multiple lines of evidence suggest that this is rarely true. We applied comparative native gel electrophoresis of chloroplast protein complexes followed by protein mass spectrometry in two species within the angiosperm genusSilene, which has highly elevated and heterogeneous rates ofclpP1 evolution. We confirmed thatclpP1 is expressed as a stable protein and forms oligomeric complexes with the nuclear-encoded Clp subunits, even in one of the most divergentSilene species. Additionally, there is a tight correlation between amino acid substitution rates inclpP1 and the nuclear-encoded Clp subunits across a broad sampling of angiosperms, suggesting continuing selection on interactions within this complex.

The Unique Evolutionary Trajectory and Dynamic Conformations of DR and IR/DR-Coexisting Plastomes of the Early Vascular Plant Selaginellaceae (Lycophyte)

Both direct repeats (DR) and inverted repeats (IR) are documented in the published plastomes of Selaginella species indicating the unusual and diverse plastome structure in the family Selaginellaceae. In this study, we newly sequenced complete plastomes of seven species from five main lineages of Selaginellaceae and also resequenced three species (Selaginella tamariscina, Selaginella uncinata, and Selaginella moellendorffii) to explore the evolutionary trajectory of Selaginellaceae plastomes. Our results showed that the plastomes of Selaginellaceae vary remarkably in size, gene contents, gene order, and GC contents. Notably, both DR and IR structures existed in the plastomes of Selaginellaceae with DR structure being an ancestral state. The occurrence of DR structure was at ∼257 Ma and remained in most subgenera of Selaginellaceae, whereas IR structure only reoccurred in Selaginella sect. Lepidophyllae (∼143 Ma) and Selaginella subg. Heterostachys (∼19 Ma). The presence of a pair of large repeats psbK-trnQ, together with DR/IR region in Selaginella bisulcata, Selaginella pennata, S. uncinata, and Selaginella hainanensis, could frequently mediate diverse homologous recombination and create approximately equal stoichiometric isomers (IR/DR-coexisting) and subgenomes. High proportion of repeats is presumably responsible for the dynamic IR/DR-coexisting plastomes, which possess a lower synonymous substitution rate (dS) compared with DR-possessing and IR-possessing plastomes. We propose that the occurrence of DR structure, together with few repeats, is possibly selected to keep the stability of plastomes and the IR/DR-coexisting plastomes also reached an equilibrium in plastome organization through highly efficient homologous recombination to maintain stability.

Lycophyte plastid genomics: extreme variation in GC, gene and intron content and multiple inversions between a direct and inverted orientation of the rRNA repeat

Lycophytes are a key group for understanding vascular plant evolution. Lycophyte plastomes are highly distinct, indicating a dynamic evolutionary history, but detailed evaluation is hindered by the limited availability of sequences. Eight diverse plastomes were sequenced to assess variation in structure and functional content across lycophytes. Lycopodiaceae plastomes have remained largely unchanged compared with the common ancestor of land plants, whereas plastome evolution in Isoetes and especially Selaginella is highly dynamic. Selaginella plastomes have the highest GC content and fewest genes and introns of any photosynthetic land plant. Uniquely, the canonical inverted repeat was converted into a direct repeat (DR) via large-scale inversion in some Selaginella species. Ancestral reconstruction identified additional putative transitions between an inverted and DR orientation in Selaginella and Isoetes plastomes. A DR orientation does not disrupt the activity of copy-dependent repair to suppress substitution rates within repeats. Lycophyte plastomes include the most archaic examples among vascular plants and the most reconfigured among land plants. These evolutionary trends correlate with the mitochondrial genome, suggesting shared underlying mechanisms. Copy-dependent repair for DR-localized genes indicates that recombination and gene conversion are not inhibited by the DR orientation. Gene relocation in lycophyte plastomes occurs via overlapping inversions rather than transposase/recombinase-mediated processes.

Plastid phylogenomic analysis of green plants: A billion years of evolutionary history

PREMISE OF THE STUDY

For the past one billion years, green plants (Viridiplantae) have dominated global ecosystems, yet many key branches in their evolutionary history remain poorly resolved. Using the largest analysis of Viridiplantae based on plastid genome sequences to date, we examined the phylogeny and implications for morphological evolution at key nodes.

METHODS

We analyzed amino acid sequences from protein-coding genes from complete (or nearly complete) plastomes for 1879 taxa, including representatives across all major clades of Viridiplantae. Much of the data used was derived from transcriptomes from the One Thousand Plants Project (1KP); other data were taken from GenBank.

KEY RESULTS

Our results largely agree with previous plastid-based analyses. Noteworthy results include (1) the position of Zygnematophyceae as sister to land plants (Embryophyta), (2) a bryophyte clade (hornworts, mosses + liverworts), (3) Equisetum + Psilotaceae as sister to Marattiales + leptosporangiate ferns, (4) cycads + Ginkgo as sister to the remaining extant gymnosperms, within which Gnetophyta are placed within conifers as sister to non-Pinaceae (Gne-Cup hypothesis), and (5) Amborella, followed by water lilies (Nymphaeales), as successive sisters to all other extant angiosperms. Within angiosperms, there is support for Mesangiospermae, a clade that comprises magnoliids, Chloranthales, monocots, Ceratophyllum, and eudicots. The placements of Ceratophyllum and Dilleniaceae remain problematic. Within Pentapetalae, two major clades (superasterids and superrosids) are recovered.

CONCLUSIONS

This plastid data set provides an important resource for elucidating morphological evolution, dating divergence times in Viridiplantae, comparisons with emerging nuclear phylogenies, and analyses of molecular evolutionary patterns and dynamics of the plastid genome.

Selaginella guihaia (Selaginellaceae): A new spikemoss species from southern China and northern Vietnam around the Gulf of Tonkin

Selaginella guihaiasp. nov. (Selaginellaceae), a new species of spikemoss from southern China and northern Vietnam around the Gulf of Tonkin (Beibu Gulf), is described and illustrated. Morphological and molecular comparisons of the new species with other similar species (S. doederleinii, S. ornata and S. trachyphylla) are provided. The morphological and molecular evidence clearly indicates S. guihaia is a distinct species. Morphologically S. guihaia differs from other species by its obviously white-margined leaves, the ventral leaves scabrous on upper surfaces throughout the basiscopic or also rarely present on upper halves, and the ovate axillary leaves.

Complete mitochondrial genomes from the ferns Ophioglossum californicum and Psilotum nudum are highly repetitive with the largest organellar introns

Currently, complete mitochondrial genomes (mitogenomes) are available from all major land plant lineages except ferns. Sequencing of fern mitogenomes could shed light on the major evolutionary transitions that established mitogenomic diversity among extant lineages. In this study, we generated complete mitogenomes from the adder's tongue fern (Ophioglossum californicum) and the whisk fern (Psilotum nudum). The Psilotum mitogenome (628 kb) contains a rich complement of genes and introns, some of which are the largest of any green plant organellar genome. In the Ophioglossum mitogenome (372 kb), gene and intron content is slightly reduced, including the loss of all four mitochondrial ccm genes. Transcripts of nuclear Ccm genes also were not detected, suggesting loss of the entire mitochondrial cytochrome c maturation pathway from Ophioglossum. Both fern mitogenomes are highly repetitive, yet they show extremely low levels of active recombination. Transcriptomic sequencing uncovered ˜1000 sites of C-to-U RNA editing in both species, plus a small number (< 60) of U-to-C edit sites. Overall, the first mitochondrial genomes of ferns show a mix of features shared with lycophytes and/or seed plants and several novel genomic features, enabling a robust reconstruction of the mitogenome in the common ancestor of vascular plants.

Unexpectedly Streamlined Mitochondrial Genome of the Euglenozoan Euglena gracilis

In this study, we describe the mitochondrial genome of the excavate flagellate Euglena gracilis. Its gene complement is reduced as compared with the well-studied sister groups Diplonemea and Kinetoplastea. We have identified seven protein-coding genes: Three subunits of respiratory complex I (nad1, nad4, and nad5), one subunit of complex III (cob), and three subunits of complex IV (cox1, cox2, and a highly divergent cox3). Moreover, fragments of ribosomal RNA genes have also been identified. Genes encoding subunits of complex V, ribosomal proteins and tRNAs were missing, and are likely located in the nuclear genome. Although mitochondrial genomes of diplonemids and kinetoplastids possess the most complex RNA processing machineries known, including trans-splicing and editing of the uridine insertion/deletion type, respectively, our transcriptomic data suggest their total absence in E. gracilis. This finding supports a scenario in which the complex mitochondrial processing machineries of both sister groups evolved relatively late in evolution from a streamlined genome and transcriptome of their common predecessor.

Comparative Bioinformatics Analysis of the Chloroplast Genomes of a Wild Diploid Gossypium and Two Cultivated Allotetraploid Species

BACKGROUND

Gossypium thurberi is a wild diploid species that has been used to improve cultivated allotetraploid cotton.G. thurberi belongs to D genome, which is an important wild bio-source for the cotton breeding and genetic research. To a certain degree, chloroplast DNA sequence information are a versatile tool for species identification and phylogenetic implications in plants. Different chloroplast loci have been utilized for evaluating phylogenetic relationships at each classification level among plant species, including at the interspecies and intraspecies levels. Present study was conducted in order to analyse the sequence of its chloroplast.

OBJECTIVES

Present study was conducted to study and compare the complete chloroplast sequence of G. thurberi, analyses of its genome structure, gene content and organization, repeat sequence and codon usage and comparison with two cultivated allotetraploid sequenced cotton species.

MATERIALS AND METHODS

The available sequence was assembled by DNAman (Version 8.1.2.378). Gene annotation was mainly performed by DOGMA. The map of genome structure and gene distribution were carried out using OGDRAW V1.1. Relative synonymous codon usage (RSCU) of different codons in each gene sample was calculated by codonW in Mobyle. To determine the repeat sequence and location, an online version of REPuter was used.

RESULTS

The G. thurberi chloroplast (cp) genome is 160264 bp in length with conserved quadripartite structure. Single copy region of cp genome is separated by the two inverted regions. The large single copy region is 88,737 bp, and the small single copy region is 20,271 bp whereas the inverted repeat is 25,628 bp each. The plastidic genome has 113 single genes and 20 duplicated genes. The singletones encode 79 proteins, 4 ribosomal RNA genes and 30 transfer RNA genes.

CONCLUSIONS

Amongst all plastidic genes only 18 genes appeared to have 1-2 introns and when compared with cpDNA of two cultivated allotetraploid, rps18 was the only duplicated gene in G.thurberi. Despite the high level of conservation in cp genome SSRs ,these are useful in analysis of genetic diversity due to their greater efficiency as opposed to genomic SSRs. Low GC content is a significant feature of plastidic genomes, which is possibly formed after endosymbiosis by DNA replication and repair.

Evolutionary Stasis in Cycad Plastomes and the First Case of Plastome GC-Biased Gene Conversion

In angiosperms, gene conversion has been known to reduce the mutational load of plastid genomes (the plastomes). Particularly, more frequent gene conversions in inverted repeat (IR) than in single copy (SC) regions result in contrasting substitution rates between these two regions. However, little has been known about the effect of gene conversion in the evolution of gymnosperm plastomes. Cycads (Cycadophyta) are the second largest gymnosperm group. Evolutionary study of their plastomes is limited to the basal cycad genus, Cycas. In this study, we addressed three questions. 1) Do the plastomes of other cycad genera evolve slowly as previously observed in the plastome of Cycas taitungensis? 2) Do substitution rates differ between their SC and IR regions? And 3) Does gene conversion occur in the cycad plastomes? If yes, is it AT-biased or GC-biased? Plastomes of eight species from other eight genera of cycads were sequenced. These plastomes are highly conserved in genome organization. Excluding ginkgo, cycad plastomes have significantly lower synonymous and nonsynonymous substitution rates than other gymnosperms, reflecting their evolutionary stasis in nucleotide mutations. In the IRs of cycad plastomes, the reduced substitution rates and GC-biased mutations are associated with a GC-biased gene conversion (gBGC) mechanism. Further investigations suggest that in cycads, gBGC is able to rectify plastome-wide mutations. Therefore, this study is the first to uncover the plastomic gBGC in seed plants. We also propose a gBGC model to interpret the dissimilar evolutionary patterns as well as the compositionally biased mutations in the SC and IR regions of cycad plastomes.

The past, present and future of mitochondrial genomics: have we sequenced enough mtDNAs?

The year 2014 saw more than a thousand new mitochondrial genome sequences deposited in GenBank—an almost 15% increase from the previous year. Hundreds of peer-reviewed articles accompanied these genomes, making mitochondrial DNAs (mtDNAs) the most sequenced and reported type of eukaryotic chromosome. These mtDNA data have advanced a wide range of scientific fields, from forensics to anthropology to medicine to molecular evolution. But for many biological lineages, mtDNAs are so well sampled that newly published genomes are arguably no longer contributing significantly to the progression of science, and in some cases they are tying up valuable resources, particularly journal editors and referees. Is it time to acknowledge that as a research community we have published enough mitochondrial genome papers? Here, I address this question, exploring the history, milestones and impacts of mitochondrial genomics, the benefits and drawbacks of continuing to publish mtDNAs at a high rate and what the future may hold for such an important and popular genetic marker. I highlight groups for which mtDNAs are still poorly sampled, thus meriting further investigation, and recommend that more energy be spent characterizing aspects of mitochondrial genomes apart from the DNA sequence, such as their chromosomal and transcriptional architectures. Ultimately, one should be mindful before writing a mitochondrial genome paper. Consider perhaps sending the sequence directly to GenBank instead, and be sure to annotate it correctly before submission.

Mitochondrial and plastid genome architecture: Reoccurring themes, but significant differences at the extremes

Mitochondrial and plastid genomes show a wide array of architectures, varying immensely in size, structure, and content. Some organelle DNAs have even developed elaborate eccentricities, such as scrambled coding regions, nonstandard genetic codes, and convoluted modes of posttranscriptional modification and editing. Here, we compare and contrast the breadth of genomic complexity between mitochondrial and plastid chromosomes. Both organelle genomes have independently evolved many of the same features and taken on similar genomic embellishments, often within the same species or lineage. This trend is most likely because the nuclear-encoded proteins mediating these processes eventually leak from one organelle into the other, leading to a high likelihood of processes appearing in both compartments in parallel. However, the complexity and intensity of genomic embellishments are consistently more pronounced for mitochondria than for plastids, even when they are found in both compartments. We explore the evolutionary forces responsible for these patterns and argue that organelle DNA repair processes, mutation rates, and population genetic landscapes are all important factors leading to the observed convergence and divergence in organelle genome architecture.