Global identification of noncoding RNAs in Saccharomyces cerevisiae by modulating an essential RNA processing pathway

Multiple structural flavors of RNase P in precursor tRNA processing

The precursor transfer RNAs (pre-tRNAs) require extensive processing to generate mature tRNAs possessing proper fold, structural stability, and functionality required to sustain cellular viability. The road to tRNA maturation follows an ordered process: 5'-processing, 3'-processing, modifications at specific sites, if any, and 3'-CCA addition before aminoacylation and recruitment to the cellular protein synthesis machinery. Ribonuclease P (RNase P) is a universally conserved endonuclease in all domains of life, performing the hydrolysis of pre-tRNA sequences at the 5' end by the removal of phosphodiester linkages between nucleotides at position -1 and +1. Except for an archaeal species: Nanoarchaeum equitans where tRNAs are transcribed from leaderless-position +1, RNase P is indispensable for life and displays fundamental variations in terms of enzyme subunit composition, mechanism of substrate recognition and active site architecture, utilizing in all cases a two metal ion-mediated conserved catalytic reaction. While the canonical RNA-based ribonucleoprotein RNase P has been well-known to occur in bacteria, archaea, and eukaryotes, the occurrence of RNA-free protein-only RNase P in eukaryotes and RNA-free homologs of Aquifex RNase P in prokaryotes has been discovered more recently. This review aims to provide a comprehensive overview of structural diversity displayed by various RNA-based and RNA-free RNase P holoenzymes towards harnessing critical RNA-protein and protein-protein interactions in achieving conserved pre-tRNA processing functionality. Furthermore, alternate roles and functional interchangeability of RNase P are discussed in the context of its employability in several clinical and biotechnological applications. This article is categorized under: RNA Processing > tRNA Processing RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes.

The specificity landscape of bacterial ribonuclease P

Developing quantitative models of substrate specificity for RNA processing enzymes is a key step toward understanding their biology and guiding applications in biotechnology and biomedicine. Optimally, models to predict relative rate constants for alternative substrates should integrate an understanding of structures of the enzyme bound to "fast" and "slow" substrates, large datasets of rate constants for alternative substrates, and transcriptomic data identifying in vivo processing sites. Such data are either available or emerging for bacterial ribonucleoprotein RNase P a widespread and essential tRNA 5' processing endonuclease, thus making it a valuable model system for investigating principles of biological specificity. Indeed, the well-established structure and kinetics of bacterial RNase P enabled the development of high throughput measurements of rate constants for tRNA variants and provided the necessary framework for quantitative specificity modeling. Several studies document the importance of conformational changes in the precursor tRNA substrate as well as the RNA and protein subunits of bacterial RNase P during binding, although the functional roles and dynamics are still being resolved. Recently, results from cryo-EM studies of E. coli RNase P with alternative precursor tRNAs are revealing prospective mechanistic relationships between conformational changes and substrate specificity. Yet, extensive uncharted territory remains, including leveraging these advances for drug discovery, achieving a complete accounting of RNase P substrates, and understanding how the cellular context contributes to RNA processing specificity in vivo.

An Overview on Identification and Regulatory Mechanisms of Long Non-coding RNAs in Fungi

For decades, more and more long non-coding RNAs (lncRNAs) have been confirmed to play important functions in key biological processes of different organisms. At present, most identified lncRNAs and those with known functional roles are from mammalian systems. However, lncRNAs have also been found in primitive eukaryotic fungi, and they have different functions in fungal development, metabolism, and pathogenicity. In this review, we highlight some recent researches on lncRNAs in the primitive eukaryotic fungi, particularly focusing on the identification of lncRNAs and their regulatory roles in diverse biological processes.

Advances in RNAi-Assisted Strain Engineering in Saccharomyces cerevisiae

Saccharomyces cerevisiae is a widely used eukaryotic model and microbial cell factory. RNA interference (RNAi) is a conserved regulatory mechanism among eukaryotes but absent from S. cerevisiae. Recent reconstitution of RNAi machinery in S. cerevisiae enables the use of this powerful tool for strain engineering. Here we first discuss the introduction of heterologous RNAi pathways in S. cerevisiae, and the design of various expression cassettes of RNAi precursor reagents for tunable, dynamic, and genome-wide regulation. We then summarize notable examples of RNAi-assisted functional genomics and metabolic engineering studies in S. cerevisiae. We conclude with the future challenges and opportunities of RNAi-based approaches, as well as the potential of other regulatory RNAs in advancing yeast engineering.

RNAseq-based transcriptome comparison of Saccharomyces cerevisiae strains isolated from diverse fermentative environments

Transcriptome analyses play a central role in unraveling the complexity of gene expression regulation in Saccharomyces cerevisiae. This species, one of the most important microorganisms for humans given its industrial applications, shows an astonishing degree of genetic and phenotypic variability among different strains adapted to specific environments. In order to gain novel insights into the Saccharomyces cerevisiae biology of strains adapted to different fermentative environments, we analyzed the whole transcriptome of three strains isolated from wine, flor wine or mezcal fermentations. An RNA-seq transcriptome comparison of the different yeasts in the samples obtained during synthetic must fermentation highlighted the differences observed in the genes that encode mannoproteins, and in those involved in aroma, sugar transport, glycerol and alcohol metabolism, which are important under alcoholic fermentation conditions. These differences were also observed in the physiology of the strains after mannoprotein and aroma determinations. This study offers an essential foundation for understanding how gene expression variations contribute to the fermentation differences of the strains adapted to unequal fermentative environments. Such knowledge is crucial to make improvements in fermentation processes and to define targets for the genetic improvement or selection of wine yeasts.

Transcriptional and Posttranscriptional Programming by Long Noncoding RNAs

Recently, several lines of evidence have suggested that noncoding RNAs, which include both small and long noncoding RNAs (ncRNAs), contribute to a significant portion of the transcriptome in eukaryotic organisms. However, the functional significance of this wide-spread occurrence of ncRNAs, and in particular, the long ncRNAs (lncRNAs), for organismal development and differentiation is unclear. The available evidence from a subset of lncRNAs suggests that certain lncRNAs, and/or the act of their transcription, are involved in important biological functions at the transcriptional and posttranscriptional level. This chapter discusses the epigenetic and nonepigenetic mechanisms by which lncRNAs and/or their transcription are involved in the programming of various biological functions in model systems, from yeast to mammals.

Sequence Analysis and Comparative Study of the Protein Subunits of Archaeal RNase P

