Origin and Evolution of the Chloroplast trnK (matK) Intron: A Model for Evolution of Group II Intron RNA Structures

GhCTSF1, a short PPR protein with a conserved role in chloroplast development and photosynthesis, participates in intron splicing of rpoC1 and ycf3-2 transcripts in cotton

Cotton is one of the most important textile fibers worldwide. As crucial agronomic traits, leaves play an essential role in the growth, disease resistance, fiber quality, and yield of cotton plants. Pentatricopeptide repeat (PPR) proteins are a large family of nuclear-encoded proteins involved in organellar or nuclear RNA metabolism. Using a virus-induced gene silencing assay, we found that cotton plants displayed variegated yellow leaf phenotypes with decreased chlorophyll content when expression of the PPR gene GhCTSF1 was silenced. GhCTSF1 encodes a chloroplast-localized protein that contains only two PPR motifs. Disruption of GhCTSF1 substantially reduces the splicing efficiency of rpoC1 intron 1 and ycf3 intron 2. Loss of function of the GhCTSF1 ortholog EMB1417 causes splicing defects in rpoC1 and ycf3-2, leading to impaired chloroplast structure and decreased photosynthetic rates in Arabidopsis. We also found that GhCTSF1 interacts with two splicing factors, GhCRS2 and GhWTF1. Defects in GhCRS2 and GhWTF1 severely affect intron splicing of rpoC1 and ycf3-2 in cotton, leading to defects in chloroplast development and a reduction in photosynthesis. Our results suggest that GhCTSF1 is specifically required for splicing rpoC1 and ycf3-2 in cooperation with GhCRS2 and GhWTF1.

Plant organellar RNA maturation

Plant organellar RNA metabolism is run by a multitude of nucleus-encoded RNA-binding proteins (RBPs) that control RNA stability, processing, and degradation. In chloroplasts and mitochondria, these post-transcriptional processes are vital for the production of a small number of essential components of the photosynthetic and respiratory machinery-and consequently for organellar biogenesis and plant survival. Many organellar RBPs have been functionally assigned to individual steps in RNA maturation, often specific to selected transcripts. While the catalog of factors identified is ever-growing, our knowledge of how they achieve their functions mechanistically is far from complete. This review summarizes the current knowledge of plant organellar RNA metabolism taking an RBP-centric approach and focusing on mechanistic aspects of RBP functions and the kinetics of the processes they are involved in.

Plastid Genome of Equisetum xylochaetum from the Atacama Desert, Chile and the Relationships of Equisetum Based on Frequently Used Plastid Genes and Network Analysis

The modern pteridophyte genus Equisetum is the only survivor of Sphenopsida, an ancient clade known from the Devonian. This genus, of nearly worldwide distribution, comprises approximately 15 extant species. However, genomic information is limited. In this study, we assembled the complete chloroplast genome of the giant species Equisetum xylochaetum from a metagenomic sequence and compared the plastid genome structure and protein-coding regions with information available for two other Equisetum species using network analysis. Equisetum chloroplast genomes showed conserved traits of quadripartite structure, gene content, and gene order. Phylogenetic analysis based on plastome protein-coding regions corroborated previous reports that Equisetum is monophyletic, and that E. xylochaetum is more closely related to E. hyemale than to E. arvense. Single-gene phylogenetic estimation and haplotype analysis showed that E. xylochaetum belonged to the subgenus Hippochaete. Single-gene haplotype analysis revealed that E. arvense, E. hyemale, E. myriochaetum, and E. variegatum resolved more than one haplotype per species, suggesting the presence of a high diversity or a high mutation rate of the corresponding nucleotide sequence. Sequences from E. bogotense appeared as a distinct group of haplotypes representing the subgenus Paramochaete that diverged from Hippochaete and Equisetum. In addition, the taxa that were frequently located at the joint region of the map were E. scirpoides and E. pratense, suggesting the presence of some plastome characters among the Equiseum subgenera.

Organellar Introns in Fungi, Algae, and Plants

Introns are ubiquitous in eukaryotic genomes and have long been considered as ‘junk RNA’ but the huge energy expenditure in their transcription, removal, and degradation indicate that they may have functional significance and can offer evolutionary advantages. In fungi, plants and algae introns make a significant contribution to the size of the organellar genomes. Organellar introns are classified as catalytic self-splicing introns that can be categorized as either Group I or Group II introns. There are some biases, with Group I introns being more frequently encountered in fungal mitochondrial genomes, whereas among plants Group II introns dominate within the mitochondrial and chloroplast genomes. Organellar introns can encode a variety of proteins, such as maturases, homing endonucleases, reverse transcriptases, and, in some cases, ribosomal proteins, along with other novel open reading frames. Although organellar introns are viewed to be ribozymes, they do interact with various intron- or nuclear genome-encoded protein factors that assist in the intron RNA to fold into competent splicing structures, or facilitate the turn-over of intron RNAs to prevent reverse splicing. Organellar introns are also known to be involved in non-canonical splicing, such as backsplicing and trans-splicing which can result in novel splicing products or, in some instances, compensate for the fragmentation of genes by recombination events. In organellar genomes, Group I and II introns may exist in nested intronic arrangements, such as introns within introns, referred to as twintrons, where splicing of the external intron may be dependent on splicing of the internal intron. These nested or complex introns, with two or three-component intron modules, are being explored as platforms for alternative splicing and their possible function as molecular switches for modulating gene expression which could be potentially applied towards heterologous gene expression. This review explores recent findings on organellar Group I and II introns, focusing on splicing and mobility mechanisms aided by associated intron/nuclear encoded proteins and their potential roles in organellar gene expression and cross talk between nuclear and organellar genomes. Potential application for these types of elements in biotechnology are also discussed.

