HIV-1 modulates the tRNA pool to improve translation efficiency

GCN2 in Viral Defence and the Subversive Tactics Employed by Viruses

The recent SARS-CoV-2 pandemic and associated COVID19 disease illustrates the important role of viral defence mechanisms in ensuring survival and recovery of the host or patient. Viruses absolutely depend on the host's protein synthesis machinery to replicate, meaning that impeding translation is a powerful way to counteract viruses. One major approach used by cells to obstruct protein synthesis is to phosphorylate the alpha subunit of eukaryotic translation initiation factor 2 (eIF2α). Mammals possess four different eIF2α-kinases: PKR, HRI, PEK/PERK, and GCN2. While PKR is currently considered the principal eIF2α-kinase involved in viral defence, the other eIF2α-kinases have also been found to play significant roles. Unsurprisingly, viruses have developed mechanisms to counteract the actions of eIF2α-kinases, or even to exploit them to their benefit. While some of these virulence factors are specific to one eIF2α-kinase, such as GCN2, others target all eIF2α-kinases. This review critically evaluates the current knowledge of viral mechanisms targeting the eIF2α-kinase GCN2. A detailed and in-depth understanding of the molecular mechanisms by which viruses evade host defence mechanisms will help to inform the development of powerful anti-viral measures.

Interferon-Stimulated Genes that Target Retrovirus Translation

The innate immune system, particularly the interferon (IFN) system, constitutes the initial line of defense against viral infections. IFN signaling induces the expression of interferon-stimulated genes (ISGs), and their products frequently restrict viral infection. Retroviruses like the human immunodeficiency viruses and the human T-lymphotropic viruses cause severe human diseases and are targeted by ISG-encoded proteins. Here, we discuss ISGs that inhibit the translation of retroviral mRNAs and thereby retrovirus propagation. The Schlafen proteins degrade cellular tRNAs and rRNAs needed for translation. Zinc Finger Antiviral Protein and RNA-activated protein kinase inhibit translation initiation factors, and Shiftless suppresses translation recoding essential for the expression of retroviral enzymes. We outline common mechanisms that underlie the antiviral activity of multifunctional ISGs and discuss potential antiretroviral therapeutic approaches based on the mode of action of these ISGs.

Proteolytic cleavage and inactivation of the TRMT1 tRNA modification enzyme by SARS-CoV-2 main protease

Nonstructural protein 5 (Nsp5) is the main protease of SARS-CoV-2 that cleaves viral polyproteins into individual polypeptides necessary for viral replication. Here, we show that Nsp5 binds and cleaves human tRNA methyltransferase 1 (TRMT1), a host enzyme required for a prevalent post-transcriptional modification in tRNAs. Human cells infected with SARS-CoV-2 exhibit a decrease in TRMT1 protein levels and TRMT1-catalyzed tRNA modifications, consistent with TRMT1 cleavage and inactivation by Nsp5. Nsp5 cleaves TRMT1 at a specific position that matches the consensus sequence of SARS-CoV-2 polyprotein cleavage sites, and a single mutation within the sequence inhibits Nsp5-dependent proteolysis of TRMT1. The TRMT1 cleavage fragments exhibit altered RNA binding activity and are unable to rescue tRNA modification in TRMT1-deficient human cells. Compared to wild-type human cells, TRMT1-deficient human cells infected with SARS-CoV-2 exhibit reduced levels of intracellular viral RNA. These findings provide evidence that Nsp5-dependent cleavage of TRMT1 and perturbation of tRNA modification patterns contribute to the cellular pathogenesis of SARS-CoV-2 infection.

Post-Transcriptional HIV-1 Latency: A Promising Target for Therapy?

Human Immunodeficiency Virus type 1 (HIV-1) latency represents a significant hurdle in finding a cure for HIV-1 infections, despite tireless research efforts. This challenge is partly attributed to the intricate nature of HIV-1 latency, wherein various host and viral factors participate in multiple physiological processes. While substantial progress has been made in discovering therapeutic targets for HIV-1 transcription, targets for the post-transcriptional regulation of HIV-1 infections have received less attention. However, cumulative evidence now suggests the pivotal contribution of post-transcriptional regulation to the viral latency in both in vitro models and infected individuals. In this review, we explore recent insights on post-transcriptional latency in HIV-1 and discuss the potential of its therapeutic targets, illustrating some host factors that restrict HIV-1 at the post-transcriptional level.

Differential Selection for Translation Efficiency Shapes Translation Machineries in Bacterial Species

Different bacterial species have dramatically different generation times, from 20-30 min in Escherichia coli to about two weeks in Mycobacterium leprae. The translation machinery in a cell needs to synthesize all proteins for a new cell in each generation. The three subprocesses of translation, i.e., initiation, elongation, and termination, are expected to be under stronger selection pressure to optimize in short-generation bacteria (SGB) such as Vibrio natriegens than in the long-generation Mycobacterium leprae. The initiation efficiency depends on the start codon decoded by the initiation tRNA, the optimal Shine-Dalgarno (SD) decoded by the anti-SD (aSD) sequence on small subunit rRNA, and the secondary structure that may embed the initiation signals and prevent them from being decoded. The elongation efficiency depends on the tRNA pool and codon usage. The termination efficiency in bacteria depends mainly on the nature of the stop codon and the nucleotide immediately downstream of the stop codon. By contrasting SGB with long-generation bacteria (LGB), we predict (1) SGB to have more ribosome RNA operons to produce ribosomes, and more tRNA genes for carrying amino acids to ribosomes, (2) SGB to have a higher percentage of genes using AUG as the start codon and UAA as the stop codon than LGB, (3) SGB to exhibit better codon and anticodon adaptation than LGB, and (4) SGB to have a weaker secondary structure near the translation initiation signals than LGB. These differences between SGB and LGB should be more pronounced in highly expressed genes than the rest of the genes. We present empirical evidence in support of these predictions.

Schlafen14 Impairs HIV-1 Expression in a Codon Usage-Dependent Manner

Schlafen (SLFN) is a family of proteins upregulated by type I interferons with a regulatory role in translation. Intriguingly, SLFN14 associates with the ribosome and can degrade rRNA, tRNA, and mRNA in vitro, but a role in translation is still unknown. Ribosomes are important regulatory hubs during translation elongation of mRNAs rich in rare codons. Therefore, we evaluated the potential role of SLFN14 in the expression of mRNAs enriched in rare codons, using HIV-1 genes as a model. We found that, in a variety of cell types, including primary immune cells, SLFN14 regulates the expression of HIV-1 and non-viral genes based on their codon adaptation index, a measurement of the synonymous codon usage bias; consequently, SLFN14 inhibits the replication of HIV-1. The potent inhibitory effect of SLFN14 on the expression of the rare codon-rich transcript HIV-1 Gag was minimized by codon optimization. Mechanistically, we found that the endoribonuclease activity of SLFN14 is required, and that ribosomal RNA degradation is involved. Therefore, we propose that SLFN14 impairs the expression of HIV-1 transcripts rich in rare codons, in a catalytic-dependent manner.

Human Schlafen 11 inhibits influenza A virus production

