Intermolecular mRNA-rRNA hybridization and the distribution of potential interaction regions in murine 18S rRNA

rRFtargetDB: a database of Ago1-mediated targets of ribosomal RNA fragments

rRNA-derived fragments (rRFs) are a class of emerging posttranscriptional regulators of gene expression likely binding to the transcripts of target genes. However, the lack of knowledge about such targets hinders our understanding of rRF functions or binding mechanisms. The paucity of resources supporting the identification of the targets of rRFs creates a bottleneck in the fast-developing field. We have previously analyzed chimeric reads in cross-linked Argonaute1-RNA complexes to help infer the guide-target pairs and binding mechanisms of multiple rRFs based on experimental data in human HEK293 cells. To efficiently disseminate these results to the research community, we designed a web-based database rRFtargetDB that preserves most of the experimental results after the removal of noise and has a user-friendly interface with flexible query options and filters allowing users to obtain comprehensive information on rRFs (or targets) of interest. rRFtargetDB is populated by ∼163,000 experimentally determined unique rRF-mRNA pairs (∼60,000 supported by ≥2 reads). Almost 30,000 rRF isoforms produced >385,000 (>156,000 with ≥2 reads) chimeras with all types of RNA targets (mRNAs and noncoding RNAs). Further analyses suggested hypothetical modes of interactions, supported by secondary structures of potential guide-target hybrids and binding motifs, essential for understanding the targeting mechanisms of rRFs. All these results (ranging from the weakest to the strongest experimental support) are presented in rRFtargetDB, whose goal is to provide a resource for building users' hypotheses on the potential roles of rRFs for experimental validation. Further, we illustrate the value/application of the database in several examples.rRFtargetDB is freely accessible at https://grigoriev-lab.camden.rutgers.edu/tardb.

Comparative analysis and classification of highly divergent mouse rDNA units based on their intergenic spacer (IGS) variability

Ribosomal DNA (rDNA) repeat units are organized into tandem clusters in eukaryotic cells. In mice, these clusters are located on at least eight chromosomes and show extensive variation in the number of repeats between mouse genomes. To analyze intra- and inter-genomic variation of mouse rDNA repeats, we selectively isolated 25 individual rDNA units using Transformation-Associated Recombination (TAR) cloning. Long-read sequencing and subsequent comparative sequence analysis revealed that each full-length unit comprises an intergenic spacer (IGS) and a ∼13.4 kb long transcribed region encoding the three rRNAs, but with substantial variability in rDNA unit size, ranging from ∼35 to ∼46 kb. Within the transcribed regions of rDNA units, we found 209 variants, 70 of which are in external transcribed spacers (ETSs); but the rDNA size differences are driven primarily by IGS size heterogeneity, due to indels containing repetitive elements and some functional signals such as enhancers. Further evolutionary analysis categorized rDNA units into distinct clusters with characteristic IGS lengths; numbers of enhancers; and presence/absence of two common SNPs in promoter regions, one of which is located within promoter (p)RNA and may influence pRNA folding stability. These characteristic features of IGSs also correlated significantly with 5'ETS variant patterns described previously and associated with differential expression of rDNA units. Our results suggest that variant rDNA units are differentially regulated and open a route to investigate the role of rDNA variation on nucleolar formation and possible associations with pathology.

Actively transcribed rDNA and distal junction (DJ) sequence are involved in association of NORs with nucleoli

Although they are organelles without a limiting membrane, nucleoli have an exclusive structure, built upon the rDNA-rich acrocentric short arms of five human chromosomes (nucleolar organizer regions or NORs). This has raised the question: what are the structural features of a chromosome required for its inclusion in a nucleolus? Previous work has suggested that sequences adjacent to the tandemly repeated rDNA repeat units (DJ, distal junction sequence) may be involved, and we have extended such studies by addressing several issues related to the requirements for the association of NORs with nucleoli. We exploited both a set of somatic cell hybrids containing individual human acrocentric chromosomes and a set of Human Artificial Chromosomes (HACs) carrying different parts of a NOR, including an rDNA unit or DJ or PJ (proximal junction) sequence. Association of NORs with nucleoli was increased when constituent rDNA was transcribed and may be also affected by the status of heterochromatin blocks formed next to the rDNA arrays. Furthermore, our data suggest that a relatively small size DJ region, highly conserved in evolution, is also involved, along with the rDNA repeats, in the localization of p-arms of acrocentric chromosomes in nucleoli. Thus, we infer a cooperative action of rDNA sequence-stimulated by its activity-and sequences distal to rDNA contributing to incorporation into nucleoli. Analysis of NOR sequences also identified LncRNA_038958 in the DJ, a candidate transcript with the region of the suggested promoter that is located close to the DJ/rDNA boundary and contains CTCF binding sites. This LncRNA may affect RNA Polymerase I and/or nucleolar activity. Our findings provide the basis for future studies to determine which RNAs and proteins interact critically with NOR sequences to organize the higher-order structure of nucleoli and their function in normal cells and pathological states.

A 3' UTR-derived small RNA connecting nitrogen and carbon metabolism in enteric bacteria

Increasing numbers of small, regulatory RNAs (sRNAs) corresponding to 3' untranslated regions (UTR) are being discovered in bacteria. One such sRNA, denoted GlnZ, corresponds to the 3' UTR of the Escherichia coli glnA mRNA encoding glutamine synthetase. Several forms of GlnZ, processed from the glnA mRNA, are detected in cells growing with limiting ammonium. GlnZ levels are regulated transcriptionally by the NtrC transcription factor and post-transcriptionally by RNase III. Consistent with the expression, E. coli cells lacking glnZ show delayed outgrowth from nitrogen starvation compared to wild type cells. Transcriptome-wide RNA-RNA interactome datasets indicated that GlnZ binds to multiple target RNAs. Immunoblots and assays of fusions confirmed GlnZ-mediated repression of glnP and sucA, encoding proteins that contribute to glutamine transport and the citric acid cycle, respectively. Although the overall sequences of GlnZ from E. coli K-12, Enterohemorrhagic E. coli and Salmonella enterica have significant differences due to various sequence insertions, all forms of the sRNA were able to regulate the two targets characterized. Together our data show that GlnZ impacts growth of E. coli under low nitrogen conditions by modulating genes that affect carbon and nitrogen flux.

