U2 and U6 snRNA genes in the microsporidian Nosema locustae: evidence for a functional spliceosome

Characterization of Pre-mRNA Splicing and Spliceosomal Machinery in Porphyridium purpureum and Evolutionary Implications for Red Algae

Pre-mRNA splicing is a highly conserved eukaryotic process, but our understanding of it is limited by a historical focus on well-studied organisms such as humans and yeast. There is considerable diversity in mechanisms and components of pre-mRNA splicing, especially in lineages that have evolved under the pressures of genome reduction. The ancestor of red algae is thought to have undergone genome reduction prior to the lineage's radiation, resulting in overall gene and intron loss in extant groups. Previous studies on the extremophilic red alga Cyanidioschyzon merolae revealed an intron-sparse genome with a highly reduced spliceosome. To determine whether these features applied to other red algae, we investigated multiple aspects of pre-mRNA splicing in the mesophilic red alga Porphyridium purpureum. Through strand-specific RNA-Seq, we observed high levels of intron retention across a large number of its introns, and nearly half of the transcripts for these genes are not spliced at all. We also discovered a relationship between variability of 5' splice site sequences and levels of splicing. To further investigate the connections between intron retention and splicing machinery, we bioinformatically assembled the P. purpureum spliceosome, and biochemically verified the presence of snRNAs. While most other core spliceosomal components are present, our results suggest highly divergent or missing U1 snRNP proteins, despite the presence of an uncharacteristically long U1 snRNA. These unusual aspects highlight the diverse nature of pre-mRNA splicing that can be seen in lesser-studied eukaryotes, raising the importance of investigating fundamental eukaryotic processes outside of model organisms.

Numerous fragmented spliceosomal introns, AT-AC splicing, and an unusual dynein gene expression pathway in Giardia lamblia

Spliceosomal introns are hallmarks of eukaryotic genomes, dividing coding regions into separate exons, which are joined during mRNA intron removal catalyzed by the spliceosome. With few known exceptions, spliceosomal introns are cis-spliced, that is, removed from one contiguous pre-mRNA transcript. The protistan intestinal parasite Giardia lamblia exhibits one of the most reduced eukaryotic genomes known, with short intergenic regions and only four known spliceosomal introns. Our genome-wide search for additional introns revealed four unusual cases of spliceosomal intron fragmentation, with consecutive exons of conserved protein-coding genes being dispersed to distant genomic sites. Independent transcripts are trans-spliced to yield contiguous mature mRNAs. Most strikingly, a dynein heavy chain subunit is both interrupted by two fragmented introns and also predicted to be assembled as two separately translated polypeptides, a remarkably complex expression pathway for a nuclear-encoded sequence. For each case, we observe extensive base-pairing potential between intron halves. This base pairing provides both a rationale for the in vivo association of independently transcribed mRNAs transcripts and the apparent specificity of splicing. Similar base-pairing potential in two cis-spliced G. lamblia introns suggests an evolutionary pathway whereby intron fragmentation of cis-spliced introns is permissible and a preliminary evolutionary step to complete gene fission. These results reveal remarkably complex genome dynamics in a severely genomically reduced parasite.

Evolution of small nuclear RNAs in S. cerevisiae, C. albicans, and other hemiascomycetous yeasts

