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

Detecting transcriptionally active regions using genomic tiling arrays

We have developed a method for interpreting genomic tiling array data, implemented as the program TranscriptionDetector. Probed loci expressed above background are identified by combining replicates in a way that makes minimal assumptions about the data. We performed medium-resolution Anopheles gambiae tiling array experiments and found extensive transcription of both coding and non-coding regions. Our method also showed improved detection of transcriptional units when applied to high-density tiling array data for ten human chromosomes.

RNase MRP RNA and human genetic diseases

RNase MRP RNA is the RNA subunit of the RNase mitochondrial RNA processing (MRP) enzyme complex that is involved in multiple cellular RNA processing events. Mutations on RNase MRP RNA gene (RMRP) cause a recessively inherited developmental disorder, cartilage-hair hypoplasia (CHH). The relationship of the genotype (RMRP mutation), RNA processing deficiency of the RNase MRP complex, and the phenotype of CHH and other skeletal dysplasias is yet to be explored.

A noncoding RNA in Saccharomyces cerevisiae is an RNase P substrate

Ribonuclease P (RNase P) is involved in regulation of noncoding RNA (ncRNA) expression in Saccharomyces cerevisiae. A hidden-in-reading-frame antisense-1 (HRA1) RNA in S. cerevisiae, which belongs to a class of ncRNAs located in the antisense strand to verified protein coding regions, was cloned for further use in RNase P assays. Escherichia coli RNase P assays in vitro of HRA1 RNA show two cleavage sites, one major and one minor in terms of rates. The same result was observed with a partially purified S. cerevisiae RNase P activity, both at 30 degrees C and 37 degrees C. These latter cells are normally grown at 30 degrees C. Predictions of the secondary structure of HRA1 RNA in silico show the cleavage sites are canonical RNase P recognition sites. A relatively small amount of endogenous HRA1 RNA was identified by RT-PCR in yeast cells. The endogenous HRA1 RNA is increased in amount in strains that are deficient in RNase P activity. A deletion of 10 nucleotides in the HRA1 gene that does not overlap with the gene coding for a protein (DRS2) in the sense strand shows no defective growth in galactose or glucose. These data indicate that HRA1 RNA is a substrate for RNase P and does not appear as a direct consequence of separate regulatory effects of the enzyme on ncRNAs.

A role for noncoding transcription in activation of the yeast PHO5 gene

Noncoding, or intergenic, transcription by RNA polymerase II (RNAPII) is remarkably widespread in eukaryotic organisms, but the effects of such transcription remain poorly understood. Here we show that noncoding transcription plays a role in activation, but not repression, of the Saccharomyces cerevisiae PHO5 gene. Histone eviction from the PHO5 promoter during activation occurs with normal kinetics even in the absence of the PHO5 TATA box, showing that transcription of the gene itself is not required for promoter remodeling. Nevertheless, we find that mutations that impair transcript elongation by RNAPII affect the kinetics of histone eviction from the PHO5 promoter. Most dramatically, inactivation of RNAPII itself abolishes eviction completely. Under repressing conditions, an approximately 2.4-kb noncoding exosome-degraded transcript is detected that originates near the PHO5 termination site and is transcribed in the antisense direction. Abrogation of this transcript delays chromatin remodeling and subsequent RNAPII recruitment to PHO5 upon activation. We propose that noncoding transcription through positioned nucleosomes can enhance chromatin plasticity so that chromatin remodeling and activation of traversed genes occur in a timely manner.

Assessing the need for sequence-based normalization in tiling microarray experiments

MOTIVATION

Increases in microarray feature density allow the construction of so-called tiling microarrays. These arrays, or sets of arrays, contain probes targeting regions of sequenced genomes at regular genomic intervals. The unbiased nature of this approach allows for the identification of novel transcribed sequences, the localization of transcription factor binding sites (ChIP-chip), and high resolution comparative genomic hybridization, among other uses. These applications are quickly growing in popularity as tiling microarrays become more affordable. To reach maximum utility, the tiling microarray platform needs be developed to the point that 1 nt resolutions are achieved and that we have confidence in individual measurements taken at this fine of resolution. Any biases in tiling array signals must be systematically removed to achieve this goal.

RESULTS

Towards this end, we investigated the importance of probe sequence composition on the efficacy of tiling microarrays for identifying novel transcription and transcription factor binding sites. We found that intensities are highly sequence dependent and can greatly influence results. We developed three metrics for assessing this sequence dependence and use them in evaluating existing sequence-based normalizations from the tiling microarray literature. In addition, we applied three new techniques for addressing this problem; one method, adapted from similar work on GeneChip brand microarrays, is based on modeling array signal as a linear function of probe sequence, the second method extends this approach by iterative weighting and re-fitting of the model, and the third technique extrapolates the popular quantile normalization algorithm for between-array normalization to probe sequence space. These three methods perform favorably to existing strategies, based on the metrics defined here.

