uORFs, reinitiation and alternative translation start sites in human mRNAs

TP53BP1, a dual-coding gene, uses promoter switching and translational reinitiation to express a smORF protein

The complexity of the metazoan proteome is significantly increased by the expression of small proteins (<100 aa) derived from smORFs within lncRNAs, uORFs, 3' UTRs and, reading frames overlapping the CDS. These smORF encoded proteins (SEPs) have diverse roles, ranging from the regulation of cellular physiological to essential developmental functions. We report the characterization of a new member of this protein family, SEP53BP1, derived from a small internal ORF that overlaps the CDS encoding 53BP1. Its expression is coupled to the utilization of an alternative, cell-type specific promoter coupled to translational reinitiation events mediated by a uORF in the alternative 5' TL of the mRNA. This uORF-mediated reinitiation at an internal ORF is also observed in zebrafish. Interactome studies indicate that the human SEP53BP1 associates with components of the protein turnover pathway including the proteasome, and the TRiC/CCT chaperonin complex, suggesting that it may play a role in cellular proteostasis.

X-CAP improves pathogenicity prediction of stopgain variants

Stopgain substitutions are the third-largest class of monogenic human disease mutations and often examined first in patient exomes. Existing computational stopgain pathogenicity predictors, however, exhibit poor performance at the high sensitivity required for clinical use. Here, we introduce a new classifier, termed X-CAP, which uses a novel training methodology and unique feature set to improve the AUROC by 18% and decrease the false-positive rate 4-fold on large variant databases. In patient exomes, X-CAP prioritizes causal stopgains better than existing methods do, further illustrating its clinical utility. X-CAP is available at https://github.com/bejerano-lab/X-CAP .

KEPI plays a negative role in the repression that accompanies translational inhibition guided by the uORF element of human CHOP transcript during stress response

Translation of the downstream coding sequence of some mRNAs may be repressed by the upstream open reading frame (uORF) at their 5'-end. The mechanism underlying this uORF-mediated translational inhibition (uORF-MTI) is not fully understood in vivo. Recently, it was found that zebrafish Endouc or its human orthologue ENDOU (Endouc/ENDOU) plays a positive role in repressing the uORF-MTI of human CHOP (uORF^chop-MTI) during stress by blocking its activity However, the repression of uORF^chop-MTI assisted by an as-yet unidentified negative effector remains to be elucidated. Compared to the upregulated CHOP transcript, we herein report that the kepi (kinase-enhanced PP1 inhibitor) transcript was downregulated in the zebrafish embryos treated with both heat shock and hypoxia. Quantitative RT-PCR also revealed that the level of kepi mRNA was noticeably decreased in both heat-shock-treated and hypoxia-exposed embryos. When kepi mRNA was microinjected into the one-celled embryos from transgenic line huORFZ, the translation of downstream GFP reporter controlled by the uORF^chop-MTI was reduced in the hypoxia-exposed embryos. In contrast, when kepi was knocked down by injection of antisense Morpholino oligonucleotide, the translation of downstream GFP reporter was induced and expressed in the brain and spinal cord of injected embryos in the absence of stress. During normal condition, overexpression of KEPI increased eIF2α phosphorylation, resulting in inducing the translation of uORF-tag mRNA, such as ATF4 and CHOP mRNAs. However, during stress condition, overexpression of KEPI decreased eIF2α phosphorylation, resulting in reducing the GFP reporter and CHOP proteins. This is the first report to demonstrate that KEPI plays a negative role in uORF^chop - mediated translation during ER stress.

Decoupling expression and editing preferences of ADAR1 p150 and p110 isoforms

Human adenosine deaminase acting on RNA 1 (ADAR1) catalyzes adenosine-to-inosine deamination reactions on double-stranded RNA molecules to regulate cellular responses to endogenous and exogenous RNA. Defective ADAR1 editing leads to disorders such as Aicardi-Goutières syndrome, an autoinflammatory disease that manifests in the brain and skin, and dyschromatosis symmetrica hereditaria, a skin pigmentation disorder. Two ADAR1 protein isoforms, p150 (150 kDa) and p110 (110 kDa), are expressed and can edit RNA, but the contribution of each isoform to the editing landscape remains unclear, largely because of the challenges in expressing p150 without p110. In this study, we demonstrate that p110 is coexpressed with p150 from the canonical p150-encoding mRNA due to leaky ribosome scanning downstream of the p150 start codon. The presence of a strong Kozak consensus context surrounding the p110 start codon suggests the p150 mRNA is optimized to leak p110 alongside expression of p150. To reduce leaky scanning and translation initiation at the p110 start codon, we introduced synonymous mutations in the coding region between the p150 and p110 start codons. Cells expressing p150 constructs with these mutations produced significantly reduced levels of p110. Editing analysis of total RNA from ADAR1 knockout cells reconstituted separately with modified p150 and p110 revealed that more than half of the A-to-I edit sites are selectively edited by p150, and the other half are edited by either p150 or p110. This method of isoform-selective editing analysis, making use of the modified p150, has the potential to be adapted for other cellular contexts.

RUFY4 exists as two translationally regulated isoforms, that localize to the mitochondrion in activated macrophages

We report here that RUFY4, a newly characterized member of the 'RUN and FYVE domain-containing' family of proteins previously associated with autophagy enhancement, is highly expressed in alveolar macrophages (AM). We show that RUFY4 interacts with mitochondria upon stimulation by microbial-associated molecular patterns of AM and dendritic cells. RUFY4 interaction with mitochondria and other organelles is dependent on a previously uncharacterized OmpH domain located immediately upstream of its C-terminal FYVE domain. Further, we demonstrate that rufy4 messenger RNA can be translated from an alternative translation initiation codon, giving rise to a N-terminally truncated form of the molecule lacking most of its RUN domain and with enhanced potential for its interaction with mitochondria. Our observations point towards a role of RUFY4 in selective mitochondria clearance in activated phagocytes.

The C99 domain of the amyloid precursor protein resides in the disordered membrane phase

Processing of the amyloid precursor protein (APP) via the amyloidogenic pathway is associated with the etiology of Alzheimer's disease. The cleavage of APP by β-secretase to generate the transmembrane 99-residue C-terminal fragment (C99) and subsequent processing of C99 by γ-secretase to yield amyloid-β (Aβ) peptides are essential steps in this pathway. Biochemical evidence suggests that amyloidogenic processing of C99 occurs in cholesterol- and sphingolipid-enriched liquid-ordered phase membrane rafts. However, direct evidence that C99 preferentially associates with these rafts has remained elusive. Here, we tested this by quantifying the affinity of C99-GFP for raft domains in cell-derived giant plasma membrane vesicles (GPMVs). We found that C99 was essentially excluded from ordered domains in vesicles from HeLa cells, undifferentiated SH-SY5Y cells, or SH-SY5Y-derived neurons; instead, ∼90% of C99 partitioned into disordered domains. The strong association of C99 with disordered domains occurred independently of its cholesterol-binding activity or homodimerization, or of the presence of the familial Alzheimer disease Arctic mutation (APP E693G). Finally, through biochemical studies we confirmed previous results, which showed that C99 is processed in the plasma membrane by α-secretase, in addition to the well-known γ-secretase. These findings suggest that C99 itself lacks an intrinsic affinity for raft domains, implying that either i) amyloidogenic processing of the protein occurs in disordered regions of the membrane, ii) processing involves a marginal subpopulation of C99 found in rafts, or iii) as-yet-unidentified protein-protein interactions with C99 in living cells drive this protein into membrane rafts to promote its cleavage therein.

Hybrid Gene Origination Creates Human-Virus Chimeric Proteins during Infection