Detecting RNA-RNA interactome

The last decade has seen a robust increase in various types of novel RNA molecules and their complexity in gene regulation. RNA molecules play a critical role in cellular events by interacting with other biomolecules, including protein, DNA, and RNA. It has been established that RNA-RNA interactions play a critical role in several biological processes by regulating the biogenesis and function of RNA molecules. Interestingly, RNA-RNA interactions regulate the biogenesis of diverse RNA molecules, including mRNAs, microRNAs, tRNAs, and circRNAs, through splicing or backsplicing. Structured RNAs like rRNA, tRNA, and snRNAs achieve their functional conformation by intramolecular RNA-RNA interactions. In addition, functional consequences of many intermolecular RNA-RNA interactions have been extensively studied in the regulation of gene expression. Hence, it is essential to understand the mechanism and functions of RNA-RNA interactions in eukaryotes. Conventionally, RNA-RNA interactions have been identified through diverse biochemical methods for decades. The advent of high-throughput RNA-sequencing technologies has revolutionized the identification of global RNA-RNA interactome in cells and their importance in RNA structure and function in gene expression regulation. Although these technologies revealed tens of thousands of intramolecular and intermolecular RNA-RNA interactions, we further look forward to future unbiased and quantitative high-throughput technologies for detecting transcriptome-wide RNA-RNA interactions. With the ability to detect RNA-RNA interactome, we expect that future studies will reveal the higher-order structures of RNA molecules and multi-RNA hybrids impacting human health and diseases. This article is categorized under: RNA Methods > RNA Analyses In Vitro and In Silico RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

Computational meta-analysis of ribosomal RNA fragments: potential targets and interaction mechanisms

The most abundant cellular RNA species, ribosomal RNA (rRNA), appears to be a source of massive amounts of non-randomly generated fragments. We found rRNA fragments (rRFs) in immunoprecipitated Argonaute (Ago-IP) complexes in human and mouse cells and in small RNA sequencing datasets. In human Ago1-IP, guanine-rich rRFs were preferentially cut in single-stranded regions of mature rRNAs between pyrimidines and adenosine, and non-randomly paired with cellular transcripts in crosslinked chimeras. Numerous identical rRFs were found in the cytoplasm and nucleus in mouse Ago2-IP. We report specific interaction motifs enriched in rRF-target pairs. Locations of such motifs on rRFs were compatible with the Ago structural features and patterns of the Ago-RNA crosslinking in both species. Strikingly, many of these motifs may bind to double-stranded regions on target RNAs, suggesting a potential pathway for regulating translation by unwinding mRNAs. Occurring on either end of rRFs and matching intronic, untranslated or coding regions in targets, such interaction sites extend the concept of microRNA seed regions. Targeting both borders of certain short introns, rRFs may be involved in their biogenesis or function, facilitated by Ago. Frequently dismissed as noise, rRFs are poised to greatly enrich the known functional spectrum of small RNA regulation.

A unique insert in the genomes of high-risk human papillomaviruses with a predicted dual role in conferring oncogenic risk

The differences between high risk and low risk human papillomaviruses (HR-HPV and LR-HPV, respectively) that contribute to the tumorigenic potential of HR-HPV are not well understood but can be expected to involve the HPV oncoproteins, E6 and E7. We combine genome comparison and machine learning techniques to identify a previously unnoticed insert near the 3’-end of the E6 oncoprotein gene that is unique to HR-HPV. Analysis of the insert sequence suggests that it exerts a dual effect, by creating a PDZ domain-binding motif at the C-terminus of E6, as well as eliminating the overlap between the E6 and E7 coding regions in HR-HPV. We show that, as a result, the insert might enable coupled termination-reinitiation of the E6 and E7 genes, supported by motifs complementary to the human 18S rRNA. We hypothesize that the added functionality of E6 and positive regulation of E7 expression jointly account for the tumorigenic potential of HR-HPV.

Stable G-quadruplex enabling sequences are selected against by the context-dependent codon bias

The distributions of secondary structural elements appear to differ between coding regions (CDS) of mRNAs compared to the untranslated regions (UTRs), presumably as a mechanism to fine-tune gene expression, including efficiency of translation. However, a systematic and comprehensive analysis of secondary structure avoidance because of potential bias in codon usage is difficult as some of the common secondary structures, such as, hairpins can be formed by numerous sequence combinations. Using G-quadruplex (GQ) as the model secondary structure we studied the impact of codon bias on GQs within the CDS. Because GQs can be predicted using specific consensus sequence motifs, they provide an excellent platform for investigation of the selectivity of such putative structures at the codon level. Using a bioinformatics approach, we calculated the frequencies of putative GQs within the CDS of a variety of species. Our results suggest that the most stable GQs appear to be significantly underrepresented within the CDS, through the use of specific synonymous codon combinations. Furthermore, we identified many peptide sequence motifs in which silent mutations can potentially alter translation via stable GQ formation. This work not only provides a comprehensive analysis on how stable secondary structures appear to be avoided within the CDS of mRNA, but also broadens the current understanding of synonymous codon usage as they relate to the structure-function relationship of RNA.

The Ribosome as a Missing Link in Prebiotic Evolution III: Over-Representation of tRNA- and rRNA-Like Sequences and Plieofunctionality of Ribosome-Related Molecules Argues for the Evolution of Primitive Genomes from Ribosomal RNA Modules

We propose that ribosomal RNA (rRNA) formed the basis of the first cellular genomes, and provide evidence from a review of relevant literature and proteonomic tests. We have proposed previously that the ribosome may represent the vestige of the first self-replicating entity in which rRNAs also functioned as genes that were transcribed into functional messenger RNAs (mRNAs) encoding ribosomal proteins. rRNAs also encoded polymerases to replicate itself and a full complement of the transfer RNAs (tRNAs) required to translate its genes. We explore here a further prediction of our “ribosome-first” theory: the ribosomal genome provided the basis for the first cellular genomes. Modern genomes should therefore contain an unexpectedly large percentage of tRNA- and rRNA-like modules derived from both sense and antisense reading frames, and these should encode non-ribosomal proteins, as well as ribosomal ones with key cell functions. Ribosomal proteins should also have been co-opted by cellular evolution to play extra-ribosomal functions. We review existing literature supporting these predictions. We provide additional, new data demonstrating that rRNA-like sequences occur at significantly higher frequencies than predicted on the basis of mRNA duplications or randomized RNA sequences. These data support our “ribosome-first” theory of cellular evolution.

