U1 small nuclear RNA and spliceosomal introns in Euglena gracilis

Comprehensive analysis of the Kinetoplastea intron landscape reveals a novel intron-containing gene and the first exclusively trans-splicing eukaryote

BACKGROUND

In trypanosomatids, a group of unicellular eukaryotes that includes numerous important human parasites, cis-splicing has been previously reported for only two genes: a poly(A) polymerase and an RNA helicase. Conversely, trans-splicing, which involves the attachment of a spliced leader sequence, is observed for nearly every protein-coding transcript. So far, our understanding of splicing in this protistan group has stemmed from the analysis of only a few medically relevant species. In this study, we used an extensive dataset encompassing all described trypanosomatid genera to investigate the distribution of intron-containing genes and the evolution of splice sites.

RESULTS

We identified a new conserved intron-containing gene encoding an RNA-binding protein that is universally present in Kinetoplastea. We show that Perkinsela sp., a kinetoplastid endosymbiont of Amoebozoa, represents the first eukaryote completely devoid of cis-splicing, yet still preserving trans-splicing. We also provided evidence for reverse transcriptase-mediated intron loss in Kinetoplastea, extensive conservation of 5' splice sites, and the presence of non-coding RNAs within a subset of retained trypanosomatid introns.

CONCLUSIONS

All three intron-containing genes identified in Kinetoplastea encode RNA-interacting proteins, with a potential to fine-tune the expression of multiple genes, thus challenging the perception of cis-splicing in these protists as a mere evolutionary relic. We suggest that there is a selective pressure to retain cis-splicing in trypanosomatids and that this is likely associated with overall control of mRNA processing. Our study provides new insights into the evolution of introns and, consequently, the regulation of gene expression in eukaryotes.

Unique features of conventional and nonconventional introns in Euglena gracilis

BACKGROUND

Nuclear introns in Euglenida have been understudied. This study aimed to investigate nuclear introns in Euglenida by identifying a large number of introns in Euglena gracilis (E. gracilis), including cis-spliced conventional and nonconventional introns, as well as trans-spliced outrons. We also examined the sequence characteristics of these introns.

RESULTS

A total of 28,337 introns and 11,921 outrons were identified. Conventional and nonconventional introns have distinct splice site features; the former harbour canonical GT/C-AG splice sites, whereas the latter are capable of forming structured motifs with their terminal sequences. We observed that short introns had a preference for canonical GT-AG introns. Notably, conventional introns and outrons in E. gracilis exhibited a distinct cytidine-rich polypyrimidine tract, in contrast to the thymidine-rich tracts observed in other organisms. Furthermore, the SL-RNAs in E. gracilis, as well as in other trans-splicing species, can form a recently discovered motif called the extended U6/5' ss duplex with the respective U6s. We also describe a novel type of alternative splicing pattern in E. gracilis. The tandem repeat sequences of introns in this protist were determined, and their contents were comparable to those in humans.

CONCLUSIONS

Our findings highlight the unique features of E. gracilis introns and provide insights into the splicing mechanism of these introns, as well as the genomics and evolution of Euglenida.

Transcriptomic and genomic identification of spliceosomal genes from Euglena gracilis

Diverse splicing types in nuclear and chloroplast genes of protist Euglena gracilis have been recognized for decades. However, the splicing machinery responsible for processing nuclear precursor messenger RNA introns, including trans-splicing of the 5' terminal outron and spliced leader (SL) RNA, remains elusive. Here, we identify 166 spliceosomal protein genes and two snRNA genes from E. gracilis by performing bioinformatics analysis from a combination of next-generation and full-length transcriptomic RNA sequencing (RNAseq) data as well as draft genomic data. With the spliceosomal proteins we identified in hand, the insensitivity of E. gracilis to some splicing modulators is revealed at the sequence level. The prevalence of SL RNA-mediated trans-splicing is estimated to be more than 70% from our full-length RNAseq data. Finally, the splicing proteomes between E. gracilis and its three evolutionary cousins within the same Excavata group are compared. In conclusion, our study characterizes the spliceosomal components in E. gracilis and provides the molecular basis for further exploration of underlying splicing mechanisms.

Heterotrophic euglenid Rhabdomonas costata resembles its phototrophic relatives in many aspects of molecular and cell biology

