Minor introns are embedded molecular switches regulated by highly unstable U6atac snRNA.

Eukaryotes have two types of spliceosomes, comprised of either major (U1, U2, U4, U5, U6) or minor (U11, U12, U4atac, U6atac; <1%) snRNPs. The high conservation of minor introns, typically one amidst many major introns in several hundred genes, despite their poor splicing, has been a long-standing enigma. Here, we discovered that the low abundance minor spliceosome's catalytic snRNP, U6atac, is strikingly unstable (t½<2 hr). We show that U6atac level depends on both RNA polymerases II and III and can be rapidly increased by cell stress-activated kinase p38MAPK, which stabilizes it, enhancing mRNA expression of hundreds of minor intron-containing genes that are otherwise suppressed by limiting U6atac. Furthermore, p38MAPK-dependent U6atac modulation can control minor intron-containing tumor suppressor PTEN expression and cytokine production. We propose that minor introns are embedded molecular switches regulated by U6atac abundance, providing a novel post-transcriptional gene expression mechanism and a rationale for the minor spliceosome's evolutionary conservation. DOI:http://dx.doi.org/10.7554/eLife.00780.001.

Vaccinia and influenza A viruses select rather than adjust tRNAs to optimize translation.

Transfer RNAs (tRNAs) are central to protein synthesis and impact translational speed and fidelity by their abundance. Here we examine the extent to which viruses manipulate tRNA populations to favor translation of their own genes. We study two very different viruses: influenza A virus (IAV), a medium-sized (13 kB genome) RNA virus; and vaccinia virus (VV), a large (200 kB genome) DNA virus. We show that the total cellular tRNA population remains unchanged following viral infection, whereas the polysome-associated tRNA population changes dramatically in a virus-specific manner. The changes in polysome-associated tRNA levels reflect the codon usage of viral genes, suggesting the existence of local tRNA pools optimized for viral translation.

Genome-wide determination of a broad ESRP-regulated posttranscriptional network by high-throughput sequencing.

Tissue-specific alternative splicing is achieved through the coordinated assembly of RNA binding proteins at specific sites to enhance or silence splicing at nearby splice sites. We used high-throughput sequencing (RNA-Seq) to investigate the complete spectrum of alternative splicing events that are regulated by the epithelium-specific splicing regulatory proteins ESRP1 and ESRP2. We also combined this analysis with direct RNA sequencing (DRS) to reveal ESRP-mediated regulation of alternative polyadenylation. To define binding motifs that mediate direct regulation of splicing and polyadenylation by ESRP, SELEX-Seq analysis was performed, coupling traditional SELEX with high-throughput sequencing. Identification and scoring of high-affinity ESRP1 binding motifs within ESRP target genes allowed the generation of RNA maps that define the position-dependent activity of the ESRPs in regulating cassette exons and alternative 3' ends. These extensive analyses provide a comprehensive picture of the functions of the ESRPs in an epithelial posttranscriptional gene expression program.

MATS: a Bayesian framework for flexible detection of differential alternative splicing from RNA-Seq data.

Ultra-deep RNA sequencing has become a powerful approach for genome-wide analysis of pre-mRNA alternative splicing. We develop MATS (multivariate analysis of transcript splicing), a bayesian statistical framework for flexible hypothesis testing of differential alternative splicing patterns on RNA-Seq data. MATS uses a multivariate uniform prior to model the between-sample correlation in exon splicing patterns, and a Markov chain Monte Carlo (MCMC) method coupled with a simulation-based adaptive sampling procedure to calculate the P-value and false discovery rate (FDR) of differential alternative splicing. Importantly, the MATS approach is applicable to almost any type of null hypotheses of interest, providing the flexibility to identify differential alternative splicing events that match a given user-defined pattern. We evaluated the performance of MATS using simulated and real RNA-Seq data sets. In the RNA-Seq analysis of alternative splicing events regulated by the epithelial-specific splicing factor ESRP1, we obtained a high RT-PCR validation rate of 86% for differential exon skipping events with a MATS FDR of <10%. Additionally, over the full list of RT-PCR tested exons, the MATS FDR estimates matched well with the experimental validation rate. Our results demonstrate that MATS is an effective and flexible approach for detecting differential alternative splicing from RNA-Seq data.