The spliceosome is a large, dynamic ribonuclear protein complex, required for the removal of intron sequences from newly synthesized eukaryotic RNAs. The spliceosome contains five essential small nuclear RNAs (snRNAs): U1, U2, U4, U5, and U6. Phylogenetic comparisons of snRNAs from protists to mammals have long demonstrated remarkable conservation in both primary sequence and secondary structure. In contrast, the snRNAs of the hemiascomycetous yeast Saccharomyces cerevisiae have highly unusual features that set them apart from the snRNAs of other eukaryotes. With an emphasis on the pathogenic yeast Candida albicans, we have now identified and compared snRNAs from newly sequenced yeast genomes, providing a perspective on spliceosome evolution within the hemiascomycetes. In addition to tracing the origins of previously identified snRNA variations present in Saccharomyces cerevisiae, we have found numerous unexpected changes occurring throughout the hemiascomycetous lineages. Our observations reveal interesting examples of RNA and protein coevolution, giving rise to altered interaction domains, losses of deeply conserved snRNA-binding proteins, and unique snRNA sequence changes within the catalytic center of the spliceosome. These same yeast lineages have experienced exceptionally high rates of intron loss, such that modern hemiascomycetous genomes contain introns in only approximately 5% of their genes. Also, the splice site sequences of those introns that remain adhere to an unusually strict consensus. Some of the snRNA variations we observe may thus reflect the altered intron landscape with which the hemiascomycetous spliceosome must contend.

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.

Free energy landscapes of RNA/RNA complexes: with applications to snRNA complexes in spliceosomes

We develop a statistical mechanical model for RNA/RNA complexes with both intramolecular and intermolecular interactions. As an application of the model, we compute the free energy landscapes, which give the full distribution for all the possible conformations, for U4/U6 and U2/U6 in major spliceosome and U4atac/U6atac and U12/U6atac in minor spliceosome. Different snRNA experiments found contrasting structures, our free energy landscape theory shows why these structures emerge and how they compete with each other. For yeast U2/U6, the model predicts that the two distinct experimental structures, the four-helix junction structure and the helix Ib-containing structure, can actually coexist and specifically compete with each other. In addition, the energy landscapes suggest possible mechanisms for the conformational switches in splicing. For instance, our calculation shows that coaxial stacking is essential for stabilizing the four-helix junction in yeast U2/U6. Therefore, inhibition of the coaxial stacking possibly by protein-binding may activate the conformational switch from the four-helix junction to the helix Ib-containing structure. Moreover, the change of the energy landscape shape gives information about the conformational changes. We find multiple (native-like and misfolded) intermediates formed through base-pairing rearrangements in snRNA complexes. For example, the unfolding of the U2/U6 undergoes a transition to a misfolded state which is functional, while in the unfolding of U12/U6atac, the functional helix Ib is found to be the last one to unfold and is thus the most stable structural component. Furthermore, the energy landscape gives the stabilities of all the possible (functional) intermediates and such information is directly related to splicing efficiency.

Molecular phylogeny of the Microsporidia: ecological, ultrastructural and taxonomic considerations

The Microsporidia are a group of obligate intracellular parasites, now thought to be derived fungi. Presented here is a comparative small subunit rDNA (ssrDNA) analysis of 125 species of Microsporidia (sequences obtained from GenBank). This analysis shows that groups or clades are formed based largely on habitat and host. This result is supported by comparative molecular analyses of the past decade, and indicates that structural and ultrastructural characters are unreliable for distinguishing among higher-level microsporidian taxa. Our findings indicate the presence of five major clades of Microsporidia which group according to habitat. We present three new classes of Microsporidia based on natural phylogenetic groupings as illustrated by the ssrDNA analysis: Aquasporidia, Marinosporidia and Terresporidia. The names of the proposed classes reflect the habitat of each group. The class Aquasporidia, found primarily in freshwater habitats, is a paraphyletic group consisting of three clades. The Marinosporidia are found in hosts of marine origin and the Terresporidia are primarily from terrestrial environments.

Rates of intron loss and gain: implications for early eukaryotic evolution

We study the intron-exon structures of 684 groups of orthologs from seven diverse eukaryotic genomes and provide maximum likelihood estimates for rates and numbers of intron losses and gains in these same genes for a variety of lineages. Rates of intron loss vary from approximately 2 x 10(-9) to 2 x 10(-10) per year. Rates of gain vary from 6 x 10(-13) to 4 x 10(-12) per possible intron insertion site per year. There is an inverse correspondence between rates of intron loss and gain, leading to a 20-fold variation among lineages in the ratio of the rates of the two processes. The observed rates of intron gain are insufficient to explain the large number of introns estimated to have been present in the plant-animal ancestor, suggesting that introns present in early eukaryotes may have been created by a fundamentally different process than more recently gained introns.

