Ribosomal frameshifting and transcriptional slippage: From genetic steganography and cryptography to adventitious use

Ribosome decision graphs for the representation of eukaryotic RNA translation complexity

The application of ribosome profiling has revealed an unexpected abundance of translation in addition to that responsible for the synthesis of previously annotated protein-coding regions. Multiple short sequences have been found to be translated within single RNA molecules, within both annotated protein-coding and noncoding regions. The biological significance of this translation is a matter of intensive investigation. However, current schematic or annotation-based representations of mRNA translation generally do not account for the apparent multitude of translated regions within the same molecules. They also do not take into account the stochasticity of the process that allows alternative translations of the same RNA molecules by different ribosomes. There is a need for formal representations of mRNA complexity that would enable the analysis of quantitative information on translation and more accurate models for predicting the phenotypic effects of genetic variants affecting translation. To address this, we developed a conceptually novel abstraction that we term ribosome decision graphs (RDGs). RDGs represent translation as multiple ribosome paths through untranslated and translated mRNA segments. We termed the latter "translons." Nondeterministic events, such as initiation, reinitiation, selenocysteine insertion, or ribosomal frameshifting, are then represented as branching points. This representation allows for an adequate representation of eukaryotic translation complexity and focuses on locations critical for translation regulation. We show how RDGs can be used for depicting translated regions and for analyzing genetic variation and quantitative genome-wide data on translation for characterization of regulatory modulators of translation.

Abolished frameshifting for predicted structure-stabilizing SARS-CoV-2 mutants: Implications to alternative conformations and their statistical structural analyses

The SARS-CoV-2 frameshifting element (FSE) has been intensely studied and explored as a therapeutic target for coronavirus diseases including COVID-19. Besides the intriguing virology, this small RNA is known to adopt many length-dependent conformations, as verified by multiple experimental and computational approaches. However, the role these alternative conformations play in the frameshifting mechanism and how to quantify this structural abundance has been an ongoing challenge. Here, we show by DMS and dual-luciferase functional assays that previously predicted FSE mutants (using the RAG graph theory approach) suppress structural transitions and abolish frameshifting. Furthermore, correlated mutation analysis of DMS data by three programs (DREEM, DRACO, and DANCE-MaP) reveals important differences in their estimation of specific RNA conformations, suggesting caution in the interpretation of such complex conformational landscapes. Overall, the abolished frameshifting in three different mutants confirms that all alternative conformations play a role in the pathways of ribosomal transition.

PRFect: a tool to predict programmed ribosomal frameshifts in prokaryotic and viral genomes

BACKGROUND

One of the stranger phenomena that can occur during gene translation is where, as a ribosome reads along the mRNA, various cellular and molecular properties contribute to stalling the ribosome on a slippery sequence and shifting the ribosome into one of the other two alternate reading frames. The alternate frame has different codons, so different amino acids are added to the peptide chain. More importantly, the original stop codon is no longer in-frame, so the ribosome can bypass the stop codon and continue to translate the codons past it. This produces a longer version of the protein, a fusion of the original in-frame amino acids, followed by all the alternate frame amino acids. There is currently no automated software to predict the occurrence of these programmed ribosomal frameshifts (PRF), and they are currently only identified by manual curation.

RESULTS

Here we present PRFect, an innovative machine-learning method for the detection and prediction of PRFs in coding genes of various types. PRFect combines advanced machine learning techniques with the integration of multiple complex cellular properties, such as secondary structure, codon usage, ribosomal binding site interference, direction, and slippery site motif. Calculating and incorporating these diverse properties posed significant challenges, but through extensive research and development, we have achieved a user-friendly approach. The code for PRFect is freely available, open-source, and can be easily installed via a single command in the terminal. Our comprehensive evaluations on diverse organisms, including bacteria, archaea, and phages, demonstrate PRFect's strong performance, achieving high sensitivity, specificity, and an accuracy exceeding 90%. The code for PRFect is freely available and installs with a single terminal command.

CONCLUSION

PRFect represents a significant advancement in the field of PRF detection and prediction, offering a powerful tool for researchers and scientists to unravel the intricacies of programmed ribosomal frameshifting in coding genes.

Alternative translation initiation produces synaptic organizer proteoforms with distinct localization and functions

While previous studies suggest that many mRNAs contain more than one translation initiation site (TIS), the biological significance of most alternative TISs and their corresponding protein isoforms (proteoforms) remains undetermined. Here we show that alternative translation initiation at a CUG and an AUG TIS in neuronal pentraxin receptor (NPR) mRNA produces two proteoforms, and their relative abundance is regulated by both neuronal activity as well as an adjacent RNA secondary structure. Downstream AUG initiation transforms the N-terminal transmembrane domain into a signal peptide, thereby converting NPR to a secreted factor sufficient to promote synaptic clustering of AMPA-type glutamate receptors. Changing the relative proteoform ratio, but not the overall NPR abundance reduces AMPA receptor in parvalbumin (PV)-positive interneurons and induces changes in learning behaviors in mice. In addition to NPR, N-terminal extensions of C1q-like synaptic organizers, mediated by upstream AUU start codons, anchor these otherwise secreted factors to the membrane. Thus, our results uncovered the plasticity of N-terminal signal sequences regulated by alternative TIS usage as a widespread mechanism to diversify protein localization and functions.

Coronavirus takeover of host cell translation and intracellular antiviral response: a molecular perspective

Coronaviruses are a group of related RNA viruses that cause respiratory diseases in humans and animals. Understanding the mechanisms of translation regulation during coronaviral infections is critical for developing antiviral therapies and preventing viral spread. Translation of the viral single-stranded RNA genome in the host cell cytoplasm is an essential step in the life cycle of coronaviruses, which affects the cellular mRNA translation landscape in many ways. Here we discuss various viral strategies of translation control, including how members of the Betacoronavirus genus shut down host cell translation and suppress host innate immune functions, as well as the role of the viral non-structural protein 1 (Nsp1) in the process. We also outline the fate of viral RNA, considering stress response mechanisms triggered in infected cells, and describe how unique viral RNA features contribute to programmed ribosomal -1 frameshifting, RNA editing, and translation shutdown evasion.

Programmable modulation of ribosomal frameshifting by mRNA targeting CRISPR-Cas12a system

Minus 1 programmed ribosomal frameshifting (-1 PRF) is a conserved translational regulation event essential for critical biological processes, including the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) replication. Efficient trans-modulation of the structured RNA element crucial to -1 PRF will endow the therapeutic application. Here, we demonstrate that CRISPR RNA can stimulate efficient -1 PRF. Assembled CRISPR-Cas12a, but not CRISPR-Cas9, complex further enhances -1 PRF efficiency through its higher capacity to stall translating ribosomes. We additionally perform CRISPR-Cas12a targeting to impair the SARS-CoV-2 frameshifting pseudoknot structure via a focused screening. We demonstrate that targeting CRISPR-Cas12a results in more than 70% suppression of -1 PRF in vitro and about 50% suppression in mammalian cells. Our results show the expanded function of the CRISPR-Cas12 system in modulating -1 PRF efficiency through stalling ribosomes and deforming frameshifting stimulatory signals, which could serve as a new strategy for future coronavirus pandemics.

ORFeus: a computational method to detect programmed ribosomal frameshifts and other non-canonical translation events

BACKGROUND

In canonical protein translation, ribosomes initiate translation at a specific start codon, maintain a single reading frame throughout elongation, and terminate at the first in-frame stop codon. However, ribosomal behavior can deviate at each of these steps, sometimes in a programmed manner. Certain mRNAs contain sequence and structural elements that cause ribosomes to begin translation at alternative start codons, shift reading frame, read through stop codons, or reinitiate on the same mRNA. These processes represent important translational control mechanisms that can allow an mRNA to encode multiple functional protein products or regulate protein expression. The prevalence of these events remains uncertain, due to the difficulty of systematic detection.

RESULTS