An ESRP-regulated splicing programme is abrogated during the epithelial-mesenchymal transition.

Alternative splicing achieves coordinated changes in post-transcriptional gene expression programmes through the activities of diverse RNA-binding proteins. Epithelial splicing regulatory proteins 1 and 2 (ESRP1 and ESRP2) are cell-type-specific regulators of transcripts that switch splicing during the epithelial-mesenchymal transition (EMT). To define a comprehensive programme of alternative splicing that is regulated during the EMT, we identified an extensive ESRP-regulated splicing network of hundreds of alternative splicing events within numerous genes with functions in cell-cell adhesion, polarity, and migration. Loss of this global ESRP-regulated epithelial splicing programme induces the phenotypic changes in cell morphology that are observed during the EMT. Components of this splicing signature provide novel molecular markers that can be used to characterize the EMT. Bioinformatics and experimental approaches revealed a high-affinity ESRP-binding motif and a predictive RNA map that governs their activity. This work establishes the ESRPs as coordinators of a complex alternative splicing network that adds an important post-transcriptional layer to the changes in gene expression that underlie epithelial-mesenchymal transitions during development and disease.

Rapid-response splicing reporter screens identify differential regulators of constitutive and alternative splicing.

Bioactive compounds have been invaluable for dissecting the mechanisms, regulation, and functions of cellular processes. However, very few such reagents have been described for pre-mRNA splicing. To facilitate their systematic discovery, we developed a high-throughput cell-based assay that measures pre-mRNA splicing by utilizing a quantitative reporter system with advantageous features. The reporter, consisting of a destabilized, intron-containing luciferase expressed from a short-lived mRNA, allows rapid screens (<4 h), thereby obviating the potential toxicity of splicing inhibitors. We describe three inhibitors (out of >23,000 screened), all pharmacologically active: clotrimazole, flunarizine, and chlorhexidine. Interestingly, none was a general splicing inhibitor. Rather, each caused distinct splicing changes of numerous genes. We further discovered the target of action of chlorhexidine and show that it is a selective inhibitor of specific Cdc2-like kinases (Clks) that phosphorylate serine-arginine-rich (SR) protein splicing factors. Our findings reveal unexpected activities of clinically used drugs in splicing and uncover differential regulation of constitutively spliced introns.

Innate immune and chemically triggered oxidative stress modifies translational fidelity.

Translational fidelity, essential for protein and cell function, requires accurate transfer RNA (tRNA) aminoacylation. Purified aminoacyl-tRNA synthetases exhibit a fidelity of one error per 10,000 to 100,000 couplings. The accuracy of tRNA aminoacylation in vivo is uncertain, however, and might be considerably lower. Here we show that in mammalian cells, approximately 1% of methionine (Met) residues used in protein synthesis are aminoacylated to non-methionyl-tRNAs. Remarkably, Met-misacylation increases up to tenfold upon exposing cells to live or non-infectious viruses, toll-like receptor ligands or chemically induced oxidative stress. Met is misacylated to specific non-methionyl-tRNA families, and these Met-misacylated tRNAs are used in translation. Met-misacylation is blocked by an inhibitor of cellular oxidases, implicating reactive oxygen species (ROS) as the misacylation trigger. Among six amino acids tested, tRNA misacylation occurs exclusively with Met. As Met residues are known to protect proteins against ROS-mediated damage, we propose that Met-misacylation functions adaptively to increase Met incorporation into proteins to protect cells against oxidative stress. In demonstrating an unexpected conditional aspect of decoding mRNA, our findings illustrate the importance of considering alternative iterations of the genetic code.

Genome-wide analysis of tRNA charging and activation of the eIF2 kinase Gcn2p.

