tRNA-derived fragments target the small ribosomal subunit to fine-tune translation

Post-transcriptional cleavage of RNA molecules to generate smaller fragments is a widespread mechanism that enlarges the structural and functional complexity of cellular RNomes. In particular, fragments deriving from both precursor and mature tRNAs represent one of the rapidly growing classes of post-transcriptional RNA pieces. Importantly, these tRNA-derived fragments (tRFs) possess distinct expression patterns, abundance, cellular localizations, or biological roles compared with their parental tRNA molecules [1]. Here we present evidence that tRFs from the halophilic archaeon Haloferax volcanii directly bind to ribosomes. In a previous genomic screen for ribosome-associated small RNAs we have identified a 26 residue long fragment originating from the 5’ part of valine tRNA (Val-tRF) to be by far the most abundant tRF in H. volcanii [2]. The Val-tRF is processed in a stress-dependent manner and was found to primarily target the small ribosomal subunit in vitro and in vivo. Translational activity was markedly reduced in the presence of Val-tRF, while control RNA fragments of similar length did not show inhibition of protein biosynthesis. Crosslinking experiments and subsequent primer extension analyses revealed the Val-tRF interaction site to surround the mRNA path in the 30S subunit. In support of this, binding experiments demonstrated that Val-tRF does compete with mRNAs for ribosome binding. Therefore this tRF represents a ribosome-associated non-coding RNA (rancRNA) capable of regulating gene expression in H. volcanii under environmental stress conditions probably by fine-tuning the rate of protein production [3].

Post-transcriptional Regulation of Mammalian Gene Expression in Non-coding Region of Target RNA

Tumor Suppressor Candidate 2 (TUSC2) is a novel tumor suppressor gene located in the human chromosome 3p21.3 region. TUSC2 mRNA transcripts could be detected on Northern blots in both normal lung and some lung cancer cell lines, but no endogenous TUSC2 protein could be detected in a majority of lung cancer cell lines. Mechanisms regulating TUSC2 protein expression and its inactivation in primary lung cancer cells are largely unknown. We investigated the role of the 5’- and 3’-untranslated regions (UTRs) of the TUSC2 gene in the regulation of TUSC2 protein expression. We found that two small upstream open-reading frames (uORFs) in the 5’UTR of TUSC2 could markedly inhibit the translational initiation of TUSC2 protein by interfering with the “scanning” of the ribosome initiation complexes. Site-specific stem-loop array reverse transcription-polymerase chain reaction (SLA-RT-PCR) verified several micoRNAs (miRNAs) targeted at 3’UTR and directed TUSC2 cleavage and degradation. In addition, we used the established let-7-targeted high mobility group A2 (Hmga2) mRNA as a model system to study the mechanism of regulation of target mRNA by miRNAs in mammalian cells under physiological conditions. There have been no evidence of direct link between mRNA downregulation and mRNA cleavages mediated by miRNAs. Here we showed that the endonucleolytic cleavages on mRNAs were initiated by mammalian miRNA in seed pairing style. Let-7 directed cleavage activities among the eight predicted potential target sites have varied efficiency, which are influenced by the positional and the structural contexts in the UTR. The 5’ cleaved RNA fragments were mostly oligouridylated at their 3’-termini and accumulated for delayed 5’–3’ degradation. RNA fragment oligouridylation played important roles in marking RNA fragments for delayed bulk degradation and in converting RNA degradation mode from 3’–5’ to 5’–3’ with cooperative efforts from both endonucleolytic and non-catalytic miRNA-induced silencing complex (miRISC). Our findings point to a mammalian miRNA-mediated mechanism for the regulation of mRNA that miRNA can decrease target mRNA through target mRNA cleavage and uridine addition

RNA-SEQUENCING APPLICATIONS: GENE EXPRESSION QUANTIFICATION AND METHYLATOR PHENOTYPE IDENTIFICATION