AVAILABILITY

http://tiling.gersteinlab.org/sequence_effects/

Antisense transcription controls cell fate in Saccharomyces cerevisiae

Entry into meiosis is a key developmental decision. We show here that meiotic entry in Saccharomyces cerevisiae is controlled by antisense-mediated regulation of IME4, a gene required for initiating meiosis. In MAT a/alpha diploids the antisense IME4 transcript is repressed by binding of the a1/alpha2 heterodimer at a conserved site located downstream of the IME4 coding sequence. MAT a/alpha diploids that produce IME4 antisense transcript have diminished sense transcription and fail to initiate meiosis. Haploids that produce the sense transcript have diminished antisense transcription and manifest several diploid phenotypes. Our data are consistent with transcription interference as a regulatory mechanism at the IME4 locus that determines cell fate.

The transcription factor Gcr1 stimulates cell growth by participating in nutrient-responsive gene expression on a global level

Transcriptomic reprogramming is critical to the coordination between growth and cell cycle progression in response to changing extracellular conditions. In Saccharomyces cerevisiae, the transcription factor Gcr1 contributes to this coordination by supporting maximum expression of G1 cyclins in addition to regulating both glucose-induced and glucose-repressed genes. We report here the comprehensive genome-wide expression profiling of gcr1Delta cells. Our data show that reduced expression of ribosomal protein genes in gcr1Delta cells is detectable both 20 min after glucose addition and in steady-state cultures of raffinose-grown cells, showing that this defect is not the result of slow growth or growth on a repressing sugar. However, the large cell phenotype of the gcr1Delta mutant occurs only in the presence of repressing sugars. GCR1 deletion also results in aberrant derepression of numerous glucose repressed loci; glucose-grown gcr1Delta cells actively respire, demonstrating that this global alteration in transcription corresponds to significant changes at the physiological level. These data offer an insight into the coordination of growth and cell division by providing an integrated view of the transcriptomic, phenotypic, and metabolic consequences of GCR1 deletion.

Termination of cryptic unstable transcripts is directed by yeast RNA-binding proteins Nrd1 and Nab3

Studies of yeast transcription have revealed the widespread distribution of intergenic RNA polymerase II transcripts. These cryptic unstable transcripts (CUTs) are rapidly degraded by the nuclear exosome. Yeast RNA binding proteins Nrd1 and Nab3 direct termination of sn/snoRNAs and recently have also been implicated in premature transcription termination of the NRD1 gene. In this paper, we show that Nrd1 and Nab3 are required for transcription termination of CUTs. In nrd1 and nab3 mutants, we observe 3'-extended transcripts originating from CUT promoters but failing to terminate through the Nrd1- and Nab3-directed pathway. Nrd1 and Nab3 colocalize to regions of the genome expressing antisense CUTs, and these transcripts require yeast nuclear exosome and TRAMP components for degradation. Dissection of a CUT terminator reveals a minimal element sufficient for Nrd1- and Nab3-directed termination. These results suggest that transcription termination of CUTs directed by Nrd1 and Nab3 is a prerequisite for rapid degradation by the nuclear exosome.

TUF love for "junk" DNA

The widespread occurrence of noncoding (nc) RNAs--unannotated eukaryotic transcripts with reduced protein coding potential--suggests that they are functionally important. Study of ncRNAs is increasing our understanding of the organization and regulation of genomes.

Transcript mapping with high-density oligonucleotide tiling arrays

MOTIVATION

High-density DNA tiling microarrays are a powerful tool for the characterization of complete transcriptomes. The two major analytical challenges are the segmentation of the hybridization signal along genomic coordinates to accurately determine transcript boundaries and the adjustment of the sequence-dependent response of the oligonucleotide probes to achieve quantitative comparability of the signal between different probes.

RESULTS

We describe a dynamic programming algorithm for finding a globally optimal fit of a piecewise constant expression profile along genomic coordinates. We developed a probe-specific background correction and scaling method that employs empirical probe response parameters determined from reference hybridizations with no need for paired mismatch probes. This combined analysis approach allows the accurate determination of dynamical changes in transcription architectures from hybridization data and will help to study the biological significance of complex transcriptional phenomena in eukaryotic genomes.

AVAILABILITY

R package tilingArray at http://www.bioconductor.org.