Large-scale genomic analysis of codon usage in dengue virus and evaluation of its phylogenetic dependence

Determination of codon pattern and evolutionary forces acting on genes linked to inflammatory bowel disease

Inflammatory bowel disease (IBD) is an inflammatory disorder of the gastrointestinal tract. The present study attempted to understand the codon usage preferences in genes associated with IBD progression. Compositional analysis, codon usage bias (CUB), Relative synonymous codon usage (RSCU), RNA structure, and expression analysis were performed to obtain a comprehensive picture of codon usage in IBD genes. Compositional analysis of 62 IBD-associated genes revealed that G and T are the most and least abundant nucleotides, respectively. ApG, CpA, and TpG dinucleotides were overrepresented or randomly used, while ApC, CpG, GpT, and TpA dinucleotides were either underrepresented or randomly used in genes related to IBD. The codons influencing the codon usage the most in IBD genes were CGC and AGG. A comparison of codon usage between IBD, and pancreatitis (non-IBD inflammatory disease) indicated that only codon CTG codon usage was significantly different between IBD and pancreatitis. At the same time, there were codons ATA, ACA, CGT, CAA, GTA, CCT, ATT, GCT, CGG, TTG, and CAG for whom codon usage was significantly different for IBD and housekeeping gene sets. The results suggest similar codon usage in at least two inflammatory disorders, IBD and pancreatitis. The analysis helps understand the codon biology, factors affecting gene expression of IBD-associated genes, and the evolution of these genes. The study helps reveal the molecular patterns associated with IBD.

Genetic background of adaptation of Crimean-Congo haemorrhagic fever virus to the different tick hosts

Crimean-Congo haemorrhagic fever orthonairovirus (CCHFV) is a negative-sense, single-stranded RNA virus with a segmented genome and the causative agent of a severe Crimean-Congo haemorrhagic fever (CCHF) disease. The virus is transmitted mainly by tick species in Hyalomma genus but other ticks such as representatives of genera Dermacentor and Rhipicephalus may also be involved in virus life cycle. To improve our understanding of CCHFV adaptation to its tick species, we compared nucleotide composition and codon usage patterns among the all CCHFV strains i) which sequences and other metadata as locality of collection and date of isolation are available in GenBank and ii) which were isolated from in-field collected tick species. These criteria fulfilled 70 sequences (24 coding for S, 23 for M, and 23 for L segment) of virus isolates originating from different representatives of Hyalomma and Rhipicephalus genera. Phylogenetic analyses confirmed that Hyalomma- and Rhipicephalus-originating CCHFV isolates belong to phylogenetically distinct CCHFV clades. Analyses of nucleotide composition among the Hyalomma- and Rhipicephalus-originating CCHFV isolates also showed significant differences, mainly in nucleotides located at the 3rd codon positions indicating changes in codon usage among these lineages. Analyses of codon adaptation index (CAI), effective number of codons (ENC), and other codon usage statistics revealed significant differences between Hyalomma- and Rhipicephalus-isolated CCHFV strains. Despite both sets of strains displayed a higher adaptation to use codons that are preferred by Hyalomma ticks than Rhipicephalus ticks, there were distinct codon usage preferences observed between the two tick species. These findings suggest that over the course of its long co-evolution with tick vectors, CCHFV has optimized its codon usage to efficiently utilize translational resources of Hyalomma species.

Comprehensive Genomics Investigation of Neboviruses Reveals Distinct Codon Usage Patterns and Host Specificity

Neboviruses (NeVs) from the Caliciviridae family have been linked to enteric diseases in bovines and have been detected worldwide. As viruses rely entirely on the cellular machinery of the host for replication, their ability to thrive in a specific host is greatly impacted by the specific codon usage preferences. Here, we systematically analyzed the codon usage bias in NeVs to explore the genetic and evolutionary patterns. Relative Synonymous Codon Usage and Effective Number of Codon analyses indicated a marginally lower codon usage bias in NeVs, predominantly influenced by the nucleotide compositional constraints. Nonetheless, NeVs showed a higher codon usage bias for codons containing G/C at the third codon position. The neutrality plot analysis revealed natural selection as the primary factor that shaped the codon usage bias in both the VP1 (82%) and VP2 (57%) genes of NeVs. Furthermore, the NeVs showed a highly comparable codon usage pattern to bovines, as reflected through Codon Adaptation Index and Relative Codon Deoptimization Index analyses. Notably, yak NeVs showed considerably different nucleotide compositional constraints and mutational pressure compared to bovine NeVs, which appear to be predominantly host-driven. This study sheds light on the genetic mechanism driving NeVs' adaptability, evolution, and fitness to their host species.