My dissertation focuses on two aspects of RNA sequencing technology. The first is the methodology for modeling the overdispersion inherent in RNA-seq data for differential expression analysis. This aspect is addressed in three sections. The second aspect is the application of RNA-seq data to identify the CpG island methylator phenotype (CIMP) by integrating datasets of mRNA expression level and DNA methylation status. Section 1: The cost of DNA sequencing has reduced dramatically in the past decade. Consequently, genomic research increasingly depends on sequencing technology. However it remains elusive how the sequencing capacity influences the accuracy of mRNA expression measurement. We observe that accuracy improves along with the increasing sequencing depth. To model the overdispersion, we use the beta-binomial distribution with a new parameter indicating the dependency between overdispersion and sequencing depth. Our modified beta-binomial model performs better than the binomial or the pure beta-binomial model with a lower false discovery rate. Section 2: Although a number of methods have been proposed in order to accurately analyze differential RNA expression on the gene level, modeling on the base pair level is required. Here, we find that the overdispersion rate decreases as the sequencing depth increases on the base pair level. Also, we propose four models and compare them with each other. As expected, our beta binomial model with a dynamic overdispersion rate is shown to be superior. Section 3: We investigate biases in RNA-seq by exploring the measurement of the external control, spike-in RNA. This study is based on two datasets with spike-in controls obtained from a recent study. We observe an undiscovered bias in the measurement of the spike-in transcripts that arises from the influence of the sample transcripts in RNA-seq. Also, we find that this influence is related to the local sequence of the random hexamer that is used in priming. We suggest a model of the inequality between samples and to correct this type of bias. Section 4: The expression of a gene can be turned off when its promoter is highly methylated. Several studies have reported that a clear threshold effect exists in gene silencing that is mediated by DNA methylation. It is reasonable to assume the thresholds are specific for each gene. It is also intriguing to investigate genes that are largely controlled by DNA methylation. These genes are called “L-shaped” genes. We develop a method to determine the DNA methylation threshold and identify a new CIMP of BRCA. In conclusion, we provide a detailed understanding of the relationship between the overdispersion rate and sequencing depth. And we reveal a new bias in RNA-seq and provide a detailed understanding of the relationship between this new bias and the local sequence. Also we develop a powerful method to dichotomize methylation status and consequently we identify a new CIMP of breast cancer with a distinct classification of molecular characteristics and clinical features.

The U5 RNA of trypanosomes deviates from the canonical U5 RNA: The Leptomonas collosoma U5 RNA and its coding gene

Fractionation of the abundant small ribonucleoproteins (RNPs) of the trypanosomatid Leptomonas collosoma revealed the existence of a group of unidentified small RNPs that were shown to fractionate differently than the well-characterized trans-spliceosomal RNPs. One of these RNAs, an 80-nt RNA, did not possess a trimethylguanosine (TMG) cap structure but did possess a 5′ phosphate terminus and an invariant consensus U5 snRNA loop 1. The gene coding for the RNA was cloned, and the coding region showed 55% sequence identity to the recently described U5 homologue of Trypanosoma brucei [Dungan, J. D., Watkins, K. P. & Agabian, N. (1996) EMBO J. 15, 4016–4029]. The L. collosoma U5 homologue exists in multiple forms of RNP complexes, a 10S monoparticle, and two subgroups of 18S particles that either contain or lack the U4 and U6 small nuclear RNAs, suggesting the existence of a U4/U6⋅U5 tri-small nuclear RNP complex. In contrast to T. brucei U5 RNA (62 nt), the L. collosoma homologue is longer (80 nt) and possesses a second stem–loop. Like the trypanosome U3, U6, and 7SL RNA genes, a tRNA gene coding for tRNACys was found 98 nt upstream to the U5 gene. A potential for base pair interaction between U5 and SL RNA in the 5′ splice site region (positions −1 and +1) and downstream from it is proposed. The presence of a U5-like RNA in trypanosomes suggests that the most essential small nuclear RNPs are ubiquitous for both cis- and trans-splicing, yet even among the trypanosomatids the U5 RNA is highly divergent.

A Toolbox for Quantitative Gene Expression in Varroa destructor : RNA Degradation in Field Samples and Systematic Analysis of Reference Gene Stability