RNA viruses are a major human health threat. The life cycles of many highly pathogenic RNA viruses like influenza A virus (IAV) and Lassa virus depends on host mRNA, because viral polymerases cleave 5'-m7G-capped host transcripts to prime viral mRNA synthesis ("cap-snatching"). We hypothesized that start codons within cap-snatched host transcripts could generate chimeric human-viral mRNAs with coding potential. We report the existence of this mechanism of gene origination, which we named "start-snatching." Depending on the reading frame, start-snatching allows the translation of host and viral "untranslated regions" (UTRs) to create N-terminally extended viral proteins or entirely novel polypeptides by genetic overprinting. We show that both types of chimeric proteins are made in IAV-infected cells, generate T cell responses, and contribute to virulence. Our results indicate that during infection with IAV, and likely a multitude of other human, animal and plant viruses, a host-dependent mechanism allows the genesis of hybrid genes.

Translational profiling of B cells infected with the Epstein-Barr virus reveals 5' leader ribosome recruitment through upstream open reading frames

The Epstein-Barr virus (EBV) genome encodes several hundred transcripts. We have used ribosome profiling to characterize viral translation in infected cells and map new translation initiation sites. We show here that EBV transcripts are translated with highly variable efficiency, owing to variable transcription and translation rates, variable ribosome recruitment to the leader region and coverage by monosomes versus polysomes. Some transcripts were hardly translated, others mainly carried monosomes, showed ribosome accumulation in leader regions and most likely represent non-coding RNAs. A similar process was visible for a subset of lytic genes including the key transactivators BZLF1 and BRLF1 in cells infected with weakly replicating EBV strains. This suggests that ribosome trapping, particularly in the leader region, represents a new checkpoint for the repression of lytic replication. We could identify 25 upstream open reading frames (uORFs) located upstream of coding transcripts that displayed 5′ leader ribosome trapping, six of which were located in the leader region shared by many latent transcripts. These uORFs repressed viral translation and are likely to play an important role in the regulation of EBV translation.

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.

Networks of ribosome flow models for modeling and analyzing intracellular traffic

The ribosome flow model with input and output (RFMIO) is a deterministic dynamical system that has been used to study the flow of ribosomes during mRNA translation. The input of the RFMIO controls its initiation rate and the output represents the ribosome exit rate (and thus the protein production rate) at the 3′ end of the mRNA molecule. The RFMIO and its variants encapsulate important properties that are relevant to modeling ribosome flow such as the possible evolution of “traffic jams” and non-homogeneous elongation rates along the mRNA molecule, and can also be used for studying additional intracellular processes such as transcription, transport, and more. Here we consider networks of interconnected RFMIOs as a fundamental tool for modeling, analyzing and re-engineering the complex mechanisms of protein production. In these networks, the output of each RFMIO may be divided, using connection weights, between several inputs of other RFMIOs. We show that under quite general feedback connections the network has two important properties: (1) it admits a unique steady-state and every trajectory converges to this steady-state; and (2) the problem of how to determine the connection weights so that the network steady-state output is maximized is a convex optimization problem. These mathematical properties make these networks highly suitable as models of various phenomena: property (1) means that the behavior is predictable and ordered, and property (2) means that determining the optimal weights is numerically tractable even for large-scale networks. For the specific case of a feed-forward network of RFMIOs we prove an additional useful property, namely, that there exists a spectral representation for the network steady-state, and thus it can be determined without any numerical simulations of the dynamics. We describe the implications of these results to several fundamental biological phenomena and biotechnological objectives.

Genome-wide maps of ribosomal occupancy provide insights into adaptive evolution and regulatory roles of uORFs during Drosophila development

Upstream open reading frames (uORFs) play important roles in regulating the main coding DNA sequences (CDSs) via translational repression. Despite their prevalence in the genomes, uORFs are overall discriminated against by natural selection. However, it remains unclear why in the genomes there are so many uORFs more conserved than expected under the assumption of neutral evolution. Here, we generated genome-wide maps of translational efficiency (TE) at the codon level throughout the life cycle of Drosophila melanogaster. We identified 35,735 uORFs that were expressed, and 32,224 (90.2%) of them showed evidence of ribosome occupancy during Drosophila development. The ribosome occupancy of uORFs is determined by genomic features, such as optimized sequence contexts around their start codons, a shorter distance to CDSs, and higher coding potentials. Our population genomic analysis suggests the segregating mutations that create or disrupt uORFs are overall deleterious in D. melanogaster. However, we found for the first time that many (68.3% of) newly fixed uORFs that are associated with ribosomes in D. melanogaster are driven by positive Darwinian selection. Our findings also suggest that uORFs play a vital role in controlling the translational program in Drosophila. Moreover, we found that many uORFs are transcribed or translated in a developmental stage-, sex-, or tissue-specific manner, suggesting that selective transcription or translation of uORFs could potentially modulate the TE of the downstream CDSs during Drosophila development.

Upstream open reading frames (uORFs) in the 5′ untranslated regions (UTRs) of messenger RNAs can potentially inhibit translation of the downstream regions that encode proteins by sequestering protein-making machinery the ribosome. Moreover, mutations that destroy existing uORFs or create new ones are known to cause human disease. Although mutations that create new uORFs are generally deleterious and are selected against, many uORFs are evolutionarily conserved across eukaryotic species. To resolve this dilemma, we used extensive mRNA-Seq and ribosome profiling to generate high-resolution genome-wide maps of ribosome occupancy and translational efficiency (TE) during the life cycle of the fruit fly D. melanogaster. This allowed us to identify the sequence features of uORFs that influence their ability to associate with ribosomes. We demonstrate for the first time that the majority of the newly fixed uORFs in D. melanogaster, especially the translated ones, are under positive Darwinian selection. We also show that uORFs exert widespread repressive effects on the translation of the downstream protein-coding region. We find that many uORFs are transcribed or translated in a developmental stage-, sex-, or tissue-specific manner. Our results suggest that during Drosophila development, changes in the TE of uORFs, as well as the inclusion/exclusion of uORFs, are frequently exploited to inversely influence the translation of the downstream protein-coding regions. Our study provides novel insights into the molecular mechanisms and functional consequences of uORF-mediated regulation.

eIF1 modulates the recognition of suboptimal translation initiation sites and steers gene expression via uORFs

Alternative translation initiation mechanisms such as leaky scanning and reinitiation potentiate the polycistronic nature of human transcripts. By allowing for reprogrammed translation, these mechanisms can mediate biological responses to stimuli. We combined proteomics with ribosome profiling and mRNA sequencing to identify the biological targets of translation control triggered by the eukaryotic translation initiation factor 1 (eIF1), a protein implicated in the stringency of start codon selection. We quantified expression changes of over 4000 proteins and 10 000 actively translated transcripts, leading to the identification of 245 transcripts undergoing translational control mediated by upstream open reading frames (uORFs) upon eIF1 deprivation. Here, the stringency of start codon selection and preference for an optimal nucleotide context were largely diminished leading to translational upregulation of uORFs with suboptimal start. Interestingly, genes affected by eIF1 deprivation were implicated in energy production and sensing of metabolic stress.

Ribosome profiling reveals translational regulation of mammalian cells in response to hypoxic stress

Background

Retinal pigment epithelium (RPE) cells transfer oxygen and nutrients from choroid to the neural retina. Reduced oxygen to RPE perturbs development and functions of blood vessels in retina. Previous efforts of genome-wide studies have been largely focused on transcriptional changes of cells in response to hypoxia. Recently developed ribosome profiling provides an opportunity to study genome-wide translational changes. To gain systemic insights into the transcriptional and translational regulation of cellular in response to hypoxic stress, we used simultaneous RNA sequencing and ribosome profiling on an RPE cells line, ARPE-19, under hypoxia condition.

Results

Both HIF-1α and EPAS1 (HIF-2α) proteins were stabilized in ARPE-19 under hypoxic stress treatment at 1 h, 2 h and 4 h. Analysis of simultaneous RNA sequencing and ribosome profiling data showed genome-wide gene expression changes at both transcriptional and translational levels. Comparative analysis of ribosome profiling and RNA-seq data revealed that hypoxia induced changes of more genes at the translational than the transcriptional levels. Ribosomes densities at 5′ untranslated region (UTR) significantly increased under hypoxic stress. Interestingly, the increase in ribosome densities at 5′ UTR is positively correlated with the presence of upstream open reading frames (uORFs) in the 5′ UTR of mRNAs.