Study on the Characteristic Codon Usage Pattern in Porcine Epidemic Diarrhea Virus Genomes and Its Host Adaptation Phenotype

Porcine epidemic diarrhea virus (PEDV), which classified in the genus Alphacoronavirus, family Coronaviridae, is one of the most important pathogens that cause heavy economic losses in pig industry. Although intensive mutation and recombination analysis of PEDV strains were provided, systematic genome analysis were needed to elucidate the evolution mechanism and codon usage adaptation profiles of the pathogen. Here, a comprehensive investigation was carried out to reveal the systematic evolutionary processes of synonymous codon usage and host-adapted evolution phenotype of PEDV genome. We found a low codon usage bias (CUB) in PEDV genome and that nucleotide compositions, natural selection, mutation pressure and geographical diversity shapes the codon usage patterns of PEDV, with natural selection dominated the overall codon usage bias in PEDV than the others. By using the relative codon deoptimization index (RCDI) and similarity index (SiD) analysis, we observed that genotype II PEDV strains showed the highest level of adaptation phenotype to Sus scrofa than another divergent clade. To the best of our knowledge, this is the first comprehensive report elaborating the codon usage and host adaptation of PEDV. The findings offer an insight into our understanding of factors involved in PEDV evolution, adaptation and fitness toward their hosts.

Evolutionary and codon usage preference insights into spike glycoprotein of SARS-CoV-2

Interaction of SARS-CoV-2 spike glycoprotein with the ACE2 cell receptor is very crucial for virus attachment to human cells. Selected mutations in SARS-CoV-2 S-protein are reported to strengthen its binding affinity to mammalian ACE2. The N501T mutation in SARS-CoV-2-CTD furnishes better support to hotspot 353 in comparison with SARS-CoV and shows higher affinity for receptor binding. Recombination analysis exhibited higher recombination events in SARS-CoV-2 strains, irrespective of their geographical origin or hosts. Investigation further supports a common origin among SARS-CoV-2 and its predecessors, SARS-CoV and bat-SARS-like-CoV. The recombination events suggest a constant exchange of genetic material among the co-infecting viruses in possible reservoirs and human hosts before SARS-CoV-2 emerged. Furthermore, a comprehensive analysis of codon usage bias (CUB) in SARS-CoV-2 revealed significant CUB among the S-genes of different beta-coronaviruses governed majorly by natural selection and mutation pressure. Various indices of codon usage of S-genes helped in quantifying its adaptability in other animal hosts. These findings might help in identifying potential experimental animal models for investigating pathogenicity for drugs and vaccine development experiments.

Analysis of codon usage patterns and influencing factors in rice tungro bacilliform virus

Rice tungro bacilliform virus (RTBV) belongs to genus Tungrovirus within the family Caulimoviridae harbors circular double-stranded DNA (dsDNA). Rice tungro disease (RTD) caused by RTBV, responsible for severe rice yield losses in South and Southeast Asia. Here, we performed a systematic evolutionary and codon usage bias (CUB) analysis of RTBV genome sequences. We analysed different bioinformatics techniques to calculate the nucleotide compositions, the relative synonymous codon usage (RSCU), and other indices. The results indicated slightly or low codon usage bias in RTBV isolates. Mutation and natural selection pressures have equally contributed to this low codon usage bias. Additionally, multiple factors such as host, geographical distribution also affect codon usage patterns in RTBV genomes. RSCU analysis revealed that RTBV shows mutation bias and prefers A and U ended codons to code amino acids. Codon usage patterns of RTBV were also found to be influenced by its host. This indicates that RTBV have evolved codon usage patterns that are specific to its host. The findings from this study are expected to increase our understanding of factors leading to viral evolution and fitness with respect to hosts and the environment.

