Aberrant mRNAs with extended 3' UTRs are substrates for rapid degradation by mRNA surveillance

The mRNA surveillance system is known to rapidly degrade aberrant mRNAs that contain premature termination codons in a process referred to as nonsense-mediated decay. A second class of aberrant mRNAs are those wherein the 3' UTR is abnormally extended due to a mutation in the polyadenylation site. We provide several observations that these abnormally 3'-extended mRNAs are degraded by the same machinery that degrades mRNAs with premature nonsense codons. First, the decay of the 3'-extended mRNAs is dependent on the same decapping enzyme and 5'-to-3' exonuclease. Second, the decay is also dependent on the proteins encoded by the UPF1, UPF2, and UPF3 genes, which are known to be specifically required for the rapid decay of mRNAs containing nonsense codons. Third, the ability of an extended 3' UTR to trigger decay is prevented by stabilizing sequences within the PGK1 coding region that are known to protect mRNAs from the rapid decay induced by premature nonsense codons. These results indicate that the mRNA surveillance system plays a role in degrading abnormally extended 3' UTRs. Based on these results, we propose a model in which the mRNA surveillance machinery degrades aberrant mRNAs due to the absence of the proper spatial arrangement of the translation-termination codon with respect to the 3' UTR element as defined by the utilization of a polyadenylation site.     

The conversion of 3' UTRs into coding regions

A possible origin of novel coding sequences is the removal of stop codons, leading to the inclusion of 3' untranslated regions (3' UTRs) within genes. We classified changes in the position of stop codons in closely related Saccharomyces species and in a mouse/rat comparison as either additions to or subtractions from coding regions. In both cases, the position of stop codons is highly labile, with more subtractions than additions found. The subtraction bias may be balanced by the input of new coding regions through gene duplication. Saccharomyces shows less stop codon lability than rodents, probably due to greater selective constraint. A higher proportion of 3' UTR incorporation events preserve frame in Saccharomyces. This higher proportion is consistent with the action of the [PSI(+)] prion as an evolutionary capacitor to facilitate 3' UTR incorporation in yeast.     

Characterization and evolution of 5' and 3' untranslated regions in eukaryotes

Untranslated regions (UTRs) in eukaryotes play a significant role in the regulation of translation and mRNA half-life, as well as interacting with specific RNA-binding proteins. However, UTRs receive less attention than more crucial elements such as genes, and the basic structural and evolutionary characteristics of UTRs of different species, and the relationship between these UTRs and the genome size and species gene number is not well understood. To address these questions, we performed a comparative analysis of 5' and 3' untranslated regions of different species by analyzing the basic characteristics of 244,976 UTRs from three eukaryote kingdoms (Plantae, Fungi, and Protista). The results showed that the UTR lengths and SSR frequencies in UTRs increased significantly with increasing species gene number while the length and G+C content in 5' UTRs and different types of repetitive sequences in 3' UTRs increased with the increase of genome size. We also found that the sequence length of 5' UTRs was significantly positively correlated with the presence of transposons and SSRs while the sequence length of 3' UTRs was significantly positively correlated with the presence of tandem repeat sequences. These results suggested that evolution of species complexity from lower organisms to higher organisms is accompanied by an increase in the regulatory complexity of UTRs, mediated by increasing UTR length, increasing G+C content of 5' UTRs, and insertion and expansion of repetitive sequences.     

Intron size, abundance, and distribution within untranslated regions of genes