Conclusion

Our results characterized translational profiles of mRNAs for a RPE cell line in response to hypoxia. In particular, uORFs play important roles in the regulation of translation efficiency by affecting ribosomes loading onto mRNAs. This study provides the first attempt to understand translational response of mammalian cells under hypoxic condition.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3996-8) contains supplementary material, which is available to authorized users.

Translational effects and coding potential of an upstream open reading frame associated with DOPA Responsive Dystonia

Upstream open reading frames (uORFs) have emerged as major post-transcriptional regulatory elements in eukaryotic species. In general, uORFs are initiated by a translation start codon within the 5' untranslated region of a gene (upstream ATG; uATG), and they are negatively correlated with translational efficiency. In addition to their translational regulatory role, some uORFs can code for biologically active short peptides. The importance of uATGs/uORFs is further underscored by human diseases associated with single nucleotide polymorphisms (SNPs), which disrupt existing uORFs or introduce novel uORFs. Although several functional proteins translated from naturally occurring uORFs have been described, the coding potential of uORFs created by SNPs has been ignored because of the a priori assumption that these proteins are short-lived with no likely impact on protein homeostasis. Thus, studies on SNP-created uORFs are limited to their translational effects, leaving unexplored the potential cellular consequences of a SNP/uORF-encoded protein. Here, we investigate functionality of a uATG/uORF introduced by a +142C>T SNP within the GCH1 gene and associated with a familial form of DOPA Responsive Dystonia. We report that the +142C>T SNP represses GCH1 translation, and introduces a short, frame shifted uORF that encodes a 73-amino acid peptide. This peptide is localized within the nucleus and compromises cell viability upon proteasome inhibition. Our work extends the list of uATG/uORF associated diseases and advances research on peptides translated from SNP-introduced uORFs, a neglected component of the proteome.

Upstream ORFs are prevalent translational repressors in vertebrates

Regulation of gene expression is fundamental in establishing cellular diversity and a target of natural selection. Untranslated mRNA regions (UTRs) are key mediators of post-transcriptional regulation. Previous studies have predicted thousands of ORFs in 5'UTRs, the vast majority of which have unknown function. Here, we present a systematic analysis of the translation and function of upstream open reading frames (uORFs) across vertebrates. Using high-resolution ribosome footprinting, we find that (i)uORFs are prevalent within vertebrate transcriptomes, (ii) the majority show signatures of active translation, and (iii)uORFs act as potent regulators of translation and RNA levels, with a similar magnitude to miRNAs. Reporter experiments reveal clear repression of downstream translation by uORFs/oORFs. uORF number, intercistronic distance, overlap with the CDS, and initiation context most strongly influence translation. Evolution has targeted these features to favor uORFs amenable to regulation over constitutively repressive uORFs/oORFs. Finally, we observe that the regulatory potential of uORFs on individual genes is conserved across species. These results provide insight into the regulatory code within mRNA leader sequences and their capacity to modulate translation across vertebrates.

Sequence and Expression Analysis of Interferon Regulatory Factor 10 (IRF10) in Three Diverse Teleost Fish Reveals Its Role in Antiviral Defense

BACKGROUND

Interferon regulatory factor (IRF) 10 was first found in birds and is present in the genome of other tetrapods (but not humans and mice), as well as in teleost fish. The functional role of IRF10 in vertebrate immunity is relatively unknown compared to IRF1-9. The target of this research was to clone and characterize the IRF10 genes in three economically important fish species that will facilitate future evaluation of this molecule in fish innate and adaptive immunity.

MOLECULAR CHARACTERIZATION OF IRF10 IN THREE FISH SPECIES

In the present study, a single IRF10 gene was cloned in grass carp Ctenopharyngodon idella and Asian swamp eel Monopterus albus, and two, named IRF10a and IRF10b, in rainbow trout Oncorhynchus mykiss. The fish IRF10 molecules share highest identities to other vertebrate IRF10s, and have a well conserved DNA binding domain, IRF-associated domain, and an 8 exon/7 intron structure with conserved intron phase. The presence of an upstream ATG or open reading frame (ORF) in the 5'-untranslated region of different fish IRF10 cDNA sequences suggests potential regulation at the translational level, and this has been verified by in vitro transcription/translation experiments of the trout IRF10a cDNA, but would still need to be validated in fish cells.

EXPRESSION ANALYSIS OF IRF10 IN VIVO AND IN VITRO

Both trout IRF10 paralogues are highly expressed in thymus, blood and spleen but are relatively low in head kidney and caudal kidney. Trout IRF10b expression is significantly higher than IRF10a in integumentary tissues i.e. gills, scales, skin, intestine, adipose fin and tail fins, suggesting that IRF10b may be more important in mucosal immunity. The expression of both trout IRF10 paralogues is up-regulated by recombinant IFN-γ. The expression of the IRF10 genes is highly induced by Poly I:C in vitro and in vivo, and by viral infection, but is less responsive to peptidoglycan and bacterial infection, suggesting an important role of fish IRF10 in antiviral defense.

Identification and expression analysis of an atypical chemokine receptor-2 (ACKR2)/CC chemokine binding protein-2 (CCBP2) in rainbow trout (Oncorhynchus mykiss)

Atypical chemokine receptors (ACKRs) have emerged as key components of the chemokine system, with an essential regulatory function in innate and adaptive immune responses and inflammation. In mammals ACKR2 is a 'scavenging' receptor for inflammatory CC chemokines and plays a central role in the resolution of in vivo inflammatory responses. An ACKR2 like gene has been identified and cloned in rainbow trout (Teleostei) in the present study, enabling the further identification of this molecule in another group of ray-finned teleost fish (Holostei), in a lobe-finned fish (Sarcopterygii-coelacanth), and in reptiles. The identity of these ACKR2 molecules is supported by their conserved structure, and by phylogenetic tree and synteny analysis. Trout ACKR2 is highly expressed in spleen and head kidney, suggesting a homeostatic role of this receptor in limiting the availability of its potential ligands. Trout ACKR2 expression can be modulated in vivo by bacterial and parasitic infections, and in vitro by PAMPs (poly I:C and peptidoglycan) and cytokines (IL-6, TNF-α, IFN-γ and IL-21) in a time dependent manner. These patterns of expression and modulation suggest that trout ACKR2 is regulated in a complex way and has an important role in control of the chemokine network in fish as in mammals.

Identification and expression analysis of two interleukin-23α (p19) isoforms, in rainbow trout Oncorhynchus mykiss and Atlantic salmon Salmo salar

Interleukin (IL)-23 is a heterodimeric IL-12 family cytokine composed of a p19 α-chain, linked to a p40 β-chain that is shared with IL-12. IL-23 is distinguished functionally from IL-12 by its ability to induce the production of IL-17, and differentiation of Th17 cells in mammals. Three isoforms of p40 (p40a, p40b and p40c) have been found in some 3R teleosts. Salmonids also possess three p40 isoforms (p40b1, p40b2 and p40c) although p40a is missing, and two copies (paralogues) of p40b are present that have presumably been retained following the 4R duplication in this fish lineage. Teleost p19 has been discovered recently in zebrafish, but to date there is limited information on expression and modulation of this molecule. In this report we have cloned two p19 paralogues (p19a and p19b) in salmonids, suggesting that a salmonid can possess six potential IL-23 isoforms. Whilst Atlantic salmon has two active p19 genes, the rainbow trout p19b gene may have been pseudogenized. The salmonid p19 translations share moderate identities (22.8-29.9%) to zebrafish and mammalian p19 molecules, but their identity was supported by structural features, a conserved 4 exon/3 intron gene organisation, and phylogenetic tree analysis. The active salmonid p19 genes are highly expressed in blood and gonad. Bacterial (Yersinia ruckeri) and viral infection in rainbow trout induces the expression of p19a, suggesting pathogen-specific induction of IL-23 isoforms. Trout p19a expression was also induced by PAMPs (poly IC and peptidoglycan) and the proinflammatory cytokine IL-1β in primary head kidney macrophages. These data may indicate diverse functional roles of trout IL-23 isoforms in regulating the immune response in fish.