We have developed a computational model to infer non-canonical translation events from ribosome profiling data.

CONCLUSION

ORFeus identifies known examples of alternative open reading frames and recoding events across different organisms and enables transcriptome-wide searches for novel events.

A highly efficient human cell-free translation system

Cell-free protein synthesis (CFPS) systems enable easy in vitro expression of proteins with many scientific, industrial, and therapeutic applications. Here we present an optimized, highly efficient human cell-free translation system that bypasses many limitations of currently used in vitro systems. This CFPS system is based on extracts from human HEK293T cells engineered to endogenously express GADD34 and K3L proteins, which suppress phosphorylation of translation initiation factor eIF2α. Overexpression of GADD34 and K3L proteins in human cells before cell lysate preparation significantly simplifies lysate preparation. We find that expression of the GADD34 and K3L accessory proteins before cell lysis maintains low levels of phosphorylation of eIF2α in the extracts. During in vitro translation reactions, eIF2α phosphorylation increases moderately in a GCN2-dependent fashion that can be inhibited by GCN2 kinase inhibitors. This new CFPS system should be useful for exploring human translation mechanisms in more physiological conditions outside the cell.

Ribosome Decision Graphs for the Representation of Eukaryotic RNA Translation Complexity

The application of ribosome profiling has revealed an unexpected abundance of translation in addition to that responsible for the synthesis of previously annotated protein-coding regions. Multiple short sequences have been found to be translated within single RNA molecules, both within annotated protein-coding and non-coding regions. The biological significance of this translation is a matter of intensive investigation. However, current schematic or annotation-based representations of mRNA translation generally do not account for the apparent multitude of translated regions within the same molecules. They also do not take into account the stochasticity of the process that allows alternative translations of the same RNA molecules by different ribosomes. There is a need for formal representations of mRNA complexity that would enable the analysis of quantitative information on translation and more accurate models for predicting the phenotypic effects of genetic variants affecting translation. To address this, we developed a conceptually novel abstraction that we term Ribosome Decision Graphs (RDGs). RDGs represent translation as multiple ribosome paths through untranslated and translated mRNA segments. We termed the later 'translons'. Non-deterministic events, such as initiation, re-initiation, selenocysteine insertion or ribosomal frameshifting are then represented as branching points. This representation allows for an adequate representation of eukaryotic translation complexity and focuses on locations critical for translation regulation. We show how RDGs can be used for depicting translated regions, analysis of genetic variation and quantitative genome-wide data on translation for characterisation of regulatory modulators of translation.

The human mitochondrial mRNA structurome reveals mechanisms of gene expression

The mammalian mitochondrial genome encodes thirteen oxidative phosphorylation system proteins, crucial in aerobic energy transduction. These proteins are translated from 9 monocistronic and 2 bicistronic transcripts, whose native structures remain unexplored, leaving fundamental molecular determinants of mitochondrial gene expression unknown. To address this gap, we developed a mitoDMS-MaPseq approach and used DREEM clustering to resolve the native human mitochondrial mt-mRNA structurome. We gained insights into mt-mRNA biology and translation regulatory mechanisms, including a unique programmed ribosomal frameshifting for the ATP8/ATP6 transcript. Furthermore, absence of the mt-mRNA maintenance factor LRPPRC led to a mitochondrial transcriptome structured differently, with specific mRNA regions exhibiting increased or decreased structuredness. This highlights the role of LRPPRC in maintaining mRNA folding to promote mt-mRNA stabilization and efficient translation. In conclusion, our mt-mRNA folding maps reveal novel mitochondrial gene expression mechanisms, serving as a detailed reference and tool for studying them in different physiological and pathological contexts.

The identification and genetic characteristics of Quang Binh virus from field-captured Culex tritaeniorhynchus (Diptera: Culicidae) from Guizhou Province, China

BACKGROUND

Mosquitoes carry a variety of viruses that can cause disease in humans, animals and livestock. Surveys for viruses carried by wild mosquitoes can significantly contribute to surveillance efforts and early detection systems. In addition to mosquito-borne viruses, mosquitoes harbor many insect-specific viruses (ISVs). Quang Binh virus (QBV) is one such example, categorized as an ISV within the Flavivirus genus (family Flaviviridae). QBV has been specifically documented in Vietnam and China, with reports limited to several mosquito species.

METHODS

The homogenate obtained from female mosquitoes was cultured on C6/36 (Aedes albopictus) and BHK-21 (baby hamster kidney) cell lines. Positive cultures were identified by reverse transcription-polymerase chain reaction (RT‒PCR) with taxon- or species-specific primers. Next-generation sequencing was employed to sequence the complete genomes of the identified positive samples. Subsequently, phylogenetic, gene homology, molecular evolutionary and genetic variation analyses were conducted.

RESULT

In 2021, a total of 32,177 adult female mosquitoes were collected from 15 counties in Guizhou Province, China. The predominant mosquito species identified were Culex tritaeniorhynchus, Armigeres subalbatus and Anopheles sinensis. Among the collected mosquitoes, three positive cultures were obtained from Cx. tritaeniorhynchus pools, revealing the presence of Quang Binh virus (QBV) RNA sequences. Phylogenetic analysis indicated that the three Guizhou isolates, along with the prototype isolate from Vietnam, formed distinct branches. These branches were primarily closely related to other QBV isolates reported in China. Comparative analysis revealed a high degree of nucleotide and amino acid homology between the Guizhou isolates and both Vietnamese and other indigenous Chinese isolates. Additionally, nonsynonymous single-nucleotide variants (SNVs) were observed in these strains compared to the QBV prototype strain.

CONCLUSION

This study represents the first report of QBV presences in Cx. tritaeniorhynchus mosquitoes in Guizhou Province, China. Phylogenetic tree analysis showed that the three Guizhou isolates were most closely related to the QBV genes found in China. In addition, the study of the genetic characteristics and variation of this virus provided a deeper understanding of QBV and enriched the baseline data of these insect-specific flaviviruses (ISFVs).

Our current understanding of the toxicity of altered mito-ribosomal fidelity during mitochondrial protein synthesis: What can it tell us about human disease?

Altered mito-ribosomal fidelity is an important and insufficiently understood causative agent of mitochondrial dysfunction. Its pathogenic effects are particularly well-known in the case of mitochondrially induced deafness, due to the existence of the, so called, ototoxic variants at positions 847C (m.1494C) and 908A (m.1555A) of 12S mitochondrial (mt-) rRNA. It was shown long ago that the deleterious effects of these variants could remain dormant until an external stimulus triggered their pathogenicity. Yet, the link from the fidelity defect at the mito-ribosomal level to its phenotypic manifestation remained obscure. Recent work with fidelity-impaired mito-ribosomes, carrying error-prone and hyper-accurate mutations in mito-ribosomal proteins, have started to reveal the complexities of the phenotypic manifestation of mito-ribosomal fidelity defects, leading to a new understanding of mtDNA disease. While much needs to be done to arrive to a clear picture of how defects at the level of mito-ribosomal translation eventually result in the complex patterns of disease observed in patients, the current evidence indicates that altered mito-ribosome function, even at very low levels, may become highly pathogenic. The aims of this review are three-fold. First, we compare the molecular details associated with mito-ribosomal fidelity to those of general ribosomal fidelity. Second, we gather information on the cellular and organismal phenotypes associated with defective translational fidelity in order to provide the necessary grounds for an understanding of the phenotypic manifestation of defective mito-ribosomal fidelity. Finally, the results of recent experiments directly tackling mito-ribosomal fidelity are reviewed and future paths of investigation are discussed.

Nontriplet feature of genetic code in Euplotes ciliates is a result of neutral evolution

