Identification of the major spliceosomal RNAs in Dictyostelium discoideum reveals developmentally regulated U2 variants and polyadenylated snRNAs

Unconventional features in the transcription and processing of spliceosomal small nuclear RNAs in the protozoan parasite Trichomonas vaginalis

Trichomonas vaginalis is a medically important protozoan parasite, and a deep-branching, evolutionarily divergent unicellular eukaryote that has conserved several key features of eukaryotic gene expression. Trichomonas vaginalis possesses a metazoan/plant-like capping apparatus, mRNAs with a cap 1 structure and spliceosomes containing the five small nuclear RNAs (snRNAs). However, in contrast to metazoan and plant snRNAs, the structurally conserved T. vaginalis snRNAs were initially identified as lacking the canonical guanosine cap nucleotide. To explain this unusual condition, we sought to investigate transcriptional and processing features of the spliceosomal snRNAs in this protist. Here, we show that T. vaginalis spliceosomal snRNA genes mostly lack typical eukaryotic promoters. In contrast to other eukaryotes, the putative TATA box in the T. vaginalis U6 snRNA gene was found to be dispensable for transcription or RNA polymerase selectivity. Moreover, U6 transcription in T. vaginalis was virtually insensitive to tagetitoxin compared with other cellular transcripts produced by the same RNA polymerase III. Most important and unexpected, snRNA transcription in T. vaginalis appears to bypass capping as we show that these transcripts retain their original 5'-triphosphate groups. In conclusion, transcription and processing of spliceosomal snRNAs in T. vaginalis deviate considerably from the conventional rules of other eukaryotes.

Abundantly expressed class of noncoding RNAs conserved through the multicellular evolution of dictyostelid social amoebas

Aggregative multicellularity has evolved multiple times in diverse groups of eukaryotes, exemplified by the well-studied development of dictyostelid social amoebas, for example, Dictyostelium discoideum However, it is still poorly understood why multicellularity emerged in these amoebas while the majority of other members of Amoebozoa are unicellular. Previously, a novel type of noncoding RNA, Class I RNAs, was identified in D. discoideum and shown to be important for normal multicellular development. Here, we investigated Class I RNA evolution and its connection to multicellular development. We identified a large number of new Class I RNA genes by constructing a covariance model combined with a scoring system based on conserved upstream sequences. Multiple genes were predicted in representatives of each major group of Dictyostelia and expression analysis confirmed that our search approach identifies expressed Class I RNA genes with high accuracy and sensitivity and that the RNAs are developmentally regulated. Further studies showed that Class I RNAs are ubiquitous in Dictyostelia and share highly conserved structure and sequence motifs. In addition, Class I RNA genes appear to be unique to dictyostelid social amoebas because they could not be identified in outgroup genomes, including their closest known relatives. Our results show that Class I RNA is an ancient class of ncRNAs, likely to have been present in the last common ancestor of Dictyostelia dating back at least 600 million years. Based on previous functional analyses and the presented evolutionary investigation, we hypothesize that Class I RNAs were involved in evolution of multicellularity in Dictyostelia.

Rapid protein evolution, organellar reductions, and invasive intronic elements in the marine aerobic parasite dinoflagellate Amoebophrya spp

BACKGROUND

Dinoflagellates are aquatic protists particularly widespread in the oceans worldwide. Some are responsible for toxic blooms while others live in symbiotic relationships, either as mutualistic symbionts in corals or as parasites infecting other protists and animals. Dinoflagellates harbor atypically large genomes (~ 3 to 250 Gb), with gene organization and gene expression patterns very different from closely related apicomplexan parasites. Here we sequenced and analyzed the genomes of two early-diverging and co-occurring parasitic dinoflagellate Amoebophrya strains, to shed light on the emergence of such atypical genomic features, dinoflagellate evolution, and host specialization.

RESULTS

We sequenced, assembled, and annotated high-quality genomes for two Amoebophrya strains (A25 and A120), using a combination of Illumina paired-end short-read and Oxford Nanopore Technology (ONT) MinION long-read sequencing approaches. We found a small number of transposable elements, along with short introns and intergenic regions, and a limited number of gene families, together contribute to the compactness of the Amoebophrya genomes, a feature potentially linked with parasitism. While the majority of Amoebophrya proteins (63.7% of A25 and 59.3% of A120) had no functional assignment, we found many orthologs shared with Dinophyceae. Our analyses revealed a strong tendency for genes encoded by unidirectional clusters and high levels of synteny conservation between the two genomes despite low interspecific protein sequence similarity, suggesting rapid protein evolution. Most strikingly, we identified a large portion of non-canonical introns, including repeated introns, displaying a broad variability of associated splicing motifs never observed among eukaryotes. Those introner elements appear to have the capacity to spread over their respective genomes in a manner similar to transposable elements. Finally, we confirmed the reduction of organelles observed in Amoebophrya spp., i.e., loss of the plastid, potential loss of a mitochondrial genome and functions.