Identification and characterisation of the IL-27 p28 subunits in fish: Cloning and comparative expression analysis of two p28 paralogues in Atlantic salmon Salmo salar

Interleukin (IL)-27 is an IL-6/IL-12 family member with pro-inflammatory and anti-inflammatory properties. It is a heterodimeric cytokine composed of an α-chain p28 and a β-chain Ebi3 (Epstein-Barr virus induce gene 3). The p28 subunit can also be secreted as a monomer and function as IL-30 that acts as an inhibitor of IL-27 signalling. At present, the p28 gene has only been described in mammals. Thus, for the first time outwith mammals, we have identified seven p28 molecules in six divergent teleost fish species, Atlantic salmon, two cichlids, two cyprinids and a yellowtail. The fish p28 molecules have higher similarities to mammalian p28 than other IL-6/12 family members. Critical residues involved in the interaction with Ebi3 and the receptor gp130 are highly conserved. The prediction that these are true orthologues is supported by phylogenetic and synteny analysis. Two p28 paralogues (p28a and p28b) sharing 72% aa identity have been identified and characterised in Atlantic salmon. There are multiple upstream ATGs in the 5'-UTR and ATTTA motifs in the 3'-UTR of both cDNA sequences, suggesting regulation at the post-transcriptional and translational level. Both salmon p28 genes are highly expressed in immune relevant tissues, such as thymus, gills, spleen and head kidney. Conversely salmon Ebi3 is highly expressed in other organs, such as liver and caudal kidney. The expression of p28 but not Ebi3 was induced by PAMPs and recombinant cytokines in head kidney cells, and in spleen by Poly I:C challenge in vivo. The dissociation of the expression and modulation of p28 and Ebi3 suggest that both p28 and Ebi3 may be secreted alone or with other partners.

N-terminal proteomics and ribosome profiling provide a comprehensive view of the alternative translation initiation landscape in mice and men

Usage of presumed 5'UTR or downstream in-frame AUG codons, next to non-AUG codons as translation start codons contributes to the diversity of a proteome as protein isoforms harboring different N-terminal extensions or truncations can serve different functions. Recent ribosome profiling data revealed a highly underestimated occurrence of database nonannotated, and thus alternative translation initiation sites (aTIS), at the mRNA level. N-terminomics data in addition showed that in higher eukaryotes around 20% of all identified protein N termini point to such aTIS, to incorrect assignments of the translation start codon, translation initiation at near-cognate start codons, or to alternative splicing. We here report on more than 1700 unique alternative protein N termini identified at the proteome level in human and murine cellular proteomes. Customized databases, created using the translation initiation mapping obtained from ribosome profiling data, additionally demonstrate the use of initiator methionine decoded near-cognate start codons besides the existence of N-terminal extended protein variants at the level of the proteome. Various newly identified aTIS were confirmed by mutagenesis, and meta-analyses demonstrated that aTIS reside in strong Kozak-like motifs and are conserved among eukaryotes, hinting to a possible biological impact. Finally, TargetP analysis predicted that the usage of aTIS often results in altered subcellular localization patterns, providing a mechanism for functional diversification.

Differential expression, modulation and bioactivity of distinct fish IL-12 isoforms: implication towards the evolution of Th1-like immune responses

IL-12 is a heterodimeric cytokine composed of an α-chain (p35) and a β-chain (p40). Primarily produced by APCs, IL-12 induces IFN-γ production in T, B and NK cells. IL-12 drives Th1-cell differentiation and IFN-γ secretion to promote cell-mediated immunity, which is essential in the defence against intracellular pathogens. The importance of IL-12 in Th1 responses is echoed by its targeted suppression by intracellular pathogens evading cell-mediated immunity. IL-12 subunits have been identified recently in fish, although reported bioactivities are limited to higher vertebrates. Here, we report the cloning of a p35 gene and two divergent p40 genes (p40b and p40c), capable of producing two functional IL-12 isoforms (p35/p40b and p35/p40c) in rainbow trout. Trout IL-12 isoforms possess distinct bioactivities with respect to the induction of IFN-γ, IL-10 and p40c expression. Trout IL-12 isoforms are differentially expressed and modulated in vivo, exhibiting specific gene expression profiles in bacterial, viral and parasitic infection models, and in vitro in stimulated macrophage and leucocyte cultures. These data imply that alternative or additional pathogen-specific Th-like cell populations may exist in fish. This study will facilitate a broader understanding of the evolutionary processes driving host-pathogen interactions and Th1-like immune responses in lower vertebrates.

Gene expression regulation by upstream open reading frames and human disease

Upstream open reading frames (uORFs) are major gene expression regulatory elements. In many eukaryotic mRNAs, one or more uORFs precede the initiation codon of the main coding region. Indeed, several studies have revealed that almost half of human transcripts present uORFs. Very interesting examples have shown that these uORFs can impact gene expression of the downstream main ORF by triggering mRNA decay or by regulating translation. Also, evidence from recent genetic and bioinformatic studies implicates disturbed uORF-mediated translational control in the etiology of many human diseases, including malignancies, metabolic or neurologic disorders, and inherited syndromes. In this review, we will briefly present the mechanisms through which uORFs regulate gene expression and how they can impact on the organism's response to different cell stress conditions. Then, we will emphasize the importance of these structures by illustrating, with specific examples, how disturbed uORF-mediated translational control can be involved in the etiology of human diseases, giving special importance to genotype-phenotype correlations. Identifying and studying more cases of uORF-altering mutations will help us to understand and establish genotype-phenotype associations, leading to advancements in diagnosis, prognosis, and treatment of many human disorders.

New universal rules of eukaryotic translation initiation fidelity

The accepted model of eukaryotic translation initiation begins with the scanning of the transcript by the pre-initiation complex from the 5′end until an ATG codon with a specific nucleotide (nt) context surrounding it is recognized (Kozak rule). According to this model, ATG codons upstream to the beginning of the ORF should affect translation. We perform for the first time, a genome-wide statistical analysis, uncovering a new, more comprehensive and quantitative, set of initiation rules for improving the cost of translation and its efficiency. Analyzing dozens of eukaryotic genomes, we find that in all frames there is a universal trend of selection for low numbers of ATG codons; specifically, 16–27 codons upstream, but also 5–11 codons downstream of the START ATG, include less ATG codons than expected. We further suggest that there is selection for anti optimal ATG contexts in the vicinity of the START ATG. Thus, the efficiency and fidelity of translation initiation is encoded in the 5′UTR as required by the scanning model, but also at the beginning of the ORF.

The observed nt patterns suggest that in all the analyzed organisms the pre-initiation complex often misses the START ATG of the ORF, and may start translation from an alternative initiation start-site. Thus, to prevent the translation of undesired proteins, there is selection for nucleotide sequences with low affinity to the pre-initiation complex near the beginning of the ORF. With the new suggested rules we were able to obtain a twice higher correlation with ribosomal density and protein levels in comparison to the Kozak rule alone (e.g. for protein levels r = 0.7 vs. r = 0.31; p<10⁻¹²).

Gene translation is an important step of the intra-cellular protein synthesis, which is a central process in all living organisms. Thus, understanding how translation efficiency is encoded in transcripts has ramifications to every biomedical discipline. The aim of the current study is to decipher the way translation initiation fidelity is encoded in eukaryotic transcripts, and how evolution shapes the beginning of transcripts. Based on the genomes of dozens of organisms we were able to derive a new, more precise, set of rules related to this process, facilitating a high resolution view of the mechanisms aiding translation initiation fidelity. Among others, we show that there is a universal trend of selection for low numbers of ATG codons upstream, but also in the 5–11 codons downstream of the START ATG, presumably to prevent translation of alternative ORFs over the main one. With the new suggested rules we were able to obtain a twice higher correlation with ribosomal density and protein levels in comparison to the previous translation initiation efficiency rule.

Nonsense-mediated mRNA decay of collagen -emerging complexity in RNA surveillance mechanisms

