Stop-codon read-through arises largely from molecular errors and is generally nonadaptive

Chimerism: A whole new perspective in gene regulation

The diversity of molecular entities emerging from a single gene are recognized. Several studies have thus established the cellular role(s) of transcript variants and protein isoforms. A step ahead in challenging the central dogma towards expanding molecular diversity is the identification of fusion genes, chimeric transcripts and chimeric proteins that harbor sequences from more than one gene. The mechanisms for generation of chimeras largely follow similar patterns across all levels of gene regulation but also have interdependence and mutual exclusivity. Whole genome and RNA-seq technologies supported by development of computational algorithms and programs for processing datasets have increasingly enabled the identification of fusion genes and chimeric transcripts, while the discovery of chimeric proteins is as yet more subtle. Earlier thought to be associated with cellular transformation, the contribution of chimeric molecules to normal physiology is also realized and found to influence the expression of their parental genes and regulate cellular pathways. This review offers a collective and comprehensive overview of cellular chimeric entities encompassing the mechanisms involved in their generation, insights on their evolution, functions in gene regulation and their current and novel clinical applications.

Most m5C Modifications in Mammalian mRNAs are Nonadaptive

5-Methylation (m5C) on mRNA molecules is a prevalent internal posttranscriptional modification in eukaryotes. Although m5C modification has been reported to regulate some biological processes, whether most mRNA m5C modifications are functional is unknown. To address this question, we analyzed the genome-wide evolutionary characteristics of m5C modifications in protein-coding genes of humans and mice. Our analysis of RNA sequencing data from 13 tissues of both species revealed that (i) the occurrence of m5C modification is exceedingly low, (ii) the fraction of m5Cs decreases with the amount of Cs across genes or tissues, (iii) m5C modifications are mostly unshared between species, and (iv) m5C sites and motifs do not exhibit greater evolutionary conservation. Additionally, we estimate that a large fraction of the observed mRNA m5C modifications may be deleterious. Together, these observations suggest that most m5C modifications in mammalian mRNAs are nonadaptive, which has important implications for understanding the biological significance of m5C and other posttranscriptional modifications.

Deciphering 3'UTR Mediated Gene Regulation Using Interpretable Deep Representation Learning

The 3' untranslated regions (3'UTRs) of messenger RNAs contain many important cis-regulatory elements that are under functional and evolutionary constraints. It is hypothesized that these constraints are similar to grammars and syntaxes in human languages and can be modeled by advanced natural language techniques such as Transformers, which has been very effective in modeling complex protein sequence and structures. Here 3UTRBERT is described, which implements an attention-based language model, i.e., Bidirectional Encoder Representations from Transformers (BERT). 3UTRBERT is pre-trained on aggregated 3'UTR sequences of human mRNAs in a task-agnostic manner; the pre-trained model is then fine-tuned for specific downstream tasks such as identifying RBP binding sites, m6A RNA modification sites, and predicting RNA sub-cellular localizations. Benchmark results show that 3UTRBERT generally outperformed other contemporary methods in each of these tasks. More importantly, the self-attention mechanism within 3UTRBERT allows direct visualization of the semantic relationship between sequence elements and effectively identifies regions with important regulatory potential. It is expected that 3UTRBERT model can serve as the foundational tool to analyze various sequence labeling tasks within the 3'UTR fields, thus enhancing the decipherability of post-transcriptional regulatory mechanisms.

Identification of nonsense variants in the genomes of 15 Murciano-Granadina bucks and analysis of their segregation in parent-offspring trios

Nonsense variants can inactivate gene function by causing the synthesis of truncated proteins or by inducing nonsense mediated decay of messenger RNAs. The occurrence of such variants in the genomes of livestock species is modulated by multiple demographic and selective factors. Even though nonsense variants can have causal effects on embryo lethality, abortions, and disease, their genomic distribution and segregation in domestic goats have not been characterized in depth yet. In this work, we have sequenced the genomes of 15 Murciano-Granadina bucks with an average coverage of 32.92 × ± 1.45 × . Bioinformatic analysis revealed 947 nonsense variants consistently detected with SnpEff and Ensembl-VEP. These variants were especially abundant in the 3'end of the protein-coding regions. Genes related to olfactory perception, ATPase activity coupled to transmembrane movement of substances, defense to virus, hormonal response, and sensory perception of taste were particularly enriched in nonsense variants. Seventeen nonsense variants expected to have harmful effects on fitness were genotyped in parent-offspring trios. We observed that several nonsense variants predicted to be lethal based on mouse knockout data did not have such effect, a finding that could be explained by the existence of multiple mechanisms counteracting lethality. These findings demonstrate that predicting the effects of putative nonsense variants on fitness is extremely challenging. As a matter of fact, such a goal could only be achieved by generating a high quality telomere-to-telomere goat reference genome combined with carefully curated annotation and functional testing of promising candidate variants.

Generating CRISPR-edited clonal lines of cultured Drosophila S2 cells

CRISPR/Cas9 genome editing is a pervasive research tool due to its relative ease of use. However, some systems are not amenable to generating edited clones due to genomic complexity and/or difficulty in establishing clonal lines. For example, Drosophila Schneider 2 (S2) cells possess a segmental aneuploid genome and are challenging to single-cell select. Here, we describe a streamlined CRISPR/Cas9 methodology for knock-in and knock-out experiments in S2 cells, whereby an antibiotic resistance gene is inserted in-frame with the coding region of a gene-of-interest. By using selectable markers, we have improved the ease and efficiency for the positive selection of null cells using antibiotic selection in feeder layers followed by cell expansion to generate clonal lines. Using this method, we generated the first acentrosomal S2 cell lines by knocking-out centriole genes Polo-like Kinase 4/Plk4 or Ana2 as proof of concept. These strategies for generating gene-edited clonal lines will add to the collection of CRISPR tools available for cultured Drosophila cells by making CRISPR more practical and therefore improving gene function studies.

Nkd1 functions downstream of Axin2 to attenuate Wnt signaling

Wnt signaling is a crucial developmental pathway involved in early development as well as stem-cell maintenance in adults and its misregulation leads to numerous diseases. Thus, understanding the regulation of this pathway becomes vitally important. Axin2 and Nkd1 are widely utilized negative feedback regulators in Wnt signaling where Axin2 functions to destabilize cytoplasmic β-catenin, and Nkd1 functions to inhibit the nuclear localization of β-catenin. Here, we set out to further understand how Axin2 and Nkd1 regulate Wnt signaling by creating axin2gh1/gh1, nkd1gh2/gh2 single mutants and axin2gh1/gh1;nkd1gh2/gh2 double mutant zebrafish using sgRNA/Cas9. All three Wnt regulator mutants were viable and had impaired heart looping, neuromast migration defects, and behavior abnormalities in common, but there were no signs of synergy in the axin2gh1/gh1;nkd1gh2/gh2 double mutants. Further, Wnt target gene expression by qRT-PCR and RNA-seq, and protein expression by mass spectrometry demonstrated that the double axin2gh1/gh1;nkd1gh2/gh2 mutant resembled the nkd1gh2/gh2 phenotype demonstrating that Nkd1 functions downstream of Axin2. In support of this, the data further demonstrates that Axin2 uniquely alters the properties of β-catenin-dependent transcription having novel readouts of Wnt activity compared with nkd1gh2/gh2 or the axin2gh1/gh1;nkd1gh2/gh2 double mutant. We also investigated the sensitivity of the Wnt regulator mutants to exacerbated Wnt signaling, where the single mutants displayed characteristic heightened Wnt sensitivity, resulting in an eyeless phenotype. Surprisingly, this phenotype was rescued in the double mutant, where we speculate that cross-talk between Wnt/β-catenin and Wnt/Planar Cell Polarity pathways could lead to altered Wnt signaling in some scenarios. Collectively, the data emphasizes both the commonality and the complexity in the feedback regulation of Wnt signaling.