RNase P, a ribozyme-based ribonucleoprotein (RNP) complex that catalyzes tRNA 5′-maturation, is ubiquitous in all domains of life, but the evolution of its protein components (RNase P proteins, RPPs) is not well understood. Archaeal RPPs may provide clues on how the complex evolved from an ancient ribozyme to an RNP with multiple archaeal and eukaryotic (homologous) RPPs, which are unrelated to the single bacterial RPP. Here, we analyzed the sequence and structure of archaeal RPPs from over 600 available genomes. All five RPPs are found in eight archaeal phyla, suggesting that these RPPs arose early in archaeal evolutionary history. The putative ancestral genomic loci of archaeal RPPs include genes encoding several members of ribosome, exosome, and proteasome complexes, which may indicate coevolution/coordinate regulation of RNase P with other core cellular machineries. Despite being ancient, RPPs generally lack sequence conservation compared to other universal proteins. By analyzing the relative frequency of residues at every position in the context of the high-resolution structures of each of the RPPs (either alone or as functional binary complexes), we suggest residues for mutational analysis that may help uncover structure-function relationships in RPPs.

Identification of novel noncoding transcripts in telomerase-negative yeast using RNA-seq

Telomerase is a ribonucleoprotein that maintains the ends of linear chromosomes in most eukaryotes. Loss of telomerase activity results in shortening of telomeric DNA and eventually a specific G2/M cell-cycle arrest known as senescence. In humans, telomere shortening occurs during aging, while inappropriate activation of telomerase is associated with approximately 90% of cancers. Previous studies have identified several classes of noncoding RNAs (ncRNA) also associated with aging-related senescence and cancer, but whether ncRNAs are also involved in short-telomere-induced senescence in yeast is unknown. Here, we report 112 putative novel lncRNAs in the yeast Saccharomyces cerevisiae, 41 of which are only expressed in telomerase-negative yeast. Expression of approximately half of the lncRNAs is strongly correlated with that of adjacent genes, suggesting this subset may influence transcription of neighboring genes. Our results reveal a new potential mechanism governing adaptive changes in senescing and post-senescent survivor yeast cells.

Time resolved DNA occupancy dynamics during the respiratory oscillation uncover a global reset point in the yeast growth program

The structural dynamics of chromatin have been implicated in the regulation of fundamental eukaryotic processes, such as DNA transcription, replication and repair. Although previous studies have revealed that the chromatin landscape, nucleosome remodeling and histone modification events are intimately tied into cellular energetics and redox state, few studies undertake defined time-resolved measurements of these state variables. Here, we use metabolically synchronous, continuously-grown yeast cultures to measure DNA occupancy and track global patterns with respect to the metabolic state of the culture. Combined with transcriptome analyses and ChIP-qPCR experiments, these paint an intriguing picture where genome-wide nucleosome focusing occurs during the recovery of energy charge, followed by clearance of the promoter regions and global transcriptional slow-down, thus indicating a nucleosome-mediated “reset point” for the cycle. The reset begins at the end of the catabolic and stress-response transcriptional programs and ends prior to the start of the anabolic and cell-growth transcriptional program, and the histones on genes from both the catabolic and anabolic superclusters are deacetylated.

Playing RNase P evolution: swapping the RNA catalyst for a protein reveals functional uniformity of highly divergent enzyme forms

The RNase P family is a diverse group of endonucleases responsible for the removal of 5′ extensions from tRNA precursors. The diversity of enzyme forms finds its extremes in the eukaryal nucleus where RNA-based catalysis by complex ribonucleoproteins in some organisms contrasts with single-polypeptide enzymes in others. Such structural contrast suggests associated functional differences, and the complexity of the ribonucleoprotein was indeed proposed to broaden the enzyme's functionality beyond tRNA processing. To explore functional overlap and differences between most divergent forms of RNase P, we replaced the nuclear RNase P of Saccharomyces cerevisiae, a 10-subunit ribonucleoprotein, with Arabidopsis thaliana PRORP3, a single monomeric protein. Surprisingly, the RNase P-swapped yeast strains were viable, displayed essentially unimpaired growth under a wide variety of conditions, and, in a certain genetic background, their fitness even slightly exceeded that of the wild type. The molecular analysis of the RNase P-swapped strains showed a minor disturbance in tRNA metabolism, but did not point to any RNase P substrates or functions beyond that. Altogether, these results indicate the full functional exchangeability of the highly dissimilar enzymes. Our study thereby establishes the RNase P family, with its combination of structural diversity and functional uniformity, as an extreme case of convergent evolution. It moreover suggests that the apparently gratuitous complexity of some RNase P forms is the result of constructive neutral evolution rather than reflecting increased functional versatility.

Many biocatalysts apparently evolved independently more than once, leading to structurally unrelated macromolecules catalyzing the same biochemical reaction. The RNase P enzyme family is an exceptional case of this phenomenon called convergent evolution. RNase P enzymes use not only unrelated, but chemically distinct macromolecules, either RNA or protein, to catalyze a specific step in the biogenesis of transfer RNAs, the ubiquitous adaptor molecules in protein synthesis. However, this fundamental difference in the identity of the actual catalyst, together with a broad variation in structural complexity of the diverse forms of RNase P, cast doubts on their functional equivalence. Here we compared two of the structurally most extreme variants of RNase P by replacing the yeast nuclear enzyme, a 10-subunit RNA-protein complex, with a single-protein from plants representing the apparently simplest form of RNase P. Surprisingly, the viability and fitness of these RNase P-swapped yeasts and their molecular analyses demonstrated the full functional exchangeability of the highly dissimilar enzymes. The RNase P family, with its combination of structural diversity and functional uniformity, thus not only truly represents an extraordinary case of convergent evolution, but also demonstrates that increased structural complexity does not necessarily entail broadened functionality, but may rather be the result of “neutral” evolutionary mechanisms.

Understanding non-coding DNA regions in yeast

Non-coding transcripts play an important role in gene expression regulation in all species, including budding and fission yeast. Such regulatory transcripts include intergenic ncRNA (non-coding RNA), 5′ and 3′ UTRs, introns and antisense transcripts. In the present review, we discuss advantages and limitations of recently developed sequencing techniques, such as ESTs, DNA microarrays, RNA-Seq (RNA sequencing), DRS (direct RNA sequencing) and TIF-Seq (transcript isoform sequencing). We provide an overview of methods applied in yeast and how each of them has contributed to our knowledge of gene expression regulation and transcription.

RNase P enzymes: divergent scaffolds for a conserved biological reaction

