Codon usage bias and recombination events for neuraminidase and hemagglutinin genes in Chinese isolates of influenza A virus subtype H9N2

Dengue virus preferentially uses human and mosquito non-optimal codons

Codon optimality refers to the effect that codon composition has on messenger RNA (mRNA) stability and translation level and implies that synonymous codons are not silent from a regulatory point of view. Here, we investigated the adaptation of virus genomes to the host optimality code using mosquito-borne dengue virus (DENV) as a model. We demonstrated that codon optimality exists in mosquito cells and showed that DENV preferentially uses nonoptimal (destabilizing) codons and avoids codons that are defined as optimal (stabilizing) in either human or mosquito cells. Human genes enriched in the codons preferentially and frequently used by DENV are upregulated during infection, and so is the tRNA decoding the nonoptimal and DENV preferentially used codon for arginine. We found that adaptation during single-host passaging in human or mosquito cells results in the selection of synonymous mutations towards DENV's preferred nonoptimal codons that increase virus fitness. Finally, our analyses revealed that hundreds of viruses preferentially use nonoptimal codons, with those infecting a single host displaying an even stronger bias, suggesting that host-pathogen interaction shapes virus-synonymous codon choice.

Novel Plastid Genome Characteristics in Fugacium kawagutii and the Trend of Accelerated Evolution of Plastid Proteins in Dinoflagellates

Typical (peridinin-containing) dinoflagellates possess plastid genomes composed of small plasmids named "minicircles". Despite the ecological importance of dinoflagellate photosynthesis in corals and marine ecosystems, the structural characteristics, replication dynamics, and evolutionary forcing of dinoflagellate plastid genomes remain poorly understood. Here, we sequenced the plastid genome of the symbiodiniacean species Fugacium kawagutii and conducted comparative analyses. We identified psbT-coding minicircles, features previously not found in Symbiodiniaceae. The copy number of F. kawagutii minicircles showed a strong diel dynamics, changing between 3.89 and 34.3 copies/cell and peaking in mid-light period. We found that F. kawagutii minicircles are the shortest among all dinoflagellates examined to date. Besides, the core regions of the minicircles are highly conserved within genus in Symbiodiniaceae. Furthermore, the codon usage bias of the plastid genomes in Heterocapsaceae, Amphidiniaceae, and Prorocentraceae species are greatly influenced by selection pressure, and in Pyrocystaceae, Symbiodiniaceae, Peridiniaceae, and Ceratiaceae species are influenced by both natural selection pressure and mutation pressure, indicating a family-level distinction in codon usage evolution in dinoflagellates. Phylogenetic analysis using 12 plastid-encoded proteins and five nucleus-encoded plastid proteins revealed accelerated evolution trend of both plastid- and nucleus-encoded plastid proteins in peridinin- and fucoxanthin-dinoflagellate plastids compared to plastid proteins of nondinoflagellate algae. These findings shed new light on the structure and evolution of plastid genomes in dinoflagellates, which will facilitate further studies on the evolutionary forcing and function of the diverse dinoflagellate plastids. The accelerated evolution documented here suggests plastid-encoded sequences are potentially useful for resolving closely related dinoflagellates.

The genetic evolution and codon usage pattern of severe fever with thrombocytopenia syndrome virus

Severe fever with thrombocytopenia syndrome (SFTS) is a newly emerging zoonotic infectious disease caused by the SFTS virus (SFTSV), which has been continuously circulating in Eastern Asia in recent years. Although the evolution of SFTSV has been investigated, the evolutionary changes associated with codon usage have not been reported. Thus, a comprehensive genetic and codon usage bias analysis of SFTSV was conducted to elucidate the genetic diversity and evolutionary relationships in a novel perspective. The study amplified and sequenced fifteen SFTSV strains from a prefecture of Zhejiang Province, Eastern China in 2020, where SFTS cases have been continuously reported in the past decade. Phylogenetic analysis was conducted based on the complete coding sequences of SFTSV segments. It suggested that all SFTSV strains circulating in Zhejiang were clustered with Japanese and Korean strains, which belonged to two different genotypes. Meanwhile, thirty-nine genetic reassortants classified into nineteen different reassortment forms were identified, while 45 recombination events in 41 SFTSV strains were found. Codon usage patterns were further analyzed to understand the evolutionary changes in relation to genotype and host. And it revealed that codon usage bias was mainly driven by natural selection rather than mutation pressure. In addition, the codon adaptation index (CAI) analysis demonstrated the strong adaptability of SFTSV to Gallus gallus and Homo sapiens. Similarity index (SiD) analysis indicated that Haemaphysalis longicornis posed a strong selection pressure to SFTSV. In conclusion, this study revealed that the genetic diversity of SFTSV is gradually increasing. The codon usage analysis suggested that codon usage bias of SFTSV was mainly driven by natural selection, and SFTSV has evolved host-specific codon usage patterns. This contributes to the development of control measures against SFTSV.