When cells are subjected to nutritional stress, uncharged tRNAs accumulate and activate Gcn2p phosphorylation of eukaryotic initiation factor-2 (eIF2) and the general amino acid control pathway. The Gcn2p regulatory domain homologous to histidyl-tRNA synthetases is proposed to bind to uncharged tRNA, directly contributing to activation of Gcn2p. Here we apply a microarray technology to analyze genome-wide changes in tRNA charging in yeast upon activation of Gcn2p in response to amino acid starvation and high salinity, a stress not directly linked to nutritional deficiency. This microarray technology is applicable for all eukaryotic cells. Strains were starved for histidine, leucine, or tryptophan and shown to rapidly induce Gcn2p phosphorylation of eIF2. The relative charging level of all tRNAs was measured before and after starvation, and Gcn2p activation and the intracellular levels of the starved amino acid correlate with the observed decrease in tRNA charging. Interestingly, in some cases, tRNAs not charged with the starved amino acid became deacylated more rapidly than tRNAs charged with the starved amino acid. This increase in uncharged tRNA levels occurred although the intracellular levels for these non-starved amino acids remained unchanged. Additionally, treatment of a wild-type strain with high salinity stress showed transient changes in the charging of several different tRNAs. These results suggest that Gcn2p can be activated by many different tRNA species in the cell. These results also depict a complex cellular relationship between tRNA charging, amino acid availability, and non-nutrient stress. These relationships are best revealed by simultaneous monitoring of the charging level of all tRNAs.

SMN deficiency causes tissue-specific perturbations in the repertoire of snRNAs and widespread defects in splicing.

The survival of motor neurons (SMN) protein is essential for the biogenesis of small nuclear RNA (snRNA)-ribonucleoproteins (snRNPs), the major components of the pre-mRNA splicing machinery. Though it is ubiquitously expressed, SMN deficiency causes the motor neuron degenerative disease spinal muscular atrophy (SMA). We show here that SMN deficiency, similar to that which occurs in severe SMA, has unexpected cell type-specific effects on the repertoire of snRNAs and mRNAs. It alters the stoichiometry of snRNAs and causes widespread pre-mRNA splicing defects in numerous transcripts of diverse genes, preferentially those containing a large number of introns, in SMN-deficient mouse tissues. These findings reveal a key role for the SMN complex in RNA metabolism and in splicing regulation and indicate that SMA is a general splicing disease that is not restricted to motor neurons.

Tissue-specific differences in human transfer RNA expression.

Over 450 transfer RNA (tRNA) genes have been annotated in the human genome. Reliable quantitation of tRNA levels in human samples using microarray methods presents a technical challenge. We have developed a microarray method to quantify tRNAs based on a fluorescent dye-labeling technique. The first-generation tRNA microarray consists of 42 probes for nuclear encoded tRNAs and 21 probes for mitochondrial encoded tRNAs. These probes cover tRNAs for all 20 amino acids and 11 isoacceptor families. Using this array, we report that the amounts of tRNA within the total cellular RNA vary widely among eight different human tissues. The brain expresses higher overall levels of nuclear encoded tRNAs than every tissue examined but one and higher levels of mitochondrial encoded tRNAs than every tissue examined. We found tissue-specific differences in the expression of individual tRNA species, and tRNAs decoding amino acids with similar chemical properties exhibited coordinated expression in distinct tissue types. Relative tRNA abundance exhibits a statistically significant correlation to the codon usage of a collection of highly expressed, tissue-specific genes in a subset of tissues or tRNA isoacceptors. Our findings demonstrate the existence of tissue-specific expression of tRNA species that strongly implicates a role for tRNA heterogeneity in regulating translation and possibly additional processes in vertebrate organisms.

Real-time RNA profiling within a single bacterium.

Characterizing the dynamics of specific RNA levels requires real-time RNA profiling in a single cell. We show that the combination of a synthetic modular genetic system with fluorescence correlation spectroscopy allows us to directly measure in real time the activity of any specific promoter in prokaryotes. Using a simple inducible gene expression system, we found that induced RNA levels within a single bacterium of Escherichia coli exhibited a pulsating profile in response to a steady input of inducer. The genetic deletion of an efflux pump system, a key determinant of antibiotic resistance, altered the pulsating transcriptional dynamics and caused overexpression of induced RNA. In contrast with population measurements, real-time RNA profiling permits identifying relationships between genotypes and transcriptional dynamics that are accessible only at the level of the single cell.

Selective charging of tRNA isoacceptors induced by amino-acid starvation.