Ribonuclease P (RNase P) catalyzes the maturation of the 5' end of precursor-tRNAs (pre-tRNA) and is conserved in all domains of life. However, the composition of RNase P varies from bacteria to archaea and eukarya, making RNase P one of the most diverse enzymes characterized. Most known RNase P enzymes contain a large catalytic RNA subunit that associates with one to 10 proteins. Recently, a protein-only form of RNase P was discovered in mitochondria and chloroplasts of many higher eukaryotes. This proteinaceous RNase P (PRORP) represents a new class of metallonucleases. Here we discuss our recent crystal structure of PRORP1 from Arabidopsis thaliana and speculate on the reasons for the replacement of catalytic RNA by a protein catalyst. We conclude, based on an analysis of the catalytic efficiencies of ribonucleoprotein (RNP) and PRORP enzymes, that the need for greater catalytic efficiency is most likely not the driving force behind the replacement of the RNA with a protein catalyst. The emergence of a protein-based RNase P more likely reflects the increasing complexity of the biological system, including difficulties in importation into organelles and vulnerability of organellar RNAs to cleavage.

Correction of spatial bias in oligonucleotide array data

Background. Oligonucleotide microarrays allow for high-throughput gene expression profiling assays. The technology relies on the fundamental assumption that observed hybridization signal intensities (HSIs) for each intended target, on average, correlate with their target's true concentration in the sample. However, systematic, nonbiological variation from several sources undermines this hypothesis. Background hybridization signal has been previously identified as one such important source, one manifestation of which appears in the form of spatial autocorrelation. Results. We propose an algorithm, pyn, for the elimination of spatial autocorrelation in HSIs, exploiting the duality of desirable mutual information shared by probes in a common probe set and undesirable mutual information shared by spatially proximate probes. We show that this correction procedure reduces spatial autocorrelation in HSIs; increases HSI reproducibility across replicate arrays; increases differentially expressed gene detection power; and performs better than previously published methods. Conclusions. The proposed algorithm increases both precision and accuracy, while requiring virtually no changes to users' current analysis pipelines: the correction consists merely of a transformation of raw HSIs (e.g., CEL files for Affymetrix arrays). A free, open-source implementation is provided as an R package, compatible with standard Bioconductor tools. The approach may also be tailored to other platform types and other sources of bias.

Analysis of tiling array expression studies with flexible designs in Bioconductor (waveTiling)

BACKGROUND

Existing statistical methods for tiling array transcriptome data either focus on transcript discovery in one biological or experimental condition or on the detection of differential expression between two conditions. Increasingly often, however, biologists are interested in time-course studies, studies with more than two conditions or even multiple-factor studies. As these studies are currently analyzed with the traditional microarray analysis techniques, they do not exploit the genome-wide nature of tiling array data to its full potential.

RESULTS

We present an R Bioconductor package, waveTiling, which implements a wavelet-based model for analyzing transcriptome data and extends it towards more complex experimental designs. With waveTiling the user is able to discover (1) group-wise expressed regions, (2) differentially expressed regions between any two groups in single-factor studies and in (3) multifactorial designs. Moreover, for time-course experiments it is also possible to detect (4) linear time effects and (5) a circadian rhythm of transcripts. By considering the expression values of the individual tiling probes as a function of genomic position, effect regions can be detected regardless of existing annotation. Three case studies with different experimental set-ups illustrate the use and the flexibility of the model-based transcriptome analysis.

CONCLUSIONS

The waveTiling package provides the user with a convenient tool for the analysis of tiling array trancriptome data for a multitude of experimental set-ups. Regardless of the study design, the probe-wise analysis allows for the detection of transcriptional effects in both exonic, intronic and intergenic regions, without prior consultation of existing annotation.

Nuclear noncoding RNA surveillance: is the end in sight?

Nuclear noncoding RNA (ncRNA) surveillance pathways play key roles in shaping the steady-state transcriptomes of eukaryotic cells. Defective and unneeded ncRNAs are primarily degraded by exoribonucleases that rely on protein cofactors to identify these RNAs. Recent studies have begun to elucidate both the mechanisms by which these cofactors recognize aberrant RNAs and the features that mark RNAs for degradation. One crucial RNA determinant is the presence of an accessible end; in addition, the failure of aberrant RNAs to fold into compact structures and assemble with specific binding proteins probably also contributes to their recognition and subsequent degradation. To date, ncRNA surveillance has been most extensively studied in budding yeast. However, mammalian cells possess nucleases and cofactors that have no known yeast counterparts, indicating that RNA surveillance pathways may be more complex in metazoans. Importantly, there is evidence that the failure of ncRNA surveillance pathways contributes to human disease.

Transcriptional regulatory circuitries in the human pathogen Candida albicans involving sense--antisense interactions

Candida albicans, a major human fungal pathogen, usually contains a diploid genome, but controls adaptation to a toxic alternative carbon source L-sorbose, by the reversible loss of one chromosome 5 (Ch5). We have previously identified multiple unique regions on Ch5 that repress the growth on sorbose. In one of the regions, the CSU51 gene determining the repressive property of the region was identified. We report here the identification of the CSU53 gene from a different region on Ch5. Most importantly, we find that CSU51 and CSU53 are associated with novel regulatory elements, ASUs, which are embedded within CSUs in an antisense configuration. ASUs act opposite to CSUs by enhancing the growth on sorbose. In respect to the CSU transcripts, the ASU long antisense transcripts are in lesser amounts, are completely overlapped, and are inversely related. ASUs interact with CSUs in natural CSU/ASU cis configurations, as well as when extra copies of ASUs are placed in trans to the CSU/ASU configurations. We suggest that ASU long embedded antisense transcripts modulate CSU sense transcripts.

Fidelity of tRNA 5'-maturation: a possible basis for the functional dependence of archaeal and eukaryal RNase P on multiple protein cofactors