Exploring the Role of AUG Triplets in Human Cap-Independent Translation Enhancing Elements

Cap-independent translation is believed to play an important role in eukaryotic protein synthesis, but the mechanisms of ribosomal recruitment and translation initiation remain largely unknown. Messenger RNA display was previously used to profile the human genome for RNA leader sequences that can enhance cap-independent translation. Surprisingly, many of the isolated sequences contain AUG triplets, suggesting a possible functional role for these motifs during translation initiation. Herein, we examine the sequence determinants of AUG triplets within a set of human translation enhancing elements (TEEs). Functional analyses performed in vitro and in cultured cells indicate that AUGs have the capacity to modulate mRNA translation either by serving as part of a larger ribosomal recruitment site or by directing the ribosome to defined initiation sites. These observations help constrain the functional role of AUG triplets in human TEEs and advance our understanding of this specific mechanism of cap-independent translation initiation.

Optimization of signal-to-noise ratio for efficient microarray probe design

MOTIVATION

Target-specific hybridization depends on oligo-probe characteristics that improve hybridization specificity and minimize genome-wide cross-hybridization. Interplay between specific hybridization and genome-wide cross-hybridization has been insufficiently studied, despite its crucial role in efficient probe design and in data analysis.

RESULTS

In this study, we defined hybridization specificity as a ratio between oligo target-specific hybridization and oligo genome-wide cross-hybridization. A microarray database, derived from the Genomic Comparison Hybridization (GCH) experiment and performed using the Affymetrix platform, contains two different types of probes. The first type of oligo-probes does not have a specific target on the genome and their hybridization signals are derived from genome-wide cross-hybridization alone. The second type includes oligonucleotides that have a specific target on the genomic DNA and their signals are derived from specific and cross-hybridization components combined together in a total signal. A comparative analysis of hybridization specificity of oligo-probes, as well as their nucleotide sequences and thermodynamic features was performed on the database. The comparison has revealed that hybridization specificity was negatively affected by low stability of the fully-paired oligo-target duplex, stable probe self-folding, G-rich content, including GGG motifs, low sequence complexity and nucleotide composition symmetry.

CONCLUSION

Filtering out the probes with defined 'negative' characteristics significantly increases specific hybridization and dramatically decreasing genome-wide cross-hybridization. Selected oligo-probes have two times higher hybridization specificity on average, compared to the probes that were filtered from the analysis by applying suggested cutoff thresholds to the described parameters. A new approach for efficient oligo-probe design is described in our study.

CONTACT

shabalin@ncbi.nlm.nih.gov or olga.matveeva@gmail.com

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The DAU cluster: a comparative analysis of 18 RHD alleles, some forming partial D antigens

BACKGROUND

The Rh system is the most complex and polymorphic blood group system in humans with more than 460 alleles known for the RHD gene. The DAU cluster of RHD alleles is characterized by the single-nucleotide change producing the p.Thr379Met amino acid substitution. It is called the DAU-0 allele and has been postulated to be the primordial allele, from which all other alleles of the DAU cluster have eventually evolved.

STUDY DESIGN AND METHODS

For two novel DAU alleles, the nucleotide sequences of all 10 exons as well as adjacent intronic regions, including the 5' and 3' untranslated regions (UTR), were determined for the RHD and RHCE genes. A phylogenetic tree for all DAU alleles was established using the neighbor-joining method with Pan troglodytes as root. Standard hemagglutination and flow cytometry tests were performed.

RESULTS

We characterized two DAU alleles, DAU-11 and DAU-5.1, closely related to DAU-3 and DAU-5, respectively. A phylogenetic analysis of the 18 known DAU alleles indicated point mutations and interallelic recombination contributing to diversification of the DAU cluster.

CONCLUSIONS

The DAU alleles encode a group of RhD protein variants, some forming partial D antigens known to permit anti-D in carriers; all are expected to cause anti-D alloimmunization in recipients of red blood cell transfusions. The DAU alleles evolved through genomic point mutations and recombination. These results suggest that the cluster of DAU alleles represent a clade, which is concordant with our previous postulate that they derived from the primordial DAU-0 allele.

Translational dysregulation in cancer: eIF4A isoforms and sequence determinants of eIF4A dependence

The malignant phenotype is largely the consequence of dysregulated gene expression. Transformed cells depend upon not just a global increase in protein synthesis but an altered translational landscape in which pro-oncogenic mRNAs are translationally up-regulated. Such mRNAs have been shown to possess longer and more structured 5'-UTRs requiring high levels of eukaryotic initiation factor 4A (eIF4A) helicase activity for efficient translation. As such there is a developing focus on targeting eIF4A as a cancer therapy. In order for such treatments to be successful, we must develop a detailed understanding of the mechanisms which make specific mRNAs more dependent on eIF4A activity than others. It is also crucial to fully characterize the potentially distinct roles of eIF4A1 and eIF4A2, which until recently were thought to be functionally interchangeable. This review will highlight the recent advances made in this field that address these issues.

A systematic computational analysis of the rRNA-3' UTR sequence complementarity suggests a regulatory mechanism influencing post-termination events in metazoan translation

Nucleic acid sequence complementarity underlies many fundamental biological processes. Although first noticed a long time ago, sequence complementarity between mRNAs and ribosomal RNAs still lacks a meaningful biological interpretation. Here we used statistical analysis of large-scale sequence data sets and high-throughput computing to explore complementarity between 18S and 28S rRNAs and mRNA 3' UTR sequences. By the analysis of 27,646 full-length 3' UTR sequences from 14 species covering both protozoans and metazoans, we show that the computed 18S rRNA complementarity creates an evolutionarily conserved localization pattern centered around the ribosomal mRNA entry channel, suggesting its biological relevance and functionality. Based on this specific pattern and earlier data showing that post-termination 80S ribosomes are not stably anchored at the stop codon and can migrate in both directions to codons that are cognate to the P-site deacylated tRNA, we propose that the 18S rRNA-mRNA complementarity selectively stabilizes post-termination ribosomal complexes to facilitate ribosome recycling. We thus demonstrate that the complementarity between 18S rRNA and 3' UTRs has a non-random nature and very likely carries information with a regulatory potential for translational control.

