Structural organization of the genes encoding the small nuclear RNAs U1 to U6 of Tetrahymena thermophila is very similar to that of plant small nuclear RNA genes

Spliceosomal snRNAs, the Essential Players in pre-mRNA Processing in Eukaryotic Nucleus: From Biogenesis to Functions and Spatiotemporal Characteristics

Spliceosomal small nuclear RNAs (snRNAs) are a fundamental class of non-coding small RNAs abundant in the nucleoplasm of eukaryotic cells, playing a crucial role in splicing precursor messenger RNAs (pre-mRNAs). They are transcribed by DNA-dependent RNA polymerase II (Pol II) or III (Pol III), and undergo subsequent processing and 3' end cleavage to become mature snRNAs. Numerous protein factors are involved in the transcription initiation, elongation, termination, splicing, cellular localization, and terminal modification processes of snRNAs. The transcription and processing of snRNAs are regulated spatiotemporally by various mechanisms, and the homeostatic balance of snRNAs within cells is of great significance for the growth and development of organisms. snRNAs assemble with specific accessory proteins to form small nuclear ribonucleoprotein particles (snRNPs) that are the basal components of spliceosomes responsible for pre-mRNA maturation. This article provides an overview of the biological functions, biosynthesis, terminal structure, and tissue-specific regulation of snRNAs.

Genome-wide identification and in-silico analysis of papain-family cysteine protease encoding genes in Tetrahymena thermophila

Tetrahymena thermophila is a promising host for recombinant protein production, but its utilization in biotechnology is mostly limited due to the presence of intracellular and extracellular papain-family cysteine proteases (PFCPs). In this study, we employed bioinformatics approaches to investigate the T. thermophila PFCP genes and their encoded proteases (TtPFCPs), the most prominent protease family in the genome. Results from the multiple sequence alignment, protein modeling, and conserved motif analyses revealed that all TtPFCPs showed considerably high homology with mammalian cysteine cathepsins and contained conserved amino acid motifs. The total of 121 TtPFCP-encoding genes, 14 of which were classified as non-peptidase homologs, were found. Remaining 107 true TtPFCPs were divided into four distinct subgroups depending on their homology with mammalian lysosomal cathepsins: cathepsin L-like (TtCATLs), cathepsin B-like (TtCATBs), cathepsin C-like (TtCATCs), and cathepsin X-like (TtCATXs) PFCPs. The majority of true TtPFCPs (96 out of the total) were in TtCATL-like peptidase subgroup. Both phylogenetic and chromosomal localization analyses of TtPFCPs supported the hypothesis that TtPFCPs likely evolved through tandem gene duplication events and predominantly accumulated on micronuclear chromosome 5. Additionally, more than half of the identified TtPFCP genes are expressed in considerably low quantities compared to the rest of the TtPFCP genes, which are expressed at a higher level. However, their expression patterns fluctuate based on the stage of the life cycle. In conclusion, this study provides the first comprehensive in-silico analysis of TtPFCP genes and encoded proteases. The results would help designing an effective strategy for protease knockout mutant cell lines to discover biological function and to improve the recombinant protein production in T. thermophila.

Evolution of plant telomerase RNAs: farther to the past, deeper to the roots

The enormous sequence heterogeneity of telomerase RNA (TR) subunits has thus far complicated their characterization in a wider phylogenetic range. Our recent finding that land plant TRs are, similarly to known ciliate TRs, transcribed by RNA polymerase III and under the control of the type-3 promoter, allowed us to design a novel strategy to characterize TRs in early diverging Viridiplantae taxa, as well as in ciliates and other Diaphoretickes lineages. Starting with the characterization of the upstream sequence element of the type 3 promoter that is conserved in a number of small nuclear RNAs, and the expected minimum TR template region as search features, we identified candidate TRs in selected Diaphoretickes genomes. Homologous TRs were then used to build covariance models to identify TRs in more distant species. Transcripts of the identified TRs were confirmed by transcriptomic data, RT-PCR and Northern hybridization. A templating role for one of our candidates was validated in Physcomitrium patens. Analysis of secondary structure demonstrated a deep conservation of motifs (pseudoknot and template boundary element) observed in all published TRs. These results elucidate the evolution of the earliest eukaryotic TRs, linking the common origin of TRs across Diaphoretickes, and underlying evolutionary transitions in telomere repeats.

Telomerase RNAs in land plants