RNase P, which catalyzes tRNA 5′-maturation, typically comprises a catalytic RNase P RNA (RPR) and a varying number of RNase P proteins (RPPs): 1 in bacteria, at least 4 in archaea and 9 in eukarya. The four archaeal RPPs have eukaryotic homologs and function as heterodimers (POP5•RPP30 and RPP21•RPP29). By studying the archaeal Methanocaldococcus jannaschii RPR's cis cleavage of precursor tRNA^Gln (pre-tRNA^Gln), which lacks certain consensus structures/sequences needed for substrate recognition, we demonstrate that RPP21•RPP29 and POP5•RPP30 can rescue the RPR's mis-cleavage tendency independently by 4-fold and together by 25-fold, suggesting that they operate by distinct mechanisms. This synergistic and preferential shift toward correct cleavage results from the ability of archaeal RPPs to selectively increase the RPR's apparent rate of correct cleavage by 11 140-fold, compared to only 480-fold for mis-cleavage. Moreover, POP5•RPP30, like the bacterial RPP, helps normalize the RPR's rates of cleavage of non-consensus and consensus pre-tRNAs. We also show that archaeal and eukaryal RNase P, compared to their bacterial relatives, exhibit higher fidelity of 5′-maturation of pre-tRNA^Gln and some of its mutant derivatives. Our results suggest that protein-rich RNase P variants might have evolved to support flexibility in substrate recognition while catalyzing efficient, high-fidelity 5′-processing.

Enrichment of unstable non-coding RNAs and their genome-wide identification

Cryptic unstable transcripts (CUTs) have been recently described as a major class of non-coding RNAs. These transcripts are, however, extremely unstable in normal cells and their analyzes pose specific technical problems. In this chapter, after a brief introduction discussing general aspects associated with the analysis of non-coding RNAs, we provide details of methods to enrich, map, and quantify this unconventional class of transcripts.

Adaptation of organisms by resonance of RNA transcription with the cellular redox cycle

Sequence variation in organisms differs across the genome and the majority of mutations are caused by oxidation, yet its origin is not fully understood. It has also been shown that the reduction-oxidation reaction cycle is the fundamental biochemical cycle that coordinates the timing of all biochemical processes in the cell, including energy production, DNA replication, and RNA transcription. We show that the temporal resonance of transcriptome biosynthesis with the oscillating binary state of the reduction-oxidation reaction cycle serves as a basis for non-random sequence variation at specific genome-wide coordinates that change faster than by accumulation of chance mutations. This work demonstrates evidence for a universal, persistent and iterative feedback mechanism between the environment and heredity, whereby acquired variation between cell divisions can outweigh inherited variation.

tRNA biology charges to the front

tRNA biology has come of age, revealing an unprecedented level of understanding and many unexpected discoveries along the way. This review highlights new findings on the diverse pathways of tRNA maturation, and on the formation and function of a number of modifications. Topics of special focus include the regulation of tRNA biosynthesis, quality control tRNA turnover mechanisms, widespread tRNA cleavage pathways activated in response to stress and other growth conditions, emerging evidence of signaling pathways involving tRNA and cleavage fragments, and the sophisticated intracellular tRNA trafficking that occurs during and after biosynthesis.

Of proteins and RNA: the RNase P/MRP family

Nuclear ribonuclease (RNase) P is a ubiquitous essential ribonucleoprotein complex, one of only two known RNA-based enzymes found in all three domains of life. The RNA component is the catalytic moiety of RNases P across all phylogenetic domains; it contains a well-conserved core, whereas peripheral structural elements are diverse. RNA components of eukaryotic RNases P tend to be less complex than their bacterial counterparts, a simplification that is accompanied by a dramatic reduction of their catalytic ability in the absence of protein. The size and complexity of the protein moieties increase dramatically from bacterial to archaeal to eukaryotic enzymes, apparently reflecting the delegation of some structural functions from RNA to proteins and, perhaps, in response to the increased complexity of the cellular environment in the more evolutionarily advanced organisms; the reasons for the increased dependence on proteins are not clear. We review current information on RNase P and the closely related universal eukaryotic enzyme RNase MRP, focusing on their functions and structural organization.

Chromatin remodeling around nucleosome-free regions leads to repression of noncoding RNA transcription

Nucleosome-free regions (NFRs) at the 5' and 3' ends of genes are general sites of transcription initiation for mRNA and noncoding RNA (ncRNA). The presence of NFRs within transcriptional regulatory regions and the conserved location of transcription start sites at NFRs strongly suggest that the regulation of NFRs profoundly affects transcription initiation. To date, multiple factors are known to facilitate transcription initiation by positively regulating the formation and/or size of NFRs in vivo. However, mechanisms to repress transcription by negatively regulating the size of NFRs have not been identified. We identified four distinct classes of NFRs located at the 5' and 3' ends of genes, within open reading frames (ORFs), and far from ORFs. The ATP-dependent chromatin-remodeling enzyme Isw2 was found enriched at all classes of NFRs. Analysis of RNA levels also demonstrated Isw2 is required to repress ncRNA transcription from many of these NFRs. Thus, by the systematic annotation of NFRs across the yeast genome and analysis of ncRNA transcription, we established, for the first time, a mechanism by which NFR size is negatively regulated to repress ncRNA transcription from NFRs. Finally, we provide evidence suggesting that one biological consequence of repression of ncRNA, by Isw2 or by the exosome, is prevention of transcriptional interference of mRNA.

Detecting novel genes with sparse arrays

Species-specific genes play an important role in defining the phenotype of an organism. However, current gene prediction methods can only efficiently find genes that share features such as sequence similarity or general sequence characteristics with previously known genes. Novel sequencing methods and tiling arrays can be used to find genes without prior information and they have demonstrated that novel genes can still be found from extensively studied model organisms. Unfortunately, these methods are expensive and thus are not easily applicable, e.g., to finding genes that are expressed only in very specific conditions. We demonstrate a method for finding novel genes with sparse arrays, applying it on the 33.9 Mb genome of the filamentous fungus Trichoderma reesei. Our computational method does not require normalisations between arrays and it takes into account the multiple-testing problem typical for analysis of microarray data. In contrast to tiling arrays, that use overlapping probes, only one 25 mer microarray oligonucleotide probe was used for every 100b. Thus, only relatively little space on a microarray slide was required to cover the intergenic regions of a genome. The analysis was done as a by-product of a conventional microarray experiment with no additional costs. We found at least 23 good candidates for novel transcripts that could code for proteins and all of which were expressed at high levels. Candidate genes were found to neighbour ire1 and cre1 and many other regulatory genes. Our simple, low-cost method can easily be applied to finding novel species-specific genes without prior knowledge of their sequence properties.

High-resolution transcription atlas of the mitotic cell cycle in budding yeast

BACKGROUND

Extensive transcription of non-coding RNAs has been detected in eukaryotic genomes and is thought to constitute an additional layer in the regulation of gene expression. Despite this role, their transcription through the cell cycle has not been studied; genome-wide approaches have only focused on protein-coding genes. To explore the complex transcriptome architecture underlying the budding yeast cell cycle, we used 8 bp tiling arrays to generate a 5 minute-resolution, strand-specific expression atlas of the whole genome.

