Analysis of synonymous codon usage in enterovirus 71

Analysis of codon usage bias of WRKY transcription factors in Helianthus annuus

Background

The phenomenon of codon usage bias is known to exist in many genomes and is mainly determined by mutation and selection. Codon usage bias analysis is a suitable strategy for identifying the principal evolutionary driving forces in different organisms. Sunflower (Helianthus annuus L.) is an annual crop that is cultivated worldwide as ornamentals, food plants and for their valuable oil. The WRKY family genes in plants play a central role in diverse regulation and multiple stress responses. Evolutionary analysis of WRKY family genes of H. annuus can provide rich genetic information for developing hybridization resources of the genus Helianthus.

Results

Bases composition analysis showed the average GC content of WRKY genes of H. annuus was 43.42%, and the average GC3 content was 39.60%, suggesting that WRKY gene family prefers A/T(U) ending codons. There were 29 codons with relative synonymous codon usage (RSCU) greater than 1 and 22 codons ending with A and U base. The effective number of codons (ENC) and codon adaptation index (CAI) in WRKY genes ranged from 43.47–61.00 and 0.14–0.26, suggesting that the codon bias was weak and WRKY genes expression level was low. Neutrality analysis found a significant correlation between GC12 and GC3. ENC-plot showed most genes on or close to the expected curve, suggesting that mutational bias played a major role in shaping codon usage. The Parity Rule 2 plot (PR2) analysis showed that the usage of AT and GC was disproportionate. A total of three codons were identified as the optimal codons.

Conclusion

Apart from natural selection effects, most of the genetic evolution in the H. annuus WRKY genome might be driven by mutation pressure. Our results provide a theoretical foundation for elaborating the genetic architecture and mechanisms of H. annuus and contributing to enrich H. annuus genetic resources.

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.

Comprehensive analysis of synonymous codon usage patterns in orf3 gene of porcine epidemic diarrhea virus in China

The ORF3 protein of porcine epidemic diarrhea virus (PEDV) is found to function as an ion channel which influences virus virulence and production. Taking consideration of the importance of PEDV orf3 gene, we have performed comprehensive analysis to investigate its synonymous codon usage patterns. In this study, the results of base composition analysis showed A/T rich and G/C poor in PEDV orf3 genes, and the most abundant base was nucleotide T. The relative synonymous codon usage value in each codon revealed that codon usage bias existed. The mean ENC value of each gene was 48.75, indicating a low codon usage bias, as well as a relatively instable change in PEDV orf3 genes. The general correlation analysis between base composition and codon usage bias indicated that mutational bias has an impact on the PEDV codon usage bias. Neutral analysis suggested that natural selection pressure takes a more important influence than mutational bias in shaping codon usage bias. Moreover, other factors including hydrophobicity and aromaticity have been also found to influence the codon usage variation among the PEDV orf3 genes. This study not only represents the most systematic analysis of codon usage patterns in PEDV orf3 genes, but also provides a basic shaping mechanism of the codon usage bias.

•

ORF3 genes are A/T rich and C/G poor, and codon usage bias existed.

•

Natural selection pressure plays a more important role than mutational bias in codon usage bias.

•

Hydrophobicity and aromaticity also influence the codon usage variation.

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.

Enterovirus 71 inhibits cytoplasmic stress granule formation during the late stage of infection

Stress granules (SGs) are host translationally silent ribonucleo-proteins formed in cells in response to multiple types of environmental stress, including viral infection. We previously showed that the nuclear protein, 68-kDa Src-associated in mitosis protein (Sam68), is recruited to cytoplasm and form the Sam68-positive SGs at 6 hpi, but the Sam68-positive SGs disassembled beyond 12 hpi, suggesting that the SGs might be inhibited during the late stage of Enterovirus 71 (EV71) infection. However, the mechanism and function of this process remains poorly understood. Thus in this study, we demonstrated that EV71 initially induced SGs formation at the early stage of EV71 infection, and confirmed that 2Apro of EV71 was the key viral component that triggered SG formation. In contrast, SGs were diminished as EV71 infection proceeding. At the same time, arsenite-induced SGs were also blocked at the late stage of EV71 infection. This disruption of SGs was caused by viral protease 3Cpro-mediated G3BP1 cleavage. Furthermore, we demonstrated that over-expression of G3BP1-SGs negatively impacted viral replication at the cytopathic effect (CPE), protein, RNA, and viral titer levels. Our novel finding may not only help us to better understand the mechanism how EV71 interacts with the SG response, but also provide mechanistic linkage between cellular stress responses and innate immune activation during EV71 infection.

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

