Synonymous codon usage analysis of hand, foot and mouth disease viruses: A comparative study on coxsackievirus A6, A10, A16, and enterovirus 71 from 2008 to 2015

Evolutionary dynamics and adaptive analysis of Seneca Valley virus

Over the past 20 years, the Seneca Valley virus (SVV) has emerged in various countries and regions around the world. Infected pigs display symptoms similar to foot-and-mouth disease and other vesicular diseases, causing severe economic losses to affected countries. In recent years, the number of SVV infections has been increasing in Brazil, China, and the United States. In this study, we comprehensively analyzed SVV genomic sequence data from the perspectives of evolutionary dynamics, phylogeography, and codon usage bias. We aimed to gain further insights into SVV's genetic diversity, spatiotemporal distribution patterns, and evolutionary adaptations. Phylogenetic analysis revealed that SVV has evolved into eight distinct lineages. Based on the results of phylogeographic analysis, it is speculated that the United States might have been the source of SVV, from where it subsequently spread to different countries and regions. Moreover, our analysis of positive selection sites in SVV capsid proteins suggests their potential importance in the process of receptor recognition. Finally, codon preference analysis indicates that natural selection has been a primary evolutionary driver influencing SVV codon usage bias. In conclusion, our in-depth investigation into SVV's origin, dissemination, evolution, and adaptation emphasizes the significance of SVV surveillance and control measures.

Enterovirus A Shows Unique Patterns of Codon Usage Bias in Conventional Versus Unconventional Clade

Enterovirus A (EV-A) species cause hand, foot and mouth disease (HFMD), threatening the health of young children. Understanding the mutual codon usage pattern of the virus and its host(s) has fundamental and applied values. Here, through examining multiple codon usage parameters, we found that the codon usage bias among EV-A strains varies and is clade-specific. EVA76, EVA89, EVA90, EVA91 and EVA92, the unconventional clade of EV-A strains, show unique codon usage pattern relative to the two conventional clades, including EVA71, CVA16, CVA6 and CVA10, etc. Analyses of Effective Number of Codon (ENC), Correspondence Analysis (COA) and Parity Rule 2 (PR2), etc., revealed that the codon usage patterns of EV-A strains are shaped by mutation pressure and natural selection. Based on the neutrality analysis, we determined the dominant role of natural selection in the formation of the codon usage bias of EV-A. In addition, we have determined the codon usage compatibility of potential hosts for EV-A strains using codon adaptation index (CAI), relative codon deoptimization index (RCDI) and similarity index (SiD) analyses, and found that EV-A showed host-specific codon adaptation patterns in different clades. Finally, we confirmed that the unique codon usage pattern of the unconventional clade affected protein expression level in human cell lines. In conclusion, we identified novel characteristics of codon usage bias in distinct EV-A clades associated with their host range, transmission and pathogenicity.

Senecavirus A Enhances Its Adaptive Evolution via Synonymous Codon Bias Evolution

Synonymous codon bias in the viral genome affects protein translation and gene expression, suggesting that the synonymous codon mutant plays an essential role in influencing virulence and evolution. However, how the recessive mutant form contributes to virus evolvability remains elusive. In this paper, we characterize how the Senecavirus A (SVA), a picornavirus, utilizes synonymous codon mutations to influence its evolution, resulting in the adaptive evolution of the virus to adverse environments. The phylogenetic tree and Median-joining (MJ)-Network of these SVA lineages worldwide were constructed to reveal SVA three-stage genetic development clusters. Furthermore, we analyzed the codon bias of the SVA genome of selected strains and found that SVA could increase the GC content of the third base of some amino acid synonymous codons to enhance the viral RNA adaptive evolution. Our results highlight the impact of recessive mutation of virus codon bias on the evolution of the SVA and uncover a previously underappreciated evolutionary strategy for SVA. They also underline the importance of understanding the genetic evolution of SVA and how SVA adapts to the adverse effects of external stress.

Host Adaptation and Evolutionary Analysis of Zaire ebolavirus: Insights From Codon Usage Based Investigations