Compositional biases in RNA viruses: Causes, consequences and applications

If each of the four nucleotides were represented equally in the genomes of viruses and the hosts they infect, each base would occur at a frequency of 25%. However, this is not observed in nature. Similarly, the order of nucleotides is not random (e.g., in the human genome, guanine follows cytosine at a frequency of ~0.0125, or a quarter the number of times predicted by random representation). Codon usage and codon order are also nonrandom. Furthermore, nucleotide and codon biases vary between species. Such biases have various drivers, including cellular proteins that recognize specific patterns in nucleic acids, that once triggered, induce mutations or invoke intrinsic or innate immune responses. In this review we examine the types of compositional biases identified in viral genomes and current understanding of the evolutionary mechanisms underpinning these trends. Finally, we consider the potential for large scale synonymous recoding strategies to engineer RNA virus vaccines, including those with pandemic potential, such as influenza A virus and Severe Acute Respiratory Syndrome Coronavirus Virus 2. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Evolution and Genomics > Computational Analyses of RNA RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition.

Phylogenetic and codon usage analysis for replicase and capsid genes of porcine circovirus 3

Porcine circovirus type 3 (PCV3) is a highly contagious virus belonging to the family Circoviridae that causes the severe dermatitis and nephropathy syndrome. To date, PCV3 has a worldwide distribution and bring huge economic losses to swine industry. Replicase (Rep) and capsid (Cap) are two major coded proteins of PCV3. Considering the large number of new PCV3 isolates were reported in the past few years and the research for the codon usage pattern of Rep and Cap genes was still a gap, phylogenetic and codon usage analysis of these two genes was performed. Phylogenetic analyses showed that Rep genes in PCV3a were dispersed with no clear clusters while corresponding sequences in PCV3b clustered into two groups and Cap genes clustered into distinct clades according to different genotypes. Relative synonymous codon usage (RSCU) analysis revealed that the codon usage bias existed and effective number of codon (ENC) analysis showed that the bias was slight low. ENC-GC3s plot indicated that mutational pressure and other factors both played a role in PCV3 codon usage and neutrality plot analysis showed that natural selection was the main force influencing the codon usage pattern. The results presented here provided the important basic data on codon usage pattern of Rep and Cap genes, and a better understanding of the evolution and potential origin of PCV3.

Genetic changes and evolutionary analysis of canine circovirus

Canine circovirus (canineCV) has been found to be associated with vasculitis, hemorrhage, hemorrhagic enteritis, and diarrhea of canines. CanineCV, like other circoviruses, may also be associated with lymphoid depletion and immunosuppression. This circovirus has been detected worldwide in different countries and species. Recombination and mutation events in the canineCV genome have been described, indicating that the virus is continuing to evolve. However, the origin, codon usage patterns, and host adaptation of canineCV remain to be studied. Here, the coding sequences of 93 canineCV sequences available in the GenBank database were used for analysis. The results showed that canineCV sequences could be classified into five genotypes, as confirmed by phylogenetic and principal component analysis (PCA). Maximum clade credibility (MCC) and maximum-likelihood (ML) trees suggested that canineCV originated from bat circovirus. G/T and A/C nucleotide biases were observed in ORF1 and ORF2, respectively, and a low codon usage bias (CUB) was found in canineCV using an effective number of codon (ENC) analysis. Correlation analysis, ENC plot analysis and neutrality plot analysis indicated that the codon usage pattern was mainly shaped by natural selection. Codon adaptation index (CAI) analysis, relative codon deoptimization index (RCDI) analysis, and similarity index (SiD) analysis revealed a better adaption to Vulpes vulpes than to Canis familiaris. Furthermore, a cross-species transmission hypothesis that canineCV may have evolved from bats (origin analysis) and subsequently adapted to wolves, arctic foxes, dogs, and red foxes, was proposed. This study contributes to our understanding of the factors related to canineCV evolution and host adaption.

Evolution and genetic diversity of SARS-CoV-2 in Africa using whole genome sequences

BACKGROUND

The ongoing SARS-CoV-2 pandemic was introduced into Africa on 14th February 2020 and has rapidly spread across the continent causing a severe public health crisis and mortality. We investigated the genetic diversity and evolution of this virus during the early outbreak months, between 14th February to 24th April 2020, using whole genome sequences.

METHODS

We performed recombination analysis against closely related CoV strains, Bayesian time scaled phylogeny, and investigation of spike protein amino acid mutations.

