Comparative analysis of orthologous eukaryotic mRNAs: potential hidden functional signals

Kinetics of programmed and spontaneous ribosome sliding along the mRNA

The ribosome can slide along mRNA without establishing codon-anticodon interactions. This movement can be regulated (programmed) by the elements encoded in the mRNA, as observed in bypassing of non-coding gap in gene 60 of bacteriophage T4, or occur spontaneously, such as during traversal by the 70S ribosome of the 3'UTRs or upon re-initiation on bacterial polycistronic genes. In this study, we investigate the kinetic mechanism underlying the programmed and spontaneous ribosome sliding. We show that the translation rate of gene 60 mRNA decreases as the ribosome approaches the take-off site, especially when the KKYK regulatory sequence in the nascent peptide reaches the constriction site in the ribosome exit tunnel. However, efficiency of bypassing increases when the ribosome traverses the gap quickly. With the non-coding gap exceeding the natural 50 nt, the processivity of sliding remains high up to 56 nt, but drops sharply beyond that due to the loss of mRNA elements support. Sliding efficiency is temperature-dependent; while temperature regulates the number of ribosomes initiating programmed bypassing, traversing the long gaps becomes increasingly unfavorable at lower temperatures. This data offers novel insights into the kinetic determinants of programmed and spontaneous ribosome sliding along the mRNA.

The regulatory landscape of 5' UTRs in translational control during zebrafish embryogenesis

The 5' UTRs of mRNAs are critical for translation regulation, but their in vivo regulatory features are poorly characterized. Here, we report the regulatory landscape of 5' UTRs during early zebrafish embryogenesis using a massively parallel reporter assay of 18,154 sequences coupled to polysome profiling. We found that the 5' UTR is sufficient to confer temporal dynamics to translation initiation, and identified 86 motifs enriched in 5' UTRs with distinct ribosome recruitment capabilities. A quantitative deep learning model, DaniO5P, revealed a combined role for 5' UTR length, translation initiation site context, upstream AUGs and sequence motifs on in vivo ribosome recruitment. DaniO5P predicts the activities of 5' UTR isoforms and indicates that modulating 5' UTR length and motif grammar contributes to translation initiation dynamics. This study provides a first quantitative model of 5' UTR-based translation regulation in early vertebrate development and lays the foundation for identifying the underlying molecular effectors.

Highlights

In vivo MPRA systematically interrogates the regulatory potential of endogenous 5' UTRs The 5' UTR alone is sufficient to regulate the dynamics of ribosome recruitment during early embryogenesis The MPRA identifies 5' UTR cis -regulatory motifs for translation initiation control 5' UTR length, upstream AUGs and motif grammar contribute to the differential regulatory capability of 5' UTR switching isoforms.

The DEAD-Box RNA Helicase Ded1 Is Associated with Translating Ribosomes

DEAD-box RNA helicases are ATP-dependent RNA binding proteins and RNA-dependent ATPases that possess weak, nonprocessive unwinding activity in vitro, but they can form long-lived complexes on RNAs when the ATPase activity is inhibited. Ded1 is a yeast DEAD-box protein, the functional ortholog of mammalian DDX3, that is considered important for the scanning efficiency of the 48S pre-initiation complex ribosomes to the AUG start codon. We used a modified PAR-CLIP technique, which we call quicktime PAR-CLIP (qtPAR-CLIP), to crosslink Ded1 to 4-thiouridine-incorporated RNAs in vivo using UV light centered at 365 nm. The irradiation conditions are largely benign to the yeast cells and to Ded1, and we are able to obtain a high efficiency of crosslinking under physiological conditions. We find that Ded1 forms crosslinks on the open reading frames of many different mRNAs, but it forms the most extensive interactions on relatively few mRNAs, and particularly on mRNAs encoding certain ribosomal proteins and translation factors. Under glucose-depletion conditions, the crosslinking pattern shifts to mRNAs encoding metabolic and stress-related proteins, which reflects the altered translation. These data are consistent with Ded1 functioning in the regulation of translation elongation, perhaps by pausing or stabilizing the ribosomes through its ATP-dependent binding.

PhyloPGM: boosting regulatory function prediction accuracy using evolutionary information

Motivation

The computational prediction of regulatory function associated with a genomic sequence is of utter importance in -omics study, which facilitates our understanding of the underlying mechanisms underpinning the vast gene regulatory network. Prominent examples in this area include the binding prediction of transcription factors in DNA regulatory regions, and predicting RNA–protein interaction in the context of post-transcriptional gene expression. However, existing computational methods have suffered from high false-positive rates and have seldom used any evolutionary information, despite the vast amount of available orthologous data across multitudes of extant and ancestral genomes, which readily present an opportunity to improve the accuracy of existing computational methods.

Results

In this study, we present a novel probabilistic approach called PhyloPGM that leverages previously trained TFBS or RNA–RBP binding predictors by aggregating their predictions from various orthologous regions, in order to boost the overall prediction accuracy on human sequences. Throughout our experiments, PhyloPGM has shown significant improvement over baselines such as the sequence-based RNA–RBP binding predictor RNATracker and the sequence-based TFBS predictor that is known as FactorNet. PhyloPGM is simple in principle, easy to implement and yet, yields impressive results.

Availability and implementation

The PhyloPGM package is available at https://github.com/BlanchetteLab/PhyloPGM

Supplementary information

Supplementary data are available at Bioinformatics online.

Pairs of mutually compensatory frameshifting mutations contribute to protein evolution

Insertions and deletions of lengths not divisible by 3 in protein-coding sequences cause frameshifts that usually induce premature stop codons and may carry a high fitness cost. However, this cost can be partially offset by a second compensatory indel restoring the reading frame. The role of such pairs of compensatory frameshifting mutations (pCFMs) in evolution has not been studied systematically. Here, we use whole-genome alignments of protein-coding genes of 100 vertebrate species, and of 122 insect species, studying the prevalence of pCFMs in their divergence. We detect a total of 624 candidate pCFM genes; six of them pass stringent quality filtering, including three human genes: RAB36, ARHGAP6, and NCR3LG1. In some instances, amino acid substitutions closely predating or following pCFMs restored the biochemical similarity of the frameshifted segment to the ancestral amino acid sequence, possibly reducing or negating the fitness cost of the pCFM. Typically, however, the biochemical similarity of the frameshifted sequence to the ancestral one was not higher than the similarity of a random sequence of a protein-coding gene to its frameshifted version, indicating that pCFMs can uncover radically novel regions of protein space. In total, pCFMs represent an appreciable and previously overlooked source of novel variation in amino acid sequences.

Identification of permissive amber suppression sites for efficient non-canonical amino acid incorporation in mammalian cells

The genetic code of mammalian cells can be expanded to allow the incorporation of non-canonical amino acids (ncAAs) by suppressing in-frame amber stop codons (UAG) with an orthogonal pyrrolysyl-tRNA synthetase (PylRS)/tRNAPylCUA (PylT) pair. However, the feasibility of this approach is substantially hampered by unpredictable variations in incorporation efficiencies at different stop codon positions within target proteins. Here, we apply a proteomics-based approach to quantify ncAA incorporation rates at hundreds of endogenous amber stop codons in mammalian cells. With these data, we compute iPASS (Identification of Permissive Amber Sites for Suppression; available at www.bultmannlab.eu/tools/iPASS), a linear regression model to predict relative ncAA incorporation efficiencies depending on the surrounding sequence context. To verify iPASS, we develop a dual-fluorescence reporter for high-throughput flow-cytometry analysis that reproducibly yields context-specific ncAA incorporation efficiencies. We show that nucleotides up- and downstream of UAG synergistically influence ncAA incorporation efficiency independent of cell line and ncAA identity. Additionally, we demonstrate iPASS-guided optimization of ncAA incorporation rates by synonymous exchange of codons flanking the amber stop codon. This combination of in silico analysis followed by validation in living mammalian cells substantially simplifies identification as well as adaptation of sites within a target protein to confer high ncAA incorporation rates.

