Bacterial origin of a chloroplast intron: conserved self-splicing group I introns in cyanobacteria

Self-splicing introns in genes of Bastillevirinae bacteriophages

Group I introns are self-splicing ribozymes that can be found in eukaryotes, prokaryotes, and quite often in their viruses. The distribution, structure, and splicing of group I introns in genes of some phage taxa like the Tevenvirinae or Twortwirinae was extensively studied. On the other hand, the prevalence of intervening sequences in most other clades of bacterial viruses remains mostly unexplored. In this paper, we describe group I autocatalytic introns in genes of phages from the Bastillevirinae subfamily. This taxon belongs to the Herelleviridae family and consists of 15 genera and 37 species, including viruses with strong antimicrobial potential. A bioinformatic search for intron-related RNA structures revealed the presence of 45 intervening sequences within 37 genes that belong to four gene families. Eight of the nine genes selected for experimental validation were spliced-four only in an infected bacteria but additional four self-spliced in vitro. Interestingly, one of the studied genes undergoes alternative splicing. To sum up, our findings expand the knowledge on the distribution and diversity of group I introns and shed new light on this neglected aspect of phage transcriptomics. Additionally, in the course of our study, we demonstrated the effectiveness of nanopore sequencing in elucidating prokaryotic splicing mechanisms.

Group I Intron as a Potential Target for Antifungal Compounds: Development of a Trans-Splicing High-Throughput Screening Strategy

The search for safe and efficient new antifungal compounds for agriculture has led to more efforts in finding new modes of action. This involves the discovery of new molecular targets, including coding and non-coding RNA. Rarely found in plants and animals but present in fungi, group I introns are of interest as their complex tertiary structure may allow selective targeting using small molecules. In this work, we demonstrate that group I introns present in phytopathogenic fungi have a self-splicing activity in vitro that can be adapted in a high-throughput screening to find new antifungal compounds. Ten candidate introns from different filamentous fungi were tested and one group ID intron found in F. oxysporum showed high self-splicing efficiency in vitro. We designed the Fusarium intron to act as a trans-acting ribozyme and used a fluorescence-based reporter system to monitor its real time splicing activity. Together, these results are opening the way to study the druggability of such introns in crop pathogen and potentially discover small molecules selectively targeting group I introns in future high-throughput screenings.

Characterizing and circumventing sequence restrictions for synthesis of circular RNA in vitro

Just as eukaryotic circular RNA (circRNA) is a product of intracellular backsplicing, custom circRNA can be synthesized in vitro using a transcription template in which transposed halves of a split group I intron flank the sequence of the RNA to be circularized. Such permuted intron-exon (PIE) constructs have been used to produce circRNA versions of ribozymes, mimics of viral RNA motifs, a streptavidin aptamer, and protein expression vectors for genetic engineering and vaccine development. One limitation of this approach is the obligatory incorporation of small RNA segments (E1 and E2) into nascent circRNA at the site of end-joining. This restriction may preclude synthesis of small circRNA therapeutics and RNA nanoparticles that are sensitive to extraneous sequence, as well as larger circRNA mimics whose sequences must precisely match those of the native species on which they are modelled. In this work, we used serial mutagenesis and in vitro selection to determine how varying E1 and E2 sequences in a thymidylate synthase (td) group I intron PIE transcription template construct affects circRNA synthesis yield. Based on our collective findings, we present guidelines for the design of custom-tailored PIE transcription templates from which synthetic circRNAs of almost any sequence may be efficiently synthesized.

The function of KptA/Tpt1 gene - a minor review

Rapid response of uni- and multicellular organisms to environmental changes and their own growth is achieved through a series of molecular mechanisms, often involving modification of macromolecules, including nucleic acids, proteins and lipids. The ADP-ribosylation process has ability to modify these different macromolecules in cells, and is closely related to the biological processes, such as DNA replication, transcription, signal transduction, cell division, stress, microbial aging and pathogenesis. In addition, tRNA plays an essential role in the regulation of gene expression, as effector molecules, no-load tRNA affects the overall gene expression level of cells under some nutritional stress. KptA/Tpt1 is an essential phosphotransferase in the process of pre-tRNA splicing, releasing mature tRNA and participating in ADP-ribose. The objective of this review is concluding the gene structure, the evolution history and the function of KptA/Tpt1 from prokaryote to eukaryote organisms. At the same time, the results of promoter elements analysis were also shown in the present study. Moreover, the problems in the function of KptA/Tpt1 that have not been clarified at the present time are summarised, and some suggestions to solve those problems are given. This review presents no only a summary of clear function of KptA/Tpt1 in the process of tRNA splicing and ADP-ribosylation of organisms, but also gives some proposals to clarify unclear problems of it in the future.

Noncanonical features and modifications on the 5'-end of bacterial sRNAs and mRNAs

Although many eukaryotic transcripts contain cap structures, it has been long thought that bacterial RNAs do not carry any special modifications on their 5'-ends. In bacteria, primary transcripts are produced by transcription initiated with a nucleoside triphosphate and are therefore triphosphorylated on 5'-ends. Some transcripts are then processed by nucleases that yield monophosphorylated RNAs for specific cellular activities. Many primary transcripts are also converted to monophosphorylated species by removal of the terminal pyrophosphate for 5'-end-dependent degradation. Recent studies surprisingly revealed an expanded repertoire of chemical groups on 5'-ends of bacterial RNAs. In addition to mono- and triphosphorylated moieties, some mRNAs and sRNAs contain cap-like structures and diphosphates on their 5'-ends. Although incorporation and removal of these groups have become better understood in recent years, the physiological significance of these modifications remain obscure. This review highlights recent studies aimed at identification and elucidation of novel modifications on the 5'-ends of bacterial RNAs and discusses possible physiological applications of the modified RNAs. This article is categorized under: RNA Turnover and Surveillance > Regulation of RNA Stability RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry RNA Processing > Capping and 5' End Modifications.

Recent Insights Into the Structure, Function, and Evolution of the RNA-Splicing Endonucleases