Tobacco BY-2 cell-free lysate: an alternative and highly-productive plant-based in vitro translation system

BACKGROUND

Cell-free protein synthesis is a rapid and efficient method for the production of recombinant proteins. Usage of prokaryotic cell-free extracts often leads to non-functional proteins. Eukaryotic counterparts such as wheat germ extract (WGE) and rabbit reticulocyte lysate (RLL) may improve solubility and promote the correct folding of eukaryotic multi-domain proteins that are difficult to express in bacteria. However, the preparation of WGEs is complex and time-consuming, whereas RLLs suffer from low yields. Here we report the development of a novel cell-free system based on tobacco Bright Yellow 2 (BY-2) cells harvested in the exponential growth phase.

RESULTS

The highly-productive BY-2 lysate (BYL) can be prepared quickly within 4-5 h, compared to 4-5 d for WGE. The efficiency of the BYL was tested using three model proteins: enhanced yellow fluorescent protein (eYFP) and two versions of luciferase. The added mRNA was optimized by testing different 5' and 3' untranslated regions (UTRs). The protein yield in batch and dialysis reactions using BYL was much higher than that of a commercial Promega WGE preparation, achieving a maximum yield of 80 μg/mL of eYFP and 100 μg/mL of luciferase, compared to only 45 μg/mL of eYFP and 35 μg/mL of luciferase in WGEs. In dialysis reactions, the BYL yielded about 400 μg/mL eYFP, representing up to 50% more of the target protein than the Promega WGE, and equivalent to the amount using 5Prime WGE system.

CONCLUSIONS

Due to the high yield and the short preparation time the BYL represents a remarkable improvement over current eukaryotic cell-free systems.

An evolutionary conserved pattern of 18S rRNA sequence complementarity to mRNA 5' UTRs and its implications for eukaryotic gene translation regulation

There are several key mechanisms regulating eukaryotic gene expression at the level of protein synthesis. Interestingly, the least explored mechanisms of translational control are those that involve the translating ribosome per se, mediated for example via predicted interactions between the ribosomal RNAs (rRNAs) and mRNAs. Here, we took advantage of robustly growing large-scale data sets of mRNA sequences for numerous organisms, solved ribosomal structures and computational power to computationally explore the mRNA–rRNA complementarity that is statistically significant across the species. Our predictions reveal highly specific sequence complementarity of 18S rRNA sequences with mRNA 5′ untranslated regions (UTRs) forming a well-defined 3D pattern on the rRNA sequence of the 40S subunit. Broader evolutionary conservation of this pattern may imply that 5′ UTRs of eukaryotic mRNAs, which have already emerged from the mRNA-binding channel, may contact several complementary spots on 18S rRNA situated near the exit of the mRNA binding channel and on the middle-to-lower body of the solvent-exposed 40S ribosome including its left foot. We discuss physiological significance of this structurally conserved pattern and, in the context of previously published experimental results, propose that it modulates scanning of the 40S subunit through 5′ UTRs of mRNAs.

Sounds of silence: synonymous nucleotides as a key to biological regulation and complexity

Messenger RNA is a key component of an intricate regulatory network of its own. It accommodates numerous nucleotide signals that overlap protein coding sequences and are responsible for multiple levels of regulation and generation of biological complexity. A wealth of structural and regulatory information, which mRNA carries in addition to the encoded amino acid sequence, raises the question of how these signals and overlapping codes are delineated along non-synonymous and synonymous positions in protein coding regions, especially in eukaryotes. Silent or synonymous codon positions, which do not determine amino acid sequences of the encoded proteins, define mRNA secondary structure and stability and affect the rate of translation, folding and post-translational modifications of nascent polypeptides. The RNA level selection is acting on synonymous sites in both prokaryotes and eukaryotes and is more common than previously thought. Selection pressure on the coding gene regions follows three-nucleotide periodic pattern of nucleotide base-pairing in mRNA, which is imposed by the genetic code. Synonymous positions of the coding regions have a higher level of hybridization potential relative to non-synonymous positions, and are multifunctional in their regulatory and structural roles. Recent experimental evidence and analysis of mRNA structure and interspecies conservation suggest that there is an evolutionary tradeoff between selective pressure acting at the RNA and protein levels. Here we provide a comprehensive overview of the studies that define the role of silent positions in regulating RNA structure and processing that exert downstream effects on proteins and their functions.

Multi-targeted priming for genome-wide gene expression assays

BACKGROUND

Complementary approaches to assaying global gene expression are needed to assess gene expression in regions that are poorly assayed by current methodologies. A key component of nearly all gene expression assays is the reverse transcription of transcribed sequences that has traditionally been performed by priming the poly-A tails on many of the transcribed genes in eukaryotes with oligo-dT, or by priming RNA indiscriminately with random hexamers. We designed an algorithm to find common sequence motifs that were present within most protein-coding genes of Saccharomyces cerevisiae and of Neurospora crassa, but that were not present within their ribosomal RNA or transfer RNA genes. We then experimentally tested whether degenerately priming these motifs with multi-targeted primers improved the accuracy and completeness of transcriptomic assays.

RESULTS

We discovered two multi-targeted primers that would prime a preponderance of genes in the genomes of Saccharomyces cerevisiae and Neurospora crassa while avoiding priming ribosomal RNA or transfer RNA. Examining the response of Saccharomyces cerevisiae to nitrogen deficiency and profiling Neurospora crassa early sexual development, we demonstrated that using multi-targeted primers in reverse transcription led to superior performance of microarray profiling and next-generation RNA tag sequencing. Priming with multi-targeted primers in addition to oligo-dT resulted in higher sensitivity, a larger number of well-measured genes and greater power to detect differences in gene expression.

CONCLUSIONS

Our results provide the most complete and detailed expression profiles of the yeast nitrogen starvation response and N. crassa early sexual development to date. Furthermore, our multi-targeting priming methodology for genome-wide gene expression assays provides selective targeting of multiple sequences and counter-selection against undesirable sequences, facilitating a more complete and precise assay of the transcribed sequences within the genome.

[Region 1112-1123 in the central domain of 18S rRNA in 40S subunits of plant ribosomes: accessibility for complementary interactions and the functional role]

