ORFik: a comprehensive R toolkit for the analysis of translation

KDM5-mediated transcriptional activation of ribosomal protein genes alters translation efficiency to regulate mitochondrial metabolism in neurons

Genes encoding the KDM5 family of transcriptional regulators are disrupted in individuals with intellectual disability (ID). To understand the link between KDM5 and ID, we characterized five Drosophila strains harboring missense alleles analogous to those observed in patients. These alleles disrupted neuroanatomical development, cognition and other behaviors, and displayed a transcriptional signature characterized by the downregulation of many ribosomal protein genes. A similar transcriptional profile was observed in KDM5C knockout iPSC-induced human glutamatergic neurons, suggesting an evolutionarily conserved role for KDM5 proteins in regulating this class of gene. In Drosophila, reducing KDM5 changed neuronal ribosome composition, lowered the translation efficiency of mRNAs required for mitochondrial function, and altered mitochondrial metabolism. These data highlight the cellular consequences of altered KDM5-regulated transcriptional programs that could contribute to cognitive and behavioral phenotypes. Moreover, they suggest that KDM5 may be part of a broader network of proteins that influence cognition by regulating protein synthesis.

De novo transcriptome assembly of mouse male germ cells reveals novel genes, stage-specific bidirectional promoter activity, and noncoding RNA expression

In mammals, the adult testis is the tissue with the highest diversity in gene expression. Much of that diversity is attributed to germ cells, primarily meiotic spermatocytes and postmeiotic haploid spermatids. Exploiting a newly developed cell purification method, we profiled the transcriptomes of such postmitotic germ cells of mice. We used a de novo transcriptome assembly approach and identified thousands of novel expressed transcripts characterized by features distinct from those of known genes. Novel loci tend to be short in length, monoexonic, and lowly expressed. Most novel genes have arisen recently in evolutionary time and possess low coding potential. Nonetheless, we identify several novel protein-coding genes harboring open reading frames that encode proteins containing matches to conserved protein domains. Analysis of mass-spectrometry data from adult mouse testes confirms protein production from several of these novel genes. We also examine overlap between transcripts and repetitive elements. We find that although distinct families of repeats are expressed with differing temporal dynamics during spermatogenesis, we do not observe a general mode of regulation wherein repeats drive expression of nonrepetitive sequences in a cell type-specific manner. Finally, we observe many fairly long antisense transcripts originating from canonical gene promoters, pointing to pervasive bidirectional promoter activity during spermatogenesis that is distinct and more frequent compared with somatic cells.

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.

Protein-coding circular RNAs - mechanism, detection, and their role in cancer and neurodegenerative diseases

Circular RNAs (circRNAs) are a unique class of non-coding RNAs and were originally thought to have no protein-coding potential due to their lack of a 5' cap and 3' poly(A) tail. However, recent studies have challenged this notion and revealed that some circRNAs have protein-coding potential. They have emerged as a key area of interest in cancer and neurodegeneration research as recent studies have identified several circRNAs that can produce functional proteins with important roles in cancer progression. The protein-coding potential of circRNAs is determined by the presence of an open reading frame (ORF) within the circular structure that can encode a protein. In some cases, the ORF can be translated into a functional protein despite the lack of traditional mRNA features. While the protein-coding potential of most circRNAs remains unclear, several studies have identified specific circRNAs that can produce functional proteins. Understanding the protein-coding potential of circRNAs is important for unravelling their biological functions and potential roles in disease. Our review provides comprehensive coverage of recent advances in the field of circRNA protein-coding capacity and its impact on cancer and neurodegenerative diseases pathogenesis and progression.

Generation of host-directed and virus-specific antivirals using targeted protein degradation promoted by small molecules and viral RNA mimics

Targeted protein degradation (TPD), as exemplified by proteolysis-targeting chimera (PROTAC), is an emerging drug discovery platform. PROTAC molecules, which typically contain a target protein ligand linked to an E3 ligase ligand, recruit a target protein to the E3 ligase to induce its ubiquitination and degradation. Here, we applied PROTAC approaches to develop broad-spectrum antivirals targeting key host factors for many viruses and virus-specific antivirals targeting unique viral proteins. For host-directed antivirals, we identified a small-molecule degrader, FM-74-103, that elicits selective degradation of human GSPT1, a translation termination factor. FM-74-103-mediated GSPT1 degradation inhibits both RNA and DNA viruses. Among virus-specific antivirals, we developed viral RNA oligonucleotide-based bifunctional molecules (Destroyers). As a proof of principle, RNA mimics of viral promoter sequences were used as heterobifunctional molecules to recruit and target influenza viral polymerase for degradation. This work highlights the broad utility of TPD to rationally design and develop next-generation antivirals.

Direct and indirect consequences of PAB1 deletion in the regulation of translation initiation, translation termination, and mRNA decay