Ebola virus (EBOV) has caused several outbreaks as the consequence of spillover events from zoonotic sources and has resulted in huge death tolls. In spite of considerable progress, a thorough know-how regarding EBOV adaptation in various host species and detailed information about the potential reservoirs of EBOV still remains obscure. The present study was executed to examine the patterns of codon usage and its associated influence in the adaptation of EBOV to potential hosts that dwell in Africa, the origin of the viral outbreaks. Correspondence analysis (CA) revealed that the codon usage signature in EBOV is a complex interplay of factors including compositional bias and natural selection, with the latter having a more pronounced impact. Low codon usage bias in EBOV indicates a flexibility of the viruses in adapting to diverse range of hosts with different codon usage architectures. EBOV adaptation in potential hosts, as estimated by codon adaptation index (CAI) and relative codon deoptimization index (RCDI), revealed that the viruses were relatively better adapted to African primates than other mammals examined, which might account for the high fatality rate of primates owing to EBOV infection. Bats have been speculated as natural reservoirs of EBOV. In the present analysis it was interesting to note that EBOV displayed lower degrees of adaptation, as estimated by CAI and RCDI, with bats in comparison to the primate hosts. Lower degrees of adaptation might contribute to long-term co-existence and circulation of the viral pathogens in bat populations. Codon usage patterns of EBOV isolates associated with different outbreaks varied significantly, with discrete patterns between the West and Central African isolates. Additional evolutionary analyses indicated that the West African Epidemic began with an initial spillover infection and there was more than one population of EBOV circulating in the natural reservoir in the Democratic Republic of the Congo. The present study yields valuable information regarding the possible circulation of EBOV in various African mammals.

The fit of codon usage of human-isolated avian influenza A viruses to human

Avian influenza A viruses (AIVs) classify into 18 hemagglutinin (HA) and 11 neuraminidase (NA) subtypes. Even though H1N1 and H3N2 subtypes usually circulate among humans leading to infection, occasionally, H5, H6, H7, H9, and H10 that circulate in poultry also infect humans, and especially H5N1 and H7N9. Efficient virus replication is a critical factor that influences infection. Codon usage of a virus must coevolve with its host for efficient viral replication, therefore, we conduct a comprehensive analysis of codon usage bias in human-isolated AIVs to test their adaptation to host expression system. The relative synonymous codon usage (RSCU) pattern, and the codon adaptation index (CAI) are calculated for this purpose. We find that all human-isolated AIVs tend to eliminate GC and CpG compositions, which may prevent activation of the host innate immune system. Although codon usage differs between AIV subtypes, our data support the conclusion that natural selection has played a major role and mutation pressure a minor role in shaping codon usage bias in all AIVs. Our efforts discover that codon usage of genes encoding surface proteins of H5N1, and the polymerase genes of H7N9 has better fit to the human expression system. This may associate with their better replication and infection in human.

Genome Analysis of Coxsackievirus A4 Isolates From Hand, Foot, and Mouth Disease Cases in Shandong, China

Coxsackievirus A4 (CVA4) is one of the most prevalent pathogens associated with hand, foot and mouth disease (HFMD), an acute febrile illness in children, and is also associated with acute localized exanthema, myocarditis, hepatitis and pancreatitis. Despite this, limited CVA4 genome sequences are currently available. Herein, complete genome sequences from CVA4 strains (n = 21), isolated from patients with HFMD in Shandong province, China between 2014 and 2016, were determined and phylogenetically characterized. Phylogenetic analysis of the VP1 gene from a larger CVA4 collection (n = 175) showed that CVA4 has evolved into four separable genotypes: A, B, C, and D; and genotype D could be further classified in to two sub-genotypes: D1 and D2. Each of the 21 newly described genomes derived from isolates that segregated with sub-genotype D2. The CVA4 genomes displayed significant intra-genotypic genetic diversity with frequent synonymous substitutions occurring at the third codon positions, particularly within the P2 region. However, VP1 was relatively stable and therefore represents a potential target for molecular diagnostics assays and also for the rational design of vaccine epitopes. The substitution rate of VP1 was estimated to be 5.12 × 10^-3 substitutions/site/year, indicative of ongoing CVA4 evolution. Mutations at amino acid residue 169 in VP1 gene may be responsible for differing virulence of CVA4 strains. Bayesian skyline plot analysis showed that the population size of CVA4 has experienced several dynamic fluctuations since 1948. In summary, we describe the phylogenetic and molecular characterization of 21 complete genomes from CVA4 isolates which greatly enriches the known genomic diversity of CVA4 and underscores the need for further surveillance of CVA4 in China.