Mechanisms and Delivery of tRNA Therapeutics

Transfer ribonucleic acid (tRNA) therapeutics will provide personalized and mutation specific medicines to treat human genetic diseases for which no cures currently exist. The tRNAs are a family of adaptor molecules that interpret the nucleic acid sequences in our genes into the amino acid sequences of proteins that dictate cell function. Humans encode more than 600 tRNA genes. Interestingly, even healthy individuals contain some mutant tRNAs that make mistakes. Missense suppressor tRNAs insert the wrong amino acid in proteins, and nonsense suppressor tRNAs read through premature stop signals to generate full length proteins. Mutations that underlie many human diseases, including neurodegenerative diseases, cancers, and diverse rare genetic disorders, result from missense or nonsense mutations. Thus, specific tRNA variants can be strategically deployed as therapeutic agents to correct genetic defects. We review the mechanisms of tRNA therapeutic activity, the nature of the therapeutic window for nonsense and missense suppression as well as wild-type tRNA supplementation. We discuss the challenges and promises of delivering tRNAs as synthetic RNAs or as gene therapies. Together, tRNA medicines will provide novel treatments for common and rare genetic diseases in humans.

Environment modulates protein heterogeneity through transcriptional and translational stop codon readthrough

Stop codon readthrough events give rise to longer proteins, which may alter the protein's function, thereby generating short-lasting phenotypic variability from a single gene. In order to systematically assess the frequency and origin of stop codon readthrough events, we designed a library of reporters. We introduced premature stop codons into mScarlet, which enabled high-throughput quantification of protein synthesis termination errors in E. coli using fluorescent microscopy. We found that under stress conditions, stop codon readthrough may occur at rates as high as 80%, depending on the nucleotide context, suggesting that evolution frequently samples stop codon readthrough events. The analysis of selected reporters by mass spectrometry and RNA-seq showed that not only translation but also transcription errors contribute to stop codon readthrough. The RNA polymerase was more likely to misincorporate a nucleotide at premature stop codons. Proteome-wide detection of stop codon readthrough by mass spectrometry revealed that temperature regulated the expression of cryptic sequences generated by stop codon readthrough in E. coli. Overall, our findings suggest that the environment affects the accuracy of protein production, which increases protein heterogeneity when the organisms need to adapt to new conditions.

Information Theory, Living Systems, and Communication Engineering

Mainstream research on information theory within the field of living systems involves the application of analytical tools to understand a broad range of life processes. This paper is dedicated to an opposite problem: it explores the information theory and communication engineering methods that have counterparts in the data transmission process by way of DNA structures and neural fibers. Considering the requirements of modern multimedia, transmission methods chosen by nature may be different, suboptimal, or even far from optimal. However, nature is known for rational resource usage, so its methods have a significant advantage: they are proven to be sustainable. Perhaps understanding the engineering aspects of methods of nature can inspire a design of alternative green, stable, and low-cost transmission.

Extended stop codon context predicts nonsense codon readthrough efficiency in human cells

Protein synthesis terminates when a stop codon enters the ribosome's A-site. Although termination is efficient, stop codon readthrough can occur when a near-cognate tRNA outcompetes release factors during decoding. Seeking to understand readthrough regulation we used a machine learning approach to analyze readthrough efficiency data from published HEK293T ribosome profiling experiments and compared it to comparable yeast experiments. We obtained evidence for the conservation of identities of the stop codon, its context, and 3'-UTR length (when termination is compromised), but not the P-site codon, suggesting a P-site tRNA role in readthrough regulation. Models trained on data from cells treated with the readthrough-promoting drug, G418, accurately predicted readthrough of premature termination codons arising from CFTR nonsense alleles that cause cystic fibrosis. This predictive ability has the potential to aid development of nonsense suppression therapies by predicting a patient's likelihood of improvement in response to drugs given their nonsense mutation sequence context.

Fitness Effects of Phenotypic Mutations at Proteome-Scale Reveal Optimality of Translation Machinery

Errors in protein translation can lead to non-genetic, phenotypic mutations, including amino acid misincorporations. While phenotypic mutations can increase protein diversity, the systematic characterization of their proteome-wide frequencies and their evolutionary impact has been lacking. Here, we developed a mechanistic model of translation errors to investigate how selection acts on protein populations produced by amino acid misincorporations. We fitted the model to empirical observations of misincorporations obtained from over a hundred mass spectrometry datasets of E. coli and S. cerevisiae. We found that on average 20% to 23% of proteins synthesized in the cell are expected to harbor at least one amino acid misincorporation, and that deleterious misincorporations are less likely to occur. Combining misincorporation probabilities and the estimated fitness effects of amino acid substitutions in a population genetics framework, we found 74% of mistranslation events in E. coli and 94% in S. cerevisiae to be neutral. We further show that the set of available synonymous tRNAs is subject to evolutionary pressure, as the presence of missing tRNAs would increase codon-anticodon cross-reactivity and misincorporation error rates. Overall, we find that the translation machinery is likely optimal in E. coli and S. cerevisiae and that both local solutions at the level of codons and a global solution such as the tRNA pool can mitigate the impact of translation errors. We provide a framework to study the evolutionary impact of codon-specific translation errors and a method for their proteome-wide detection across organisms and conditions.

Novel MFSD7-ATP5I fusion promotes migration and invasion of human sarcoma

Fusion genes have been implicated in the development and progression of several types of sarcomas, serving as valuable diagnostic and prognostic markers, as well as potential therapeutic targets. We discovered a novel major facilitator superfamily domain-containing 7 (MFSD7) and adenosine triphosphate 5I (ATP5I) gene fusion from sarcomas. In this study, the MFSD7-ATP5I fusion transcript was screened using RNA sequencing in 55 sarcoma samples and sixteen normal samples. The MFSD7-ATP5I fusion transcript was detected in 58% of sarcoma samples. The correlation between the expression of MFSD7-ATP5I fusion transcript and clinicopathological information was analyzed, and MFSD7-ATP5I expression is associated with marked pleomorphism and lower tumor necrosis. Cell migration and invasion was significantly reduced by knock-down of MFSD7-ATP5I. Cell migration and invasion was increased by overexpression of MFSD7-ATP5I. A phosphokinase assay demonstrated that MFSD7-ATP5I is involved in the GSK-3 pathway. The current study found that MFSD7-ATP5I is associated with increasing pleomorphism and decreasing necrosis of tumors. And our gain and loss of function experiments prove that MFSD7-ATP5I promotes the invasiveness of tumor cells.

A roadmap to cure CHD2-related disorders

Coalition to Cure CHD2 (CCC) is a patient advocacy group dedicated to improving the lives of those affected by CHD2-related disorders (CHD2-RD) by increasing education, building community, and accelerating research to uncover a cure. CHD2 is a chromatin remodeler that was identified in 2013 as being a genetic cause for developmental and epileptic encephalopathies. Pathogenic changes in CHD2 can cause treatment-resistant epilepsy, intellectual and developmental delays, and autism, and some individuals experience neurodevelopmental regression. There are currently no targeted therapies available for CHD2-related disorders. Haploinsufficiency of CHD2 is a causative mechanism of disease for individuals with pathogenic variants (primarily truncating) in CHD2. Recently, identification of individuals with deletion of nearby gene CHASERR, a regulator of CHD2 gene expression, has established dosage sensitivity in CHD2 and solidified the CHASERR gene as a potential therapeutic target for CHD2 levels. Through collaboration with our community and our scientific advisory board, CCC has created a Roadmap to Cure CHD2 as our guide toward a targeted cure that can benefit our community, with steps including (1) identifying and defining patients, (2) developing models of CHD2, (3) studying models of CHD2, (4) testing therapies, (5) involving patients, and (6) reaching a cure. Despite some of the challenges inherent in CHD2 research including establishing animal and cellular models that recapitulate the CHD2 clinical phenotype, identifying measurable outcomes and reliable biomarkers, or testing emerging therapeutic approaches, CCC continues to engage with our community to support ongoing research that aligns with our priorities. CCC sees new and exciting opportunities for additional research that can move our community toward our common goal of a cure that will improve the lives of individuals and their families now and in the future.