Euglenids represent a group of protists with diverse modes of feeding. To date, only a partial genomic sequence of Euglena gracilis and transcriptomes of several phototrophic and secondarily osmotrophic species are available, while primarily heterotrophic euglenids are seriously undersampled. In this work, we begin to fill this gap by presenting genomic and transcriptomic drafts of a primary osmotroph, Rhabdomonas costata. The current genomic assembly length of 100 Mbp is 14× smaller than that of E. gracilis. Despite being too fragmented for comprehensive gene prediction it provided fragments of the mitochondrial genome and comparison of the transcriptomic and genomic data revealed features of its introns, including several candidates for nonconventional types. A set of 39,456 putative R. costata proteins was predicted from the transcriptome. Annotation of the mitochondrial core metabolism provides the first data on the facultatively anaerobic mitochondrion of R. costata, which in most respects resembles the mitochondrion of E. gracilis with a certain level of streamlining. R. costata can synthetise thiamine by enzymes of heterogenous provenances and haem by a mitochondrial-cytoplasmic C4 pathway with enzymes orthologous to those found in E. gracilis. The low percentage of green algae-affiliated genes supports the ancestrally osmotrophic status of this species.

A New Type of Circular RNA derived from Nonconventional Introns in Nuclear Genes of Euglenids

Nuclear protein-coding genes of euglenids (Discoba, Euglenozoa, Euglenida) contain conventional (spliceosomal) and nonconventional introns. The latter have been found only in euglenozoans. A unique feature of nonconventional introns is the ability to form a stable and slightly conserved RNA secondary structure bringing together intron ends and placing adjacent exons in proximity. To date, little is known about the mechanism of their excision (e.g. whether it involves the spliceosome or not). The tubA gene of Euglena gracilis harbors three conventional and three nonconventional introns. While the conventional introns are excised as lariats, nonconventional introns are present in the cell solely as circular RNAs with full-length ends. Based on this discovery as well as on previous observations indicating that nonconventional introns are observed frequently at unique positions of genes, we suggest that this new type of intronic circRNA might play a role in intron mobility.

Comparative molecular cell biology of phototrophic euglenids and parasitic trypanosomatids sheds light on the ancestor of Euglenozoa

Parasitic trypanosomatids and phototrophic euglenids are among the most extensively studied euglenozoans. The phototrophic euglenid lineage arose relatively recently through secondary endosymbiosis between a phagotrophic euglenid and a prasinophyte green alga that evolved into the euglenid secondary chloroplast. The parasitic trypanosomatids (i.e. Trypanosoma spp. and Leishmania spp.) and the freshwater phototrophic euglenids (i.e. Euglena gracilis) are the most evolutionary distant lineages in the Euglenozoa phylogenetic tree. The molecular and cell biological traits they share can thus be considered as ancestral traits originating in the common euglenozoan ancestor. These euglenozoan ancestral traits include common mitochondrial presequence motifs, respiratory chain complexes containing various unique subunits, a unique ATP synthase structure, the absence of mitochondria-encoded transfer RNAs (tRNAs), a nucleus with a centrally positioned nucleolus, closed mitosis without dissolution of the nuclear membrane and nucleoli, a nuclear genome containing the unusual 'J' base (β-D-glucosyl-hydroxymethyluracil), processing of nucleus-encoded precursor messenger RNAs (pre-mRNAs) via spliced-leader RNA (SL-RNA) trans-splicing, post-transcriptional gene silencing by the RNA interference (RNAi) pathway and the absence of transcriptional regulation of nuclear gene expression. Mitochondrial uridine insertion/deletion RNA editing directed by guide RNAs (gRNAs) evolved in the ancestor of the kinetoplastid lineage. The evolutionary origin of other molecular features known to be present only in either kinetoplastids (i.e. polycistronic transcripts, compaction of nuclear genomes) or euglenids (i.e. monocistronic transcripts, huge genomes, many nuclear cis-spliced introns, polyproteins) is unclear.

Ploeotids Represent Much of the Phylogenetic Diversity of Euglenids

Ploeotids are an assemblage of rigid phagotrophic euglenids that have 10-12 pellicular strips and glide on their posterior flagellum. Molecular phylogenies place them as a poorly resolved, likely paraphyletic assemblage outside the Spirocuta clade of flexible euglenids, which includes the well-known phototrophs and primary osmotrophs. Here, we report SSU rRNA gene sequences from 38 ploeotids, using both single-cell and culture-based methods. Several contain group I or non-canonical introns. Our phylogenetic analyses place ploeotids in 8 distinct clades: Olkasia n. gen., Hemiolia n. gen., Liburna n. gen., Lentomonas, Decastava, Keelungia, Ploeotiidae, and Entosiphon. Ploeotia vitrea, the type of Ploeotia, is closely related to P. oblonga and Serpenomonas costata, but not to Lentomonas. Ploeotia cf. vitrea sensu Lax and Simpson 2013 is not related to P. vitrea and has a different pellicle strip architecture (as imaged by scanning electron microscopy): it instead represents a novel genus and species, Olkasia polycarbonata. We also describe new genera, Hemiolia and Liburna, for the morphospecies Anisonema trepidum and A. glaciale. A recent system proposing 13 suprafamilial taxa that include ploeotids is not supported by our phylogenies. The exact relationships between ploeotid groups remain unresolved and multigene phylogenetics or phylogenomics are needed to address this uncertainty.