Comparative Chloroplast Genomes of Four Lycoris Species (Amaryllidaceae) Provides New Insight into Interspecific Relationship and Phylogeny

The genus Lycoris (Amaryllidaceae) consists of about 20 species, which is endemic to East Asia. Although the Lycoris species is of great horticultural and medical importance, challenges in accurate species identification persist due to frequent natural hybridization and large-scale intraspecific variation. In this study, we sequenced chloroplast genomes of four Lycoris species and retrieved seven published chloroplast (cp) genome sequences in this genus for comparative genomic and phylogenetic analyses. The cp genomes of these four newly sequenced species were found to be 158,405-158,498 bp with the same GC content of 37.8%. The structure of the genomes exhibited the typical quadripartite structure with conserved gene order and content. A total of 113 genes (20 duplicated) were identified, including 79 protein-coding genes (PCGs), 30 tRNAs, and 4 rRNAs. Phylogenetic analysis showed that the 11 species were clustered into three main groups, and L. sprengeri locate at the base of Lycoriss. The L. radiata was suggested to be the female donor of the L. incarnata, L. shaanxiensis, and L. squamigera. The L. straminea and L. houdyshelii may be derived from L. anhuiensis, L. chinensis, or L. longituba. These results could not only offer a genome-scale platform for identification and utilization of Lycoris but also provide a phylogenomic framework for future studies in this genus.

nMAT3 is an essential maturase splicing factor required for holo-complex I biogenesis and embryo development in Arabidopsis thaliana plants

Group-II introns are self-splicing mobile genetic elements consisting of catalytic intron-RNA and its related intron-encoded splicing maturase protein cofactor. Group-II sequences are particularly plentiful within the mitochondria of land plants, where they reside within many critical gene loci. During evolution, the plant organellar introns have degenerated, such as they lack regions that are are required for splicing, and also lost their evolutionary related maturase proteins. Instead, for their splicing the organellar introns in plants rely on different host-acting protein cofactors, which may also provide a means to link cellular signals with respiratory functions. The nuclear genome of Arabidopsis thaliana encodes four maturase-related factors. Previously, we showed that three of the maturases, nMAT1, nMAT2 and nMAT4, function in the excision of different group-II introns in Arabidopsis mitochondria. The function of nMAT3 (encoded by the At5g04050 gene locus) was found to be essential during early embryogenesis. Using a modified embryo-rescue method, we show that nMAT3-knockout plants are strongly affected in the splicing of nad1 introns 1, 3 and 4 in Arabidopsis mitochondria, resulting in complex-I biogenesis defects and altered respiratory activities. Functional complementation of nMAT3 restored the organellar defects and embryo-arrested phenotypes associated with the nmat3 mutant line. Notably, nMAT3 and nMA4 were found to act on the same RNA targets but have no redundant functions in the splicing of nad1 transcripts. The two maturases, nMAT3 and nMAT4 are likely to cooperate together in the maturation of nad1 pre-RNAs. Our results provide important insights into the roles of maturases in mitochondria gene expression and the biogenesis of the respiratory system during early plant life.

A comparative analysis of the complete chloroplast genomes of three Chrysanthemum boreale strains

BACKGROUND

Chrysanthemum boreale Makino (Anthemideae, Asteraceae) is a plant of economic, ornamental and medicinal importance. We characterized and compared the chloroplast genomes of three C. boreale strains. These were collected from different geographic regions of Korea and varied in floral morphology.

METHODS

The chloroplast genomes were obtained by next-generation sequencing techniques, assembled de novo, annotated, and compared with one another. Phylogenetic analysis placed them within the Anthemideae tribe.

RESULTS

The sizes of the complete chloroplast genomes of the C. boreale strains were 151,012 bp (strain 121002), 151,098 bp (strain IT232531) and 151,010 bp (strain IT301358). Each genome contained 80 unique protein-coding genes, 4 rRNA genes and 29 tRNA genes. Comparative analyses revealed a high degree of conservation in the overall sequence, gene content, gene order and GC content among the strains. We identified 298 single nucleotide polymorphisms (SNPs) and 106 insertions/deletions (indels) in the chloroplast genomes. These variations were more abundant in non-coding regions than in coding regions. Long dispersed repeats and simple sequence repeats were present in both coding and noncoding regions, with greater frequency in the latter. Regardless of their location, these repeats can be used for molecular marker development. Phylogenetic analysis revealed the evolutionary relationship of the species in the Anthemideae tribe. The three complete chloroplast genomes will be valuable genetic resources for studying the population genetics and evolutionary relationships of Asteraceae species.

Zygnema circumcarinatum UTEX 1559 chloroplast and mitochondrial genomes provide insight into land plant evolution

The complete chloroplast and mitochondrial genomes of Charophyta have shed new light on land plant terrestrialization. Here, we report the organellar genomes of the Zygnema circumcarinatum strain UTEX 1559, and a comparative genomics investigation of 33 plastomes and 18 mitogenomes of Chlorophyta, Charophyta (including UTEX 1559 and its conspecific relative SAG 698-1a), and Embryophyta. Gene presence/absence was determined across these plastomes and mitogenomes. A comparison between the plastomes of UTEX 1559 (157 548 bp) and SAG 698-1a (165 372 bp) revealed very similar gene contents, but substantial genome rearrangements. Surprisingly, the two plastomes share only 85.69% nucleotide sequence identity. The UTEX 1559 mitogenome size is 215 954 bp, the largest among all sequenced Charophyta. Interestingly, this large mitogenome contains a 50 kb region without homology to any other organellar genomes, which is flanked by two 86 bp direct repeats and contains 15 ORFs. These ORFs have significant homology to proteins from bacteria and plants with functions such as primase, RNA polymerase, and DNA polymerase. We conclude that (i) the previously published SAG 698-1a plastome is probably from a different Zygnema species, and (ii) the 50 kb region in the UTEX 1559 mitogenome might be recently acquired as a mobile element.

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.