Review of phage display: A jack-of-all-trades and master of most biomolecule display

Phage display was first described by George P. Smith when it was shown that virus particles were capable of presenting foreign proteins on their surface. The technology has paved the way for the evolution of various biomolecules presentation and diverse selection strategies. This unique feature has been applied as a versatile platform for numerous applications in drug discovery, protein engineering, diagnostics, and vaccine development. Over the decades, the limits of biomolecules displayed on phage particles have expanded from peptides to proteomes and even alternative scaffolds. This has allowed phage display to be viewed as a versatile display platform to accommodate various biomolecules ranging from small peptides to larger proteomes which has significantly impacted advancements in the biomedical industry. This review will explore the vast array of biomolecules that have been successfully employed in phage display technology in biomedical research.

Genes for highly abundant proteins in Escherichia coli avoid 5' codons that promote ribosomal initiation

In many species highly expressed genes (HEGs) over-employ the synonymous codons that match the more abundant iso-acceptor tRNAs. Bacterial transgene codon randomization experiments report, however, that enrichment with such "translationally optimal" codons has little to no effect on the resultant protein level. By contrast, consistent with the view that ribosomal initiation is rate limiting, synonymous codon usage following the 5' ATG greatly influences protein levels, at least in part by modifying RNA stability. For the design of bacterial transgenes, for simple codon based in silico inference of protein levels and for understanding selection on synonymous mutations, it would be valuable to computationally determine initiation optimality (IO) scores for codons for any given species. One attractive approach is to characterize the 5' codon enrichment of HEGs compared with the most lowly expressed genes, just as translational optimality scores of codons have been similarly defined employing the full gene body. Here we determine the viability of this approach employing a unique opportunity: for Escherichia coli there is both the most extensive protein abundance data for native genes and a unique large-scale transgene codon randomization experiment enabling objective definition of the 5' codons that cause, rather than just correlate with, high protein abundance (that we equate with initiation optimality, broadly defined). Surprisingly, the 5' ends of native genes that specify highly abundant proteins avoid such initiation optimal codons. We find that this is probably owing to conflicting selection pressures particular to native HEGs, including selection favouring low initiation rates, this potentially enabling high efficiency of ribosomal usage and low noise. While the classical HEG enrichment approach does not work, rendering simple prediction of native protein abundance from 5' codon content futile, we report evidence that initiation optimality scores derived from the transgene experiment may hold relevance for in silico transgene design for a broad spectrum of bacteria.

Evolution of termination codons of proteins and the TAG-TGA paradox

In most eukaryotes and prokaryotes TGA is used at a significantly higher frequency than TAG as termination codon of protein-coding genes. Although this phenomenon has been recognized several years ago, there is no generally accepted explanation for the TAG-TGA paradox. Our analyses of human mutation data revealed that out of the eighteen sense codons that can give rise to a nonsense codon by single base substitution, the CGA codon is exceptional: it gives rise to the TGA stop codon at an order of magnitude higher rate than the other codons. Here we propose that the TAG-TGA paradox is due to methylation and hypermutabilty of CpG dinucleotides. In harmony with this explanation, we show that the coding genomes of organisms with strong CpG methylation have a significant bias for TGA whereas those from organisms that lack CpG methylation use TGA and TAG termination codons with similar probability.

Gene product diversity: adaptive or not?

One gene does not equal one RNA or protein. The genomic revolution has revealed numerous different RNA and protein molecules that can be produced from one gene, such as circular RNAs generated by back-splicing, proteins with residues mismatching the genomic encoding because of RNA editing, and proteins extended in the C terminus via stop codon readthrough in translation. Are these diverse products results of exquisite gene regulations or imprecise biological processes? While there are cases where the gene product diversity appears beneficial, genome-scale patterns suggest that much of this diversity arises from nonadaptive, molecular errors. This finding has important implications for studying the functions of diverse gene products and for understanding the fundamental properties and evolution of cellular life.

Unraveling usnic acid: a comparison of biosynthetic gene clusters between two reindeer lichen (Cladonia rangiferina and C. uncialis)

Lichenized fungi are known for their production of a diversity of secondary metabolites, many of which have broad biological and pharmacological applications. By far the most well-studied of these metabolites is usnic acid. While this metabolite has been well-known and researched for decades, the gene cluster responsible for its production was only recently identified from the species Cladonia uncialis. Usnic acid production varies considerably in the genus Cladonia, even among closely related taxa, and many species, such as C. rangiferina, have been inferred to be incapable of producing the metabolite based on analysis by thin-layer chromatography (TLC). We sequenced and examined the usnic acid biosynthetic gene clusters, or lack thereof, from four closely related Cladonia species (C. oricola, C. rangiferina, C. stygia, and C. subtenuis), and compare them against those of C. uncialis. We complement this comparison with tiered chemical profile analyses to confirm the presence or absence of usnic acid in select samples, using both HPLC and LC-MS. Despite long-standing reporting that C. rangiferina lacks the ability to produce usnic acid, we observed functional gene clusters from the species and detected usnic acid when extracts were examined by LC-MS. By contrast, C. stygia and C. oricola, have been previously described as lacking the ability to produce usnic acid, lacked the gene cluster entirely, and no usnic acid could be detected in C. oricola extracts via HPLC or LC-MS. This work suggests that chemical profiles attained through inexpensive and low-sensitivity methods like TLC may fail to detect low abundance metabolites that can be taxonomically informative. This study also bolsters understanding of the usnic acid gene cluster in lichens, revealing differences among domains of the polyketide synthase which may explain observed differences in expression. These results reinforce the need for comprehensive characterization of lichen secondary metabolite profiles with sensitive LC-MS methods.

Stop codon usage as a window into genome evolution: mutation, selection, biased gene conversion and the TAG paradox

Protein coding genes terminate with one of three stop codons (TAA, TGA, or TAG) that, like synonymous codons, are not employed equally. With TGA and TAG having identical nucleotide content, analysis of their differential usage provides an unusual window into the forces operating on what are ostensibly functionally identical residues. Across genomes and between isochores within the human genome, TGA usage increases with G + C content but, with a common G + C → A + T mutation bias, this cannot be explained by mutation bias-drift equilibrium. Increased usage of TGA in G + C-rich genomes or genomic regions is also unlikely to reflect selection for the optimal stop codon, as TAA appears to be universally optimal, probably because it has the lowest read-through rate. Despite TAA being favored by selection and mutation bias, as with codon usage bias G + C pressure is the prime determinant of between-species TGA usage trends. In species with strong G + C-biased gene conversion (gBGC), such as mammals and birds, the high usage and conservation of TGA is best explained by an A + T → G + C repair bias. How to explain TGA enrichment in other G + C-rich genomes is less clear. Enigmatically, across bacterial and archaeal species and between human isochores TAG usage is mostly unresponsive to G + C pressure. This unresponsiveness we dub the TAG paradox as currently no mutational, selective, or gBGC model provides a well-supported explanation. That TAG does increase with G + C usage across eukaryotes makes the usage elsewhere yet more enigmatic. We suggest resolution of the TAG paradox may provide insights into either an unknown but common selective preference (probably at the DNA/RNA level) or an unrecognized complexity to the action of gBGC.