The triplet nature of the genetic code is considered a universal feature of known organisms. However, frequent stop codons at internal mRNA positions in Euplotes ciliates ultimately specify ribosomal frameshifting by one or two nucleotides depending on the context, thus posing a nontriplet feature of the genetic code of these organisms. Here, we sequenced transcriptomes of eight Euplotes species and assessed evolutionary patterns arising at frameshift sites. We show that frameshift sites are currently accumulating more rapidly by genetic drift than they are removed by weak selection. The time needed to reach the mutational equilibrium is several times longer than the age of Euplotes and is expected to occur after a several-fold increase in the frequency of frameshift sites. This suggests that Euplotes are at an early stage of the spread of frameshifting in expression of their genome. In addition, we find the net fitness burden of frameshift sites to be noncritical for the survival of Euplotes. Our results suggest that fundamental genome-wide changes such as a violation of the triplet character of genetic code can be introduced and maintained solely by neutral evolution.

A highly efficient human cell-free translation system

Cell-free protein synthesis (CFPS) systems enable easy in vitro expression of proteins with many scientific, industrial, and therapeutic applications. Here we present an optimized, highly efficient human cell-free translation system that bypasses many limitations of currently used in vitro systems. This CFPS system is based on extracts from human HEK293T cells engineered to endogenously express GADD34 and K3L proteins, which suppress phosphorylation of translation initiation factor eIF2α. Overexpression of GADD34 and K3L proteins in human cells significantly simplifies cell lysate preparation. The new CFPS system improves the translation of 5' cap-dependent mRNAs as well as those that use internal ribosome entry site (IRES) mediated translation initiation. We find that expression of the GADD34 and K3L accessory proteins before cell lysis maintains low levels of phosphorylation of eIF2α in the extracts. During in vitro translation reactions, eIF2α phosphorylation increases moderately in a GCN2-dependent fashion that can be inhibited by GCN2 kinase inhibitors. We also find evidence for activation of regulatory pathways related to eukaryotic elongation factor 2 (eEF2) phosphorylation and ribosome quality control in the extracts. This new CFPS system should be useful for exploring human translation mechanisms in more physiological conditions outside the cell.

CRISPR screening reveals a dependency on ribosome recycling for efficient SARS-CoV-2 programmed ribosomal frameshifting and viral replication

During translation of the genomic RNA of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative virus in the COVID-19 pandemic, host ribosomes undergo programmed ribosomal frameshifting (PRF) at a conserved structural element. Although PRF is essential for coronavirus replication, host factors that regulate this process have not yet been identified. Here we perform genome-wide CRISPR-Cas9 knockout screens to identify regulators of SARS-CoV-2 PRF. These screens reveal that loss of ribosome recycling factors markedly decreases frameshifting efficiency and impairs SARS-CoV-2 viral replication. Mutational studies support a model wherein efficient removal of ribosomal subunits at the ORF1a stop codon is required for frameshifting of trailing ribosomes. This dependency upon ribosome recycling is not observed with other non-pathogenic human betacoronaviruses and is likely due to the unique position of the ORF1a stop codon in the SARS clade of coronaviruses. These findings therefore uncover host factors that support efficient SARS-CoV-2 translation and replication.

How much (evil) intelligence can be encoded by 30 kb?

Genomes of most RNA viruses are rarely larger than the size of an average human gene (10-15 kb) and still code for a number of biologically active polypeptides that modify the immune system and metabolism of the host organism in an amazingly complex way. Prolonged coevolution developed tricks by which viruses can dodge many protective mechanisms of the host and lead to the formation of molecular mimicry patterns. Some viruses inhibit the interferon response, interfere with the membrane destroying effects of the activated complement cascade. They can replicate in cellular compartments formed by inner membranes of the cell hiding their characteristic features from diverse pattern recognition receptors. In many cases-and in this respect, the new coronavirus is a champion-they can exploit our own defensive mechanisms to cause serious harm, severe symptoms and frequently deadly disease.

Restriction of Flaviviruses by an Interferon-Stimulated Gene SHFL/C19orf66

Flaviviruses (the genus Flavivirus of the Flaviviridae family) include many arthropod-borne viruses, often causing life-threatening diseases in humans, such as hemorrhaging and encephalitis. Although the flaviviruses have a significant clinical impact, it has become apparent that flavivirus replication is restricted by cellular factors induced by the interferon (IFN) response, which are called IFN-stimulated genes (ISGs). SHFL (shiftless antiviral inhibitor of ribosomal frameshifting) is a novel ISG that inhibits dengue virus (DENV), West Nile virus (WNV), Zika virus (ZIKV), and Japanese encephalitis virus (JEV) infections. Interestingly, SHFL functions as a broad-spectrum antiviral factor exhibiting suppressive activity against various types of RNA and DNA viruses. In this review, we summarize the current understanding of the molecular mechanisms by which SHFL inhibits flavivirus infection and discuss the molecular basis of the inhibitory mechanism using a predicted tertiary structure of SHFL generated by the program AlphaFold2.

Comparative genome analysis of mycobacteria focusing on tRNA and non-coding RNA

BACKGROUND

The Mycobacterium genus encompasses at least 192 named species, many of which cause severe diseases such as tuberculosis. Non-tuberculosis mycobacteria (NTM) can also infect humans and animals. Some are of emerging concern because they show high resistance to commonly used antibiotics while others are used and evaluated in bioremediation or included in anticancer vaccines.

RESULTS

We provide the genome sequences for 114 mycobacterial type strains and together with 130 available mycobacterial genomes we generated a phylogenetic tree based on 387 core genes and supported by average nucleotide identity (ANI) data. The 244 genome sequences cover most of the species constituting the Mycobacterium genus. The genome sizes ranged from 3.2 to 8.1 Mb with an average of 5.7 Mb, and we identified 14 new plasmids. Moreover, mycobacterial genomes consisted of phage-like sequences ranging between 0 and 4.64% dependent on mycobacteria while the number of IS elements varied between 1 and 290. Our data also revealed that, depending on the mycobacteria, the number of tRNA and non-coding (nc) RNA genes differ and that their positions on the chromosome varied. We identified a conserved core set of 12 ncRNAs, 43 tRNAs and 18 aminoacyl-tRNA synthetases among mycobacteria.

CONCLUSIONS

Phages, IS elements, tRNA and ncRNAs appear to have contributed to the evolution of the Mycobacterium genus where several tRNA and ncRNA genes have been horizontally transferred. On the basis of our phylogenetic analysis, we identified several isolates of unnamed species as new mycobacterial species or strains of known mycobacteria. The predicted number of coding sequences correlates with genome size while the number of tRNA, rRNA and ncRNA genes does not. Together these findings expand our insight into the evolution of the Mycobacterium genus and as such they establish a platform to understand mycobacterial pathogenicity, their evolution, antibiotic resistance/tolerance as well as the function and evolution of ncRNA among mycobacteria.

Requirement of a functional ion channel for Sindbis virus glycoprotein transport, CPV-II formation, and efficient virus budding

Many viruses encode ion channel proteins that oligomerize to form hydrophilic pores in membranes of virus-infected cells and the viral membrane in some enveloped viruses. Alphavirus 6K, human immunodeficiency virus type 1 Vpu (HIV-Vpu), influenza A virus M2 (IAV-M2), and hepatitis C virus P7 (HCV-P7) are transmembrane ion channel proteins that play essential roles in virus assembly, budding, and entry. While the oligomeric structures and mechanisms of ion channel activity are well-established for M2 and P7, these remain unknown for 6K. Here we investigated the functional role of the ion channel activity of 6K in alphavirus assembly by utilizing a series of Sindbis virus (SINV) ion channel chimeras expressing the ion channel helix from Vpu or M2 or substituting the entire 6K protein with full-length P7, in cis. We demonstrate that the Vpu helix efficiently complements 6K, whereas M2 and P7 are less efficient. Our results indicate that while SINV is primarily insensitive to the M2 ion channel inhibitor amantadine, the Vpu inhibitor 5-N, N-Hexamethylene amiloride (HMA), significantly reduces SINV release, suggesting that the ion channel activity of 6K similar to Vpu, promotes virus budding. Using live-cell imaging of SINV with a miniSOG-tagged 6K and mCherry-tagged E2, we further demonstrate that 6K and E2 colocalize with the Golgi apparatus in the secretory pathway. To contextualize the localization of 6K in the Golgi, we analyzed cells infected with SINV and SINV-ion channel chimeras using transmission electron microscopy. Our results provide evidence for the first time for the functional role of 6K in type II cytopathic vacuoles (CPV-II) formation. We demonstrate that in the absence of 6K, CPV-II, which originates from the Golgi apparatus, is not detected in infected cells, with a concomitant reduction in the glycoprotein transport to the plasma membrane. Substituting a functional ion channel, M2 or Vpu localizing to Golgi, restores CPV-II production, whereas P7, retained in the ER, is inadequate to induce CPV-II formation. Altogether our results indicate that ion channel activity of 6K is required for the formation of CPV-II from the Golgi apparatus, promoting glycoprotein spike transport to the plasma membrane and efficient virus budding.