The binding of the 18S RNA of the 40S subunits of wheat germ ribosomes to an oligodeoxyribonucleotide complementary to the 1112-1123 region of the central domain of this RNA molecule has been studied. The selective binding of this oligomer to the complementary RNA fragment and the inhibition of the translation of uncapped chimeric RNA containing enhancer sequences in the 5'-untranslated region upstream of the reporter sequence coding for beta-glucuronidase has been shown in a cell-free protein-synthesizing system. The use of a derivative of the aforementioned oligomer containing an alkylating group at the 5' end allowed for the demonstration that the 1112-1123 region of 18S RNA can form a heteroduplex with the complementary sequence of the oligomer. The data obtained show that the 1112-1123 region in loop 27 of the central domain of 18S RNA of 40S ribosomal subunits is exposed on the subunit surface and probably participates in the cap-independent binding of the subunits to mRNA due to the complementary interaction with the enhancer sequences.

Evolution of alternative and constitutive regions of mammalian 5'UTRs

BACKGROUND

Alternative splicing (AS) in protein-coding sequences has emerged as an important mechanism of regulation and diversification of animal gene function. By contrast, the extent and roles of alternative events including AS and alternative transcription initiation (ATI) within the 5'-untranslated regions (5'UTRs) of mammalian genes are not well characterized.

RESULTS

We evaluated the abundance, conservation and evolution of putative regulatory control elements, namely, upstream start codons (uAUGs) and open reading frames (uORFs), in the 5'UTRs of human and mouse genes impacted by alternative events. For genes with alternative 5'UTRs, the fraction of alternative sequences (those present in a subset of the transcripts) is much greater than that in the corresponding coding sequence, conceivably, because 5'UTRs are not bound by constraints on protein structure that limit AS in coding regions. Alternative regions of mammalian 5'UTRs evolve faster and are subject to a weaker purifying selection than constitutive portions. This relatively weak selection results in over-abundance of uAUGs and uORFs in the alternative regions of 5'UTRs compared to constitutive regions. Nevertheless, even in alternative regions, uORFs evolve under a stronger selection than the rest of the sequences, indicating that some of the uORFs are conserved regulatory elements; some of the non-conserved uORFs could be involved in species-specific regulation.

CONCLUSION

The findings on the evolution and selection in alternative and constitutive regions presented here are consistent with the hypothesis that alternative events, namely, AS and ATI, in 5'UTRs of mammalian genes are likely to contribute to the regulation of translation.

Low enzymatic activity haplotypes of the human catechol-O-methyltransferase gene: enrichment for marker SNPs

Catechol-O-methyltransferase (COMT) is an enzyme that plays a key role in the modulation of catechol-dependent functions such as cognition, cardiovascular function, and pain processing. Three common haplotypes of the human COMT gene, divergent in two synonymous and one nonsynonymous (val(158)met) position, designated as low (LPS), average (APS), and high pain sensitive (HPS), are associated with experimental pain sensitivity and risk of developing chronic musculoskeletal pain conditions. APS and HPS haplotypes produce significant functional effects, coding for 3- and 20-fold reductions in COMT enzymatic activity, respectively. In the present study, we investigated whether additional minor single nucleotide polymorphisms (SNPs), accruing in 1 to 5% of the population, situated in the COMT transcript region contribute to haplotype-dependent enzymatic activity. Computer analysis of COMT ESTs showed that one synonymous minor SNP (rs769224) is linked to the APS haplotype and three minor SNPs (two synonymous: rs6267, rs740602 and one nonsynonymous: rs8192488) are linked to the HPS haplotype. Results from in silico and in vitro experiments revealed that inclusion of allelic variants of these minor SNPs in APS or HPS haplotypes did not modify COMT function at the level of mRNA folding, RNA transcription, protein translation, or enzymatic activity. These data suggest that neutral variants are carried with APS and HPS haplotypes, while the high activity LPS haplotype displays less linked variation. Thus, both minor synonymous and nonsynonymous SNPs in the coding region are markers of functional APS and HPS haplotypes rather than independent contributors to COMT activity.

Identification of protein-coding sequences using the hybridization of 18S rRNA and mRNA during translation

We introduce a new approach in this article to distinguish protein-coding sequences from non-coding sequences utilizing a period-3, free energy signal that arises from the interactions of the 3′-terminal nucleotides of the 18S rRNA with mRNA. We extracted the special features of the amplitude and the phase of the period-3 signal in protein-coding regions, which is not found in non-coding regions, and used them to distinguish protein-coding sequences from non-coding sequences. We tested on all the experimental genes from Saccharomyces cerevisiae and Schizosaccharomyces pombe. The identification was consistent with the corresponding information from GenBank, and produced better performance compared to existing methods that use a period-3 signal. The primary tests on some fly, mouse and human genes suggests that our method is applicable to higher eukaryotic genes. The tests on pseudogenes indicated that most pseudogenes have no period-3 signal. Some exploration of the 3′-tail of 18S rRNA and pattern analysis of protein-coding sequences supported further our assumption that the 3′-tail of 18S rRNA has a role of synchronization throughout translation elongation process. This, in turn, can be utilized for the identification of protein-coding sequences.

Expression patterns of protein kinases correlate with gene architecture and evolutionary rates

Background

Protein kinase (PK) genes comprise the third largest superfamily that occupy ∼2% of the human genome. They encode regulatory enzymes that control a vast variety of cellular processes through phosphorylation of their protein substrates. Expression of PK genes is subject to complex transcriptional regulation which is not fully understood.

Principal Findings

Our comparative analysis demonstrates that genomic organization of regulatory PK genes differs from organization of other protein coding genes. PK genes occupy larger genomic loci, have longer introns, spacer regions, and encode larger proteins. The primary transcript length of PK genes, similar to other protein coding genes, inversely correlates with gene expression level and expression breadth, which is likely due to the necessity to reduce metabolic costs of transcription for abundant messages. On average, PK genes evolve slower than other protein coding genes. Breadth of PK expression negatively correlates with rate of non-synonymous substitutions in protein coding regions. This rate is lower for high expression and ubiquitous PKs, relative to low expression PKs, and correlates with divergence in untranslated regions. Conversely, rate of silent mutations is uniform in different PK groups, indicating that differing rates of non-synonymous substitutions reflect variations in selective pressure. Brain and testis employ a considerable number of tissue-specific PKs, indicating high complexity of phosphorylation-dependent regulatory network in these organs. There are considerable differences in genomic organization between PKs up-regulated in the testis and brain. PK genes up-regulated in the highly proliferative testicular tissue are fast evolving and small, with short introns and transcribed regions. In contrast, genes up-regulated in the minimally proliferative nervous tissue carry long introns, extended transcribed regions, and evolve slowly.