Complex early genes

We use the pattern of intron conservation in 684 groups of orthologs from seven fully sequenced eukaryotic genomes to provide maximum likelihood estimates of the number of introns present in the same orthologs in various eukaryotic ancestors. We find: (i) intron density in the plant-animal ancestor was high, perhaps two-thirds that of humans and three times that of Drosophila; and (ii) intron density in the ancestral bilateran was also high, equaling that of humans and four times that of Drosophila. We further find that modern introns are generally very old, with two-thirds of modern bilateran introns dating to the ancestral bilateran and two-fifths of modern plant, animal, and fungus introns dating to the plant-animal ancestor. Intron losses outnumber gains over a large range of eukaryotic lineages. These results show that early eukaryotic gene structures were very complex, and that simplification, not embellishment, has dominated subsequent evolution.

Transfer of Nosema locustae (Microsporidia) to Antonospora locustae n. comb. based on molecular and ultrastructural data

Nosema locustae is a microsporidian parasite of grasshopper pests that is used as a biological control agent, and is one of the emerging model systems for microsporidia. Due largely to its diplokaryotic nuclei, N. locustae has been classified in the genus Nosema, a large genus with members that infect a wide variety of insects. However, some molecular studies have cast doubt on the validity of certain Nosema species, and on the taxonomic position of N. locustae. To clarify the affinities of this important insect parasite we sequenced part of the rRNA operon of N. locustae and conducted a phylogenetic analysis using the complete small subunit rRNA gene. Nosema locustae is only distantly related to the nominotypic N. bombycis, and is instead closely related to Antonospora scoticae, a recently described parasite of bees. We examined the ultrastructure of mature N. locustae spores, and found the spore wall to differ from true Nosema species in having a multi-layered exospore resembling that of Antonospora (one of the distinguishing features of that genus). Based on both molecular and morphological evidence, therefore, we propose transferring N. locustae to the genus Antonospora, as Antonospora locustae n. comb.

Mystery of intron gain

For nearly 15 years, it has been widely believed that many introns were recently acquired by the genes of multicellular organisms. However, the mechanism of acquisition has yet to be described for a single animal intron. Here, we report a large-scale computational analysis of the human, Drosophila melanogaster, Caenorhabditis elegans, and Arabidopsis thaliana genomes. We divided 147,796 human intron sequences into batches of similar lengths and aligned them with each other. Different types of homologies between introns were found, but none showed evidence of simple intron transposition. Also, 106,902 plant, 39,624 Drosophila, and 6021 C. elegans introns were examined. No single case of homologous introns in nonhomologous genes was detected. Thus, we found no example of transposition of introns in the last 50 million years in humans, in 3 million years in Drosophila and C. elegans, or in 5 million years in Arabidopsis. Either new introns do not arise via transposition of other introns or intron transposition must have occurred so early in evolution that all traces of homology have been lost.

Microsporidia: biology and evolution of highly reduced intracellular parasites

Microsporidia are a large group of microbial eukaryotes composed exclusively of obligate intracellular parasites of other eukaryotes. Almost 150 years of microsporidian research has led to a basic understanding of many aspects of microsporidian biology, especially their unique and highly specialized mode of infection, where the parasite enters its host through a projectile tube that is expelled at high velocity. Molecular biology and genomic studies on microsporidia have also drawn attention to many other unusual features, including a unique core carbon metabolism and genomes in the size range of bacteria. These seemingly simple parasites were once thought to be the most primitive eukaryotes; however, we now know from molecular phylogeny that they are highly specialized fungi. The fungal nature of microsporidia indicates that microsporidia have undergone severe selective reduction permeating every level of their biology: From cell structures to metabolism, and from genomics to gene structure, microsporidia are reduced.