Epistasis-Driven Evolution of the SARS-CoV-2 Secondary Structure

Epistasis is an evolutionary phenomenon whereby the fitness effect of a mutation depends on the genetic background in which it arises. A key source of epistasis in an RNA molecule is its secondary structure, which contains functionally important topological motifs held together by hydrogen bonds between Watson–Crick (WC) base pairs. Here we study epistasis in the secondary structure of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) by examining properties of derived alleles arising from substitution mutations at ancestral WC base-paired and unpaired (UP) sites in 15 conserved topological motifs across the genome. We uncover fewer derived alleles and lower derived allele frequencies at WC than at UP sites, supporting the hypothesis that modifications to the secondary structure are often deleterious. At WC sites, we also find lower derived allele frequencies for mutations that abolish base pairing than for those that yield G·U “wobbles,” illustrating that weak base pairing can partially preserve the integrity of the secondary structure. Last, we show that WC sites under the strongest epistatic constraint reside in a three-stemmed pseudoknot motif that plays an essential role in programmed ribosomal frameshifting, whereas those under the weakest epistatic constraint are located in 3’ UTR motifs that regulate viral replication and pathogenicity. Our findings demonstrate the importance of epistasis in the evolution of the SARS-CoV-2 secondary structure, as well as highlight putative structural and functional targets of different forms of natural selection.

Three-Layered Complex Interactions among Capsidless (+)ssRNA Yadokariviruses, dsRNA Viruses, and a Fungus

We have previously discovered a virus neo-lifestyle exhibited by a capsidless positive-sense (+), single-stranded (ss) RNA virus YkV1 (family Yadokariviridae) and an unrelated double-stranded (ds) RNA virus YnV1 (proposed family "Yadonushiviridae") in a phytopathogenic ascomycete, Rosellinia necatrix. YkV1 has been proposed to replicate in the capsid provided by YnV1 as if it were a dsRNA virus and enhance YnV1 replication in return. Recently, viruses related to YkV1 (yadokariviruses) have been isolated from diverse ascomycetous fungi. However, it remains obscure whether such viruses generally show the YkV1-like lifestyle. Here, we identified partner viruses for three distinct yadokariviruses, YkV3, YkV4a, and YkV4b, isolated from R. necatrix that were coinfected with multiple dsRNA viruses phylogenetically distantly related to YnV1. We first established transformants of R. necatrix carrying single yadokarivirus cDNAs and fused them with infectants by single partner candidate dsRNA viruses. Consequently, YkV3 and YkV4s replicated only in the presence of RnMBV3 (family Megabirnaviridae) and RnMTV1 (proposed family "Megatotiviridae"), respectively. The partners were mutually interchangeable between the two YkV4 strains and three RnMTV1 strains but not between other combinations involving YkV1 or YkV3. In contrast to YkV1 enhancing YnV1 accumulation, YkV4s reduced RnMTV1 accumulation to different degrees according to strains. Interestingly, YkV4 rescued the host R. necatrix from impaired growth induced by RnMTV1. YkV3 exerted no apparent effect on its partner (RnMBV3) or host fungus. Overall, we revealed that while yadokariviruses generally require partner dsRNA viruses for replication, each yadokarivirus partners with a different dsRNA virus species in the three diverse families and shows a distinct symbiotic relation in a fungus. IMPORTANCE A capsidless (+)ssRNA virus YkV1 (family Yadokariviridae) highjacks the capsid of an unrelated dsRNA virus YnV1 (proposed family "Yadonushiviridae") in a phytopathogenic ascomycete, while YkV1 trans-enhances YnV1 replication. Herein, we identified the dsRNA virus partners of three yadokariviruses (YkV3, YkV4a, and YkV4b) with genome organization different from YkV1 as being different from YnV1 at the suborder level. Their partners were mutually interchangeable between the two YkV4 strains and three strains of the partner virus RnMTV1 (proposed family "Megatotiviridae") but not between other combinations involving YkV1 or YkV3. Unlike YkV1, YkV4s reduced RnMTV1 accumulation and rescued the host fungus from impaired growth induced by RnMTV1. YkV3 exerted no apparent effect on its partner (RnMBV3, family Megabirnaviridae) or host fungus. These revealed that while each yadokarivirus has a species-specific partnership with a dsRNA virus, yadokariviruses collectively partner extremely diverse dsRNA viruses and show three-layered complex mutualistic/antagonistic interactions in a fungus.

The complexity landscape of viral genomes

BACKGROUND

Viruses are among the shortest yet highly abundant species that harbor minimal instructions to infect cells, adapt, multiply, and exist. However, with the current substantial availability of viral genome sequences, the scientific repertory lacks a complexity landscape that automatically enlights viral genomes' organization, relation, and fundamental characteristics.

RESULTS

This work provides a comprehensive landscape of the viral genome's complexity (or quantity of information), identifying the most redundant and complex groups regarding their genome sequence while providing their distribution and characteristics at a large and local scale. Moreover, we identify and quantify inverted repeats abundance in viral genomes. For this purpose, we measure the sequence complexity of each available viral genome using data compression, demonstrating that adequate data compressors can efficiently quantify the complexity of viral genome sequences, including subsequences better represented by algorithmic sources (e.g., inverted repeats). Using a state-of-the-art genomic compressor on an extensive viral genomes database, we show that double-stranded DNA viruses are, on average, the most redundant viruses while single-stranded DNA viruses are the least. Contrarily, double-stranded RNA viruses show a lower redundancy relative to single-stranded RNA. Furthermore, we extend the ability of data compressors to quantify local complexity (or information content) in viral genomes using complexity profiles, unprecedently providing a direct complexity analysis of human herpesviruses. We also conceive a features-based classification methodology that can accurately distinguish viral genomes at different taxonomic levels without direct comparisons between sequences. This methodology combines data compression with simple measures such as GC-content percentage and sequence length, followed by machine learning classifiers.

CONCLUSIONS

This article presents methodologies and findings that are highly relevant for understanding the patterns of similarity and singularity between viral groups, opening new frontiers for studying viral genomes' organization while depicting the complexity trends and classification components of these genomes at different taxonomic levels. The whole study is supported by an extensive website (https://asilab.github.io/canvas/) for comprehending the viral genome characterization using dynamic and interactive approaches.

A Transfectable Fusagravirus from a Japanese Strain of Cryphonectria carpinicola with Spherical Particles