Conclusions/Significance

PK genomic architecture, the size of gene functional domains and evolutionary rates correlate with the pattern of gene expression. Structure and evolutionary divergence of tissue-specific PK genes is related to the proliferative activity of the tissue where these genes are predominantly expressed. Our data provide evidence that physiological requirements for transcription intensity, ubiquitous expression, and tissue-specific regulation shape gene structure and affect rates of evolution.

The presence of rRNA sequences in polyadenylated RNA and its potential functions

Accumulating evidence has shown that various lengths of ribosomal RNA (rRNA) sequences are widely present in polyadenylated RNA. This review article will discuss these polyadenylated rRNA containing transcripts (PART). PART are highly abundant and widely expressed in various tissues. It appears that there may be two types of PART. One type, type I, contains the rRNA segments (from approximately 10 nucleotides up to several hundred nucleotides) located within the transcripts. It has been demonstrated that short rRNA sequences within type I PART may function as cis-regulatory elements that regulate translational efficiency. The other type, type II, contains large portions or almost entire sequences of rRNA with a cap at the 5' end and poly(A) at 3' end. Recent work has shown that some type II PART have functional significance for some neurodegenerative disease processes and may play an important role in the pathogenesis of diseases. Further investigation in this area is critical to understanding the basic biology of PART and the potential role of PART in diseases.

Antisense masking reveals contributions of mRNA-rRNA base pairing to translation of Gtx and FGF2 mRNAs

We previously showed that a 9-nucleotide sequence from the 5' leader of the Gtx homeodomain mRNA facilitates translation initiation by base pairing to 18S rRNA. These earlier studies tested the Gtx element in isolation; we now assess the physiological relevance of this element in the context of two natural mRNAs that contain this sequence in their 5' leaders, Gtx itself and FGF2 (fibroblast growth factor 2). 2'-O-Methyl-modified RNA oligonucleotides were employed to block mRNA-rRNA base pairing by targeting either the Gtx-binding site in 18S rRNA or Gtx elements in recombinant mRNAs containing the Gtx or FGF2 5' leaders linked to a reporter cistron. Studies in cell-free lysates and transfected COS-7 cells showed that translation of mRNAs containing the Gtx or FGF2 5' leaders was decreased by > 50% when oligonucleotides targeting either the rRNA or mRNA were used. Specificity was demonstrated by showing that translation of the recombinant mRNAs was unaffected by control oligonucleotides. In addition, the specific oligonucleotides did not affect the translation of recombinant mRNAs in which the Gtx elements were mutated. Experiments performed using constructs containing Gtx and FGF2 5' leader and coding sequences ruled out possible effects of the reporter cistron. Furthermore, two-dimensional gel electrophoresis revealed that the oligonucleotides used in this study had little overall effect on the proteomes of cells transfected with these oligonucleotides. This study demonstrates that mRNA-rRNA base pairing affects the expression of two cellular mRNAs and describes a new approach for investigating putative mRNA-rRNA base pairing interactions in mammalian cells.

The ribosome filter redux

The ribosome filter hypothesis postulates that ribosomes are not simply translation machines but also function as regulatory elements that differentially affect or filter the translation of particular mRNAs. On the basis of new information, we take the opportunity here to review the ribosome filter hypothesis, suggest specific mechanisms of action, and discuss recent examples from the literature that support it.

Ribosomal frameshifting in decoding antizyme mRNAs from yeast and protists to humans: close to 300 cases reveal remarkable diversity despite underlying conservation

The protein antizyme is a negative regulator of intracellular polyamine levels. Ribosomes synthesizing antizyme start in one ORF and at the codon 5′ adjacent to its stop codon, shift +1 to a second and partially overlapping ORF which encodes most of the protein. The ribosomal frameshifting is a sensor and effector of an autoregulatory circuit which is conserved in animals, fungi and protists. Stimulatory signals encoded 5′ and 3′ of the shift site act to program the frameshifting. Despite overall conservation, many individual branches have evolved specific features surrounding the frameshift site. Among these are RNA pseudoknots, RNA stem-loops, conserved primary RNA sequences, nascent peptide sequences and branch-specific ‘shifty’ codons.

An mRNA-rRNA base-pairing mechanism for translation initiation in eukaryotes

Base-pairing of messenger RNA to ribosomal RNA is a mechanism of translation initiation in prokaryotes. Although analogous base-pairing has been suggested to affect the translation of various eukaryotic mRNAs, direct evidence has been lacking. To test such base-pairing, we developed a yeast system that uses ribosomes containing a mouse-yeast hybrid 18S rRNA. Using this system, we demonstrate that a 9-nucleotide element found in the mouse Gtx homeodomain mRNA facilitates translation initiation by base-pairing to 18S rRNA. Various point mutations in the Gtx element and in either the hybrid or wild-type yeast 18S rRNAs confirmed the requirement for an intact complementary match. The presence of the Gtx element in various mRNAs suggests that this element affects the translation of groups of mRNAs. We discuss the possibility that other mRNA elements affect translation by base-pairing to different sites in the 18S rRNA.

Cap-independent translation of tobacco etch virus is conferred by an RNA pseudoknot in the 5'-leader

The tobacco etch virus (TEV) 5'-leader promotes cap-independent translation in a 5'-proximal position and promotes internal initiation when present in the intercistronic region of a dicistronic mRNA, indicating that the leader contains an internal ribosome entry site. The TEV 143-nucleotide 5'-leader folds into a structure that contains two domains, each of which contains an RNA pseudoknot. Mutational analysis of the TEV 5'-leader identified pseudoknot (PK) 1 within the 5'-proximal domain and an upstream single-stranded region flanking PK1 as necessary to promote cap-independent translation. Mutations to either stem or to loops 2 or 3 of PK1 substantially disrupted cap-independent translation. The sequence of loop 3 in PK1 is complementary to a region in 18 S rRNA that is conserved throughout eukaryotes. Mutations within L3 that disrupted its potential base pairing with 18 S rRNA reduced cap-independent translation, whereas mutations that maintained the potential for base pairing with 18 S rRNA had little effect. These results indicated that the TEV 5'-leader functionally substitutes for a 5'-cap and promotes cap-independent translation through a 45-nucleotide pseudoknot-containing domain.