Reprogramming of translation in yeast cells impaired for ribosome recycling favors short, efficiently translated mRNAs

In eukaryotes, 43S preinitiation complex (PIC) formation is a rate-determining step of translation. Ribosome recycling following translation termination produces free 40S subunits for re-assembly of 43S PICs. Yeast mutants lacking orthologs of mammalian eIF2D (Tma64), and either MCT-1 (Tma20) or DENR (Tma22), are broadly impaired for 40S recycling; however, it was unknown whether this defect alters the translational efficiencies (TEs) of particular mRNAs. Here, we conducted ribosome profiling of a yeast tma64∆/tma20∆ double mutant and observed a marked reprogramming of translation, wherein the TEs of the most efficiently translated ('strong') mRNAs increase, while those of 'weak' mRNAs generally decline. Remarkably, similar reprogramming was seen on reducing 43S PIC assembly by inducing phosphorylation of eIF2α or by decreasing total 40S subunit levels by depleting Rps26. Our findings suggest that strong mRNAs outcompete weak mRNAs in response to 43S PIC limitation achieved in various ways, in accordance with previous mathematical modeling.

Uncovering deeply conserved motif combinations in rapidly evolving noncoding sequences

BACKGROUND

Animal genomes contain thousands of long noncoding RNA (lncRNA) genes, a growing subset of which are thought to be functionally important. This functionality is often mediated by short sequence elements scattered throughout the RNA sequence that correspond to binding sites for small RNAs and RNA binding proteins. Throughout vertebrate evolution, the sequences of lncRNA genes changed extensively, so that it is often impossible to obtain significant alignments between sequences of lncRNAs from evolutionary distant species, even when synteny is evident. This often prohibits identifying conserved lncRNAs that are likely to be functional or prioritizing constrained regions for experimental interrogation.

RESULTS

We introduce here LncLOOM, a novel algorithmic framework for the discovery and evaluation of syntenic combinations of short motifs. LncLOOM is based on a graph representation of the input sequences and uses integer linear programming to efficiently compare dozens of sequences that have thousands of bases each and to evaluate the significance of the recovered motifs. We show that LncLOOM is capable of identifying specific, biologically relevant motifs which are conserved throughout vertebrates and beyond in lncRNAs and 3'UTRs, including novel functional RNA elements in the CHASERR lncRNA that are required for regulation of CHD2 expression.

CONCLUSIONS

We expect that LncLOOM will become a broadly used approach for the discovery of functionally relevant elements in the noncoding genome.

eIF1 discriminates against suboptimal initiation sites to prevent excessive uORF translation genome-wide

The translation preinitiation complex (PIC) scans the mRNA for an AUG codon in a favorable context. Previous findings suggest that the factor eIF1 discriminates against non-AUG start codons by impeding full accommodation of Met-tRNA_i in the P site of the 40S ribosomal subunit, necessitating eIF1 dissociation for start codon selection. Consistent with this, yeast eIF1 substitutions that weaken its binding to the PIC increase initiation at UUG codons on a mutant his4 mRNA and particular synthetic mRNA reporters; and also at the AUG start codon of the mRNA for eIF1 itself owing to its poor Kozak context. It was not known however whether such eIF1 mutants increase initiation at suboptimal start codons genome-wide. By ribosome profiling, we show that the eIF1-L96P variant confers increased translation of numerous upstream open reading frames (uORFs) initiating with either near-cognate codons (NCCs) or AUGs in poor context. The increased uORF translation is frequently associated with the reduced translation of the downstream main coding sequences (CDS). Initiation is also elevated at certain NCCs initiating amino-terminal extensions, including those that direct mitochondrial localization of the GRS1 and ALA1 products, and at a small set of main CDS AUG codons with especially poor context, including that of eIF1 itself. Thus, eIF1 acts throughout the yeast translatome to discriminate against NCC start codons and AUGs in poor context; and impairing this function enhances the repressive effects of uORFs on CDS translation and alters the ratios of protein isoforms translated from near-cognate versus AUG start codons.

Tuning of NK-Specific HLA-C Expression by Alternative mRNA Splicing

A complex system regulating HLA-C expression in NK cells, driven by an NK-specific promoter that produces alternatively spliced variants of the 5'-UTR has been recently identified. Exon content of the NK-specific 5'-UTR varies strikingly across HLA-C alleles, with some exons being allele specific. In order to investigate the possibility that allelic variation in the 5'-UTR modulates HLA-C expression levels, cDNAs containing several distinct classes of 5'-UTR were compared. Subtle changes in 5'-UTR content had a significant effect on the expression of HLA-C ^* 03 and HLA-C ^* 12 cDNA clones, suggesting that alternative splicing can fine-tune the level of protein expression. The HLA-C ^* 06 allele was found to be highly expressed in relation to the other alleles studied. However, its increased expression was primarily associated with differences in the peptide-binding groove. Although the impact of allele-specific alternative splicing of NK-Pro transcripts on protein levels can be modest when compared with the effect of changes in peptide-loading, alternative splicing may represent an additional regulatory mechanism to fine-tune HLA-C levels within NK cells in distinct tissue environments or at different stages of maturation in order to achieve optimal levels of missing-self recognition.

Temperature-dependent regulation of upstream open reading frame translation in S. cerevisiae

BACKGROUND

Translation of an mRNA in eukaryotes starts at an AUG codon in most cases, but near-cognate codons (NCCs) such as UUG, ACG, and AUU can also be used as start sites at low levels in Saccharomyces cerevisiae. Initiation from NCCs or AUGs in the 5'-untranslated regions (UTRs) of mRNAs can lead to translation of upstream open reading frames (uORFs) that might regulate expression of the main ORF (mORF). Although there is some circumstantial evidence that the translation of uORFs can be affected by environmental conditions, little is known about how it is affected by changes in growth temperature.

RESULTS

Using reporter assays, we found that changes in growth temperature can affect translation from NCC start sites in yeast cells, suggesting the possibility that gene expression could be regulated by temperature by altering use of different uORF start codons. Using ribosome profiling, we provide evidence that growth temperature regulates the efficiency of translation of nearly 200 uORFs in S. cerevisiae. Of these uORFs, most that start with an AUG codon have increased translational efficiency at 37 °C relative to 30 °C and decreased efficiency at 20 °C. For translationally regulated uORFs starting with NCCs, we did not observe a general trend for the direction of regulation as a function of temperature, suggesting mRNA-specific features can determine the mode of temperature-dependent regulation. Consistent with this conclusion, the position of the uORFs in the 5'-leader relative to the 5'-cap and the start codon of the main ORF correlates with the direction of temperature-dependent regulation of uORF translation. We have identified several novel cases in which changes in uORF translation are inversely correlated with changes in the translational efficiency of the downstream main ORF. Our data suggest that translation of these mRNAs is subject to temperature-dependent, uORF-mediated regulation.

CONCLUSIONS

Our data suggest that alterations in the translation of specific uORFs by temperature can regulate gene expression in S. cerevisiae.

Functionally Significant Features in the 5' Untranslated Region of the ABCA1 Gene and Their Comparison in Vertebrates