Retroelement origins of pre-mRNA splicing

Recent cryo-EM structures of a group II intron caught in the process of invading DNA have given new insight into the mechanisms of both splicing and retrotransposition. Conformational dynamics involving the branch-site helix domain VI are responsible for substrate exchange between the two steps of splicing. These structural rearrangements have strong parallels with the movement of the branch-site helix in the spliceosome during catalysis. This is strong evidence for the spliceosome evolving from a group II intron ancestor. We observe other topological changes in the overall structure of the catalytic domain V that may occur in the spliceosome as well. Therefore, studying group II introns not only provides us with insight into the evolutionary origins of the spliceosome, but also may inform the design of experiments to further probe structure-function relationships in this eukaryotic splicing apparatus. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution.

Sequencing the Plastid Genome of Giant Ragweed (Ambrosia trifida, Asteraceae) From a Herbarium Specimen

We report the first plastome sequence of giant ragweed (Ambrosia trifida); with this new genome information, we assessed the phylogeny of Asteraceae and the transcriptional profiling against glyphosate resistance in giant ragweed. Assembly and genic features show a normal angiosperm quadripartite plastome structure with no signatures of deviation in gene directionality. Comparative analysis revealed large inversions across the plastome of giant ragweed and the previously sequenced members of the plant family. Asteraceae plastid genomes contain two inversions of 22.8 and 3.3 kb; the former is located between trnS-GCU and trnG-UCC genes, and the latter between trnE-UUC and trnT-GGU genes. The plastid genome sequences of A. trifida and the related species, Ambrosia artemisiifolia, are identical in gene content and arrangement, but they differ in length. The phylogeny is well-resolved and congruent with previous hypotheses about the phylogenetic relationship of Asteraceae. Transcriptomic analysis revealed divergence in the relative expressions at the exonic and intronic levels, providing hints toward the ecological adaptation of the genus. Giant ragweed shows various levels of glyphosate resistance, with introns displaying higher expression patterns at resistant time points after the assumed herbicide treatment.

Plastome-Wide Rearrangements and Gene Losses in Carnivorous Droseraceae

The plastid genomes of four related carnivorous plants (Drosera regia, Drosera erythrorhiza, Aldrovanda vesiculosa, and Dionaea muscipula) were sequenced to examine changes potentially induced by the transition to carnivory. The plastid genomes of the Droseraceae show multiple rearrangements, gene losses, and large expansions or contractions of the inverted repeat. All the ndh genes are lost or nonfunctional, as well as in some of the species, clpP1, ycf1, ycf2 and some tRNA genes. Uniquely, among land plants, the trnK gene has no intron. Carnivory in the Droseraceae coincides with changes in plastid gene content similar to those induced by parasitism and mycoheterotrophy, suggesting parallel changes in chloroplast function due to the similar switch from autotrophy to (mixo-) heterotrophy. A molecular phylogeny of the taxa based on all shared plastid genes indicates that the "snap-traps" of Aldrovanda and Dionaea have a common origin.

Full plastome sequence of the fern Vandenboschia speciosa (Hymenophyllales): structural singularities and evolutionary insights

We provide here the first full chloroplast genome sequence, i.e., the plastome, for a species belonging to the fern order Hymenophyllales. The phylogenetic position of this order within leptosporangiate ferns, together with the general scarcity of information about fern plastomes, places this research as a valuable study on the analysis of the diversity of plastomes throughout fern evolution. Gene content of V. speciosa plastome was similar to that in most ferns, although there were some characteristic gene losses and lineage-specific differences. In addition, an important number of genes required U to C RNA editing for proper protein translation and two genes showed start codons alternative to the canonical AUG (AUA). Concerning gene order, V. speciosa shared the specific 30-kb inversion of euphyllophytes plastomes and the 3.3-kb inversion of fern plastomes, keeping the ancestral gene order shared by eusporangiate and early leptosporangiate ferns. Conversely, V. speciosa has expanded IR regions comprising the rps7, rps12, ndhB and trnL genes in addition to rRNA and other tRNA genes, a condition shared with several eusporangiate ferns, lycophytes and hornworts, as well as most seed plants.

Characterization of chloroplast genomes of Alnus rubra and Betula cordifolia, and their use in phylogenetic analyses in Betulaceae

OBJECTIVE

Betulaceae is a relatively small birch family that comprises about 160 deciduous trees and shrubs. Chloroplast (cp) genome sequencing of Alnus rubra and Betula cordifolia was carried out to elucidate their molecular features and phylogenetic relationship among species in Betulaceae family.

METHODS

Chloroplast genome sequencing was carried out using next generation sequencing method. Molecular and genomic features of the two cp genomes were characterized with other cp genomes in Betulaceae. Also, molecular phylogenetic analysis was performed using the whole cp genome sequences.

RESULTS

The average cp genome length was 160,136 bp among the Betulaceae species. Base compositions of the cp genomes were skewed toward a high AT ratio, with an average of 63.4%. We identified 117 different genes 83 with protein coding, 4 with ribosomal RNA, and 30 with tRNA. Eighteen genes contained introns which were conserved among the cp genomes of all Betulaceae. We mined 82 SSRs from the cp genomes of A. rubra, A. cordifolia, and A. nana. The SSRs were variable in motif repeat numbers and presence/absence among the cp genomes.

CONCLUSION

Chloroplast genome-wide sequence comparison from 11 Betulaceae species and one cp genome of evergreen oak revealed that the patterns of sequence variations were congruent with two subfamily classification Betuloideae (Alnus and Betula) and Corylaceae (Corylus, Ostrya, and Carpinus). Subsequent phylogenetic analysis also supports the sub-classifications of these species.