RESULTS

Recombination signals were observed between the Afr-SARS-CoV-2 sequences and reference sequences within the RdRPs and S genes. The evolutionary rate of the Afr-SARS-CoV-2 was 4.133 × 10-4 Highest Posterior Density (HPD 4.132 × 10-4 to 4.134 × 10-4) substitutions/site/year. The time to most recent common ancestor (TMRCA) of the African strains was December 7th 2019, (95% HPD November 12th 2019-December 29th 2019). The Afr-SARCoV-2 sequences diversified into two lineages A and B, with B being more diverse with multiple sub-lineages confirmed by both maximum clade credibility (MCC) tree and PANGOLIN software. There was a high prevalence of the D614G spike protein amino acid mutation 59/69 (82.61%) among the African strains.

CONCLUSION

This study has revealed a rapidly diversifying viral population with the G614G spike protein variant dominatinge advocate for up scaling NGS sequencing platforms across Africa to enhance surveillance and aid control effort of SARS-CoV-2 in Africa.

Comprehensive analysis of synonymous codon usage patterns and influencing factors of porcine epidemic diarrhea virus

Porcine epidemic diarrhea virus (PEDV) is an enteric pathogen belonging to the family Coronaviridae that causes the porcine epidemic diarrhea, a highly contagious disease with high mortality in piglets and symptoms that include dehydration and severe diarrhea. Considering the high frequency of genetic mutations in PEDV and its potential for interspecies transmission, as it can infect and replicate in bat and human cells, a comprehensive analysis of its codon usage bias was performed. The effective number of codons (ENC) and the relative synonymous codon usage (RSCU) were determined, revealing codon usage bias in the PEDV genome. Principal component analysis (PCA), an ENC plot, and a parity rule 2 (PR2) plot showed that mutation pressure and natural selection have influenced the codon usage bias of the PEDV genomes. Correlation analysis with GRAVY and aromaticity values and neutrality plot analysis indicated that natural selection was the main force influencing the codon usage pattern, while mutation pressure played a minor role. This study provides valuable basic data for further fundamental research on evolution of PEDV.

Genetic and evolutionary analysis of emerging H3N2 canine influenza virus

The H3N2 canine influenza virus (CIV) originated from an avian species. Since its emergence, it has circulated in multiple states and has caused pandemics among dog populations; however, no comprehensive studies have explored the causes driving these ongoing cases. The study of the codon usage patterns of viruses can reveal the genetic changes required for the viruses to adapt to new hosts and the external environment. Here we performed a thorough genetic, evolutionary, and codon usage analysis. We identified three evolutionary H3N2 CIV clades from a timescaled phylogenetic tree, namely, Origin, China, and Korea/USA, by principal component analysis (PCA). Additionally, we found a low codon usage bias and that mutation pressure, natural selection, and dinucleotide abundance shape the codon usage bias of H3N2 CIVs, with natural selection being more crucial than the others. Moreover, the human codon adaptation index was similar to that of dogs (the natural host) and cats. In addition, the H3N2 CIV similarity index values were higher than those of the avian influenza virus (AIV), suggesting viral adaptation to the host. Therefore, H3N2 CIVs may pose a potential risk to public health in the future, and further epidemiologic, evolutionary, and pathogenetic studies are required.

Comparative analysis of codon usage pattern and its influencing factors in Schistosoma japonicum and Ascaris suum

Schistosoma japonicum and Ascaris suum are considered as the major parasites of human which cause various life threatening diseases such as schistomiasis and ascariasis. The codon usage bias (CUB) is known as the phenomenon of more usage of a specific codon than the other synonymous codons for an amino acid. The factors that influence the codon usage bias are mutation pressure, natural selection, gene expression, gene length, GC content, RNA stability, recombination rates, codon position etc. Here we had used various bioinformatic tools and statistical analyses to understand the compositional features, expression level and codon usage bias in the genes of these two species.After estimating the effective number of codon (ENC) in both the species, codon usage bias was found to be low and gene expression was high. The nucleobase A and T were used most often than C and G. From neutrality plot and correspondence analysis it was found that both natural selection and mutation pressure played an important role in shaping the codon usage pattern of both species. Moreover, natural selection played a major role while mutation pressure played a minor role in shaping the codon usage bias in S. japonicum and A.suum. This is the first report on the codon usage biology in S. japonicum and A.suum, and the factors influencing their codon usage bias. These results are expected to be useful for genetic engineering and evolutionary studies.

Comprehensive analysis of the codon usage patterns in the envelope glycoprotein E2 gene of the classical swine fever virus