Single nucleotide polymorphisms located in 5′ untranslated regions (5′UTRs) can regulate gene expression and have clinical impact. Recognition of functionally significant sequences within 5′UTRs is crucial in next-generation sequencing applications. Furthermore, information about the behavior of 5′UTRs during gene evolution is scarce. Using the example of the ATP-binding cassette transporter A1 (ABCA1) gene (Tangier disease), we describe our algorithm for functionally significant sequence finding. 5′UTR features (upstream start and stop codons, open reading frames (ORFs), GC content, motifs, and secondary structures) were studied using freely available bioinformatics tools in 55 vertebrate orthologous genes obtained from Ensembl and UCSC. The most conserved sequences were suggested as hot spots. Exon and intron enhancers and silencers (sc35, ighg2 cgamma2, ctnt, gh-1, and fibronectin eda exon), transcription factors (TFIIA, TATA, NFAT1, NFAT4, and HOXA13), some of them cancer related, and microRNA (hsa-miR-4474-3p) were localized to these regions. An upstream ORF, overlapping with the main ORF in primates and possibly coding for a small bioactive peptide, was also detected. Moreover, we showed several features of 5′UTRs, such as GC content variation, hairpin structure conservation or 5′UTR segmentation, which are interesting from a phylogenetic point of view and can stimulate further evolutionary oriented research.

Annotation and Expression of IDN2-like and FDM-like Genes in Sexual and Aposporous Hypericum perforatum L. accessions

The protein IDN2, together with the highly similar interactors FDM1 and FDM2, is required for RNA-directed DNA methylation (RdDM) and siRNA production. Epigenetic regulation of gene expression is required to restrict cell fate determination in A. thaliana ovules. Recently, three transcripts sharing high similarity with the A. thaliana IDN2 and FDM1-2 were found to be differentially expressed in ovules of apomictic Hypericum perforatum L. accessions. To gain further insight into the expression and regulation of these genes in the context of apomixis, we investigated genomic, transcriptional and functional aspects of the gene family in this species. The H. perforatum genome encodes for two IDN2-like and 7 FDM-like genes. Differential and heterochronic expression of FDM4-like genes was found in H. perforatum pistils. The involvement of these genes in reproduction and seed development is consistent with the observed reduction of the seed set and high variability in seed size in A. thaliana IDN2 and FDM-like knockout lines. Differential expression of IDN2-like and FDM-like genes in H. perforatum was predicted to affect the network of potential interactions between these proteins. Furthermore, pistil transcript levels are modulated by cytokinin and auxin but the effect operated by the two hormones depends on the reproductive phenotype.

Mutations in MAGT1 lead to a glycosylation disorder with a variable phenotype

Congenital disorders of glycosylation (CDG) are a group of rare metabolic diseases, due to impaired protein and lipid glycosylation. We identified two patients with defective serum transferrin glycosylation and mutations in the MAGT1 gene. These patients present with a phenotype that is mainly characterized by intellectual and developmental disability. MAGT1 has been described to be a subunit of the oligosaccharyltransferase (OST) complex and more specifically of the STT3B complex. However, it was also claimed that MAGT1 is a magnesium (Mg²⁺) transporter. So far, patients with mutations in MAGT1 were linked to a primary immunodeficiency, characterized by chronic EBV infections attributed to a Mg²⁺ homeostasis defect (XMEN). We compared the clinical and cellular phenotype of our two patients to that of an XMEN patient that we recently identified. All three patients have an N-glycosylation defect, as was shown by the study of different substrates, such as GLUT1 and SHBG, demonstrating that the posttranslational glycosylation carried out by the STT3B complex is dysfunctional in all three patients. Moreover, MAGT1 deficiency is associated with an enhanced expression of TUSC3, the homolog protein of MAGT1, pointing toward a compensatory mechanism. Hence, we delineate MAGT1-CDG as a disorder associated with two different clinical phenotypes caused by defects in glycosylation.

Tma64/eIF2D, Tma20/MCT-1, and Tma22/DENR Recycle Post-termination 40S Subunits In Vivo

The recycling of ribosomal subunits after translation termination is critical for efficient gene expression. Tma64 (eIF2D), Tma20 (MCT-1), and Tma22 (DENR) function as 40S recycling factors in vitro, but it is unknown whether they perform this function in vivo. Ribosome profiling of tma deletion strains revealed 80S ribosomes queued behind the stop codon, consistent with a block in 40S recycling. We found that unrecycled ribosomes could reinitiate translation at AUG codons in the 3' UTR, as evidenced by peaks in the footprint data and 3' UTR reporter analysis. In vitro translation experiments using reporter mRNAs containing upstream open reading frames (uORFs) further established that reinitiation increased in the absence of these proteins. In some cases, 40S ribosomes appeared to rejoin with 60S subunits and undergo an 80S reinitiation process in 3' UTRs. These results support a crucial role for Tma64, Tma20, and Tma22 in recycling 40S ribosomal subunits at stop codons and translation reinitiation.

Eukaryotic translational termination efficiency is influenced by the 3' nucleotides within the ribosomal mRNA channel

When a stop codon is at the 80S ribosomal A site, there are six nucleotides (+4 to +9) downstream that are inferred to be occupying the mRNA channel. We examined the influence of these downstream nucleotides on translation termination success or failure in mammalian cells at the three stop codons. The expected hierarchy in the intrinsic fidelity of the stop codons (UAA>UAG>>UGA) was observed, with highly influential effects on termination readthrough mediated by nucleotides at position +4 and position +8. A more complex influence was observed from the nucleotides at positions +5 and +6. The weakest termination contexts were most affected by increases or decreases in the concentration of the decoding release factor (eRF1), indicating that eRF1 binding to these signals was rate-limiting. When termination efficiency was significantly reduced by cognate suppressor tRNAs, the observed influence of downstream nucleotides was maintained. There was a positive correlation between experimentally measured signal strength and frequency of the signal in eukaryotic genomes, particularly in Saccharomyces cerevisiae and Drosophila melanogaster. We propose that termination efficiency is not only influenced by interrogation of the stop signal directly by the release factor, but also by downstream ribosomal interactions with the mRNA nucleotides in the entry channel.

eIF1A residues implicated in cancer stabilize translation preinitiation complexes and favor suboptimal initiation sites in yeast

The translation pre-initiation complex (PIC) scans the mRNA for an AUG codon in favorable context, and AUG recognition stabilizes a closed PIC conformation. The unstructured N-terminal tail (NTT) of yeast eIF1A deploys five basic residues to contact tRNA_i, mRNA, or 18S rRNA exclusively in the closed state. Interestingly, EIF1AX mutations altering the human eIF1A NTT are associated with uveal melanoma (UM). We found that substituting all five basic residues, and seven UM-associated substitutions, in yeast eIF1A suppresses initiation at near-cognate UUG codons and AUGs in poor context. Ribosome profiling of NTT substitution R13P reveals heightened discrimination against unfavorable AUG context genome-wide. Both R13P and K16D substitutions destabilize the closed complex at UUG codons in reconstituted PICs. Thus, electrostatic interactions involving the eIF1A NTT stabilize the closed conformation and promote utilization of suboptimal start codons. We predict UM-associated mutations alter human gene expression by increasing discrimination against poor initiation sites.

Non-AUG translation: a new start for protein synthesis in eukaryotes

This review by Kearse and Wilusz discusses the profound impact of non-AUG start codons in eukaryotic translation. It describes how misregulation of non-AUG initiation events contributes to multiple human diseases, including cancer and neurodegeneration, and how modulation of non-AUG usage may represent a novel therapeutic strategy.

Although it was long thought that eukaryotic translation almost always initiates at an AUG start codon, recent advancements in ribosome footprint mapping have revealed that non-AUG start codons are used at an astonishing frequency. These non-AUG initiation events are not simply errors but instead are used to generate or regulate proteins with key cellular functions; for example, during development or stress. Misregulation of non-AUG initiation events contributes to multiple human diseases, including cancer and neurodegeneration, and modulation of non-AUG usage may represent a novel therapeutic strategy. It is thus becoming increasingly clear that start codon selection is regulated by many trans-acting initiation factors as well as sequence/structural elements within messenger RNAs and that non-AUG translation has a profound impact on cellular states.

Efficient and Accurate Translation Initiation Directed by TISU Involves RPS3 and RPS10e Binding and Differential Eukaryotic Initiation Factor 1A Regulation