Dynamic evolution of inverted repeats in Euglenophyta plastid genomes

Photosynthetic euglenids (Euglenophyta) are a monophyletic group of unicellular eukaryotes characterized by the presence of plastids, which arose as the result of the secondary endosymbiosis. Many Euglenophyta plastid (pt) genomes have been characterized recently, but they represented mainly one family - Euglenaceae. Here, we report a comparative analysis of plastid genomes from eight representatives of the family Phacaceae. Newly sequenced plastid genomes share a number of features including synteny and gene content, except for genes mat2 and mat5 encoding maturases. The observed diversity of intron number and presence/absence of maturases corroborated previously suggested correlation between the number of maturases in the pt genome and intron proliferation. Surprisingly, pt genomes of taxa belonging to Discoplastis and Lepocinclis encode two inverted repeat (IR) regions containing the rDNA operon, which are absent from the Euglenaceae. By mapping the presence/absence of IR region on the obtained phylogenomic tree, we reconstructed the most probable events in the evolution of IRs in the Euglenophyta. Our study highlights the dynamic nature of the Euglenophyta plastid genome, in particular with regards to the IR regions that underwent losses repeatedly.

The Complete Plastome Sequence of an Antarctic Bryophyte Sanionia uncinata (Hedw.) Loeske

Organellar genomes of bryophytes are poorly represented with chloroplast genomes of only four mosses, four liverworts and two hornworts having been sequenced and annotated. Moreover, while Antarctic vegetation is dominated by the bryophytes, there are few reports on the plastid genomes for the Antarctic bryophytes. Sanionia uncinata (Hedw.) Loeske is one of the most dominant moss species in the maritime Antarctic. It has been researched as an important marker for ecological studies and as an extremophile plant for studies on stress tolerance. Here, we report the complete plastome sequence of S. uncinata, which can be exploited in comparative studies to identify the lineage-specific divergence across different species. The complete plastome of S. uncinata is 124,374 bp in length with a typical quadripartite structure of 114 unique genes including 82 unique protein-coding genes, 37 tRNA genes and four rRNA genes. However, two genes encoding the α subunit of RNA polymerase (rpoA) and encoding the cytochrome b_6/f complex subunit VIII (petN) were absent. We could identify nuclear genes homologous to those genes, which suggests that rpoA and petN might have been relocated from the chloroplast genome to the nuclear genome.

Molecular Mechanism and Evolution of Nuclear Pre-mRNA and Group II Intron Splicing: Insights from Cryo-Electron Microscopy Structures

Nuclear pre-mRNA splicing and group II intron self-splicing both proceed by two-step transesterification reactions via a lariat intron intermediate. Recently determined cryo-electron microscopy (cryo-EM) structures of catalytically active spliceosomes revealed the RNA-based catalytic core and showed how pre-mRNA substrates and reaction products are positioned in the active site. These findings highlight a strong structural similarity to the group II intron active site, strengthening the notion that group II introns and spliceosomes evolved from a common ancestor. Prp8, the largest and most conserved protein in the spliceosome, cradles the active site RNA. Prp8 and group II intron maturase have a similar domain architecture, suggesting that they also share a common evolutionary origin. The interactions between maturase and key group II intron RNA elements, such as the exon-binding loop and domains V and VI, are recapitulated in the interactions between Prp8 and key elements in the spliceosome's catalytic RNA core. Structural comparisons suggest that the extensive RNA scaffold of the group II intron was gradually replaced by proteins as the spliceosome evolved. A plausible model of spliceosome evolution is discussed.

Choosing and Using Introns in Molecular Phylogenetics

Introns are now commonly used in molecular phylogenetics in an attempt to recover gene trees that are concordant with species trees, but there are a range of genomic, logistical and analytical considerations that are infrequently discussed in empirical studies that utilize intron data. This review outlines expedient approaches for locus selection, overcoming paralogy problems, recombination detection methods and the identification and incorporation of LVHs in molecular systematics. A range of parsimony and Bayesian analytical approaches are also described in order to highlight the methods that can currently be employed to align sequences and treat indels in subsequent analyses. By covering the main points associated with the generation and analysis of intron data, this review aims to provide a comprehensive introduction to using introns (or any non-coding nuclear data partition) in contemporary phylogenetics.

Mobile Group II Introns as Ancestral Eukaryotic Elements

The duality of group II introns, capable of carrying out both self-splicing and retromobility reactions, is hypothesized to have played a profound role in the evolution of eukaryotes. These introns likely provided the framework for the emergence of eukaryotic retroelements, spliceosomal introns and other key components of the spliceosome. Group II introns are found in all three domains of life and are therefore considered to be exceptionally successful mobile genetic elements. Initially identified in organellar genomes, group II introns are found in bacteria, chloroplasts, and mitochondria of plants and fungi, but not in nuclear genomes. Although there is no doubt that prokaryotic and organellar group II introns are evolutionary related, there are remarkable differences in survival strategies between them. Furthermore, an evolutionary relationship of group II introns to eukaryotic retroelements, including telomeres, and spliceosomes is unmistakable.

Group II introns: versatile ribozymes and retroelements

Group II introns are catalytic RNAs (ribozymes) and retroelements found in the genomes of bacteria, archaebacteria, and organelles of some eukaryotes. The prototypical retroelement form consists of a structurally conserved RNA and a multidomain reverse transcriptase protein, which interact with each other to mediate splicing and mobility reactions. A wealth of biochemical, cross-linking, and X-ray crystal structure studies have helped to reveal how the two components cooperate to carry out the splicing and mobility reactions. In addition to the standard retroelement form, group II introns have evolved into derivative forms by either losing specific splicing or mobility characteristics, or becoming functionally specialized. Of particular interest are the eukaryotic derivatives-the spliceosome, spliceosomal introns, and non-LTR retroelements-which together make up approximately half of the human genome. On a practical level, the properties of group II introns have been exploited to develop group II intron-based biotechnological tools. WIREs RNA 2016, 7:341-355. doi: 10.1002/wrna.1339 For further resources related to this article, please visit the WIREs website.