RESULTS

We discovered 523 antisense transcripts, of which 80 cycle or are located opposite periodically expressed mRNAs, 135 unannotated intergenic non-coding RNAs, of which 11 cycle, and 109 cell-cycle-regulated protein-coding genes that had not previously been shown to cycle. We detected periodic expression coupling of sense and antisense transcript pairs, including antisense transcripts opposite of key cell-cycle regulators, like FAR1 and TAF2.

CONCLUSIONS

Our dataset presents the most comprehensive resource to date on gene expression during the budding yeast cell cycle. It reveals periodic expression of both protein-coding and non-coding RNA and profiles the expression of non-annotated RNAs throughout the cell cycle for the first time. This data enables hypothesis-driven mechanistic studies concerning the functions of non-coding RNAs.

Life without RNAi: noncoding RNAs and their functions in Saccharomyces cerevisiae

There are a number of well-characterized and fundamental roles for noncoding RNAs (ncRNAs) in gene regulation in all kingdoms of life. ncRNAs, such as ribosomal RNAs, transfer RNAs, small nuclear RNAs, small nucleolar RNAs, and small interfering RNAs, can serve catalytic and scaffolding functions in transcription, messenger RNA processing, translation, and RNA degradation. Recently, our understanding of gene expression has been dramatically challenged by the identification of large and diverse populations of novel ncRNAs in the eukaryotic genomes surveyed thus far. Studies carried out using the budding yeast Saccharomyces cerevisiae indicate that at least some coding genes are regulated by these novel ncRNAs. S. cerevisiae lacks RNA interference (RNAi) and, thus, provides an ideal system for studying the RNAi-independent mechanisms of ncRNA-based gene regulation. The current picture of gene regulation is one of great unknowns, in which the transcriptional environment surrounding a given locus may have as much to do with its regulation as its DNA sequence or local chromatin structure. Drawing on the recent research in S. cerevisiae and other organisms, this review will discuss the identification of ncRNAs, their origins and processing, and several models that incorporate ncRNAs into the regulation of gene expression and chromatin structure.

Unexpected diversity of RNase P, an ancient tRNA processing enzyme: challenges and prospects

For an enzyme functioning predominantly in a seemingly housekeeping role of 5' tRNA maturation, RNase P displays a remarkable diversity in subunit make-up across the three domains of life. Despite the protein complexity of this ribonucleoprotein enzyme increasing dramatically from bacteria to eukarya, the catalytic function rests with the RNA subunit during evolution. However, the recent demonstration of a protein-only human mitochondrial RNase P has added further intrigue to the compositional variability of this enzyme. In this review, we discuss some possible reasons underlying the structural diversity of the active sites, and use them as thematic bases for elaborating new directions to understand how functional variations might have contributed to the complex evolution of RNase P.

Broadening the mission of an RNA enzyme

The "RNA World" hypothesis suggests that life developed from RNA enzymes termed ribozymes, which carry out reactions without assistance from proteins. Ribonuclease (RNase) P is one ribozyme that appears to have adapted these origins to modern cellular life by adding protein to the RNA core in order to broaden the potential functions. This RNA-protein complex plays diverse roles in processing RNA, but its best-understood reaction is pre-tRNA maturation, resulting in mature 5' ends of tRNAs. The core catalytic activity resides in the RNA subunit of almost all RNase P enzymes but broader substrate tolerance is required for recognizing not only the diverse sequences of tRNAs, but also additional cellular RNA substrates. This broader substrate tolerance is provided by the addition of protein to the RNA core and allows RNase P to selectively recognize different RNAs, and possibly ribonucleoprotein (RNP) substrates. Thus, increased protein content correlated with evolution from bacteria to eukaryotes has further enhanced substrate potential enabling the enzyme to function in a complex cellular environment.

Toggle involving cis-interfering noncoding RNAs controls variegated gene expression in yeast

The identification of specific functional roles for the numerous long noncoding (nc)RNAs found in eukaryotic transcriptomes is currently a matter of intense study amid speculation that these ncRNAs have key regulatory roles. We have identified a pair of cis-interfering ncRNAs in yeast that contribute to the control of variegated gene expression at the FLO11 locus by implementing a regulatory circuit that toggles between two stable states. These capped, polyadenylated ncRNAs are transcribed across the large intergenic region upstream of the FLO11 ORF. As with mammalian long intervening (li)ncRNAs, these yeast ncRNAs (ICR1 and PWR1) are themselves regulated by transcription factors (Sfl1 and Flo8) and chromatin remodelers (Rpd3L) that are key elements in phenotypic transitions in yeast. The mechanism that we describe explains the unanticipated role of a histone deacetylase complex in activating gene expression, because Rpd3L mutants force the ncRNA circuit into a state that silences the expression of the adjacent variegating gene.

Genome reprogramming during sporulation

When environmental conditions compromise survival, single celled organisms, such as the budding yeast S. cerevisiae, induce and complete a differentiation program called sporulation. The first step consists of meiosis, which generates genetic diversity within the eventual haploid cells. The post-meiotic maturation stage reinforces protective barriers, such as the spore wall, against deleterious external conditions. In later stages of sporulation, the spore nucleus becomes highly compacted, likely sharing certain characteristics with the metazoan male gamete, the spermatozoon. The sporulation differentiation program involves many chromatin-related events, including execution of a precise transcription program involving more than one thousand genes. Here, we review how chromatin structure and genome reprogramming regulate the sporulation transcription program, and how post-meiotic events reorganize spore chromatin.

Unravelling cis-regulatory elements in the genome of the smallest photosynthetic eukaryote: phylogenetic footprinting in Ostreococcus

We used a phylogenetic footprinting approach, adapted to high levels of divergence, to estimate the level of constraint in intergenic regions of the extremely gene dense Ostreococcus algae genomes (Chlorophyta, Prasinophyceae). We first benchmarked our method against the Saccharomyces sensu stricto genome data and found that the proportion of conserved non-coding sites was consistent with those obtained with methods using calibration by the neutral substitution rate. We then applied our method to the complete genomes of Ostreococcus tauri and O. lucimarinus, which are the most divergent species from the same genus sequenced so far. We found that 77% of intergenic regions in Ostreococcus still contain some phylogenetic footprints, as compared to 88% for Saccharomyces, corresponding to an average rate of constraint on intergenic region of 17% and 30%, respectively. A comparison with some known functional cis-regulatory elements enabled us to investigate whether some transcriptional regulatory pathways were conserved throughout the green lineage. Strikingly, the size of the phylogenetic footprints depends on gene orientation of neighboring genes, and appears to be genus-specific. In Ostreococcus, 5' intergenic regions contain four times more conserved sites than 3' intergenic regions, whereas in yeast a higher frequency of constrained sites in intergenic regions between genes on the same DNA strand suggests a higher frequency of bidirectional regulatory elements. The phylogenetic footprinting approach can be used despite high levels of divergence in the ultrasmall Ostreococcus algae, to decipher structure of constrained regulatory motifs, and identify putative regulatory pathways conserved within the green lineage.