Order of removal of conventional and nonconventional introns from nuclear transcripts of Euglena gracilis

Nuclear genes of euglenids and marine diplonemids harbor atypical, nonconventional introns which are not observed in the genomes of other eukaryotes. Nonconventional introns do not have the conserved borders characteristic for spliceosomal introns or the sequence complementary to U1 snRNA at the 5' end. They form a stable secondary structure bringing together both exon/intron junctions, nevertheless, this conformation does not resemble the form of self-splicing or tRNA introns. In the genes studied so far, frequent nonconventional introns insertions at new positions have been observed, whereas conventional introns have been either found at the conserved positions, or simply lost. In this work, we examined the order of intron removal from Euglena gracilis transcripts of the tubA and gapC genes, which contain two types of introns: nonconventional and spliceosomal. The relative order of intron excision was compared for pairs of introns belonging to different types. Furthermore, intermediate products of splicing were analyzed using the PacBio Next Generation Sequencing system. The analysis led to the main conclusion that nonconventional introns are removed in a rapid way but later than spliceosomal introns. Moreover, the observed accumulation of transcripts with conventional introns removed and nonconventional present may suggest the existence of a time gap between the two types of splicing.

The existence of conventional spliceosomal introns in genes of eukaryotic organisms is a well-known theorem. However, genes of the unicellular algae group, euglenids, contain also another type of introns, so-called nonconventional ones. They lack canonical borders, a feature most characteristic for conventional introns and form a stable secondary structure bringing together their ends. Along with the increasing popularity of whole genome studies, nonconventional introns were also disclosed in the genes of other protists, diplonemids. In this study we were particularly interested which introns–conventional or nonconventional–are removed earlier from euglenids’ pre-mRNA. To track this process we analyzed transcripts of the two Euglena gracilis genes. The relative order of intron excision was compared for pairs of introns belonging to different types. We also surveyed thousands of intermediate products of splicing using the Next-Generation Sequencing system. Summarizing the results of both experiments, we proved that spliceosomal introns are removed at an earlier stage of pre-mRNA maturation than nonconventional ones.

Intermediate introns in nuclear genes of euglenids - are they a distinct type?

Background

Nuclear genes of euglenids contain two major types of introns: conventional spliceosomal and nonconventional introns. The latter are characterized by variable non-canonical borders, RNA secondary structure that brings intron ends together, and an unknown mechanism of removal. Some researchers also distinguish intermediate introns, which combine features of both types. They form a stable RNA secondary structure and are classified into two subtypes depending on whether they contain one (intermediate/nonconventional subtype) or both (conventional/intermediate subtype) canonical spliceosomal borders. However, it has been also postulated that most introns classified as intermediate could simply be special cases of conventional or nonconventional introns.

Results

Sequences of tubB, hsp90 and gapC genes from six strains of Euglena agilis were obtained. They contain four, six, and two or three introns, respectively (the third intron in the gapC gene is unique for just one strain). Conventional introns were present at three positions: two in the tubB gene (at one position conventional/intermediate introns were also found) and one in the gapC gene. Nonconventional introns are present at ten positions: two in the tubB gene (at one position intermediate/nonconventional introns were also found), six in hsp90 (at four positions intermediate/nonconventional introns were also found), and two in the gapC gene.

Conclusions

Sequence and RNA secondary structure analyses of nonconventional introns confirmed that their most strongly conserved elements are base pairing nucleotides at positions +4, +5 and +6/ -8, −7 and −6 (in most introns CAG/CTG nucleotides were observed). It was also confirmed that the presence of the 5' GT/C end in intermediate/nonconventional introns is not the result of kinship with conventional introns, but is due to evolutionary pressure to preserve the purine at the 5' end. However, an example of a nonconventional intron with GC-AG ends was shown, suggesting the possibility of intron type conversion between nonconventional and conventional. Furthermore, an analysis of conventional introns revealed that the ability to form a stable RNA secondary structure by some introns is probably not a result of their relationship with nonconventional introns. It was also shown that acquisition of new nonconventional introns is an ongoing process and can be observed at the level of a single species. In the recently acquired intron in the gapC gene an extended direct repeats at the intron-exon junctions are present, suggesting that double-strand break repair process could be the source of new nonconventional introns.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0620-5) contains supplementary material, which is available to authorized users.

Distribution of conventional and nonconventional introns in tubulin (α and β) genes of euglenids