The classical swine fever virus (CSFV), circulating worldwide, is a highly contagious virus. Since the emergence of CSFV, it has caused great economic loss in swine industry. The envelope glycoprotein E2 gene of the CSFV is an immunoprotective antigen that induces the immune system to produce neutralizing antibodies. Therefore, it is essential to study the codon usage of the E2 gene of the CSFV. In this study, 140 coding sequences of the E2 gene were analyzed. The value of effective number of codons (ENC) showed low codon usage bias in the E2 gene. Our study showed that codon usage could be described mainly by mutation pressure ENC plot analysis combined with principal component analysis (PCA) and translational selection-correlation analysis between the general average hydropathicity (Gravy) and aromaticity (Aroma), and nucleotides at the third position of codons (A3s, T3s, G3s, C3s and GC3s). Furthermore, the neutrality analysis, which explained the relationship between GC12s and GC3s, revealed that natural selection had a key role compared with mutational bias during the evolution of the E2 gene. These results lay a foundation for further research on the molecular evolution of CSFV.

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.

The ecology and adaptive evolution of influenza A interspecies transmission

Since 2013, there have been several alarming influenza-related events; the spread of highly pathogenic avian influenza H5 viruses into North America, the detection of H10N8 and H5N6 zoonotic infections, the ongoing H7N9 infections in China and the continued zoonosis of H5N1 viruses in parts of Asia and the Middle East. The risk of a new influenza pandemic increases with the repeated interspecies transmission events that facilitate reassortment between animal influenza strains; thus, it is of utmost importance to understand the factors involved that promote or become a barrier to cross-species transmission of Influenza A viruses (IAVs). Here, we provide an overview of the ecology and evolutionary adaptations of IAVs, with a focus on a review of the molecular factors that enable interspecies transmission of the various virus gene segments.

Evolution of codon usage in Zika virus genomes is host and vector specific

The codon usage patterns of viruses reflect the evolutionary changes that allow them to optimize their survival and adapt their fitness to the external environment and, most importantly, their hosts. Here we report the genotype-specific codon usage patterns of Zika virus (ZIKV) strains from the current and previous outbreaks. Several genotype-specific and common codon usage traits were noted in the ZIKV coding sequences, indicating their independent evolutionary origins from a common ancestor. The overall influence of natural selection was more profound than that of mutation pressure, acting on a specific set of viral genes in the Asian-genotype ZIKV strains from the recent outbreak. An interplay between codon adaptation and deoptimization may have allowed the virus to adapt to multiple host and vectors and is reported for the first time in ZIKV genomes. Combining our codon analysis with geographical data on Aedes populations in the Americas suggested that ZIKV has evolved host- and vector-specific codon usage patterns to maintain successful replication and transmission chains within multiple hosts and vectors.

Molecular and Phenotypic Characterization of a Highly Evolved Type 2 Vaccine-Derived Poliovirus Isolated from Seawater in Brazil, 2014

A type 2 vaccine-derived poliovirus (VDPV), differing from the Sabin 2 strain at 8.6% (78/903) of VP1 nucleotide positions, was isolated from seawater collected from a seaport in São Paulo State, Brazil. The P1/capsid region is related to the Sabin 2 strain, but sequences within the 5'-untranslated region and downstream of the P1 region were derived from recombination with other members of Human Enterovirus Species C (HEV-C). The two known attenuating mutations had reverted to wild-type (A481G in the 5'-UTR and Ile143Thr in VP1). The VDPV isolate had lost the temperature sensitive phenotype and had accumulated amino acid substitutions in neutralizing antigenic (NAg) sites 3a and 3b. The date of the initiating OPV dose, estimated from the number of synonymous substitutions in the capsid region, was approximately 8.5 years before seawater sampling, a finding consistent with a long time of virus replication and possible transmission among several individuals. Although no closely related type 2 VDPVs were detected in Brazil or elsewhere, this VDPV was found in an area with a mobile population, where conditions may favor both viral infection and spread. Environmental surveillance serves as an important tool for sensitive and early detection of circulating poliovirus in the final stages of global polio eradication.

Genomic analysis of codon usage shows influence of mutation pressure, natural selection, and host features on Marburg virus evolution

BACKGROUND

The Marburg virus (MARV) has a negative-sense single-stranded RNA genome, belongs to the family Filoviridae, and is responsible for several outbreaks of highly fatal hemorrhagic fever. Codon usage patterns of viruses reflect a series of evolutionary changes that enable viruses to shape their survival rates and fitness toward the external environment and, most importantly, their hosts. To understand the evolution of MARV at the codon level, we report a comprehensive analysis of synonymous codon usage patterns in MARV genomes. Multiple codon analysis approaches and statistical methods were performed to determine overall codon usage patterns, biases in codon usage, and influence of various factors, including mutation pressure, natural selection, and its two hosts, Homo sapiens and Rousettus aegyptiacus.