A novel dsRNA virus (Cryphonectria carpinicola fusagravirus 1, CcFGV1), isolated from a Japanese strain (JS13) of Cryphonectria carpinicola, was thoroughly characterized. The biological comparison of a set of isogenic CcFGV1-infected and -free (JS13VF) strains indicated asymptomatic infection by CcFGV1. The sequence analysis showed that the virus has a two open reading frame (ORF) genome of 9.6 kbp with the RNA-directed RNA polymerase domain encoded by ORF2. The N-terminal sequencing and peptide mass fingerprinting showed an N-terminally processed or degraded product (150 kDa) of the 5′-proximal ORF1-encoded protein (1462 amino acids) to make up the CcFGV1 spherical particles of ~40 nm in diameter. Interestingly, a portion of CcFGV1 dsRNA co-fractionated with a host protein of 70 kDa. The purified CcFGV1 particles were used to transfect protoplasts of JS13VF as well as the standard strain of an experimental model filamentous fungal host Cryphonectria parasitica. CcFGV1 was confirmed to be associated with asymptomatic infection of both fungi. RNA silencing was shown to target the virus in C. parasitica, resulting in reduced CcFGV1 accumulation by comparing the CcFGV1 content between RNA silencing-competent and -deficient strains. These results indicate the transfectability of spherical particles of a fusagravirus associated with asymptomatic infection.

Mechanisms of protein evolution

How do proteins evolve? How do changes in sequence mediate changes in protein structure, and in turn in function? This question has multiple angles, ranging from biochemistry and biophysics to evolutionary biology. This review provides a brief integrated view of some key mechanistic aspects of protein evolution. First, we explain how protein evolution is primarily driven by randomly acquired genetic mutations and selection for function, and how these mutations can even give rise to completely new folds. Then, we also comment on how phenotypic protein variability, including promiscuity, transcriptional and translational errors, may also accelerate this process, possibly via "plasticity-first" mechanisms. Finally, we highlight open questions in the field of protein evolution, with respect to the emergence of more sophisticated protein systems such as protein complexes, pathways, and the emergence of pre-LUCA enzymes.

Mutagenic Analysis of the HIV Restriction Factor Shiftless

The interferon-induced host cell protein shiftless (SFL) was reported to inhibit human immunodeficiency virus (HIV) infection by blocking the –1 programmed ribosomal frameshifting (–1PRF) required for expression of the Gag-Pol polyprotein. However, it is not clear how SFL inhibits –1PRF. To address this question, we focused on a 36 amino acids comprising region (termed required for antiviral activity (RAA)) that is essential for suppression of –1PRF and HIV infection and is missing from SFL short (SFLS), a splice variant of SFL with unknown function. Here, we confirm that SFL, but not SFLS, inhibits HIV –1PRF and show that inhibition is cell-type-independent. Mutagenic and biochemical analyses demonstrated that the RAA region is required for SFL self-interactions and confirmed that it is necessary for ribosome association and binding to the HIV RNA. Analysis of SFL mutants with six consecutive amino-acids-comprising deletions in the RAA region suggests effects on binding to the HIV RNA, complete inhibition of –1PRF, inhibition of Gag-Pol expression, and antiviral activity. In contrast, these amino acids did not affect SFL expression and were partially dispensable for SFL self-interactions and binding to the ribosome. Collectively, our results support the notion that SFL binds to the ribosome and the HIV RNA in order to block –1PRF and HIV infection, and suggest that the multimerization of SFL may be functionally important.

Translational Frameshifting in the chlD Gene Gives a Clue to the Coevolution of the Chlorophyll and Cobalamin Biosyntheses

Today, hundreds of prokaryotic species are able to synthesize chlorophyll and cobalamin (vitamin B12). An important step in the biosynthesis of these coenzymes is the insertion of a metal ion into a porphyrin ring. Namely, Mg-chelatase ChlIDH and aerobic Co-chelatase CobNST are utilized in the chlorophyll and vitamin B12 pathways, respectively. The corresponding subunits of these enzymes have common evolutionary origin. Recently, we have identified a highly conserved frameshifting signal in the chlD gene. This unusual regulatory mechanism allowed production of both the small and the medium chelatase subunits from the same gene. Moreover, the chlD gene appeared early in the evolution and could be at the starting point in the development of the chlorophyll and B12 pathways. Here, we studied the possible coevolution of these two pathways through the analysis of the chelatase genes. To do that, we developed a specialized Web database with comprehensive information about more than 1200 prokaryotic genomes. Further analysis allowed us to split the coevolution of the chlorophyll and B12 pathway into eight distinct stages.

Maximum depth sequencing reveals an ON/OFF replication slippage switch and apparent in vivo selection for bifidobacterial pilus expression

The human gut microbiome, of which the genus Bifidobacterium is a prevalent and abundant member, is thought to sustain and enhance human health. Several surface-exposed structures, including so-called sortase-dependent pili, represent important bifidobacterial gut colonization factors. Here we show that expression of two sortase-dependent pilus clusters of the prototype Bifidobacterium breve UCC2003 depends on replication slippage at an intragenic G-tract, equivalents of which are present in various members of the Bifidobacterium genus. The nature and extent of this slippage is modulated by the host environment. Involvement of such sortase-dependent pilus clusters in microbe-host interactions, including bacterial attachment to the gut epithelial cells, has been shown previously and is corroborated here for one case. Using a Maximum Depth Sequencing strategy aimed at excluding PCR and sequencing errors introduced by DNA polymerase reagents, specific G-tract sequences in B. breve UCC2003 reveal a range of G-tract lengths whose plasticity within the population is functionally utilized. Interestingly, replication slippage is shown to be modulated under in vivo conditions in a murine model. This in vivo modulation causes an enrichment of a G-tract length which appears to allow biosynthesis of these sortase-dependent pili. This work provides the first example of productive replication slippage influenced by in vivo conditions. It highlights the potential for microdiversity generation in “beneficial” gut commensals.

Role of Polyamine-Induced Dimerization of Antizyme in Its Cellular Functions

The polyamines, spermine (Spm) and spermidine (Spd), are important for cell growth and function. Their homeostasis is strictly controlled, and a key downregulator of the polyamine pool is the polyamine-inducible protein, antizyme 1 (OAZ1). OAZ1 inhibits polyamine uptake and targets ornithine decarboxylase (ODC), the rate-limiting enzyme of polyamine biosynthesis, for proteasomal degradation. Here we report, for the first time, that polyamines induce dimerization of mouse recombinant full-length OAZ1, forming an (OAZ1)₂-Polyamine complex. Dimerization could be modulated by functionally active C-methylated spermidine mimetics (MeSpds) by changing the position of the methyl group along the Spd backbone—2-MeSpd was a poor inducer as opposed to 1-MeSpd, 3-MeSpd, and Spd, which were good inducers. Importantly, the ability of compounds to inhibit polyamine uptake correlated with the efficiency of the (OAZ1)₂-Polyamine complex formation. Thus, the (OAZ1)₂-Polyamine complex may be needed to inhibit polyamine uptake. The efficiency of polyamine-induced ribosomal +1 frameshifting of OAZ1 mRNA could also be differentially modulated by MeSpds—2-MeSpd was a poor inducer of OAZ1 biosynthesis and hence a poor downregulator of ODC activity unlike the other MeSpds. These findings offer new insight into the OAZ1-mediated regulation of polyamine homeostasis and provide the chemical tools to study it.

A frameshift in time

The efficiency with which ribosomes shift reading frames when decoding viral RNA may change over the course of an infection.

Ribosome profiling of porcine reproductive and respiratory syndrome virus reveals novel features of viral gene expression

The arterivirus porcine reproductive and respiratory syndrome virus (PRRSV) causes significant economic losses to the swine industry worldwide. Here we apply ribosome profiling (RiboSeq) and parallel RNA sequencing (RNASeq) to characterise the transcriptome and translatome of both species of PRRSV and to analyse the host response to infection. We calculated programmed ribosomal frameshift (PRF) efficiency at both sites on the viral genome. This revealed the nsp2 PRF site as the second known example where temporally regulated frameshifting occurs, with increasing -2 PRF efficiency likely facilitated by accumulation of the PRF-stimulatory viral protein, nsp1β. Surprisingly, we find that PRF efficiency at the canonical ORF1ab frameshift site also increases over time, in contradiction of the common assumption that RNA structure-directed frameshift sites operate at a fixed efficiency. This has potential implications for the numerous other viruses with canonical PRF sites. Furthermore, we discovered several highly translated additional viral ORFs, the translation of which may be facilitated by multiple novel viral transcripts. For example, we found a highly expressed 125-codon ORF overlapping nsp12, which is likely translated from novel subgenomic RNA transcripts that overlap the 3' end of ORF1b. Similar transcripts were discovered for both PRRSV-1 and PRRSV-2, suggesting a potential conserved mechanism for temporally regulating expression of the 3'-proximal region of ORF1b. We also identified a highly translated, short upstream ORF in the 5' UTR, the presence of which is highly conserved amongst PRRSV-2 isolates. These findings reveal hidden complexity in the gene expression programmes of these important nidoviruses.