The nuclear genomes of euglenids contain three types of introns: conventional spliceosomal introns, nonconventional introns for which a splicing mechanism is unknown (variable noncanonical borders, RNA secondary structure bringing together intron ends), and so-called intermediate introns, which combine features of conventional and nonconventional introns. Analysis of two genes, tubA and tubB, from 20 species of euglenids reveals contrasting distribution patterns of conventional and nonconventional introns--positions of conventional introns are conserved, whereas those of the nonconventional ones are unique to individual species or small groups of closely related taxa. Moreover, in the group of phototrophic euglenids, 11 events of conventional intron loss versus 15 events of nonconventional intron gain were identified. A comparison of all nonconventional intron sequences highlighted the most conserved elements in their sequence and secondary structure. Our results led us to put forward two hypotheses. 1) The first one posits that mutational changes in intron sequence could lead to a change in their excision mechanism--intermediate introns would then be a transitional form between the conventional and nonconventional introns. 2) The second hypothesis concerns the origin of nonconventional introns--because of the presence of inverted repeats near their ends, insertion of MITE-like transposon elements is proposed as a possible source of new introns.

Genomic characterization of Neoparamoeba pemaquidensis (Amoebozoa) and its kinetoplastid endosymbiont

We have performed a genomic characterization of a kinetoplastid protist living within the amoebozoan Neoparamoeba pemaquidensis. The genome of this "Ichthyobodo-related organism" was found to be unexpectedly large, with at least 11 chromosomes between 1.0 and 3.5 Mbp and a total genome size of at least 25 Mbp.

A possible role for short introns in the acquisition of stroma-targeting peptides in the flagellate Euglena gracilis

The chloroplasts of Euglena gracilis bounded by three membranes arose via secondary endosymbiosis of a green alga in a heterotrophic euglenozoan host. Many genes were transferred from symbiont to the host nucleus. A subset of Euglena nuclear genes of predominately symbiont, but also host, or other origin have obtained complex presequences required for chloroplast targeting. This study has revealed the presence of short introns (41–93 bp) either in the second half of presequence-encoding regions or shortly downstream of them in nine nucleus-encoded E. gracilis genes for chloroplast proteins (Eno29, GapA, PetA, PetF, PetJ, PsaF, PsbM, PsbO, and PsbW). In addition, the E. gracilis Pbgd gene contains two introns in the second half of presequence-encoding region and one at the border of presequence-mature peptide-encoding region. Ten of 12 introns present within presequence-encoding regions or shortly downstream of them identified in this study have typical eukaryotic GT/AG borders, are T-rich, 45–50 bp long, and pairwise sequence identities range from 27 to 61%. Thus single recombination events might have been mediated via these cis-spliced introns. A double crossing over between these cis-spliced introns and trans-spliced introns present in 5′-UTRs of Euglena nuclear genes is also likely to have occurred. Thus introns and exon-shuffling could have had an important role in the acquisition of chloroplast targeting signals in E. gracilis. The results are consistent with a late origin of photosynthetic euglenids.

U3 snoRNA genes are multi-copy and frequently linked to U5 snRNA genes in Euglena gracilis

Background

U3 snoRNA is a box C/D small nucleolar RNA (snoRNA) involved in the processing events that liberate 18S rRNA from the ribosomal RNA precursor (pre-rRNA). Although U3 snoRNA is present in all eukaryotic organisms, most investigations of it have focused on fungi (particularly yeasts), animals and plants. Relatively little is known about U3 snoRNA and its gene(s) in the phylogenetically broad assemblage of protists (mostly unicellular eukaryotes). In the euglenozoon Euglena gracilis, a distant relative of the kinetoplastid protozoa, Southern analysis had previously revealed at least 13 bands hybridizing with U3 snoRNA, suggesting the existence of multiple copies of U3 snoRNA genes.

Results

Through screening of a λ genomic library and PCR amplification, we recovered 14 U3 snoRNA gene variants, defined by sequence heterogeneities that are mostly located in the U3 3'-stem-loop domain. We identified three different genomic arrangements of Euglena U3 snoRNA genes: i) stand-alone, ii) linked to tRNA^Arg genes, and iii) linked to a U5 snRNA gene. In arrangement ii), the U3 snoRNA gene is positioned upstream of two identical tRNA^Arg genes that are convergently transcribed relative to the U3 gene. This scenario is reminiscent of a U3 snoRNA-tRNA gene linkage previously described in trypanosomatids. We document here twelve different U3 snoRNA-U5 snRNA gene arrangements in Euglena; in each case, the U3 gene is linked to a downstream and convergently oriented U5 gene, with the intergenic region differing in length and sequence among the variants.

Conclusion

The multiple U3 snoRNA-U5 snRNA gene linkages, which cluster into distinct families based on sequence similarities within the intergenic spacer, presumably arose by genome, chromosome, and/or locus duplications. We discuss possible reasons for the existence of the unusually large number of U3 snoRNA genes in the Euglena genome. Variability in the signal intensities of the multiple Southern hybridization bands raises the possibility that Euglena contains a naturally aneuploid chromosome complement.

Phylogeny of Phagotrophic Euglenids (Euglenozoa) as Inferred from Hsp90 Gene Sequences

