Molecular characterization at the RNA and gene levels of U3 snoRNA from a unicellular green alga, Chlamydomonas reinhardtii

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.

High-throughput sequencing of the chloroplast and mitochondrion of Chlamydomonas reinhardtii to generate improved de novo assemblies, analyze expression patterns and transcript speciation, and evaluate diversity among laboratory strains and wild isolates

Chlamydomonas reinhardtii is a unicellular chlorophyte alga that is widely studied as a reference organism for understanding photosynthesis, sensory and motile cilia, and for development of an algal-based platform for producing biofuels and bio-products. Its highly repetitive, ~205-kbp circular chloroplast genome and ~15.8-kbp linear mitochondrial genome were sequenced prior to the advent of high-throughput sequencing technologies. Here, high coverage shotgun sequencing was used to assemble both organellar genomes de novo. These new genomes correct dozens of errors in the prior genome sequences and annotations. Genome sequencing coverage indicates that each cell contains on average 83 copies of the chloroplast genome and 130 copies of the mitochondrial genome. Using protocols and analyses optimized for organellar transcripts, RNA-Seq was used to quantify their relative abundances across 12 different growth conditions. Forty-six percent of total cellular mRNA is attributable to high expression from a few dozen chloroplast genes. RNA-Seq data were used to guide gene annotation, to demonstrate polycistronic gene expression, and to quantify splicing of psaA and psbA introns. In contrast to a conclusion from a recent study, we found that chloroplast transcripts are not edited. Unexpectedly, cytosine-rich polynucleotide tails were observed at the 3'-end of all mitochondrial transcripts. A comparative genomics analysis of eight laboratory strains and 11 wild isolates of C. reinhardtii identified 2658 variants in the organellar genomes, which is 1/10th as much genetic diversity as is found in the nucleus.

Prasinovirus Attack of Ostreococcus Is Furtive by Day but Savage by Night

Prasinoviruses are large DNA viruses that infect diverse genera of green microalgae worldwide in aquatic ecosystems, but molecular knowledge of their life cycles is lacking. Several complete genomes of both these viruses and their marine algal hosts are now available and have been used to show the pervasive presence of these species in microbial metagenomes. We have analyzed the life cycle of Ostreococcus tauri virus 5 (OtV5), a lytic virus, using transcriptome sequencing (RNA-Seq) from 12 time points of healthy or infected Ostreococcus tauri cells over a day/night cycle in culture. In the day, viral gene transcription remained low while host nitrogen metabolism gene transcription was initially strongly repressed for two successive time points before being induced for 8 h, but during the night, viral transcription increased steeply while host nitrogen metabolism genes were repressed and many host functions that are normally reduced in the dark appeared to be compensated either by genes expressed from the virus or by increased expression of a subset of 4.4% of the host's genes. Some host cells underwent lysis progressively during the night, but a larger proportion were lysed the following morning. Our data suggest that the life cycles of algal viruses mirror the diurnal rhythms of their hosts.IMPORTANCE Prasinoviruses are common in marine environments, and although several complete genomes of these viruses and their hosts have been characterized, little is known about their life cycles. Here we analyze in detail the transcriptional changes occurring over a 27-h-long experiment in a natural diurnal rhythm, in which the growth of host cells is to some extent synchronized, so that host DNA replication occurs late in the day or early in the night and cell division occurs during the night. Surprisingly, viral transcription remains quiescent over the daytime, when the most energy (from light) is available, but during the night viral transcription activates, accompanied by expression of a few host genes that are probably required by the virus. Although our experiment was accomplished in the lab, cyclical changes have been documented in host transcription in the ocean. Our observations may thus be relevant for eukaryotic phytoplankton in natural environments.

RNAiFOLD: a constraint programming algorithm for RNA inverse folding and molecular design