Origin and evolution of emerging Liao ning Virus (genus Seadornavirus, family Reoviridae)

Background

Liao ning virus (LNV) is a member of the genus Seadornavirus, family Reoviridae and has been isolated from kinds of vectors in Asia and Australia. However, there are no systematic studies describe the molecular genetic evolution and migration of LNVs. With the development of bioinformatics, viral genetic data combining the information of virus isolation time and locations could be integrated to infer the virus evolution and spread in nature.

Methods

Here, a phylogenetic and phylogeographic analysis using Bayesian Markov chain Monte Carlo simulations was conducted on the LNVs isolated from a variety of vectors during 1990–2014 to identify the evolution and migration patterns of LNVs.

Results

The results demonstrated that the LNV could be divided into 3 genotypes, of which genotype 1 mainly composed of LNVs isolated from Australia during 1990 to 2014 and the original LNV strain (LNV-NE97–31) isolated from Liaoning province in northern China in 1997, genotype 2 comprised of the isolates all from Xinjiang province in western China and genotype 3 consisted the isolates from Qinghai and Shanxi province of central China. LNVs emerged about 272 years ago and gradually evolved into three lineages in the order genotype 1, genotype 2 and genotype 3. Following phylogeographic analysis, it shows genotype 1 LNVs transmitted from Australia (113°E-153°E,10°S-42°S) to Liaoning province (118°E-125°E,38°N-43°N) in Northeast Asian continent then further spread across the central part of China to western China (75°E-95°E,35°N-50°N).

Conclusion

LNVs were initially isolated from Liaoning province of China in the Northeast Asia, however, the present study revealed that LNVs were first appeared in Australia in the South Pacific region and transmitted to mainland China then rapidly spread across China and evolved three different genotypes. The above results suggested that LNV had the characteristics of long-distance transmission and there were great genetic diversity existed in the LNV population. Notably, current information of 80 strains of LNVs are limited. It is of great importance to strengthen the surveillance of LNVs to explore its real origin in nature and monitoring of the LNVs’ population variation and maintain vigilance to avoid LNV breaking through the species barrier and further clarify its relationship to human and animal infection.

Analysis of Synonymous Codon Usage Bias in Flaviviridae Virus

Background

Flaviviridae viruses are single-stranded, positive-sense RNA viruses, which threat human constantly mediated by mosquitoes, ticks, and sandflies. Considering the recent increase in the prevalence of the family virus and its risk potential, we investigated the codon usage pattern to understand its evolutionary processes and provide some useful data to develop the medications for most of Flaviviridae viruses.

Results

The overall extent of codon usage bias in 65 Flaviviridae viruses is low with the average value of GC contents being 50.5% and the highest value being 55.9%; the lowest value is 40.2%. ENC values of Flaviviridae virus genes vary from 48.75 to 57.83 with a mean value of 55.56. U- and A-ended codons are preferred in the Flaviviridae virus. Correlation analysis shows that the positive correlation between ENC value and GC content at the third nucleotide positions was significant in this family virus. The result of analysis of ENC, neutrality plot analysis, and correlation analysis revealed that codon usage bias of all the viruses was affected mainly by natural selection. Meanwhile, according to correspondence analysis (CoA) based on RSCU and phylogenetic analysis, the Flaviviridae viruses mainly are made up of two groups, Group I (Yellow fever virus, Apoi virus, Tembusu virus, Dengue virus 1, and others) and Group II (West Nile virus lineage 2, Japanese encephalitis virus, Usutu virus, Kedougou virus, and others).

Conclusions

All in, the bias of codon usage pattern is affected not only by compositional constraints but also by natural selection. Phylogenetic analysis also illustrates that codon usage bias of virus can serve as an effective means of evolutionary classification in Flaviviridae virus.

Computational screening of medicinal plant phytochemicals to discover potent pan-serotype inhibitors against dengue virus