Enterovirus 71 (EV71) and coxsackievirus A16 (CVA16) have been considered major pathogens of hand, foot and mouth disease (HFMD) throughout the world for decades. In recent years, coxsackievirus A6 (CVA6) and coxsackievirus A10 (CVA10) have raised attention as two other serious pathogens of HFMD. The present study focused on the synonymous codon usage of four viruses isolated from 2008 to 2015, with particular attention on P1 (encoding capsid proteins) and P2-P3 regions (both encoding non-structural proteins) in the genomic RNA. Relative synonymous codon usage, effective number of codons, neutrality and correspondence were analyzed. The results indicated that these viruses prefer A/T at the third position in codons rather than G/C. The most frequent codons of 4 essential and 2 semi-essential amino acids, as well as a key amino acid of metabolic junctions (Glu) used in the four viruses are also the most frequently used in humans. Effective number of codons (ENC) values indicated weak codon usage bias in all the viruses. Relatively, the force of mutation pressure in the P1 region was found to be stronger than that in the P2-P3 region, and this force in the P1 region of CVA6 and EV71 was stronger than that of CVA10 and A16. The neutrality analysis results implied that mutation pressure plays a minor role in shaping codon bias of these viruses. Correspondence analysis indicated that the codon usage of EV71 strains varied much more than that of other viruses. In conclusion, the present study provides novel and comparative insight into the evolution of HFMD pathogens at the codon level.

Comparative characterization analysis of synonymous codon usage bias in classical swine fever virus

Classical swine fever virus (CSFV) is responsible for the highly contagious viral disease of swine, and causes great economic loss in the swine-raising industry. Considering the significance of CSFV, a systemic analysis was performed to study its codon usage patterns. In this study, using the complete genome sequences of 76 CSFV representing three genotypes, we firstly analyzed the relative nucleotide composition, effective number of codon (ENC) and synonymous codon usage in CSFV genomes. The results showed that CSFV is GC-moderate genome and the third-ended codons are not preferentially used. Every ENC values in CSFV genomes are >50, indicating that the codon usage bias is comparatively slight. Subsequently, we performed the correspondence analysis (COA) to investigate synonymous codon usage variation among all of the CSFV genomes. We found that codon usage bias in these CSFV genomes is greatly influenced by G + C mutation, which suggests that mutational pressure may be the main factor determining the codon usage biases. Moreover, most of the codon usage bias among different CSFV ORFs is directly related to the nucleotide composition. Other factors, such as hydrophobicity and aromaticity, also influence the codon usage variation among CSFV genomes. Our study represents the most comprehensive analysis of codon usage patterns in CSFV genome and provides a basic understanding of the mechanisms for its codon usage bias.

Enterovirus 71 inhibits cellular type I interferon signaling by inhibiting host RIG-I ubiquitination

Enterovirus 71 (EV71) is a human pathogen that induces hand, foot, and mouth disease (HFMD) and fatal neurological diseases in young children and infants. Pathogenicity of EV71 is likely related to its ability to evade host innate immunity through inhibiting cellular type I interferon signaling. However, it is less well understood the molecular events governing this process. In this study, we found that EV71 infection suppressed the induction of antiviral immunity by inhibiting the expression levels of IFN-β and IFN-stimulated genes (ISGs), such as ISG54 and ISG56, at the late stage of viral infection. At the same time, our results showed that EV71 infection significantly inhibited ubiquitination of RIG-I. In contrast, up-regulation of RIG-I ubiquitination promoted expression of IFN-β and ISGs, suggesting that inhibition of cellular type I interferon signaling was caused by down-regulation of RIG-I ubiquitination during EV71 infection. These results suggest that inhibition of RIG-I-mediated type I IFN responses by EV71 may contribute to the pathogenesis of viral infection.

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

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

The nuclear protein Sam68 is recruited to the cytoplasmic stress granules during enterovirus 71 infection

Our previous study found that the nuclear protein, 68-kDa Src-associated in mitosis protein (Sam68), is translocated to the cytoplasm and forms punctate pattern during enterovirus 71 (EV71) infection [Virus Research, 180 (2014), 1-11]. However, the exact function of this punctate pattern in cytoplasm during EV71 infection remains unknown. In this study, we firstly have examined this punctate pattern of Sam68 re-localization in the cytoplasm, and observed the obvious recruitments of Sam68 to the EV71-induced stress granules (SGs). Sam68, belongs to the KH domain family of RNA binding proteins (RBPs), was then confirmed that its KH domain was essential for this recruitment. Nevertheless, Knockdown of Sam68 expression using ShRNA had no effects on SGs assembly, indicating that Sam68 is not a constitutive component of the SGs during EV71 infection. Lastly, we investigated the importance of microtubulin transport to SGs aggregation, and revealed that microtubule depolymerization inhibited SGs formation, suggesting that EV71-induced SGs move throughout the cytoplasm in a microtubule-dependent manner. Taken together, these results illuminated that EV71 infections can induce SGs formation, and Sam68, as a SGs component, migrates alone with SGs dependent on intact microtubule upon the viral infections. These findings may provide novel underlying mechanism for delineating the role of SGs during EV71 infection.