Cytoplasmic poly(A)-binding protein (PABPC; Pab1 in yeast) is thought to be involved in multiple steps of post-transcriptional control, including translation initiation, translation termination, and mRNA decay. To understand these roles of PABPC in more detail for endogenous mRNAs, and to distinguish its direct effects from indirect effects, we have employed RNA-Seq and Ribo-Seq to analyze changes in the abundance and translation of the yeast transcriptome, as well as mass spectrometry to assess the abundance of the components of the yeast proteome, in cells lacking the PAB1 gene. We observed drastic changes in the transcriptome and proteome, as well as defects in translation initiation and termination, in pab1Δ cells. Defects in translation initiation and the stabilization of specific classes of mRNAs in pab1Δ cells appear to be partly indirect consequences of reduced levels of specific initiation factors, decapping activators, and components of the deadenylation complex in addition to the general loss of Pab1's direct role in these processes. Cells devoid of Pab1 also manifested a nonsense codon readthrough phenotype indicative of a defect in translation termination, but this defect may be a direct effect of the loss of Pab1 as it could not be attributed to significant reductions in the levels of release factors.

Annotation of uORFs in the OMIM genes allows to reveal pathogenic variants in 5'UTRs

An increasing number of studies emphasize the role of non-coding variants in the development of hereditary diseases. However, the interpretation of such variants in clinical genetic testing still remains a critical challenge due to poor knowledge of their pathogenicity mechanisms. It was previously shown that variants in 5'-untranslated regions (5'UTRs) can lead to hereditary diseases due to disruption of upstream open reading frames (uORFs). Here, we performed a manual annotation of upstream translation initiation sites (TISs) in human disease-associated genes from the OMIM database and revealed ∼4.7 thousand of TISs related to uORFs. We compared our TISs with the previous studies and provided a list of 'high confidence' uORFs. Using a luciferase assay, we experimentally validated the translation of uORFs in the ETFDH, PAX9, MAST1, HTT, TTN,GLI2 and COL2A1 genes, as well as existence of N-terminal CDS extension in the ZIC2 gene. Besides, we created a tool to annotate the effects of genetic variants located in uORFs. We revealed the variants from the HGMD and ClinVar databases that disrupt uORFs and thereby could lead to Mendelian disorders. We also showed that the distribution of uORFs-affecting variants differs between pathogenic and population variants. Finally, drawing on manually curated data, we developed a machine-learning algorithm that allows us to predict the TISs in other human genes.

Gene Expression Analysis in gla-Mutant Zebrafish Reveals Enhanced Ca2+ Signaling Similar to Fabry Disease

Fabry disease (FD) is an X-linked inborn metabolic disorder due to partial or complete lysosomal α-galactosidase A deficiency. FD is characterized by progressive renal insufficiency and cardio- and cerebrovascular involvement. Restricted access on Gb3-independent tissue injury experimental models has limited the understanding of FD pathophysiology and delayed the development of new therapies. Accumulating glycosphingolipids, mainly Gb3 and lysoGb3, are Fabry specific markers used in clinical follow up. However, recent studies suggest there is a need for additional markers to monitor FD clinical course or response to treatment. We used a gla-knockout zebrafish (ZF) to investigate alternative biomarkers in Gb3-free-conditions. RNA sequencing was used to identify transcriptomic signatures in kidney tissues discriminating gla-mutant (M) from wild type (WT) ZF. Gene Ontology (GO) and KEGG pathways analysis showed upregulation of immune system activation and downregulation of oxidative phosphorylation pathways in kidneys from M ZF. In addition, upregulation of the Ca²⁺ signaling pathway was also detectable in M ZF kidneys. Importantly, disruption of mitochondrial and lysosome-related pathways observed in M ZF was validated by immunohistochemistry. Thus, this ZF model expands the pathophysiological understanding of FD, the Gb3-independent effects of gla mutations could be used to explore new therapeutic targets for FD.

Isoform-resolved mRNA profiling of ribosome load defines interplay of HIF and mTOR dysregulation in kidney cancer

Hypoxia inducible factor (HIF) and mammalian target of rapamycin (mTOR) pathways orchestrate responses to oxygen and nutrient availability. These pathways are frequently dysregulated in cancer, but their interplay is poorly understood, in part because of difficulties in simultaneous measurement of global and mRNA-specific translation. Here, we describe a workflow for measurement of ribosome load of mRNAs resolved by their transcription start sites (TSSs). Its application to kidney cancer cells reveals extensive translational reprogramming by mTOR, strongly affecting many metabolic enzymes and pathways. By contrast, global effects of HIF on translation are limited, and we do not observe reported translational activation by HIF2A. In contrast, HIF-dependent alterations in TSS usage are associated with robust changes in translational efficiency in a subset of genes. Analyses of the interplay of HIF and mTOR reveal that specific classes of HIF1A and HIF2A transcriptional target gene manifest different sensitivity to mTOR, in a manner that supports combined use of HIF2A and mTOR inhibitors in treatment of kidney cancer.

A high-resolution map of human RNA translation