Characterization of an inactivated whole-virus bivalent vaccine that induces balanced protective immunity against coxsackievirus A6 and A10 in mice

Coxsackievirus A6 (CVA6) and CVA10 are two of the major pathogens associated with hand, foot and mouth disease (HFMD) in children. The majority of CVA6 and CVA10 infections result in mild, self-limiting episodes (fever and herpangina) in pediatric populations; however, in some cases, can proceed to severe neurological disease and death. Efforts to mitigate viral transmission to decrease the morbidity and mortality associated with infection would be greatly strengthened by the availability of an efficacious CVA6 and CVA10 bivalent vaccine. Here we report the immunogenicity and protective efficacy of a bivalent combination vaccine comprised of formaldehyde-inactivated, whole-virus CVA6 and CVA10. We demonstrate that subcutaneous delivery of the bivalent vaccine can induce antigen-specific systemic immune responses, particularly the induction of polyfunctional T cells, which elicit active immunization to achieve a protection rate of >80% in the infected neonatal mice. Furthermore, passive transfer of the antisera from vaccinated mice potently protected recipient mice against CVA6 and CVA10 challenge. Importantly, the bivalent vaccine could induce high levels of IgG and neutralizing antibodies in adult female mice and the maternal antibody transmitted to the recipient mice played an important role in controlling homotypic and heterotypic CVA6 and CVA10 infections and viral replication in vivo. Collectively, these findings indicate that there is no immunological interference between the two antigens with respect to their ability to induce virus-specific immune responses, and thus provides proof-of-concept for further development of multivalent vaccines for broad protection against HFMD.

The adaptation of codon usage of +ssRNA viruses to their hosts

Viruses depend on their host's cellular structure to survive. Most of them do not have tRNAs, their translation relies on hosts' tRNA pools. Over the course of evolution, viruses needed to optimally exploit cellular processes of their host. Thus, codon usage of a virus should coevolve with its host to efficiently and rapidly replicate. Some viruses can invade a broad spectrum of hosts (BSTVs), while others can invade a narrow spectrum only (NSTVs). Consequently, we test the hypothesis that similarity of codon usage preference and the degree of matching between BSTVs and their hosts will be lower than that of NSTVs, which only need to coevolve with few hosts. We compare the patterns of codon usage in 255 virus genomes to test this hypothesis. Our results show that NSTVs have a higher degree of matching to their hosts' tRNA pools than BSTVs. Further, analysis of the effective number of codons (ENC) infers that codon usage bias of NSTVs is relatively stronger than that of BSTVs. Thus, codon usage of NSTVs tends to better match their host than that of BSTVs. This supports the hypothesis that viruses adapt to the expression system of their host(s).

Rapid detection of hand, foot and mouth disease enterovirus genotypes by multiplex PCR

Hand, foot and mouth disease (HFMD) is a pediatric disease associated with infection by enterovirus (EV) genotypes. The major HFMD EV pathogens are enterovirus A71 (EVA71) and coxsackievirus A16 (CVA16); however, recently, coxsackievirus A6 (CVA6) and coxsackievirus A10 (CVA10) have also emerged. EV genotypes cannot be distinguished on clinical grounds and a new methodology for the rapid detection of the four major HFMD EV genotypes is urgently required. In the present study, a multiplex real-time PCR assay was established for the simultaneous detection of CVA6, CVA10, CVA16 and EVA71. The specificity and sensitivity of the assay was determined on a validation panel of clinical samples, comprising cerebrospinal fluid (n = 51), blood (n = 39), feces (n = 58) and throat swabs (n = 29). The results showed that the multiplex real-time PCR exhibited high specificity, no cross-reactivity with other EV genotypes, lower limits of detection for CVA6, CVA10, CVA16 and EVA71 were 4 × 10³, 4 × 10², 5 × 10², and 3 × 10³ copies/μL, respectively and had comparable sensitivity to singleplex assays testing clinical samples. The multiplex real-time PCR methodology established in this study can be employed for the rapid detection of the four most prevalent HFMD-associated EVs, for epidemiologic surveillance of circulating EV genotypes and for assessing treatment responses and vaccine studies.