Synthetic biology is a rapidly emerging discipline with long-term ramifications that range from single-molecule detection within cells to the creation of synthetic genomes and novel life forms. Truly phenomenal results have been obtained by pioneering groups--for instance, the combinatorial synthesis of genetic networks, genome synthesis using BioBricks, and hybridization chain reaction (HCR), in which stable DNA monomers assemble only upon exposure to a target DNA fragment, biomolecular self-assembly pathways, etc. Such work strongly suggests that nanotechnology and synthetic biology together seem poised to constitute the most transformative development of the 21st century. In this paper, we present a Constraint Programming (CP) approach to solve the RNA inverse folding problem. Given a target RNA secondary structure, we determine an RNA sequence which folds into the target structure; i.e. whose minimum free energy structure is the target structure. Our approach represents a step forward in RNA design--we produce the first complete RNA inverse folding approach which allows for the specification of a wide range of design constraints. We also introduce a Large Neighborhood Search approach which allows us to tackle larger instances at the cost of losing completeness, while retaining the advantages of meeting design constraints (motif, GC-content, etc.). Results demonstrate that our software, RNAiFold, performs as well or better than all state-of-the-art approaches; nevertheless, our approach is unique in terms of completeness, flexibility, and the support of various design constraints. The algorithms presented in this paper are publicly available via the interactive webserver http://bioinformatics.bc.edu/clotelab/RNAiFold; additionally, the source code can be downloaded from that site.

The small subunit processome in ribosome biogenesis—progress and prospects

The small subunit (SSU) processome is a 2.2-MDa ribonucleoprotein complex involved in the processing, assembly, and maturation of the SSU of eukaryotic ribosomes. The identities of many of the factors involved in SSU biogenesis have been elucidated over the past 40 years. However, as our understanding increases, so do the number of questions about the nature of this complicated process. Cataloging the components is the first step toward understanding the molecular workings of a system. This review will focus on how identifying components of ribosome biogenesis has led to the knowledge of how these factors, protein and RNA alike, associate with one another into subcomplexes, with a concentration on the small ribosomal subunit. We will also explore how this knowledge of subcomplex assembly has informed our understanding of the workings of the ribosome synthesis system as a whole.

Features of a unique intronless cluster of class I small heat shock protein genes in tandem with box C/D snoRNA genes on chromosome 6 in tomato (Solanum lycopersicum)

Physical clustering of genes has been shown in plants; however, little is known about gene clusters that have different functions, particularly those expressed in the tomato fruit. A class I 17.6 small heat shock protein (Sl17.6 shsp) gene was cloned and used as a probe to screen a tomato (Solanum lycopersicum) genomic library. An 8.3-kb genomic fragment was isolated and its DNA sequence determined. Analysis of the genomic fragment identified intronless open reading frames of three class I shsp genes (Sl17.6, Sl20.0, and Sl20.1), the Sl17.6 gene flanked by Sl20.1 and Sl20.0, with complete 5' and 3' UTRs. Upstream of the Sl20.0 shsp, and within the shsp gene cluster, resides a box C/D snoRNA cluster made of SlsnoR12.1 and SlU24a. Characteristic C and D, and C' and D', boxes are conserved in SlsnoR12.1 and SlU24a while the upstream flanking region of SlsnoR12.1 carries TATA box 1, homol-E and homol-D box-like cis sequences, TM6 promoter, and an uncharacterized tomato EST. Molecular phylogenetic analysis revealed that this particular arrangement of shsps is conserved in tomato genome but is distinct from other species. The intronless genomic sequence is decorated with cis elements previously shown to be responsive to cues from plant hormones, dehydration, cold, heat, and MYC/MYB and WRKY71 transcription factors. Chromosomal mapping localized the tomato genomic sequence on the short arm of chromosome 6 in the introgression line (IL) 6-3. Quantitative polymerase chain reaction analysis of gene cluster members revealed differential expression during ripening of tomato fruit, and relatively different abundances in other plant parts.

Human box C/D snoRNAs with miRNA like functions: expanding the range of regulatory RNAs

Small nucleolar RNAs (snoRNAs) and microRNAs are two classes of non-protein-coding RNAs with distinct functions in RNA modification or post-transcriptional gene silencing. In this study, we introduce novel insights to RNA-induced gene activity adjustments in human cells by identifying numerous snoRNA-derived molecules with miRNA-like function, including H/ACA box snoRNAs and C/D box snoRNAs. In particular, we demonstrate that several C/D box snoRNAs give rise to gene regulatory RNAs, named sno-miRNAs here. Our data are complementing the increasing number of studies in the field of small RNAs with regulatory functions. In massively deep sequencing of small RNA fractions we identified high copy numbers of sub-sequences from >30 snoRNAs with lengths of ≥18 nt. RNA secondary structure prediction indicated for a majority of candidates a location in predicted stem regions. Experimental analysis revealed efficient gene silencing for 11 box C/D sno-miRNAs, indicating cytoplasmic processing and recruitment to the RNA silencing machinery. Assays in four different human cell lines indicated variations in both the snoRNA levels and their processing to active sno-miRNAs. In addition we show that box D elements are predominantly flanking at least one of the sno-miRNA strands, while the box C element locates within the sequence of the sno-miRNA guide strand.

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.