Whole-genome tiling array analysis of Mycobacterium leprae RNA reveals high expression of pseudogenes and noncoding regions

Whole-genome sequence analysis of Mycobacterium leprae has revealed a limited number of protein-coding genes, with half of the genome composed of pseudogenes and noncoding regions. We previously showed that some M. leprae pseudogenes are transcribed at high levels and that their expression levels change following infection. In order to clarify the RNA expression profile of the M. leprae genome, a tiling array in which overlapping 60-mer probes cover the entire 3.3-Mbp genome was designed. The array was hybridized with M. leprae RNA from the SHR/NCrj-rnu nude rat, and the results were compared to results from an open reading frame array and confirmed by reverse transcription-PCR. RNA expression was detected from genes, pseudogenes, and noncoding regions. The signal intensities obtained from noncoding regions were higher than those from pseudogenes. Expressed noncoding regions include the M. leprae unique repetitive sequence RLEP and other sequences without any homology to known functional noncoding RNAs. Although the biological functions of RNA transcribed from M. leprae pseudogenes and noncoding regions are not known, RNA expression analysis will provide insights into the bacteriological significance of the species. In addition, our study suggests that M. leprae will be a useful model organism for the study of the molecular mechanism underlying the creation of pseudogenes and the role of microRNAs derived from noncoding regions.

Chromatin- and transcription-related factors repress transcription from within coding regions throughout the Saccharomyces cerevisiae genome

Previous studies in Saccharomyces cerevisiae have demonstrated that cryptic promoters within coding regions activate transcription in particular mutants. We have performed a comprehensive analysis of cryptic transcription in order to identify factors that normally repress cryptic promoters, to determine the amount of cryptic transcription genome-wide, and to study the potential for expression of genetic information by cryptic transcription. Our results show that a large number of factors that control chromatin structure and transcription are required to repress cryptic transcription from at least 1,000 locations across the S. cerevisiae genome. Two results suggest that some cryptic transcripts are translated. First, as expected, many cryptic transcripts contain an ATG and an open reading frame of at least 100 codons. Second, several cryptic transcripts are translated into proteins. Furthermore, a subset of cryptic transcripts tested is transiently induced in wild-type cells following a nutritional shift, suggesting a possible physiological role in response to a change in growth conditions. Taken together, our results demonstrate that, during normal growth, the global integrity of gene expression is maintained by a wide range of factors and suggest that, under altered genetic or physiological conditions, the expression of alternative genetic information may occur.

Non-coding RNA prediction and verification in Saccharomyces cerevisiae

Non-coding RNA (ncRNA) play an important and varied role in cellular function. A significant amount of research has been devoted to computational prediction of these genes from genomic sequence, but the ability to do so has remained elusive due to a lack of apparent genomic features. In this work, thermodynamic stability of ncRNA structural elements, as summarized in a Z-score, is used to predict ncRNA in the yeast Saccharomyces cerevisiae. This analysis was coupled with comparative genomics to search for ncRNA genes on chromosome six of S. cerevisiae and S. bayanus. Sets of positive and negative control genes were evaluated to determine the efficacy of thermodynamic stability for discriminating ncRNA from background sequence. The effect of window sizes and step sizes on the sensitivity of ncRNA identification was also explored. Non-coding RNA gene candidates, common to both S. cerevisiae and S. bayanus, were verified using northern blot analysis, rapid amplification of cDNA ends (RACE), and publicly available cDNA library data. Four ncRNA transcripts are well supported by experimental data (RUF10, RUF11, RUF12, RUF13), while one additional putative ncRNA transcript is well supported but the data are not entirely conclusive. Six candidates appear to be structural elements in 5′ or 3′ untranslated regions of annotated protein-coding genes. This work shows that thermodynamic stability, coupled with comparative genomics, can be used to predict ncRNA with significant structural elements.

Recent advances in DNA sequence technology have made it possible to sequence entire genomes. Once a genome is sequenced, it becomes necessary to identify the set of genes and other functional elements within the genome. This is particularly challenging as much of the genomic sequence does not appear to perform any function and is loosely referred to as “junk.” Identifying functional elements among the “junk” is difficult. Experimental methods have been developed for this purpose but they are time-consuming, expensive, and often provide an incomplete picture. Thus, it is important to develop the ability to identify these functional elements using computational methods. Protein-coding genes are relatively easy to identify computationally, but other categories of functional elements present a significantly greater challenge. In this work, we used a computational approach to identify genes that do not encode for a protein but rather function as an RNA molecule. We then used experimental methods to verify our predictions and thereby validate the computational method.

Whole-genome microarrays: applications and technical issues

DNA microarrays have become a mainstream tool in experimental plant biology. The constant improvements in the technological platforms have enabled the development of the tiling DNA microarrays that cover the whole genome, which in turn catalyzed the wide variety of creative applications of such microarrays in the areas as diverse as global studies of genetic variation, DNA-binding proteins, DNA methylation, and chromatin and transcriptome dynamics. This chapter attempts to summarize such applications as well as discusses some technical and strategic issues that are particular to the use of tiling microarrays.

RNase P: increased versatility through protein complexity?

Ribonuclease P (RNase P) is an essential enzyme that catalyzes the 5' endonucleolytic cleavage of precursor transfer RNAs (pretRNAs). It is found in all phylogenetic domains: bacteria, archaea and eukaryotes. The bacterial enzyme consists of a single, catalytic RNA subunit and one small protein, while the archaeal and eukaryotic enzymes have 4-10 proteins in addition to a similar RNA subunit. The bacterial RNA acts as a ribozyme at high salt in vitro; however the added protein optimizes kinetics and makes specific contacts with the pre-tRNA substrate. The bacterial protein subunit also appears to be required for the processing of non-tRNA substrates by broadening recognition tolerance. In addition, the immense increase in protein content in the eukaryotic enzymes suggests substantially enlarged capacity for recognition of additional substrates. Recently intron-encoded box C/D snoRNAs were shown to be likely substrates for RNase P, with several lines of evidence suggesting that the nuclear holoenzyme binds tightly to, and can cleave single-stranded RNA in a sequence dependent fashion. The possible involvement of RNase P in additional RNA processing or turnover pathways would be consistent with previous findings that RNase MRP, a variant of RNase P that has evolved to participate in ribosomal RNA processing, is also involved in turnover of specific messenger RNAs. Here, involvement of RNase P in multiple RNA processing pathways is discussed.