RNA-splicing endonuclease (EndA) cleaves out introns from archaeal and eukaryotic precursor (pre)-tRNA and is essential for tRNA maturation. In archaeal EndA, the molecular mechanisms underlying complex assembly, substrate recognition, and catalysis have been well understood. Recently, certain studies have reported novel findings including the identification of new subunit types in archaeal EndA structures, providing insights into the mechanism underlying broad substrate specificity. Further, metagenomics analyses have enabled the acquisition of numerous DNA sequences of EndAs and intron-containing pre-tRNAs from various species, providing information regarding the co-evolution of substrate specificity of archaeal EndAs and tRNA genetic diversity, and the evolutionary pathway of archaeal and eukaryotic EndAs. Although the complex structure of the heterothermic form of eukaryotic EndAs is unknown, previous reports regarding their functions indicated that mutations in human EndA cause neurological disorders including pontocerebellar hypoplasia and progressive microcephaly, and yeast EndA significantly cleaves mitochondria-localized mRNA encoding cytochrome b mRNA processing 1 (Cpb1) for mRNA maturation. This mini-review summarizes the aforementioned results, discusses their implications, and offers my personal opinion regarding future directions for the analysis of the structure and function of EndAs.

The RNA-splicing endonuclease from the euryarchaeaon Methanopyrus kandleri is a heterotetramer with constrained substrate specificity