Nonsense-mediated mRNA decay (NMD) is an evolutionarily conserved mRNA surveillance system that degrades mRNA transcripts that harbour a premature translation-termination codon (PTC), thus reducing the synthesis of truncated proteins that would otherwise have deleterious effects. Although extensive research has identified a conserved repertoire of NMD factors, these studies have been performed with a restricted set of genes and gene constructs with relatively few exons. As a consequence, NMD mechanisms are poorly understood for genes with large 3' terminal exons, and the applicability of the current models to large multi-exon genes is not clear. In this Commentary, we present an overview of the current understanding of NMD and discuss how analysis of nonsense mutations in the collagen gene family has provided new mechanistic insights into this process. Although NMD of the collagen genes with numerous small exons is consistent with the widely accepted exon-junction complex (EJC)-dependent model, the degradation of Col10a1 transcripts with nonsense mutations cannot be explained by any of the current NMD models. Col10a1 NMD might represent a fail-safe mechanism for genes that have large 3' terminal exons. Defining the mechanistic complexity of NMD is important to allow us to understand the pathophysiology of the numerous genetic disorders caused by PTC mutations.

5' UTR control of native ERG and of Tmprss2:ERG variants activity in prostate cancer

ERG, a member of the ETS transcription factor family, is frequently overexpressed in prostate cancer as a result of its fusion to the androgen-responsive Tmprss2 gene. Different genomic rearrangements and alternative splicing events around the junction region lead to multiple combination of Tmprss2:ERG fusion transcripts that correlate with different tumor aggressiveness, but their specific functions and biological activities are still unclear. The complexity of ERG expression pattern is compounded by the use of alternative promoters, splice sites, polyadenylation sites and translation initiation sites in both the native and fusion contexts. Our systematic characterization of native ERG and Tmprss2:ERG variants reveals that their different oncogenic potential is impacted by the status of the Ets domain and the configuration of the 5′ UTR region. In particular, expression and activity of functional ERG and Tmprss2:ERG variants are influenced both by translation initiation signals within the different isoforms and by inhibitory upstream Open Reading Frames (uORF) in their 5′ UTRs. Stable expression of ERG and Tmprss2:ERG variants promoted cell migration/invasion, induced a block of proliferation and induced a senescence-like state, suggesting a role for these variants in the prostate tumorigenesis process. In addition to Tmprss2:ERG fusion products, a group of related native ERG isoforms is also highly over-expressed in fusion-carrying prostate cancers, and share the same translation initiation site (in ERG exon 4) with the commonly observed Tmprss2 exon1 joined to ERG exon 4 (T1:E4) fusion-derived variant. Usage of this ATG can be preferentially down-regulated by directed antisense-based compounds, possibly representing the basis of a targeted approach that distinguishes between tumor–associated and normal ERG.

Whole human genome proteogenomic mapping for ENCODE cell line data: identifying protein-coding regions

BACKGROUND

Proteogenomic mapping is an approach that uses mass spectrometry data from proteins to directly map protein-coding genes and could aid in locating translational regions in the human genome. In concert with the ENcyclopedia of DNA Elements (ENCODE) project, we applied proteogenomic mapping to produce proteogenomic tracks for the UCSC Genome Browser, to explore which putative translational regions may be missing from the human genome.

RESULTS

We generated ~1 million high-resolution tandem mass (MS/MS) spectra for Tier 1 ENCODE cell lines K562 and GM12878 and mapped them against the UCSC hg19 human genome, and the GENCODE V7 annotated protein and transcript sets. We then compared the results from the three searches to identify the best-matching peptide for each MS/MS spectrum, thereby increasing the confidence of the putative new protein-coding regions found via the whole genome search. At a 1% false discovery rate, we identified 26,472, 24,406, and 13,128 peptides from the protein, transcript, and whole genome searches, respectively; of these, 481 were found solely via the whole genome search. The proteogenomic mapping data are available on the UCSC Genome Browser at http://genome.ucsc.edu/cgi-bin/hgTrackUi?db=hg19&g=wgEncodeUncBsuProt.

CONCLUSIONS

The whole genome search revealed that ~4% of the uniquely mapping identified peptides were located outside GENCODE V7 annotated exons. The comparison of the results from the disparate searches also identified 15% more spectra than would have been found solely from a protein database search. Therefore, whole genome proteogenomic mapping is a complementary method for genome annotation when performed in conjunction with other searches.

An upstream open reading frame modulates ebola virus polymerase translation and virus replication

Ebolaviruses, highly lethal zoonotic pathogens, possess longer genomes than most other non-segmented negative-strand RNA viruses due in part to long 5′ and 3′ untranslated regions (UTRs) present in the seven viral transcriptional units. To date, specific functions have not been assigned to these UTRs. With reporter assays, we demonstrated that the Zaire ebolavirus (EBOV) 5′-UTRs lack internal ribosomal entry site function. However, the 5′-UTRs do differentially regulate cap-dependent translation when placed upstream of a GFP reporter gene. Most dramatically, the 5′-UTR derived from the viral polymerase (L) mRNA strongly suppressed translation of GFP compared to a β-actin 5′-UTR. The L 5′-UTR is one of four viral genes to possess upstream AUGs (uAUGs), and ablation of each uAUG enhanced translation of the primary ORF (pORF), most dramatically in the case of the L 5′-UTR. The L uAUG was sufficient to initiate translation, is surrounded by a “weak” Kozak sequence and suppressed pORF translation in a position-dependent manner. Under conditions where eIF2α was phosphorylated, the presence of the uORF maintained translation of the L pORF, indicating that the uORF modulates L translation in response to cellular stress. To directly address the role of the L uAUG in virus replication, a recombinant EBOV was generated in which the L uAUG was mutated to UCG. Strikingly, mutating two nucleotides outside of previously-defined protein coding and cis-acting regulatory sequences attenuated virus growth to titers 10–100-fold lower than a wild-type virus in Vero and A549 cells. The mutant virus also exhibited decreased viral RNA synthesis as early as 6 hours post-infection and enhanced sensitivity to the stress inducer thapsigargin. Cumulatively, these data identify novel mechanisms by which EBOV regulates its polymerase expression, demonstrate their relevance to virus replication and identify a potential therapeutic target.

Filoviruses (Ebola and Marburg viruses) are emerging zoonotic pathogens that cause lethal hemorrhagic fever in humans and have the potential to be employed as bioterrorism agents. Currently, approved therapeutics to treat filovirus infections are not available and new treatment strategies could be facilitated by improved mechanistic insight into the virus replication cycle. Compared to other related viruses, filovirus messenger RNAs have unusually long 5′ untranslated regions (UTRs) with undefined functions. In the Zaire ebolavirus (EBOV) genome, four of its seven messenger RNAs have 5′-UTRs with a small upstream open reading frame (uORF). We found that a uORF present in the EBOV polymerase (L) 5′-UTR suppresses L protein production and established a reporter assay to demonstrate that this uORF maintains L translation following the induction of an innate immune response; a phenomenon observed with several uORF-containing cellular messenger RNAs. The presence of the uORF is important for optimal virus replication, because a mutant virus lacking the upstream reading frame replicates less efficiently than a wildtype virus, an attenuation which is more pronounced following the induction of cellular stress. These studies define a novel mechanism by which filovirus upstream open reading frames modulate virus protein translation in the face of an innate immune response and highlight their importance in filovirus replication.

Selective constraint on the upstream open reading frames that overlap with coding sequences in animals