Emergence of Dengue as one of the deadliest viral diseases prompts the need for development of effective therapeutic agents. Dengue virus (DV) exists in four different serotypes and infection caused by one serotype predisposes its host to another DV serotype heterotypic re-infection. We undertook virtual ligand screening (VLS) to filter compounds against DV that may inhibit inclusively all of its serotypes. Conserved non-structural DV protein targets such as NS1, NS3/NS2B and NS5, which play crucial role in viral replication, infection cycle and host interaction, were selected for screening of vital antiviral drug leads. A dataset of plant based natural antiviral derivatives was developed. Molecular docking was performed to estimate the spatial affinity of target compounds for the active sites of DV’s NS1, NS3/NS2B and NS5 proteins. The drug likeliness of the screened compounds was followed by ADMET analysis whereas the binding behaviors were further elucidated through molecular dynamics (MD) simulation experiments. VLS screened three potential compounds including Canthin-6-one 9-O-beta-glucopyranoside, Kushenol W and Kushenol K which exhibited optimal binding with all the three conserved DV proteins. This study brings forth novel scaffolds against DV serotypes to serve as lead molecules for further optimization and drug development against all DV serotypes with equal effect against multiple disease causing DV proteins. We therefore anticipate that the insights given in the current study could be regarded valuable towards exploration and development of a broad-spectrum natural anti-dengue therapy.

Analysis of codon usage patterns and influencing factors in Nipah virus

•

Nipah virus (NiV) genes are AT-rich.

•

Codon usage bias of NiV genes is low.

•

Patterns of codon usage bias differ across the genomes of NiV.

•

Both mutation pressure and natural selection influenced codon usage bias of NiV genes.

Codon usage bias (CUB) is the unequal usage of synonymous codons of an amino acid in which some codons are used more often than others and is widely used in understanding molecular biology, genetics, and functional regulation of gene expression. Nipah virus (NiV) is an emerging zoonotic paramyxovirus that causes fatal disease in both humans and animals. NiV was first identified during an outbreak of a disease in Malaysia in 1998 and then occurred periodically since 2001 in India, Bangladesh, and the Philippines. We used bioinformatics tools to analyze the codon usage patterns in a genome-wide manner among 11 genomes of NiV as no work was reported yet. The compositional properties revealed that the overall GC and AT contents were 41.96 and 58.04%, respectively i.e. Nipah virus genes were AT-rich. Correlation analysis between overall nucleotide composition and its 3^rd codon position suggested that both mutation pressure and natural selection might influence the CUB across Nipah genomes. Neutrality plot revealed natural selection might have played a major role while mutation pressure had a minor role in shaping the codon usage bias in NiV genomes.

Codon adaptation biases among sylvatic and urban genotypes of Dengue virus type 2

Dengue virus (DENV) emerged from the sylvatic environment and colonized urban settings, being sustained in a human-Aedes-human transmission chain, mainly by the bites of females of the anthropophilic species Aedes aegypti. Herein, we sought evidence for fine-tuning in viral codon usage, possibly due to viral adaptation to human transmission. We compared the codon adaptation of DENV serotype 2 (DENV-2) genotypes from urban and sylvatic habitats and tried to correlate the findings with key evolutionary determinants. We found that DENV-2 codons of urban and sylvatic genotypes had a higher CAI to humans than to Ae. aegypti. Remarkably, we found no significant differences in codon adaptation to human between urban American/Asian and sylvatic DENV-2 genotypes. Moreover, CAI values were significantly different, when comparing all genotypes to Ae. aegypti codon preferences, with lower values for sylvatic than urban genotypes. In summary, our findings suggest the presence of a molecular signature among the genotypes that circulate in sylvatic and urban environments, and may help explain the trafficking of DENV-2 strains to an urban cycle.

•

DENV-2 codons of all genotypes had a higher CAI to humans than to Ae. Aegypti.

•

CAI values for the sylvatic genotype were the lowest in humans and Ae. Aegypti.

•

Similar CAI values were recovered for the American/Asian and the Sylvatic genotypes.

Analysis of codon usage bias of Crimean-Congo hemorrhagic fever virus and its adaptation to hosts