Translated small open reading frames (smORFs) can have important regulatory roles and encode microproteins, yet their genome-wide identification has been challenging. We determined the ribosome locations across six primary human cell types and five tissues and detected 7,767 smORFs with translational profiles matching those of known proteins. The human genome was found to contain highly cell-type- and tissue-specific smORFs and a subset that encodes highly conserved amino acid sequences. Changes in the translational efficiency of upstream-encoded smORFs (uORFs) and the corresponding main ORFs predominantly occur in the same direction. Integration with 456 mass-spectrometry datasets confirms the presence of 603 small peptides at the protein level in humans and provides insights into the subcellular localization of these small proteins. This study provides a comprehensive atlas of high-confidence translated smORFs derived from primary human cells and tissues in order to provide a more complete understanding of the translated human genome.

Single-cell RNA-seq-based proteogenomics identifies glioblastoma-specific transposable elements encoding HLA-I-presented peptides

We analyze transposable elements (TEs) in glioblastoma (GBM) patients using a proteogenomic pipeline that combines single-cell transcriptomics, bulk RNA sequencing (RNA-seq) samples from tumors and healthy-tissue cohorts, and immunopeptidomic samples. We thus identify 370 human leukocyte antigen (HLA)-I-bound peptides encoded by TEs differentially expressed in GBM. Some of the peptides are encoded by repeat sequences from intact open reading frames (ORFs) present in up to several hundred TEs from recent long interspersed nuclear element (LINE)-1, long terminal repeat (LTR), and SVA subfamilies. Other HLA-I-bound peptides are encoded by single copies of TEs from old subfamilies that are expressed recurrently in GBM tumors and not expressed, or very infrequently and at low levels, in healthy tissues (including brain). These peptide-coding, GBM-specific, highly recurrent TEs represent potential tumor-specific targets for cancer immunotherapies.

Variation in upstream open reading frames contributes to allelic diversity in maize protein abundance

SignificanceProteins are the machinery which execute essential cellular functions. However, measuring their abundance within an organism can be difficult and resource-intensive. Cells use a variety of mechanisms to control protein synthesis from mRNA, including short open reading frames (uORFs) that lie upstream of the main coding sequence. Ribosomes can preferentially translate uORFs instead of the main coding sequence, leading to reduced translation of the main protein. In this study, we show that uORF sequence variation between individuals can lead to different rates of protein translation and thus variable protein abundances. We also demonstrate that natural variation in uORFs occurs frequently and can be linked to whole-plant phenotypes, indicating that uORF sequence variation likely contributes to plant adaptation.

Small Open Reading Frames, How to Find Them and Determine Their Function

Advances in genomics and molecular biology have revealed an abundance of small open reading frames (sORFs) across all types of transcripts. While these sORFs are often assumed to be non-functional, many have been implicated in physiological functions and a significant number of sORFs have been described in human diseases. Thus, sORFs may represent a hidden repository of functional elements that could serve as therapeutic targets. Unlike protein-coding genes, it is not necessarily the encoded peptide of an sORF that enacts its function, sometimes simply the act of translating an sORF might have a regulatory role. Indeed, the most studied sORFs are located in the 5′UTRs of coding transcripts and can have a regulatory impact on the translation of the downstream protein-coding sequence. However, sORFs have also been abundantly identified in non-coding RNAs including lncRNAs, circular RNAs and ribosomal RNAs suggesting that sORFs may be diverse in function. Of the many different experimental methods used to discover sORFs, the most commonly used are ribosome profiling and mass spectrometry. These can confirm interactions between transcripts and ribosomes and the production of a peptide, respectively. Extensions to ribosome profiling, which also capture scanning ribosomes, have further made it possible to see how sORFs impact the translation initiation of mRNAs. While high-throughput techniques have made the identification of sORFs less difficult, defining their function, if any, is typically more challenging. Together, the abundance and potential function of many of these sORFs argues for the necessity of including sORFs in gene annotations and systematically characterizing these to understand their potential functional roles. In this review, we will focus on the high-throughput methods used in the detection and characterization of sORFs and discuss techniques for validation and functional characterization.

Combined nanopore and single-molecule real-time sequencing survey of human betaherpesvirus 5 transcriptome

Long-read sequencing (LRS), a powerful novel approach, is able to read full-length transcripts and confers a major advantage over the earlier gold standard short-read sequencing in the efficiency of identifying for example polycistronic transcripts and transcript isoforms, including transcript length- and splice variants. In this work, we profile the human cytomegalovirus transcriptome using two third-generation LRS platforms: the Sequel from Pacific BioSciences, and MinION from Oxford Nanopore Technologies. We carried out both cDNA and direct RNA sequencing, and applied the LoRTIA software, developed in our laboratory, for the transcript annotations. This study identified a large number of novel transcript variants, including splice isoforms and transcript start and end site isoforms, as well as putative mRNAs with truncated in-frame ORFs (located within the larger ORFs of the canonical mRNAs), which potentially encode N-terminally truncated polypeptides. Our work also disclosed a highly complex meshwork of transcriptional read-throughs and overlaps.