ARC-1, a sequence element complementary to an internal 18S rRNA segment, enhances translation efficiency in plants when present in the leader or intercistronic region of mRNAs

The sequences of different plant viral leaders with known translation enhancer ability show partial complementarity to the central region of 18S rRNA. Such complementarity might serve as a means to attract 40S ribosomal subunits and explain in part the translation-enhancing property of these sequences. To verify this notion, we designed beta-glucuronidase (GUS) mRNAs differing only in the nature of 10 nt inserts in the center of their 41 base leaders. These were complementary to consecutive domains of plant 18S rRNA. Sucrose gradient analysis revealed that leaders with inserts complementary to regions 1105-1114 and 1115-1124 ('ARC-1') of plant 18S rRNA bound most efficiently to the 40S ribosomal subunit after dissociation from 80S ribosomes under conditions of high ionic strength, a treatment known to remove translation initiation factors. Using wheat germ cell-free extracts, we could demonstrate that mRNAs with these leaders were translated more than three times more efficiently than a control lacking such a complementarity. Three linked copies of the insert enhanced translation of reporter mRNA to levels comparable with those directed by the natural translation enhancing leaders of tobacco mosaic virus and potato virus Y RNAs. Moreover, inserting the same leaders as intercistronic sequences in dicistronic mRNAs substantially increased translation of the second cistron, thereby revealing internal ribosome entry site activity. Thus, for plant systems, the complementary interaction between mRNA leader and the central region of 18S rRNA allows cap-independent binding of mRNA to the 43S pre-initiation complex without assistance of translation initiation factors.

Patterns in interspecies similarity correlate with nucleotide composition in mammalian 3'UTRs

Post-transcriptional regulation and the formation of mRNA 3' ends are crucial for gene expression in eukaryotes. Interspecies conservation of many sequences within 3'UTRs reveals selective constraint due to similar function. To study the pattern of conservation within 3'UTRs, we compiled and aligned 50 sets of complete orthologous 3'UTRs from four orders of mammals. We observed a mosaic pattern of conservation, with alternating regions of high (phylogenetic footprints) and low similarity. Conservation in 3'UTRs correlates with their base composition and also with the synonymous substitution rate in corresponding coding regions. The non-uniform distribution of conservation is more pronounced for 3'UTRs with a moderate or low level of overall conservation, where invariant nucleotides are more numerous, and their runs of lengths 4-7 occur more frequently than if conservation were random. Many runs of invariant nucleotides are AU-rich or pyrimidine-rich. Some of these runs coincide with known functional cis- elements of eukaryotic mRNAs, such as the U-rich upstream element, polyadenylation signal and DICE regulatory signal. More divergent regions of multiple alignments of 3'UTRs are often more G- and/or C-rich. Our results provide evidence on the importance of moderately conserved regions in 3'UTRs and suggest that regulatory functions of 3'UTRs might utilize gene-specific information in these regions.

Synergism of the 3'-untranslated region and an internal ribosome entry site differentially enhances the translation of a plant virus coat protein

The use of internal ribosome entry sites (IRESs) is one of the unorthodox mechanisms exploited by viruses to initiate the translation of internal genes. Herein, we report a plant virus exploiting an IRES and its 3'-untranslated region (UTR) to express its internal genes, notably the 3'-proximal viral coat protein gene. Hibiscus chlorotic ringspot virus (HCRSV), a positive-strand non-polyadenylated RNA virus, was demonstrated to harbor a unique 100-nucleotide (nt) IRES, located 124 nt upstream of the coat protein gene, that could function in wheat germ extract, rabbit reticulocyte lysate, and mammalian cells. In comparison with other known IRESs of picornaviruses and eukaryotic mRNAs, this 100-nt IRES is distinctively short and simple. The IRES activity was tested in homologous and heterologous bicistronic constructs, and the expression of the 3'-proximal gene was enhanced when the 3'-UTR was present. When the IRES element was bisected, each half still possessed IRES activity and could initiate internal translation on its own. Site-directed mutagenesis and deletion analyses revealed that the primary sequence within the 5' half was crucial for IRES activity, whereas the primary sequence of the second half and a GNRA motif were non-essential. To our knowledge, this is the first report describing a mechanism whereby an IRES, located in the 3' portion of the virus genome, co-operates with the 3'-UTR to enhance gene expression differentially.

The ribosome filter hypothesis

A variety of posttranscriptional mechanisms affects the processing, subcellular localization, and translation of messenger RNAs (mRNAs). Translational control appears to occur primarily at the initiation rather than the elongation stage. It has been suggested that translation is mediated largely by means of a cap-binding/scanning mechanism. On the basis of recent findings, we propose here that differential binding of particular mRNAs to eukaryotic 40S ribosomal subunits before translation may also selectively affect rates of polypeptide chain production. In this view, ribosomal subunits themselves are considered to be regulatory elements or filters that mediate interactions between particular mRNAs and components of the translation machinery. Differences in these interactions affect how efficiently individual mRNAs compete for ribosomal subunits. These competitive interactions would depend in part on the complementarity between sequences in mRNA and rRNA, as well as on structural differences among ribosomes in different cell types. By these means, translation may either be enhanced through increased recruitment of ribosomes or inhibited through strong interactions that sequester mRNAs. We propose that ribosomal filters may be important in cell differentiation and describe experimental tests for the filter hypothesis.

Complementarity between the mRNA 5' untranslated region and 18S ribosomal RNA can inhibit translation