Crimean-Congo hemorrhagic fever virus (CCHFV) is a negative-sense, single stranded RNA virus with a three-segmented genome that belongs to the genus Nairovirus within the family Bunyaviridae. CCHFV uses Hyalomma ticks as a vector to infect humans with a wide range of clinical signs, from asymptomatic to Zika-like syndrome. Despite significant progress in genomic analyses, the influences of viral relationships with different hosts on overall viral fitness, survival, and evading the host's immune systems remain unknown. To better understand the evolutionary characteristics of CCHFV, we performed a comprehensive analysis of the codon usage pattern in 179 CCHFV strains by calculating the relative synonymous codon usage (RSCU), effective number of codons (ENC), codon adaptation index (CAI), and other indicators. The results indicate that the codon usage bias of CCHFV is relatively low. Several lines of evidence support the hypothesis that a translation selection factor is shaping codon usage pattern in this virus. A correspondence analysis (CA) showed that other factors, such as base composition, aromaticity, and hydrophobicity may also be involved in shaping the codon usage pattern of CCHFV. Additionally, the results from a comparative analysis of RSCU between CCHFV and its hosts suggest that CCHFV tends to evolve codon usage patterns that are comparable to those of its hosts. Furthermore, the selection pressures from Homo sapiens, Bos taurus, and Ovis aries on the CCHFV RSCU patterns were dominant when compared with selection pressure from Hyalomma spp. vectors. Taken together, both natural selection and mutation pressure are important for shaping the codon usage pattern of CCHFV. We believe that such findings will assist researchers in understanding the evolution of CCHFV and its adaptation to its hosts.

Flavivirus Infection Uncouples Translation Suppression from Cellular Stress Responses

As obligate parasites, viruses strictly depend on host cell translation for the production of new progeny, yet infected cells also synthesize antiviral proteins to limit virus infection. Modulation of host cell translation therefore represents a frequent strategy by which viruses optimize their replication and spread. Here we sought to define how host cell translation is regulated during infection of human cells with dengue virus (DENV) and Zika virus (ZIKV), two positive-strand RNA flaviviruses. Polysome profiling and analysis of de novo protein synthesis revealed that flavivirus infection causes potent repression of host cell translation, while synthesis of viral proteins remains efficient. Selective repression of host cell translation was mediated by the DENV polyprotein at the level of translation initiation. In addition, DENV and ZIKV infection suppressed host cell stress responses such as the formation of stress granules and phosphorylation of the translation initiation factor eIF2α (α subunit of eukaryotic initiation factor 2). Mechanistic analyses revealed that translation repression was uncoupled from the disruption of stress granule formation and eIF2α signaling. Rather, DENV infection induced p38-Mnk1 signaling that resulted in the phosphorylation of the eukaryotic translation initiation factor eIF4E and was essential for the efficient production of virus particles. Together, these results identify the uncoupling of translation suppression from the cellular stress responses as a conserved strategy by which flaviviruses ensure efficient replication in human cells.

For efficient production of new progeny, viruses need to balance their dependency on the host cell translation machinery with potentially adverse effects of antiviral proteins produced by the infected cell. To achieve this, many viruses evolved mechanisms to manipulate host cell translation. Here we find that infection of human cells with two major human pathogens, dengue virus (DENV) and Zika virus (ZIKV), leads to the potent repression of host cell translation initiation, while the synthesis of viral protein remains unaffected. Unlike other RNA viruses, these flaviviruses concomitantly suppress host cell stress responses, thereby uncoupling translation suppression from stress granule formation. We identified that the p38-Mnk1 cascade regulating phosphorylation of eIF4E is a target of DENV infection and plays an important role in virus production. Our results define several molecular interfaces by which flaviviruses hijack host cell translation and interfere with stress responses to optimize the production of new virus particles.

Analysis of Synonymous Codon Usage Bias of Zika Virus and Its Adaption to the Hosts