Alternative translation initiation codons for the plastid maturase MatK: unraveling the pseudogene misconception in the Orchidaceae

BACKGROUND

The plastid maturase MatK has been implicated as a possible model for the evolutionary "missing link" between prokaryotic and eukaryotic splicing machinery. This evolutionary implication has sparked investigations concerning the function of this unusual maturase. Intron targets of MatK activity suggest that this is an essential enzyme for plastid function. The matK gene, however, is described as a pseudogene in many photosynthetic orchid species due to presence of premature stop codons in translations, and its high rate of nucleotide and amino acid substitution.

RESULTS

Sequence analysis of the matK gene from orchids identified an out-of-frame alternative AUG initiation codon upstream from the consensus initiation codon used for translation in other angiosperms. We demonstrate translation from the alternative initiation codon generates a conserved MatK reading frame. We confirm that MatK protein is expressed and functions in sample orchids currently described as having a matK pseudogene using immunodetection and reverse-transcription methods. We demonstrate using phylogenetic analysis that this alternative initiation codon emerged de novo within the Orchidaceae, with several reversal events at the basal lineage and deep in orchid history.

CONCLUSION

These findings suggest a novel evolutionary shift for expression of matK in the Orchidaceae and support the function of MatK as a group II intron maturase in the plastid genome of land plants including the orchids.

Alternative translation initiation codons for the plastid maturase MatK: unraveling the pseudogene misconception in the Orchidaceae

BACKGROUND

The plastid maturase MatK has been implicated as a possible model for the evolutionary "missing link" between prokaryotic and eukaryotic splicing machinery. This evolutionary implication has sparked investigations concerning the function of this unusual maturase. Intron targets of MatK activity suggest that this is an essential enzyme for plastid function. The matK gene, however, is described as a pseudogene in many photosynthetic orchid species due to presence of premature stop codons in translations, and its high rate of nucleotide and amino acid substitution.

RESULTS

Sequence analysis of the matK gene from orchids identified an out-of-frame alternative AUG initiation codon upstream from the consensus initiation codon used for translation in other angiosperms. We demonstrate translation from the alternative initiation codon generates a conserved MatK reading frame. We confirm that MatK protein is expressed and functions in sample orchids currently described as having a matK pseudogene using immunodetection and reverse-transcription methods. We demonstrate using phylogenetic analysis that this alternative initiation codon emerged de novo within the Orchidaceae, with several reversal events at the basal lineage and deep in orchid history.

CONCLUSION

These findings suggest a novel evolutionary shift for expression of matK in the Orchidaceae and support the function of MatK as a group II intron maturase in the plastid genome of land plants including the orchids.

Complete chloroplast genome sequence of MD-2 pineapple and its comparative analysis among nine other plants from the subclass Commelinidae

BACKGROUND

Pineapple (Ananas comosus var. comosus) is known as the king of fruits for its crown and is the third most important tropical fruit after banana and citrus. The plant, which is indigenous to South America, is the most important species in the Bromeliaceae family and is largely traded for fresh fruit consumption. Here, we report the complete chloroplast sequence of the MD-2 pineapple that was sequenced using the PacBio sequencing technology.

RESULTS

In this study, the high error rate of PacBio long sequence reads of A. comosus's total genomic DNA were improved by leveraging on the high accuracy but short Illumina reads for error-correction via the latest error correction module from Novocraft. Error corrected long PacBio reads were assembled by using a single tool to produce a contig representing the pineapple chloroplast genome. The genome of 159,636 bp in length is featured with the conserved quadripartite structure of chloroplast containing a large single copy region (LSC) with a size of 87,482 bp, a small single copy region (SSC) with a size of 18,622 bp and two inverted repeat regions (IRA and IRB) each with the size of 26,766 bp. Overall, the genome contained 117 unique coding regions and 30 were repeated in the IR region with its genes contents, structure and arrangement similar to its sister taxon, Typha latifolia. A total of 35 repeats structure were detected in both the coding and non-coding regions with a majority being tandem repeats. In addition, 205 SSRs were detected in the genome with six protein-coding genes contained more than two SSRs. Comparative chloroplast genomes from the subclass Commelinidae revealed a conservative protein coding gene albeit located in a highly divergence region. Analysis of selection pressure on protein-coding genes using Ka/Ks ratio showed significant positive selection exerted on the rps7 gene of the pineapple chloroplast with P less than 0.05. Phylogenetic analysis confirmed the recent taxonomical relation among the member of commelinids which support the monophyly relationship between Arecales and Dasypogonaceae and between Zingiberales to the Poales, which includes the A. comosus.

CONCLUSIONS

The complete sequence of the chloroplast of pineapple provides insights to the divergence of genic chloroplast sequences from the members of the subclass Commelinidae. The complete pineapple chloroplast will serve as a reference for in-depth taxonomical studies in the Bromeliaceae family when more species under the family are sequenced in the future. The genetic sequence information will also make feasible other molecular applications of the pineapple chloroplast for plant genetic improvement.

Complete Chloroplast Genome of the Wollemi Pine (Wollemia nobilis): Structure and Evolution