RESULTS

Nucleotide composition and relative synonymous codon usage (RSCU) analysis revealed that MARV shows mutation bias and prefers U- and A-ended codons to code amino acids. Effective number of codons analysis indicated that overall codon usage among MARV genomes is slightly biased. The Parity Rule 2 plot analysis showed that GC and AU nucleotides were not used proportionally which accounts for the presence of natural selection. Codon usage patterns of MARV were also found to be influenced by its hosts. This indicates that MARV have evolved codon usage patterns that are specific to both of its hosts. Moreover, selection pressure from R. aegyptiacus on the MARV RSCU patterns was found to be dominant compared with that from H. sapiens. Overall, mutation pressure was found to be the most important and dominant force that shapes codon usage patterns in MARV.

CONCLUSIONS

To our knowledge, this is the first detailed codon usage analysis of MARV and extends our understanding of the mechanisms that contribute to codon usage and evolution of MARV.

Published sequences do not support transfer of oseltamivir resistance mutations from avian to human influenza A virus strains

Background

Tamiflu (oseltamivir phosphate ester, OE) is a widely used antiviral active against influenza A virus. Its active metabolite, oseltamivir carboxylate (OC), is chemically stable and secreted into wastewater treatment plants. OC contamination of natural habitats of waterfowl might induce OC resistance in influenza viruses persistently infecting waterfowl, and lead to transfer of OC-resistance from avian to human influenza. The aim of this study was to evaluate whether such has occurred.

Methods

A genomics approach including phylogenetic analysis and probability calculations for homologous recombination was applied on altogether 19,755 neuraminidase (N1 and N2) genes from virus sampled in humans and birds, with and without resistance mutations.

Results

No evidence for transfer of OE resistance mutations from avian to human N genes was obtained, and events suggesting recombination between human and avian influenza virus variants could not be traced in the sequence material studied.

Conclusions

The results indicate that resistance in influenza viruses infecting humans is due to the selection pressure posed by the global OE administration in humans rather than transfer from avian influenza A virus strains carrying mutations induced by environmental exposure to OC.

Electronic supplementary material

The online version of this article (doi:10.1186/s12879-015-0860-9) contains supplementary material, which is available to authorized users.

Species based synonymous codon usage in fusion protein gene of Newcastle disease virus

Newcastle disease is highly pathogenic to poultry and many other avian species. However, the Newcastle disease virus (NDV) has also been reported from many non-avian species. The NDV fusion protein (F) is a major determinant of its pathogenicity and virulence. The functionalities of F gene have been explored for the development of vaccine and diagnostics against NDV. Although the F protein is well studied but the codon usage and its nucleotide composition from NDV isolated from different species have not yet been explored. In present study, we have analyzed the factors responsible for the determination of codon usage in NDV isolated from four major avian host species. The F gene of NDV is analyzed for its base composition and its correlation with the bias in codon usage. Our result showed that random mutational pressure is responsible for codon usage bias in F protein of NDV isolates. Aromaticity, GC3s, and aliphatic index were not found responsible for species based synonymous codon usage bias in F gene of NDV. Moreover, the low amount of codon usage bias and expression level was further confirmed by a low CAI value. The phylogenetic analysis of isolates was found in corroboration with the relatedness of species based on codon usage bias. The relationship between the host species and the NDV isolates from the host does not represent a significant correlation in our study. The present study provides a basic understanding of the mechanism involved in codon usage among species.

Characterization of the porcine epidemic diarrhea virus codon usage bias

•

Codon usage bias of porcine epidemic diarrhea virus is low.

•

Mutational bias and natural selection pressure influence codon usage bias of porcine epidemic diarrhea virus.

•

Natural selection plays an increasingly significant role during evolution of porcine epidemic diarrhea virus.

Porcine epidemic diarrhea virus (PEDV) has been responsible for several recent outbreaks of porcine epidemic diarrhea (PED) and has caused great economic loss in the swine-raising industry. Considering the significance of PEDV, a systemic analysis was performed to study its codon usage patterns. The relative synonymous codon usage value of each codon revealed that codon usage bias exists and that PEDV tends to use codons that end in T. The mean ENC value of 47.91 indicates that the codon usage bias is low. However, we still wanted to identify the cause of this codon usage bias. A correlation analysis between the codon compositions (A3s, T3s, G3s, C3s, and GC3s), the ENC values, and the nucleotide contents (A%, T%, G%, C%, and GC%) indicated that mutational bias plays role in shaping the PEDV codon usage bias. This was further confirmed by a principal component analysis between the codon compositions and the axis values. Using the Gravy, Aroma, and CAI values, a role of natural selection in the PEDV codon usage pattern was also identified. Neutral analysis indicated that natural selection pressure plays a more important role than mutational bias in codon usage bias. Natural selection also plays an increasingly significant role during PEDV evolution. Additionally, gene function and geographic distribution also influence the codon usage bias to a degree.