To elucidate the molecular nature of evolutionary changes of telomeres in the plant order Asparagales, we aimed to characterize telomerase RNA subunits (TRs) in these plants. The unusually long telomere repeat unit in Allium plants (12 nt) allowed us to identify TRs in transcriptomic data of representative species of the Allium genus. Orthologous TRs were then identified in Asparagales plants harbouring telomere DNA composed of TTAGGG (human type) or TTTAGGG (Arabidopsis-type) repeats. Further, we identified TRs across the land plant phylogeny, including common model plants, crop plants, and plants with unusual telomeres. Several lines of functional testing demonstrate the templating telomerase function of the identified TRs and disprove a functionality of the only previously reported plant telomerase RNA in Arabidopsis thaliana. Importantly, our results change the existing paradigm in plant telomere biology which has been based on the existence of a relatively conserved telomerase reverse transcriptase subunit (TERT) associating with highly divergent TRs even between closely related plant taxa. The finding of a monophyletic origin of genuine TRs across land plants opens the possibility to identify TRs directly in transcriptomic or genomic data and/or predict telomere sequences synthesized according to the respective TR template region.

Knock-Down of a Novel snoRNA in Tetrahymena Reveals a Dual Role in 5.8S rRNA Processing and Generation of a 26S rRNA Fragment

In eukaryotes, 18S, 5.8S, and 28S rRNAs are transcribed as precursor molecules that undergo extensive modification and nucleolytic processing to form the mature rRNA species. Central in the process are the small nucleolar RNAs (snoRNAs). The majority of snoRNAs guide site specific chemical modifications but a few are involved in defining pre-rRNA cleavages. Here, we describe an unusual snoRNA (TtnuCD32) belonging to the box C/D subgroup from the ciliate Tetrahymena thermophila. We show that TtnuCD32 is unlikely to function as a modification guide snoRNA and that it is critical for cell viability. Cell lines with genetic knock-down of TtnuCD32 were impaired in growth and displayed two novel and apparently unrelated phenotypes. The most prominent phenotype is the accumulation of processing intermediates of 5.8S rRNA. The second phenotype is the decrease in abundance of a ~100 nt 26S rRNA fragment of unknown function. Sequence analysis demonstrated that TtnuCD32 share features with the essential snoRNA U14 but an alternative candidate (TtnuCD25) was more closely related to other U14 sequences. This, together with the fact that the observed rRNA processing phenotypes were not similar to what has been observed in U14 depleted cells, suggests that TtnuCD32 is a U14 homolog that has gained novel functions.

Budding yeast telomerase RNA transcription termination is dictated by the Nrd1/Nab3 non-coding RNA termination pathway

The RNA component of budding yeast telomerase (Tlc1) occurs in two forms, a non-polyadenylated form found in functional telomerase and a rare polyadenylated version with unknown function. Previous work suggested that the functional Tlc1 polyA- RNA is processed from the polyA+ form, but the mechanisms regulating its transcription termination and 3'-end formation remained unclear. Here we examined transcription termination of Tlc1 RNA in the sequences 3' of the TLC1 gene and relate it to telomere maintenance. Strikingly, disruption of all probable or cryptic polyadenylation signals near the 3'-end blocked the accumulation of the previously reported polyA+ RNA without affecting the level, function or specific 3' nucleotide of the mature polyA- form. A genetic approach analysing TLC1 3'-end sequences revealed that transcription terminates upstream of the polyadenylation sites. Furthermore, the results also demonstrate that the function of this Tlc1 terminator depends on the Nrd1/Nab3 transcription termination pathway. The data thus show that transcription termination of the budding yeast telomerase RNA occurs as that of snRNAs and Tlc1 functions in telomere maintenance are not strictly dependent on a polyadenylated precursor, even if the polyA+ form can serve as intermediate in a redundant termination/maturation pathway.

Experimental identification and analysis of macronuclear non-coding RNAs from the ciliate Tetrahymena thermophila