The Wollemi pine (Wollemia nobilis) is a rare Southern conifer with striking morphological similarity to fossil pines. A small population of W. nobilis was discovered in 1994 in a remote canyon system in the Wollemi National Park (near Sydney, Australia). This population contains fewer than 100 individuals and is critically endangered. Previous genetic studies of the Wollemi pine have investigated its evolutionary relationship with other pines in the family Araucariaceae, and have suggested that the Wollemi pine genome contains little or no variation. However, these studies were performed prior to the widespread use of genome sequencing, and their conclusions were based on a limited fraction of the Wollemi pine genome. In this study, we address this problem by determining the entire sequence of the W. nobilis chloroplast genome. A detailed analysis of the structure of the genome is presented, and the evolution of the genome is inferred by comparison with the chloroplast sequences of other members of the Araucariaceae and the related family Podocarpaceae. Pairwise alignments of whole genome sequences, and the presence of unique pseudogenes, gene duplications and insertions in W. nobilis and Araucariaceae, indicate that the W. nobilis chloroplast genome is most similar to that of its sister taxon Agathis. However, the W. nobilis genome contains an unusually high number of repetitive sequences, and these could be used in future studies to investigate and conserve any remnant genetic diversity in the Wollemi pine.

Organellar maturases: A window into the evolution of the spliceosome

During the evolution of eukaryotic genomes, many genes have been interrupted by intervening sequences (introns) that must be removed post-transcriptionally from RNA precursors to form mRNAs ready for translation. The origin of nuclear introns is still under debate, but one hypothesis is that the spliceosome and the intron-exon structure of genes have evolved from bacterial-type group II introns that invaded the eukaryotic genomes. The group II introns were most likely introduced into the eukaryotic genome from an α-proteobacterial predecessor of mitochondria early during the endosymbiosis event. These self-splicing and mobile introns spread through the eukaryotic genome and later degenerated. Pieces of introns became part of the general splicing machinery we know today as the spliceosome. In addition, group II introns likely brought intron maturases with them to the nucleus. Maturases are found in most bacterial introns, where they act as highly specific splicing factors for group II introns. In the spliceosome, the core protein Prp8 shows homology to group II intron-encoded maturases. While maturases are entirely intron specific, their descendant of the spliceosomal machinery, the Prp8 protein, is an extremely versatile splicing factor with multiple interacting proteins and RNAs. How could such a general player in spliceosomal splicing evolve from the monospecific bacterial maturases? Analysis of the organellar splicing machinery in plants may give clues on the evolution of nuclear splicing. Plants encode various proteins which are closely related to bacterial maturases. The organellar genomes contain one maturase each, named MatK in chloroplasts and MatR in mitochondria. In addition, several maturase genes have been found in the nucleus as well, which are acting on mitochondrial pre-RNAs. All plant maturases show sequence deviation from their progenitor bacterial maturases, and interestingly are all acting on multiple organellar group II intron targets. Moreover, they seem to function in the splicing of group II introns together with a number of additional nuclear-encoded splicing factors, possibly acting as an organellar proto-spliceosome. Together, this makes them interesting models for the early evolution of nuclear spliceosomal splicing. In this review, we summarize recent advances in our understanding of the role of plant maturases and their accessory factors in plants. This article is part of a Special Issue entitled: Chloroplast Biogenesis.

Phylogeny of North American Tolypella (Charophyceae, Charophyta) based on plastid DNA sequences with a description of Tolypella ramosissima sp. nov

Characeae (Charophyceae, Charophyta) contains two tribes with six genera: tribe Chareae with four genera and tribe Nitelleae, which includes Tolypella and Nitella. This paper uses molecular and morphological data to elucidate the phylogeny of Tolypella species in North America. In the most comprehensive taxonomic treatment of Characeae, 16 Tolypella species worldwide were subsumed into two species, T. intricata and T. nidifica, in two sections, Rothia and Tolypella respectively. It was further suggested that Tolypella might be a derived group within Nitella. In this investigation into species diversity and relationships in North American Tolypella, sequence data from the plastid genes atpB, psbC, and rbcL were assembled for a broad range of charophycean and land plant taxa. Molecular data were used in conjunction with morphology to test monophyly of the genus and species within it. Phylogenetic analyses of the sequence data showed that Characeae is monophyletic but that Nitelleae is paraphyletic with Tolypella sister to a monophyletic Nitella + Chareae. The results also supported the monophyly of Tolypella and the sections Rothia and Tolypella. Morphologically defined species were supported as clades with little or no DNA sequence differences. In addition, molecular data revealed several lineages and a new species (T. ramosissima sp. nov.), which suggests greater species diversity in Tolypella than previously recognized.

Root microbiome relates to plant host evolution in maize and other Poaceae

Prokaryote-eukaryote interactions are primordial, but host selection of its bacterial community remains poorly understood. Because eukaryote evolution affects numerous traits shaping the ecology of their microbiome, we can expect that many evolutionary changes in the former will have the potential to impact on the composition of the latter. Consequently, the more phylogenetically distant the eukaryotic hosts, the more distinct their associated bacterial communities should be. We tested this with plants, by comparing the bacterial communities associated with maize genotypes or other Poaceae. 16S rRNA taxonomic microarray analysis showed that the genetic distance between rhizobacterial communities correlated significantly with the phylogenetic distance (derived from chloroplastic sequences) between Poaceae genotypes. This correlation was also significant when considering specific bacterial populations from all main bacterial divisions, instead of the whole rhizobacterial community. These results indicate that eukaryotic host's evolutionary history can be a significant factor shaping directly the assembly and composition of its associated bacterial compartment.

A comparison of 16 DNA regions for use as phylogenetic markers in the pleurocarpous moss genus Plagiothecium (Hypnales)

PREMISE OF THE STUDY

Within the Hypnales-the most derived and species-rich order of pleurocarpous mosses - phylogenies at or below the family level often show poor resolution. In preparation for a phylogeny of the genus Plagiothecium, we wished to identify the DNA markers best suited for evolutionary reconstruction in this group of hypnalean pleurocarps.

METHODS