Molecular phylogenies of euglenids are usually based on ribosomal RNA genes that do not resolve the branching order among the deeper lineages. We addressed deep euglenid phylogeny using the cytosolic form of the heat-shock protein 90 gene (hsp90), which has already been employed with some success in other groups of euglenozoans and eukaryotes in general. Hsp90 sequences were generated from three taxa of euglenids representing different degrees of ultrastructural complexity, namely Petalomonas cantuscygni and wild isolates of Entosiphon sulcatum, and Peranema trichophorum. The hsp90 gene sequence of P. trichophorum contained three short introns (ranging from 27 to 31 bp), two of which had non-canonical borders GG-GG and GG-TG and two 10-bp inverted repeats, suggesting a structure similar to that of the non-canonical introns described in Euglena gracilis. Phylogenetic analyses confirmed a closer relationship between kinetoplastids and diplonemids than to euglenids, and supported previous views regarding the branching order among primarily bacteriovorous, primarily eukaryovorous, and photosynthetic euglenids. The position of P. cantuscygni within Euglenozoa, as well as the relative support for the nodes including it were strongly dependent on outgroup selection. The results were most consistent when the jakobid Reclinomonas americana was used as the outgroup. The most robust phylogenies place P. cantuscygni as the most basal branch within the euglenid clade. However, the presence of a kinetoplast-like mitochondrial inclusion in P. cantuscygni deviates from the currently accepted apomorphy-based definition of the kinetoplastid clade and highlights the necessity of detailed studies addressing the molecular nature of the euglenid and diplonemid mitochondrial genome.

Studies of DNA dumbbells VIII. Melting analysis of DNA dumbbells with dinucleotide repeat stem sequences

Melting curves and circular dichroism spectra were measured for a number of DNA dumbbell and linear molecules containing dinucleotide repeat sequences of different lengths. To study effects of different sequences on the melting and spectroscopic properties, six DNA dumbbells whose stems contain the central sequences (AA)(10), (AC)(10), (AG)(10), (AT)(10), (GC)(10), and (GG)(10) were prepared. These represent the minimal set of 10 possible dinucleotide repeats. To study effects of dinucleotide repeat length, dumbbells with the central sequences (AG)(n), n = 5 and 20, were prepared. Control molecules, dumbbells with a random central sequence, (RN)(n), n = 5, 10, and 20, were also prepared. The central sequence of each dumbbell was flanked on both sides by the same 12 base pairs and T(4) end-loops. Melting curves were measured by optical absorbance and differential scanning calorimetry in solvents containing 25, 55, 85, and 115 mM Na(+). CD spectra were collected from 20 to 45 degrees C and [Na(+)] from 25 to 115 mM. The spectral database did not reveal any apparent temperature dependence in the pretransition region. Analysis of the melting thermodynamics evaluated as a function of Na(+) provided a means for quantitatively estimating the counterion release with melting for the different sequences. Results show a very definite sequence dependence, indicating the salt-dependent properties of duplex DNA are also sequence dependent. Linear DNA molecules containing the (AG)(n) and (RN)(n), sequences, n = 5, 10, 20, and 30, were also prepared and studied. The linear DNA molecules had the exact sequences of the dumbbell stems. That is, the central repeat sequence in each linear duplex was flanked on both sides by the same 12-bp sequence. Melting and CD studies were also performed on the linear DNA molecules. Comparison of results obtained for the same sequences in dumbbell and linear molecular environments reveals several interesting features of the interplay between sequence-dependent structural variability, sequence length, and the unconstrained (linear) or constrained (dumbbell) molecular environments.

A high density of ancient spliceosomal introns in oxymonad excavates

Background

Certain eukaryotic genomes, such as those of the amitochondriate parasites Giardia and Trichomonas, have very low intron densities, so low that canonical spliceosomal introns have only recently been discovered through genome sequencing. These organisms were formerly thought to be ancient eukaryotes that diverged before introns originated, or at least became common. Now however, they are thought to be members of a supergroup known as excavates, whose members generally appear to have low densities of canonical introns. Here we have used environmental expressed sequence tag (EST) sequencing to identify 17 genes from the uncultivable oxymonad Streblomastix strix, to survey intron densities in this most poorly studied excavate group.

Results

We find that Streblomastix genes contain an unexpectedly high intron density of about 1.1 introns per gene. Moreover, over 50% of these are at positions shared between a broad spectrum of eukaryotes, suggesting theyare very ancient introns, potentially present in the last common ancestor of eukaryotes.

Conclusion

The Streblomastix data show that the genome of the ancestor of excavates likely contained many introns and the subsequent evolution of introns has proceeded very differently in different excavate lineages: in Streblomastix there has been much stasis while in Trichomonas and Giardia most introns have been lost.

The evolution of spliceosomal introns: patterns, puzzles and progress