Eukaryotic snoRNAs: a paradigm for gene expression flexibility

Small nucleolar RNAs (snoRNAs) are one of the most ancient and numerous families of non-protein-coding RNAs (ncRNAs). The main function of snoRNAs - to guide site-specific rRNA modification - is the same in Archaea and all eukaryotic lineages. In contrast, as revealed by recent genomic and RNomic studies, their genomic organization and expression strategies are the most varied. Seemingly snoRNA coding units have adopted, in the course of evolution, all the possible ways of being transcribed, thus providing a unique paradigm of gene expression flexibility. By focusing on representative fungal, plant and animal genomes, we review here all the documented types of snoRNA gene organization and expression, and we provide a comprehensive account of snoRNA expressional freedom by precisely estimating the frequency, in each genome, of each type of genomic organization. We finally discuss the relevance of snoRNA genomic studies for our general understanding of ncRNA family evolution and expression in eukaryotes.

In vivo footprint of a picornavirus internal ribosome entry site reveals differences in accessibility to specific RNA structural elements

Internal ribosome entry site (IRES) elements were described in picornaviruses as an essential region of the viral RNA. Understanding of IRES function requires a detailed knowledge of each step involved in the internal initiation process, from RNA folding and IRES-protein interaction to ribosome recruitment. Thus, deciphering IRES accessibility to external agents due to RNA structural features, as well as RNA-protein protection within living cells, is of primary importance. In this study, two chemical reagents, dimethylsulfate (DMS) and aminomethylpsoralen, have been used to footprint the entire IRES of foot-and-mouth disease virus (FMDV) in living cells; these reagents enter the cell membrane and interact with nucleic acids in a structure-dependent manner. For FMDV, as in other picornaviruses, viral infection is dependent on the correct function of the IRES; therefore, the IRES region itself constitutes a useful target of antiviral drugs. Here, the in vivo footprint of a picornavirus IRES element in the context of a biologically active mRNA is shown for the first time. The accessibility of unpaired adenosine and cytosine nucleotides in the entire FMDV IRES was first obtained in vitro by DMS probing; subsequently, this information was used to interpret the footprint data obtained in vivo for the mRNA encompassing the IRES element in the intercistronic space. The results of DMS accessibility and UV-psoralen cross-linking studies in the competitive cellular environment provided evidence for differences in RNA structure from data obtained in vitro, and provided essential information to identify appropriate targets within the FMDV IRES aimed at combating this important pathogen.

Novel TRF1/BRF target genes revealed by genome-wide analysis of Drosophila Pol III transcription

Metazoans have evolved multiple paralogues of the TATA binding protein (TBP), adding another tunable level of gene control at core promoters. While TBP-related factor 1 (TRF1) shares extensive homology with TBP and can direct both Pol II and Pol III transcription in vitro, TRF1 target sites in vivo have remained elusive. Here, we report the genome-wide identification of TRF1-binding sites using high-resolution genome tiling microarrays. We found 354 TRF1-binding sites genome-wide with approximately 78% of these sites displaying colocalization with BRF. Strikingly, the majority of TRF1 target genes are Pol III-dependent small noncoding RNAs such as tRNAs and small nonmessenger RNAs. We provide direct evidence that the TRF1/BRF complex is functionally required for the activity of two novel TRF1 targets (7SL RNA and small nucleolar RNAs). Our studies suggest that unlike most other eukaryotic organisms that rely on TBP for Pol III transcription, in Drosophila and possibly other insects the alternative TRF1/BRF complex appears responsible for the initiation of all known classes of Pol III transcription.

INFO-RNA--a fast approach to inverse RNA folding

MOTIVATION

The structure of RNA molecules is often crucial for their function. Therefore, secondary structure prediction has gained much interest. Here, we consider the inverse RNA folding problem, which means designing RNA sequences that fold into a given structure.

RESULTS

We introduce a new algorithm for the inverse folding problem (INFO-RNA) that consists of two parts; a dynamic programming method for good initial sequences and a following improved stochastic local search that uses an effective neighbor selection method. During the initialization, we design a sequence that among all sequences adopts the given structure with the lowest possible energy. For the selection of neighbors during the search, we use a kind of look-ahead of one selection step applying an additional energy-based criterion. Afterwards, the pre-ordered neighbors are tested using the actual optimization criterion of minimizing the structure distance between the target structure and the mfe structure of the considered neighbor. We compared our algorithm to RNAinverse and RNA-SSD for artificial and biological test sets. Using INFO-RNA, we performed better than RNAinverse and in most cases, we gained better results than RNA-SSD, the probably best inverse RNA folding tool on the market.