Thinking Outside the Frame: Impacting Genomes Capacity by Programmed Ribosomal Frameshifting

Translation facilitates the transfer of the genetic information stored in the genome via messenger RNAs to a functional protein and is therefore one of the most fundamental cellular processes. Programmed ribosomal frameshifting is a ubiquitous alternative translation event that is extensively used by viruses to regulate gene expression from overlapping open reading frames in a controlled manner. Recent technical advances in the translation field enabled the identification of precise mechanisms as to how and when ribosomes change the reading frame on mRNAs containing cis-acting signals. Several studies began also to illustrate that trans-acting RNA modulators can adjust the timing and efficiency of frameshifting illuminating that frameshifting can be a dynamically regulated process in cells. Here, we intend to summarize these new findings and emphasize how it fits in our current understanding of PRF mechanisms as previously described.

Minding the message: tactics controlling RNA decay, modification, and translation in virus-infected cells

With their categorical requirement for host ribosomes to translate mRNA, viruses provide a wealth of genetically tractable models to investigate how gene expression is remodeled post-transcriptionally by infection-triggered biological stress. By co-opting and subverting cellular pathways that control mRNA decay, modification, and translation, the global landscape of post-transcriptional processes is swiftly reshaped by virus-encoded factors. Concurrent host cell-intrinsic countermeasures likewise conscript post-transcriptional strategies to mobilize critical innate immune defenses. Here we review strategies and mechanisms that control mRNA decay, modification, and translation in animal virus-infected cells. Besides settling infection outcomes, post-transcriptional gene regulation in virus-infected cells epitomizes fundamental physiological stress responses in health and disease.

Identifying Inhibitors of −1 Programmed Ribosomal Frameshifting in a Broad Spectrum of Coronaviruses

Recurrent outbreaks of novel zoonotic coronavirus (CoV) diseases in recent years have highlighted the importance of developing therapeutics with broad-spectrum activity against CoVs. Because all CoVs use −1 programmed ribosomal frameshifting (−1 PRF) to control expression of key viral proteins, the frameshift signal in viral mRNA that stimulates −1 PRF provides a promising potential target for such therapeutics. To test the viability of this strategy, we explored whether small-molecule inhibitors of −1 PRF in SARS-CoV-2 also inhibited −1 PRF in a range of bat CoVs—the most likely source of future zoonoses. Six inhibitors identified in new and previous screens against SARS-CoV-2 were evaluated against the frameshift signals from a panel of representative bat CoVs as well as MERS-CoV. Some drugs had strong activity against subsets of these CoV-derived frameshift signals, while having limited to no effect on −1 PRF caused by frameshift signals from other viruses used as negative controls. Notably, the serine protease inhibitor nafamostat suppressed −1 PRF significantly for multiple CoV-derived frameshift signals. These results suggest it is possible to find small-molecule ligands that inhibit −1 PRF specifically in a broad spectrum of CoVs, establishing frameshift signals as a viable target for developing pan-coronaviral therapeutics.

Insights from structural studies of the Cardiovirus 2A protein

Cardioviruses are single-stranded RNA viruses of the family Picornaviridae. In addition to being the first example of internal ribosome entry site (IRES) utilization, cardioviruses also employ a series of alternative translation strategies, such as Stop-Go translation and programmed ribosome frameshifting. Here, we focus on cardiovirus 2A protein, which is not only a primary virulence factor, but also exerts crucial regulatory functions during translation, including activation of viral ribosome frameshifting and inhibition of host cap-dependent translation. Only recently, biochemical and structural studies have allowed us to close the gaps in our knowledge of how cardiovirus 2A is able to act in diverse translation-related processes as a novel RNA-binding protein. This review will summarize these findings, which ultimately may lead to the discovery of other RNA-mediated gene expression strategies across a broad range of RNA viruses.

The SARS-CoV2 and mitochondria: the impact on cell fate

Coronavirus infection causes endoplasmic reticulum stress inside the cells, which inhibits protein folding. Prolonged endoplasmic reticulum stress causes an apoptotic process of unfolded protein response-induced cell death. Endoplasmic reticulum stress rapidly induces the activation of mTORC1, responsible for the induction of the IRE1-JNK pathway. IRE1-JNK stands out for its dual nature: pro-apoptotic in the first stage of infection, anti-apoptotic in persistently infected cells. Once penetrated the cells, the virus can deflect the mitochondrial function by implementing both waterfalls pro-apoptotic and anti-apoptotic response. The virus prevents, through Open Reading Frame 9b (ORF-9b) interacting with mitochondria, the response of the type I interferon of the cells affected by the infection and is fundamental for generating an antiviral cellular state. ORF-9b has effects on mitochondrial dynamics, inducing fusion and autophagy and promoting cell survival. The recognition of ORF-9b has made it possible to identify it as a molecular target of some existing potentially effective drugs (Midostaurin and Ruxolitinib). Other drugs, with the same target, are currently being tested. Given the great importance of mitochondria in virus-host interaction, in-depth knowledge of the actors and pathways involved is essential to continue developing new therapeutic strategies against SARS CoV2.

Detection and Characterization of RNA Viruses in Red Macroalgae (Bangiaceae) and Their Food Product (Nori Sheets)

Persistent RNA viruses, which have been suggested to form symbiotic relationships with their hosts, have been reported to occur in eukaryotes, such as plants, fungi, and algae. Based on empirical findings, these viruses may also be present in commercially cultivated macroalgae. Accordingly, the present study aimed to screen red macroalgae (family Bangiaceae conchocelis and Neopyropia yezoensis thallus) and processed nori sheets (N. yezoensis) for persistent RNA viruses using fragmented and primer-ligated dsRNA sequencing (FLDS) and targeted reverse transcription PCR (RT-PCR). A Totiviridae-related virus was detected in the conchocelis of Neoporphyra haitanensis, which is widely cultivated in China, while two Mitoviridae-related viruses were found in several conchocelis samples and all N. yezoensis-derived samples (thallus and nori sheets). Mitoviridae-related viruses in N. yezoensis are widespread among cultivated species and not expected to inhibit host growth. Mitoviridae-related viruses were also detected in several phylogenetically distant species in the family Bangiaceae, which suggests that these viruses persisted and coexist in the family Bangiaceae over a long period of time. The present study is the first to report persistent RNA viruses in nori sheets and their raw materials.

Coordination of -1 programmed ribosomal frameshifting by transcript and nascent chain features revealed by deep mutational scanning

Programmed ribosomal frameshifting (PRF) is a translational recoding mechanism that enables the synthesis of multiple polypeptides from a single transcript. During translation of the alphavirus structural polyprotein, the efficiency of -1PRF is coordinated by a 'slippery' sequence in the transcript, an adjacent RNA stem-loop, and a conformational transition in the nascent polypeptide chain. To characterize each of these effectors, we measured the effects of 4530 mutations on -1PRF by deep mutational scanning. While most mutations within the slip-site and stem-loop reduce the efficiency of -1PRF, the effects of mutations upstream of the slip-site are far more variable. We identify several regions where modifications of the amino acid sequence of the nascent polypeptide impact the efficiency of -1PRF. Molecular dynamics simulations of polyprotein biogenesis suggest the effects of these mutations primarily arise from their impacts on the mechanical forces that are generated by the translocon-mediated cotranslational folding of the nascent polypeptide chain. Finally, we provide evidence suggesting that the coupling between cotranslational folding and -1PRF depends on the translation kinetics upstream of the slip-site. These findings demonstrate how -1PRF is coordinated by features within both the transcript and nascent chain.