Novel homozygous nonsense mutation of MLIP and compensatory alternative splicing

Despite the growing accessibility of clinical sequencing, functional interpretation of variants remains a major hurdle to molecular diagnostics of Mendelian diseases. We aimed to describe a new adult-onset myopathy with muscle weakness and hyperCKemia caused by a nonsense variant in muscular LMNA-interacting protein (MLIP). Following RNA-sequencing, differential expression analysis uncovered a significant downregulation of this gene, which had a surprisingly mild effect on MLIP protein expression. RT-PCR and long-read sequencing (LRS) both support an important transcriptome shift in the patient, where decreased MLIP levels are seemingly due to nonsense-mediated decay of transcripts containing the exon 5 mutation. Moreover, a compensatory mechanism upregulates the functionally lacking isoforms and generates novel transcripts. These results support the recently discovered clinical implications of MLIP variants in myopathies, highlighting for the first time its relevance in adult-onset cases. These results also underline the power of LRS as a tool for the functional assessment of variants of unknown significance (VUS), as well as the definition of accurate isoform profile annotations in a tissue-specific manner.

Unusual mammalian usage of TGA stop codons reveals that sequence conservation need not imply purifying selection

The assumption that conservation of sequence implies the action of purifying selection is central to diverse methodologies to infer functional importance. GC-biased gene conversion (gBGC), a meiotic mismatch repair bias strongly favouring GC over AT, can in principle mimic the action of selection, this being thought to be especially important in mammals. As mutation is GC→AT biased, to demonstrate that gBGC does indeed cause false signals requires evidence that an AT-rich residue is selectively optimal compared to its more GC-rich allele, while showing also that the GC-rich alternative is conserved. We propose that mammalian stop codon evolution provides a robust test case. Although in most taxa TAA is the optimal stop codon, TGA is both abundant and conserved in mammalian genomes. We show that this mammalian exceptionalism is well explained by gBGC mimicking purifying selection and that TAA is the selectively optimal codon. Supportive of gBGC, we observe (i) TGA usage trends are consistent at the focal stop codon and elsewhere (in UTR sequences); (ii) that higher TGA usage and higher TAA→TGA substitution rates are predicted by a high recombination rate; and (iii) across species the difference in TAA <-> TGA substitution rates between GC-rich and GC-poor genes is largest in genomes that possess higher between-gene GC variation. TAA optimality is supported both by enrichment in highly expressed genes and trends associated with effective population size. High TGA usage and high TAA→TGA rates in mammals are thus consistent with gBGC’s predicted ability to “drive” deleterious mutations and supports the hypothesis that sequence conservation need not be indicative of purifying selection. A general trend for GC-rich trinucleotides to reside at frequencies far above their mutational equilibrium in high recombining domains supports the generality of these results.

Why Is the UAG (Amber) Stop Codon Almost Absent in Highly Expressed Bacterial Genes?

The genome hypothesis postulates that genes in a genome tend to conform to their species’ usage of the codon catalog and the GC content of the DNA. Thus, codon frequencies differ across organisms, including the three termination codons in the standard genetic code. Here, we analyze the frequencies of stop codons in a group of highly expressed genes from 196 prokaryotes under strong translational selection. The occurrence of the three translation termination codons is highly biased, with UAA (ochre) being the most prevalent in almost all bacteria. In contrast, UAG (amber) is the least frequent termination codon, e.g., only 321 occurrences (7.4%) in E. coli K-12 substr. W3110. Of the 253 highly expressed genes, only two end with an UAG codon. The strength of the selective bias against UAG in highly expressed genes varies among bacterial genomes, but it is not affected by the GC content of these genomes. In contrast, increased GC content results in a decrease in UAA abundance with a concomitant increase in UGA abundance. We propose that readthrough efficiency and context effects could explain the prevalence of UAA over UAG, particularly in highly expressed genes. Findings from this communication can be utilized for the optimization of gene expression.

Variation in Release Factor Abundance Is Not Needed to Explain Trends in Bacterial Stop Codon Usage

In bacteria stop codons are recognized by one of two class I release factors (RF1) recognizing TAG, RF2 recognizing TGA, and TAA being recognized by both. Variation across bacteria in the relative abundance of RF1 and RF2 is thus hypothesized to select for different TGA/TAG usage. This has been supported by correlations between TAG:TGA ratios and RF1:RF2 ratios across multiple bacterial species, potentially also explaining why TAG usage is approximately constant despite extensive variation in GC content. It is, however, possible that stop codon trends are determined by other forces and that RF ratios adapt to stop codon usage, rather than vice versa. Here, we determine which direction of the causal arrow is the more parsimonious. Our results support the notion that RF1/RF2 ratios become adapted to stop codon usage as the same trends, notably the anomalous TAG behavior, are seen in contexts where RF1:RF2 ratios cannot be, or are unlikely to be, causative, that is, at 3'untranslated sites never used for translation termination, in intragenomic analyses, and across archaeal species (that possess only one RF1). We conclude that specifics of RF biology are unlikely to fully explain TGA/TAG relative usage. We discuss why the causal relationships for the evolution of synonymous stop codon usage might be different from those affecting synonymous sense codon usage, noting that transitions between TGA and TAG require two-point mutations one of which is likely to be deleterious.

Mammalian splicing divergence is shaped by drift, buffering in trans, and a scaling law

This study globally investigates the allelic splicing pattern in multiple tissues of an F1 hybrid mouse and reveals the underlying driving forces shaping such tissue-dependent splicing divergence.

Alternative splicing is ubiquitous, but the mechanisms underlying its pattern of evolutionary divergence across mammalian tissues are still underexplored. Here, we investigated the cis-regulatory divergences and their relationship with tissue-dependent trans-regulation in multiple tissues of an F1 hybrid between two mouse species. Large splicing changes between tissues are highly conserved and likely reflect functional tissue-dependent regulation. In particular, micro-exons frequently exhibit this pattern with high inclusion levels in the brain. Cis-divergence of splicing appears to be largely non-adaptive. Although divergence is in general associated with higher densities of sequence variants in regulatory regions, events with high usage of the dominant isoform apparently tolerate more mutations, explaining why their exon sequences are highly conserved but their intronic splicing site flanking regions are not. Moreover, we demonstrate that non-adaptive mutations are often masked in tissues where accurate splicing likely is more important, and experimentally attribute such buffering effect to trans-regulatory splicing efficiency.

Accurate Prediction of Protein Sequences for Proteogenomics Data Integration

This book chapter discusses proteogenomics data integration and provides an overview into the different omics layer involved in defining the proteome of a living organism. Various aspects of genome variability affecting either the sequence or abundance level of proteins are discussed in this book chapter, such as the effect of single-nucleotide variants or larger genomic structural variants on the proteome. Next, various sequencing technologies are introduced and discussed from a proteogenomics data integration perspective such as those providing short- and long-read sequencing and listing their respective advantages and shortcomings for accurate protein variant prediction using genomic/transcriptomics sequencing data. Finally, the various bioinformatics tools used to process and analyze DNA/RNA sequencing data are discussed with the ultimate goal of obtaining accurately predicted sample-specific protein sequences that can be used as a drop-in replacement in existing approaches for peptide and protein identification using popular database search engines such as MSFragger, SearchGUI/PeptideShaker.

Mutations That Affect the Surface Expression of TRPV6 Are Associated with the Upregulation of Serine Proteases in the Placenta of an Infant