The ciliate Tetrahymena thermophila is an important eukaryotic model organism that has been used in pioneering studies of general phenomena, such as ribozymes, telomeres, chromatin structure and genome reorganization. Recent work has shown that Tetrahymena has many classes of small RNA molecules expressed during vegetative growth or sexual reorganization. In order to get an overview of medium-sized (40-500 nt) RNAs expressed from the Tetrahymena genome, we created a size-fractionated cDNA library from macronuclear RNA and analyzed 80 RNAs, most of which were previously unknown. The most abundant class was small nucleolar RNAs (snoRNAs), many of which are formed by an unusual maturation pathway. The modifications guided by the snoRNAs were analyzed bioinformatically and experimentally and many Tetrahymena-specific modifications were found, including several in an essential, but not conserved domain of ribosomal RNA. Of particular interest, we detected two methylations in the 5'-end of U6 small nuclear RNA (snRNA) that has an unusual structure in Tetrahymena. Further, we found a candidate for the first U8 outside metazoans, and an unusual U14 candidate. In addition, a number of candidates for new non-coding RNAs were characterized by expression analysis at different growth conditions.

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.

Genome-wide evolutionary analysis of the noncoding RNA genes and noncoding DNA of Paramecium tetraurelia

The compact genome of the unicellular eukaryote Paramecium tetraurelia contains noncoding DNA (ncDNA) distributed into >39,000 intergenic sequences and >90,000 introns of 390 base pairs (bp) and 25 bp on average, respectively. Here we analyzed the molecular features of the ncRNA genes, introns, and intergenic sequences of this genome. We mainly used computational programs and comparative genomics possible because the P. tetraurelia genome had formed throughout whole-genome duplications (WGDs). We characterized 417 5S rRNA, snRNA, snoRNA, SRP RNA, and tRNA putative genes, 415 of which map within intergenic sequences, and two, within introns. The evolution of these ncRNA genes appears to have mainly involved purifying selection and gene deletion. We then compared the introns that interrupt the protein-coding gene duplicates arisen from the recent WGD and identified a population of a few thousands of introns having evolved under most stringent constraints (>95% of identity). We also showed that low nucleotide substitution levels characterize the 50 and 80-115 base pairs flanking, respectively, the stop and start codons of the protein-coding genes. Lower substitution levels mark the base pairs flanking the highly transcribed genes, or the start codons of the genes of the sets with a high number of WGD-related sequences. Finally, adjacent to protein-coding genes, we characterized 32 DNA motifs able to encode stable and evolutionary conserved RNA secondary structures and defining putative expression controlling elements. Fourteen DNA motifs with similar properties map distant from protein-coding genes and may encode regulatory 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.

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

Most eukaryotic mRNAs depend upon precise removal of introns by the spliceosome, a complex of RNAs and proteins. Splicing of pre-mRNA is known to take place in Dictyostelium discoideum, and we previously isolated the U2 spliceosomal RNA experimentally. In this study, we identified the remaining major spliceosomal RNAs in Dictyostelium by a bioinformatical approach. Expression was verified from 17 small nuclear RNA (snRNA) genes. All these genes are preceded by a putative noncoding RNA gene promoter. Immunoprecipitation showed that snRNAs U1, U2, U4, and U5, but not U6, carry the conserved trimethylated 5' cap structure. A number of divergent U2 species are expressed in Dictyostelium. These RNAs carry the U2 RNA hallmark sequence and structure motifs but have an additional predicted stem-loop structure at the 5' end. Surprisingly, and in contrast to the other spliceosomal RNAs in this study, the new U2 variants were enriched in the cytoplasm and were developmentally regulated. Furthermore, all of the snRNAs could also be detected as polyadenylated species, and polyadenylated U1 RNA was demonstrated to be located in the cytoplasm.

Identification of an essential proximal sequence element in the promoter of the telomerase RNA gene of Tetrahymena thermophila

Telomerase is a ribonucleoprotein reverse transcriptase that synthesizes and maintains telomeric DNA. Studies of telomeres and telomerase are facilitated by the large number of linear DNA molecules found in ciliated protozoa, such as Tetrahymena thermophila. To examine the expression of telomerase, we investigated the transcription of the RNA polymerase III-directed gene encoding the RNA subunit (TER1) of this enzyme. A chimeric gene containing the Glaucoma chattoni TER1 transcribed region flanked by 5' and 3' Tetrahymena regions was used to identify promoter elements following transformation of Tetrahymena cells. Disruption of a conserved proximal sequence element (PSE) located at -55 in the Tetrahymena TER1 5' flanking region eliminated expression of the chimeric gene. In addition, mutation of an A/T-rich element at -25 decreased expression markedly. A gel mobility shift assay and protein-DNA cross-linking identified a PSE-binding polypeptide of 50-60 kDa in Tetrahymena extracts. Gel filtration analysis revealed a native molecular mass of approximately 160 kDa for this binding activity. Our results point to a similar architecture between ciliate telomerase RNA and metazoan U6 small nuclear RNA promoters.