Zika virus (ZIKV) is a mosquito-borne virus (arbovirus) in the family Flaviviridae, and the symptoms caused by ZIKV infection in humans include rash, fever, arthralgia, myalgia, asthenia and conjunctivitis. Codon usage bias analysis can reveal much about the molecular evolution and host adaption of ZIKV. To gain insight into the evolutionary characteristics of ZIKV, we performed a comprehensive analysis on the codon usage pattern in 46 ZIKV strains by calculating the effective number of codons (ENc), codon adaptation index (CAI), relative synonymous codon usage (RSCU), and other indicators. The results indicate that the codon usage bias of ZIKV is relatively low. Several lines of evidence support the hypothesis that translational selection plays a role in shaping the codon usage pattern of ZIKV. The results from a correspondence analysis (CA) indicate that other factors, such as base composition, aromaticity, and hydrophobicity may also be involved in shaping the codon usage pattern of ZIKV. Additionally, the results from a comparative analysis of RSCU between ZIKV and its hosts suggest that ZIKV tends to evolve codon usage patterns that are comparable to those of its hosts. Moreover, selection pressure from Homo sapiens on the ZIKV RSCU patterns was found to be dominant compared with that from Aedes aegypti and Aedes albopictus. Taken together, both natural translational selection and mutation pressure are important for shaping the codon usage pattern of ZIKV. Our findings contribute to understanding the evolution of ZIKV and its adaption to its hosts.

Genome Investigations of Vector Competence in Aedes aegypti to Inform Novel Arbovirus Disease Control Approaches

Dengue (DENV), yellow fever, chikungunya, and Zika virus transmission to humans by a mosquito host is confounded by both intrinsic and extrinsic variables. Besides virulence factors of the individual arboviruses, likelihood of virus transmission is subject to variability in the genome of the primary mosquito vector, Aedes aegypti. The “vectorial capacity” of A. aegypti varies depending upon its density, biting rate, and survival rate, as well as its intrinsic ability to acquire, host and transmit a given arbovirus. This intrinsic ability is known as “vector competence”. Based on whole transcriptome analysis, several genes and pathways have been predicated to have an association with a susceptible or refractory response in A. aegypti to DENV infection. However, the functional genomics of vector competence of A. aegypti is not well understood, primarily due to lack of integrative approaches in genomic or transcriptomic studies. In this review, we focus on the present status of genomics studies of DENV vector competence in A. aegypti as limited information is available relative to the other arboviruses. We propose future areas of research needed to facilitate the integration of vector and virus genomics and environmental factors to work towards better understanding of vector competence and vectorial capacity in natural conditions.

Analysis of genotype diversity and evolution of Dengue virus serotype 2 using complete genomes

BACKGROUND

Dengue is one of the most common arboviral diseases prevalent worldwide and is caused by Dengue viruses (genus Flavivirus, family Flaviviridae). There are four serotypes of Dengue Virus (DENV-1 to DENV-4), each of which is further subdivided into distinct genotypes. DENV-2 is frequently associated with severe dengue infections and epidemics. DENV-2 consists of six genotypes such as Asian/American, Asian I, Asian II, Cosmopolitan, American and sylvatic. Comparative genomic study was carried out to infer population structure of DENV-2 and to analyze the role of evolutionary and spatiotemporal factors in emergence of diversifying lineages.

METHODS

Complete genome sequences of 990 strains of DENV-2 were analyzed using Bayesian-based population genetics and phylogenetic approaches to infer genetically distinct lineages. The role of spatiotemporal factors, genetic recombination and selection pressure in the evolution of DENV-2 is examined using the sequence-based bioinformatics approaches.

RESULTS

DENV-2 genetic structure is complex and consists of fifteen subpopulations/lineages. The Asian/American genotype is observed to be diversified into seven lineages. The Asian I, Cosmopolitan and sylvatic genotypes were found to be subdivided into two lineages, each. The populations of American and Asian II genotypes were observed to be homogeneous. Significant evidence of episodic positive selection was observed in all the genes, except NS4A. Positive selection operational on a few codons in envelope gene confers antigenic and lineage diversity in the American strains of Asian/American genotype. Selection on codons of non-structural genes was observed to impact diversification of lineages in Asian I, cosmopolitan and sylvatic genotypes. Evidence of intra/inter-genotype recombination was obtained and the uncertainty in classification of recombinant strains was resolved using the population genetics approach.

DISCUSSION

Complete genome-based analysis revealed that the worldwide population of DENV-2 strains is subdivided into fifteen lineages. The population structure of DENV-2 is spatiotemporal and is shaped by episodic positive selection and recombination. Intra-genotype diversity was observed in four genotypes (Asian/American, Asian I, cosmopolitan and sylvatic). Episodic positive selection on envelope and non-structural genes translates into antigenic diversity and appears to be responsible for emergence of strains/lineages in DENV-2 genotypes. Understanding of the genotype diversity and emerging lineages will be useful to design strategies for epidemiological surveillance and vaccine design.

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.