AVAILABILITY

www.bioinf.uni-freiburg.de?Subpages/software.html.

Imp3p and Imp4p mediate formation of essential U3-precursor rRNA (pre-rRNA) duplexes, possibly to recruit the small subunit processome to the pre-rRNA

In eukaryotes, formation of short duplexes between the U3 small nucleolar RNA (snoRNA) and the precursor rRNA (pre-rRNA) at multiple sites is a prerequisite for three endonucleolytic cleavages that initiate small subunit biogenesis by releasing the 18S rRNA precursor from the pre-rRNA. The most likely role of these RNA duplexes is to guide the U3 snoRNA and its associated proteins, designated the small subunit processome, to the target cleavage sites on the pre-rRNA. Studies by others in Saccharomyces cerevisiae have identified the proteins Mpp10p, Imp3p, and Imp4p as candidates to mediate U3-pre-rRNA interactions. We report here that Imp3p and Imp4p appear to stabilize an otherwise unstable duplex between the U3 snoRNA hinge region and complementary bases in the external transcribed spacer of the pre-rRNA. In addition, Imp4p, but not Imp3p, seems to rearrange the U3 box A stem structure to expose the site that base-pairs with the 5' end of the 18S rRNA, thereby mediating duplex formation at a second site. By mediating formation of both essential U3-pre-rRNA duplexes, Imp3p and Imp4p may help the small subunit processome to dock onto the pre-rRNA, an event indispensable for ribosome biogenesis and hence for cell growth.

A new algorithm for RNA secondary structure design

The function of many RNAs depends crucially on their structure. Therefore, the design of RNA molecules with specific structural properties has many potential applications, e.g. in the context of investigating the function of biological RNAs, of creating new ribozymes, or of designing artificial RNA nanostructures. Here, we present a new algorithm for solving the following RNA secondary structure design problem: given a secondary structure, find an RNA sequence (if any) that is predicted to fold to that structure. Unlike the (pseudoknot-free) secondary structure prediction problem, this problem appears to be hard computationally. Our new algorithm, "RNA Secondary Structure Designer (RNA-SSD)", is based on stochastic local search, a prominent general approach for solving hard combinatorial problems. A thorough empirical evaluation on computationally predicted structures of biological sequences and artificially generated RNA structures as well as on empirically modelled structures from the biological literature shows that RNA-SSD substantially out-performs the best known algorithm for this problem, RNAinverse from the Vienna RNA Package. In particular, the new algorithm is able to solve structures, consistently, for which RNAinverse is unable to find solutions. The RNA-SSD software is publically available under the name of RNA Designer at the RNASoft website (www.rnasoft.ca).

Plant dicistronic tRNA-snoRNA genes: a new mode of expression of the small nucleolar RNAs processed by RNase Z

Small nucleolar RNAs (snoRNAs) guiding modifications of ribosomal RNAs and other RNAs display diverse modes of gene organization and expression depending on the eukaryotic system: in animals most are intron encoded, in yeast many are monocistronic genes and in plants most are polycistronic (independent or intronic) genes. Here we report an unprecedented organization: plant dicistronic tRNA-snoRNA genes. In Arabidopsis thaliana we identified a gene family encoding 12 novel box C/D snoRNAs (snoR43) located just downstream from tRNA(Gly) genes. We confirmed that they are transcribed, probably from the tRNA gene promoter, producing dicistronic tRNA(Gly)-snoR43 precursors. Using transgenic lines expressing a tagged tRNA-snoR43.1 gene we show that the dicistronic precursor is accurately processed to both snoR43.1 and tRNA(Gly). In addition, we show that a recombinant RNase Z, the plant tRNA 3' processing enzyme, efficiently cleaves the dicistronic precursor in vitro releasing the snoR43.1 from the tRNA(Gly). Finally, we describe a similar case in rice implicating a tRNA(Met-e) expressed in fusion with a novel C/D snoRNA, showing that this mode of snoRNA expression is found in distant plant species.

Trypanosoma brucei 5'ETS A'-cleavage is directed by 3'-adjacent sequences, but not two U3 snoRNA-binding elements, which are all required for subsequent pre-small subunit rRNA processing events