An efficient genetic algorithm for structural RNA pairwise alignment and its application to non-coding RNA discovery in yeast

Background

Aligning RNA sequences with low sequence identity has been a challenging problem since such a computation essentially needs an algorithm with high complexities for taking structural conservation into account. Although many sophisticated algorithms for the purpose have been proposed to date, further improvement in efficiency is necessary to accelerate its large-scale applications including non-coding RNA (ncRNA) discovery.

Results

We developed a new genetic algorithm, Cofolga2, for simultaneously computing pairwise RNA sequence alignment and consensus folding, and benchmarked it using BRAliBase 2.1. The benchmark results showed that our new algorithm is accurate and efficient in both time and memory usage. Then, combining with the originally trained SVM, we applied the new algorithm to novel ncRNA discovery where we compared S. cerevisiae genome with six related genomes in a pairwise manner. By focusing our search to the relatively short regions (50 bp to 2,000 bp) sandwiched by conserved sequences, we successfully predict 714 intergenic and 1,311 sense or antisense ncRNA candidates, which were found in the pairwise alignments with stable consensus secondary structure and low sequence identity (≤ 50%). By comparing with the previous predictions, we found that > 92% of the candidates is novel candidates. The estimated rate of false positives in the predicted candidates is 51%. Twenty-five percent of the intergenic candidates has supports for expression in cell, i.e. their genomic positions overlap those of the experimentally determined transcripts in literature. By manual inspection of the results, moreover, we obtained four multiple alignments with low sequence identity which reveal consensus structures shared by three species/sequences.

Conclusion

The present method gives an efficient tool complementary to sequence-alignment-based ncRNA finders.

Genome-wide search for yeast RNase P substrates reveals role in maturation of intron-encoded box C/D small nucleolar RNAs

Ribonuclease P (RNase P) is an essential endonuclease responsible for the 5'-end maturation of precursor tRNAs. Bacterial RNase P also processes precursor 4.5S RNA, tmRNA, 30S preribosomal RNA, and several reported protein-coding RNAs. Eukaryotic nuclear RNase P is far more complex than in the bacterial form, employing multiple essential protein subunits in addition to the catalytic RNA subunit. RNomic studies have shown that RNase P binds other RNAs in addition to tRNAs, but no non-tRNA substrates have previously been identified. Additional substrates were identified by using a multipronged approach in the budding yeast Saccharomyces cerevisiae. First, RNase P-dependant changes in RNA abundance were examined on whole-genome microarrays by using strains containing temperature sensitive (TS) mutations in two of the essential RNase P subunits, Pop1p and Rpr1r. Second, RNase P was rapidly affinity-purified, and copurified RNAs were identified by using a genome-wide microarray. Third, to identify RNAs that do not change abundance when RNase P is depleted but accumulate as larger precursors, >80 potential small RNA substrates were probed directly by Northern blot analysis with RNA from the RNase P TS mutants. Numerous potential substrates were identified, of which we characterized the box C/D intron-encoded small nucleolar RNAs (snoRNAs), because these both copurify with RNase P and accumulate larger forms in the RNase P temperature-sensitive mutants. It was previously known that two pathways existed for excising these snoRNAs, one using the pre-mRNA splicing path and the other that was independent of splicing. RNase P appears to participate in the splicing-independent path for the box C/D intron-encoded snoRNAs.

The Nrd1-Nab3-Sen1 termination complex interacts with the Ser5-phosphorylated RNA polymerase II C-terminal domain

RNA polymerase II (Pol II) in Saccharomyces cerevisiae can terminate transcription via several pathways. To study how a mechanism is chosen, we analyzed recruitment of Nrd1, which cooperates with Nab3 and Sen1 to terminate small nucleolar RNAs and other short RNAs. Budding yeast contains three C-terminal domain (CTD) interaction domain (CID) proteins, which bind the CTD of the Pol II largest subunit. Rtt103 and Pcf11 act in mRNA termination, and both preferentially interact with CTD phosphorylated at Ser2. The crystal structure of the Nrd1 CID shows a fold similar to that of Pcf11, but Nrd1 preferentially binds to CTD phosphorylated at Ser5, the form found proximal to promoters. This indicates why Nrd1 cross-links near 5' ends of genes and why the Nrd1-Nab3-Sen1 termination pathway acts specifically at short Pol II-transcribed genes. Nrd1 recruitment to genes involves a combination of interactions with CTD and Nab3.

At-TAX: a whole genome tiling array resource for developmental expression analysis and transcript identification in Arabidopsis thaliana

Gene expression maps for model organisms, including Arabidopsis thaliana, have typically been created using gene-centric expression arrays. Here, we describe a comprehensive expression atlas, Arabidopsis thaliana Tiling Array Express (At-TAX), which is based on whole-genome tiling arrays. We demonstrate that tiling arrays are accurate tools for gene expression analysis and identified more than 1,000 unannotated transcribed regions. Visualizations of gene expression estimates, transcribed regions, and tiling probe measurements are accessible online at the At-TAX homepage.

Whole-genome expression profiling of the marine diatom Thalassiosira pseudonana identifies genes involved in silicon bioprocesses

Formation of complex inorganic structures is widespread in nature. Diatoms create intricately patterned cell walls of inorganic silicon that are a biomimetic model for design and generation of three-dimensional silica nanostructures. To date, only relatively simple silica structures can be generated in vitro through manipulation of known diatom phosphoproteins (silaffins) and long-chain polyamines. Here, we report the use of genome-wide transcriptome analyses of the marine diatom Thalassiosira pseudonana to identify additional candidate gene products involved in the biological manipulation of silicon. Whole-genome oligonucleotide tiling arrays and tandem mass spectrometry identified transcripts for >8,000 genes, approximately 3,000 of which were not previously described and included noncoding and antisense RNAs. Gene-specific expression profiles detected a set of 75 genes induced only under low concentrations of silicon but not under low concentrations of nitrogen or iron, alkaline pH, or low temperatures. Most of these induced gene products were predicted to contain secretory signals and/or transmembrane domains but displayed no homology to known proteins. Over half of these genes were newly discovered, identified only through the use of tiling arrays. Unexpectedly, a common set of 84 genes were induced by both silicon and iron limitations, suggesting that biological manipulation of silicon may share pathways in common with iron or, alternatively, that iron may serve as a required cofactor for silicon processes. These results provide insights into the transcriptional and translational basis for the biological generation of elaborate silicon nanostructures by these ecologically important microbes.