In eubacteria, base pairing between the 3' end of 16S rRNA and the ribosome-binding site of mRNA is required for efficient initiation of translation. An interaction between the 18S rRNA and the mRNA was also proposed for translation initiation in eukaryotes. Here, we used an antisense RNA approach in vivo to identify the regions of 18S rRNA that might interact with the mRNA 5' untranslated region (5' UTR). Various fragments covering the entire mouse 18S rRNA gene were cloned 5' of a cat reporter gene in a eukaryotic vector, and translation products were analyzed after transient expression in human cells. For the largest part of 18S rRNA, we show that the insertion of complementary fragments in the mRNA 5' UTR do not impair translation of the downstream open reading frame (ORF). When translation inhibition is observed, reduction of the size of the complementary sequence to less than 200 nt alleviates the inhibitory effect. A single fragment complementary to the 18S rRNA 3' domain retains its inhibitory potential when reduced to 100 nt. Deletion analyses show that two distinct sequences of approximately 25 nt separated by a spacer sequence of 50 nt are required for the inhibitory effect. Sucrose gradient fractionation of polysomes reveals that mRNAs containing the inhibitory sequences accumulate in the fractions with 40S ribosomal subunits, suggesting that translation is blocked due to stalling of initiation complexes. Our results support an mRNA-rRNA base pairing to explain the translation inhibition observed and suggest that this region of 18S rRNA is properly located for interacting with mRNA.

Internal ribosome entry sites (IRESs): reality and use

IRESs are known to recruit ribosomes directly, without a previous scanning of untranslated region of mRNA by the ribosomes. IRESs have been found in a number of viral and cellular mRNAs. Experimentally, IRESs are commonly used to direct the expression of the second cistrons of bicistronic mRNAs. The mechanism of action of IRESs is not fully understood and a certain number of laboratories were not successful in using them in a reliable manner. Three observations done in our laboratory suggested that IRESs might not work as functionally as it was generally believed. Stem loops added before IRESs inhibited mRNA translation. When added into bicistronic mRNAs, IRESs initiated translation of the second cistrons efficiently only when the intercistronic region contained about 80 nucleotides, and they did not work any more effectively with intercistronic regions containing at least 300-400 nucleotides. Conversely, IRESs inserted at any position into the coding region of a cistron interrupted its translation and initiated translation of the following cistron. The first two data are hardly compatible with the idea that IRESs are able to recruit ribosomes without using the classical scanning mechanism. IRESs are highly structured and cannot be scanned by the 40S ribosomal subunit. We suggest that IRESs are short-circuited and are essentially potent stimulators favoring translation in particular physiological situations.

rRNA-complementarity in the 5' untranslated region of mRNA specifying the Gtx homeodomain protein: evidence that base- pairing to 18S rRNA affects translational efficiency

Numerous eukaryotic mRNAs contain sequences complementary to segments of the 18S and 28S rRNAs. Previous studies raised the possibility that these complementarities might allow mRNA-rRNA interactions and affect rates of translation. In the present study, we investigated the mRNA encoding the mouse Gtx homeodomain protein. This mRNA contains within its 5' untranslated region (UTR) a segment that is complementary to two regions of the 18S rRNA, located at nucleotides 701-741 and 1104-1136. A Gtx RNA probe containing this complementarity could be photochemically cross-linked to ribosomal subunits through a linkage to 18S rRNA but not to 28S rRNA. Oligonucleotide-directed RNase H digestion of the rRNA and a reverse transcription analysis localized the cross-linked probe to the complementary segment of 18S rRNA at nucleotides 1104-1136 but not at nucleotides 701-741. To determine whether complementarity in the Gtx mRNA affected translation, a mutational analysis was performed with a Gtx-luciferase fusion construct and four related constructs with altered complementarity to the 18S rRNA. These constructs were examined for their ability to be translated in cell-free lysates prepared from P19 embryonal carcinoma and C6 glioma cell lines and after cellular transfection into these same cell lines. In both cell-free translation and transfection studies, the rate of translation decreased more than 9-fold as the degree of complementarity to nucleotides 1104-1136 of the 18S rRNA increased. We hypothesize that segments complementary to rRNA, such as those contained within the Gtx mRNA, form a category of cis-acting regulatory elements in mRNAs that affect translation by base pairing to rRNA within ribosomes.

rRNA complementarity within mRNAs: a possible basis for mRNA-ribosome interactions and translational control

Our recent demonstration that many eukaryotic mRNAs contain sequences complementary to rRNA led to the hypothesis that these sequences might mediate specific interactions between mRNAs and ribosomes and thereby affect translation. In the present experiments, the ability of complementary sequences to bind to rRNA was investigated by using photochemical cross-linking. RNA probes with perfect complementarity to 18S or 28S rRNA were shown to cross-link specifically to the corresponding rRNA within intact ribosomal subunits. Similar results were obtained by using probes based on natural mRNA sequences with varying degrees of complementarity to the 18S rRNA. RNase H cleavage localized four such probes to complementary regions of the 18S rRNA. The effects of complementarity on translation were assessed by using the mRNA encoding ribosomal protein S15. This mRNA contains a sequence within its coding region that is complementary to the 18S rRNA at 20 of 22 nucleotides. RNA from an S15-luciferase fusion construct was translated in a cell-free lysate and compared with the translation of four related constructs that were mutated to decrease complementarity to the 18S rRNA. These mutations did not alter the amino acid sequence or the codon bias. A correlation between complementarity and translation was observed; constructs with less complementarity increased the amount of translation up to 54%. These findings raised the possibility that direct base-pairing of particular mRNAs to rRNAs within ribosomes may function as a mechanism of translational control.

rRNA-like sequences occur in diverse primary transcripts: implications for the control of gene expression

Many eukaryotic mRNAs contain sequences that resemble segments of 28S and 18S rRNAs, and these rRNA-like sequences are present in both the sense and antisense orientations. Some are similar to highly conserved regions of the rRNAs, whereas others have sequence similarities to expansion segments. In particular, four 18S rRNA-like sequences are found in several hundred different genes, and the location of these four sequences within the various genes is not random. One of these rRNA-like sequences is preferentially located within protein coding regions immediately upstream of the termination codon of a number of genes. Northern blot analysis of poly(A)+ RNA from different vertebrates (chicken, cattle, rat, mouse, and human) revealed that a large number of discrete RNA molecules hybridize at high stringency to cloned probes prepared from the 28S or 18S rRNA sequences that were found to match those in mRNAs. Inhibition of polymerase II activity, which prevents the synthesis of most mRNAs, abolished most of the hybridization to the rRNA probes. We consider the hypotheses that rRNA-like sequences may have spread throughout eukaryotic genomes and that their presence in primary transcripts may differentially affect gene expression.