Analysis of synonymous codon usage patterns in duck hepatitis A virus: a comparison on the roles of mutual pressure and natural selection

Codon usage patterns of duck hepatitis A virus (HAV) were studied in the present study. The major trends of codon usage patterns were analyzed using principal component analysis on the basis of the relative synonymous codon usage values. Correlation analysis was utilized to reveal the associations of the first two major axes of PCA and nucleotide- or amino acid-relevant indices. Our results showed that compositional constraint and/or mutational pressure are major factors influencing codon usage bias patterns of HAV. However, the influence of natural selection is also prevalent, as indicated by strongly significant correlations between the hydrophobicity, aromaticity, aliphaticity and ionization and the first axis of PCA. Also, ionization could characterize the second axis of PCA. At last, maximum likelihood phylogenetic analysis shows that there are no remarkable geographic clustering patterns of HAV strains in the phylogenetic tree. However, through MANOVA test, there are significant differences on the codon usage patterns among HAV strains from different countries. In conclusion, both mutational pressure and natural selection are of equally great importance to codon usage patterns of duck HAV genomes.

Extensive homologous recombination in classical swine fever virus: A re-evaluation of homologous recombination events in the strain AF407339

In this short report, the genome-wide homologous recombination events were re-evaluated for classical swine fever virus (CSFV) strain AF407339. We challenged a previous study which suggested only one recombination event in AF407339 based on 25 CSFV genomes. Through our re-analysis on the 25 genomes in the previous study and the 41 genomes used in the present study, we argued that there should be possibly at least two clear recombination events happening in AF407339 through genome-wide scanning. The reasons for identifying only one recombination event in the previous study might be due to the limited number of available CSFV genome sequences at that time and the limited usage of detection methods. In contrast, as identified by most detection methods using all available CSFV genome sequences, two major recombination events were found at the starting and ending zones of the genome AF407339, respectively. The first one has two parents AF333000 (minor) and AY554397 (major) with beginning and ending breakpoints located at 19 and 607 nt of the genome respectively. The second one has two parents AF531433 (minor) and GQ902941 (major) with beginning and ending breakpoints at 8397 and 11,078 nt of the genome respectively. Phylogenetic incongruence analysis using neighbor-joining algorithm with 1000 bootstrapping replicates further supported the existence of these two recombination events. In addition, we also identified additional 18 recombination events on the available CSFV strains. Some of them may be trivial and can be ignored. In conclusion, CSFV might have relatively high frequency of homologous recombination events. Genome-wide scanning of identifying recombination events should utilize multiple detection methods so as to reduce the risk of misidentification.

Natural Selection Determines Synonymous Codon Usage Patterns of Neuraminidase (NA) Gene of the Different Subtypes of Influenza A Virus in Canada

Synonymous codon usage patterns of neuraminidase (NA) gene of 64 subtypes (one is a mixed subtype) of influenza A virus found in Canada were analyzed. In total, 1422 NA sequences were analyzed. Among the subtypes, H1N1 is the prevailing one with 516 NCBI accession records, followed by H3N2, H3N8, and H4N6. The year of 2009 has the highest report records for the NA sequences in Canada, corresponding to the 2009 pandemic event. Correspondence analysis on the RSCU values of the four major subtypes showed that they had distinct clustering patterns in the two-dimensional scatter plot, indicating that different subtypes of IAV utilized different preferential codons. This subtype clustering pattern implied the important influence of natural selection, which could be further evidenced by an extremely flattened regression line in the neutrality plot (GC12 versus G3s plot) and a significant phylogenetic signal on the distribution of different subtypes in the clades of the phylogenetic tree (λ statistic). In conclusion, different subtypes of IAV showed an evolutionary differentiation on choosing different optimal codons. Natural selection played a deterministic role to structure IAV codon usage patterns in Canada.

Evidence of selection pressures of neuraminidase gene (NA) of influenza A virus subtype H5N1 on different hosts in Guangxi Province of China

In the present study, the possible evidence of positive selection was analyzed for the neuraminidase (NA) sequences of Guangxi H5N1 strains of China. Based on an overall site-specific positive selection analysis, it was found that NA gene of H5N1 Guangxi strains underwent purifying selection and no significant positively selected sites were identified. For the branch-specific positive selection analysis, there was no positive selection evidence for the branches leading to different poultry hosts (chicken, duck and goose). Conclusively, positive selection seems not possible (if not rare) for the NA gene in influenza H5N1 subtype, at least for the samples found in Guangxi Province of China.