Telomerase reverse transcriptase genes identified in Tetrahymena thermophila and Oxytricha trifallax

Telomerase reverse transcriptase (TERT) has been identified as the catalytic subunit of the chromosome end-replicating enzyme in Euplotes, yeasts, and mammals. However, it was not reported among the protein components of purified Tetrahymena telomerase, the first telomerase identified and the most thoroughly studied. It therefore seemed possible that Tetrahymena used an alternative telomerase that lacked a TERT protein. We now report the cloning and sequencing of a Tetrahymena thermophila gene whose encoded protein has the properties expected for a TERT, including large size (133 kDa), basicity (calculated pI = 10.0), and reverse transcriptase sequence motifs with telomerase-specific features. The expression of mRNA from the Tetrahymena TERT gene increases dramatically at 2-5 h after conjugation, preceding de novo addition of telomeres to macronuclear DNA molecules. We also report the cloning and sequencing of the ortholog from Oxytricha trifallax. The Oxytricha macronuclear TERT gene has no introns, whereas that of Tetrahymena has 18 introns. Sequence comparisons reveal a new amino acid sequence motif (CP), conserved among the ciliated protozoan TERTs, and allow refinement of previously identified motifs. A phylogenetic tree of the known TERTs follows the phylogeny of the organisms in which they are found, consistent with an ancient origin rather than recent transposition. The conservation of TERTs among eukaryotes supports the model that telomerase has a conserved core (TERT plus the RNA subunit), with other subunits of the holoenzyme being more variable among species.

The mouse telomerase RNA 5"-end lies just upstream of the telomerase template sequence

Telomerase is a ribonucleoprotein enzyme with an essential RNA component. Embedded within the telomerase RNA is a template sequence for telomere synthesis. We have characterized the structure of the 5' regions of the human and mouse telomerase-RNA genes, and have found a striking difference in the location of the template sequence: Whereas the 5'-end of the human telomerase RNA lies 45 nt from the telomerase-RNA template sequence, the 5'-end of the mouse telomerase RNA lies just 2 nt from the telomerase-RNA template sequence. Analysis of genomic sequences flanking the 5'-end of the human and mouse telomerase RNA-coding sequences reveals similar promoter-element arrangements typical of mRNA-type promoters: a TATA box-like element and an upstream region containing a consensus CCAAT box. This putative promoter structure contrasts with that of the ciliate telomerase-RNA genes whose structure resembles RNA polymerase III U6 small nuclear RNA (snRNA) promoters. These and other comparisons suggest that, during evolution, both the RNA-polymerase specificity of telomerase RNA-gene promoters and, more recently, the position of the template sequence in the telomerase RNA changed.

Expression of Tetrahymena snRNA gene variants including a U1 gene with mutations in the 5' splice site recognition sequence

The expression of U1, U2 and U5 snRNA gene variants has been studied under different physiological states of Tetrahymena. Variants of all three snRNA genes are expressed. Among the snRNAs detected is U1-3, a variant with 66 mutations compared to the normal U1 snRNA. Three of these mutations affect the 5' splice site recognition sequence. The U1-3 snRNA is present in a few hundred copies per cell. The expression of Tetrahymena snRNA genes is independent of the physiological state of the cell.

Identification of developmentally regulated sea urchin U5 snRNA genes

A PCR approach was used to isolate repeated U5 small nuclear RNA (snRNA) genes from the sea urchin Lytechinus variegatus. A 1.3 kb repeat, LvU5.0, and three other variants, LvU5.1-U5.3, that differ in the coding region and in the proximal sequence element (PSE) region were isolated. Southern Blot analysis indicate that the U5 snRNA genes, unlike other embryonically expressed snRNA genes (U1, U2 and U6), are not found in a simple tandem repeat, but instead, exist in several heterogeneous clusters each with a small number of genes. The U5 PSE has limited sequence similarity with the other sea urchin PSEs. However, when used in a mobility shift assay the U5 PSE forms a protein/DNA complex that is very similar to the complex formed with the U6 PSE. An RNase protection assay used to monitor the accumulation of U5 snRNA during development shows that at least two U5 variants are coordinately expressed during embryogenesis.

A single telomerase RNA is sufficient for the synthesis of variable telomeric DNA repeats in ciliates of the genus Paramecium