CONCLUSION

These results expand the range of atypical genome features found in basal dinoflagellates and raise questions regarding speciation and the evolutionary mechanisms at play while parastitism was selected for in this particular unicellular lineage.

Nucleomorph Small RNAs in Cryptophyte and Chlorarachniophyte Algae

The regulation of gene expression and RNA maturation underlies fundamental processes such as cell homeostasis, development, and stress acclimation. The biogenesis and modification of RNA is tightly controlled by an array of regulatory RNAs and nucleic acid-binding proteins. While the role of small RNAs (sRNAs) in gene expression has been studied in-depth in select model organisms, little is known about sRNA biology across the eukaryotic tree of life. We used deep sequencing to explore the repertoires of sRNAs encoded by the miniaturized, endosymbiotically derived “nucleomorph” genomes of two single-celled algae, the cryptophyte Guillardia theta and the chlorarachniophyte Bigelowiella natans. A total of 32.3 and 35.3 million reads were generated from G. theta and B. natans, respectively. In G. theta, we identified nucleomorph U1, U2, and U4 spliceosomal small nuclear RNAs (snRNAs) as well as 11 C/D box small nucleolar RNAs (snoRNAs), five of which have potential plant and animal homologs. The snoRNAs are predicted to perform 2′-O methylation of rRNA (but not snRNA). In B. natans, we found the previously undetected 5S rRNA as well as six orphan sRNAs. Analysis of chlorarachniophyte snRNAs shed light on the removal of the miniature 18–21 nt introns found in B. natans nucleomorph genes. Neither of the nucleomorph genomes appears to encode RNA pseudouridylation machinery, and U5 snRNA cannot be found in the cryptophyte G. theta. Considering the central roles of U5 snRNA and RNA modifications in other organisms, cytoplasm-to-nucleomorph RNA shuttling in cryptophyte algae is a distinct possibility.

Spliceosomal introns in Trichomonas vaginalis revisited

Background

The human protozoan parasite Trichomonas vaginalis is an organism of interest for understanding eukaryotic evolution. Despite having an unusually large genome and a rich gene repertoire among protists, spliceosomal introns in T. vaginalis appear rare: only 62 putative introns have been annotated in this genome, and little or no experimental evidence exists to back up these predictions.

Results

This study revisited the 62 annotated introns of T. vaginalis derived from the genome sequencing plus previous publications. After experimental validation and a new genome-wide search, we confirmed the presence of introns in 32 genes and 18 others were concluded to be intronless. Sequence analyses classified the validated introns into two types, based on distinctive features such as length and conservation of splice site motifs.

Conclusions

Our study provides an updated list of intron-containing genes in the genome of T. vaginalis. Our findings suggests the existence of two intron ‘families’ spread among T. vaginalis protein-coding genes. Additional studies are needed to understand the functional separation of these two classes of introns and to assess the existence of further introns in the T. vaginalis genome.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-3196-7) contains supplementary material, which is available to authorized users.

Global characterization of the Dicer-like protein DrnB roles in miRNA biogenesis in the social amoeba Dictyostelium discoideum

Micro (mi)RNAs regulate gene expression in many eukaryotic organisms where they control diverse biological processes. Their biogenesis, from primary transcripts to mature miRNAs, have been extensively characterized in animals and plants, showing distinct differences between these phylogenetically distant groups of organisms. However, comparably little is known about miRNA biogenesis in organisms whose evolutionary position is placed in between plants and animals and/or in unicellular organisms. Here, we investigate miRNA maturation in the unicellular amoeba Dictyostelium discoideum, belonging to Amoebozoa, which branched out after plants but before animals. High-throughput sequencing of small RNAs and poly(A)-selected RNAs demonstrated that the Dicer-like protein DrnB is required, and essentially specific, for global miRNA maturation in D. discoideum. Our RNA-seq data also showed that longer miRNA transcripts, generally preceded by a T-rich putative promoter motif, accumulate in a drnB knock-out strain. For two model miRNAs we defined the transcriptional start sites (TSSs) of primary (pri)-miRNAs and showed that they carry the RNA polymerase II specific m⁷G-cap. The generation of the 3ʹ-ends of these pri-miRNAs differs, with pri-mir-1177 reading into the downstream gene, and pri-mir-1176 displaying a distinct end. This 3´-end is processed to shorter intermediates, stabilized in DrnB-depleted cells, of which some carry a short oligo(A)-tail. Furthermore, we identified 10 new miRNAs, all DrnB dependent and developmentally regulated. Thus, the miRNA machinery in D. discoideum shares features with both plants and animals, which is in agreement with its evolutionary position and perhaps also an adaptation to its complex lifestyle: unicellular growth and multicellular development.

CRISPR/Cas9 mediated targeting of multiple genes in Dictyostelium