Recently, we reported a case of an infant with neonatal severe under-mineralizing skeletal dysplasia caused by mutations within both alleles of the TRPV6 gene. One mutation results in an in frame stop codon (R510stop) that leads to a truncated, nonfunctional TRPV6 channel, and the second in a point mutation (G660R) that, surprisingly, does not affect the Ca2+ permeability of TRPV6. We mimicked the subunit composition of the unaffected heterozygous parent and child by coexpressing the TRPV6 G660R and R510stop mutants and combinations with wild type TRPV6. We show that both the G660R and R510stop mutant subunits are expressed and result in decreased calcium uptake, which is the result of the reduced abundancy of functional TRPV6 channels within the plasma membrane. We compared the proteomic profiles of a healthy placenta with that of the diseased infant and detected, exclusively in the latter two proteases, HTRA1 and cathepsin G. Our results implicate that the combination of the two mutant TRPV6 subunits, which are expressed in the placenta of the diseased child, is responsible for the decreased calcium uptake, which could explain the skeletal dysplasia. In addition, placental calcium deficiency also appears to be associated with an increase in the expression of proteases.

Stop codon readthrough alters the activity of a POU/Oct transcription factor during Drosophila development

Background

A number of cellular processes have evolved in metazoans that increase the proteome repertoire in relation to the genome, such as alternative splicing and translation recoding. Another such process, translational stop codon readthrough (SCR), generates C-terminally extended protein isoforms in many eukaryotes, including yeast, plants, insects, and humans. While comparative genome analyses have predicted the existence of programmed SCR in many species including humans, experimental proof of its functional consequences are scarce.

Results

We show that SCR of the Drosophila POU/Oct transcription factor Ventral veins lacking/Drifter (Vvl/Dfr) mRNA is prevalent in certain tissues in vivo, reaching a rate of 50% in the larval prothoracic gland. Phylogenetically, the C-terminal extension is conserved and harbors intrinsically disordered regions and amino acid stretches implied in transcriptional activation. Elimination of Vvl/Dfr translational readthrough by CRISPR/Cas9 mutagenesis changed the expression of a large number of downstream genes involved in processes such as chromatin regulation, neurogenesis, development, and immune response. As a proof-of-principle, we demonstrate that the C-terminal extension of Vvl/Dfr is necessary for correct timing of pupariation, by increasing the capacity to regulate its target genes. The extended Vvl/Dfr isoform acts in synergy with the transcription factor Molting defective (Mld) to increase the expression and biosynthesis of the steroid hormone ecdysone, thereby advancing pupariation. Consequently, late-stage larval development was prolonged and metamorphosis delayed in vvl/dfr readthrough mutants.

Conclusions

We demonstrate that translational recoding of a POU/Oct transcription factor takes place in a highly tissue-specific and temporally controlled manner. This dynamic and regulated recoding is necessary for normal expression of a large number of genes involved in many cellular and developmental processes. Loss of Vvl/Dfr translational readthrough negatively affects steroid hormone biosynthesis and delays larval development and progression into metamorphosis. Thus, this study demonstrates how SCR of a transcription factor can act as a developmental switch in a spatiotemporal manner, feeding into the timing of developmental transitions between different life-cycle stages.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01106-0.

Mammalian circular RNAs result largely from splicing errors

Ubiquitous in eukaryotes, circular RNAs (circRNAs) comprise a large class of mostly non-coding RNAs produced by back-splicing. Although some circRNAs have demonstrated biochemical activities, whether most circRNAs are functional is unknown. Here, we test the hypothesis that circRNA production primarily results from splicing error and so is deleterious instead of beneficial. In support of the error hypothesis, our analysis of RNA sequencing data from 11 shared tissues of humans, macaques, and mice finds that (1) back-splicing is much rarer than linear-splicing, (2) the rate of back-splicing diminishes with the splicing amount, (3) the overall prevalence of back-splicing in a species declines with its effective population size, and (4) circRNAs are overall evolutionarily unconserved. We estimate that more than 97% of the observed circRNA production is deleterious. We identify a small number of functional circRNA candidates, and the genome-wide trend strongly suggests that circRNAs are largely non-functional products of splicing errors.

Effective Population Size Predicts Local Rates but Not Local Mitigation of Read-through Errors

In correctly predicting that selection efficiency is positively correlated with the effective population size (Ne), the nearly neutral theory provides a coherent understanding of between-species variation in numerous genomic parameters, including heritable error (germline mutation) rates. Does the same theory also explain variation in phenotypic error rates and in abundance of error mitigation mechanisms? Translational read-through provides a model to investigate both issues as it is common, mostly nonadaptive, and has good proxy for rate (TAA being the least leaky stop codon) and potential error mitigation via "fail-safe" 3' additional stop codons (ASCs). Prior theory of translational read-through has suggested that when population sizes are high, weak selection for local mitigation can be effective thus predicting a positive correlation between ASC enrichment and Ne. Contra to prediction, we find that ASC enrichment is not correlated with Ne. ASC enrichment, although highly phylogenetically patchy, is, however, more common both in unicellular species and in genes expressed in unicellular modes in multicellular species. By contrast, Ne does positively correlate with TAA enrichment. These results imply that local phenotypic error rates, not local mitigation rates, are consistent with a drift barrier/nearly neutral model.

Context-aware synthetic biology by controller design: Engineering the mammalian cell

The rise of systems biology has ushered a new paradigm: the view of the cell as a system that processes environmental inputs to drive phenotypic outputs. Synthetic biology provides a complementary approach, allowing us to program cell behavior through the addition of synthetic genetic devices into the cellular processor. These devices, and the complex genetic circuits they compose, are engineered using a design-prototype-test cycle, allowing for predictable device performance to be achieved in a context-dependent manner. Within mammalian cells, context effects impact device performance at multiple scales, including the genetic, cellular, and extracellular levels. In order for synthetic genetic devices to achieve predictable behaviors, approaches to overcome context dependence are necessary. Here, we describe control systems approaches for achieving context-aware devices that are robust to context effects. We then consider cell fate programing as a case study to explore the potential impact of context-aware devices for regenerative medicine applications.

The Role of N6 -Methyladenosine Modified Circular RNA in Pathophysiological Processes

Circular RNA (circRNA) is a type of covalently closed and endogenous non-coding RNA (ncRNA) with tissue- and cell-specific expression patterns generated by a non-canonical splicing event. Previous reports have indicated that circRNAs exert their functions in different ways, thereby participating in various pathophysiological processes. N6 -methyladenosine (m6A) methylation occurs in the N6-position, which is the most abundant and conserved internal transcriptional modification in eukaryotes, including mRNA and ncRNAs. Accumulating evidences confirm that m6A modification also exists in the circRNA and greatly affects the biological functions of circRNA. Their dysregulated expression can be a cause of various pathophysiological processes, such as spermatogenesis, myoblast differentiation, cancer, cardiovascular disease, mental illness and so on. Understanding the role of m6A-modified circRNAs in pathophysiological processes may contribute to better understanding the physiological mechanisms and develop new biomarkers. This review summarizes the regulatory mechanism of m6A modification on circRNA metabolism and the role of m6A-modified circRNAs in pathophysiological processes. This article may pave the way for a better understanding of the role of epigenetically modified circRNAs in pathophysiological process.

Transcriptome-wide investigation of stop codon readthrough in Saccharomyces cerevisiae