Upstream open reading frames (uORFs) are translational regulatory elements located in 5′ untranslated regions. They can significantly repress the translation of the downstream coding sequences (CDS), and participate in the spatio-temporal regulations of protein translation. Notwithstanding this biological significance, the selective constraint on uORFs remains underexplored. Particularly, the uORFs that partially overlap with CDS with a different reading frame (overlapping uORFs, or “VuORFs”) may lead to strong translational inhibition or N-terminal truncation of the peptides encoded by the affected CDS. By analyzing VuORF-containing transcripts (designated as “VuORF transcripts”) in human, mouse, and fruit fly, we demonstrate that VuORFs are in general slightly deleterious - the proportion of genes that encode at least one VuORF transcript is significantly smaller than expected in all of the three examined species. In addition, this proportion is significantly smaller in fruit fly than in mammals, indicating a higher efficiency of removing VuORFs in the former species because of its larger effective population size. Furthermore, the deleterious effect of a VuORF depends on the sequence context of its start codon (VuAUG). VuORFs with an optimal VuAUG context are more strongly disfavored than those with a suboptimal context in all of the three examined species. And the propensity to remove optimal-context VuAUGs is stronger in fruit fly than in mammals. Intriguingly, however, the currently observable optimal-context VuAUGs (but not suboptimal-context VuAUGs) are more conserved than expected. These observations suggest that the regulatory functions of VuORFs may have been gained fortuitously in organisms with a small effective population size because the slightly deleterious effect of these elements can be better tolerated in these organisms, thus allowing opportunities for the development of novel biological functions. Nevertheless, once the functions of VuORFs were established, they became subject to negative selection.

Possible link between the synthesis of GR alpha isoforms and eIF2 alpha phosphorylation state

Glucocorticoid hormones regulate numerous physiological processes and are widely used in the treatment of inflammation, autoimmune disease and cancer. Glucocorticoid receptor (GR) - a transcription factor, derived from a single gene, is responsible for the diverse actions of glucocorticoids. It was shown that GR gene gives rise a variety of mRNA species that produces several protein isoforms, among them GRα is the most abundant. In addition, GRα N-end-truncated protein isoforms (A, B, C, D) are generated by translational mechanisms. As it was found that the ratio between the translational isoforms amounts varied in different tissues and cell lines and distinct isoforms could control transcription of different sets of genes, molecular mechanisms underlining the synthesis of translational GRα isoforms are of great interest. It was considered that GRα isoform A is translated by a conventional linear scanning, isoform B is translated by leaky scanning, isoform C is translated by leaky scanning and ribosomal shunt whereas translation of isoform D occurs through ribosomal shunt only. Since the sequence organization of GRα mRNA strongly resembles the cases of ATF4 or ATF5, the well-known examples of reinitiation-dependent synthesis of functional isoforms, we hypothesize that translation of isoform C could be controlled by reinitiation mechanism also. If this assumption is correct, the ratio between GRα N-end isoforms could depend on the eIF2α phosphorylation state that could provide an additional connection between the GR and cellular stresses. We believe that this hypothesis could be of interest to plan more robust experiments or for better interpretation of available data.

Hidden coding potential of eukaryotic genomes: nonAUG started ORFs

It is widely considered that the vast majority of eukaryotic mRNAs contain only one open reading frame (ORF) and encode single protein. However, eukaryotic ribosomes can initiate translation at alternative start codons due to leaky scanning or reinitiation mechanisms that provides an opportunity to synthesize several protein products. Recent investigations also demonstrated that alternative translation from nonAUG start codons and AUG codons in a weak nucleotide context could make an important contribution to eukaryotic proteomes. However, accurate prediction of alternative start codons demands detailed investigation of mRNA features influencing their recognition by eukaryotic ribosomes. In this work, we present the results of computational analysis of characteristics of yeast and mammalian mRNAs potentially involved in the recognition of nonAUG start codons. It was found that sequence features of nonAUG started Saccharomyces cerevisiae upstream ORFs (uORFs) were adjusted for efficient translation and these uORFs could frequently encode functional polypeptides. In particular, our initial studies revealed that predicted tertiary structures downstream of nonAUG start sites in mammalian mRNAs were energetically more stable than those predicted for AUG start sites with strong Kozak context. We hypothesize that presence of such stable tertiary structure downstream of nonAUG start sites could be an important factor for the ribosome to recognize noncanonical start codons.

Extensive translatome remodeling during ER stress response in mammalian cells

In this work we have described the translatome of two mammalian cell lines, NIH3T3 and Jurkat, by scoring the relative polysome association of ∼10,000 mRNA under normal and ER stress conditions. We have found that translation efficiencies of mRNA correlated poorly with transcript abundance, although a general tendency was observed so that the highest translation efficiencies were found in abundant mRNA. Despite the differences found between mouse (NIH3T3) and human (Jurkat) cells, both cell types share a common translatome composed by ∼800–900 mRNA that encode proteins involved in basic cellular functions. Upon stress, an extensive remodeling in translatomes was observed so that translation of ∼50% of mRNA was inhibited in both cell types, this effect being more dramatic for those mRNA that accounted for most of the cell translation. Interestingly, we found two subsets comprising 1000–1500 mRNA whose translation resisted or was induced by stress. Translation arrest resistant class includes many mRNA encoding aminoacyl tRNA synthetases, ATPases and enzymes involved in DNA replication and stress response such as BiP. This class of mRNA is characterized by high translation rates in both control and stress conditions. Translation inducible class includes mRNA whose translation was relieved after stress, showing a high enrichment in early response transcription factors of bZIP and zinc finger C2H2 classes. Unlike yeast, a general coordination between changes in translation and transcription upon stress (potentiation) was not observed in mammalian cells. Among the different features of mRNA analyzed, we found a relevant association of translation efficiency with the presence of upstream ATG in the 5′UTR and with the length of coding sequence of mRNA, and a looser association with other parameters such as the length and the G+C content of 5′UTR. A model for translatome remodeling during the acute phase of stress response in mammalian cells is proposed.

Alternative splicing within the elk-1 5' untranslated region serves to modulate initiation events downstream of the highly conserved upstream open reading frame 2

The 5' untranslated region (UTR) plays a central role in the regulation of mammalian translation initiation. Key components include RNA structure, upstream AUGs (uAUGs), upstream open reading frames (uORFs), and internal ribosome entry site elements that can interact to modulate the readout. We previously reported the characterization of two alternatively spliced 5' UTR isoforms of the human elk-1 gene. Both contain two uAUGs and a stable RNA stem-loop, but the long form (5' UTR(L)) was more repressive than the short form (5' UTR(S)) for initiation at the ELK-1 AUG. We now demonstrate that ELK-1 expression arises by a combination of leaky scanning and reinitiation, with the latter mediated by the small uORF2 conserved in both spliced isoforms. In HEK293T cells, a considerable fraction of ribosomes scans beyond the ELK-1 AUG in a reinitiation mode. These are sequestered by a series of out-of-frame AUG codons that serve to prevent access to a second in-frame AUG start site used to express short ELK-1 (sELK-1), an N-terminally truncated form of ELK-1 that has been observed only in neuronal cells. We present evidence that all these events are fine-tuned by the nature of the 5' UTR and the activity of the α subunit of eukaryotic initiation factor 2 and provide insights into the neuronal specificity of sELK-1 expression.

The genomic features that affect the lengths of 5' untranslated regions in multicellular eukaryotes

Background

The lengths of 5’UTRs of multicellular eukaryotes have been suggested to be subject to stochastic changes, with upstream start codons (uAUGs) as the major constraint to suppress 5’UTR elongation. However, this stochastic model cannot fully explain the variations in 5’UTR length. We hypothesize that the selection pressure on a combination of genomic features is also important for 5’UTR evolution. The ignorance of these features may have limited the explanatory power of the stochastic model. Furthermore, different selective constraints between vertebrates and invertebrates may lead to differences in the determinants of 5’UTR length, which have not been systematically analyzed.

Methods

Here we use a multiple linear regression model to delineate the correlation between 5’UTR length and the combination of a series of genomic features (G+C content, observed-to-expected (OE) ratios of uAUGs, upstream stop codons (uSTOPs), methylation-related CG/UG dinucleotides, and mRNA-destabilizing UU/UA dinucleotides) in six vertebrates (human, mouse, rat, chicken, African clawed frog, and zebrafish) and four invertebrates (fruit fly, mosquito, sea squirt, and nematode). The relative contributions of each feature to the variation of 5’UTR length were also evaluated.

Results