Four different types (α4, α'2, (αβ)2 and ϵ2) of RNA-splicing endonucleases (EndAs) for RNA processing are known to exist in the Archaea. Only the (αβ)2 and ϵ2 types can cleave non-canonical introns in precursor (pre)-tRNA. Both enzyme types possess an insert associated with a specific loop, allowing broad substrate specificity in the catalytic α units. Here, the hyperthermophilic euryarchaeon Methanopyrus kandleri (MKA) was predicted to harbor an (αβ)2-type EndA lacking the specific loop. To characterize MKA EndA enzymatic activity, we constructed a fusion protein derived from MKA α and β subunits (fMKA EndA). In vitro assessment demonstrated complete removal of the canonical bulge-helix-bulge (BHB) intron structure from MKA pre-tRNAAsn. However, removal of the relaxed BHB structure in MKA pre-tRNAGlu was inefficient compared to crenarchaeal (αβ)2 EndA, and the ability to process the relaxed intron within mini-helix RNA was not detected. fMKA EndA X-ray structure revealed a shape similar to that of other EndA types, with no specific loop. Mapping of EndA types and their specific loops and the tRNA gene diversity among various Archaea suggest that MKA EndA is evolutionarily related to other (αβ)2-type EndAs found in the Thaumarchaeota, Crenarchaeota and Aigarchaeota but uniquely represents constrained substrate specificity.

Transcription of the apicoplast genome

Many members of the Apicomplexa contain a remnant chloroplast, known as an apicoplast. The apicoplast encodes numerous genes, and loss of the organelle is lethal. Here, we present a summary of what is known about apicoplast transcription. Unlike plant chloroplasts, there is a single RNA polymerase, and initial transcription is polycistronic. RNA is then cleaved into tRNA, mRNA and rRNA molecules. Significant levels of antisense transcription have been reported, together with a single case of RNA editing. Polycistronic transcription is also observed in the related algae Chromera and Vitrella, which retain a photosynthetic chloroplast. Surprisingly, a polyU tail is added to Chromera and Vitrella transcripts which encode proteins involved in photosynthesis. No such tail is added to Plasmodium transcripts. Transcription in the Apicomplexa is remarkably similar to that seen in the chloroplast of the related peridinin dinoflagellate algae, reflecting the common evolutionary origins of the organelle.

The essential gene set of a photosynthetic organism

Synechococcus elongatus PCC 7942 is a model organism used for studying photosynthesis and the circadian clock, and it is being developed for the production of fuel, industrial chemicals, and pharmaceuticals. To identify a comprehensive set of genes and intergenic regions that impacts fitness in S. elongatus, we created a pooled library of ∼ 250,000 transposon mutants and used sequencing to identify the insertion locations. By analyzing the distribution and survival of these mutants, we identified 718 of the organism's 2,723 genes as essential for survival under laboratory conditions. The validity of the essential gene set is supported by its tight overlap with well-conserved genes and its enrichment for core biological processes. The differences noted between our dataset and these predictors of essentiality, however, have led to surprising biological insights. One such finding is that genes in a large portion of the TCA cycle are dispensable, suggesting that S. elongatus does not require a cyclic TCA process. Furthermore, the density of the transposon mutant library enabled individual and global statements about the essentiality of noncoding RNAs, regulatory elements, and other intergenic regions. In this way, a group I intron located in tRNA(Leu), which has been used extensively for phylogenetic studies, was shown here to be essential for the survival of S. elongatus. Our survey of essentiality for every locus in the S. elongatus genome serves as a powerful resource for understanding the organism's physiology and defines the essential gene set required for the growth of a photosynthetic organism.

Circularly permuted tRNA genes: their expression and implications for their physiological relevance and development

A number of genome analyses and searches using programs that focus on the RNA-specific bulge-helix-bulge (BHB) motif have uncovered a wide variety of disrupted tRNA genes. The results of these analyses have shown that genetic information encoding functional RNAs is described in the genome cryptically and is retrieved using various strategies. One such strategy is represented by circularly permuted tRNA genes, in which the sequences encoding the 5′-half and 3′-half of the specific tRNA are separated and inverted on the genome. Biochemical analyses have defined a processing pathway in which the termini of tRNA precursors (pre-tRNAs) are ligated to form a characteristic circular RNA intermediate, which is then cleaved at the acceptor-stem to generate the typical cloverleaf structure with functional termini. The sequences adjacent to the processing site located between the 3′-half and the 5′-half of pre-tRNAs potentially form a BHB motif, which is the dominant recognition site for the tRNA-intron splicing endonuclease, suggesting that circularization of pre-tRNAs depends on the splicing machinery. Some permuted tRNAs contain a BHB-mediated intron in their 5′- or 3′-half, meaning that removal of an intron, as well as swapping of the 5′- and 3′-halves, are required during maturation of their pre-tRNAs. To date, 34 permuted tRNA genes have been identified from six species of unicellular algae and one archaeon. Although their physiological significance and mechanism of development remain unclear, the splicing system of BHB motifs seems to have played a key role in the formation of permuted tRNA genes. In this review, current knowledge of circularly permuted tRNA genes is presented and some unanswered questions regarding these species are discussed.

Bacterial group I introns: mobile RNA catalysts

Group I introns are intervening sequences that have invaded tRNA, rRNA and protein coding genes in bacteria and their phages. The ability of group I introns to self-splice from their host transcripts, by acting as ribozymes, potentially renders their insertion into genes phenotypically neutral. Some group I introns are mobile genetic elements due to encoded homing endonuclease genes that function in DNA-based mobility pathways to promote spread to intronless alleles. Group I introns have a limited distribution among bacteria and the current assumption is that they are benign selfish elements, although some introns and homing endonucleases are a source of genetic novelty as they have been co-opted by host genomes to provide regulatory functions. Questions regarding the origin and maintenance of group I introns among the bacteria and phages are also addressed.

Mobile elements in a single-filament orange Guaymas Basin Beggiatoa ("Candidatus Maribeggiatoa") sp. draft genome: evidence for genetic exchange with cyanobacteria

The draft genome sequence of a single orange Beggiatoa ("Candidatus Maribeggiatoa") filament collected from a microbial mat at a hydrothermal site in Guaymas Basin (Gulf of California, Mexico) shows evidence of extensive genetic exchange with cyanobacteria, in particular for sensory and signal transduction genes. A putative homing endonuclease gene and group I intron within the 23S rRNA gene; several group II catalytic introns; GyrB and DnaE inteins, also encoding homing endonucleases; multiple copies of sequences similar to the fdxN excision elements XisH and XisI (required for heterocyst differentiation in some cyanobacteria); and multiple sequences related to an open reading frame (ORF) (00024_0693) of unknown function all have close non-Beggiatoaceae matches with cyanobacterial sequences. Sequences similar to the uncharacterized ORF and Xis elements are found in other Beggiatoaceae genomes, a variety of cyanobacteria, and a few phylogenetically dispersed pleiomorphic or filamentous bacteria. We speculate that elements shared among filamentous bacterial species may have been exchanged in microbial mats and that some of them may be involved in cell differentiation.

Insights into the physiology and ecology of the brackish-water-adapted Cyanobacterium Nodularia spumigena CCY9414 based on a genome-transcriptome analysis

Nodularia spumigena is a filamentous diazotrophic cyanobacterium that dominates the annual late summer cyanobacterial blooms in the Baltic Sea. But N. spumigena also is common in brackish water bodies worldwide, suggesting special adaptation allowing it to thrive at moderate salinities. A draft genome analysis of N. spumigena sp. CCY9414 yielded a single scaffold of 5,462,271 nucleotides in length on which genes for 5,294 proteins were annotated. A subsequent strand-specific transcriptome analysis identified more than 6,000 putative transcriptional start sites (TSS). Orphan TSSs located in intergenic regions led us to predict 764 non-coding RNAs, among them 70 copies of a possible retrotransposon and several potential RNA regulators, some of which are also present in other N2-fixing cyanobacteria. Approximately 4% of the total coding capacity is devoted to the production of secondary metabolites, among them the potent hepatotoxin nodularin, the linear spumigin and the cyclic nodulapeptin. The transcriptional complexity associated with genes involved in nitrogen fixation and heterocyst differentiation is considerably smaller compared to other Nostocales. In contrast, sophisticated systems exist for the uptake and assimilation of iron and phosphorus compounds, for the synthesis of compatible solutes, and for the formation of gas vesicles, required for the active control of buoyancy. Hence, the annotation and interpretation of this sequence provides a vast array of clues into the genomic underpinnings of the physiology of this cyanobacterium and indicates in particular a competitive edge of N. spumigena in nutrient-limited brackish water ecosystems.

Invasion of protein coding genes by green algal ribosomal group I introns

The spread of group I introns depends on their association with intron-encoded homing endonucleases. Introns that encode functional homing endonuclease genes (HEGs) are highly invasive, whereas introns that only encode the group I ribozyme responsible for self-splicing are generally stably inherited (i.e., vertical inheritance). A number of recent case studies have provided new knowledge on the evolution of group I introns, however, there are still large gaps in understanding of their distribution on the tree of life, and how they have spread into new hosts and genic sites. During a larger phylogenetic survey of chlorophyceaen green algae, we found that 23 isolates contain at least one group I intron in the rbcL chloroplast gene. Structural analyses show that the introns belong to one of two intron lineages, group IA2 intron-HEG (GIY-YIG family) elements inserted after position 462 in the rbcL gene, and group IA1 introns inserted after position 699. The latter intron type sometimes encodes HNH homing endonucleases. The distribution of introns was analyzed on an exon phylogeny and patterns were recovered that are consistent with vertical inheritance and possible horizontal transfer. The rbcL 462 introns are thus far reported only within the Volvocales, Hydrodictyaceae and Bracteacoccus, and closely related isolates of algae differ in the presence of rbcL introns. Phylogenetic analysis of the intron conserved regions indicates that the rbcL699 and rbcL462 introns have distinct evolutionary origins. The rbcL699 introns were likely derived from ribosomal RNA L2449 introns, whereas the rbcL462 introns form a close relationship with psbA introns.

Circular RNA interference effector molecules (WO10084371)

Just over a decade after the discovery of RNA interference (RNAi), the RNAi field has begun travelling the bumpy road towards the clinic. Notwithstanding this extraordinarily rapid progress, and despite holding great promise for substantial future clinical impact, RNAi-inducing agents unfortunately exhibit certain physico-chemical and pharmacokinetic drawbacks. The patent application WO10084371 utilizes the exquisite biological mechanism of auto-catalytic intron splicing to generate circularized interfering RNAs (ciRNAs) in order to remedy the issue of rapid nuclease-mediated degradation of RNAi-inducers. The present patent evaluation assesses the utility of the ciRNAs in light of commonly used nucleotide modification strategies for modulating the properties of siRNAs, as well as scrutinizes the experimental data substantiating the patent application in question. The ciRNAs disclosed in WO10084371 exhibit potency on par with exogenously introduced short hairpin RNAs, albeit displaying increased exonuclease resistance. However, experimental validation as to the exact silencing mechanism is sparse, although the potential Dicer-substrate function of the ciRNAs is indeed a promising aspect. Despite the novel, elegant approach of WO10084371, the current gold standard nucleotide modifications appear to deliver sufficient stability for therapeutic applications, meaning that solving the issue of siRNA delivery still remains the major hurdle for clinical success.

A conserved lysine residue in the crenarchaea-specific loop is important for the crenarchaeal splicing endonuclease activity

In Archaea, splicing endonuclease (EndA) recognizes and cleaves precursor RNAs to remove introns. Currently, EndAs are classified into three families according to their subunit structures: homotetramer, homodimer, and heterotetramer. The crenarchaeal heterotetrameric EndAs can be further classified into two subfamilies based on the size of the structural subunit. Subfamily A possesses a structural subunit similar in size to the catalytic subunit, whereas subfamily B possesses a structural subunit significantly smaller than the catalytic subunit. Previously, we solved the crystal structure of an EndA from Pyrobaculum aerophilum. The endonuclease was classified into subfamily B, and the structure revealed that the enzyme lacks an N-terminal subdomain in the structural subunit. However, no structural information is available for crenarchaeal heterotetrameric EndAs that are predicted to belong to subfamily A. Here, we report the crystal structure of the EndA from Aeropyrum pernix, which is predicted to belong to subfamily A. The enzyme possesses the N-terminal subdomain in the structural subunit, revealing that the two subfamilies of heterotetrameric EndAs are structurally distinct. EndA from A. pernix also possesses an extra loop region that is characteristic of crenarchaeal EndAs. Our mutational study revealed that the conserved lysine residue in the loop is important for endonuclease activity. Furthermore, the sequence characteristics of the loops and the positions towards the substrate RNA according to a docking model prompted us to propose that crenarchaea-specific loops and an extra amino acid sequence at the catalytic loop of nanoarchaeal EndA are derived by independent convergent evolution and function for recognizing noncanonical bulge-helix-bulge motif RNAs as substrates.

Preparation of group I introns for biochemical studies and crystallization assays by native affinity purification

The study of functional RNAs of various sizes and structures requires efficient methods for their synthesis and purification. Here, 23 group I intron variants ranging in length from 246 to 341 nucleotides—some containing exons—were subjected to a native purification technique previously applied only to shorter RNAs (<160 nucleotides). For the RNAs containing both exons, we adjusted the original purification protocol to allow for purification of radiolabeled molecules. The resulting RNAs were used in folding assays on native gel electrophoresis and in self-splicing assays. The intron-only RNAs were subjected to the regular native purification scheme, assayed for folding and employed in crystallization screens. All RNAs that contained a 3′ overhang of one nucleotide were efficiently cleaved off from the support and were at least 90% pure after the non-denaturing purification. A representative subset of these RNAs was shown to be folded and self-splicing after purification. Additionally, crystals were grown for a 286 nucleotide long variant of the Clostridium botulinum intron. These results demonstrate the suitability of the native affinity purification method for the preparation of group I introns. We hope these findings will stimulate a broader application of this strategy to the preparation of other large RNA molecules.

A natural ribozyme with 3',5' RNA ligase activity

Using electrophoresis, sequencing and enzymatic digestion, we show that the group I intron from the cyanobacterium Anabaena sp. PCC 7120 catalyzes phosphodiester bond formation using a triphosphate on the 5'-terminal nucleotide, much like protein polymerases and engineered ribozymes. In the process, this ribozyme forms a unique circular RNA that incorporates the exogenous guanosine cofactor added during self-splicing. This finding may have relevance to a prebiotic RNA world and to modern biology.

Evidence from glycine transfer RNA of a frozen accident at the dawn of the genetic code

BACKGROUND

Transfer RNA (tRNA) is the means by which the cell translates DNA sequence into protein according to the rules of the genetic code. A credible proposition is that tRNA was formed from the duplication of an RNA hairpin half the length of the contemporary tRNA molecule, with the point at which the hairpins were joined marked by the canonical intron insertion position found today within tRNA genes. If these hairpins possessed a 3'-CCA terminus with different combinations of stem nucleotides (the ancestral operational RNA code), specific aminoacylation and perhaps participation in some form of noncoded protein synthesis might have occurred. However, the identity of the first tRNA and the initial steps in the origin of the genetic code remain elusive.

RESULTS

Here we show evidence that glycine tRNA was the first tRNA, as revealed by a vestigial imprint in the anticodon loop sequences of contemporary descendents. This provides a plausible mechanism for the missing first step in the origin of the genetic code. In 448 of 466 glycine tRNA gene sequences from bacteria, archaea and eukaryote cytoplasm analyzed, CCA occurs immediately upstream of the canonical intron insertion position, suggesting the first anticodon (NCC for glycine) has been captured from the 3'-terminal CCA of one of the interacting hairpins as a result of an ancestral ligation.

CONCLUSION

That this imprint (including the second and third nucleotides of the glycine tRNA anticodon) has been retained through billions of years of evolution suggests Crick's 'frozen accident' hypothesis has validity for at least this very first step at the dawn of the genetic code.

REVIEWERS

This article was reviewed by Dr Eugene V. Koonin, Dr Rob Knight and Dr David H Ardell.

The complete nucleotide sequence of the cassava (Manihot esculenta) chloroplast genome and the evolution of atpF in Malpighiales: RNA editing and multiple losses of a group II intron

The complete sequence of the chloroplast genome of cassava (Manihot esculenta, Euphorbiaceae) has been determined. The genome is 161,453 bp in length and includes a pair of inverted repeats (IR) of 26,954 bp. The genome includes 128 genes; 96 are single copy and 16 are duplicated in the IR. There are four rRNA genes and 30 distinct tRNAs, seven of which are duplicated in the IR. The infA gene is absent; expansion of IRb has duplicated 62 amino acids at the 3' end of rps19 and a number of coding regions have large insertions or deletions, including insertions within the 23S rRNA gene. There are 17 intron-containing genes in cassava, 15 of which have a single intron while two (clpP, ycf3) have two introns. The usually conserved atpF group II intron is absent and this is the first report of its loss from land plant chloroplast genomes. The phylogenetic distribution of the atpF intron loss was determined by a PCR survey of 251 taxa representing 34 families of Malpighiales and 16 taxa from closely related rosids. The atpF intron is not only missing in cassava but also from closely related Euphorbiaceae and other Malpighiales, suggesting that there have been at least seven independent losses. In cassava and all other sequenced Malphigiales, atpF gene sequences showed a strong association between C-to-T substitutions at nucleotide position 92 and the loss of the intron, suggesting that recombination between an edited mRNA and the atpF gene may be a possible mechanism for the intron loss.

A CRM domain protein functions dually in group I and group II intron splicing in land plant chloroplasts

The CRM domain is a recently recognized RNA binding domain found in three group II intron splicing factors in chloroplasts, in a bacterial protein that associates with ribosome precursors, and in a family of uncharacterized proteins in plants. To elucidate the functional repertoire of proteins with CRM domains, we studied CFM2 (for CRM Family Member 2), which harbors four CRM domains. RNA coimmunoprecipitation assays showed that CFM2 in maize (Zea mays) chloroplasts is associated with the group I intron in pre-trnL-UAA and group II introns in the ndhA and ycf3 pre-mRNAs. T-DNA insertions in the Arabidopsis thaliana ortholog condition a defective-seed phenotype (strong allele) or chlorophyll-deficient seedlings with impaired splicing of the trnL group I intron and the ndhA, ycf3-int1, and clpP-int2 group II introns (weak alleles). CFM2 and two previously described CRM proteins are bound simultaneously to the ndhA and ycf3-int1 introns and act in a nonredundant fashion to promote their splicing. With these findings, CRM domain proteins are implicated in the activities of three classes of catalytic RNA: group I introns, group II introns, and 23S rRNA.

The molecular evolution and structural organization of group I introns at position 1389 in nuclear small subunit rDNA of myxomycetes

The number of nuclear group I introns from myxomycetes is rapidly increasing in GenBank as more rDNA sequences from these organisms are being sequenced. They represent an interesting and complex group of intervening sequences because several introns are mobile (or inferred to be mobile) and many contain large and unusual insertions in peripheral loops. Here we describe related group I introns at position 1389 in the small subunit rDNA of representatives from the myxomycete family Didymiaceae. Phylogenetic analyses support a common origin and mainly vertical inheritance of the intron. All S1389 introns from the Didymiaceae belong to the IC1 subclass of nuclear group I introns. The central catalytic core region of about 100 nt appears divergent in sequence composition even though the introns reside in closely related species. Furthermore, unlike the majority of group I introns from myxomycetes the S1389 introns do not self-splice as naked RNA in vitro under standard conditions, consistent with a dependence on host factors for folding or activity. Finally, the myxomycete S1389 introns are exclusively found within the family Didymiaceae, which suggests that this group I intron was acquired after the split between the families Didymiaceae and Physaraceae.

Group I introns and gnra tetraloops: remnants of ‘The RNA world’?

GNRA tetraloops, found in high frequency in natural RNAs, make loop-receptor interactions, stabilizing the tertiary structure of Group I introns, a class of small RNAs. Analyzing 230 Group I introns, to study the distribution and sequence pattern of the GNRA tetraloops, we suggest that these features reflect the ancestral nature of these catalytic molecules, in a prebiotic RNA world. The adenosine rich GNRA tetraloops would have interacted with each other through long range RNA-RNA interactions to form higher order structures forming potential sites that render the propensity for the short RNAs to bind to metal ions from the prebiotic pool, aiding them to act as metalloenzymes.

Secretory phospholipase A2 hydrolysis of phospholipid analogues is dependent on water accessibility to the active site

A new and unnatural type of phospholipids with the head group attached to the 2-position of the glycerol backbone has been synthesized and shown to be a good substrate for secretory phospholipase A2 (sPLA2). To investigate the unexpected sPLA2 activity, we have compared three different phospholipids by using fluorescence techniques and HPLC, namely: (R)-1,2-dipalmitoyl-glycero-3-phosphocholine (hereafter referred to as 1R), (R)-1-O-hexadecyl-2-palmitoyl-glycero-3-phoshocholine (2R), and (S)-1-O-hexadecyl-3-palmitoyl-glycero-2-phosphocholine (3S). Furthermore, to understand the underlying mechanisms for the observed differences, we have performed molecular dynamics simulations to clarify on a structural level the substrate specificity of sPLA2 toward phospholipid analogues with their head groups in the 2-position of the glycerol backbone. We have studied the lipids above 1R, 2R, and 3S as well as their enantiomers 1S, 2S, and 3R. In the simulations of sPLA2-1S and sPLA2-3R, structural distortion in the binding cleft induced by the phospholipids showed that these are not substrates for sPLA2. In the case of the phospholipids 1R, 2R, and 3S, our simulations revealed that the difference observed experimentally in sPLA2 activity might be caused by reduced access of water molecules to the active site. We have monitored the number of water molecules that enter the active site region for the different sPLA2-phospholipid complexes and found that the probability of a water molecule reaching the correct position such that hydrolysis can occur is reduced for the unnatural lipids. The relative water count follows 1R > 2R > 3S. This is in good agreement with experimental data that indicate the same trend for sPLA2 activity: 1R > 2R > 3S.

The restriction fold turns to the dark side: a bacterial homing endonuclease with a PD-(D/E)-XK motif

The homing endonuclease I-Ssp6803I causes the insertion of a group I intron into a bacterial tRNA gene-the only example of an invasive mobile intron within a bacterial genome. Using a computational fold prediction, mutagenic screen and crystal structure determination, we demonstrate that this protein is a tetrameric PD-(D/E)-XK endonuclease - a fold normally used to protect a bacterial genome from invading DNA through the action of restriction endonucleases. I-Ssp6803I uses its tetrameric assembly to promote recognition of a single long target site, whereas restriction endonuclease tetramers facilitate cooperative binding and cleavage of two short sites. The limited use of the PD-(D/E)-XK nucleases by mobile introns stands in contrast to their frequent use of LAGLIDADG and HNH endonucleases - which in turn, are rarely incorporated into restriction/modification systems.

Pseudocyphellaria crocata, P. neglecta and P. perpetua from the Northern and Southern Hemispheres are a phylogenetic species and share cyanobionts

Pseudocyphellaria crocata, P. neglecta and P. perpetua specimens were examined to investigate links between genetic variation and morphology, geographical distribution and cyanobiont specificity. Fungal internal transcribed spacer (ITS), beta-tubulin and cyanobacterial tRNA(Leu) (UAA) intron sequences were used to investigate symbiont diversity in these lichens. Specimens were morphologically distinct but could not be distinguished by ITS sequences. Phylogenetic analyses split the P. crocata specimens into two clades, the larger of which contained P. neglecta and P. perpetua. Five cyanobionts were identified; two of these were in a number of specimens, while three were each restricted to a single lichen thallus. Fungus-specific molecular markers indicated that all specimens belonged to a single phylogenetic species. However, this may contain a cryptic species. Geography was linked to genetic diversity with Canadian specimens forming a monophyletic group, and most Southern Hemisphere specimens grouping together, although Chile represented a hot spot of genetic diversity. There was no connection between fungal genetic diversity and cyanobiont choice, consistent with the presence of a common pool of cyanobionts.

Molecular evolution of the trnL(UAA) intron in bryophytes

Structure, variability, and molecular evolution of the trnL(UAA) intron in bryophytes (mosses and liverworts) is analyzed based on more than 1000 sequences representing all classes, including comparisons of lengths and GC-contents, sequence similarities, evolutionary rates and ti/tv ratios of the major lineages and selected genera. Secondary structure analyses of the more variable stem-loop regions facilitated recognition of sequence repeats and minute inversions that often occurred independently in non-related lineages, thus supporting alignment construction and homology assessment. The most length-variable stem-loop region P8 does not share a common evolutionary history across all major bryophyte lineages. Independent nucleotide additions such as internally repeated sequence segments resulted in non-homologous P8 sequences that cannot be folded into a common P8 secondary structure, neither for all bryophytes nor for liverworts or mosses. To address evolutionary patterns, separate analyses of P6/P8 and the remaining intron (core) have to be performed, as overall values of the complete intron are misleading. It is argued that a transition bias observed above the genus level in the core structure is caused by structural constraints, not by its higher GC-content in comparison to the more AT-rich P6 and P8. Compensating base pair changes detected in highly conserved elements are often characteristic of the major bryophyte lineages (classes). Sequence divergence and evolutionary rates are generally higher in liverworts than in mosses, resulting in ambiguous alignments of P6 and P8 even within classes. In mosses, trends towards length reduction of P8 and lower evolutionary rates of the intron are observed. Average intraspecific variation is less than 1%, corresponding to 2-3 mutations in the complete intron.

The chloroplast trnT-trnF region in the seed plant lineage Gnetales

The trnT-trnF region is located in the large single-copy region of the chloroplast genome. It consists of the trnL intron, a group I intron, and the trnT-trnL and trnL-trnF intergenic spacers. We analyzed the evolution of the region in the three genera of the gymnosperm lineage Gnetales (Gnetum, Welwitschia, and Ephedra), with especially dense sampling in Gnetum for which we sequenced 41 accessions, representing most of the 25-35 species. The trnL intron has a conserved secondary structure and contains elements that are homologous across land plants, while the spacers are so variable in length and composition that homology cannot be found even among the three genera. Palindromic sequences that form hairpin structures were detected in the trnL-trnF spacer, but neither spacer contained promoter elements for the tRNA genes. The absence of promoters, presence of hairpin structures in the trnL-trnF spacer, and high sequence variation in both spacers together suggest that trnT and trnF are independently transcribed. Our model for the expression and processing of the genes tRNA(Thr)(UGU), tRNA(Leu)(UAA), and tRNA(Phe) (GAA) therefore attributes the seemingly neutral evolution of the two spacers to their escape from functional constraints.

Evolution of the trnF(GAA) Gene in Arabidopsis Relatives and the Brassicaceae Family: Monophyletic Origin and Subsequent Diversification of a Plastidic Pseudogene

Recently, we used the 5'-trnL(UAA)-trnF(GAA) region of the chloroplast DNA for phylogeographic reconstructions and phylogenetic analysis among the genera Arabidopsis, Boechera, Rorippa, Nasturtium, and Cardamine. Despite the fact that extensive gene duplications are rare among the chloroplast genome of higher plants, within these taxa the anticodon domain of the trnF(GAA) gene exhibit extensive gene duplications with one to eight tandemly repeated copies in close 5' proximity of the functional gene. Interestingly, even in Arabidopsis thaliana we found six putative pseudogenic copies of the functional trnF gene within the 5'-intergenic trnL-trnF spacer. A reexamination of trnL(UAA)-trnF(GAA) regions from numerous published phylogenetic studies among halimolobine, cardaminoid, and other cruciferous taxa revealed not only extensive trnF gene duplications but also favor the hypothesis about a single origin of trnF pseudogene formation during evolution of the Brassicaceae family 16-21 MYA. Conserved sequence motifs from this tandemly repeated region are codistributed nonrandomly throughout the plastome, and we found some similarities with a DNA sequence duplication in the rps7 gene and its adjacent spacer. Our results demonstrate the potential evolutionary dynamics of a plastidic region generally regarded as highly conserved and probably cotranscribed and, as shown here for several genera among cruciferous plants, greatly characterized by parallel gains and losses of duplicated trnF copies.

Analysis of 2'-phosphotransferase (Tpt1p) from Saccharomyces cerevisiae: evidence for a conserved two-step reaction mechanism

Tpt1p is an essential protein responsible for the 2'-phosphotransferase step of tRNA splicing in Saccharomyces cerevisiae, in which the splice junction 2'-phosphate of ligated tRNA is transferred to NAD to form mature tRNA and ADP-ribose 1''-2'' cyclic phosphate. We showed previously that Tpt1p is a member of a family of functional 2'-phosphotransferases found in eukaryotes, eubacteria, and archaea, that the Escherichia coli protein (KptA) is highly specific for 2'-phosphorylated RNAs despite the lack of obvious natural substrates, and that KptA acts on a trinucleotide substrate through an intermediate in which RNA is ADP-ribosylated at the 2'-phosphate. This mechanism is similar to a proposed mechanism of NAD-dependent histone deacetylases. We present evidence here that this mechanism is conserved in S. cerevisiae, and we identify residues important for the second step of the reaction, during which the intermediate is resolved into products. We examined 21 Tpt1 protein variants mutated in conserved residues or blocks of residues and show that one of them, Tpt1 K69A/R71S protein, accumulates large amounts of intermediate with trinucleotide substrate due to a very slow second step. This intermediate can be trapped on beads when formed with biotin-NAD. We also show that Tpt1 K69A/R71S protein forms an intermediate with the natural ligated tRNA substrate and demonstrate that, as expected, this mutation is lethal in yeast. The high degree of conservation of these residues suggests that the entire Tpt1p family is involved in a similar two-step chemical reaction.

A self-splicing group I intron in DNA polymerase genes of T7-like bacteriophages

Group I introns are inserted into genes of a wide variety of bacteriophages of gram-positive bacteria. However, among the phages of enteric and other gram-negative proteobacteria, introns have been encountered only in phage T4 and several of its close relatives. Here we report the insertion of a self-splicing group I intron in the coding sequence of the DNA polymerase genes of PhiI and W31, phages that are closely related to T7. The introns belong to subgroup IA2 and both contain an open reading frame, inserted into structural element P6a, encoding a protein belonging to the HNH family of homing endonucleases. The introns splice efficiently in vivo and self-splice in vitro under mild conditions of ionic strength and temperature. We conclude that there is no barrier for maintenance of group I introns in phages of proteobacteria.

Repeat-type distribution in trnL intron does not correspond with species phylogeny: comparison of the genetic markers 16S rRNA and trnL intron in heterocystous cyanobacteria

tRNA(Leu) UAA (trnL) intron sequences are used as genetic markers for differentiating cyanobacteria and for constructing phylogenies, since the introns are thought to be more variable among close relatives than is the 16S rRNA gene, the conventional phylogenetic marker. The evolution of trnL intron sequences and their utility as a phylogenetic marker were analysed among heterocystous cyanobacteria with maximum-parsimony, maximum-likelihood and Bayesian inference by comparing their evolutionary information to that of the 16S rRNA gene. Trees inferred from the 16S rRNA gene and the distribution of two repeat classes in the P6b stem-loop of the trnL intron were in clear conflict. The results show that, while similar heptanucleotide repeat classes I and II in the P6b stem-loop of the trnL intron could be found among distant relatives, some close relatives harboured different repeat classes with a high sequence difference. Moreover, heptanucleotide repeat class II and other sequences from the P6b stem-loop of the trnL intron interrupted several other intergenic regions in the genomes of heterocystous cyanobacteria. Cluster analyses based on conserved intron sequences without loops P6b, P9 and parts of P5 corresponded in most clades to the 16S rRNA gene phylogeny, although the relationships were not resolved well, according to low bootstrap support. Thus, the hypervariable loop sequences of the trnL intron, especially the P6b stem-loop, cannot be used for phylogenetic analysis and conclusions cannot be drawn about species relationships on the basis of these elements. Evolutionary scenarios are discussed considering the origin of the repeats.

Distinct sites of phosphorothioate substitution interfere with folding and splicing of the Anabaena group I intron

Although the active site of group I introns is phylogenetically conserved, subclasses of introns have evolved different mechanisms of stabilizing the catalytic core. Large introns contain weakly conserved 'peripheral' domains that buttress the core through predicted interhelical contacts, while smaller introns use loop-helix interactions for stability. In all cases, specific and non-specific magnesium ion binding accompanies folding into the active structure. Whether similar RNA-RNA and RNA-magnesium ion contacts play related functional roles in different introns is not clear, particularly since it can be difficult to distinguish interactions directly involved in catalysis from those important for RNA folding. Using phosphorothioate interference with RNA activity and structure in the small (249 nt) group I intron from Anabaena, we used two independent assays to detect backbone phosphates important for catalysis and those involved in intron folding. Comparison of the interference sites identified in each assay shows that positions affecting catalysis cluster primarily in the conserved core of the intron, consistent with conservation of functionally important phosphates, many of which are magnesium ion binding sites, in diverse group I introns, including those from Azoarcus and Tetrahymena. However, unique sites of folding interference located outside the catalytic core imply that different group I introns, even within the same subclass, use distinct sets of tertiary interactions to stabilize the structure of the catalytic core.

Recent horizontal intron transfer to a chloroplast genome

Evidence is presented for the recent, horizontal transfer of a self-splicing, homing group II intron from a cyanobacteria to the chloroplast genome of Euglena myxocylindracea. The psbA gene of E.myxocylindracea was found to contain a single 2566 nt group II intron with a gene in domain 4 for a 575 amino acid maturase. The predicted secondary structure and tertiary interactions of the group II intron, as well as the derived maturase primary sequence, most closely resemble the homing intron of the cyanobacterium Calothrix and the rnl introns of Porphyra purpurea mitochondria, while being only distantly related to all other Euglena plastid introns and maturases. All main functional domains of the intron-encoded proteins of known homing introns are conserved, including reverse transcriptase domains 1-7, the zinc finger domain and domain X. The close relationship with cyanobacterial introns was confirmed by phylogenetic analysis. Both the full-length psbA intron and a Delta-maturase variant self-splice in vitro in two independent assays. The psbA intron is the first example of a self-splicing chloroplast group II intron from any organism. These results support the conclusion that the psbA intron is the result of a recent horizontal transfer into the E.myxocylindracea chloroplast genome from a cyanobacterial donor and should prompt a reconsideration of horizontal transfer mechanisms to account for the origin of other chloroplast genetic elements.

Chloroplast RNA processing and stability

Primary chloroplast transcripts are processed in a number of ways, including intron splicing, internal cleavage of polycistronic RNAs, and endonucleolytic or exonucleolytic cleavages at the transcript termini. All chloroplast RNAs are also subject to degradation, although a curious feature of many chloroplast mRNAs is their relative longevity. Some of these processes, e.g., psbA splicing and stability of a number of chloroplast mRNAs, are regulated in response to light-dark cycles or nutrient availability. This review highlights recent advances in our understanding of these processes in the model organism Chlamydomonas reinhardtii, focusing on results since the extensive reviews published in 1998 [Herrin DL et al. 1998 (pp. 183-195), Nickelsen Y 1998 (pp. 151-163), Stern DB and Drager RG 1998 (pp. 164-182), in Rochaix JD et al. (eds) The Molecular Biology of Chloroplasts and Mitochondria in Chlamydomonas. Kluwer Academic Publishers, Dordrecht, The Netherlands]. We also allude to studies with other organisms, and to the potential impact of the Chlamydomonas genome project where appropriate.

Lariat formation and a hydrolytic pathway in plant chloroplast group II intron splicing

Lariat formation has been studied intensively only with a few self-splicing group II introns, and little is known about how the numerous diverse introns in plant organelles are excised. Several of these introns have branch-points that are not a single bulge but are adjoined by A:A, A:C, A:G and G:G pairs. Using a highly sensitive in vivo approach, we demonstrate that all but one of the barley chloroplast introns splice via the common pathway that produces a branched product. RNA editing does not improve domain 5 and 6 structures of these introns. The conserved branch-point in tobacco rpl16 is chosen even if an adjacent unpaired adenosine is available, suggesting that spatial arrangements in domain 6 determine correct branch-point selection. Lariats were not detected for the chloroplast trnV intron, which lacks an unpaired adenosine in domain 6. Instead, this intron is released as linear molecules that undergo further polyadenylation. trnV, which is conserved throughout plant evolution, constitutes the first example of naturally occurring hydrolytic group II intron splicing in vivo.

Species of cyanolichens from Pseudocyphellaria with indistinguishable ITS sequences have different photobionts

•  Cyanobacteria were isolated from bipartite cyanolichen species of Pseudocyphellaria and the identity of the major photobionts established. The specificity of the cyanobacterial-fungal association was also examined. •  Comparison of 16S rRNA gene sequences distinguished cyanobacterial and green algal isolates, and both 16S rRNA gene and tRNA^Leu  (UAA) intron sequences of isolates and lichen thalli identified candidate photobionts. In addition, the genetic diversity of both the cyanobiont and mycobiont was investigated using the comparison of tRNA^Leu  (UAA) intron sequences and ITS sequences, respectively. •  The 16S rRNA gene sequences identified two species-specific photobionts with similar sequences; however, the tRNA^Leu  (UAA) intron sequences unambiguously discriminated between the two symbiotic cyanobacterial strains. Moreover, the fungal ITS sequences of the two corresponding lichens, Pseudocyphellaria crocata and Pseudocyphellaria neglecta, showed little variation. •  The cyanobacterial-fungal associations of P. crocata and P. neglecta were specific for all samples. However, the similarity of the ITS sequences raised the possibility that they represent the same species and that their different morphology is influenced by the cyanobacterial symbiont.

Group I self-splicing intron in the recA gene of Bacillus anthracis

Self-splicing introns are rarely found in bacteria and bacteriophages. They are classified into group I and II according to their structural features and splicing mechanisms. While the group I introns are occasionally found in protein-coding regions of phage genomes and in several tRNA genes of cyanobacteria and proteobacteria, they had not been found in protein-coding regions of bacterial genomes. Here we report a group I intron in the recA gene of Bacillus anthracis which was initially found by DNA sequencing as an intervening sequence (IVS). By using reverse transcriptase PCR, the IVS was shown to be removable from the recA precursor mRNA for RecA that was being translated in E. coli. The splicing was visualized in vitro with labeled free GTP, indicating that it is a group I intron, which is also implied by its predicted secondary structure. The RecA protein of B. anthracis expressed in E. coli was functional in its ability to complement a recA defect. When recA-negative E. coli cells were irradiated with UV, the Bacillus RecA reduced the UV susceptibility of the recA mutant, regardless of the presence of intron.

The cyanobacterial tRNA(Leu) (UAA) intron: evolutionary patterns in a genetic marker

Cyanobacterial tRNA(Leu) (UAA) intron sequences from natural populations of Nostoc and other cyanobacteria were compared. Variation between the different introns was not randomly distributed but strongly restricted by the secondary and tertiary structure of the intron. Although all Nostoc sequences examined shared high similarity, differences were observed in one stem-loop. This stem-loop could be divided into two classes, both built up from two base pairing heptanucleotide repeats. Size variation was primarily caused by different numbers of repeats, but some strains also contained additional sequences in this stem-loop not following the heptanucleotide repeat motif. Several sequences showing similarity with these additional sequences were identified in the Nostoc punctiforme genome. Furthermore, the regions flanking these sequences contained the same, or similar, heptanucleotide repeats as those flanking the corresponding sequences in the intron. It is proposed that both slipped strand mispairing during replication and homologous recombination among different loci in the genome are important processes causing variation between introns.

A study of the phylogeny of Brassica rapa, B. nigra, Raphanus sativus, and their related genera using noncoding regions of chloroplast DNA

There are two evolutionary lineages in the genus Brassica: the rapa/oleracea lineage and the nigra lineage. Using nuclear DNA sequences such as the intergenic spacer between 5S rRNA genes and the internal transcribed spacer between 18S and 25S rRNA genes, we and others had previously demonstrated that Raphanus sativus is closely related to the nigra lineage. In the present study, we sequenced the chloroplast noncoding region between trnT and trnF and that between trnD and trnT in seven species and showed that R. sativus is more closely related to the rapa/oleracea lineage than to the nigra lineage. The conflicting results from nuclear DNA and chloroplast DNA support the hypothesis that Raphanus was derived from a hybridization between the rapa/oleracea and the nigra lineages. We estimated the date of this hybridization event to be 60% of the divergence time between the two Brassica lineages. In addition, the pattern and rate of nucleotide substitution were studied. There were more transversions than transitions in these noncoding regions, which have a high AT content. Furthermore, the proportion of transversions among the substitutions at a site increases with increasing A + T content of its two adjacent nucleotides. An influence of immediate 5(') pyrimidine on substitution pattern is also observed when both adjacent bases in the two DNA strands are A or T. The rate of nucleotide substitution in the trnL group I intron is only about one third of the rate in the nearby intergenic spacers in the trnT-trnF fragment. The rate of nucleotide substitution in the rapa/oleracea lineage is at least 1.5 times that in the nigra lineage.

Novel group I intron in the tRNA(Leu)(UAA) gene of a gamma-proteobacterium isolated from a deep subsurface environment

A group I intron has been found to interrupt the anticodon loop of the tRNA(Leu)(UAA) gene in a bacterium belonging to the gamma-subdivision of Proteobacteria and isolated from a deep subsurface environment. The subsurface isolate SMCC D0715 was identified as belonging to the genus Pseudomonas. The group I intron from this isolate is the first to be reported for gamma-proteobacteria, and the first instance of a tRNA(Leu)(UAA) group I intron to be found in a group of bacteria other than cyanobacteria. The 231-nucleotide (nt) intron's sequence has group I conserved elements and folds into a bona fide group I secondary structure with canonical base-paired segments P1 to P9 and a paired region, P10. The D0715 intron possesses the 11-nt motif CCUACG. UAUGG in its P8 region, a feature not common in bacterial introns. To date, phylogenetic analysis has shown that bacterial introns form two distinct families, and their complex distribution suggests that both lateral transfer and common ancestry have taken part in the evolutionary history of these elements.