Translation of mRNA into a polypeptide is terminated when the release factor eRF1 recognizes a UAA, UAG, or UGA stop codon in the ribosomal A site and stimulates nascent peptide release. However, stop codon readthrough can occur when a near-cognate tRNA outcompetes eRF1 in decoding the stop codon, resulting in the continuation of the elongation phase of protein synthesis. At the end of a conventional mRNA coding region, readthrough allows translation into the mRNA 3'-UTR. Previous studies with reporter systems have shown that the efficiency of termination or readthrough is modulated by cis-acting elements other than stop codon identity, including two nucleotides 5' of the stop codon, six nucleotides 3' of the stop codon in the ribosomal mRNA channel, and stem-loop structures in the mRNA 3'-UTR. It is unknown whether these elements are important at a genome-wide level and whether other mRNA features proximal to the stop codon significantly affect termination and readthrough efficiencies in vivo. Accordingly, we carried out ribosome profiling analyses of yeast cells expressing wild-type or temperature-sensitive eRF1 and developed bioinformatics strategies to calculate readthrough efficiency, and to identify mRNA and peptide features which influence that efficiency. We found that the stop codon (nt +1 to +3), the nucleotide after it (nt +4), the codon in the P site (nt -3 to -1), and 3'-UTR length are the most influential features in the control of readthrough efficiency, while nts +5 to +9 had milder effects. Additionally, we found low readthrough genes to have shorter 3'-UTRs compared to high readthrough genes in cells with thermally inactivated eRF1, while this trend was reversed in wild-type cells. Together, our results demonstrated the general roles of known regulatory elements in genome-wide regulation and identified several new mRNA or peptide features affecting the efficiency of translation termination and readthrough.

Readthrough Errors Purge Deleterious Cryptic Sequences, Facilitating the Birth of Coding Sequences

De novo protein-coding innovations sometimes emerge from ancestrally noncoding DNA, despite the expectation that translating random sequences is overwhelmingly likely to be deleterious. The "preadapting selection" hypothesis claims that emergence is facilitated by prior, low-level translation of noncoding sequences via molecular errors. It predicts that selection on polypeptides translated only in error is strong enough to matter and is strongest when erroneous expression is high. To test this hypothesis, we examined noncoding sequences located downstream of stop codons (i.e., those potentially translated by readthrough errors) in Saccharomyces cerevisiae genes. We identified a class of "fragile" proteins under strong selection to reduce readthrough, which are unlikely substrates for co-option. Among the remainder, sequences showing evidence of readthrough translation, as assessed by ribosome profiling, encoded C-terminal extensions with higher intrinsic structural disorder, supporting the preadapting selection hypothesis. The cryptic sequences beyond the stop codon, rather than spillover effects from the regular C-termini, are primarily responsible for the higher disorder. Results are robust to controlling for the fact that stronger selection also reduces the length of C-terminal extensions. These findings indicate that selection acts on 3' UTRs in Saccharomyces cerevisiae to purge potentially deleterious variants of cryptic polypeptides, acting more strongly in genes that experience more readthrough errors.

Mammalian Alternative Translation Initiation Is Mostly Nonadaptive

Alternative translation initiation (ATLI) refers to the existence of multiple translation initiation sites per gene and is a widespread phenomenon in eukaryotes. ATLI is commonly assumed to be advantageous through creating proteome diversity or regulating protein synthesis. We here propose an alternative hypothesis that ATLI arises primarily from nonadaptive initiation errors presumably due to the limited ability of ribosomes to distinguish sequence motifs truly signaling translation initiation from similar sequences. Our hypothesis, but not the adaptive hypothesis, predicts a series of global patterns of ATLI, all of which are confirmed at the genomic scale by quantitative translation initiation sequencing in multiple human and mouse cell lines and tissues. Similarly, although many codons differing from AUG by one nucleotide can serve as start codons, our analysis suggests that using non-AUG start codons is mostly disadvantageous. These and other findings strongly suggest that ATLI predominantly results from molecular error, requiring a major revision of our understanding of the precision and regulation of translation initiation.

Characterization of Bicistronic Transcription in Budding Yeast

Bicistronic transcripts (operon-like transcripts) have occasionally been reported in eukaryotes, including unicellular yeasts, plants, and humans, despite the fact that they lack trans-splice mechanisms. However, the characteristics of eukaryotic bicistronic transcripts are poorly understood, except for those in nematodes. Here, we describe the genomic, transcriptomic, and ribosome profiling features of bicistronic transcripts in unicellular yeasts. By comparing the expression level of bicistronic transcripts with their monocistronic equivalents, we identify two main categories of bicistronic transcripts: highly and lowly expressed. These two categories exhibit quite different features. First, highly expressed bicistronic transcripts have higher conservation within and between strains and shorter intergenic spacers with higher GC content and less stable secondary structure. Second, genes in highly expressed bicistronic transcripts have lower translation efficiency, with the second gene showing statistically significant lower translation efficiency than the first. Finally, the genes found in these highly expressed bicistronic transcripts tend to be younger, with more recent origins. Together, these results suggest that bicistronic transcripts in yeast are heterogeneous. We further propose that at least some highly expressed bicistronic transcripts appear to play a role in modulating monocistronic translation.IMPORTANCE Operons, where a single mRNA transcript encodes multiple adjacent proteins, are a widespread feature of bacteria and archaea. In contrast, the genes of eukaryotes are generally considered monocistronic. However, a number of studies have revealed the presence of bicistronic transcripts in eukaryotes, including humans. The basic features of these transcripts are largely unknown in eukaryotes, especially in organisms lacking trans-splice mechanisms. Our analyses characterize bicistronic transcripts in one such eukaryotic group, yeasts. We show that highly expressed bicistronic transcripts have unusual features compared to lowly expressed bicistronic transcripts, with several features influencing translational modulation.

In silico identification of novel open reading frames in Plasmodium falciparum oocyte and salivary gland sporozoites using proteogenomics framework

Background

Plasmodium falciparum causes the deadliest form of malaria, which remains one of the most prevalent infectious diseases. Unfortunately, the only licensed vaccine showed limited protection and resistance to anti-malarial drug is increasing, which can be largely attributed to the biological complexity of the parasite’s life cycle. The progression from one developmental stage to another in P. falciparum involves drastic changes in gene expressions, where its infectivity to human hosts varies greatly depending on the stage. Approaches to identify candidate genes that are responsible for the development of infectivity to human hosts typically involve differential gene expression analysis between stages. However, the detection may be limited to annotated proteins and open reading frames (ORFs) predicted using restrictive criteria.

Methods

The above problem is particularly relevant for P. falciparum; whose genome annotation is relatively incomplete given its clinical significance. In this work, systems proteogenomics approach was used to address this challenge, as it allows computational detection of unannotated, novel Open Reading Frames (nORFs), which are neglected by conventional analyses. Two pairs of transcriptome/proteome were obtained from a previous study where one was collected in the mosquito-infectious oocyst sporozoite stage, and the other in the salivary gland sporozoite stage with human infectivity. They were then re-analysed using the proteogenomics framework to identify nORFs in each stage.

Results

Translational products of nORFs that map to antisense, intergenic, intronic, 3′ UTR and 5′ UTR regions, as well as alternative reading frames of canonical proteins were detected. Some of these nORFs also showed differential expression between the two life cycle stages studied. Their regulatory roles were explored through further bioinformatics analyses including the expression regulation on the parent reference genes, in silico structure prediction, and gene ontology term enrichment analysis.

Conclusion

The identification of nORFs in P. falciparum sporozoites highlights the biological complexity of the parasite. Although the analyses are solely computational, these results provide a starting point for further experimental validation of the existence and functional roles of these nORFs,

Sequence- and structure-specific cytosine-5 mRNA methylation by NSUN6