The origins and importance of spliceosomal introns comprise one of the longest-abiding mysteries of molecular evolution. Considerable debate remains over several aspects of the evolution of spliceosomal introns, including the timing of intron origin and proliferation, the mechanisms by which introns are lost and gained, and the forces that have shaped intron evolution. Recent important progress has been made in each of these areas. Patterns of intron-position correspondence between widely diverged eukaryotic species have provided insights into the origins of the vast differences in intron number between eukaryotic species, and studies of specific cases of intron loss and gain have led to progress in understanding the underlying molecular mechanisms and the forces that control intron evolution.

An ancient spliceosomal intron in the ribosomal protein L7a gene (Rpl7a) of Giardia lamblia

Background

Only one spliceosomal-type intron has previously been identified in the unicellular eukaryotic parasite, Giardia lamblia (a diplomonad). This intron is only 35 nucleotides in length and is unusual in possessing a non-canonical 5' intron boundary sequence, CT, instead of GT.

Results

We have identified a second spliceosomal-type intron in G. lamblia, in the ribosomal protein L7a gene (Rpl7a), that possesses a canonical GT 5' intron boundary sequence. A comparison of the two known Giardia intron sequences revealed extensive nucleotide identity at both the 5' and 3' intron boundaries, similar to the conserved sequence motifs recently identified at the boundaries of spliceosomal-type introns in Trichomonas vaginalis (a parabasalid). Based on these observations, we searched the partial G. lamblia genome sequence for these conserved features and identified a third spliceosomal intron, in an unassigned open reading frame. Our comprehensive analysis of the Rpl7a intron in other eukaryotic taxa demonstrates that it is evolutionarily conserved and is an ancient eukaryotic intron.

Conclusion

An analysis of the phylogenetic distribution and properties of the Rpl7a intron suggests its utility as a phylogenetic marker to evaluate particular eukaryotic groupings. Additionally, analysis of the G. lamblia introns has provided further insight into some of the conserved and unique features possessed by the recently identified spliceosomal introns in related organisms such as T. vaginalis and Carpediemonas membranifera.

Unusual features of fibrillarin cDNA and gene structure in Euglena gracilis: evolutionary conservation of core proteins and structural predictions for methylation-guide box C/D snoRNPs throughout the domain Eucarya

Box C/D ribonucleoprotein (RNP) particles mediate O^2′-methylation of rRNA and other cellular RNA species. In higher eukaryotic taxa, these RNPs are more complex than their archaeal counterparts, containing four core protein components (Snu13p, Nop56p, Nop58p and fibrillarin) compared with three in Archaea. This increase in complexity raises questions about the evolutionary emergence of the eukaryote-specific proteins and structural conservation in these RNPs throughout the eukaryotic domain. In protists, the primarily unicellular organisms comprising the bulk of eukaryotic diversity, the protein composition of box C/D RNPs has not yet been extensively explored. This study describes the complete gene, cDNA and protein sequences of the fibrillarin homolog from the protozoon Euglena gracilis, the first such information to be obtained for a nucleolus-localized protein in this organism. The E.gracilis fibrillarin gene contains a mixture of intron types exhibiting markedly different sizes. In contrast to most other E.gracilis mRNAs characterized to date, the fibrillarin mRNA lacks a spliced leader (SL) sequence. The predicted fibrillarin protein sequence itself is unusual in that it contains a glycine-lysine (GK)-rich domain at its N-terminus rather than the glycine-arginine-rich (GAR) domain found in most other eukaryotic fibrillarins. In an evolutionarily diverse collection of protists that includes E.gracilis, we have also identified putative homologs of the other core protein components of box C/D RNPs, thereby providing evidence that the protein composition seen in the higher eukaryotic complexes was established very early in eukaryotic cell evolution.

A role for the Psi-U mismatch in the recognition of the 5' splice site of yeast introns by the U1 small nuclear ribonucleoprotein particle

The U1 small nuclear ribonucleoprotein particle (snRNP)/5' splice site (5'SS) interaction in yeast is essential for the splicing process and depends on the formation of a short RNA duplex between the 5' arm of U1 snRNA and the 1st intronic nucleotides. This RNA/RNA interaction is characterized by the presence of a mismatch that occurs with almost all yeast introns and concerns nucleotides 4 on the pre-mRNA (a U) and 5 on U1 snRNA (a Psi). The latter nucleotide is well conserved from yeast to vertebrates, but its role in yeast and the significance of the associated mismatch in the U1 snRNA/5'SS interaction have never been fully explained. We report here that the presence of this mismatch is a determinant of stability that mainly affects the off rate of the interaction. To our knowledge this is the first report assigning a function to this noncanonical interaction. We also performed SELEX (systematic evolution of ligands by exponential enrichment) experiments by immunoprecipitating U1 snRNP and the associated RNA. The artificial phylogeny derived from these experiments allows the isolation of the selective pressure due to U1 snRNP binding on the 5'SS of yeast introns.