Trypanosoma brucei pre-rRNA processing commences by cleavage near the 5' end of 5.8 S sequences. The 5' external transcribed spacer (5'ETS) is removed from pre-small subunit (SSU) rRNAs by sequential cleavages at internal A' and A0 sites, and A1 at the 5' end of SSU rRNA. The A' and A0 sites positionally resemble the U3 small nucleolar RNA-dependent, primary pre-rRNA cleavages of vertebrates and yeast, respectively. Uniquely in T. brucei, two U3-crosslinkable 5'ETS sites are essential for SSU rRNA production: site1b is novel in its 3' location to the A' site, and site3 lies upstream of A0 in a position analogous to the yeast U3-binding site. Here, in vivo analysis of mutated 5'ETS sequences shows that sequences 5' to the A' site are not needed for A' cleavage or SSU rRNA production. A' cleavage is linked to, but is not sufficient to trigger, downstream pre-SSU rRNA processing events. These events require an intact 11 nt sequence, 3'-adjacent to A', which directs efficient and accurate A' cleavage. Neither the A' nearby site1b nor the site3 U3-binding elements affect A' processing, yet each is required for A0 and A1 cleavage, and SSU rRNA production. The same U3 3' hinge bases evidently bind a core element, UGUu/gGGU, within site1a and site3; the U3-site1b interaction is less reliant on base-pairing than the U3-site3 interaction. As yeast U3 5' hinge bases pair to 5'ETS sequences, it is clear that distinct U3 hinge regions can interact at both novel and related 5'ETS sites to promote 3'-proximal 5'ETS processing events in diverse organisms. The T. brucei data fit a model wherein processing factors assemble at the 5'ETS site1a to affect A' cleavage and stabilize a U3-site1b complex, which may work in concert with the downstream U3-site3 complex to assist processing events leading to ribosomal SSU production.

An unexpected, conserved element of the U3 snoRNA is required for Mpp10p association

The U3 small nucleolar ribonucleoprotein (snoRNP) is composed of a small nucleolar RNA (snoRNA) and at least 10 proteins. The U3 snoRNA base pairs with the pre-rRNA to carry out the A0, A1, and A2 processing reactions that lead to the release of the 18S rRNA from the nascent pre-rRNA transcript. The yeast U3 snoRNA can be divided into a short 5' domain (nt 1-39) and a larger 3' domain (73 to the 3' end) separated by a stretch of nucleotides called the hinge region (nt 40-72). The sequences required for pre-rRNA base pairing are found in the 5' domain and hinge region whereas the 3' domain is largely covered with proteins. Mpp10p, one of the protein components unique to the U3 snoRNP, plays a role in processing at the A1 and A2 sites. Because of its critical role in U3 snoRNP function, we determined which sequences in the U3 snoRNA are required for Mpp10p association. Unlike fibrillarin and all the previous U3 snoRNP components studied in this manner, sequences in the 3' domain are not sufficient for Mpp10p association. Instead, a conserved sequence element in the U3 snoRNA hinge region is required, placing Mpp10p near the 5' domain that carries out the pre-rRNA base-pairing interactions in the functional center of the U3 snoRNP.

Analysis of intergenic spacer transcripts suggests 'read-around' transcription of the extrachromosomal circular rDNA in Euglena gracilis

We report here the sequence of the 1743 bp intergenic spacer (IGS) that separates the 3'-end of the large subunit ribosomal RNA (rRNA) gene from the 5'-end of the small subunit (SSU) rRNA gene in the circular, extrachromosomal ribosomal DNA (rDNA) of Euglena gracilis. The IGS contains a 277 nt stretch of sequence that is related to a sequence found in ITS 1, an internal transcribed spacer between the SSU and 5.8S rRNA genes. Primer extension analysis of IGS transcripts identified three abundant reverse transcriptase stops that may be analogous to the transcription initiation site (TIS) and two processing sites (A' and A0) that are found in this region in other eukaryotes. Features that could influence processing at these sites include an imperfect palindrome near site A0 and a sequence near site A' that could potentially base pair with U3 small nucleolar RNA. Our identification of the TIS (verified by mung bean nuclease analysis) is considered tentative because we also detected low-abundance transcripts upstream of this site throughout the entire IGS. This result suggests the possibility of 'read-around' transcription, i.e. transcription that proceeds multiple times around the rDNA circle without termination.