Funding: This project has received funding from the European Union’s Seventh Framework Programme for research, technological development and demonstration under grant agreement no. 613960 (SMARTBEES) (http://www.smartbees-fp7.eu/) and Veterinary Medicines Directorate, Department for Environment Food & Rural Affairs (Project # VM0517) (https://www.gov.uk/government/organisations/veterinary-medicines-directorate). CHM was supported by a Biosciences Knowledge Transfer Network Biotechnology and Biological Sciences Research Council (KTN-BBSRC CASE) Studentship (BB/L502467/1) (http://www.bbsrc.ac.uk/). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Acknowledgments We gratefully acknowledge Mr Sebastian Bacz’s expert help and advice with beekeeping.

RNA folding kinetics regulates translation of phage MS2 maturation gene

The gene for the maturation protein of the single-stranded RNA coliphage MS2 is preceded by an untranslated leader of 130 nt, which folds into a cloverleaf, i.e., three stem–loop structures enclosed by a long distance interaction (LDI). This LDI prevents translation because its 3′ moiety contains the Shine–Dalgarno sequence of the maturation gene. Previously, several observations suggested that folding of the cloverleaf is kinetically delayed, providing a time window for ribosomes to access the RNA. Here we present direct evidence for this model. In vitro experiments show that ribosome binding to the maturation gene is faster than refolding of the denatured cloverleaf. This folding delay appears related to special properties of the leader sequence. We have replaced the three stem–loop structures by a single five nt loop. This change does not affect the equilibrium structure of the LDI. Nevertheless, in this construct, the folding delay has virtually disappeared, suggesting that now the RNA folds faster than ribosomes can bind. Perturbation of the cloverleaf by an insertion makes the maturation start permanently accessible. A pseudorevertant that evolved from an infectious clone carrying the insertion had overcome this defect. It showed a wild-type folding delay before closing down the maturation gene. This experiment reveals the biological significance of retarded cloverleaf formation.

Noncytopathic Sindbis virus RNA vectors for heterologous gene expression

Infection of vertebrate cells with alphaviruses normally leads to prodigious expression of virus-encoded genes and a dramatic inhibition of host protein synthesis. Recombinant Sindbis viruses and replicons have been useful as vectors for high level foreign gene expression, but the cytopathic effects of viral replication have limited their use to transient studies. We recently selected Sindbis replicons capable of persistent, noncytopathic growth in BHK cells and describe here a new generation of Sindbis vectors useful for long-term foreign gene expression based on such replicons. Foreign genes of interest as well as the dominant selectable marker puromycin N-acteyltransferase, which confers resistance to the drug puromycin, were expressed as subgenomic transcripts of noncytopathic replicons or defective-interfering genomes complemented in trans by a replicon. Based on these strategies, we developed vectors that can be initiated via either RNA or DNA transfection and analyzed them for their level and stability of foreign gene expression. Noncytopathic Sindbis vectors express reasonably high levels of protein in nearly every cell. These vectors should prove to be flexible tools for the rapid expression of heterologous genes under conditions in which cellular metabolism is not perturbed, and we illustrate their utility with a number of foreign proteins.

Inhibition of gene expression in human cells through small molecule-RNA interactions

Small molecules that bind their biological receptors with high affinity and selectivity can be isolated from randomized pools of combinatorial libraries. RNA-protein interactions are important in many cellular functions, including transcription, RNA splicing, and translation. One example of such interactions is the mechanism of trans-activation of HIV-1 gene expression that requires the interaction of Tat protein with the trans-activation responsive region (TAR) RNA, a 59-base stem-loop structure located at the 5′ end of all nascent HIV-1 transcripts. Here we demonstrate the isolation of small TAR RNA-binding molecules from an encoded combinatorial library. We have made an encoded combinatorial tripeptide library of 24,389 possible members from d-and l-alpha amino acids on TentaGel resin. Using on-bead screening we have identified a small family of mostly heterochiral tripeptides capable of structure-specific binding to the bulge loop of TAR RNA. In vitro binding studies reveal stereospecific discrimination when the best tripeptide ligand is compared with diastereomeric peptide sequences. In addition, the most strongly binding tripeptide was shown to suppress transcriptional activation by Tat protein in human cells with an IC50 of ≈50 nM. Our results indicate that tripeptide RNA ligands are cell permeable, nontoxic to cells, and capable of inhibiting expression of specific genes by interfering with RNA-protein interactions.

Suppression of gene silencing: A general strategy used by diverse DNA and RNA viruses of plants

In transgenic and nontransgenic plants, viruses are both initiators and targets of a defense mechanism that is similar to posttranscriptional gene silencing (PTGS). Recently, it was found that potyviruses and cucumoviruses encode pathogenicity determinants that suppress this defense mechanism. Here, we test diverse virus types for the ability to suppress PTGS. Nicotiana benthamiana exhibiting PTGS of a green fluorescent protein transgene were infected with a range of unrelated viruses and various potato virus X vectors producing viral pathogenicity factors. Upon infection, suppression of PTGS was assessed in planta through reactivation of green fluorescence and confirmed by molecular analysis. These experiments led to the identification of three suppressors of PTGS and showed that suppression of PTGS is widely used as a counter-defense strategy by DNA and RNA viruses. However, the spatial pattern and degree of suppression varied extensively between viruses. At one extreme, there are viruses that suppress in all tissues of all infected leaves, whereas others are able to suppress only in the veins of new emerging leaves. This variation existed even between closely related members of the potexvirus group. Collectively, these results suggest that virus-encoded suppressors of gene silencing have distinct modes of action, are targeted against distinct components of the host gene-silencing machinery, and that there is dynamic evolution of the host and viral components associated with the gene-silencing mechanism.

Specific mutations in a viral RNA pseudoknot drastically change ribosomal frameshifting efficiency

Many viruses regulate protein synthesis by −1 ribosomal frameshifting using an RNA pseudoknot. Frameshifting is vital for viral reproduction. Using the information gained from the recent high-resolution crystal structure of the beet western yellow virus pseudoknot, a systematic mutational analysis has been carried out in vitro and in vivo. We find that specific nucleotide tertiary interactions at the junction between the two stems of the pseudoknot are crucial. A triplex is found between stem 1 and loop 2, and triplex interactions are required for frameshifting function. For some mutations, loss of one hydrogen bond is sufficient to abolish frameshifting. Furthermore, mutations near the 5′ end of the pseudoknot can increase frameshifting by nearly 300%, possibly by modifying ribosomal contacts. It is likely that the selection of suitable mutations can thus allow viruses to adjust frameshifting efficiencies and thereby regulate protein synthesis in response to environmental change.

Reciprocal subtraction differential RNA display: An efficient and rapid procedure for isolating differentially expressed gene sequences

A reciprocal subtraction differential RNA display (RSDD) approach has been developed that permits the rapid and efficient identification and cloning of both abundant and rare differentially expressed genes. RSDD comprises reciprocal subtraction of cDNA libraries followed by differential RNA display. The RSDD strategy was applied to analyze the gene expression alterations resulting during cancer progression as adenovirus-transformed rodent cells developed an aggressive transformed state, as documented by elevated anchorage-independence and enhanced in vivo oncogenesis in nude mice. This approach resulted in the identification and cloning of both known and a high proportion (>65%) of unknown sequences, including cDNAs displaying elevated expression as a function of progression (progression-elevated gene) and cDNAs displaying suppressed expression as a function of progression (progression-suppressed gene). Sixteen differentially expressed genes, including five unknown progression-elevated genes and six unknown progression-suppressed genes, have been characterized. The RSDD scheme should find wide application for the effective detection and isolation of differentially expressed genes.

Virus resistance and gene silencing in plants can be induced by simultaneous expression of sense and antisense RNA

Many examples of extreme virus resistance and posttranscriptional gene silencing of endogenous or reporter genes have been described in transgenic plants containing sense or antisense transgenes. In these cases of either cosuppression or antisense suppression, there appears to be induction of a surveillance system within the plant that specifically degrades both the transgene and target RNAs. We show that transforming plants with virus or reporter gene constructs that produce RNAs capable of duplex formation confer virus immunity or gene silencing on the plants. This was accomplished by using transcripts from one sense gene and one antisense gene colocated in the plant genome, a single transcript that has self-complementarity, or sense and antisense transcripts from genes brought together by crossing. A model is presented that is consistent with our data and those of other workers, describing the processes of induction and execution of posttranscriptional gene silencing.

A nucleolar protein related to ribosomal protein L7 is required for an early step in large ribosomal subunit biogenesis

The Saccharomyces cerevisiae Rlp7 protein has extensive identity and similarity to the large ribosomal subunit L7 proteins and shares an RNA-binding domain with them. Rlp7p is not a ribosomal protein; however, it is encoded by an essential gene and therefore must perform a function essential for cell growth. In this report, we show that Rlp7p is a nucleolar protein that plays a critical role in processing of precursors to the large ribosomal subunit RNAs. Pulse–chase labeling experiments with Rlp7p-depleted cells reveal that neither 5.8SS, 5.8SL, nor 25S is produced, indicating that both the major and minor processing pathways are affected. Analysis of processing intermediates by primer extension indicates that Rlp7p-depleted cells accumulate the 27SA3 precursor RNA, which is normally the major substrate (85%) used to produce the 5.8S and 25S rRNAs, and the ratio of 27SBL to 27SBS precursors changes from approximately 1:8 to 8:1 (depleted cells). Because 27SA3 is the direct precursor to 27SBS, we conclude that Rlp7p is specifically required for the 5′ to 3′ exonucleolytic trimming of the 27SA3 into the 27SBS precursor. As it is essential for processing in both the major and minor pathways, we propose that Rlp7p may act as a specificity factor that binds precursor rRNAs and tethers the enzymes that carry out the early 5′ to 3′ exonucleolytic reactions that generate the mature rRNAs. Rlp7p may also be required for the endonucleolytic cleavage in internal transcribed spacer 2 that separates the 5.8S rRNA from the 25S rRNA.

Yeast mutations in multiple complementation groups inhibit brome mosaic virus RNA replication and transcription and perturb regulated expression of the viral polymerase-like gene

Brome mosaic virus (BMV), a member of the alphavirus-like superfamily of positive-strand RNA viruses, encodes two proteins, 1a and 2a, that interact with each other, with unidentified host proteins, and with host membranes to form the viral RNA replication complex. Yeast expressing 1a and 2a support replication and subgenomic mRNA synthesis by BMV RNA3 derivatives. Using a multistep selection and screening process, we have isolated yeast mutants in multiple complementation groups that inhibit BMV-directed gene expression. Three complementation groups, represented by mutants mab1–1, mab2–1, and mab3–1 (for maintenance of BMV functions), were selected for initial study. Each of these mutants has a single, recessive, chromosomal mutation that inhibits accumulation of positive- and negative-strand RNA3 and subgenomic mRNA. BMV-directed gene expression was inhibited when the RNA replication template was introduced by in vivo transcription from DNA or by transfection of yeast with in vitro transcripts, confirming that cytoplasmic RNA replication steps were defective. mab1–1, mab2–1, and mab3–1 slowed yeast growth to varying degrees and were temperature-sensitive, showing that the affected genes contribute to normal cell growth. In wild-type yeast, expression of the helicase-like 1a protein increased the accumulation of 2a mRNA and the polymerase-like 2a protein, revealing a new level of viral regulation. In association with their other effects, mab1–1 and mab2–1 blocked the ability of 1a to stimulate 2a mRNA and protein accumulation, whereas mab3–1 had elevated 2a protein accumulation. Together, these results show that BMV RNA replication in yeast depends on multiple host genes, some of which directly or indirectly affect the regulated expression and accumulation of 2a.

Comparative mutational analysis of cis-acting RNA signals for translational frameshifting in HIV-1 and HTLV-2

Human immunodeficiency virus type 1 (HIV-1) and human T cell leukemia virus type II (HTLV-2) use a similar mechanism for –1 translational frameshifting to overcome the termination codon in viral RNA at the end of the gag gene. Previous studies have identified two important RNA signals for frameshifting, the slippery sequence and a downstream stem–loop structure. However, there have been somewhat conflicting reports concerning the individual contributions of these sequences. In this study we have performed a comprehensive mutational analysis of the cis-acting RNA sequences involved in HIV-1 gag–pol and HTLV-2 gag–pro frameshifting. Using an in vitro translation system we determined frameshifting efficiencies for shuffled HIV-1/HTLV-2 RNA elements in a background of HIV-1 or HTLV-2 sequences. We show that the ability of the slippery sequence and stem–loop to promote ribosomal frameshifting is influenced by the flanking upstream sequence and the nucleotides in the spacer element. A wide range of frameshift efficiency rates was observed for both viruses when shuffling single sequence elements. The results for HIV-1/HTLV-2 chimeric constructs represent strong evidence supporting the notion that the viral wild-type sequences are not designed for maximal frameshifting activity but are optimized to a level suited to efficient viral replication.