Schlafen (SLFN) proteins are a subset of interferon-stimulated early response genes with antiviral properties. An antiviral mechanism of SLFN11 was previously demonstrated in human immunodeficiency virus type 1 (HIV-1)-infected cells, and it was shown that SLFN11 inhibited HIV-1 virus production in a codon usage-specific manner. The codon usage patterns of many viruses are vastly different from those of their hosts. The codon usage-specific inhibition of HIV-1 expression by SLFN11 suggests that SLFN11 may be able to inhibit other viruses with a suboptimal codon usage pattern. However, the effect of SLFN11 on the replication of influenza A virus (IAV) has never been reported. The induction of SLFN11 expression was observed upon IAV infection. The reduction of SLFN11 expression also promotes influenza virus replication. Moreover, we found that overexpression of SLFN11 could reduce the expression of a reporter gene with a viral codon usage pattern, and the inhibition of viral hemagglutinin (HA) gene was codon-specific as the expression of codon optimized HA was not affected. These results indicate that SLFN11 inhibits the influenza A virus in a codon-specific manner and that SLFN11 may contribute to innate defense against influenza A viruses.

Different Patterns of Codon Usage and Amino Acid Composition across Primate Lentiviruses

A common feature of the mammalian Lentiviruses (family Retroviridae) is an RNA genome that contains an extremely high frequency of adenine (31.7-38.2%) while being extremely poor in cytosine (13.9-21.2%). Such a biased nucleotide composition has implications for codon usage, causing a striking difference between the frequency of synonymous codons in Lentiviruses and that in their hosts. To test whether primate Lentiviruses present differences in codon and amino acid composition, we assembled a dataset of genome sequences that includes SIV species infecting Old-World monkeys and African apes, HIV-2, and the four groups of HIV-1. Using principal component analysis, we found that HIV-1 shows a significant enrichment in adenine plus thymine in the third synonymous codon position and in adenine and guanine in the first and second nonsynonymous codon positions. Similarly, we observed an enrichment in adenine and in guanine in nonsynonymous first and second codon positions, which affects the amino acid composition of the proteins Gag, Pol, Vif, Vpr, Tat, Rev, Env, and Nef. This result suggests an effect of natural selection in shaping codon usage. Under the hypothesis that the use of synonyms in HIV-1 could reflect adaptation to that of genes expressed in specific cell types, we found a highly significant correlation between codon usage in HIV-1 and monocytes, which was remarkably higher than that with B and T lymphocytes. This finding is in line with the notion that monocytes represent an HIV-1 reservoir in infected patients, and it could help understand how this reservoir is established and maintained.

Schlafen 12 restricts HIV-1 latency reversal by a codon-usage dependent post-transcriptional block in CD4+ T cells

Latency is a major barrier towards virus elimination in HIV-1-infected individuals. Yet, the mechanisms that contribute to the maintenance of HIV-1 latency are incompletely understood. Here we describe the Schlafen 12 protein (SLFN12) as an HIV-1 restriction factor that establishes a post-transcriptional block in HIV-1-infected cells and thereby inhibits HIV-1 replication and virus reactivation from latently infected cells. The inhibitory activity is dependent on the HIV-1 codon usage and on the SLFN12 RNase active sites. Within HIV-1-infected individuals, SLFN12 expression in PBMCs correlated with HIV-1 plasma viral loads and proviral loads suggesting a link with the general activation of the immune system. Using an RNA FISH-Flow HIV-1 reactivation assay, we demonstrate that SLFN12 expression is enriched in infected cells positive for HIV-1 transcripts but negative for HIV-1 proteins. Thus, codon-usage dependent translation inhibition of HIV-1 proteins participates in HIV-1 latency and can restrict the amount of virus release after latency reversal.

Murine Gammaherpesvirus 68 ORF45 Stimulates B2 Retrotransposon and Pre-tRNA Activation in a Manner Dependent on Mitogen-Activated Protein Kinase (MAPK) Signaling

RNA polymerase III (RNAPIII) transcribes a variety of noncoding RNAs, including tRNA (tRNA) and the B2 family of short interspersed nuclear elements (SINEs). B2 SINEs are noncoding retrotransposons that possess tRNA-like promoters and are normally silenced in healthy somatic tissue. Infection with the murine gammaherpesvirus MHV68 induces transcription of both SINEs and tRNAs, in part through the activity of the viral protein kinase ORF36. Here, we identify the conserved MHV68 tegument protein ORF45 as an additional activator of these RNAPIII loci. MHV68 ORF45 and ORF36 form a complex, resulting in an additive induction RNAPIII and increased ORF45 expression. ORF45-induced RNAPIII transcription is dependent on its activation of the extracellular signal-regulated kinase (ERK) mitogen-activated protein kinase (MAPK) signaling pathway, which in turn increases the abundance of the RNAPIII transcription factor Brf1. Other viral and nonviral activators of MAPK/ERK signaling also increase the levels of Brf1 protein, B2 SINE RNA, and tRNA, suggesting that this is a common strategy to increase RNAPIII activity. IMPORTANCE Gammaherpesviral infection alters the gene expression landscape of a host cell, including through the induction of noncoding RNAs transcribed by RNA polymerase III (RNAPIII). Among these are a class of repetitive genes known as retrotransposons, which are normally silenced elements and can copy and spread throughout the genome, and transfer RNAs (tRNAs), which are fundamental components of protein translation machinery. How these loci are activated during infection is not well understood. Here, we identify ORF45 from the model murine gammaherpesvirus MHV68 as a novel activator of RNAPIII transcription. To do so, it engages the MAPK/ERK signaling pathway, which is a central regulator of cellular response to environmental stimuli. Activation of this pathway leads to the upregulation of a key factor required for RNAPIII activity, Brf1. These findings expand our understanding of the regulation and dysregulation of RNAPIII transcription and highlight how viral cooption of key signaling pathways can impact host gene expression.

Utilizing the codon adaptation index to evaluate the susceptibility to HIV-1 and SARS-CoV-2 related coronaviruses in possible target cells in humans

Comprehensive identification of possible target cells for viruses is crucial for understanding the pathological mechanism of virosis. The susceptibility of cells to viruses depends on many factors. Besides the existence of receptors at the cell surface, effective expression of viral genes is also pivotal for viral infection. The regulation of viral gene expression is a multilevel process including transcription, translational initiation and translational elongation. At the translational elongation level, the translational efficiency of viral mRNAs mainly depends on the match between their codon composition and cellular translational machinery (usually referred to as codon adaptation). Thus, codon adaptation for viral ORFs in different cell types may be related to their susceptibility to viruses. In this study, we selected the codon adaptation index (CAI) which is a common codon adaptation-based indicator for assessing the translational efficiency at the translational elongation level to evaluate the susceptibility to two-pandemic viruses (HIV-1 and SARS-CoV-2) of different human cell types. Compared with previous studies that evaluated the infectivity of viruses based on codon adaptation, the main advantage of our study is that our analysis is refined to the cell-type level. At first, we verified the positive correlation between CAI and translational efficiency and strengthened the rationality of our research method. Then we calculated CAI for ORFs of two viruses in various human cell types. We found that compared to high-expression endogenous genes, the CAIs of viral ORFs are relatively low. This phenomenon implied that two kinds of viruses have not been well adapted to translational regulatory machinery in human cells. Also, we indicated that presumptive susceptibility to viruses according to CAI is usually consistent with the results of experimental research. However, there are still some exceptions. Finally, we found that two viruses have different effects on cellular translational mechanisms. HIV-1 decouples CAI and translational efficiency of endogenous genes in host cells and SARS-CoV-2 exhibits increased CAI for its ORFs in infected cells. Our results implied that at least in cases of HIV-1 and SARS-CoV-2, CAI can be regarded as an auxiliary index to assess cells' susceptibility to viruses but cannot be used as the only evidence to identify viral target cells.