We found that 14%~33% of the 5’UTR length variations can be explained by a linear combination of the analyzed genomic features. The most important genomic features are the OE ratios of uSTOPs and G+C content. The surprisingly large weightings of uSTOPs highlight the importance of selection on upstream open reading frames (which include both uAUGs and uSTOPs), rather than on uAUGs per se. Furthermore, G+C content is the most important determinants for most invertebrates, but for vertebrates its effect is second to uSTOPs. We also found that shorter 5’UTRs are affected more by the stochastic process, whereas longer 5’UTRs are affected more by selection pressure on genomic features.

Conclusions

Our results suggest that upstream open reading frames may be the real target of selection, rather than uAUGs. We also show that the selective constraints on genomic features of 5’UTRs differ between vertebrates and invertebrates, and between longer and shorter 5’UTRs. A more comprehensive model that takes these findings into consideration is needed to better explain 5’UTR length evolution.

Translational regulation mechanisms of aquaporin-4 supramolecular organization in astrocytes

The two predominant isoforms of Aquaporin-4 (AQP4), AQP4-M23 and AQP4-M1, assemble in the plasma membrane to form supramolecular structures called Orthogonal Array of Particles (OAPs) whose dimension is tightly associated to the M1/M23 ratio. Here, we explore translational regulation contribution to M1/M23 expression in primary cultures of rat astrocytes, and analyze the role of M1 mRNA 5'untranslated region (5'UTR) in this mechanism. Using isoform-specific RNAi we found that in rat astrocytes primary cultures a large proportion of M23 protein derives from M1 mRNA translation. Furthermore, site-specific mutagenesis of the 5'UTR sequence of AQP4-M1 mRNA indicates that a multiple-site leaky scanning mechanism, an out-of-frame upstream ORF (uORF), and a reinitiation mechanism are able to modulate the M1/M23 ratio and consequently, OAPs formation. These mechanisms are likely to be shared by different species, including human, and they can also be assumed to play a role in those pathophysiological situations where the organization of AQP4 in supramolecular structures (OAPs) is involved. Finally, we report that, when transfected in Hela cells, the longer rat AQP4 isoform, called Mz, which is not present in human impairs OAPs formation.

Secondary structure and the role in translation initiation of the 5'-terminal region of p53 mRNA

The p53 protein is one of the major factors involved in cell cycle control, DNA repair, and induction of apoptosis. We determined the secondary structure of the 5'-terminal region of p53 mRNA that includes two major translation initiation codons AUG1 and AUG2, responsible for the synthesis of p53 and its N-truncated isoform ΔN-p53. It turned out that a part of the coding sequence was involved in the folding of the 5' untranslated region for p53. The most characteristic structural elements in the 5'-terminal region of p53 mRNA were two hairpin motifs. In one of them, the initiation codon AUG1 was embedded while the other hairpin has been earlier shown to bind the Mdm2 protein. Alternative mechanisms of p53 mRNA translation initiation were investigated in vitro using model mRNA templates. The results confirmed that initiation from AUG1 was mostly cap-dependent. The process was stimulated by a cap structure and strongly inhibited by a stable hairpin at the template 5' end. Upon inhibition, the remaining protein fraction was synthesized in a cap-independent process, which was strongly stimulated by the addition of a cap analogue. The translation initiation from AUG2 showed a largely cap-independent character. The 5' cap structure actually decreased initiation from this site which argues against a leaky scanning mechanism but might suggest the presence of an IRES. Moreover, blocking cap-dependent translation from AUG1 by the stable hairpin did not change the level of initiation from AUG2. Upon addition of the cap analogue, translation from this site was even increased.

Alternative splicing of the human gene SYBL1 modulates protein domain architecture of Longin VAMP7/TI-VAMP, showing both non-SNARE and synaptobrevin-like isoforms

Background

The control of intracellular vesicle trafficking is an ideal target to weigh the role of alternative splicing in shaping genomes to make cells. Alternative splicing has been reported for several Soluble N-ethylmaleimide-sensitive factor Attachment protein REceptors of the vesicle (v-SNAREs) or of the target membrane (t-SNARES), which are crucial to intracellular membrane fusion and protein and lipid traffic in Eukaryotes. However, splicing has not yet been investigated in Longins, i.e. the most widespread v-SNAREs. Longins are essential in Eukaryotes and prototyped by VAMP7, Sec22b and Ykt6, sharing a conserved N-terminal Longin domain which regulates membrane fusion and subcellular targeting. Human VAMP7/TI-VAMP, encoded by gene SYBL1, is involved in multiple cell pathways, including control of neurite outgrowth.

Results

Alternative splicing of SYBL1 by exon skipping events results in the production of a number of VAMP7 isoforms. In-frame or frameshift coding sequence modifications modulate domain architecture of VAMP7 isoforms, which can lack whole domains or domain fragments and show variant or extra domains. Intriguingly, two main types of VAMP7 isoforms either share the inhibitory Longin domain and lack the fusion-promoting SNARE motif, or vice versa. Expression analysis in different tissues and cell lines, quantitative real time RT-PCR and confocal microscopy analysis of fluorescent protein-tagged isoforms demonstrate that VAMP7 variants have different tissue specificities and subcellular localizations. Moreover, design and use of isoform-specific antibodies provided preliminary evidence for the existence of splice variants at the protein level.

Conclusions

Previous evidence on VAMP7 suggests inhibitory functions for the Longin domain and fusion/growth promoting activity for the Δ-longin molecule. Thus, non-SNARE isoforms with Longin domain and non-longin SNARE isoforms might have somehow opposite regulatory functions. When considering splice variants as "natural mutants", evidence on modulation of subcellular localization by variation in domain combination can shed further light on targeting determinants. Although further work will be needed to characterize identified variants, our data might open the route to unravel novel molecular partners and mechanisms, accounting for the multiplicity of functions carried out by the different members of the Longin proteins family.

Tandem termination signal in plant mRNAs

It was proposed that if some mRNA characteristics resulted in a low efficiency of termination signal, an additional closely located stop codon (tandem stop codons) could be used to prevent the harmful readthrough. However, the role of tandem terminators in higher eukaryotes was not verified and remains hypothetical. In this work the sequence features of Arabidopsis thaliana and Oryza sativa mRNAs were analyzed. It was found that plant mRNAs with UGA terminator were characterized by a higher frequency of nonsense codons in the first triplet position of 3'-UTR that could result from a weak natural selection for "reserve" stop signal. Interestingly, the presence of tandem stop codons positively correlated with a specific amino acid composition in the C-terminal position of the encoded proteins. In particular, C-terminal glycine positively correlated with significantly higher frequencies of reserve terminators at the beginning positions of 3'-UTR in UGA-containing mRNAs. This finding coincides with some earlier observations concerning the role of glycine and its codons in inefficient termination of translation and recoding (e.g., 2A oligopeptide).

Ribosomal scanning on the 5'-untranslated region of the human immunodeficiency virus RNA genome

Translation initiation on most eukaryotic mRNAs occurs via a cap-dependent scanning mechanism and its efficiency is modulated by their 5'-untranslated regions (5'-UTR). The human immunodeficiency virus type 1 (HIV-1) 5'-UTR contains a stable TAR hairpin directly at its 5'-end, which possibly masks the cap structure. In addition, the 5'-UTR is relatively long and contains several stable RNA structures that are essential for viral replication. These characteristics may interfere with ribosomal scanning and suggest that translation is initiated via internal entry of ribosomes. Literature on the HIV-1 5'-UTR-driven translation initiation mechanism is controversial. Both scanning and internal initiation have been shown to occur in various experimental systems. To gain further insight in the translation initiation process, we determined which part of the 5'-UTR is scanned. To do so, we introduced upstream AUGs at various positions across the 5'-UTR and determined the effect on expression of a downstream reporter gene that was placed under control of the gag start codon. This strategy allowed us to determine the window of ribosomal scanning on the HIV-1 5'-UTR.

Stress-inducible alternative translation initiation of human cytomegalovirus latency protein pUL138