A comparison of synonymous codon usage bias patterns in DNA and RNA virus genomes: quantifying the relative importance of mutational pressure and natural selection

Codon usage bias patterns have been broadly explored for many viruses. However, the relative importance of mutation pressure and natural selection is still under debate. In the present study, I tried to resolve controversial issues on determining the principal factors of codon usage patterns for DNA and RNA viruses, respectively, by examining over 38000 ORFs. By utilizing variation partitioning technique, the results showed that 27% and 21% of total variation could be attributed to mutational pressure, while 5% and 6% of total variation could be explained by natural selection for DNA and RNA viruses, respectively, in codon usage patterns. Furthermore, the combined effect of mutational pressure and natural selection on influencing codon usage patterns of viruses is substantial (explaining 10% and 8% of total variation of codon usage patterns). With respect to GC variation, GC content is always negatively and significantly correlated with aromaticity. Interestingly, the signs for the significant correlations between GC, gene lengths, and hydrophobicity are completely opposite between DNA and RNA viruses, being positive for DNA viruses while being negative for RNA viruses. At last, GC12 versus G3s plot suggests that natural selection is more important than mutational pressure on influencing the GC content in the first and second codon positions.

Genetic analysis of the PB1-F2 gene of equine influenza virus

Amino acid and phylogenetic analysis showed that PB1-F2 of H7N7 and H3N8 equine influenza virus (EIV) evolved into two and nine clades, respectively. The PB1-F2 gene of H7N7 EIV strains isolated after 1973 was identified as being replaced by that of H3N8 EIV circulating simultaneously. All H7N7 EIV strains before 1973 contained PB1-F2 of 34 aa, while 90 aa after 1973. Meanwhile, most H3N8 EIV strains contained PB1-F2 of 90 or 81 aa. Evolutionary rate of PB1-F2 of H3N8 EIV was similar with that of PB1 in previous study. Meanwhile, PB1-F2 of H7N7 EIV evolved in significantly higher rate when compared with PB1-F2 of H3N8 EIV. Codon usage analysis revealed that PB1-F2 gene of EIV was less biased, which was identified as being determined by three main factors: mutational bias, selection pressure, and gene length. Our studies first in details report the genetic evolution, evolutionary rate, and the factors influencing codon usage bias of PB1-F2 of EIV.

Connecting the study of wild influenza with the potential for pandemic disease

Continuing outbreaks of pathogenic (H5N1) and pandemic (SOIVH1N1) influenza have underscored the need to understand the origin, characteristics, and evolution of novel influenza A virus (IAV) variants that pose a threat to human health. In the last 4-5years, focus has been placed on the organization of large-scale surveillance programs to examine the phylogenetics of avian influenza virus (AIV) and host-virus relationships in domestic and wild animals. Here we review the current gaps in wild animal and environmental surveillance and the current understanding of genetic signatures in potentially pandemic strains.

Base composition and translational selection are insufficient to explain codon usage bias in plant viruses

Viral codon usage bias may be the product of a number of synergistic or antagonistic factors, including genomic nucleotide composition, translational selection, genomic architecture, and mutational or repair biases. Most studies of viral codon bias evaluate only the relative importance of genomic base composition and translational selection, ignoring other possible factors. We analyzed the codon preferences of ssRNA (luteoviruses and potyviruses) and ssDNA (geminiviruses) plant viruses that infect translationally distinct monocot and dicot hosts. We found that neither genomic base composition nor translational selection satisfactorily explains their codon usage biases. Furthermore, we observed a strong relationship between the codon preferences of viruses in the same family or genus, regardless of host or genomic nucleotide content. Our results suggest that analyzing codon bias as either due to base composition or translational selection is a false dichotomy that obscures the role of other factors. Constraints such as genomic architecture and secondary structure can and do influence codon usage in plant viruses, and likely in viruses of other hosts.

Genome dynamics in three different geographical isolates of white spot syndrome virus (WSSV)

White spot syndrome virus (WSSV), the sole member of the monotypic family Nimaviridae, is considered an extremely lethal shrimp pathogen. Despite its impact, some essential biological characteristics related to WSSV genome dynamics, such as the synonymous codon usage pattern and selection pressure in genes, remain to be elucidated. The results show that compositional limitations and mutational pressure determine the codon usage bias and base composition in WSSV. Furthermore, different forces of selective pressure are acting across various regions of the WSSV genome. Finally, this study points out the possible occurrence of two major recombination events.

Homologous recombination in negative sense RNA viruses