Structural and molecular basis for Cardiovirus 2A protein as a viral gene expression switch

Programmed –1 ribosomal frameshifting (PRF) in cardioviruses is activated by the 2A protein, a multi-functional virulence factor that also inhibits cap-dependent translational initiation. Here we present the X-ray crystal structure of 2A and show that it selectively binds to a pseudoknot-like conformation of the PRF stimulatory RNA element in the viral genome. Using optical tweezers, we demonstrate that 2A stabilises this RNA element, likely explaining the increase in PRF efficiency in the presence of 2A. Next, we demonstrate a strong interaction between 2A and the small ribosomal subunit and present a cryo-EM structure of 2A bound to initiated 70S ribosomes. Multiple copies of 2A bind to the 16S rRNA where they may compete for binding with initiation and elongation factors. Together, these results define the structural basis for RNA recognition by 2A, show how 2A-mediated stabilisation of an RNA pseudoknot promotes PRF, and reveal how 2A accumulation may shut down translation during virus infection.

Many RNA viruses employ programmed –1 ribosomal frameshifting (PRF) to expand their coding capacity and optimize production of viral proteins. Here, the authors report structural and biophysical analysis of protein 2A from a cardiovirus, with insights into the mechanism of its PRF-stimulatory function.

Shiftless inhibits flavivirus replication in vitro and is neuroprotective in a mouse model of Zika virus pathogenesis

Flaviviruses such as Zika virus and West Nile virus have the potential to cause severe neuropathology if they invade the central nervous system. The type I interferon response is well characterized as contributing to control of flavivirus-induced neuropathogenesis. However, the interferon-stimulated gene (ISG) effectors that confer these neuroprotective effects are less well studied. Here, we used an ISG expression screen to identify Shiftless (SHFL, C19orf66) as a potent inhibitor of diverse positive-stranded RNA viruses, including multiple members of the Flaviviridae (Zika, West Nile, dengue, yellow fever, and hepatitis C viruses). In cultured cells, SHFL functions as a viral RNA-binding protein that inhibits viral replication at a step after primary translation of the incoming genome. The murine ortholog, Shfl, is expressed constitutively in multiple tissues, including the central nervous system. In a mouse model of Zika virus infection, Shfl ^-/- knockout mice exhibit reduced survival, exacerbated neuropathological outcomes, and increased viral replication in the brain and spinal cord. These studies demonstrate that Shfl is an important antiviral effector that contributes to host protection from Zika virus infection and virus-induced neuropathological disease.

Codon usage bias

Codon usage bias is the preferential or non-random use of synonymous codons, a ubiquitous phenomenon observed in bacteria, plants and animals. Different species have consistent and characteristic codon biases. Codon bias varies not only with species, family or group within kingdom, but also between the genes within an organism. Codon usage bias has evolved through mutation, natural selection, and genetic drift in various organisms. Genome composition, GC content, expression level and length of genes, position and context of codons in the genes, recombination rates, mRNA folding, and tRNA abundance and interactions are some factors influencing codon bias. The factors shaping codon bias may also be involved in evolution of the universal genetic code. Codon-usage bias is critical factor determining gene expression and cellular function by influencing diverse processes such as RNA processing, protein translation and protein folding. Codon usage bias reflects the origin, mutation patterns and evolution of the species or genes. Investigations of codon bias patterns in genomes can reveal phylogenetic relationships between organisms, horizontal gene transfers, molecular evolution of genes and identify selective forces that drive their evolution. Most important application of codon bias analysis is in the design of transgenes, to increase gene expression levels through codon optimization, for development of transgenic crops. The review gives an overview of deviations of genetic code, factors influencing codon usage or bias, codon usage bias of nuclear and organellar genes, computational methods to determine codon usage and the significance as well as applications of codon usage analysis in biological research, with emphasis on plants.

Characterization of Four Novel dsRNA Viruses Isolated from Mucor hiemalis Strains

We previously screened the total nucleic acid extracts of 123 Mucor strains for the presence of dsRNA molecules without further molecular analyses. Here, we characterized five novel dsRNA genomes isolated from four different Mucor hiemalis strains with next-generation sequencing (NGS), namely Mucor hiemalis virus 1a (MhV1a) from WRL CN(M) 122; Mucor hiemalis virus 1b (MhV1b) from NRRL 3624; Mucor hiemalis virus 2 (MhV2) from NRRL 3616; and Mucor hiemalis virus 3 (MhV3) and Mucor hiemalis virus (MhV4) from NRRL 3617 strains. Genomes contain two open reading frames (ORF), which encode the coat protein (CP) and the RNA dependent RNA polymerase (RdRp), respectively. In MhV1a and MhV1b, it is predicted to be translated as a fusion protein via -1 ribosomal frameshift, while in MhV4 via a rare +1 (or-2) ribosomal frameshift. In MhV2 and MhV3, the presence of specific UAAUG pentanucleotide motif points to the fact for coupled translation termination and reinitialization. MhV1a, MhV2, and MhV3 are part of the clade representing the genus Victorivirus, while MhV4 is seated in Totivirus genus clade. The detected VLPs in Mucor strains were from 33 to 36 nm in diameter. Hybridization analysis revealed that the dsRNA molecules of MhV1a-MhV4 hybridized to the corresponding molecules.

Investigating molecular mechanisms of 2A-stimulated ribosomal pausing and frameshifting in Theilovirus

The 2A protein of Theiler's murine encephalomyelitis virus (TMEV) acts as a switch to stimulate programmed -1 ribosomal frameshifting (PRF) during infection. Here, we present the X-ray crystal structure of TMEV 2A and define how it recognises the stimulatory RNA element. We demonstrate a critical role for bases upstream of the originally predicted stem-loop, providing evidence for a pseudoknot-like conformation and suggesting that the recognition of this pseudoknot by beta-shell proteins is a conserved feature in cardioviruses. Through examination of PRF in TMEV-infected cells by ribosome profiling, we identify a series of ribosomal pauses around the site of PRF induced by the 2A-pseudoknot complex. Careful normalisation of ribosomal profiling data with a 2A knockout virus facilitated the identification, through disome analysis, of ribosome stacking at the TMEV frameshifting signal. These experiments provide unparalleled detail of the molecular mechanisms underpinning Theilovirus protein-stimulated frameshifting.

Temperature-Dependent Conformational Evolution of SARS CoV-2 RNA Genome Using Network Analysis

Understanding the dynamics of the SARS CoV-2 RNA genome and its dependence on temperature is necessary to fight the current COVID-19 crisis. Computationally, the handling of large data is a major challenge in the elucidation of the structures of RNA. This work presents network analysis as an important tool to see the conformational evolution and the most dominant structures of the RNA genome at six different temperatures. It effectively distinguished different communities of RNA having structural variation. It is found that at higher temperatures (348 K and above), 80% of the RNA structure is destroyed in both the SPC/E and mTIP3P water models. The thermal denaturation free energy change ΔΔG value calculated for the long-lived structure at higher temperatures of 348 and 363 K ranges from 2.58 to 2.78 kcal/mol for the SPC/E water model, which agrees well with the experimentally reported thermal denaturation free energy range of 2.874 kcal/mol of SARS CoV-NP at normal pH. At higher temperatures, the stability of RNA conformation is found to be due to the existence of non-native base pairs in the SPC/E water model.

Twice exploration of tRNA +1 frameshifting in an elongation cycle of protein synthesis