The highly abundant N6-methyladenosine (m6A) RNA modification affects most aspects of mRNA function, yet the precise function of the rarer 5-methylcytidine (m5C) remains largely unknown. Here, we map m5C in the human transcriptome using methylation-dependent individual-nucleotide resolution cross-linking and immunoprecipitation (miCLIP) combined with RNA bisulfite sequencing. We identify NSUN6 as a methyltransferase with strong substrate specificity towards mRNA. NSUN6 primarily targeted three prime untranslated regions (3'UTR) at the consensus sequence motif CTCCA, located in loops of hairpin structures. Knockout and rescue experiments revealed enhanced mRNA and translation levels when NSUN6-targeted mRNAs were methylated. Ribosome profiling further demonstrated that NSUN6-specific methylation correlated with translation termination. While NSUN6 was dispensable for mouse embryonic development, it was down-regulated in human tumours and high expression of NSUN6 indicated better patient outcome of certain cancer types. In summary, our study identifies NSUN6 as a methyltransferase targeting mRNA, potentially as part of a quality control mechanism involved in translation termination fidelity.

Insights Into Translatomics in the Nervous System

Most neurological disorders are caused by abnormal gene translation. Generally, dysregulation of elements involved in the translational process disrupts homeostasis in neurons and neuroglia. Better understanding of how the gene translation process occurs requires detailed analysis of transcriptomic and proteomic profile data. However, a lack of strictly direct correlations between mRNA and protein levels limits translational investigation by combining transcriptomic and proteomic profiling. The much better correlation between proteins and translated mRNAs than total mRNAs in abundance and insufficiently sensitive proteomics approach promote the requirement of advances in translatomics technology. Translatomics which capture and sequence the mRNAs associated with ribosomes has been effective in identifying translational changes by genetics or projections, ribosome stalling, local translation, and transcript isoforms in the nervous system. Here, we place emphasis on the main three translatomics methods currently used to profile mRNAs attached to ribosome-nascent chain complex (RNC-mRNA). Their prominent applications in neurological diseases including glioma, neuropathic pain, depression, fragile X syndrome (FXS), neurodegenerative disorders are outlined. The content reviewed here expands our understanding on the contributions of aberrant translation to neurological disease development.

Revisiting the phosphotyrosine binding pocket of Fyn SH2 domain led to the identification of novel SH2 superbinders

Protein engineering through directed evolution is an effective way to obtain proteins with novel functions with the potential applications as tools for diagnosis or therapeutics. Many natural proteins have undergone directed evolution in vitro in the test tubes in the laboratories worldwide, resulting in the numerous protein variants with novel or enhanced functions. we constructed here an SH2 variant library by randomizing 8 variable residues in its phosphotyrosine (pTyr) binding pocket. Selection of this library by a pTyr peptide led to the identification of SH2 variants with enhanced affinities measured by EC50. Fluorescent polarization was then applied to quantify the binding affinities of the newly identified SH2 variants. As a result, three SH2 variants, named V3, V13 and V24, have comparable binding affinities with the previously identified SH2 triple-mutant superbinder. Biolayer Interferometry assay was employed to disclose the kinetics of the binding of these SH2 superbinders to the phosphotyrosine peptide. The results indicated that all the SH2 superbinders have two-orders increase of the dissociation rate when binding the pTyr peptide while there was no significant change in their associate rates. Intriguingly, though binding the pTyr peptide with comparable affinity with other SH2 superbinders, the V3 does not bind to the sTyr peptide. However, variant V13 and V24 have cross-reactivity with both pTyr and sTyr peptides. The newly identified superbinders could be utilized as tools for the identification of pTyr-containing proteins from tissues under different physiological or pathophysiological conditions and may have the potential in the therapeutics.

Introduction of a leaky stop codon as molecular tool in Chlamydomonas reinhardtii

Expression of proteins in the chloroplast or mitochondria of the model green alga Chlamydomonas reinhardtii can be achieved by directly inserting transgenes into organellar genomes, or through nuclear expression and post-translational import. A number of tools have been developed in the literature for achieving high expression levels from the nuclear genome despite messy genomic integration and widespread silencing of transgenes. Here, recent advances in the field are combined and two systems of bicistronic expression, based on ribosome reinitiation or ribosomal skip induced by a viral 2A sequence, are compared side-by-side. Further, the small subunit of Rubisco (RBCS) was developed as a functional nuclear reporter for successful chloroplast import and restoration of photosynthesis: To be able to combine RBCS with a Venus fluorescent reporter without compromising photosynthetic activity, a leaky stop codon is introduced as a novel molecular tool that allows the simultaneous expression of functional and fluorescently tagged versions of the protein from a single construct.

A few good peptides: MHC class I-based cancer immunosurveillance and immunoevasion

The remarkable success of immune checkpoint inhibitors demonstrates the potential of tumour-specific CD8⁺ T cells to prevent and treat cancer. Although the number of lives saved by immunotherapy mounts, only a relatively small fraction of patients are cured. Here, we review two of the factors that limit the application of CD8⁺ T cell immunotherapies: difficulties in identifying tumour-specific peptides presented by MHC class I molecules and the ability of tumour cells to impair antigen presentation as they evolve under T cell selection. We describe recent advances in understanding how peptides are generated from non-canonical translation of defective ribosomal products, relate this to the dysregulated translation that is a feature of carcinogenesis and propose dysregulated translation as an important new source of tumour-specific peptides. We discuss how the synthesis and function of components of the antigen-processing and presentation pathway, including the recently described immunoribosome, are manipulated by tumours for immunoevasion and point to common druggable targets that may enhance immunotherapy.

Selection Shapes Synonymous Stop Codon Use in Mammals

Phylogenetic models of the evolution of protein-coding sequences can provide insights into the selection pressures that have shaped them. In the application of these models synonymous nucleotide substitutions, which do not alter the encoded amino acid, are often assumed to have limited functional consequences and used as a proxy for the neutral rate of evolution. The ratio of nonsynonymous to synonymous substitution rates is then used to categorize the selective regime that applies to the protein (e.g., purifying selection, neutral evolution, diversifying selection). Here, we extend the Muse and Gaut model of codon evolution to explore the extent of purifying selection acting on substitutions between synonymous stop codons. Using a large collection of coding sequence alignments, we estimate that a high proportion (approximately 57%) of mammalian genes are affected by selection acting on stop codon preference. This proportion varies substantially by codon, with UGA stop codons far more likely to be conserved. Genes with evidence of selection acting on synonymous stop codons have distinctive characteristics, compared to unconserved genes with the same stop codon, including longer [Formula: see text] untranslated regions (UTRs) and shorter mRNA half-life. The coding regions of these genes are also much more likely to be under strong purifying selection pressure. Our results suggest that the preference for UGA stop codons found in many multicellular eukaryotes is selective rather than mutational in origin.

Premature termination codon readthrough in Drosophila varies in a developmental and tissue-specific manner

Despite their essential function in terminating translation, readthrough of stop codons occurs more frequently than previously supposed. However, little is known about the regulation of stop codon readthrough by anatomical site and over the life cycle of animals. Here, we developed a set of reporters to measure readthrough in Drosophila melanogaster. A focused RNAi screen in whole animals identified upf1 as a mediator of readthrough, suggesting that the stop codons in the reporters were recognized as premature termination codons (PTCs). We found readthrough rates of PTCs varied significantly throughout the life cycle of flies, being highest in older adult flies. Furthermore, readthrough rates varied dramatically by tissue and, intriguingly, were highest in fly brains, specifically neurons and not glia. This was not due to differences in reporter abundance or nonsense-mediated mRNA decay (NMD) surveillance between these tissues. Readthrough rates also varied within neurons, with cholinergic neurons having highest readthrough compared with lowest readthrough rates in dopaminergic neurons. Overall, our data reveal temporal and spatial variation of PTC-mediated readthrough in animals, and suggest that readthrough may be a potential rescue mechanism for PTC-harboring transcripts when the NMD surveillance pathway is inhibited.

Non-AUG Translation Initiation Generates Peroxisomal Isoforms of 6-Phosphogluconate Dehydrogenase in Fungi

Proteins destined for transport to specific organelles usually contain targeting information, which are embedded in their sequence. Many enzymes are required in more than one cellular compartment and different molecular mechanisms are used to achieve dual localization. Here we report a cryptic type 2 peroxisomal targeting signal encoded in the 5′ untranslated region of fungal genes coding for 6-phosphogluconate dehydrogenase (PGD), a key enzyme of the oxidative pentose phosphate pathway. The conservation of the cryptic PTS2 motif suggests a biological function. We observed that translation from a non-AUG start codon generates an N-terminally extended peroxisomal isoform of Ustilago maydis PGD. Non-canonical initiation occurred at the sequence AGG AUU, consisting of two near-cognate start codons in tandem. Taken together, our data reveal non-AUG translation initiation as an additional mechanism to achieve the dual localization of a protein required both in the cytosol and the peroxisomes.

Selective breeding modifies mef2ca mutant incomplete penetrance by tuning the opposing Notch pathway

Deleterious genetic mutations allow developmental biologists to understand how genes control development. However, not all loss of function genetic mutants develop phenotypic changes. Many deleterious mutations only produce a phenotype in a subset of mutant individuals, a phenomenon known as incomplete penetrance. Incomplete penetrance can confound analyses of gene function and our understanding of this widespread phenomenon remains inadequate. To better understand what controls penetrance, we capitalized on the zebrafish mef2ca mutant which produces craniofacial phenotypes with variable penetrance. Starting with a characterized mef2ca loss of function mutant allele, we used classical selective breeding methods to generate zebrafish strains in which mutant-associated phenotypes consistently appear with low or high penetrance. Strikingly, our selective breeding for low penetrance converted the mef2ca mutant allele behavior from homozygous lethal to homozygous viable. Meanwhile, selective breeding for high penetrance converted the mef2ca mutant allele from fully recessive to partially dominant. Comparing the selectively-bred low- and high-penetrance strains revealed that the strains initially respond similarly to the mutation, but then gene expression differences between strains emerge during development. Thus, altered temporal genetic circuitry can manifest through selective pressure to modify mutant penetrance. Specifically, we demonstrate differences in Notch signaling between strains, and further show that experimental manipulation of the Notch pathway phenocopies penetrance changes occurring through selective breeding. This study provides evidence that penetrance is inherited as a liability-threshold trait. Our finding that vertebrate animals can overcome a deleterious mutation by tuning genetic circuitry complements other reported mechanisms of overcoming deleterious mutations such as transcriptional adaptation of compensatory genes, alternative mRNA splicing, and maternal deposition of wild-type transcripts, which are not observed in our system. The selective breeding approach and the resultant genetic circuitry change we uncovered advances and expands our current understanding of genetic and developmental resilience.

Some deleterious gene mutations only affect a subset of genetically mutant animals. This widespread phenomenon, known as mutant incomplete penetrance, complicates discovery of causative gene mutations in both model organisms and human disease. This study utilized the zebrafish mef2ca transcription factor mutant that produces craniofacial skeleton defects with incomplete penetrance. Selectively breeding zebrafish families for low- or high-penetrance mutants for many generations created different zebrafish strains with consistently low or high penetrance. Comparing these strains allowed us to gain insight into the mechanisms that control penetrance. Specifically, genes under the control of mef2ca are initially similarly expressed between the two strains, but differences between strains emerge during development. We found that genetic manipulation of these downstream genes mimics the effects of our selective breeding. Thus, selective breeding for penetrance can change the genetic circuitry downstream of the mutated gene. We propose that small differences in gene circuitry between individuals is one mechanism underlying susceptibility or resilience to genetic mutations.

The preponderance of nonsynonymous A-to-I RNA editing in coleoids is nonadaptive

A-to-I editing enzymatically converts the base adenosine (A) in RNA molecules to inosine (I), which is recognized as guanine (G) in translation. Exceptionally abundant A-to-I editing was recently discovered in the neural tissues of coleoids (octopuses, squids, and cuttlefishes), with a greater fraction of nonsynonymous sites than synonymous sites subject to high levels of editing. Although this phenomenon is thought to indicate widespread adaptive editing, its potential advantage is unknown. Here we propose an alternative, nonadaptive explanation. Specifically, increasing the cellular editing activity permits some otherwise harmful G-to-A nonsynonymous substitutions, because the As are edited to Is at sufficiently high levels. These high editing levels are constrained upon substitutions, resulting in the predominance of nonsynonymous editing at highly edited sites. Our evidence for this explanation suggests that the prevalent nonsynonymous editing in coleoids is generally nonadaptive, as in species with much lower editing activities.

The neural tissues of coleoids have a greater fraction of nonsynonymous sites than synonymous sites subject to high levels of A-to-I RNA editing, a pattern thought to indicate widespread adaptive editing. Here the authors propose and provide evidence for an alternative, nonadaptive explanation.

In eubacteria, unlike eukaryotes, there is no evidence for selection favouring fail-safe 3' additional stop codons

Errors throughout gene expression are likely deleterious, hence genomes are under selection to ameliorate their consequences. Additional stop codons (ASCs) are in-frame nonsense ‘codons’ downstream of the primary stop which may be read by translational machinery should the primary stop have been accidentally read through. Prior evidence in several eukaryotes suggests that ASCs are selected to prevent potentially-deleterious consequences of read-through. We extend this evidence showing that enrichment of ASCs is common but not universal for single cell eukaryotes. By contrast, there is limited evidence as to whether the same is true in other taxa. Here, we provide the first systematic test of the hypothesis that ASCs act as a fail-safe mechanism in eubacteria, a group with high read-through rates. Contra to the predictions of the hypothesis we find: there is paucity, not enrichment, of ASCs downstream; substitutions that degrade stops are more frequent in-frame than out-of-frame in 3’ sequence; highly expressed genes are no more likely to have ASCs than lowly expressed genes; usage of the leakiest primary stop (TGA) in highly expressed genes does not predict ASC enrichment even compared to usage of non-leaky stops (TAA) in lowly expressed genes, beyond downstream codon +1. Any effect at the codon immediately proximal to the primary stop can be accounted for by a preference for a T/U residue immediately following the stop, although if anything, TT- and TC- starting codons are preferred. We conclude that there is no compelling evidence for ASC selection in eubacteria. This presents an unusual case in which the same error could be solved by the same mechanism in eukaryotes and prokaryotes but is not. We discuss two possible explanations: that, owing to the absence of nonsense mediated decay, bacteria may solve read-through via gene truncation and in eukaryotes certain prion states cause raised read-through rates.

In all organisms, gene expression is error-prone. One such error, translational read-through, occurs where the primary stop codon of an expressed gene is missed by the translational machinery. Failure to terminate is likely to be costly, hence genomes are under selection to prevent this from happening. One proposed error-proofing strategy involves in-frame proximal additional stop codons (ASCs) which may act as a ‘fail-safe’ mechanism by providing another opportunity for translation to terminate. There is evidence for ASC enrichment in several eukaryotes. We extend this evidence showing it to be common but not universal in single celled eukaryotes. However, the situation in bacteria is poorly understood, despite bacteria having high read-through rates. Here, we test the fail-safe hypothesis within a broad range of bacteria. To our surprise, we find that not only are ASCs not enriched, but they may even be selected against. This provides evidence for an unusual circumstance where eukaryotes and prokaryotes could solve the same problem the same way but don’t. What are we to make of this? We suggest that if read-through is the problem, ASCs are not necessarily the expected solution. Owing to the absence of nonsense-mediated decay, a process that makes gene truncation in eukaryotes less viable, we propose bacteria may rescue a leaky stop by mutation that creates a new stop upstream. Alternatively, raised read-through rates in some particular conditions in eukaryotes might explain the difference.