Recombination is an important process that influences biological evolution at many different levels. More and more homologous recombination events have been reported among negative sense RNA viruses recently. While sporadic authentic examples indicate that homologous recombination does occur, recombination seems to be generally rare or even absent in most negative sense RNA viruses, and most of the homologous recombination events reported in the literature were likely generated artificially due to lab contamination or inappropriate bioinformatics methods. Homologous recombination in negative sense RNA viruses should be reported with caution in the future, and only after stringent quality control efforts. Moreover, co-infection experiments should be performed to confirm whether recombination can occur.

Analysis of codon usage in type 1 and the new genotypes of duck hepatitis virus

In this study, an abundant (A + U)% and low codon bias were revealed in duck hepatitis virus type 1 (DHV-1) and the new serotype strains isolated from Taiwan, South Korea and Mainland China (DHV-N). The general correlation between base composition and codon usage bias suggests that mutational pressure rather than natural selection is the main factor that determines the codon usage bias in these samples. By comparative analysis of the codon usage patterns of 40 ORFs of DHV, we found that all of DHV-1 strains grouped in genotype C; the DHV-N strains isolated in South Korea and China clustered into genotypes B; and the DHV-N strains isolated from Taiwan clustered into genotypes A. The findings revealed that more than one subtype of DHV-1 circulated in East Asia. Furthermore, the results of phylogenetic analyses based on RSCU values and Clustal W method indicated obvious phylogenetic congruities. This suggested that better genome consistency of DHV may exist in nature and phylogenetic analyses based on RSCU values maybe a good method in classifying genotypes of the virus. Our work might give some clues to the features and some evolutionary information of DHV.

Analysis of codon usage and nucleotide composition bias in polioviruses

Background

Poliovirus, the causative agent of poliomyelitis, is a human enterovirus and a member of the family of Picornaviridae and among the most rapidly evolving viruses known. Analysis of codon usage can reveal much about the molecular evolution of the viruses. However, little information about synonymous codon usage pattern of polioviruses genome has been acquired to date.

Methods

The relative synonymous codon usage (RSCU) values, effective number of codon (ENC) values, nucleotide contents and dinucleotides were investigated and a comparative analysis of codon usage pattern for open reading frames (ORFs) among 48 polioviruses isolates including 31 of genotype 1, 13 of genotype 2 and 4 of genotype 3.

Results

The result shows that the overall extent of codon usage bias in poliovirus samples is low (mean ENC = 53.754 > 40). The general correlation between base composition and codon usage bias suggests that mutational pressure rather than natural selection is the main factor that determines the codon usage bias in those polioviruses. Depending on the RSCU data, it was found that there was a significant variation in bias of codon usage among three genotypes. Geographic factor also has some effect on the codon usage pattern (exists in the genotype-1 of polioviruses). No significant effect in gene length or vaccine derived polioviruses (DVPVs), wild viruses and live attenuated virus was observed on the variations of synonymous codon usage in the virus genes. The relative abundance of dinucleotide (CpG) in the ORFs of polioviruses are far below expected values especially in DVPVs and attenuated virus of polioviruses genotype 1.

Conclusion

The information from this study may not only have theoretical value in understanding poliovirus evolution, especially for DVPVs genotype 1, but also have potential value for the development of poliovirus vaccines.

Analysis of codon usage in Newcastle disease virus

In this study, the relative synonymous codon usage (RSCU) values, effective number of codon (ENC) values, nucleotide contents, and dinucleotide were used to investigate codon usage pattern of each protein-coding gene and genome among 31 Newcastle disease virus (NDV) isolates. The result shows that the overall extent of codon usage bias in NDV is low (mean ENC = 56.15 > 40). The good correlation between the (C + G)₁₂% and (G + C)₃% suggests that the mutational pressure, rather than natural selection, is the main factor that determines the codon usage bias and base component in NDV. It is observed that synonymous codon usage pattern in NDV genes is gene function and geography specific, but not host specific. By contrasting synonymous codon usage patterns of different NDV isolates, we suggest that more than one genotype of NDV circulates in waterfowl in USA; and gene length has no significant effect on the variations of synonymous codon usage in these virus genes. CpG under-represented is a characteristic for NDV to fit in its host. These results not only provide an insight into the variation of codon usage pattern among the genomes of NDV, but also may help in understanding the processes governing the evolution of NDV.

No observed effect of homologous recombination on influenza C virus evolution

The occurrence of homologous recombination in influenza viruses has been under some debate recently. To determine the extent of homologous recombination in influenza C virus, recombination analyses of all available gene sequences of influenza C virus were carried out. No recombination signal was found. With the previous evidence in influenza A and B viruses, it seems that homologous recombination has minimal or no effect on influenza virus evolution.