Inducing tRNA +1 frameshifting to read a quadruplet codon has the potential to incorporate a non-natural amino acid into the polypeptide chain. While this strategy is being considered for genome expansion in biotechnology and bioengineering endeavors, a major limitation is a lack of understanding of where the shift occurs in an elongation cycle of protein synthesis. Here, we use the high-efficiency +1-frameshifting SufB2 tRNA, containing an extra nucleotide in the anticodon loop, to address this question. Physical and kinetic measurements of the ribosome reading frame of SufB2 identify twice exploration of +1 frameshifting in one elongation cycle, with the major fraction making the shift during translocation from the aminoacyl-tRNA binding (A) site to the peptidyl-tRNA binding (P) site and the remaining fraction making the shift within the P site upon occupancy of the A site in the +1-frame. We demonstrate that the twice exploration of +1 frameshifting occurs during active protein synthesis and that each exploration is consistent with ribosomal conformational dynamics that permits changes of the reading frame. This work indicates that the ribosome itself is a determinant of changes of the reading frame and reveals a mechanistic parallel of +1 frameshifting with –1 frameshifting.

Targeting Ribosomal Frameshifting as an Antiviral Strategy: From HIV-1 to SARS-CoV-2

Treatment of HIV-1 has largely involved targeting viral enzymes using a cocktail of inhibitors. However, resistance to these inhibitors and toxicity in the long term have pushed the field to identify new therapeutic targets. To that end, -1 programmed ribosomal frameshifting (-1 PRF) has gained attention as a potential node for therapeutic intervention. In this process, a ribosome moves one nucleotide backward in the course of translating a mRNA, revealing a new reading frame for protein synthesis. In HIV-1, -1 PRF allows the virus to regulate the ratios of enzymatic and structural proteins as needed for correct viral particle assembly. Two RNA structural elements are central to -1 PRF in HIV: a slippery sequence and a highly conserved stable hairpin called the HIV-1 frameshifting stimulatory signal (FSS). Dysregulation of -1 PRF is deleterious for the virus. Thus, -1 PRF is an attractive target for new antiviral development. It is important to note that HIV-1 is not the only virus exploiting -1 PRF for regulating production of its proteins. Coronaviruses, including the COVID-19 pandemic virus SARS-CoV-2, also rely on -1 PRF. In SARS-CoV-2 and other coronaviruses, -1 PRF is required for synthesis of RNA-dependent RNA polymerase and several other nonstructural proteins. Coronaviruses employ a more complex RNA structural element for regulating -1 PRF called a pseudoknot. The purpose of this Account is primarily to review the development of molecules targeting HIV-1 -1 PRF. These approaches are case studies illustrating how the entire pipeline from screening to the generation of high-affinity leads might be implemented. We consider both target-based and function-based screening, with a particular focus on our group's approach beginning with a resin-bound dynamic combinatorial library (RBDCL) screen. We then used rational design approaches to optimize binding affinity, selectivity, and cellular bioavailability. Our tactic is, to the best of our knowledge, the only study resulting in compounds that bind specifically to the HIV-1 FSS RNA and reduce infectivity of laboratory and drug-resistant strains of HIV-1 in human cells. Lessons learned from strategies targeting -1 PRF HIV-1 might provide solutions in the development of antivirals in areas of unmet medical need. This includes the development of new frameshift-altering therapies for SARS-CoV-2, approaches to which are very recently beginning to appear.

Genes with 5' terminal oligopyrimidine tracts preferentially escape global suppression of translation by the SARS-CoV-2 Nsp1 protein

Viruses rely on the host translation machinery to synthesize their own proteins. Consequently, they have evolved varied mechanisms to co-opt host translation for their survival. SARS-CoV-2 relies on a nonstructural protein, Nsp1, for shutting down host translation. However, it is currently unknown how viral proteins and host factors critical for viral replication can escape a global shutdown of host translation. Here, using a novel FACS-based assay called MeTAFlow, we report a dose-dependent reduction in both nascent protein synthesis and mRNA abundance in cells expressing Nsp1. We perform RNA-seq and matched ribosome profiling experiments to identify gene-specific changes both at the mRNA expression and translation levels. We discover that a functionally coherent subset of human genes is preferentially translated in the context of Nsp1 expression. These genes include the translation machinery components, RNA binding proteins, and others important for viral pathogenicity. Importantly, we uncovered a remarkable enrichment of 5' terminal oligo-pyrimidine (TOP) tracts among preferentially translated genes. Using reporter assays, we validated that 5' UTRs from TOP transcripts can drive preferential expression in the presence of Nsp1. Finally, we found that LARP1, a key effector protein in the mTOR pathway, may contribute to preferential translation of TOP transcripts in response to Nsp1 expression. Collectively, our study suggests fine-tuning of host gene expression and translation by Nsp1 despite its global repressive effect on host protein synthesis.

The SARS-CoV-2 Programmed -1 Ribosomal Frameshifting Element Crystal Structure Solved to 2.09 Å Using Chaperone-Assisted RNA Crystallography

The programmed -1 ribosomal frameshifting element (PFSE) of SARS-CoV-2 is a well conserved structured RNA found in all coronaviruses' genomes. By adopting a pseudoknot structure in the presence of the ribosome, the PFSE promotes a ribosomal frameshifting event near the stop codon of the first open reading frame Orf1a during translation of the polyprotein pp1a. Frameshifting results in continuation of pp1a via a new open reading frame, Orf1b, that produces the longer pp1ab polyprotein. Polyproteins pp1a and pp1ab produce nonstructural proteins NSPs 1-10 and NSPs 1-16, respectively, which contribute vital functions during the viral life cycle and must be present in the proper stoichiometry. Both drugs and sequence alterations that affect the stability of the -1 programmed ribosomal frameshifting element disrupt the stoichiometry of the NSPs produced, which compromise viral replication. For this reason, the -1 programmed frameshifting element is considered a promising drug target. Using chaperone assisted RNA crystallography, we successfully crystallized and solved the three-dimensional structure of the PFSE. We observe a three-stem H-type pseudoknot structure with the three stems stacked in a vertical orientation stabilized by two triple base pairs at the stem 1/stem 2 and stem 1/stem 3 junctions. This structure provides a new conformation of PFSE distinct from the bent conformations inferred from midresolution cryo-EM models and provides a high-resolution framework for mechanistic investigations and structure-based drug design.

A Novel Frameshifting Inhibitor Having Antiviral Activity against Zoonotic Coronaviruses

Recent outbreaks of zoonotic coronaviruses, such as Middle East respiratory syndrome coronavirus (MERS-CoV) and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), have caused tremendous casualties and great economic shock. Although some repurposed drugs have shown potential therapeutic efficacy in clinical trials, specific therapeutic agents targeting coronaviruses have not yet been developed. During coronavirus replication, a replicase gene cluster, including RNA-dependent RNA polymerase (RdRp), is alternatively translated via a process called -1 programmed ribosomal frameshift (−1 PRF) by an RNA pseudoknot structure encoded in viral RNAs. The coronavirus frameshifting has been identified previously as a target for antiviral therapy. In this study, the frameshifting efficiencies of MERS-CoV, SARS-CoV and SARS-CoV-2 were determined using an in vitro −1 PRF assay system. Our group has searched approximately 9689 small molecules to identify potential −1 PRF inhibitors. Herein, we found that a novel compound, 2-(5-acetylthiophen-2yl)furo[2,3-b]quinoline (KCB261770), inhibits the frameshifting of MERS-CoV and effectively suppresses viral propagation in MERS-CoV-infected cells. The inhibitory effects of 87 derivatives of furo[2,3-b]quinolines were also examined showing less prominent inhibitory effect when compared to compound KCB261770. We demonstrated that KCB261770 inhibits the frameshifting without suppressing cap-dependent translation. Furthermore, this compound was able to inhibit the frameshifting, to some extent, of SARS-CoV and SARS-CoV-2. Therefore, the novel compound 2-(5-acetylthiophen-2yl)furo[2,3-b]quinoline may serve as a promising drug candidate to interfere with pan-coronavirus frameshifting.