We have previously characterized a 21-kDa protein encoded by UL138 (pUL138) as a viral factor inherent to low-passage strains of human cytomegalovirus (HCMV) that is required for latent infection in vitro. pUL138 is encoded on 3.6-, 2.7-, and 1.4-kb 3' coterminal transcripts that are produced during productive and latent infections. pUL138 is encoded at the 3' end of each transcript and is preceded by an extensive 5' sequence (approximately 0.5 to 2.5 kb) containing several putative open reading frames (ORFs). We determined that three putative ORFs upstream of UL138 (UL133, UL135, and UL136) encode proteins. The UL138 transcripts are polycistronic, such that each transcript expresses pUL138 in addition to the most-5' ORF. The upstream coding sequences (CDS) present a significant challenge for the translation of pUL138 in mammalian cells. We hypothesized that sequences 5' of UL138 mediate translation initiation of pUL138 by alternative strategies. Accordingly, a 663-nucloetide (nt) sequence overlapping the UL136 CDS supported expression of a downstream cistron in a bicistronic reporter system. We did not detect cryptic promoter activity or RNA splicing events that could account for downstream cistron expression. These data are consistent with the sequence element functioning as an internal ribosome entry site (IRES). Interestingly, pUL138 expression from the 3.6- and 2.7-kb transcripts was induced by serum stress, which concomitantly inhibited normal cap-dependent translation. Our work suggests that an alternative and stress-inducible strategy of translation initiation ensures expression of pUL138 under a variety of cellular contexts. The UL138 polycistronic transcripts serve to coordinate the expression of multiple proteins, including a viral determinant of HCMV latency.

Interrelations between the nucleotide context of human start AUG codon, N-end amino acids of the encoded protein and initiation of translation

It is known that the recognition of AUG triplet by eukaryotic ribosomes as a translation start site strongly depends on its nucleotide context. However, the relative significance of different context positions is not fully clear. In particular, it concerns the role of 3'-end part of the context located at the beginning of the protein-coding sequence. The significant bias observed in nucleotide frequencies in positions +4, +5, +6 (corresponding to the second codon of CDS) could result from different reasons and their contribution to start codon recognition and initiation of translation is under discussion. In this study, we conducted a comparative computational analysis of the human mRNA samples containing different nucleotides (adenine, guanine or pyrimidine) in the essential context position -3. It was found that the presence of G in position +4 could be important for the context variant GnnAUG but not for AnnAUG. Interestingly, the second position of proteins encoded by mRNAs with AnnAUG context variant was specifically and significantly enriched with serine whereas the presence of GnnAUG context also correlated with a higher occurrence of alanine and glycine. It is likely that the efficiency of translation initiation process can depend on the interplay between 5'-context part, 3'-context part and the type of amino acid in the second position of the encoded protein.

Transcripts from a novel BMPR2 termination mutation escape nonsense mediated decay by downstream translation re-initiation: implications for treating pulmonary hypertension

Bone morphogenetic protein receptor type 2 (BMPR2) gene mutations are a major risk factor for heritable pulmonary arterial hypertension (HPAH), an autosomal dominant fatal disease. We have previously shown that BMPR2 transcripts that contain premature termination codon (PTC) mutations are rapidly and nearly completely degraded through nonsense mediated decay (NMD). Here we report a unique PTC mutation (W13X) that did not behave in the predicted manner. We found that patient-derived cultured lymphocytes (CLs) contained readily detectable levels of the PTC-containing transcript. Further analysis suggested that this transcript escaped NMD by translational re-initiation at a downstream Kozak sequence, resulting in the omission of 173 amino acids. Treatment of CLs containing the PTC with an aminoglycoside decreased the truncated protein levels, with a reciprocal increase in full-length BMPR2 protein and, importantly, BMPR-II signaling. This is the first demonstration of aminoglycoside-mediated 'repair' of a BMPR2 mutation at the protein level in patient-derived cells and has obvious implications for treatment of HPAH where no disease-specific treatment options are available. Our data also suggest the need for a more thorough characterization of mutations prior to labeling them as haploinsufficient or dominant negative based simply on sequencing data.

Upstream ORF affects MYCN translation depending on exon 1b alternative splicing

Background

The MYCN gene is transcribed into two major mRNAs: one full-length (MYCN) and one exon 1b-spliced (MYCN^Δ1b) mRNA. But nothing is known about their respective ability to translate the MYCN protein.

Methods

Plasmids were prepared to enable translation from the upstream (uORF) and major ORF of the two MYCN transcripts. Translation was studied after transfection in neuroblastoma SH-EP cell line. Impact of the upstream AUG on translation was evaluated after directed mutagenesis. Functional study with the two MYCN mRNAs was conducted by a cell viability assay. Existence of a new protein encoded by the MYCN^Δ1b uORF was explored by designing a rabbit polyclonal antibody against a specific epitope of this protein.

Results

Both are translated, but higher levels of protein were seen with MYCN^Δ1b mRNA. An upstream ORF was shown to have positive cis-regulatory activity on translation from MYCN but not from MYCN^Δ1b mRNA. In transfected SH-EP neuroblastoma cells, high MYCN dosage obtained with MYCN^Δ1b mRNA translation induces an antiapoptotic effect after serum deprivation that was not observed with low MYCN expression obtained with MYCN mRNA. Here, we showed that MYCNOT: MYCN Overlap Transcript, a new protein of unknown function is translated from the upstream AUG of MYCN^Δ1b mRNA.

Conclusions

Existence of upstream ORF in MYCN transcripts leads to a new level of MYCN regulation. The resulting MYCN dosage has a weak but significant anti-apoptotic activity after intrinsic apoptosis induction.

Differential use of an in-frame translation initiation codon regulates human mu opioid receptor (OPRM1)

The pharmacological effects of morphine and morphine-like drugs are mediated primarily through the micro opioid receptor. Here we show that differential use of an in-frame translational start codon in the 5'-untranslated region of the OPRM1 generates different translational products in vivo and in vitro. The 5'-end of the OPRM1 gene is necessary for initiating the alternate form and for subsequent degradation of the protein. Initiation of OPRM1 at the upstream site decreases the initiation at the main AUG site. However, alternative initiation of the long form of OPRM1 produces a protein with a short half-life, resulting from degradation mediated by the ubiquitin-proteasome pathway. Reporter and degradation assays showed that mutations of this long form at the second and third lysines reduce ubiquitin-dependent proteasome degradation, stabilizing the protein. The data suggest that MOP expression is controlled in part by initiation of the long form of MOP at the alternate site.

Alternative translation start sites and hidden coding potential of eukaryotic mRNAs

It is widely suggested that a eukaryotic mRNA typically contains one translation start site and encodes a single functional protein product. However, according to current points of view on translation initiation mechanisms, eukaryotic ribosomes can recognize several alternative translation start sites and the number of experimentally verified examples of alternative translation is growing rapidly. Also, the frequent occurrence of alternative translation events and their functional significance are supported by the results of computational evaluations. The functional role of alternative translation and its contribution to eukaryotic proteome complexity are discussed.

Conserved upstream open reading frames in higher plants

BACKGROUND

Upstream open reading frames (uORFs) can down-regulate the translation of the main open reading frame (mORF) through two broad mechanisms: ribosomal stalling and reducing reinitiation efficiency. In distantly related plants, such as rice and Arabidopsis, it has been found that conserved uORFs are rare in these transcriptomes with approximately 100 loci. It is unclear how prevalent conserved uORFs are in closely related plants.

RESULTS

We used a homology-based approach to identify conserved uORFs in five cereals (monocots) that could potentially regulate translation. Our approach used a modified reciprocal best hit method to identify putative orthologous sequences that were then analysed by a comparative R-nomics program called uORFSCAN to find conserved uORFs.

CONCLUSION

This research identified new genes that may be controlled at the level of translation by conserved uORFs. We report that conserved uORFs are rare (<150 loci contain them) in cereal transcriptomes, are generally short (less than 100 nt), highly conserved (50% median amino acid sequence similarity), position independent in their 5'-UTRs, and their start codon context and the usage of rare codons for translation does not appear to be important.