Canonical translation initiation involves ribosomal scanning, but short 5' untranslated region (5'UTR) mRNAs are translated in a scanning-independent manner. The extent and mechanism of scanning-independent translation are not fully understood. Here we report that short 5'UTR mRNAs constitute a substantial fraction of the translatome. Short 5'UTR mRNAs are enriched with TISU (translation initiator of short 5'UTR), a 12-nucleotide element directing efficient scanning-independent translation. Comprehensive mutagenesis revealed that each AUG codon-flanking nucleotide of TISU contributes to translational strength, but only a few are important for accuracy. Using site-specific UV cross-linking of ribosomal complexes assembled on TISU mRNA, we demonstrate specific binding of TISU to ribosomal proteins at the E and A sites. We identified RPS3 as the major TISU binding protein in the 48S complex A site. Upon 80S complex formation, RPS3 interaction is weakened and switched to RPS10e (formerly called RPS10). We further demonstrate that TISU is particularly dependent on eukaryotic initiation factor 1A (eIF1A) which interacts with both RPS3 and RPS10e. Our findings suggest that the cap-recruited ribosome specifically binds the TISU nucleotides at the A and E sites in cooperation with eIF1A to promote scanning arrest.

Nucleotide diversity analysis highlights functionally important genomic regions

We analyzed functionality and relative distribution of genetic variants across the complete Oryza sativa genome, using the 40 million single nucleotide polymorphisms (SNPs) dataset from the 3,000 Rice Genomes Project (http://snp-seek.irri.org), the largest and highest density SNP collection for any higher plant. We have shown that the DNA-binding transcription factors (TFs) are the most conserved group of genes, whereas kinases and membrane-localized transporters are the most variable ones. TFs may be conserved because they belong to some of the most connected regulatory hubs that modulate transcription of vast downstream gene networks, whereas signaling kinases and transporters need to adapt rapidly to changing environmental conditions. In general, the observed profound patterns of nucleotide variability reveal functionally important genomic regions. As expected, nucleotide diversity is much higher in intergenic regions than within gene bodies (regions spanning gene models), and protein-coding sequences are more conserved than untranslated gene regions. We have observed a sharp decline in nucleotide diversity that begins at about 250 nucleotides upstream of the transcription start and reaches minimal diversity exactly at the transcription start. We found the transcription termination sites to have remarkably symmetrical patterns of SNP density, implying presence of functional sites near transcription termination. Also, nucleotide diversity was significantly lower near 3′ UTRs, the area rich with regulatory regions.

Systematic analysis of the PTEN 5' leader identifies a major AUU initiated proteoform

Abundant evidence for translation within the 5' leaders of many human genes is rapidly emerging, especially, because of the advent of ribosome profiling. In most cases, it is believed that the act of translation rather than the encoded peptide is important. However, the wealth of available sequencing data in recent years allows phylogenetic detection of sequences within 5' leaders that have emerged under coding constraint and therefore allow for the prediction of functional 5' leader translation. Using this approach, we previously predicted a CUG-initiated, 173 amino acid N-terminal extension to the human tumour suppressor PTEN. Here, a systematic experimental analysis of translation events in the PTEN 5' leader identifies at least two additional non-AUG-initiated PTEN proteoforms that are expressed in most human cell lines tested. The most abundant extended PTEN proteoform initiates at a conserved AUU codon and extends the canonical AUG-initiated PTEN by 146 amino acids. All N-terminally extended PTEN proteoforms tested retain the ability to downregulate the PI3K pathway. We also provide evidence for the translation of two conserved AUG-initiated upstream open reading frames within the PTEN 5' leader that control the ratio of PTEN proteoforms.

High temperature attenuates the gravitropism of inflorescence stems by inducing SHOOT GRAVITROPISM 5 alternative splicing in Arabidopsis

In higher plants, gravitropism proceeds through three sequential steps in the responding organs: perception of gravity signals, signal transduction and asymmetric cell elongation. Light and temperature also influence the gravitropic orientation of plant organs. A series of Arabidopsis shoot gravitropism (sgr) mutants has been shown to exhibit disturbed shoot gravitropism. SGR5 is functionally distinct from other SGR members in that it mediates the early events of gravitropic responses in inflorescence stems. Here, we demonstrated that SGR5 alternative splicing produces two protein variants (SGR5α and SGR5β) in modulating the gravitropic response of inflorescence stems at high temperatures. SGR5β inhibits SGR5α function by forming non-DNA-binding heterodimers. Transgenic plants overexpressing SGR5β (35S:SGR5β) exhibit reduced gravitropic growth of inflorescence stems, as observed in the SGR5-deficient sgr5-5 mutant. Interestingly, SGR5 alternative splicing is accelerated at high temperatures, resulting in the high-level accumulation of SGR5β transcripts. When plants were exposed to high temperatures, whereas gravitropic curvature was reduced in Col-0 inflorescence stems, it was uninfluenced in the inflorescence stems of 35S:SGR5β transgenic plants and sgr5-5 mutant. We propose that the thermoresponsive alternative splicing of SGR5 provides an adaptation strategy by which plants protect the shoots from hot air under high temperature stress in natural habitats.

Extraordinary Genetic Diversity in a Wood Decay Mushroom

Populations of different species vary in the amounts of genetic diversity they possess. Nucleotide diversity π, the fraction of nucleotides that are different between two randomly chosen genotypes, has been known to range in eukaryotes between 0.0001 in Lynx lynx and 0.16 in Caenorhabditis brenneri. Here, we report the results of a comparative analysis of 24 haploid genotypes (12 from the United States and 12 from European Russia) of a split-gill fungus Schizophyllum commune. The diversity at synonymous sites is 0.20 in the American population of S. commune and 0.13 in the Russian population. This exceptionally high level of nucleotide diversity also leads to extreme amino acid diversity of protein-coding genes. Using whole-genome resequencing of 2 parental and 17 offspring haploid genotypes, we estimate that the mutation rate in S. commune is high, at 2.0 × 10⁻⁸ (95% CI: 1.1 × 10⁻⁸ to 4.1 × 10⁻⁸) per nucleotide per generation. Therefore, the high diversity of S. commune is primarily determined by its elevated mutation rate, although high effective population size likely also plays a role. Small genome size, ease of cultivation and completion of the life cycle in the laboratory, free-living haploid life stages and exceptionally high variability of S. commune make it a promising model organism for population, quantitative, and evolutionary genetics.

Divergent Alanyl-tRNA Synthetase Genes of Vanderwaltozyma polyspora Descended from a Common Ancestor through Whole-Genome Duplication Followed by Asymmetric Evolution

Cytoplasmic and mitochondrial forms of a eukaryotic aminoacyl-tRNA synthetase (aaRS) are generally encoded by two distinct nuclear genes, one of eukaryotic origin and the other of mitochondrial origin. However, in most known yeasts, only the mitochondrial-origin alanyl-tRNA synthetase (AlaRS) gene is retained and plays a dual-functional role. Here, we present a novel scenario of AlaRS evolution in the yeast Vanderwaltozyma polyspora. V. polyspora possesses two significantly diverged AlaRS gene homologues, one encoding the cytoplasmic form and the other its mitochondrial counterpart. Clever selection of transcription and translation initiation sites enables the two isoforms to be localized and thus functional in their respective cellular compartments. However, the two isoforms can also be stably expressed and function in the reciprocal compartments by insertion or removal of a mitochondrial targeting signal. Synteny and phylogeny analyses revealed that the AlaRS homologues of V. polyspora arose from a dual-functional common ancestor through whole-genome duplication (WGD). Moreover, the mitochondrial form had higher synonymous (1.6-fold) and nonsynonymous (2.8-fold) substitution rates than did its cytoplasmic counterpart, presumably due to a lesser constraint imposed on components of the mitochondrial translational apparatus. Our study suggests that asymmetric evolution confers the divergence between the AlaRS paralogues of V. polyspora.

Decoding mechanisms by which silent codon changes influence protein biogenesis and function

SCOPE

Synonymous codon usage has been a focus of investigation since the discovery of the genetic code and its redundancy. The occurrences of synonymous codons vary between species and within genes of the same genome, known as codon usage bias. Today, bioinformatics and experimental data allow us to compose a global view of the mechanisms by which the redundancy of the genetic code contributes to the complexity of biological systems from affecting survival in prokaryotes, to fine tuning the structure and function of proteins in higher eukaryotes. Studies analyzing the consequences of synonymous codon changes in different organisms have revealed that they impact nucleic acid stability, protein levels, structure and function without altering amino acid sequence. As such, synonymous mutations inevitably contribute to the pathogenesis of complex human diseases. Yet, fundamental questions remain unresolved regarding the impact of silent mutations in human disorders. In the present review we describe developments in this area concentrating on mechanisms by which synonymous mutations may affect protein function and human health.

PURPOSE

This synopsis illustrates the significance of synonymous mutations in disease pathogenesis. We review the different steps of gene expression affected by silent mutations, and assess the benefits and possible harmful effects of codon optimization applied in the development of therapeutic biologics.

PHYSIOLOGICAL AND MEDICAL RELEVANCE

Understanding mechanisms by which synonymous mutations contribute to complex diseases such as cancer, neurodegeneration and genetic disorders, including the limitations of codon-optimized biologics, provides insight concerning interpretation of silent variants and future molecular therapies.

Evolution at protein ends: major contribution of alternative transcription initiation and termination to the transcriptome and proteome diversity in mammals

Alternative splicing (AS), alternative transcription initiation (ATI) and alternative transcription termination (ATT) create the extraordinary complexity of transcriptomes and make key contributions to the structural and functional diversity of mammalian proteomes. Analysis of mammalian genomic and transcriptomic data shows that contrary to the traditional view, the joint contribution of ATI and ATT to the transcriptome and proteome diversity is quantitatively greater than the contribution of AS. Although the mean numbers of protein-coding constitutive and alternative nucleotides in gene loci are nearly identical, their distribution along the transcripts is highly non-uniform. On average, coding exons in the variable 5' and 3' transcript ends that are created by ATI and ATT contain approximately four times more alternative nucleotides than core protein-coding regions that diversify exclusively via AS. Short upstream exons that encompass alternative 5'-untranslated regions and N-termini of proteins evolve under strong nucleotide-level selection whereas in 3'-terminal exons that encode protein C-termini, protein-level selection is significantly stronger. The groups of genes that are subject to ATI and ATT show major differences in biological roles, expression and selection patterns.

Environmental manipulations generate bidirectional shifts in both behavior and gene regulation in a crossbred mouse model of extremes in trait anxiety

Although gene-environment interactions are known to significantly influence psychopathology-related disease states, only few animal models cover both the genetic background and environmental manipulations. Therefore, we have taken advantage of the bidirectionally inbred high (HAB) and low (LAB) anxiety-related behavior mouse lines to generate HAB × LAB F1 hybrids that intrinsically carry both lines’ genetic characteristics, and subsequently raised them in three different environments—standard, enriched (EE) and chronic mild stress (CMS). Assessing genetic correlates of trait anxiety, we focused on two genes already known to play a role in HAB vs. LAB mice, corticotropin releasing hormone receptor type 1 (Crhr1) and high mobility group nucleosomal binding domain 3 (Hmgn3). While EE F1 mice showed decreased anxiety-related and increased explorative behaviors compared to controls, CMS sparked effects in the opposite direction. However, environmental treatments affected the expression of the two genes in distinct ways. Thus, while expression ratios of Hmgn3 between the HAB- and LAB-specific alleles remained equal, total expression resembled the one observed in HAB vs. LAB mice, i.e., decreased after EE and increased after CMS treatment. On the other hand, while total expression of Crhr1 remained unchanged between the groups, the relative expression of HAB- and LAB-specific alleles showed a clear effect following the environmental modifications. Thus, the environmentally driven bidirectional shift of trait anxiety in this F1 model strongly correlated with Hmgn3 expression, irrespective of allele-specific expression patterns that retained the proportions of basic differential HAB vs. LAB expression, making this gene a match for environment-induced modifications. An involvement of Crhr1 in the bidirectional behavioral shift could, however, rather be due to different effects of the HAB- and LAB-specific alleles described here. Both candidate genes therefore deserve attention in the complex regulation of anxiety-related phenotypes including environment-mediated effects.

Translation initiation in eukaryotes: versatility of the scanning model

It is generally accepted that the initiation of translation in eukaryotes involves the binding of the 40S ribosomal subunit to the capped 5' end of an mRNA and subsequent scanning of 5' UTR in search of an initiation codon. However, until recently this has remained a mere hypothesis. This review describes the novel experimental evidence in support of this classical model. Data on the participation of various factors in the eukaryotic initiation process are summarized. The sequence of initiation events is described in light of the latest experimental data. The existing physical models of scanning are presented. Special attention is paid to discussion of alternative models of eukaryotic initiation of translation. It is demonstrated that the canonical mechanism of initiation is more versatile than previously thought.

Reinitiation and other unconventional posttermination events during eukaryotic translation

During ribosome recycling, posttermination complexes are dissociated by ABCE1 and eRF1 into 60S and tRNA/mRNA-associated 40S subunits, after which tRNA and mRNA are released by eIF1/eIF1A, Ligatin, or MCT-1/DENR. Occasionally, 40S subunits remain associated with mRNA and reinitiate at nearby AUGs. We recapitulated reinitiation using a reconstituted mammalian translation system. The presence of eIF2, eIF3, eIF1, eIF1A, and Met-tRNAi(Met) was sufficient for recycled 40S subunits to remain on mRNA, scan bidirectionally, and reinitiate at upstream and downstream AUGs if mRNA regions flanking the stop codon were unstructured. Imposition of 3' directionality additionally required eIF4F. Strikingly, posttermination ribosomes were not stably anchored on mRNA and migrated bidirectionally to codons cognate to the P site tRNA. Migration depended on the mode of peptide release (puromycin > eRF1⋅eRF3) and nature of tRNA and was enhanced by eEF2. The mobility of posttermination ribosomes suggests that some reinitiation events could involve 80S ribosomes rather than 40S subunits.

Hyper conserved elements in vertebrate mRNA 3'-UTRs reveal a translational network of RNA-binding proteins controlled by HuR

Little is known regarding the post-transcriptional networks that control gene expression in eukaryotes. Additionally, we still need to understand how these networks evolve, and the relative role played in them by their sequence-dependent regulatory factors, non-coding RNAs (ncRNAs) and RNA-binding proteins (RBPs). Here, we used an approach that relied on both phylogenetic sequence sharing and conservation in the whole mapped 3'-untranslated regions (3'-UTRs) of vertebrate species to gain knowledge on core post-transcriptional networks. The identified human hyper conserved elements (HCEs) were predicted to be preferred binding sites for RBPs and not for ncRNAs, namely microRNAs and long ncRNAs. We found that the HCE map identified a well-known network that post-transcriptionally regulates histone mRNAs. We were then able to discover and experimentally confirm a translational network composed of RNA Recognition Motif (RRM)-type RBP mRNAs that are positively controlled by HuR, another RRM-type RBP. HuR shows a preference for these RBP mRNAs bound in stem-loop motifs, confirming its role as a 'regulator of regulators'. Analysis of the transcriptome-wide HCE distribution revealed a profile of prevalently small clusters separated by unconserved intercluster RNA stretches, which predicts the formation of discrete small ribonucleoprotein complexes in the 3'-UTRs.

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.

Evolutionary patterns of non-coding RNAs

A plethora of new functions of non-coding RNAs (ncRNAs) have been discovered in past few years. In fact, RNA is emerging as the central player in cellular regulation, taking on active roles in multiple regulatory layers from transcription, RNA maturation, and RNA modification to translational regulation. Nevertheless, very little is known about the evolution of this "Modern RNA World" and its components. In this contribution, we attempt to provide at least a cursory overview of the diversity of ncRNAs and functional RNA motifs in non-translated regions of regular messenger RNAs (mRNAs) with an emphasis on evolutionary questions. This survey is complemented by an in-depth analysis of examples from different classes of RNAs focusing mostly on their evolution in the vertebrate lineage. We present a survey of Y RNA genes in vertebrates and study the molecular evolution of the U7 snRNA, the snoRNAs E1/U17, E2, and E3, the Y RNA family, the let-7 microRNA (miRNA) family, and the mRNA-like evf-1 gene. We furthermore discuss the statistical distribution of miRNAs in metazoans, which suggests an explosive increase in the miRNA repertoire in vertebrates. The analysis of the transcription of ncRNAs suggests that small RNAs in general are genetically mobile in the sense that their association with a hostgene (e.g. when transcribed from introns of a mRNA) can change on evolutionary time scales. The let-7 family demonstrates, that even the mode of transcription (as intron or as exon) can change among paralogous ncRNA.

The mechanism of eukaryotic translation initiation: new insights and challenges

Translation initiation in eukaryotes is a highly regulated and complex stage of gene expression. It requires the action of at least 12 initiation factors, many of which are known to be the targets of regulatory pathways. Here we review our current understanding of the molecular mechanics of eukaryotic translation initiation, focusing on recent breakthroughs from in vitro and in vivo studies. We also identify important unanswered questions that will require new ideas and techniques to solve.

In vitro activity of human translation initiation factor eIF4B is not affected by phosphomimetic amino acid substitutions S422D and S422E

Eukaryotic translation initiation factor eIF4B is necessary for ribosomal scanning through structured mRNA leaders. In higher eukaryotes, eIF4B serves as a downstream effector of several signaling pathways. In response to mitogenic stimuli, eIF4B undergoes multiple phosphorylations which are thought to regulate its activity. Recently, Ser422 was identified as a predominant site for human eIF4B phosphorylation via several signaling pathways, and phosphomimetic amino acid substitutions S422D or S422E were shown to activate eIF4B in living cells. However, stimulatory role of these modifications has never been analyzed directly. Here, using both mammalian reconstituted translation initiation assay and complete cell-free translation system, we perform a comparison of recombinant eIF4B derivatives with the wild type recombinant protein, and do not find any difference in their activities. On the contrary, native eIF4B purified from HeLa cells reveals significantly higher activity in both assays. Thus, the effects of S422D and S422E substitutions on eIF4B activity in living cells observed previously either require some other protein modification(s), or may only be manifested in an intact cell. Our study raises the question on whether the phosphorylation of Ser422 is sufficient for eIF4B activation observed upon mitogenic stimulation.

eIF4E-binding protein regulation of mRNAs with differential 5'-UTR secondary structure: a polyelectrostatic model for a component of protein-mRNA interactions

Control of translation in eukaryotes is complex, depending on the binding of various factors to mRNAs. Available data for subsets of mRNAs that are translationally up- and down-regulated in yeast eIF4E-binding protein (4E-BP) deletion mutants are coupled with reported mRNA secondary structure measurements to investigate whether 5′-UTR secondary structure varies between the subsets. Genes with up-regulated translational efficiencies in the caf20Δ mutant have relatively high averaged 5′-UTR secondary structure. There is no apparent wide-scale correlation of RNA-binding protein preferences with the increased 5′-UTR secondary structure, leading us to speculate that the secondary structure itself may play a role in differential partitioning of mRNAs between eIF4E/4E-BP repression and eIF4E/eIF4G translation initiation. Both Caf20p and Eap1p contain stretches of positive charge in regions of predicted disorder. Such regions are also present in eIF4G and have been reported to associate with mRNA binding. The pattern of these segments, around the canonical eIF4E-binding motif, varies between each 4E-BP and eIF4G. Analysis of gene ontology shows that yeast proteins containing predicted disordered segments, with positive charge runs, are enriched for nucleic acid binding. We propose that the 4E-BPs act, in part, as differential, flexible, polyelectrostatic scaffolds for mRNAs.

Molecular mechanism of scanning and start codon selection in eukaryotes

The correct translation of mRNA depends critically on the ability to initiate at the right AUG codon. For most mRNAs in eukaryotic cells, this is accomplished by the scanning mechanism, wherein the small (40S) ribosomal subunit attaches to the 5' end of the mRNA and then inspects the leader base by base for an AUG in a suitable context, using complementarity with the anticodon of methionyl initiator tRNA (Met-tRNAiMet) as the key means of identifying AUG. Over the past decade, a combination of yeast genetics, biochemical analysis in reconstituted systems, and structural biology has enabled great progress in deciphering the mechanism of ribosomal scanning. A robust molecular model now exists, describing the roles of initiation factors, notably eukaryotic initiation factor 1 (eIF1) and eIF1A, in stabilizing an "open" conformation of the 40S subunit with Met-tRNAiMet bound in a low-affinity state conducive to scanning and in triggering rearrangement into a "closed" conformation incompatible with scanning, which features Met-tRNAiMet more tightly bound to the "P" site and base paired with AUG. It has also emerged that multiple DEAD-box RNA helicases participate in producing a single-stranded "landing pad" for the 40S subunit and in removing the secondary structure to enable the mRNA to traverse the 40S mRNA-binding channel in the single-stranded form for base-by-base inspection in the P site.

Functional elements in initiation factors 1, 1A, and 2β discriminate against poor AUG context and non-AUG start codons

Yeast eIF1 inhibits initiation at non-AUG triplets, but it was unknown whether it also discriminates against AUGs in suboptimal context. As in other eukaryotes, the yeast gene encoding eIF1 (SUI1) contains an AUG in poor context, which could underlie translational autoregulation. Previously, eIF1 mutations were identified that increase initiation at UUG codons (Sui(-) phenotype), and we obtained mutations with the opposite phenotype of suppressing UUG initiation (Ssu(-) phenotype). Remarkably, Sui(-) mutations in eukaryotic translation initiation factor 1 (eIF1), eIF1A, and eIF2β all increase SUI1 expression in a manner diminished by introducing the optimal context at the SUI1 AUG, whereas Ssu(-) mutations in eIF1 and eIF1A decrease SUI1 expression with the native, but not optimal, context present. Therefore, discrimination against weak context depends on specific residues in eIFs 1, 1A, and 2β that also impede selection of non-AUGs, suggesting that context nucleotides and AUG act coordinately to stabilize the preinitiation complex. Although eIF1 autoregulates by discriminating against poor context in yeast and mammals, this mechanism does not prevent eIF1 overproduction in yeast, accounting for the hyperaccuracy phenotype afforded by SUI1 overexpression.

Translation initiation: a regulatory role for poly(A) tracts in front of the AUG codon in Saccharomyces cerevisiae

The 5'-UTR serves as the loading dock for ribosomes during translation initiation and is the key site for translation regulation. Many genes in the yeast Saccharomyces cerevisiae contain poly(A) tracts in their 5'-UTRs. We studied these pre-AUG poly(A) tracts in a set of 3274 recently identified 5'-UTRs in the yeast to characterize their effect on in vivo protein abundance, ribosomal density, and protein synthesis rate in the yeast. The protein abundance and the protein synthesis rate increase with the length of the poly(A), but exhibit a dramatic decrease when the poly(A) length is ≥12. The ribosomal density also reaches the lowest level when the poly(A) length is ≥12. This supports the hypothesis that a pre-AUG poly(A) tract can bind to translation initiation factors to enhance translation initiation, but a long (≥12) pre-AUG poly(A) tract will bind to Pab1p, whose binding size is 12 consecutive A residues in yeast, resulting in repression of translation. The hypothesis explains why a long pre-AUG poly(A) leads to more efficient translation initiation than a short one when PABP is absent, and why pre-AUG poly(A) is short in the early genes but long in the late genes of vaccinia virus.

Constraints and plasticity in genome and molecular-phenome evolution

Multiple constraints variously affect different parts of the genomes of diverse life forms. The selective pressures that shape the evolution of viral, archaeal, bacterial and eukaryotic genomes differ markedly, even among relatively closely related animal and bacterial lineages; by contrast, constraints affecting protein evolution seem to be more universal. The constraints that shape the evolution of genomes and phenomes are complemented by the plasticity and robustness of genome architecture, expression and regulation. Taken together, these findings are starting to reveal complex networks of evolutionary processes that must be integrated to attain a new synthesis of evolutionary biology.

Disruptive mRNA folding increases translational efficiency of catechol-O-methyltransferase variant

Catechol-O-methyltransferase (COMT) is a major enzyme controlling catecholamine levels that plays a central role in cognition, affective mood and pain perception. There are three common COMT haplotypes in the human population reported to have functional effects, divergent in two synonymous and one nonsynonymous position. We demonstrate that one of the haplotypes, carrying the non-synonymous variation known to code for a less stable protein, exhibits increased protein expression in vitro. This increased protein expression, which would compensate for lower protein stability, is solely produced by a synonymous variation (C¹⁶⁶T) situated within the haplotype and located in the 5′ region of the RNA transcript. Based on mRNA secondary structure predictions, we suggest that structural destabilization near the start codon caused by the T allele could be related to the observed increase in COMT expression. Our folding simulations of the tertiary mRNA structures demonstrate that destabilization by the T allele lowers the folding transition barrier, thus decreasing the probability of occupying its native state. These data suggest a novel structural mechanism whereby functional synonymous variations near the translation initiation codon affect the translation efficiency via entropy-driven changes in mRNA dynamics and present another example of stable compensatory genetic variations in the human population.

Reduced stability of mRNA secondary structure near the translation-initiation site in dsDNA viruses

Background

Recent studies have demonstrated a selection pressure for reduced mRNA secondary-structure stability near the start codon of coding sequences. This selection pressure can be observed in bacteria, archaea, and eukaryotes, and is likely caused by the requirement of efficient translation initiation in cellular organism.

Results

Here, we surveyed the complete genomes of 650 dsDNA virus strains for signals of reduced stability of mRNA secondary structure near the start codon. Our analysis included viruses infecting eukaryotic, prokaryotic, and archaeic hosts. We found that many viruses showed evidence for reduced mRNA secondary-structure stability near the start codon. The effect was most pronounced in viruses infecting prokaryotes, but was also observed in viruses infecting eukaryotes and archaea. The reduction in stability generally increased with increasing genomic GC content. For bacteriophage, the reduction was correlated with a corresponding reduction of stability in the phage hosts.

Conclusions

We conclude that reduced stability of the mRNA secondary structure near the start codon is a common feature for dsDNA viruses, likely driven by the same selective pressures that cause it in cellular organisms.

Identification of compounds that decrease the fidelity of start codon recognition by the eukaryotic translational machinery

Translation initiation in eukaryotes involves more than a dozen protein factors. Alterations in six factors have been found to reduce the fidelity of start codon recognition by the ribosomal preinitiation complex in yeast, a phenotype referred to as Sui(-). No small molecules are known that affect the fidelity of start codon recognition. Such compounds would be useful tools for probing the molecular mechanics of translation initiation and its regulation. To find compounds with this effect, we set up a high-throughput screen using a dual luciferase assay in S. cerevisiae. Screening of over 55,000 compounds revealed two structurally related molecules that decrease the fidelity of start codon selection by approximately twofold in the dual luciferase assay. This effect was confirmed using additional in vivo assays that monitor translation from non-AUG start codons. Both compounds increase translation of a natural upstream open reading frame previously shown to initiate translation at a UUG. The compounds were also found to exacerbate increased use of UUG as a start codon (Sui(-) phenotype) conferred by haploinsufficiency of wild-type eukaryotic initiation factor (eIF) 1, or by mutation in eIF1. Furthermore, the effects of the compounds are suppressed by overexpressing eIF1, which is known to restore the fidelity of start codon selection in strains harboring Sui(-) mutations in various other initiation factors. Together, these data strongly suggest that the compounds affect the translational machinery itself to reduce the accuracy of selecting AUG as the start codon.

Identification of evolutionarily conserved non-AUG-initiated N-terminal extensions in human coding sequences

In eukaryotes, it is generally assumed that translation initiation occurs at the AUG codon closest to the messenger RNA 5' cap. However, in certain cases, initiation can occur at codons differing from AUG by a single nucleotide, especially the codons CUG, UUG, GUG, ACG, AUA and AUU. While non-AUG initiation has been experimentally verified for a handful of human genes, the full extent to which this phenomenon is utilized--both for increased coding capacity and potentially also for novel regulatory mechanisms--remains unclear. To address this issue, and hence to improve the quality of existing coding sequence annotations, we developed a methodology based on phylogenetic analysis of predicted 5' untranslated regions from orthologous genes. We use evolutionary signatures of protein-coding sequences as an indicator of translation initiation upstream of annotated coding sequences. Our search identified novel conserved potential non-AUG-initiated N-terminal extensions in 42 human genes including VANGL2, FGFR1, KCNN4, TRPV6, HDGF, CITED2, EIF4G3 and NTF3, and also affirmed the conservation of known non-AUG-initiated extensions in 17 other genes. In several instances, we have been able to obtain independent experimental evidence of the expression of non-AUG-initiated products from the previously published literature and ribosome profiling data.

Defining potentially conserved RNA regulons of homologous zinc-finger RNA-binding proteins

Background

Glucose inhibition of gluconeogenic growth suppressor 2 protein (Gis2p) and zinc-finger protein 9 (ZNF9) are conserved yeast and human zinc-finger proteins. The function of yeast Gis2p is unknown, but human ZNF9 has been reported to bind nucleic acids, and mutations in the ZNF9 gene cause the neuromuscular disease myotonic dystrophy type 2. To explore the impact of these proteins on RNA regulation, we undertook a systematic analysis of the RNA targets and of the global implications for gene expression.

Results

Hundreds of mRNAs were associated with Gis2p, mainly coding for RNA processing factors, chromatin modifiers and GTPases. Target mRNAs contained stretches of G(A/U)(A/U) trinucleotide repeats located in coding sequences, which are sufficient for binding to both Gis2p and ZNF9, thus implying strong structural conservation. Predicted ZNF9 targets belong to the same functional categories as seen in yeast, indicating functional conservation, which is further supported by complementation of the large cell-size phenotype of gis2 mutants with ZNF9. We further applied a matched-sample proteome-transcriptome analysis suggesting that Gis2p differentially coordinates expression of RNA regulons, primarily by reducing mRNA and protein levels of genes required for ribosome assembly and by selectively up-regulating protein levels of myosins.

Conclusions

This integrated systematic exploration of RNA targets for homologous RNA-binding proteins indicates an unexpectedly high conservation of the RNA-binding properties and of potential targets, thus predicting conserved RNA regulons. We also predict regulation of muscle-specific genes by ZNF9, adding a potential link to the myotonic dystrophy related phenotypes seen in ZNF9 mouse models.

Evolution of prokaryotic genes by shift of stop codons

De novo origin of coding sequence remains an obscure issue in molecular evolution. One of the possible paths for addition (subtraction) of DNA segments to (from) a gene is stop codon shift. Single nucleotide substitutions can destroy the existing stop codon, leading to uninterrupted translation up to the next stop codon in the gene's reading frame, or create a premature stop codon via a nonsense mutation. Furthermore, short indels-caused frameshifts near gene's end may lead to premature stop codons or to translation past the existing stop codon. Here, we describe the evolution of the length of coding sequence of prokaryotic genes by change of positions of stop codons. We observed cases of addition of regions of 3'UTR to genes due to mutations at the existing stop codon, and cases of subtraction of C-terminal coding segments due to nonsense mutations upstream of the stop codon. Many of the observed stop codon shifts cannot be attributed to sequencing errors or rare deleterious variants segregating within bacterial populations. The additions of regions of 3'UTR tend to occur in those genes in which they are facilitated by nearby downstream in-frame triplets which may serve as new stop codons. Conversely, subtractions of coding sequence often give rise to in-frame stop codons located nearby. The amino acid composition of the added region is significantly biased, compared to the overall amino acid composition of the genes. Our results show that in prokaryotes, shift of stop codon is an underappreciated contributor to functional evolution of gene length.

Alternative translation start sites are conserved in eukaryotic genomes

Alternative start AUG codons within a single transcript can contribute to diversity of the proteome; however, their functional significance remains controversial. Here, we provide comparative genomics evidence that alternative start codons are under negative selection in vertebrates, insects and yeast. In genes where the annotated start codon (sAUG) resides within the suboptimal nucleotide context, the downstream in-frame AUG codons (dAUG) among the first ∼30 codon sites are significantly more conserved between species than in genes where the sAUG resides within the optimal context. Proteomics data show that this difference is not an annotation artifact and that dAUGs are in fact under selection as alternative start sites. The key optimal, and sometimes suboptimal, context-determining nucleotides of both the sAUG and dAUGs are conserved. Selection for secondary start sites is stronger in genes with the weak primary start site. Genes with multiple conserved start sites are enriched for transcription factors, and tend to have longer 5'UTRs and higher degree of alternative splicing. Together, these results imply that the use of alternative start sites by means of leaky mRNA scanning is a functional mechanism under selection for increased efficiency of translation and/or for translation of different N-terminal protein variants.

Identifying statistical dependence in genomic sequences via mutual information estimates

Questions of understanding and quantifying the representation and amount of information in organisms have become a central part of biological research, as they potentially hold the key to fundamental advances. In this paper, we demonstrate the use of information-theoretic tools for the task of identifying segments of biomolecules (DNA or RNA) that are statistically correlated. We develop a precise and reliable methodology, based on the notion of mutual information, for finding and extracting statistical as well as structural dependencies. A simple threshold function is defined, and its use in quantifying the level of significance of dependencies between biological segments is explored. These tools are used in two specific applications. First, they are used for the identification of correlations between different parts of the maize zmSRp32 gene. There, we find significant dependencies between the 5' untranslated region in zmSRp32 and its alternatively spliced exons. This observation may indicate the presence of as-yet unknown alternative splicing mechanisms or structural scaffolds. Second, using data from the FBI's combined DNA index system (CODIS), we demonstrate that our approach is particularly well suited for the problem of discovering short tandem repeats-an application of importance in genetic profiling.

A universal trend of reduced mRNA stability near the translation-initiation site in prokaryotes and eukaryotes

Recent studies have suggested that the thermodynamic stability of mRNA secondary structure near the start codon can regulate translation efficiency in Escherichia coli, and that translation is more efficient the less stable the secondary structure. We survey the complete genomes of 340 species for signals of reduced mRNA secondary structure near the start codon. Our analysis includes bacteria, archaea, fungi, plants, insects, fishes, birds, and mammals. We find that nearly all species show evidence for reduced mRNA stability near the start codon. The reduction in stability generally increases with increasing genomic GC content. In prokaryotes, the reduction also increases with decreasing optimal growth temperature. Within genomes, there is variation in the stability among genes, and this variation correlates with gene GC content, codon bias, and gene expression level. For birds and mammals, however, we do not find a genome-wide trend of reduced mRNA stability near the start codon. Yet the most GC rich genes in these organisms do show such a signal. We conclude that reduced stability of the mRNA secondary structure near the start codon is a universal feature of all cellular life. We suggest that the origin of this reduction is selection for efficient recognition of the start codon by initiator-tRNA.

Synonymous mutations are mutations that change the nucleotide sequence of a gene without changing the amino-acid sequence. Because these mutations don't alter the expressed protein, they are frequently also called silent mutations. Yet increasing evidence demonstrates that synonymous mutations are not that silent. In particular, experimental work in Escherichia coli has shown that the choice of synonymous codons near the start codon can greatly influence protein production. Codons that allow the mRNA to fold into a stable secondary structure seem to inhibit efficient translation initiation. This observation suggests that selection should prefer reduced mRNA stability near the start codon in many organisms. Here, we show that this prediction generally holds true in most organisms, including bacteria, archaea, fungi, plants, insects, and fishes. In birds and mammals it doesn't hold true genome-wide, but it does hold true in the most GC-rich genes. In all organisms, the extent to which mRNA stability is reduced increases with increasing GC content. In prokaryotes, it also increases with decreasing optimal growing temperature. Thus, it seems that all organisms have to optimize their synonymous sites near the start codon to guarantee efficient protein translation.

De Novo metabolic engineering and the promise of synthetic DNA

The uncertain price and tight supply of crude oil and the ever-increasing demand for clean energy have prompted heightened attention to the development of sustainable fuel technologies that ensure continued economic development while maintaining stewardship of the environment. In the face of these enormous challenges, biomass has emerged as a viable alternative to petroleum for the production of energy, chemicals, and materials owing to its abundance, inexpensiveness, and carbon-neutrality. Moreover, the immense ease and efficiency of biological systems at converting biomass-derived feedstocks into fuels, chemicals, and materials has generated renewed interest in biotechnology as a replacement for traditional chemical processes. Aided by the ever-expanding repertoire of microbial genetics and plant biotechnology, improved understanding of gene regulation and cellular metabolism, and incessantly accumulating gene and protein data, scientists are now contemplating engineering microbial cell factories to produce fuels, chemical feedstocks, polymers and pharmaceuticals in an economically and environmentally sustainable way. This goal resonates with that of metabolic engineering - the improvement of cellular properties through the intelligent design, rational modification, or directed evolution of biochemical pathways, and arguably, metabolic engineering seems best positioned to achieve the concomittant goals of environmental stewardship and economic prolificity.Improving a host organism's cellular traits and the potential design of new phenotypes is strongly dependent on the ability to effectively control the organism's genetic machinery. In fact, finely-tuned gene expression is imperative for achieving an optimal balance between pathway expression and cell viability, while avoiding cytotoxicity due to accumulation of certain gene products or metabolites. Early attempts to engineer a cell's metabolism almost exclusively relied on merely deleting or over-expressing single or multiple genes using recombinant DNA, and intervention targets were predominantly selected based on knowledge of the stoichiometry, kinetics, and regulation of the pathway of interest. However, the distributive nature of metabolic control, as opposed to the existence of a single rate-limiting step, predicates the controlled expression of multiple enzymes in several coordinated pathways to achieve the desired flux, and, as such, simple strategies involving either deleting or over-expressing genes are greatly limited in this context. On the other hand, the use of synthetic or modified promoters, riboswitches, tunable intergenic regions, and translation modulators such as internal ribosome entry sequences, upstream open reading frames, optimized mRNA secondary structures, and RNA silencing have been shown to be enormously conducive to achieving the fine-tuning of gene expression. These modifications to the genetic machinery of the host organism can be best achieved via the use of synthetic DNA technology, and the constant improvement in the affordability and quality of oligonucleotide synthesis suggests that these might well become the mainstay of the metabolic engineering toolbox in the years to come. The possibilities that arise with the use of synthetic oligonucleotides will be delineated herein.