CRISPR/Cas9 has emerged in various organisms as a powerful technology for targeted gene knockout; however, no reports of editing the Dictyostelium genome efficiently using this system are available. We describe here the application of CRISPR/Cas9-mediated gene modification in Dictyostelium. The endogenous tRNA-processing system for expressing sgRNA was approximately 10 times more effective than the commonly used U6 promoter. The resulting sgRNA affected the sub-nuclear localisation of Cas9, indicating that the expression level of sgRNA was sufficiently high to form Cas9 and sgRNA complexes within the nucleus. The all-in-one vector containing Cas9 and sgRNA was transiently expressed to generate mutants in five PI3K genes. Mutation detective PCR revealed the mutagenesis frequency of the individual genes to be between 72.9% and 100%. We confirmed that all five targeting loci in the four independent clones had insertion/deletion mutations in their target sites. Thus, we show that the CRISPR/Cas9 system can be used in Dictyostelium cells to enable efficient genome editing of multiple genes. Since this system utilises transient expression of the all-in-one vector, it has the advantage that the drug resistance cassette is not integrated into the genome and simple vector construction, involving annealing two oligo-DNAs.

The Long Noncoding RNA Transcriptome of Dictyostelium discoideum Development

Dictyostelium discoideum live in the soil as single cells, engulfing bacteria and growing vegetatively. Upon starvation, tens of thousands of amoebae enter a developmental program that includes aggregation, multicellular differentiation, and sporulation. Major shifts across the protein-coding transcriptome accompany these developmental changes. However, no study has presented a global survey of long noncoding RNAs (ncRNAs) in D. discoideum To characterize the antisense and long intergenic noncoding RNA (lncRNA) transcriptome, we analyzed previously published developmental time course samples using an RNA-sequencing (RNA-seq) library preparation method that selectively depletes ribosomal RNAs (rRNAs). We detected the accumulation of transcripts for 9833 protein-coding messenger RNAs (mRNAs), 621 lncRNAs, and 162 putative antisense RNAs (asRNAs). The noncoding RNAs were interspersed throughout the genome, and were distinct in expression level, length, and nucleotide composition. The noncoding transcriptome displayed a temporal profile similar to the coding transcriptome, with stages of gradual change interspersed with larger leaps. The transcription profiles of some noncoding RNAs were strongly correlated with known differentially expressed coding RNAs, hinting at a functional role for these molecules during development. Examining the mitochondrial transcriptome, we modeled two novel antisense transcripts. We applied yet another ribosomal depletion method to a subset of the samples to better retain transfer RNA (tRNA) transcripts. We observed polymorphisms in tRNA anticodons that suggested a post-transcriptional means by which D. discoideum compensates for codons missing in the genomic complement of tRNAs. We concluded that the prevalence and characteristics of long ncRNAs indicate that these molecules are relevant to the progression of molecular and cellular phenotypes during development.

The C-Terminal SynMuv/DdDUF926 Domain Regulates the Function of the N-Terminal Domain of DdNKAP

NKAP (NF-κB activating protein) is a highly conserved SR (serine/arginine-rich) protein involved in transcriptional control and splicing in mammals. We identified DdNKAP, the Dictyostelium discoideum ortholog of mammalian NKAP, as interacting partner of the nuclear envelope protein SUN-1. DdNKAP harbors a number of basic RDR/RDRS repeats in its N-terminal domain and the SynMuv/DUF926 domain at its C-terminus. We describe a novel and direct interaction between DdNKAP and Prp19 (Pre mRNA processing factor 19) which might be relevant for the observed DdNKAP ubiquitination. Genome wide analysis using cross-linking immunoprecipitation-high-throughput sequencing (CLIP-seq) revealed DdNKAP association with intergenic regions, exons, introns and non-coding RNAs. Ectopic expression of DdNKAP and its domains affects several developmental aspects like stream formation, aggregation, and chemotaxis. We conclude that DdNKAP is a multifunctional protein, which might influence Dictyostelium development through its interaction with RNA and RNA binding proteins. Mutants overexpressing full length DdNKAP and the N-terminal domain alone (DdN-NKAP) showed opposite phenotypes in development and opposite expression profiles of several genes and rRNAs. The observed interaction between DdN-NKAP and the DdDUF926 domain indicates that the DdDUF926 domain acts as negative regulator of the N-terminus.

Variant U1 snRNAs are implicated in human pluripotent stem cell maintenance and neuromuscular disease