For each of 25 collections of Plagiothecium and associated taxa, 16 DNA regions were sequenced: nuclear ITS and 26S, and plastid rps4, rps4-trnL, trnL-F, trnK (matK)-psbA, psbA-trnH, trnM-V, trnD-T, rbcL, atpB-rbcL, psbT-H, rpoC1 exon 2 (partial), the trnG intron, the rpl16 intron and the plastid ribosomal spacer DNA (cpITS). Each region was evaluated on the basis of its ability to resolve clades, the amount of homoplasy present in the data set, and the relative ease of obtaining the data. Descriptive statistics for each region are given.

KEY RESULTS

Under-utilized plastid markers for bryophytes such as trnK-psbA, rps4-trnL, and trnD-T outperformed more traditional markers such as trnL-F and rps4. Individual plastid topologies were similar, suggesting that only a limited amount of plastid data are needed to recover a backbone phylogeny. Adding a small amount of nuclear ribosomal data to a large plastid matrix restructured the recovered topology, emphasizing the importance of sampling multiple genomes and the need for new low-copy nuclear markers in bryophyte systematics.

CONCLUSIONS

Future genus-level phylogenies of pleurocarpous mosses should target under-utilized plastid markers such as trnK-psbA and rps4-trnL in conjunction with low-copy nuclear markers.

Group II intron splicing factors in plant mitochondria

Group II introns are large catalytic RNAs (ribozymes) which are found in bacteria and organellar genomes of several lower eukaryotes, but are particularly prevalent within the mitochondrial genomes (mtDNA) in plants, where they reside in numerous critical genes. Their excision is therefore essential for mitochondria biogenesis and respiratory functions, and is facilitated in vivo by various protein cofactors. Typical group II introns are classified as mobile genetic elements, consisting of the self-splicing ribozyme and its intron-encoded maturase protein. A hallmark of maturases is that they are intron specific, acting as cofactors which bind their own cognate containing pre-mRNAs to facilitate splicing. However, the plant organellar introns have diverged considerably from their bacterial ancestors, such as they lack many regions which are necessary for splicing and also lost their evolutionary related maturase ORFs. In fact, only a single maturase has been retained in the mtDNA of various angiosperms: the matR gene encoded in the fourth intron of the NADH-dehydrogenase subunit 1 (nad1 intron 4). Their degeneracy and the absence of cognate ORFs suggest that the splicing of plant mitochondria introns is assisted by trans-acting cofactors. Interestingly, in addition to MatR, the nuclear genomes of angiosperms also harbor four genes (nMat 1-4), which are closely related to maturases and contain N-terminal mitochondrial localization signals. Recently, we established the roles of two of these paralogs in Arabidopsis, nMAT1 and nMAT2, in the splicing of mitochondrial introns. In addition to the nMATs, genetic screens led to the identification of other genes encoding various factors, which are required for the splicing and processing of mitochondrial introns in plants. In this review we will summarize recent data on the splicing and processing of mitochondrial introns and their implication in plant development and physiology, with a focus on maturases and their accessory splicing cofactors.

The evolution of the plastid chromosome in land plants: gene content, gene order, gene function

This review bridges functional and evolutionary aspects of plastid chromosome architecture in land plants and their putative ancestors. We provide an overview on the structure and composition of the plastid genome of land plants as well as the functions of its genes in an explicit phylogenetic and evolutionary context. We will discuss the architecture of land plant plastid chromosomes, including gene content and synteny across land plants. Moreover, we will explore the functions and roles of plastid encoded genes in metabolism and their evolutionary importance regarding gene retention and conservation. We suggest that the slow mode at which the plastome typically evolves is likely to be influenced by a combination of different molecular mechanisms. These include the organization of plastid genes in operons, the usually uniparental mode of plastid inheritance, the activity of highly effective repair mechanisms as well as the rarity of plastid fusion. Nevertheless, structurally rearranged plastomes can be found in several unrelated lineages (e.g. ferns, Pinaceae, multiple angiosperm families). Rearrangements and gene losses seem to correlate with an unusual mode of plastid transmission, abundance of repeats, or a heterotrophic lifestyle (parasites or myco-heterotrophs). While only a few functional gene gains and more frequent gene losses have been inferred for land plants, the plastid Ndh complex is one example of multiple independent gene losses and will be discussed in detail. Patterns of ndh-gene loss and functional analyses indicate that these losses are usually found in plant groups with a certain degree of heterotrophy, might rendering plastid encoded Ndh1 subunits dispensable.

The evolution of chloroplast genome structure in ferns

The plastid genome (plastome) is a rich source of phylogenetic and other comparative data in plants. Most land plants possess a plastome of similar structure. However, in a major group of plants, the ferns, a unique plastome structure has evolved. The gene order in ferns has been explained by a series of genomic inversions relative to the plastome organization of seed plants. Here, we examine for the first time the structure of the plastome across fern phylogeny. We used a PCR-based strategy to map and partially sequence plastomes. We found that a pair of partially overlapping inversions in the region of the inverted repeat occurred in the common ancestor of most ferns. However, the ancestral (seed plant) structure is still found in early diverging branches leading to the osmundoid and filmy fern lineages. We found that a second pair of overlapping inversions occurred on a branch leading to the core leptosporangiates. We also found that the unique placement of the gene matK in ferns (lacking a flanking intron) is not a result of a large-scale inversion, as previously thought. This is because the intron loss maps to an earlier point on the phylogeny than the nearby inversion. We speculate on why inversions may occur in pairs and what this may mean for the dynamics of plastome evolution.

Molecular evolution and positive Darwinian selection of the chloroplast maturase matK