Most research concerning the evolution of introns has largely considered introns within coding sequences (CDSs), without regard for introns located within untranslated regions (UTRs) of genes. Here, we directly determined intron size, abundance, and distribution in UTRs of genes using full-length cDNA libraries and complete genome sequences for four species, Arabidopsis thaliana, Drosophila melanogaster, human, and mouse. Overall intron occupancy (introns/exon kbp) is lower in 5' UTRs than CDSs, but intron density (intron occupancy in regions containing introns) tends to be higher in 5' UTRs than in CDSs. Introns in 5' UTRs are roughly twice as large as introns in CDSs, and there is a sharp drop in intron size at the 5' UTR-CDS boundary. We propose a mechanistic explanation for the existence of selection for larger intron size in 5' UTRs, and outline several implications of this hypothesis. We found introns to be randomly distributed within 5' UTRs, so long as a minimum required exon size was assumed. Introns in 3' UTRs were much less abundant than in 5' UTRs. Though this was expected for human and mouse that have intron-dependent nonsense-mediated decay (NMD) pathways that discourage the presence of introns within the 3' UTR, it was also true for A. thaliana and D. melanogaster, which may lack intron-dependent NMD. Our findings have several implications for theories of intron evolution and genome evolution in general.     

Presence of ATG triplets in 5' untranslated regions of eukaryotic cDNAs correlates with a 'weak' context of the start codon.

MOTIVATION: The context of the start codon (typically, AUG) and the features of the 5' Untranslated Regions (5' UTRs) are important for understanding translation regulation in eukaryotic mRNAs and for accurate prediction of the coding region in genomic and cDNA sequences. The presence of AUG triplets in 5' UTRs (upstream AUGs) might effect the initiation rate and, in the context of gene prediction, could reduce the accuracy of the identification of the authentic start. To reveal potential connections between the presence of upstream AUGs and other features of 5' UTRs, such as their length and the start codon context, we undertook a systematic analysis of the available eukaryotic 5' UTR sequences. RESULTS: We show that a large fraction of 5' UTRs in the available cDNA sequences, 15-53% depending on the organism, contain upstream ATGs. A negative correlation was observed between the information content of the translation start signal and the length of the 5' UTR. Similarly, a negative correlation exists between the 'strength' of the start context and the number of upstream ATGs. Typically, cDNAs containing long 5' UTRs with multiple upstream ATGs have a 'weak' start context, and in contrast, cDNAs containing short 5' UTRs without ATGs have 'strong' starts. These counter-intuitive results may be interpreted in terms of upstream AUGs having an important role in the regulation of translation efficiency by ensuring low basal translation level via double negative control and creating the potential for additional regulatory mechanisms. One of such mechanisms, supported by experimental studies of some mRNAs, includes removal of the AUG-containing portion of the 5' UTR by alternative splicing. AVAILABILITY: An ATG_ EVALUATOR program is available upon request or at www.itba.mi.cnr.it/webgene. CONTACT: rogozin@ncbi.nlm.nih.gov, milanesi@itba.mi.cnr.it.     

Control of translation by the 5'- and 3'-terminal regions of the dengue virus genome

The genomic RNAs of flaviviruses such as dengue virus (DEN) have a 5' m7GpppN cap like those of cellular mRNAs but lack a 3' poly(A) tail. We have studied the contributions to translational expression of 5'- and 3'-terminal regions of the DEN serotype 2 genome by using luciferase reporter mRNAs transfected into Vero cells. DCLD RNA contained the entire DEN 5' and 3' untranslated regions (UTRs), as well as the first 36 codons of the capsid coding region fused to the luciferase reporter gene. Capped DCLD RNA was as efficiently translated in Vero cells as capped GLGpA RNA, a reporter with UTRs from the highly expressed alpha-globin mRNA and a 72-residue poly(A) tail. Analogous reporter RNAs with regulatory sequences from West Nile and Sindbis viruses were also strongly expressed. Although capped DCLD RNA was expressed much more efficiently than its uncapped form, uncapped DCLD RNA was translated 6 to 12 times more efficiently than uncapped RNAs with UTRs from globin mRNA. The 5' cap and DEN 3' UTR were the main sources of the translational efficiency of DCLD RNA, and they acted synergistically in enhancing translation. The DEN 3' UTR increased mRNA stability, although this effect was considerably weaker than the enhancement of translational efficiency. The DEN 3' UTR thus has translational regulatory properties similar to those of a poly(A) tail. Its translation-enhancing effect was observed for RNAs with globin or DEN 5' sequences, indicating no codependency between viral 5' and 3' sequences. Deletion studies showed that translational enhancement provided by the DEN 3' UTR is attributable to the cumulative contributions of several conserved elements, as well as a nonconserved domain adjacent to the stop codon. One of the conserved elements was the conserved sequence (CS) CS1 that is complementary to cCS1 present in the 5' end of the DEN polyprotein open reading frame. Complementarity between CS1 and cCS1 was not required for efficient translation.     