The U1 small nuclear (sn)RNA (U1) is a multifunctional ncRNA, known for its pivotal role in pre-mRNA splicing and regulation of RNA 3' end processing events. We recently demonstrated that a new class of human U1-like snRNAs, the variant (v)U1 snRNAs (vU1s), also participate in pre-mRNA processing events. In this study, we show that several human vU1 genes are specifically upregulated in stem cells and participate in the regulation of cell fate decisions. Significantly, ectopic expression of vU1 genes in human skin fibroblasts leads to increases in levels of key pluripotent stem cell mRNA markers, including NANOG and SOX2. These results reveal an important role for vU1s in the control of key regulatory networks orchestrating the transitions between stem cell maintenance and differentiation. Moreover, vU1 expression varies inversely with U1 expression during differentiation and cell re-programming and this pattern of expression is specifically de-regulated in iPSC-derived motor neurons from Spinal Muscular Atrophy (SMA) type 1 patient's. Accordingly, we suggest that an imbalance in the vU1/U1 ratio, rather than an overall reduction in Uridyl-rich (U)-snRNAs, may contribute to the specific neuromuscular disease phenotype associated with SMA.

Developmental analysis of spliceosomal snRNA isoform expression

Pre-mRNA splicing is a critical step in eukaryotic gene expression that contributes to proteomic, cellular, and developmental complexity. Small nuclear (sn)RNAs are core spliceosomal components; however, the extent to which differential expression of snRNA isoforms regulates splicing is completely unknown. This is partly due to difficulties in the accurate analysis of the spatial and temporal expression patterns of snRNAs. Here, we use high-throughput RNA-sequencing (RNA-seq) data to profile expression of four major snRNAs throughout Drosophila development. This analysis shows that individual isoforms of each snRNA have distinct expression patterns in the embryo, larva, and pharate adult stages. Expression of these isoforms is more heterogeneous during embryogenesis; as development progresses, a single isoform from each snRNA subtype gradually dominates expression. Despite the lack of stable snRNA orthologous groups during evolution, this developmental switching of snRNA isoforms also occurs in distantly related vertebrate species, such as Xenopus, mouse, and human. Our results indicate that expression of snRNA isoforms is regulated and lays the foundation for functional studies of individual snRNA isoforms.

U1 small nuclear RNA variants differentially form ribonucleoprotein particles in vitro

The U1 small nuclear (sn)RNA participates in splicing of pre-mRNAs by recognizing and binding to 5' splice sites at exon/intron boundaries. U1 snRNAs associate with 5' splice sites in the form of ribonucleoprotein particles (snRNPs) that are comprised of the U1 snRNA and 10 core components, including U1A, U1-70K, U1C and the 'Smith antigen', or Sm, heptamer. The U1 snRNA is highly conserved across a wide range of taxa; however, a number of reports have identified the presence of expressed U1-like snRNAs in multiple species, including humans. While numerous U1-like molecules have been shown to be expressed, it is unclear whether these variant snRNAs have the capacity to form snRNPs and participate in splicing. The purpose of the present study was to further characterize biochemically the ability of previously identified human U1-like variants to form snRNPs and bind to U1 snRNP proteins. A bioinformatics analysis provided support for the existence of multiple expressed variants. In vitro gel shift assays, competition assays, and immunoprecipitations (IPs) revealed that the variants formed high molecular weight assemblies to varying degrees and associated with core U1 snRNP proteins to a lesser extent than the canonical U1 snRNA. Together, these data suggest that the human U1 snRNA variants analyzed here are unable to efficiently bind U1 snRNP proteins. The current work provides additional biochemical insights into the ability of the variants to assemble into snRNPs.

A rare combination of ribonucleotide reductases in the social amoeba Dictyostelium discoideum

Ribonucleotide reductases (RNRs) catalyze the only pathway for de novo synthesis of deoxyribonucleotides needed for DNA replication and repair. The vast majority of eukaryotes encodes only a class I RNR, but interestingly some eukaryotes, including the social amoeba Dictyostelium discoideum, encode both a class I and a class II RNR. The amino acid sequence of the D. discoideum class I RNR is similar to other eukaryotic RNRs, whereas that of its class II RNR is most similar to the monomeric class II RNRs found in Lactobacillus spp. and a few other bacteria. Here we report the first study of RNRs in a eukaryotic organism that encodes class I and class II RNRs. Both classes of RNR genes were expressed in D. discoideum cells, although the class I transcripts were more abundant and strongly enriched during mid-development compared with the class II transcript. The quaternary structure, allosteric regulation, and properties of the diiron-oxo/radical cofactor of D. discoideum class I RNR are similar to those of the mammalian RNRs. Inhibition of D. discoideum class I RNR by hydroxyurea resulted in a 90% reduction in spore formation and decreased the germination viability of the surviving spores by 75%. Class II RNR could not compensate for class I inhibition during development, and an excess of vitamin B₁₂ coenzyme, which is essential for class II activity, did not improve spore formation. We suggest that class I is the principal RNR during D. discoideum development and growth and is important for spore formation, possibly by providing dNTPs for mitochondrial replication.

Identification and verification of microRNAs by high-throughput sequencing

High-throughput sequencing methods have become invaluable for detection and analysis of small RNAs. The results are millions of sequences that need to be carefully analyzed by computational methods and preferentially verified by different experimental techniques. Here we describe how to use high-throughput sequencing followed by bioinformatics and northern blot to identify one particular class of small RNA, microRNAs.