"Transfer" of power: The intersection of DNA virus infection and tRNA biology

Transfer RNAs (tRNAs) are at the heart of the molecular biology central dogma, functioning to decode messenger RNAs into proteins. As obligate intracellular parasites, viruses depend on the host translation machinery, including host tRNAs. Thus, the ability of a virus to fine-tune tRNA expression elicits the power to impact the outcome of infection. DNA viruses commonly upregulate the output of RNA polymerase III (Pol III)-dependent transcripts, including tRNAs. Decades after these initial discoveries we know very little about how mature tRNA pools change during viral infection, as tRNA sequencing methodology has only recently reached proficiency. Here, we review perturbation of tRNA biogenesis by DNA virus infection, including an emerging player called tRNA-derived fragments (tRFs). We discuss how tRNA dysregulation shifts the power landscape between the host and virus, highlighting the potential for tRNA-based antivirals as a future therapeutic.

The Ramp Atlas: facilitating tissue and cell-specific ramp sequence analyses through an intuitive web interface

Ramp sequences occur when the average translational efficiency of codons near the 5′ end of highly expressed genes is significantly lower than the rest of the gene sequence, which counterintuitively increases translational efficiency by decreasing downstream ribosomal collisions. Here, we show that the relative codon adaptiveness within different tissues changes the existence of a ramp sequence without altering the underlying genetic code. We present the first comprehensive analysis of tissue and cell type-specific ramp sequences and report 3108 genes with ramp sequences that change between tissues and cell types, which corresponds with increased gene expression within those tissues and cells. The Ramp Atlas (https://ramps.byu.edu/) allows researchers to query precomputed ramp sequences in 18 388 genes across 62 tissues and 66 cell types and calculate tissue-specific ramp sequences from user-uploaded FASTA files through an intuitive web interface. We used The Ramp Atlas to identify seven SARS-CoV-2 genes and seven human SARS-CoV-2 entry factor genes with tissue-specific ramp sequences that may help explain viral proliferation within those tissues. We anticipate that The Ramp Atlas will facilitate personalized and creative tissue-specific ramp sequence analyses for both human and viral genes that will increase our ability to utilize this often-overlooked regulatory region.

Profiling Selective Packaging of Host RNA and Viral RNA Modification in SARS-CoV-2 Viral Preparations

Viruses package host RNAs in their virions which are associated with a range of functions in the viral life cycle. Previous transcriptomic profiling of host RNA packaging mostly focused on retroviruses. Which host RNAs are packaged in other viruses at the transcriptome level has not been thoroughly examined. Here we perform proof-of-concept studies using both small RNA and large RNA sequencing of six different SARS-CoV-2 viral isolates grown on VeroE6 cells to profile host RNAs present in cell free viral preparations and to explore SARS-CoV-2 genomic RNA modifications. We find selective enrichment of specific host transfer RNAs (tRNAs), tRNA fragments and signal recognition particle (SRP) RNA in SARS-CoV-2 viral preparations. Different viral preparations contain the same set of host RNAs, suggesting a common mechanism of packaging. We estimate that a single SARS-CoV-2 particle likely contains up to one SRP RNA and four tRNA molecules. We identify tRNA modification differences between the tRNAs present in viral preparations and those in the uninfected VeroE6 host cells. Furthermore, we find uncharacterized candidate modifications in the SARS-CoV-2 genomic RNA. Our results reveal an under-studied aspect of viral-host interactions that may be explored for viral therapeutics.

Manipulation of RNA polymerase III by Herpes Simplex Virus-1

RNA polymerase III (Pol III) transcribes noncoding RNA, including transfer RNA (tRNA), and is commonly targeted during cancer and viral infection. We find that Herpes Simplex Virus-1 (HSV-1) stimulates tRNA expression 10-fold. Perturbation of host tRNA synthesis requires nuclear viral entry, but not synthesis of specific viral transcripts. tRNA with a specific codon bias were not targeted—rather increased transcription was observed from euchromatic, actively transcribed loci. tRNA upregulation is linked to unique crosstalk between the Pol II and III transcriptional machinery. While viral infection results in depletion of Pol II on host mRNA promoters, we find that Pol II binding to tRNA loci increases. Finally, we report Pol III and associated factors bind the viral genome, which suggests a previously unrecognized role in HSV-1 gene expression. These findings provide insight into mechanisms by which HSV-1 alters the host nuclear environment, shifting key processes in favor of the pathogen.

RNA Polymerase III (Pol III) transcribes non-coding RNA, including tRNAs. Applying different RNA-Seq techniques, Dremel et al. provide the Pol III transcriptional landscape of Herpes simplex virus 1 (HSV-1) infected cells. Infection leads to an increase in tRNA expression from host euchromatin and Pol II re-localization to tRNA loci. They also find that Pol III – associated factors bind to the viral genome.

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.

Cooperative methylation of human tRNA3Lys at positions A58 and U54 drives the early and late steps of HIV-1 replication

Retroviral infection requires reverse transcription, and the reverse transcriptase (RT) uses cellular tRNA as its primer. In humans, the TRMT6-TRMT61A methyltransferase complex incorporates N¹-methyladenosine modification at tRNA position 58 (m¹A58); however, the role of m¹A58 as an RT-stop site during retroviral infection has remained questionable. Here, we constructed TRMT6 mutant cells to determine the roles of m¹A in HIV-1 infection. We confirmed that tRNA₃^Lys m¹A58 was required for in vitro plus-strand strong-stop by RT. Accordingly, infectivity of VSV-G pseudotyped HIV-1 decreased when the virus contained m¹A58-deficient tRNA₃^Lys instead of m¹A58-modified tRNA₃^Lys. In TRMT6 mutant cells, the global protein synthesis rate was equivalent to that of wild-type cells. However, unexpectedly, plasmid-derived HIV-1 expression showed that TRMT6 mutant cells decreased accumulation of HIV-1 capsid, integrase, Tat, Gag, and GagPol proteins without reduction of HIV-1 RNAs in cells, and fewer viruses were produced. Moreover, the importance of 5,2′-O-dimethyluridine at U54 of tRNA₃^Lys as a second RT-stop site was supported by conservation of retroviral genome-tRNA^Lys sequence-complementarity, and TRMT6 was required for efficient 5-methylation of U54. These findings illuminate the fundamental importance of tRNA m¹A58 modification in both the early and late steps of HIV-1 replication, as well as in the cellular tRNA modification network.

Structure of human cytomegalovirus virion reveals host tRNA binding to capsid-associated tegument protein pp150

Under the Baltimore nucleic acid-based virus classification scheme, the herpesvirus human cytomegalovirus (HCMV) is a Class I virus, meaning that it contains a double-stranded DNA genome-and no RNA. Here, we report sub-particle cryoEM reconstructions of HCMV virions at 2.9 Å resolution revealing structures resembling non-coding transfer RNAs (tRNAs) associated with the virion's capsid-bound tegument protein, pp150. Through deep sequencing, we show that these RNA sequences match human tRNAs, and we built atomic models using the most abundant tRNA species. Based on our models, tRNA recruitment is mediated by the electrostatic interactions between tRNA phosphate groups and the helix-loop-helix motif of HCMV pp150. The specificity of these interactions may explain the absence of such tRNA densities in murine cytomegalovirus and other human herpesviruses.

The Functional Consequences of the Novel Ribosomal Pausing Site in SARS-CoV-2 Spike Glycoprotein RNA

The SARS-CoV-2 Spike glycoprotein (S protein) acquired a unique new 4 amino acid -PRRA- insertion sequence at amino acid residues (aa) 681-684 that forms a new furin cleavage site in S protein as well as several new glycosylation sites. We studied various statistical properties of the -PRRA- insertion at the RNA level (CCUCGGCGGGCA). The nucleotide composition and codon usage of this sequence are different from the rest of the SARS-CoV-2 genome. One of such features is two tandem CGG codons, although the CGG codon is the rarest codon in the SARS-CoV-2 genome. This suggests that the insertion sequence could cause ribosome pausing as the result of these rare codons. Due to population variants, the Nextstrain divergence measure of the CCU codon is extremely large. We cannot exclude that this divergence might affect host immune responses/effectiveness of SARS-CoV-2 vaccines, possibilities awaiting further investigation. Our experimental studies show that the expression level of original RNA sequence "wildtype" spike protein is much lower than for codon-optimized spike protein in all studied cell lines. Interestingly, the original spike sequence produces a higher titer of pseudoviral particles and a higher level of infection. Further mutagenesis experiments suggest that this dual-effect insert, comprised of a combination of overlapping translation pausing and furin sites, has allowed SARS-CoV-2 to infect its new host (human) more readily. This underlines the importance of ribosome pausing to allow efficient regulation of protein expression and also of cotranslational subdomain folding.

Translational adaptation of human viruses to the tissues they infect

Viruses need to hijack the translational machinery of the host cell for a productive infection to happen. However, given the dynamic landscape of tRNA pools among tissues, it is unclear whether different viruses infecting different tissues have adapted their codon usage toward their tropism. Here, we collect the coding sequences of 502 human-infecting viruses and determine that tropism explains changes in codon usage. Using the tRNA abundances across 23 human tissues from The Cancer Genome Atlas (TCGA), we build an in silico model of translational efficiency that validates the correspondence of the viral codon usage with the translational machinery of their tropism. For instance, we detect that severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is specifically adapted to the upper respiratory tract and alveoli. Furthermore, this correspondence is specifically defined in early viral proteins. The observed tissue-specific translational efficiency could be useful for the development of antiviral therapies and vaccines.

Non-optimal codon usage preferences of coronaviruses determine their promiscuity for infecting multiple hosts

Circulating animal coronaviruses occasionally infect humans. The SARS-CoV-2 is responsible for the current worldwide outbreak of COVID-19 that has resulted in 2 112 844 deaths as of late January 2021. We compared genetic code preferences in 496 viruses, including 34 coronaviruses and 242 corresponding hosts, to uncover patterns that distinguish single- and 'promiscuous' multiple-host-infecting viruses. Based on a codon usage preference score, promiscuous viruses were shown to significantly employ nonoptimal codons, namely codons that involve 'wobble' binding to anticodons, as compared to single-host viruses. The codon adaptation index (CAI) and the effective number of codons (ENC) were calculated for all viruses and hosts. Promiscuous viruses were less adapted hosts vs single-host viruses (P-value = 4.392e-11). All coronaviruses exploit nonoptimal codons to infect multiple hosts. We found that nonoptimal codon preferences at the beginning of viral coding sequences enhance the translational efficiency of viral proteins within the host. Finally, coronaviruses lack endogenous RNA degradation motifs to a significant degree, thereby increasing viral mRNA burden and infection load. To conclude, we found that promiscuously infecting coronaviruses prefer nonoptimal codon usage to remove degradation motifs from their RNAs and to dramatically increase their viral RNA production rates.

Codon usage analysis of zoonotic coronaviruses reveals lower adaptation to humans by SARS-CoV-2

Since 2002, the world has witnessed major outbreaks of acute respiratory illness by three zoonotic coronaviruses (CoVs), which differ from each other in pathogenicity. Reasons for the lower pathogenicity of SARS-CoV-2 than the other two zoonotic coronaviruses, SARS-CoV and MERS-CoV, are not well understood. We herein compared the codon usage patterns of the three zoonotic CoVs causing severe acute respiratory syndromes and four human-specific CoVs (NL63, 229E, OC43, and HKU1) causing mild diseases. We found that the seven viruses have different codon usages, with SARS-CoV-2 having the lowest effective number of codons (ENC) among the zoonotic CoVs. Human codon adaptation index (CAI) analysis revealed that the CAI value of SARS-CoV-2 is the lowest among the zoonotic CoVs. The ENC and CAI values of SARS-CoV-2 were more similar to those of the less-pathogenic human-specific CoVs. To further investigate adaptive evolution within SARS-CoV-2, we examined codon usage patterns in 3573 genomes of SARS-CoV-2 collected over the initial 4 months of the pandemic. We showed that the ENC values and the CAI values of SARS-CoV-2 were decreasing over the period. The low ENC and CAI values could be responsible for the lower pathogenicity of SARS-CoV-2. While mutational pressure appears to shape codon adaptation in the overall genomes of SARS-CoV-2 and other zoonotic CoVs, the E gene of SARS-CoV-2, which has the highest codon usage bias, appears to be under strong natural selection. Data from the study contribute to our understanding of the pathogenicity and evolution of SARS-CoV-2 in humans.

Human Schlafen 11 exploits codon preference discrimination to attenuate viral protein synthesis of prototype foamy virus (PFV)

Recently, the Schlafen (SLFN) proteins have been identified as a novel interferon-stimulated family with antiviral properties. In this study, we reported that SLFN11 inhibited prototype foamy virus (PFV) replication. Over-expression of human SLFN11 reduced viral production, while knockdown of SLFN11 enhanced viral infectivity. In addition, SLFN11 from cattle and African green monkey also suppressed PFV production. Both the ATPase activity and helicase activity of SLFN11 were required for its inhibitory function. Dephosphorylation activated the antiviral activity of SLFN11. More importantly, SLFN11 inhibited the expression of viral protein, which was rescued by viral gene codon optimization. Together, our results demonstrated that SLFN11 impaired PFV viral protein synthesis by exploiting the distinct codon usage between the virus and the host. These findings further broaden our understanding of the antiviral properties of the SLFN family and the molecular mechanism of PFV latent infection.

Retroviral Restriction Factors and Their Viral Targets: Restriction Strategies and Evolutionary Adaptations

The evolutionary conflict between retroviruses and their vertebrate hosts over millions of years has led to the emergence of cellular innate immune proteins termed restriction factors as well as their viral antagonists. Evidence accumulated in the last two decades has substantially increased our understanding of the elaborate mechanisms utilized by these restriction factors to inhibit retroviral replication, mechanisms that either directly block viral proteins or interfere with the cellular pathways hijacked by the viruses. Analyses of these complex interactions describe patterns of accelerated evolution for these restriction factors as well as the acquisition and evolution of their virus-encoded antagonists. Evidence is also mounting that many restriction factors identified for their inhibition of specific retroviruses have broader antiviral activity against additional retroviruses as well as against other viruses, and that exposure to these multiple virus challenges has shaped their adaptive evolution. In this review, we provide an overview of the restriction factors that interfere with different steps of the retroviral life cycle, describing their mechanisms of action, adaptive evolution, viral targets and the viral antagonists that evolved to counter these factors.

Metabolic Specialization and Codon Preference of Lignocellulolytic Genes in the White Rot Basidiomycete Ceriporiopsis subvermispora

Ceriporiopsis subvermispora is a white-rot fungus with a high specificity towards lignin mineralization when colonizing dead wood or lignocellulosic compounds. Its lignocellulose degrading system is formed by cellulose hydrolytic enzymes, manganese peroxidases, and laccases that catalyze the efficient depolymerization and mineralization of lignocellulose. To determine if this metabolic specialization has modified codon usage of the lignocellulolytic system, improving its adaptation to the fungal translational machine, we analyzed the adaptation to host codon usage (CAI), tRNA pool (tAI, and AAtAI), codon pair bias (CPB), and the number of effective codons (Nc). These indexes were correlated with gene expression of C. subvermispora, in the presence of glucose and Aspen wood. General gene expression was not correlated with the index values. However, in media containing Aspen wood, the induction of expression of lignocellulose-degrading genes, showed significantly (p < 0.001) higher values of CAI, AAtAI, CPB, tAI, and lower values of Nc than non-induced genes. Cellulose-binding proteins and manganese peroxidases presented the highest adaptation values. We also identified an expansion of genes encoding glycine and glutamic acid tRNAs. Our results suggest that the metabolic specialization to use wood as the sole carbon source has introduced a bias in the codon usage of genes involved in lignocellulose degradation. This bias reduces codon diversity and increases codon usage adaptation to the tRNA pool available in C. subvermispora. To our knowledge, this is the first study showing that codon usage is modified to improve the translation efficiency of a group of genes involved in a particular metabolic process.

Evolutionarily missing and conserved tRNA genes in human and avian

Viral infection heavily relies on host transfer RNA (tRNA) for viral RNA decoding. Counterintuitively, not all tRNA species based on anticodon are matched to all 64-triplet codons during evolution. Life solves this problem by cognate tRNA species via wobbling decoding. We found that 14 out of 64 tRNA genes in humans and the main avian species (chicken and duck) were parallelly missing, including 8 tRNA-A₃₄NN and 6 tRNA-G₃₄NN species. By analyzing the conservation of key motifs in tRNA genes, we found that box A and B served as intragenic tRNA promoters were evolutionally conserved among human, chicken, and duck. Thus, decoding viral RNA by similar wobbling strategies and tRNA transcripts may be parallelly used by human, chicken, and duck. We envisioned that many basic mechanisms regarding viral RNA decoding were possibly conserved in these hosts and may consequently promote cross-species infection. Transfer RNAs (tRNAs) are essentially required for gene decoding. Despite the universal nature of genetic codon, not all tRNA genes are common to all organisms. Here, we would like to discuss fundamental problems and possible effects arising from the evolutionarily missing and conserved tRNA genes in human, chicken, and duck (Alkatib et al., 2012; Ou et al., 2019; Rogalski et al., 2008). Among these three organisms, viruses especially the avian influenza virus can crossly infect (Pepin et al., 2010). For multi-host viruses, similar viral RNA decoding strategies may be parallelly used by different hosts. Because viral cross-species infection heavily relies on host tRNAs of different species for viral RNA decoding (Ou et al., 2020; van Weringh et al., 2011). We envisioned that many basic mechanisms regarding viral RNA decoding were possibly conserved in these three hosts and may consequently promote cross-species infection.

High-throughput interrogation of programmed ribosomal frameshifting in human cells

Programmed ribosomal frameshifting (PRF) is the controlled slippage of the translating ribosome to an alternative frame. This process is widely employed by human viruses such as HIV and SARS coronavirus and is critical for their replication. Here, we developed a high-throughput approach to assess the frameshifting potential of a sequence. We designed and tested >12,000 sequences based on 15 viral and human PRF events, allowing us to systematically dissect the rules governing ribosomal frameshifting and discover novel regulatory inputs based on amino acid properties and tRNA availability. We assessed the natural variation in HIV gag-pol frameshifting rates by testing >500 clinical isolates and identified subtype-specific differences and associations between viral load in patients and the optimality of PRF rates. We devised computational models that accurately predict frameshifting potential and frameshifting rates, including subtle differences between HIV isolates. This approach can contribute to the development of antiviral agents targeting PRF.

Emerging Roles of tRNAs in RNA Virus Infections

Viruses rely on the host cell translation machinery for efficient synthesis of their own proteins. Emerging evidence highlights different roles for host transfer RNAs (tRNAs) in the process of virus replication. For instance, different RNA viruses manipulate host tRNA pools to favor viral protein translation. Interestingly, specific host tRNAs are used as reverse transcription primers and are packaged into retroviral virions. Recent data also demonstrate the formation of tRNA-derived fragments (tRFs) upon infection to facilitate viral replication. Here, we comprehensively discuss how RNA viruses exploit distinct aspects of the host tRNA biology for their benefit. In light of the recent advances in the field, we propose that host tRNA-related pathways and mechanisms represent promising cellular targets for the development of novel antiviral strategies.

Genetic and codon usage bias analyses of major capsid protein gene in Ranavirus

The Ranavirus (one genus of Iridovidae family) is an emerging pathogen that infects fish, amphibian, and reptiles, and causes great economical loss and ecological threat to farmed and wild animals globally. The major capsid protein (MCP) has been used as genetic typing marker and as target to design vaccines. Herein, the codon usage pattern of 73 MCP genes of Ranavirus and Lymphocystivirus are studied by calculating effective number of codons (ENC), relative synonymous codon usage (RSCU), codon adaptation index (CAI), and relative codon deoptimization index (RCDI), and similarity index (SiD). The Ranavirus are confirmed to be classified into five groups by using phylogenetic analysis, and varied nucleotide compositions and hierarchical cluster analysis based on RSCU. The results revealed different codon usage patterns among Lymphocystivirus and five groups of Ranavirus. Ranavirus had six over-represented codons ended with G/C nucleotide, while Lymphocystivirus had six over-represented codons ended with A/T nucleotide. A comparative analysis of parameters that define virus and host relatedness in terms of codon usage were analyzed indicated that Amphibian-like ranaviruses (ALRVs) seem to possess lower ENC values and higher CAIs in contrast to other ranaviruses isolated from fishes, and two groups (FV3-like and CMTV-like group) of them had received higher selection pressure from their hosts as having higher relative codon deoptimization index (RCDI) and similarity index (SiD). The correspondence analysis (COA) and Spearman's rank correlation analyses revealed that nucleotide compositions, relative dinucleotide frequency, mutation pressure, and natural translational selection shape the codon usage pattern in MCP genes and the ENC-GC_3S and neutrality plots indicated that the natural selection is the predominant factor. These results contribute to understanding the evolution of Ranavirus and their adaptions to their hosts.

Codon usage similarity between viral and some host genes suggests a codon-specific translational regulation

The codon usage pattern is a specific characteristic of each species; however, the codon usage of all of the genes in a genome is not uniform. Intriguingly, most viruses have codon usage patterns that are vastly different from the optimal codon usage of their hosts. How viral genes with different codon usage patterns are efficiently expressed during a viral infection is unclear. An analysis of the similarity between viral codon usage and the codon usage of the individual genes of a host genome has never been performed. In this study, we demonstrated that the codon usage of human RNA viruses is similar to that of some human genes, especially those involved in the cell cycle. This finding was substantiated by its concordance with previous reports of an upregulation at the protein level of some of these biological processes. It therefore suggests that some suboptimal viral codon usage patterns may actually be compatible with cellular translational machineries in infected conditions.

Translational recoding: canonical translation mechanisms reinterpreted

During canonical translation, the ribosome moves along an mRNA from the start to the stop codon in exact steps of one codon at a time. The collinearity of the mRNA and the protein sequence is essential for the quality of the cellular proteome. Spontaneous errors in decoding or translocation are rare and result in a deficient protein. However, dedicated recoding signals in the mRNA can reprogram the ribosome to read the message in alternative ways. This review summarizes the recent advances in understanding the mechanisms of three types of recoding events: stop-codon readthrough, –1 ribosome frameshifting and translational bypassing. Recoding events provide insights into alternative modes of ribosome dynamics that are potentially applicable to other non-canonical modes of prokaryotic and eukaryotic translation.

The tRNA-associated dysregulation in immune responses and immune diseases

Transfer RNA (tRNA), often considered as a housekeeping molecule, mainly participates in protein translation by transporting amino acids to the ribosome. Nevertheless, accumulating evidence has shown that tRNAs are closely related to various physiological and pathological processes. The proper functioning of the immune system is the key to human health. The aim of this review is to investigate the relationships between tRNAs and the immune system. We detail the biogenesis and structure of tRNAs and summarize the pathogen tRNA-mediated infection and host responses. In addition, we address recent advances in different aspects of tRNA-associated dysregulation in immune responses and immune diseases, such as tRNA molecules, tRNA modifications, tRNA derivatives and tRNA aminoacylation. Therefore, tRNAs play an important role in immune regulation. Although our knowledge of tRNAs in the context of immunity remains, for the most part, unknown, this field deserves in-depth research to provide new ideas for the treatment of immune diseases.

Modulation of HIV-1 Gag/Gag-Pol frameshifting by tRNA abundance

A hallmark of translation in human immunodeficiency virus type 1 (HIV-1) is a –1 programmed ribosome frameshifting event that produces the Gag-Pol fusion polyprotein. The constant Gag to Gag-Pol ratio is essential for the virion structure and infectivity. Here we show that the frameshifting efficiency is modulated by Leu-tRNA^Leu that reads the UUA codon at the mRNA slippery site. This tRNA^Leu isoacceptor is particularly rare in human cell lines derived from T-lymphocytes, the cells that are targeted by HIV-1. When UUA decoding is delayed, the frameshifting follows an alternative route, which maintains the Gag to Gag-Pol ratio constant. A second potential slippery site downstream of the first one is normally inefficient but can also support –1-frameshifting when altered by a compensatory resistance mutation in response to current antiviral drug therapy. Together these different regimes allow the virus to maintain a constant –1-frameshifting efficiency to ensure successful virus propagation.

Interplay between Intrinsic and Innate Immunity during HIV Infection

Restriction factors are antiviral components of intrinsic immunity which constitute a first line of defense by blocking different steps of the human immunodeficiency virus (HIV) replication cycle. In immune cells, HIV infection is also sensed by several pattern recognition receptors (PRRs), leading to type I interferon (IFN-I) and inflammatory cytokines production that upregulate antiviral interferon-stimulated genes (ISGs). Several studies suggest a link between these two types of immunity. Indeed, restriction factors, that are generally interferon-inducible, are able to modulate immune responses. This review highlights recent knowledge of the interplay between restriction factors and immunity inducing antiviral defenses. Counteraction of this intrinsic and innate immunity by HIV viral proteins will also be discussed.

Mechanisms and biomedical implications of -1 programmed ribosome frameshifting on viral and bacterial mRNAs

Some proteins are expressed as a result of a ribosome frameshifting event that is facilitated by a slippery site and downstream secondary structure elements in the mRNA. This review summarizes recent progress in understanding mechanisms of –1 frameshifting in several viral genes, including IBV 1a/1b, HIV‐1 gag‐pol, and SFV 6K, and in Escherichia coli dnaX. The exact frameshifting route depends on the availability of aminoacyl‐tRNAs: the ribosome normally slips into the –1‐frame during tRNA translocation, but can also frameshift during decoding at condition when aminoacyl‐tRNA is in limited supply. Different frameshifting routes and additional slippery sites allow viruses to maintain a constant production of their key proteins. The emerging idea that tRNA pools are important for frameshifting provides new direction for developing antiviral therapies.

Host Restriction Factors and Human Immunodeficiency Virus (HIV-1): A Dynamic Interplay Involving All Phases of the Viral Life Cycle

Mammalian cells have evolved several mechanisms to prevent or block lentiviral infection and spread. Among the innate immune mechanisms, the signaling cascade triggered by type I interferon (IFN) plays a pivotal role in limiting the burden of HIV-1. In the presence of IFN, human cells upregulate the expression of a number of genes, referred to as IFN-stimulated genes (ISGs), many of them acting as antiviral restriction factors (RFs). RFs are dominant proteins that target different essential steps of the viral cycle, thereby providing an early line of defense against the virus. The identification and characterization of RFs have provided unique insights into the molecular biology of HIV-1, further revealing the complex host-pathogen interplay that characterizes the infection. The presence of RFs drove viral evolution, forcing the virus to develop specific proteins to counteract their activity. The knowledge of the mechanisms that prevent viral infection and their viral counterparts may offer new insights to improve current antiviral strategies. This review provides an overview of the RFs targeting HIV-1 replication and the mechanisms that regulate their expression as well as their impact on viral replication and the clinical course of the disease.

The human ortholog of archaeal Pus10 produces pseudouridine 54 in select tRNAs where its recognition sequence contains a modified residue

The nearly conserved U54 of tRNA is mostly converted to a version of ribothymidine (T) in Bacteria and eukaryotes and to a version of pseudouridine (Ψ) in Archaea. Conserved U55 is nearly always modified to Ψ55 in all organisms. Orthologs of TrmA and TruB that produce T54 and Ψ55, respectively, in Bacteria and eukaryotes are absent in Archaea. Pus10 produces both Ψ54 and Ψ55 in Archaea. Pus10 orthologs are found in nearly all sequenced archaeal and most eukaryal genomes, but not in yeast and bacteria. This coincides with the presence of Ψ54 in most archaeal tRNAs and some animal tRNAs, but its absence from yeast and bacteria. Moreover, Ψ54 is found in several tRNAs that function as primers for retroviral DNA synthesis. Previously, no eukaryotic tRNA Ψ54 synthase had been identified. We show here that human Pus10 can produce Ψ54 in select tRNAs, including tRNA^Lys3, the primer for HIV reverse transcriptase. This synthase activity of Pus10 is restricted to the cytoplasm and is distinct from nuclear Pus10, which is known to be involved in apoptosis. The sequence GUUCAm¹AAUC (m¹A is 1-methyladenosine) at position 53-61 of tRNA along with a stable acceptor stem results in maximum Ψ54 synthase activity. This recognition sequence is unique for a Ψ synthase in that it contains another modification. In addition to Ψ54, SF9 cells-derived recombinant human Pus10 can also generate Ψ55, even in tRNAs that do not contain the Ψ54 synthase recognition sequence. This activity may be redundant with that of TruB.

An improved estimation of tRNA expression to better elucidate the coevolution between tRNA abundance and codon usage in bacteria

The degree to which codon usage can be explained by tRNA abundance in bacterial species is often inadequate, partly because differential tRNA abundance is often approximated by tRNA copy numbers. To better understand the coevolution between tRNA abundance and codon usage, we provide a better estimate of tRNA abundance by profiling tRNA mapped reads (tRNA tpm) using publicly available RNA Sequencing data. To emphasize the feasibility of our approach, we demonstrate that tRNA tpm is consistent with tRNA abundances derived from RNA fingerprinting experiments in Escherichia coli, Bacillus subtilis, and Salmonella enterica. Furthermore, we do not observe an appreciable reduction in tRNA sequencing efficiency due to post-transcriptional methylations in the seven bacteria studied. To determine optimal codons, we calculate codon usage in highly and lowly expressed genes determined by protein per transcript. We found that tRNA tpm is sensitive to identify more translationally optimal codons than gene copy number and early tRNA fingerprinting abundances. Additionally, tRNA tpm improves the predictive power of tRNA adaptation index over codon preference. Our results suggest that dependence of codon usage on tRNA availability is not always associated with species growth-rate. Conversely, tRNA availability is better optimized to codon usage in fast-growing than slow-growing species.

Role of host tRNAs and aminoacyl-tRNA synthetases in retroviral replication

The lifecycle of retroviruses and retrotransposons includes a reverse transcription step, wherein dsDNA is synthesized from genomic RNA for subsequent insertion into the host genome. Retroviruses and retrotransposons commonly appropriate major components of the host cell translational machinery, including cellular tRNAs, which are exploited as reverse transcription primers. Nonpriming functions of tRNAs have also been proposed, such as in HIV-1 virion assembly, and tRNA-derived fragments may also be involved in retrovirus and retrotransposon replication. Moreover, host cellular proteins regulate retroviral replication by binding to tRNAs and thereby affecting various steps in the viral lifecycle. For example, in some cases, tRNA primer selection is facilitated by cognate aminoacyl-tRNA synthetases (ARSs), which bind tRNAs and ligate them to their corresponding amino acids, but also have many known nontranslational functions. Multi-omic studies have revealed that ARSs interact with both viral proteins and RNAs and potentially regulate retroviral replication. Here, we review the currently known roles of tRNAs and their derivatives in retroviral and retrotransposon replication and shed light on the roles of tRNA-binding proteins such as ARSs in this process.

Bioinformatics and Translation Elongation

Codon usage depends on mutation bias, tRNA-mediated selection, and the need for high efficiency and accuracy in translation. One codon in a synonymous codon family is often strongly over-used, especially in highly expressed genes, which often leads to a high dN/dS ratio because dS is very small. Many different codon usage indices have been proposed to measure codon usage and codon adaptation. Sense codon could be misread by release factors and stop codons misread by tRNAs, which also contribute to codon usage in rare cases. This chapter outlines the conceptual framework on codon evolution, illustrates codon-specific and gene-specific codon usage indices, and presents their applications. A new index for codon adaptation that accounts for background mutation bias (Index of Translation Elongation) is presented and contrasted with codon adaptation index (CAI) which does not consider background mutation bias. They are used to re-analyze data from a recent paper claiming that translation elongation efficiency matters little in protein production. The reanalysis disproves the claim.

Hijacking of the Ubiquitin/Proteasome Pathway by the HIV Auxiliary Proteins

The ubiquitin-proteasome system (UPS) ensures regulation of the protein pool in the cell by ubiquitination of proteins followed by their degradation by the proteasome. It plays a central role in the cell under normal physiological conditions as well as during viral infections. On the one hand, the UPS can be used by the cell to degrade viral proteins, thereby restricting the viral infection. On the other hand, it can also be subverted by the virus to its own advantage, notably to induce degradation of cellular restriction factors. This makes the UPS a central player in viral restriction and counter-restriction. In this respect, the human immunodeficiency viruses (HIV-1 and 2) represent excellent examples. Indeed, many steps of the HIV life cycle are restricted by cellular proteins, some of which are themselves components of the UPS. However, HIV itself hijacks the UPS to mediate defense against several cellular restriction factors. For example, the HIV auxiliary proteins Vif, Vpx and Vpu counteract specific restriction factors by the recruitment of cellular UPS components. In this review, we describe the interplay between HIV and the UPS to illustrate its role in the restriction of viral infections and its hijacking by viral proteins for counter-restriction.

Codon usage of HIV regulatory genes is not determined by nucleotide composition

Codon usage bias can be a result of either mutational bias or selection for translational efficiency and/or accuracy. Previous data has suggested that nucleotide composition constraint was the main determinant of HIV codon usage, and that nucleotide composition and codon usage were different between the regulatory genes, tat and rev, and other viral genes. It is not clear whether translational selection contributed to the codon usage difference and how nucleotide composition and translational selection interact to determine HIV codon usage. In this study, a model of codon bias due to GC composition with modification for the A-rich third codon position was used to calculate predicted HIV codon frequencies based on its nucleotide composition. The predicted codon usage of each gene was compared with the actual codon frequency. The predicted codon usage based on GC composition matched well with the actual codon frequencies for the structural genes (gag, pol and env). However, the codon usage of the regulatory genes (tat and rev) could not be predicted. Codon usage of the regulatory genes was also relatively unbiased showing the highest effective number of codons (ENC). Moreover, the codon adaptation index (CAI) of the regulatory genes showed better adaptation to human codons when compared to other HIV genes. Therefore, the early expressed genes responsible for regulation of the replication cycle, tat and rev, were more similar to humans in terms of codon usage and GC content than other HIV genes. This may help these genes to be expressed efficiently during the early stages of infection.

Genetic and codon usage bias analyses of polymerase genes of equine influenza virus and its relation to evolution

BACKGROUND

Equine influenza is a major health problem of equines worldwide. The polymerase genes of influenza virus have key roles in virus replication, transcription, transmission between hosts and pathogenesis. Hence, the comprehensive genetic and codon usage bias of polymerase genes of equine influenza virus (EIV) were analyzed to elucidate the genetic and evolutionary relationships in a novel perspective.

RESULTS

The group - specific consensus amino acid substitutions were identified in all polymerase genes of EIVs that led to divergence of EIVs into various clades. The consistent amino acid changes were also detected in the Florida clade 2 EIVs circulating in Europe and Asia since 2007. To study the codon usage patterns, a total of 281,324 codons of polymerase genes of EIV H3N8 isolates from 1963 to 2015 were systemically analyzed. The polymerase genes of EIVs exhibit a weak codon usage bias. The ENc-GC3s and Neutrality plots indicated that natural selection is the major influencing factor of codon usage bias, and that the impact of mutation pressure is comparatively minor. The methods for estimating host imposed translation pressure suggested that the polymerase acidic (PA) gene seems to be under less translational pressure compared to polymerase basic 1 (PB1) and polymerase basic 2 (PB2) genes. The multivariate statistical analysis of polymerase genes divided EIVs into four evolutionary diverged clusters - Pre-divergent, Eurasian, Florida sub-lineage 1 and 2.

CONCLUSIONS

Various lineage specific amino acid substitutions observed in all polymerase genes of EIVs and especially, clade 2 EIVs underwent major variations which led to the emergence of a phylogenetically distinct group of EIVs originating from Richmond/1/07. The codon usage bias was low in all the polymerase genes of EIVs that was influenced by the multiple factors such as the nucleotide compositions, mutation pressure, aromaticity and hydropathicity. However, natural selection was the major influencing factor in defining the codon usage patterns and evolution of polymerase genes of EIVs.

Regulation of human immunodeficiency virus type 1 (HIV-1) mRNA translation

Human immunodeficiency virus type 1 (HIV-1) mRNA translation is a complex process that uses the host translation machinery to synthesise viral proteins. Several mechanisms for HIV-1 mRNA translation initiation have been proposed including (1) cap-dependent, eIF4E-dependent, (2) cap-dependent, cap-binding complex-dependent, (3) internal ribosome entry sites, and (4) ribosome shunting. While these mechanisms promote HIV-1 mRNA translation in the context of in vitro systems and subgenomic constructs, there are substantial knowledge gaps in understanding how they regulate viral protein production in the context of full-length virus infection. In this review, we will summarise the different translation mechanisms used by HIV-1 mRNAs and the challenges in understanding how they regulate protein synthesis during viral infection.

Biased Nucleotide Composition and Differential Codon Usage Pattern in HIV-1 and HIV-2

HIV-1 and HIV-2 are closely related retroviruses with differences in pathogenicity and geographic distribution. HIV-2 infection is associated with slower disease progression and transmission, longer latency period, low or undetectable plasmatic viral loads, and reduced likelihood of progression to AIDS, compared to HIV-1. In this investigation, we analyzed HIV-2 genes and genomes and compared them with that of HIV-1 belonging to various subtypes. Comparative analysis of the effective number of codons (ENC) for each of the nine genes of the two viruses revealed that the tat gene of HIV-2 had a higher ENC value compared to HIV-1 tat, reflecting lower levels of expression of HIV-2 tat. Lower levels of tat protein particularly during the early stages of infection could result in a lower viral load, lower viral set point, and delayed progression of disease in HIV-2-infected individuals compared to HIV-1-infected subjects. Furthermore, the GC3 composition of the regulatory genes of HIV-2 was ≥50%, suggesting a firm effort by these viruses to adapt themselves to evolutionary survival. We hypothesize that differential codon usage could be one of the possible factors that could contribute to the diminished pathogenicity of HIV-2 in the host as compared to HIV-1.

Predictive biophysical modeling and understanding of the dynamics of mRNA translation and its evolution

mRNA translation is the fundamental process of decoding the information encoded in mRNA molecules by the ribosome for the synthesis of proteins. The centrality of this process in various biomedical disciplines such as cell biology, evolution and biotechnology, encouraged the development of dozens of mathematical and computational models of translation in recent years. These models aimed at capturing various biophysical aspects of the process. The objective of this review is to survey these models, focusing on those based and/or validated on real large-scale genomic data. We consider aspects such as the complexity of the models, the biophysical aspects they regard and the predictions they may provide. Furthermore, we survey the central systems biology discoveries reported on their basis. This review demonstrates the fundamental advantages of employing computational biophysical translation models in general, and discusses the relative advantages of the different approaches and the challenges in the field.

Impact of the biased nucleotide composition of viral RNA genomes on RNA structure and codon usage

We are interested in the influence of nucleotide composition on the fundamental characteristics of the virus RNA genome. Most RNA viruses have genomes with a distinct nucleotide composition, e.g. ranging from minimally 12.9 % to maximally 40.3 % (C- and U-count, respectively, in coronavirus HKU). We present a global analysis of diverse virus types, including plus-strand, minus-strand and double-strand RNA viruses, for the impact of this nucleotide preference on the predicted structure of the RNA genome that is packaged in virion particles and on the codon usage in the viral open reading frames. Several virus-specific features will be described, but also some general conclusions were drawn. Without exception, the virus-specific nucleotide bias was enriched in the unpaired, single-stranded regions of the RNA genome, thus creating an even more striking virus-specific signature. We present a simple mechanism that is based on elementary aspects of RNA structure folding to explain this general trend. In general, the nucleotide bias was the major determinant of the virus-specific codon usages, thus limiting a role for codon selection and translational control. We will discuss molecular and evolutionary scenarios that may be responsible for the diverse nucleotide biases of RNA viruses.

Analysis of the synonymous codon usage bias in recently emerged enterovirus D68 strains

Understanding the codon usage pattern of a pathogen and relationship between pathogen and host's codon usage patterns has fundamental and applied interests. Enterovirus D68 (EV-D68) is an emerging pathogen with a potentially high public health significance. In the present study, the synonymous codon usage bias of 27 recently emerged, and historical EV-D68 strains was analyzed. In contrast to previously studied enteroviruses (enterovirus 71 and poliovirus), EV-D68 and human host have a high discrepancy between favored codons. Analysis of viral synonymous codon usage bias metrics, viral nucleotide/dinucleotide compositional parameters, and viral protein properties showed that mutational pressure is more involved in shaping the synonymous codon usage bias of EV-D68 than translation selection. Computation of codon adaptation indices allowed to estimate expression potential of the EV-D68 genome in several commonly used laboratory animals. This approach requires experimental validation and may provide an auxiliary tool for the rational selection of laboratory animals to model emerging viral diseases. Enterovirus D68 genome compositional and codon usage data can be useful for further pathogenesis, animal model, and vaccine design studies.