A spliceosomal intron in Giardia lamblia

Short introns occur in numerous protist lineages, but there are no reports of intervening sequences in the protists Giardia lamblia and Trichomonas vaginalis, which may represent the deepest known branches in the eukaryotic line of descent. We have discovered a 35-bp spliceosomal intron in a gene encoding a putative [2Fe-2S] ferredoxin of G. lamblia. The Giardia intron contains a canonical splice site at its 3' end (AG), a noncanonical splice site at its 5' end (CT), and a branch point sequence that fits the yeast consensus sequence of TACTAAC except for the first nucleotide (AACTAAC). We have also identified several G. lamblia genes with spliceosomal peptides, including homologues of eukaryote-specific spliceosomal peptides (Prp8 and Prp11), several DExH-box RNA-helicases that have homologues in eubacteria, but serve essential functions in the splicing of introns in eukaryotes, and 11 predicted archaebacteria-like Sm and like-Sm core peptides, which coat small nuclear RNAs. Phylogenetic analyses show the Giardia Sm core peptides are the products of multiple, ancestral gene duplications followed by divergence, but they retain strong similarity to Sm and like-Sm peptides of other eukaryotes. Although we have documented only a single intron in Giardia, it likely has other introns and fully functional, spliceosomal machinery. If introns were added during eukaryotic evolution (the introns-late hypothesis), then these results push back the date of this event before the branching of G. lamblia.

The new phylogeny of eukaryotes

Molecular phylogeny has been regarded as the ultimate tool for the reconstruction of relationships among eukaryotes-especially the different protist groups-given the difficulty in interpreting morphological data from an evolutionary point of view. In fact, the use of ribosomal RNA as a marker has provided the first well resolved eukaryotic phylogenies, leading to several important evolutionary hypotheses. The most significant is that several early-emerging, amitochondriate lineages, are living relics from the early times of eukaryotic evolution. The use of alternative protein markers and the recognition of several molecular phylogeny reconstruction artefacts, however, have strongly challenged these ideas. The putative early emerging lineages have been demonstrated as late-emerging ones, artefactually misplaced to the base of the tree. The present state of eukaryotic evolution is best described by a multifurcation, in agreement with the 'big bang' hypothesis that assumes a rapid diversification of the major eukaryotic phyla. For further resolution, the analysis of genomic data through improved phylogenetic methods will be required.

Evolutionary appearance of genes encoding proteins associated with box H/ACA snoRNAs: cbf5p in Euglena gracilis, an early diverging eukaryote, and candidate Gar1p and Nop10p homologs in archaebacteria

A reverse transcription-polymerase chain reaction (RT-PCR) approach was used to clone a cDNA encoding the Euglena gracilis homolog of yeast Cbf5p, a protein component of the box H/ACA class of snoRNPs that mediate pseudouridine formation in eukaryotic rRNA. Cbf5p is a putative pseudouridine synthase, and the Euglena homolog is the first full-length Cbf5p sequence to be reported for an early diverging unicellular eukaryote (protist). Phylogenetic analysis of putative pseudouridine synthase sequences confirms that archaebacterial and eukaryotic (including Euglena) Cbf5p proteins are specifically related and are distinct from the TruB/Pus4p clade that is responsible for formation of pseudouridine at position 55 in eubacterial (TruB) and eukaryotic (Pus4p) tRNAs. Using a bioinformatics approach, we also identified archaebacterial genes encoding candidate homologs of yeast Gar1p and Nop10p, two additional proteins known to be associated with eukaryotic box H/ACA snoRNPs. These observations raise the possibility that pseudouridine formation in archaebacterial rRNA may be dependent on analogs of the eukaryotic box H/ACA snoRNPs, whose evolutionary origin may therefore predate the split between Archaea (archaebacteria) and Eucarya (eukaryotes). Database searches further revealed, in archaebacterial and some eukaryotic genomes, two previously unrecognized groups of genes (here designated 'PsuX' and 'PsuY') distantly related to the Cbf5p/TruB gene family.

Membrane heredity and early chloroplast evolution

Membrane heredity was central to the unique symbiogenetic origin from cyanobacteria of chloroplasts in the ancestor of Plantae (green plants, red algae, glaucophytes) and to subsequent lateral transfers of plastids to form even more complex photosynthetic chimeras. Each symbiogenesis integrated disparate genomes and several radically different genetic membranes into a more complex cell. The common ancestor of Plantae evolved transit machinery for plastid protein import. In later secondary symbiogeneses, signal sequences were added to target proteins across host perialgal membranes: independently into green algal plastids (euglenoids, chlorarachneans) and red algal plastids (alveolates, chromists). Conservatism and innovation during early plastid diversification are discussed.