A phylogenetic study of Drosophila splicing assembly chaperone RNP-4F associated U4-/U6-snRNA secondary structure

The rnp-4f gene in Drosophila melanogaster encodes nuclear protein RNP-4F. This encoded protein is represented by homologs in other eukaryotic species, where it has been shown to function as an intron splicing assembly factor. Here, RNP-4F is believed to initially bind to a recognition sequence on U6-snRNA, serving as a chaperone to facilitate its association with U4-snRNA by intermolecular hydrogen bonding. RNA conformations are a key factor in spliceosome function, so that elucidation of changing secondary structures for interacting snRNAs is a subject of considerable interest and importance. Among the five snRNAs which participate in removal of spliceosomal introns, there is a growing consensus that U6-snRNA is the most structurally dynamic and may constitute the catalytic core. Previous studies by others have generated potential secondary structures for free U4- and U6-snRNAs, including the Y-shaped U4-/U6-snRNA model. These models were based on study of RNAs from relatively few species, and the popular Y-shaped model remains to be systematically re-examined with reference to the many new sequences generated by recent genomic sequencing projects. We have utilized a comparative phylogenetic approach on 60 diverse eukaryotic species, which resulted in a revised and improved U4-/U6-snRNA secondary structure. This general model is supported by observation of abundant compensatory base mutations in every stem, and incorporates more of the nucleotides into base-paired associations than in previous models, thus being more energetically stable. We have extensively sampled the eukaryotic phylogenetic tree to its deepest roots, but did not find genes potentially encoding either U4- or U6-snRNA in the Giardia and Trichomonas data-bases. Our results support the hypothesis that nuclear introns in these most deeply rooted eukaryotes may represent evolutionary intermediates, sharing characteristics of both group II and spliceosomal introns. An unexpected result of this study was discovery of a potential competitive binding site for Drosophila splicing assembly factor RNP-4F to a 5’-UTR regulatory region within its own premRNA, which may play a role in negative feedback control.

MicroRNAs in Amoebozoa: deep sequencing of the small RNA population in the social amoeba Dictyostelium discoideum reveals developmentally regulated microRNAs

The RNA interference machinery has served as a guardian of eukaryotic genomes since the divergence from prokaryotes. Although the basic components have a shared origin, silencing pathways directed by small RNAs have evolved in diverse directions in different eukaryotic lineages. Micro (mi)RNAs regulate protein-coding genes and play vital roles in plants and animals, but less is known about their functions in other organisms. Here, we report, for the first time, deep sequencing of small RNAs from the social amoeba Dictyostelium discoideum. RNA from growing single-cell amoebae as well as from two multicellular developmental stages was sequenced. Computational analyses combined with experimental data reveal the expression of miRNAs, several of them exhibiting distinct expression patterns during development. To our knowledge, this is the first report of miRNAs in the Amoebozoa supergroup. We also show that overexpressed miRNA precursors generate miRNAs and, in most cases, miRNA* sequences, whose biogenesis is dependent on the Dicer-like protein DrnB, further supporting the presence of miRNAs in D. discoideum. In addition, we find miRNAs processed from hairpin structures originating from an intron as well as from a class of repetitive elements. We believe that these repetitive elements are sources for newly evolved miRNAs.

Spliceosomal genes in the D. discoideum genome: a comparison with those in H. sapiens, D. melanogaster, A. thaliana and S. cerevisiae

Little is known about pre-mRNA splicing in Dictyostelium discoideum although its genome has been completely sequenced. Our analysis suggests that pre-mRNA splicing plays an important role in D. discoideum gene expression as two thirds of its genes contain at least one intron. Ongoing curation of the genome to date has revealed 40 genes in D. discoideum with clear evidence of alternative splicing, supporting the existence of alternative splicing in this unicellular organism. We identified 160 candidate U2-type spliceosomal proteins and related factors in D. discoideum based on 264 known human genes involved in splicing. Spliceosomal small ribonucleoproteins (snRNPs), PRP19 complex proteins and late-acting proteins are highly conserved in D. discoideum and throughout the metazoa. In non-snRNP and hnRNP families, D. discoideum orthologs are closer to those in A. thaliana, D. melanogaster and H. sapiens than to their counterparts in S. cerevisiae. Several splicing regulators, including SR proteins and CUG-binding proteins, were found in D. discoideum, but not in yeast. Our comprehensive catalog of spliceosomal proteins provides useful information for future studies of splicing in D. discoideum where the efficient genetic and biochemical manipulation will also further our general understanding of pre-mRNA splicing.

Sequence and generation of mature ribosomal RNA transcripts in Dictyostelium discoideum