Computer analyses of complete genomes suggest that some archaebacteria employ both eukaryotic and eubacterial mechanisms in translation initiation

The translation initiation mechanism of archaebacteria is still not clearly understood. Our previous work showed that ATG triplets before start codons have been strongly depleted in eukaryotic genomes, presumably because ribosome of eukaryotes scans mRNA from the 5' to 3' direction to find proper start codons. Extra ATG triplets before start codons would confuse the process and thus they have been negatively selected in eukaryotic genomes. In eubacterial genomes, on the other hand, ribosome binds to the Shine-Dalgarno (SD) sequence at once without mRNA scanning, and the characteristic patterns of ATG triplet depletion were not observed (Saito, R., Tomita, M., 1999. On negative selection against ATG triplets near start codons in eukaryotic and procaryotic genomes. J. Mol. Evol. 48, 213-217). The ATG triplet analysis on archaebacterial genomes revealed that Methanococcus jannaschii and Pyrococcus horikoshii show patterns similar to eukaryotes, implying that these species employ scanning of mRNA from the 5' to 3' direction in the process of translation initiation. On the other hand, our earlier study found that these archaea have SD-like sequences, which are complementary to the 3' end sequence of 16S rRNA, as in eubacterial translation initiation (Osada, Y., Saito, R., Tomita, M. Analysis of base-pairing potentials between 16S rRNA and 5' UTR for translation initiation in various procaryotes. Bioinformatics, in press). These two results combined lead us to conclude that these archaea probably use a hybrid mechanism; their ribosome scans mRNAs from the 5' to 3' direction and then 16S rRNA binds to the SD-like sequence of the 5' UTR.     

Eukaryotic mRNAs encoding abundant and scarce proteins are statistically dissimilar in many structural features

It is well known that non-coding mRNA sequences are dissimilar in many structural features. For individual mRNAs correlations were found for some of these features and their translational efficiency. However, no systematic statistical analysis was undertaken to relate protein abundance and structural characteristics of mRNA encoding the given protein. We have demonstrated that structural and contextual features of eukaryotic mRNAs encoding high- and low-abundant proteins differ in the 5′ untranslated regions (UTR). Statistically, 5′ UTRs of low-expression mRNAs are longer, their guanine plus cytosine content is higher, they have a less optimal context of the translation initiation codons of the main open reading frames and contain more frequently upstream AUG than 5′ UTRs of high-expression mRNAs. Apart from the differences in 5′ UTRs, high-expression mRNAs contain stronger termination signals. Structural features of low- and high-expression mRNAs are likely to contribute to the yield of their protein products.     

A test of the null model for 5' UTR evolution based on GC content

Eukaryotic mRNAs are headed by a stretch of noncoding sequence,the 5' untranslated region (UTR). It has been proposed thatthe length of 5' UTRs is selectively neutral and evolves undera process of stochastic destruction and recruitment of corepromoter elements, combined with selection against the prematureinitiation of translation. We test this null model by investigatingwhether 5' UTR length varies with genomic GC content, an implicitprediction of the model. Using simulations, we show that thenull model predicts a positive relationship between GC contentand UTR length for genes regulated by a TATA box. Although thisprediction is borne out qualitatively in genomic data from yeast,fruit flies, and humans, we find marked quantitative discrepancies.We conclude that UTR length may be shaped to some degree bythe forces considered in the null model but that the model failsto provide a complete explanation for UTR length evolution.