Aminoacylated (charged) transfer RNA isoacceptors read different messenger RNA codons for the same amino acid. The concentration of an isoacceptor and its charged fraction are principal determinants of the translation rate of its codons. A recent theoretical model predicts that amino-acid starvation results in 'selective charging' where the charging levels of some tRNA isoacceptors will be low and those of others will remain high. Here, we developed a microarray for the analysis of charged fractions of tRNAs and measured charging for all Escherichia coli tRNAs before and during leucine, threonine or arginine starvation. Before starvation, most tRNAs were fully charged. During starvation, the isoacceptors in the leucine, threonine or arginine families showed selective charging when cells were starved for their cognate amino acid, directly confirming the theoretical prediction. Codons read by isoacceptors that retain high charging can be used for efficient translation of genes that are essential during amino-acid starvation. Selective charging can explain anomalous patterns of codon usage in the genes for different families of proteins.

Exploring the regulation of tRNA distribution on the genomic scale.

Though up to 20% of the total RNA in bacterial cells is tRNA, the regulation of tRNA distribution on the genomic level remains unclear. tRNA distribution is governed by four processes: transcription, processing of precursor tRNA, degradation of precursor tRNA and degradation of mature tRNA. To elucidate the relationship between these processes in the regulation of tRNA production, the relative tRNA distribution was measured using a microarray specifically designed for tRNA. We developed a procedure that selectively labels 3'-CCA-containing RNAs with the fluorophores Cy3 or Cy5. The labeled tRNAs were then hybridized to microarrays printed with complementary DNA probes. The regulation of tRNA distribution in Bacillus subtilis was explored for a wild-type strain and a mutant strain with significantly decreased levels of RNase P, the enzyme required for the 5' maturation of all tRNA. The strains were either grown under a variety of conditions at doubling times ranging from 0.1 to 2.2 doublings per hour to investigate growth-related changes in the tRNA abundance or treated with the transcriptional inhibitor rifampicin to analyze mature tRNA degradation. Our results confirm that transcription and processing contribute significantly to the distribution of the 35 tRNA species in B.subtilis, and suggest a role for the degradation of precursor tRNA. Mature tRNA degradation occurs with little specificity for individual tRNA species and on the hour time-scale, indicating that degradation of mature tRNA plays only a minor role in the regulation of tRNA distribution. Aside from transcription, the final tRNA distribution appears to be derived from a balance between processing and precursor degradation activities.

Molecular characterization of the mouse involuted thymus: aberrations in expression of transcription regulators in thymocyte and epithelial compartments.

Despite playing a critical role in the development of naive T cells, the thymus is involuted with age. Whether a single age-associated defect or multiple aberrations contribute to thymic involution remains controversial. Here, we determined molecular aberrations in the thymocyte and epithelium compartments of the aging thymus. We demonstrated that total thymocyte numbers declined with a stepwise kinetics; clear demarcations occurred at 1.5, 3, 12 and 22 months of age. By quantitative PCR, a 2.4-fold reduction in the copies of signal joint TCR-excised circle (sjTREC)/10(5) thymocytes was first detected at 3 months; no further reduction observed thereafter. Nevertheless, the combined reductions in thymocyte numbers and sjTREC/10(5) cells caused a 7-fold decrease in sjTREC/thymus by 3 months, 21-fold by 18 months and 72-fold by 22 months as compared to 1 month. We showed aberration in expression of E2A, a transcription regulator critical for TCR beta rearrangement. While E2A expression declined 3-fold by 3 months and 18-fold by 7 months, expression of LMO2, a negative regulator of E2A activities, increased 5-fold by 18 months. Interestingly, expression of pre-T alpha and its transcriptional regulator HEB were not reduced with age. Furthermore, keratin-8 expression, specific for cortical thymic epithelium, declined 3-fold by 7 months and remained stable thereafter. In contrast, Foxn1 expression was reduced 3-fold by 3 months, 16-fold by 12 months and 37-fold by 18 months. IL-7 expression was not reduced until 7 months and reached 15-fold reduction by 22 months. Thus, the data demonstrate that thymic involution results not from a single defect, but culminates from an array of molecular aberrations in both the developing thymocytes and thymic epithelials.