The amoeba Dictyostelium discoideum is a well established model organism for studying numerous aspects of cellular and developmental functions. Its ribosomal RNA (rRNA) is encoded in an extrachromosomal palindrome that exists in ∼100 copies in the cell. In this study, we have set out to investigate the sequence of the expressed rRNA. For this, we have ligated the rRNA ends and performed RT-PCR on these circular RNAs. Sequencing revealed that the mature 26 S, 17 S, 5.8 S, and 5 S rRNAs have sizes of 3741, 1871, 162, and 112 nucleotides, respectively. Unlike the published data, all mature rRNAs of the same type uniformly display the same start and end nucleotides in the analyzed AX2 strain. We show the existence of a short lived primary transcript covering the rRNA transcription unit of 17 S, 5.8 S, and 26 S rRNA. Northern blots and RT-PCR reveal that from this primary transcript two precursor molecules of the 17 S and two precursors of the 26 S rRNA are generated. We have also determined the sequences of these precursor molecules, and based on these data, we propose a model for the maturation of the rRNAs in Dictyostelium discoideum that we compare with the processing of the rRNA transcription unit of Saccharomyces cerevisiae.

Comparative genomics of the social amoebae Dictyostelium discoideum and Dictyostelium purpureum

Background

The social amoebae (Dictyostelia) are a diverse group of Amoebozoa that achieve multicellularity by aggregation and undergo morphogenesis into fruiting bodies with terminally differentiated spores and stalk cells. There are four groups of dictyostelids, with the most derived being a group that contains the model species Dictyostelium discoideum.

Results

We have produced a draft genome sequence of another group dictyostelid, Dictyostelium purpureum, and compare it to the D. discoideum genome. The assembly (8.41 × coverage) comprises 799 scaffolds totaling 33.0 Mb, comparable to the D. discoideum genome size. Sequence comparisons suggest that these two dictyostelids shared a common ancestor approximately 400 million years ago. In spite of this divergence, most orthologs reside in small clusters of conserved synteny. Comparative analyses revealed a core set of orthologous genes that illuminate dictyostelid physiology, as well as differences in gene family content. Interesting patterns of gene conservation and divergence are also evident, suggesting function differences; some protein families, such as the histidine kinases, have undergone little functional change, whereas others, such as the polyketide synthases, have undergone extensive diversification. The abundant amino acid homopolymers encoded in both genomes are generally not found in homologous positions within proteins, so they are unlikely to derive from ancestral DNA triplet repeats. Genes involved in the social stage evolved more rapidly than others, consistent with either relaxed selection or accelerated evolution due to social conflict.

Conclusions

The findings from this new genome sequence and comparative analysis shed light on the biology and evolution of the Dictyostelia.

Plasmodium falciparum spliceosomal RNAs: 3' and 5' end processing

The major spliceosomal small nuclear ribonucleoproteins (snRNPs) consist of snRNA (U1, U2, U4 or U5) and several proteins which can be unique or common to each snRNP particle. The common proteins are known as Sm proteins; they are crucial for RNP assembly and nuclear import of spliceosomal RNPs. This paper reports detecting the interaction between Plasmodium falciparum snRNAs and Sm proteins, and the usual 5' trimethylated caps on the snRNAs, by immunoprecipitation with specific antibodies. Furthermore, an unusual poly(A) tail was detected on these non-coding RNAs.

Evolution of spliceosomal snRNA genes in metazoan animals

While studies of the evolutionary histories of protein families are commonplace, little is known on noncoding RNAs beyond microRNAs and some snoRNAs. Here we investigate in detail the evolutionary history of the nine spliceosomal snRNA families (U1, U2, U4, U5, U6, U11, U12, U4atac, and U6atac) across the completely or partially sequenced genomes of metazoan animals. Representatives of the five major spliceosomal snRNAs were found in all genomes. None of the minor splicesomal snRNAs were detected in nematodes or in the shotgun traces of Oikopleura dioica, while in all other animal genomes at most one of them is missing. Although snRNAs are present in multiple copies in most genomes, distinguishable paralogue groups are not stable over long evolutionary times, although they appear independently in several clades. In general, animal snRNA secondary structures are highly conserved, albeit, in particular, U11 and U12 in insects exhibit dramatic variations. An analysis of genomic context of snRNAs reveals that they behave like mobile elements, exhibiting very little syntenic conservation.

Computational identification of four spliceosomal snRNAs from the deep-branching eukaryote Giardia intestinalis

RNAs processing other RNAs is very general in eukaryotes, but is not clear to what extent it is ancestral to eukaryotes. Here we focus on pre-mRNA splicing, one of the most important RNA-processing mechanisms in eukaryotes. In most eukaryotes splicing is predominantly catalysed by the major spliceosome complex, which consists of five uridine-rich small nuclear RNAs (U-snRNAs) and over 200 proteins in humans. Three major spliceosomal introns have been found experimentally in Giardia; one Giardia U-snRNA (U5) and a number of spliceosomal proteins have also been identified. However, because of the low sequence similarity between the Giardia ncRNAs and those of other eukaryotes, the other U-snRNAs of Giardia had not been found. Using two computational methods, candidates for Giardia U1, U2, U4 and U6 snRNAs were identified in this study and shown by RT-PCR to be expressed. We found that identifying a U2 candidate helped identify U6 and U4 based on interactions between them. Secondary structural modelling of the Giardia U-snRNA candidates revealed typical features of eukaryotic U-snRNAs. We demonstrate a successful approach to combine computational and experimental methods to identify expected ncRNAs in a highly divergent protist genome. Our findings reinforce the conclusion that spliceosomal small-nuclear RNAs existed in the last common ancestor of eukaryotes.