The current status of cDNA cloning

The cloning of cDNAs, copies of cellular RNA, is one of the classical technologies in molecular biology. Over the past 30 years cDNA cloning technologies have been improved to enable the cloning of large cDNA collections, which are fundamental to today's understanding of the utilization of genetic information. With the discovery of noncoding RNAs, additional new approaches to the cloning of short RNAs have been developed. However, with the realization that much larger portions of genomes are transcribed than anticipated from genome annotations, cDNA cloning faces new challenges to uncover rare transcripts and to make the corresponding cDNAs available for functional studies. This review provides an overview on the current status of cDNA cloning and possibilities for the discovery and characterization of new RNA families.

Transcriptional interference: an unexpected layer of complexity in gene regulation

Much of the genome is transcribed into long untranslated RNAs, mostly of unknown function. Growing evidence suggests that transcription of sense and antisense untranslated RNAs in eukaryotes can repress a neighboring gene by a phenomenon termed transcriptional interference. Transcriptional interference by the untranslated RNA may prevent recruitment of the initiation complex or prevent transcriptional elongation. Recent work in yeast, mammals, and Drosophila highlights the diverse roles that untranslated RNAs play in development. Previously, untranslated RNAs of the bithorax complex of Drosophila were proposed to be required for its activation. Recent studies show that these untranslated RNAs in fact silence Ultrabithorax in early embryos, probably by transcriptional interference.

Antisense artifacts in transcriptome microarray experiments are resolved by actinomycin D

Recent transcription profiling studies have revealed an unanticipatedly large proportion of antisense transcription across eukaryotic and bacterial genomes. However, the extent and significance of antisense transcripts is controversial partly because experimental artifacts are suspected. Here, we present a method to generate clean genome-wide transcriptome profiles, using actinomycin D (ActD) during reverse transcription. We show that antisense artifacts appear to be triggered by spurious synthesis of second-strand cDNA during reverse transcription reactions. Strand-specific hybridization signals obtained from Saccharomyces cerevisiae tiling arrays were compared between samples prepared with and without ActD. Use of ActD removed about half of the detectable antisense transcripts, consistent with their being artifacts, while sense expression levels and about 200 antisense transcripts were not affected. Our findings thus facilitate a more accurate assessment of the extent and position of antisense transcription, towards a better understanding of its role in cells.

A view of RNase P

Major progress in the study of RNase P has resulted from crystallography of bacterial catalytic subunits and the discovery of catalytic activity in eukaryotes. Several new substrates have also been identified, primarily in bacteria but also in yeast. Our current world should be called the "RNA-protein world" rather than the "protein world".

Identification of an in vitro transcription-based artifact affecting oligonucleotide microarrays

This study identified the widely used T7 in vitro transcription system as a major source of artifact in the tiling array data from nine eukaryotic genomes. The most affected probes contained a sequence motif complementary to the +1 to +9 initial transcribed sequence (ITS) of the T7-(dT)(24) primer. The abundance of 5' ITS cRNA fragments produced during target preparation was sufficient to drive undesirable hybridization. A new T7-(dT)(24) primer with a modified ITS was designed that shifts the artifactual motifs as predicted and reduces the effect of the artifact. A computational algorithm was generated to filter out the likely artifactual probes from existing whole-genome tiling array data and improve probe selection. Further studies of Arabidopsis thaliana were conducted using both T7-(dT)(24) primers. While the artifact affected transcript discovery with tiling arrays, it showed only a minor impact on measurements of gene expression using commercially available 'gene-only' expression arrays.

Turnover and function of noncoding RNA polymerase II transcripts

In the past few years, especially since the discovery of RNA interference (RNAi), our understanding of the role of RNA in gene expression has undergone a significant transformation. This change has been brought about by growing evidence that RNA is more complex and transcription more promiscuous than has previously been thought. Many of the new transcripts are of so-called noncoding RNA (ncRNA); i.e., RNA that does not code for proteins such as mRNA, or intrinsic parts of the cellular machinery such as the highly structured RNA components of ribosomes (rRNA) and the small nuclear RNA (snRNA) components of the splicing machinery. It is becoming increasingly apparent that ncRNAs have very important roles in gene expression. This paper focuses on work from our laboratory in which we have investigated the roles and turnover of ncRNA located within the gene pre-mRNA, which we refer to as intragenic ncRNA. Also discussed are some investigations of intergenic ncRNA transcription and how these two classes of ncRNA may interrelate.

Comparative analysis of structured RNAs in S. cerevisiae indicates a multitude of different functions

Background

Non-coding RNAs (ncRNAs) are an emerging focus for both computational analysis and experimental research, resulting in a growing number of novel, non-protein coding transcripts with often unknown functions. Whole genome screens in higher eukaryotes, for example, provided evidence for a surprisingly large number of ncRNAs. To supplement these searches, we performed a computational analysis of seven yeast species and searched for new ncRNAs and RNA motifs.

Results

A comparative analysis of the genomes of seven yeast species yielded roughly 2800 genomic loci that showed the hallmarks of evolutionary conserved RNA secondary structures. A total of 74% of these regions overlapped with annotated non-coding or coding genes in yeast. Coding sequences that carry predicted structured RNA elements belong to a limited number of groups with common functions, suggesting that these RNA elements are involved in post-transcriptional regulation and/or cellular localization. About 700 conserved RNA structures were found outside annotated coding sequences and known ncRNA genes. Many of these predicted elements overlapped with UTR regions of particular classes of protein coding genes. In addition, a number of RNA elements overlapped with previously characterized antisense transcripts. Transcription of about 120 predicted elements located in promoter regions and other, previously un-annotated, intergenic regions was supported by tiling array experiments, ESTs, or SAGE data.

Conclusion

Our computational predictions strongly suggest that yeasts harbor a substantial pool of several hundred novel ncRNAs. In addition, we describe a large number of RNA structures in coding sequences and also within antisense transcripts that were previously characterized using tiling arrays.