Large-scale comparison of intron positions among animal, plant, and fungal genes

We purge large databases of animal, plant, and fungal intron-containing genes to a 20% similarity level and then identify the most similar animal-plant, animal-fungal, and plant-fungal protein pairs. We identify the introns in each BLAST 2.0 alignment and score matched intron positions and slid (near-matched, within six nucleotides) intron positions automatically. Overall we find that 10% of the animal introns match plant positions, and a further 7% are "slides." Fifteen percent of fungal introns match animal positions, and 13% match plant positions. Furthermore, the number of alignments with high numbers of matches deviates greatly from the Poisson expectation. The 30 animal-plant alignments with the highest matches (for which 44% of animal introns match plant positions) when aligned with fungal genes are also highly enriched for triple matches: 39% of the fungal introns match both animal and plant positions. This is strong evidence for ancestral introns predating the animal-plant-fungal divergence, and in complete opposition to any expectations based on random insertion. In examining the slid introns, we show that at least half are caused by imperfections in the alignments, and are most likely to be actual matches at common positions. Thus, our final estimates are that approximately equal 14% of animal introns match plant positions, and that approximately equal 17-18% of fungal introns match animal or plant positions, all of these being likely to be ancestral in the eukaryotes.

Multiple forms of U2 snRNA coexist in the silk moth Bombyx mori

Eight U2 snRNA variants were isolated from several Bombyx mori U2-specific RT-PCR libraries. U2 sequences and secondary structures were generated and examined in terms of potential RNA and protein interactions. Analysis indicated that nucleotide changes occurred in both stem/loop and single-stranded areas. Changes in the double stranded areas were either compensatory, single substitutions (e.g. C <--> U) or prevented the double-stranded formation of one or two base pairs. The polymorphisms were clustered in moderately conserved regions. Some of the changes observed generated stronger base pairing. Inter-species conserved protein or RNA-binding sites were relatively unaffected. No polymorphic sites were found in known functional sequences. Bombyx mori and Drosophila melanogaster U2 sequences are 95% and 70% similar at the 5'- and the 3'-ends of the molecule, respectively. Phylogenetic analysis of the U2 sequences demonstrates remarkable conservation across species.

Reconstructing/deconstructing the earliest eukaryotes: how comparative genomics can help

We could reconstruct the evolution of eukaryote-specific molecular and cellular machinery if some living eukaryotes retained primitive cellular structures and we knew which eukaryotes these were. It's not clear that either is the case, but the expanding protist genomic database could help us in several ways.

Alpha and beta subunits of pyruvate dehydrogenase E1 from the microsporidian Nosema locustae: mitochondrion-derived carbon metabolism in microsporidia

Microsporidia are highly adapted eukaryotic intracellular parasites that infect a variety of animals. Microsporidia contain no recognisable mitochondrion, but recently have been shown to have evolved from fungi and to possess heat shock protein genes derived from mitochondria. These findings make it clear that microsporidian ancestors were mitochondrial, yet it remains unknown whether they still contain the organelle, and if so what its role in microsporidian metabolism might be. Here we have characterised genes encoding the alpha and beta subunits of pyruvate dehydrogenase complex E1 (PDH, EC 1.2.4.1) from the microsporidian Nosema locustae. All other amitochondriate eukaryotes studied to date have lost the PDH complex and replaced it with pyruvate:ferredoxin oxidoreductase (PFOR). Nevertheless, molecular phylogeny shows that these Nosema enzymes are most closely related to mitochondrial PDH from other eukaryotes, demonstrating that elements of mitochondrial metabolism have been retained in microsporidia, and that PDH has not been wholly lost. However, there is still no evidence for a mitochondrion in microsporidia, and neither PDH subunit is predicted to encode an amino terminal leader sequence that could function as a mitochondrion-targeting transit peptide, raising questions as to whether these proteins function in a relic organelle or in the cytosol. Moreover, it is also unclear whether these proteins remain part of the PDH complex, or whether they have been retained for another purpose. We propose that microsporidia may utilise a unique pyruvate decarboxylation pathway involving PDH, demonstrating once again the diversity of core metabolism in amitochondriate eukaryotes.