Spliceosomal snRNAs in the unicellular eukaryote Trichomonas vaginalis are structurally conserved but lack a 5'-cap structure

Few genes in the divergent eukaryote Trichomonas vaginalis have introns, despite the unusually large gene repertoire of this human-infective parasite. These introns are characterized by extended conserved regulatory motifs at the 5' and 3' boundaries, a feature shared with another divergent eukaryote, Giardia lamblia, but not with metazoan introns. This unusual characteristic of T. vaginalis introns led us to examine spliceosomal small nuclear RNAs (snRNAs) predicted to mediate splicing reactions via interaction with intron motifs. Here we identify T. vaginalis U1, U2, U4, U5, and U6 snRNAs, present predictions of their secondary structures, and provide evidence for interaction between the U2/U6 snRNA complex and a T. vaginalis intron. Structural models predict that T. vaginalis snRNAs contain conserved sequences and motifs similar to those found in other examined eukaryotes. These data indicate that mechanisms of intron recognition as well as coordination of the two catalytic steps of splicing have been conserved throughout eukaryotic evolution. Unexpectedly, we found that T. vaginalis spliceosomal snRNAs lack the 5' trimethylguanosine cap typical of snRNAs and appear to possess unmodified 5' ends. Despite the lack of a cap structure, U1, U2, U4, and U5 genes are transcribed by RNA polymerase II, whereas the U6 gene is transcribed by RNA polymerase III.

De novo search for non-coding RNA genes in the AT-rich genome of Dictyostelium discoideum: performance of Markov-dependent genome feature scoring

Genome data are increasingly important in the computational identification of novel regulatory non-coding RNAs (ncRNAs). However, most ncRNA gene-finders are either specialized to well-characterized ncRNA gene families or require comparisons of closely related genomes. We developed a method for de novo screening for ncRNA genes with a nucleotide composition that stands out against the background genome based on a partial sum process. We compared the performance when assuming independent and first-order Markov-dependent nucleotides, respectively, and used Karlin-Altschul and Karlin-Dembo statistics to evaluate the significance of hits. We hypothesized that a first-order Markov-dependent process might have better power to detect ncRNA genes since nearest-neighbor models have been shown to be successful in predicting RNA structures. A model based on a first-order partial sum process (analyzing overlapping dinucleotides) had better sensitivity and specificity than a zeroth-order model when applied to the AT-rich genome of the amoeba Dictyostelium discoideum. In this genome, we detected 94% of previously known ncRNA genes (at this sensitivity, the false positive rate was estimated to be 25% in a simulated background). The predictions were further refined by clustering candidate genes according to sequence similarity and/or searching for an ncRNA-associated upstream element. We experimentally verified six out of 10 tested ncRNA gene predictions. We conclude that higher-order models, in combination with other information, are useful for identification of novel ncRNA gene families in single-genome analysis of D. discoideum. Our generalizable approach extends the range of genomic data that can be searched for novel ncRNA genes using well-grounded statistical methods.

Systematic identification of non-coding RNA 2,2,7-trimethylguanosine cap structures in Caenorhabditis elegans

Background

The 2,2,7-trimethylguanosine (TMG) cap structure is an important functional characteristic of ncRNAs with critical cellular roles, such as some snRNAs. Here we used immunoprecipitation with both K121 and R1131 anti-TMG antibodies to systematically identify the TMG cap structures for all presently characterized ncRNAs in C. elegans.

Results

The two anti-TMG antibodies precipitated a similar group of the C. elegans ncRNAs. All snRNAs known to have a TMG cap structure were found in the precipitate, indicating that our identification system was efficient. Other ncRNA families related to splicing, such as SL RNAs and Sm Y RNAs, were also found in the precipitate, as were 7 C/D box snoRNAs. Further analysis showed that the SL RNAs and the Sm Y RNAs shared a very similar Sm binding site element (AAU_4–5GGA), which sequence composition differed somewhat from those of other U snRNAs. There were also 16 ncRNAs without an Sm binding site element in the precipitate, suggesting that for these ncRNAs, TMG formation may occur independently of Sm proteins.

Conclusion

Our results showed that most ncRNAs predicted to be transcribed by RNA polymerase II had a TMG cap, while those predicted to be transcribed by RNA plymerase III or located in introns did not have a TMG cap structure. Compared to ncRNAs without a TMG cap, TMG-capped ncRNAs tended to have higher expression levels. Five functionally non-annotated ncRNAs also have a TMG cap structure, which might be helpful for identifying the cellular roles of these ncRNAs.