Structural conservation and variation among U5 small nuclear RNAs from trypanosomatid protozoa

U5 snRNAs in trypanosomatid protozoa do not contain the trimethylguanosine cap structures that are often targeted in snRNA isolation procedures. As a result, the trypanosomatids are not well represented in the database of available U5 snRNA sequences. We have isolated and determined the sequence of the U5 snRNA from Crithidia fasciculata. Comparison with previously published trypanosomatid U5 snRNA sequences allows us to deduce the pattern of structural conservation and variation among these very divergent snRNA molecules.

A new twist in trypanosome RNA metabolism: cis-splicing of pre-mRNA

It has been known for almost a decade and a half that in trypanosomes all mRNAs are trans-spliced by addition to the 5' end of the spliced leader (SL) sequence. During the same time period the conviction developed that classical cis-splicing introns are not present in the trypanosome genome and that the trypanosome gene arrangement is highly compact with small intergenic regions separating one gene from the next. We have now discovered that these tenets are no longer true. Poly(A) polymerase (PAP) genes in Trypanosoma brucei and Trypanosoma cruzi are split by intervening sequences of 653 and 302 nt, respectively. The intervening sequences occur at identical positions in both organisms and obey the GT/AG rule of cis-splicing introns. PAP mRNAs are trans-spliced at the very 5' end as well as internally at the 3' splice site of the intervening sequence. Interestingly, 11 nucleotide positions past the actual 5' splice site are conserved between the T. bruceiand T. cruzi introns. Point mutations in these conserved positions, as well as in the AG dinucleotide of the 3' splice site, abolish intron removal in vivo. Our results, together with the recent discovery of cis-splicing introns in Euglena gracilis, suggest that both trans- and cis-splicing are ancient acquisitions of the eukaryotic cell.

Viral and cellular mRNA capping: past and prospects

This chapter focuses on the history of the discovery of cap and an update of research on viral and cellular-messenger RNA (mRNA) capping. Cap structures of the type m⁷ GpppN(m)pN(m)p are present at the 5′ ends of nearly all eukaryotic cellular and viral mRNAs. A cap is added to cellular mRNA precursors and to the transcripts of viruses that replicate in the nucleus during the initial phases of transcription and before other processing events, including internal N⁶A methylation, 3′-poly (A) addition, and exon splicing. Despite the variations on the methylation theme, the important biological consequences of a cap structure appear to correlate with the N⁷-methyl on the 5′-terminal G and the two pyrophosphoryl bonds that connect m⁷G in a 5′–5′ configuration to the first nucleotide of mRNA. In addition to elucidating the biochemical mechanisms of capping and the downstream effects of this 5′- modification on gene expression, the advent of gene cloning has made available an ever-increasing amount of information on the proteins responsible for producing caps and the functional effects of other cap-related interactions. Genetic approaches have demonstrated the lethal consequences of cap failure in yeasts, and complementation studies have shown the evolutionary functional conservation of capping from unicellular to metazoan organisms.

Spliced leader-associated RNA from Crithidia fasciculata contains a structure resembling stem/loop II of U1 snRNA

In contrast to earlier proposals, recent evidence suggests that trans-spliceosomes in trypanosomatid protozoa may contain a homolog of U1 small nuclear (sn) RNA (Schnare, M.N. and Gray, M.W. (1999) J. Biol. Chem. 274, 23,691-23,694). However, the candidate trypanosomatid U1 snRNA is unconventional because it lacks the highly conserved stem/loop II present in all other U1 snRNAs. Trypanosomatids also possess a unique spliced leader-associated (SLA) RNA of unknown function. We present the complete sequence of the SLA RNA from Crithidia fasciculata and propose that it may contribute a U1 snRNA-like stem/loop II to the trans-spliceosome.

A candidate U1 small nuclear RNA for trypanosomatid protozoa

In trypanosomatid protozoa, all mRNAs obtain identical 5'-ends by trans-splicing of the 5'-terminal 39 nucleotides of a small spliced leader RNA to appropriate acceptor sites in pre-mRNA. Although this process involves spliceosomal small nuclear (sn) RNAs, it is thought that trypanosomatids do not contain a homolog of the cis-spliceosomal U1 snRNA. We show here that a trypanosomatid protozoon, Crithidia fasciculata, contains a novel small RNA that displays several features characteristic of a U1 snRNA, including (i) a methylguanosine cap and additional 5'-terminal modifications, (ii) a potential binding site for common core proteins that are present in other trans-spliceosomal ribonucleoproteins, (iii) a U1-like 5'-terminal sequence, and (iv) a U1-like stem/loop I structure. Because trypanosomatid pre-mRNAs do not appear to contain cis-spliced introns, we argue that this previously unrecognized RNA species is a good candidate to be a trans-spliceosomal U1 snRNA.