Paramecium telomeric DNA consists largely of a random distribution of TTGGGG and TTTGGG repeats. Given the precise nature of other ciliate telomerases, it has been postulated that there are two distinct types of the Paramecium enzyme, each synthesizing perfect telomeric repeats: one with a template RNA that specifies the addition of TTTGGG and the second dictating the synthesis of TTGGGG repeats. We have cloned and sequenced telomerase RNA genes from Paramecium tetraurelia, P. primaurelia, P. multimicronucleatum, and P. caudatum. Surprisingly, a single gene encodes telomerase RNA in all four species, although an apparently nontranscribed pseudogene is also present in the genome of P. primaurelia. The overall lengths of the telomerase RNAs range between 202 and 209 nucleotides, and they can be folded into a conserved secondary structure similar to that derived for other ciliate RNAs. All Paramecium telomerase RNAs examined include a template specific for the synthesis of TTGGGG telomeric repeats, which has not been posttranscriptionally edited to account for the conventional synthesis of TTTGGG repeats. On the basis of these results, possible mechanisms for the synthesis of variable telomeric repeats by Paramecium telomerase are discussed.

A unique U5-->A substitution in the Physarum polycephalum U1 snRNA: evidence at the RNA and gene levels

The 5' terminal sequence of U1 snRNA that base-pairs with the intron 5' splice site in the course of spliceosome assembly was considered to be universally conserved. A study of the P polycephalum U1 snRNA at both RNA and gene levels shows that there are exceptions to this rule: the P polycephalum U1 snRNA has a U to A substitution at position 5, that is partially compensated by a high frequency of T residue at position +4 of introns. In contrast to the yeast genome, the P polycephalum genome contains several U1 snRNA coding sequences (about 20). They either encode the U1A snRNA expressed in microplasmodia or correspond to the previously cloned U1B coding sequence. Both coding sequences show the U5A substitution. The ratio of U1A versus U1B coding sequences is of about 3. A U1A gene was cloned. The 60 nt region upstream of the coding sequence has the same sequence as in the U1B gene. The U1B gene is probably expressed at another stage of the P polycephalum life cycle.

Ciliate telomerase RNA structural features

Telomerase RNA is an integral part of telomerase, the ribonucleoprotein enzyme that catalyzes the synthesis of telomeric DNA. The RNA moiety contains a templating domain that directs the synthesis of a species-specific telomeric repeat and may also be important for enzyme structure and/or catalysis. Phylogenetic comparisons of telomerase RNA sequences from various Tetrahymena spp. and hypotrich ciliates have revealed two conserved secondary structure models that share many features. We have cloned and sequenced the telomerase RNA genes from an additional six Tetrahymena spp. (T. vorax, T. borealis, T. australis, T. silvana, T. capricornis and T. paravorax). Inclusion of these sequences, most notably that from T. paravorax, in a phylogenetic comparative analysis allowed us to more narrowly define structural elements that may be necessary for a minimal telomerase RNA. A primary sequence element, positioned 5' of the template and conserved between all previously known ciliate telomerase RNAs, has been reduced from 5'-(C)UGUCA-3' to the 4 nt sequence 5'-GUCA-3'. Conserved secondary structural features and the impact they have on the general organization of ciliate telomerase RNAs is discussed.

Telomerase RNAs of different ciliates have a common secondary structure and a permuted template

Telomerase is composed of protein and RNA. The RNA serves as a template for telomere DNA synthesis and may also be important for enzyme structure or catalysis. We have used the presence of conserved sequence elements in the promoter and template regions to amplify by PCR the telomerase RNA genes from six different hypotrichous ciliates: Oxytricha nova, Oxytricha trifallax, Stylonychia mytilis, Stylonychia lemnae, Euplotes aediculatus, and Euplotes eurystomus. RNaseH cleavage of the O. nova RNA in extracts by use of a complementary oligonucleotide leads to loss of telomerase activity, supporting the identification. Primary sequence and biochemical experiments suggest that the templates of Oxytricha and Stylonychia are circularly permuted relative to that of E. aediculatus. On the basis of the pause sites, the former two add G4T4 during a single primer elongation cycle, whereas E. aediculatus adds G3T4G. The only primary sequence element outside the template that is conserved between these phylogenetically distant telomerase RNAs is the sequence 5'-(C)UGUCA-3', which precedes the template regions by exactly two bases. We propose a common secondary structure model that is based on nucleotide covariations, a model which resembles that proposed previously for tetrahymenine telomerase RNAs.