Detection of RNA structures in porcine EST data and related mammals

Background

Non-coding RNAs (ncRNAs) are involved in a wide spectrum of regulatory functions. Within recent years, there have been increasing reports of observed polyadenylated ncRNAs and mRNA like ncRNAs in eukaryotes. To investigate this further, we examined the large data set in the Sino-Danish PigEST resource which also contains expression information distributed on 97 non-normalized cDNA libraries.

Results

We constructed a pipeline, EST2ncRNA, to search for known and novel ncRNAs. The pipeline utilises sequence similarity to ncRNA databases (blast), structure similarity to Rfam (RaveNnA) as well as multiple alignments to predict conserved novel putative RNA structures (RNAz). EST2ncRNA was fed with 48,000 contigs and 73,000 singletons available from the PigEST resource. Using the pipeline we identified known RNA structures in 137 contigs and single reads (conreads), and predicted high confidence RNA structures in non-protein coding regions of additional 1,262 conreads. Of these, structures in 270 conreads overlap with existing predictions in human. To sum up, the PigEST resource comprises trans-acting elements (ncRNAs) in 715 contigs and 340 singletons as well as cis-acting elements (inside UTRs) in 311 contigs and 51 singletons, of which 18 conreads contain both predictions of trans- and cis-acting elements. The predicted RNAz candidates were compared with the PigEST expression information and we identify 114 contigs with an RNAz prediction and expression in at least ten of the non-normalised cDNA libraries. We conclude that the contigs with RNAz and known predictions are in general expressed at a much lower level than protein coding transcripts. In addition, we also observe that our ncRNA candidates constitute about one to two percent of the genes expressed in the cDNA libraries. Intriguingly, the cDNA libraries from developmental (brain) tissues contain the highest amount of ncRNA candidates, about two percent. These observations are related to existing knowledge and hypotheses about the role of ncRNAs in higher organisms. Furthermore, about 80% porcine coding transcripts (of 18,600 identified) as well as less than one-third ORF-free transcripts are conserved at least in the closely related bovine genome. Approximately one percent of the coding and 10% of the remaining matches are unique between the PigEST data and cow genome. Based on the pig-cow alignments, we searched for similarities to 16 other organisms by UCSC available alignments, which resulted in a 87% coverage by the human genome for instance.

Conclusion

Besides recovering several of the already annotated functional RNA structures, we predicted a large number of high confidence conserved secondary structures in polyadenylated porcine transcripts. Our observations of relatively low expression levels of predicted ncRNA candidates together with the observations of higher relative amount in cDNA libraries from developmental stages are in agreement with the current paradigm of ncRNA roles in higher organisms and supports the idea of polyadenylated ncRNAs.

Treasure hunt in an amoeba: non-coding RNAs in Dictyostelium discoideum

The traditional view of RNA being merely an intermediate in the transfer of genetic information, as mRNA, spliceosomal RNA, tRNA, and rRNA, has become outdated. The recent discovery of numerous regulatory RNAs with a plethora of functions in biological processes has truly revolutionized our understanding of gene regulation. Tiny RNAs such as microRNAs and small interfering RNAs play vital roles at different levels of gene control. Small nucleolar RNAs are much more abundant than previously recognized, and new functions beyond processing and modification of rRNA have recently emerged. Longer non-coding RNAs (ncRNAs) can also have important regulatory roles in the cell, e.g., antisense RNAs that control their target mRNAs. The majority of these important findings arose from analyses in various model organisms. In this review, we focus on ncRNAs in the social amoeba Dictyostelium discoideum. This important genetically tractable model organism has recently received renewed attention in terms of discovery, regulation and functional studies of ncRNAs. Old and recent findings are discussed and put in context of what we today know about ncRNAs in other organisms.

Evidence that noncoding RNA dutA is a multicopy suppressor of Dictyostelium discoideum STAT protein Dd-STATa

Dd-STATa, a Dictyostelium discoideum homologue of metazoan STAT transcription factors, is necessary for culmination. We created a mutant strain with partial Dd-STATa activity and used it to screen for unlinked suppressor genes. We screened approximately 450,000 clones from a slug-stage cDNA library for their ability to rescue the culmination defect when overexpressed. There were 12 multicopy suppressors of Dd-STATa, of which 4 encoded segments of a known noncoding RNA, dutA. Expression of dutA is specific to the pstA zone, the region where Dd-STATa is activated. In suppressed strains the expression patterns of several putative Dd-STATa target genes become similar to the wild-type strain. In addition, the amount of the tyrosine-phosphorylated form of Dd-STATa is significantly increased in the suppressed strain. These results indicate that partial copies of dutA may act upstream of Dd-STATa to regulate tyrosine phosphorylation by an unknown mechanism.