It is not clear whether matK evolves under Darwinian selection. In this study, 70 plant groups, representing 2,279 species at various evolutionary levels, were used to illustrate the molecular adaptation and evolutionary dynamics of gene divergence in matKs. Selective influences were investigated using standard dN/dS ratio methods. Analyses revealed the presence of positive selection in matKs of 32 plant groups. More positively selected sites were detected in the N-terminal region than in the RT domain and domain X of matK. Moreover, removing amino acid sites that are under positive selection has a significant effect on the bootstrap values of phylogenetic reconstruction. Our results suggest that the rapidly evolving matK evolves under positive selection in some lineages of land plants. Several regions of matK have experienced molecular adaptation, which fine-tunes maturase performance.

Chloroplast RNA metabolism

The chloroplast genome encodes proteins required for photosynthesis, gene expression, and other essential organellar functions. Derived from a cyanobacterial ancestor, the chloroplast combines prokaryotic and eukaryotic features of gene expression and is regulated by many nucleus-encoded proteins. This review covers four major chloroplast posttranscriptional processes: RNA processing, editing, splicing, and turnover. RNA processing includes the generation of transcript 5' and 3' termini, as well as the cleavage of polycistronic transcripts. Editing converts specific C residues to U and often changes the amino acid that is specified by the edited codon. Chloroplasts feature introns of groups I and II, which undergo protein-facilitated cis- or trans-splicing in vivo. Each of these RNA-based processes involves proteins of the pentatricopeptide motif-containing family, which does not occur in prokaryotes. Plant-specific RNA-binding proteins may underpin the adaptation of the chloroplast to the eukaryotic context.

Malpighiales phylogenetics: Gaining ground on one of the most recalcitrant clades in the angiosperm tree of life

The eudicot order Malpighiales contains ∼16000 species and is the most poorly resolved large rosid clade. To clarify phylogenetic relationships in the order, we used maximum likelihood, Bayesian, and parsimony analyses of DNA sequence data from 13 gene regions, totaling 15604 bp, and representing all three genomic compartments (i.e., plastid: atpB, matK, ndhF, and rbcL; mitochondrial: ccmB, cob, matR, nad1B-C, nad6, and rps3; and nuclear: 18S rDNA, PHYC, and newly developed low-copy EMB2765). Our sampling of 190 taxa includes representatives from all families of Malpighiales. These data provide greatly increased support for the recent additions of Aneulophus, Bhesa, Centroplacus, Ploiarium, and Rafflesiaceae to Malpighiales; sister relations of Phyllanthaceae + Picrodendraceae, monophyly of Hypericaceae, and polyphyly of Clusiaceae. Oxalidales + Huaceae, followed by Celastrales are successive sisters to Malpighiales. Parasitic Rafflesiaceae, which produce the world's largest flowers, are confirmed as embedded within a paraphyletic Euphorbiaceae. Novel findings show a well-supported placement of Ctenolophonaceae with Erythroxylaceae + Rhizophoraceae, sister-group relationships of Bhesa + Centroplacus, and the exclusion of Medusandra from Malpighiales. New taxonomic circumscriptions include the addition of Bhesa to Centroplacaceae, Medusandra to Peridiscaceae (Saxifragales), Calophyllaceae applied to Clusiaceae subfamily Kielmeyeroideae, Peraceae applied to Euphorbiaceae subfamily Peroideae, and Huaceae included in Oxalidales.

Conservation of selection on matK following an ancient loss of its flanking intron

The chloroplast gene trnK and its associated group II intron appear to be absent in a large and ancient clade that includes nearly 90% of fern species. However, the maturase protein encoded within the intron (matK) is still present and located on the boundary of a large-scale inversion. We surveyed the chloroplast genome sequence of clade-member Adiantum capillus-veneris for evidence of a still present but fragmented trnK intron. Lack of signature structural domains and sequence motifs in the genome indicate loss of the trnK intron through degradation in an ancestor of the clade. In plants, matK preferentially catalyzes splicing of the trnK intron, but may also have a generalist function, splicing other group II introns in the chloroplast genome. We therefore tested whether a shift in selective constraint has occurred after loss of the trnK intron. Using previously unavailable sequences for several ferns, we compared matK sequences of the intron-less fern clade to sequences from seed plants and ferns with the intron and found no significant differences in selection among lineages using multiple methods. We conclude that matK in ferns has maintained its apparently ancient and generalized function in chloroplasts, even after the loss of its co-evolved group II intron. Finally, we also present primers that will allow amplification and nucleotide sequencing of the phylogenetically useful matK gene in additional fern taxa.

Relationships and evolution of matK in a group of leafless orchids (Corallorhiza and Corallorhizinae; Orchidaceae: Epidendroideae)

Corallorhizinae are a small group of Old and New World temperate orchids of which a core monophyletic group comprises Govenia, Cremastra, Aplectrum, Oreorchis and the leafless Corallorhiza, and which according to phylogenetic analysis of nuclear ITS and plastid matK sequences, are related in this way: (Govenia (Cremastra (Aplectrum (Oreorchis (Corallorhiza))))). This hypothesis is consistent with the progressive deletion of the trnK intron and matK ORF. Frameshift-resulting indels yield a predicted loss of translation for the critical "domain X" region of matK and are evidence that matK is a probable pseudogene in Aplectrum, Oreorchis, and Corallorhiza. Within Corallorhiza, a previous hypothesis based on plastid DNA restriction site analysis is confirmed, with the thickened-labellum C. striata group being sister to the thin-labellum remainder of the genus, within which the circumboreal C. trifida is sister to the remainder, which then comprise two further sister groups: C. maculata + C. bulbosa + C. mertensiana and C. odontorhiza + C. wisteriana. A close relationship between C. striata and the recently described Appalachian C. bentleyi is shown; in particular, C. bentleyi is more closely allied to a southern Mexican population of C. striata than it is to northern North American C. striata populations, suggesting that two lineages, each with Mexican and northern North American populations, exist within the C. striata group.