Molecular characteristics and physiology of microsporidia

A survey of the molecular features of microsporidia is presented which attempts to comment on unresolved questions concerning the physiology of these amitochondrial intracellular parasites. Various transports of host-derived molecules can be predicted and trehalose appears as a potential reserve of glucose for energy metabolism. Significant insights into membrane lipids, polyamine metabolism and sporogony-specific proteins have been gained. Some species, such as Encephalitozoon cuniculi, are heterogeneous entities and harbor a small genome. Although showing a variation in genome size of 8.5-fold, microsporidia share reduced rDNA genes. Finally, data on gene organization and a possible evolutionary relationship with fungi are considered.

Sequence and analysis of chromosome I of the amitochondriate intracellular parasite Encephalitozoon cuniculi (Microspora)

A DNA sequencing program was applied to the small (<3 Mb) genome of the microsporidian Encephalitozoon cuniculi, an amitochondriate eukaryotic parasite of mammals, and the sequence of the smallest chromosome was determined. The approximately 224-kb E. cuniculi chromosome I exhibits a dyad symmetry characterized by two identical 37-kb subtelomeric regions which are divergently oriented and extend just downstream of the inverted copies of an 8-kb duplicated cluster of six genes. Each subtelomeric region comprises a single 16S-23S rDNA transcription unit, flanked by various tandemly repeated sequences, and ends with approximately 1 kb of heterogeneous telomeric repeats. The central (or core) region of the chromosome harbors a highly compact arrangement of 132 potential protein-coding genes plus two tRNA genes (one gene per 1.14 kb). Most genes occur as single copies with no identified introns. Of these putative genes, only 53 could be assigned to known functions. A number of genes from the transcription and translation machineries as well as from other cellular processes display characteristic eukaryotic signatures or are clearly eukaryote-specific.

The microsporidian Encephalitozoon

Encephalitozoon cuniculi is an attractive model system for amitochondriate intracellular eukaryotic parasites. It is characterized by a very small genome (below 3 Mbp) and a unique invasion apparatus. Furthermore, as an infectious agent, it is important in human and veterinary medicine. The compactness of its genome involves the reduction of rDNA sequences as well as of some protein-coding genes and intergenic regions. Its highly differentiated apparatus to penetrate the host cell, an extrusome-like polar tube, is composed of novel proteins and may permit various pathways of infestation. Completion of the systematic E. cuniculi sequencing project should provide an important reference system for the comparative genomics of amitochondriate and mitochondriate parasites. Further analysis of orphan genes should help to identify factors that are responsible for its intracellular parasitic way of life.

Coupled nucleotide covariations reveal dynamic RNA interaction patterns

Evolutionarily conserved structures in related RNA molecules contain coordinated variations (covariations) of paired nucleotides. Analysis of covariations is a very powerful approach to deduce phylogenetically conserved (i.e., functional) conformations, including tertiary interactions. Here we discuss conserved RNA folding pathways that are revealed by covariation patterns. In such pathways, structural requirements for alternative pairings cause some nucleotides to covary with two different partners. Such "coupled" covariations between three or more nucleotides were found in various types of RNAs. The analysis of coupled covariations can unravel important features of RNA folding dynamics and improve phylogeny reconstruction in some cases. Importantly, it is necessary to distinguish between multiple covariations determined by mutually exclusive structures and those determined by tertiary contacts.

Towards the minimal eukaryotic parasitic genome