Nucleotide composition of the Zika virus RNA genome and its codon usage

BACKGROUND

RNA viruses have genomes with a distinct nucleotide composition and codon usage. We present the global characteristics of the RNA genome of Zika virus (ZIKV), an emerging pathogen within the Flavivirus genus. ZIKV was first isolated in 1947 in Uganda, caused a widespread epidemic in South and Central America and the Caribbean in 2015 and has recently been associated with microcephaly in newborns.

METHODS

The nearly 11 kb positive-stranded RNA genome of ZIKV was analyzed for its nucleotide composition, also in the context of the folded RNA molecule. Nucleotide trends were investigated along the genome length by skew analyses and we analyzed the codons used for translation of the ZIKV proteins.

RESULTS

ZIKV RNA has a biased nucleotide composition in being purine-rich and pyrimidine-poor. This preference for purines is a general characteristic of the mosquito-borne and tick-borne flaviviruses. The virus-specific nucleotide bias is further enriched in the unpaired, single-stranded regions of the structured ZIKV RNA genome, thus further imposing this ZIKV-specific signature. The codons used for translation of the ZIKV proteins is also unusual, but we show that it is the underlying bias in nucleotide composition of the viral RNA that largely dictates these codon preferences.

CONCLUSIONS

The ZIKV RNA genome has a biased nucleotide composition that dictates the codon usage of this flavivirus. We discuss the evolutionary scenarios and molecular mechanisms that may be responsible for these distinctive ZIKV RNA genome features.

Characterization of codon usage pattern and influencing factors in Japanese encephalitis virus

Recently, several outbreaks of Japanese encephalitis (JE), caused by Japanese encephalitis virus (JEV), have been reported and it has become cause of concern across the world. In this study, detailed analysis of JEV codon usage pattern was performed. The relative synonymous codon usage (RSCU) values along with mean effective number of codons (ENC) value of 55.30 indicated the presence of low codon usages bias in JEV. The effect of mutational pressure on codon usage bias was confirmed by significant correlations of A3s, U3s, G3s, C3s, GC3s, ENC values, with overall nucleotide contents (A%, U%, G%, C%, and GC%). The correlation analysis of A3s, U3s, G3s, C3s, GC3s, with axis values of correspondence analysis (CoA) further confirmed the role of mutational pressure. However, the correlation analysis of Gravy values and Aroma values with A3s, U3s, G3s, C3s, and GC3s, indicated the presence of natural selection on codon usage bias in addition to mutational pressure. The natural selection was further confirmed by codon adaptation index (CAI) analysis. Additionally, relative dinucleotide frequencies, geographical distribution, and evolutionary processes also influenced the codon usage pattern to some extent.

Genome-wide analysis of synonymous codon usage in Huaiyangshan virus and other bunyaviruses

Huaiyangshan virus (HYSV) is a newly discovered bunyavirus, which is transmitted by ticks and causes hemorrhagic fever-like illness in human. The interplay of codon usage among viruses and their hosts is expected to affect viral survival, evasion from host's immune system and evolution. However, little is known about the codon usage in HYSV genome. In the present study, we analyzed synonymous codon usage in 120 available full-length HYSV sequences and performed a comparative analysis of synonymous codon usage patterns in HYSV and 42 other bunyaviruses. The relative synonymous codon usage (RSCU) analysis showed that the preferred synonymous codons were G/C-ended. A comparative analysis of RSCU between HYSV and its hosts reflected that codon usage patterns of HYSV were mostly coincident with that of its hosts. Our data suggested that although mutational bias dominated codon usage, patterns of codon usage in HYSV were also under the influence of nature selection. Phylogenetic analysis based on RSCU values across different HYSV strains and 42 other bunyaviruses suggested that codon usage pattern in HYSV was the most similar with that of Uukuniemi virus among these bunyaviruses and that viruses belonged to Phlebovirus showed a diversity of codon usage patterns.