Microsporidia are well-known to infect immunocompromised patients and are also responsible for clinical syndromes in immunocompetent individuals. In recent years, evidence has been obtained in support of a very close relationship between Microsporidia and Fungi. In some species, the compaction of the genome and genes is remarkable. Thus, a systematic sequencing project has been initiated for the 2.9 Mbp genome of Encephalitozoon cuniculi, which will be useful for future comparative genomic studies.

Microsporidia: emerging advances in understanding the basic biology of these unique organisms

Microsporidia are long-known parasites of a wide variety of invertebrate and vertebrate hosts. The emergence of these obligate intracellular organisms as important opportunistic pathogens during the AIDS pandemic and the discovery of new species in humans renewed interest in this unique group of organisms. This review summarises recent advances in the field of molecular biology of microsporidia which (i) contributed to the understanding of the natural origin of human-infecting microsporidia, (ii) revealed unique genetic features of their dramatically reduced genome and (iii) resulted in the correction of their phylogenetic placement among eukaryotes from primitive protozoans to highly evolved organisms related to fungi. Microsporidia might serve as new intracellular model organisms in the future given that gene transfer systems will be developed.

Molecular techniques for detection, species differentiation, and phylogenetic analysis of microsporidia

Microsporidia are obligate intracellular protozoan parasites that infect a broad range of vertebrates and invertebrates. These parasites are now recognized as one of the most common pathogens in human immunodeficiency virus-infected patients. For most patients with infectious diseases, microbiological isolation and identification techniques offer the most rapid and specific determination of the etiologic agent. This is not a suitable procedure for microsporidia, which are obligate intracellular parasites requiring cell culture systems for growth. Therefore, the diagnosis of microsporidiosis currently depends on morphological demonstration of the organisms themselves. Although the diagnosis of microsporidiosis and identification of microsporidia by light microscopy have greatly improved during the last few years, species differentiation by these techniques is usually impossible and transmission electron microscopy may be necessary. Immunfluorescent-staining techniques have been developed for species differentiation of microsporidia, but the antibodies used in these procedures are available only at research laboratories at present. During the last 10 years, the detection of infectious disease agents has begun to include the use of nucleic acid-based technologies. Diagnosis of infection caused by parasitic organisms is the last field of clinical microbiology to incorporate these techniques and molecular techniques (e.g., PCR and hybridization assays) have recently been developed for the detection, species differentiation, and phylogenetic analysis of microsporidia. In this paper we review human microsporidial infections and describe and discuss these newly developed molecular techniques.

Microsporidia are related to Fungi: evidence from the largest subunit of RNA polymerase II and other proteins

We have determined complete gene sequences encoding the largest subunit of the RNA polymerase II (RBP1) from two Microsporidia, Vairimorpha necatrix and Nosema locustae. Phylogenetic analyses of these and other RPB1 sequences strongly support the notion that Microsporidia are not early-diverging eukaryotes but instead are specifically related to Fungi. Our reexamination of elongation factors EF-1alpha and EF-2 sequence data that had previously been taken as support for an early (Archezoan) divergence of these amitochondriate protists show such support to be weak and likely caused by artifacts in phylogenetic analyses. These EF data sets are, in fact, not inconsistent with a Microsporidia + Fungi relationship. In addition, we show that none of these proteins strongly support a deep divergence of Parabasalia and Metamonada, the other amitochondriate protist groups currently thought to compose early branches. Thus, the phylogenetic placement among eukaryotes for these protist taxa is in need of further critical examination.

The recent origins of spliceosomal introns revisited

Does the intron/exon structure of eukaryotic genes belie their ancient assembly by exon-shuffling or have introns been inserted into preformed genes during eukaryotic evolution? These are the central questions in the ongoing 'introns-early' versus 'introns-late' controversy. The phylogenetic distribution of spliceosomal introns continues to strongly favor the intronslate theory. The introns-early theory, however, has claimed support from intron phase and protein structure correlations.