Characterization of the porcine epidemic diarrhea virus codon usage bias

RNA editing sites and triplet usage in exomes of bat RNA virus genomes of the family Paramyxoviridae

Bats contain a diverse spectrum of viral species in their bodies. The RNA virus family Paramyxoviridae tends to infect several vertebrate species, which are accountable for a variety of devastating infections in both humans and animals. Viruses of this kind include measles, mumps, and Hendra. Some synonymous codons are favoured over others in mRNAs during gene-to-protein synthesis process. Such phenomenon is termed as codon usage bias (CUB). Our research emphasized many aspects that shape the CUB of genes in the Paramyxoviridae family found in bats. Here, the nitrogenous base A occurred the most. AT was found to be abundant in the coding sequences of the Paramyxoviridae family. RSCU data revealed that A or T ending codons occurred more frequently than predicted. Furthermore, 3 overrepresented codons (CAT, AGA, and GCA) and 7 underrepresented codons (CCG, TCG, CGC, CGG, CGT, GCG and ACG) were detected in the viral genomes. Correspondence analysis, neutrality plot, and parity plots highlight the combined impact of mutational pressure and natural selection on CUB. The neutrality plot of GC12 against GC3 yielded a regression coefficient value of 0.366, indicating that natural selection had a significant (63.4 %) impact. Moreover, RNA editing analysis was done, which revealed the highest frequency of C to T mutations. The results of our research revealed the pattern of codon usage and RNA editing sites in Paramyxoviridae genomes.

Comparative genomics revealed new insights into the plastome evolution of Ludwigia (Onagraceae, Myrtales)

The primrose-willow (Ludwigia L.), a well-defined genus of the Onagraceae family, comprises 87 species widely distributed worldwide. In this study, we sequenced and characterized the complete chloroplast (cp) genomes of three species in the genus, including Ludwigia adscendens, Ludwigia hyssopifolia, and Ludwigia prostrata. Three Ludwigia cp genomes ranged from 158,354 to 159,592 bp in size, and each contained 113 genes, including 79 unique protein-coding genes (PCGs), four rRNA genes, and 30 tRNA genes. A comparison of the Ludwigia cp genomes revealed that they were highly conserved in gene composition, gene orientation, and GC content. Moreover, we compared the structure of cp genomes and reconstructed phylogenetic relationships with related species in the Onagraceae family. Regarding contraction/expansion of inverted repeat (IR) region, two kinds of expansion IR region structures were found in Oenothera, Chamaenerion, and Epilobium genera, with primitive IR structures in Ludwigia and Circeae genera. The regions clpP, ycf2, and ycf1 genes possessed highly divergent nucleotides among all available cp genomes of the Onagraceae family. The phylogenetic reconstruction using 79 PCGs from 39 Onagraceae cp genomes inferred that Ludwigia (including L. adscendens, L. hyssopifolia, L. prostrata, and Ludwigia octovalvis) clade was monophyletic and well-supported by the bootstrap and posterior probability values. This study provides the reference cp genomes of three Ludwigia species, which can be used for species identification and phylogenetic reconstruction of Ludwigia and Onagraceae taxa.

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.

Relative synonymous codon usage and codon pair analysis of depression associated genes

Depression negatively impacts mood, behavior, and mental and physical health. It is the third leading cause of suicides worldwide and leads to decreased quality of life. We examined 18 genes available at the genetic testing registry (GTR) from the National Center for Biotechnological Information to investigate molecular patterns present in depression-associated genes. Different genotypes and differential expression of the genes are responsible for ensuing depression. The present study, investigated codon pattern analysis, which might play imperative roles in modulating gene expression of depression-associated genes. Of the 18 genes, seven and two genes tended to up- and down-regulate, respectively, and, for the remaining genes, different genotypes, an outcome of SNPs were responsible alone or in combination with differential expression for different conditions associated with depression. Codon context analysis revealed the abundance of identical GTG-GTG and CTG-CTG pairs, and the rarity of methionine-initiated codon pairs. Information based on codon usage, preferred codons, rare, and codon context might be used in constructing a deliverable synthetic construct to correct the gene expression level of the human body, which is altered in the depressive state. Other molecular signatures also revealed the role of evolutionary forces in shaping codon usage.

Analysis of codon usage patterns of porcine enteric alphacoronavirus and its host adaptability

Porcine enteric alphacoronavirus (PEAV) is a newly emerging swine enteropathogen that poses a threat to the swine industry. To understand the PEAV genome evolution, we performed a comprehensive analysis of the codon usage patterns in fifty-nine PEAV strains currently available. Phylogenetic analysis showed that PEAV can be divided into six lineages. Effective number of codons analysis demonstrated that the PEAV genome exhibits a low codon usage bias (CUB). Nucleotide composition analysis indicated that the PEAV genome has the most abundant nucleotide U content, with GC content (39.37% ± 0.08%) much lower than AU content (60.63% ± 0.08%). Neutrality and effective number of codons plot analyses suggested that natural selection rather than mutation pressure dominates the CUB of PEAV. Host adaptation analysis revealed that PEAV fits the codon usage pattern of non-human primates, humans and mice better than that of pigs. Our data enriches information on PEAV evolution, host adaptability, and cross-species transmission.

Description and phylogenetic analysis of the complete mitochondrial genome in Eulaelaps silvestris provides new insights into the molecular classification of the family Haemogamasidae

In this study, the mitochondrial genome of Eulaelaps silvestris, which parasitizes Apodemus chevrieri, was sequenced and assembled to fill the gap in understanding the molecular evolution of the genus Eulaelaps. The E. silvestris mitochondrial genome is a double-stranded DNA molecule with a length of 14 882 bp, with a distinct AT preference for base composition and a notably higher AT content than GC content. The arrangement between genes is relatively compact, with a total of 10 gene intergenic regions and 12 gene overlap regions. All protein-coding genes had a typical ATN initiation codon, and only 2 protein-coding genes had an incomplete termination codon T. Out of the 13 protein-coding genes, the 5 most frequently used codons ended in A/U, with only 1 codon ending in G/C had an relative synonymous codon usage value >1. Except for trnS1 and trnS2, which lacked the D arm, all other tRNAs were able to form a typical cloverleaf structure; and there were a total of 38 mismatches in the folding process of tRNA genes. Unlike the gene arrangement order of the arthropod hypothetical ancestor, the E. silvestris mitochondrial genome underwent fewer rearrangements, mainly near tRNA genes and control regions. Both the maximum likelihood tree and the Bayesian tree showed that the family Haemogamasidae is most closely related to the family Dermanyssidae. The results not only provide a theoretical basis for studying the phylogenetic relationships of the genus Eulaelaps, but also provide molecular evidence that the family Haemogamasidae does not belong to the subfamily Laelapidae.

Coronavirus accessory protein ORF3 biology and its contribution to viral behavior and pathogenesis

Coronavirus porcine epidemic diarrhea virus (PEDV) is classified in the genus Alphacoronavirus, family Coronaviridae that encodes the only accessory protein, ORF3 protein. However, how ORF3 contributes to viral pathogenicity, adaptability, and replication is obscure. In this review, we summarize current knowledge and identify gaps in many aspects of ORF3 protein in PEDV, with emphasis on its unique biological features, including membrane topology, Golgi retention mechanism, potential intrinsic disordered property, functional motifs, protein glycosylation, and codon usage phenotypes related to genetic evolution and gene expression. In addition, we propose intriguing questions related to ORF3 protein that we hope to stimulate further studies and encourage collaboration among virologists worldwide to provide constructive knowledge about the unique characteristics and biological functions of ORF3 protein, by which their potential role in clarifying viral behavior and pathogenesis can be possible.

Molecular and Structural Evolution of Porcine Epidemic Diarrhea Virus

To analyze the evolutionary characteristics of the highly contagious porcine epidemic diarrhea virus (PEDV) at the molecular and structural levels, we analyzed the complete genomes of 647 strains retrieved from the GenBank database. The results showed that the spike (S) gene exhibited larger dS (synonymous substitutions per synonymous site) values than other PEDV genes. In the selective pressure analysis, eight amino acid (aa) sites of the S protein showed strong signals of positive selection, and seven of them were located on the surface of the S protein (S1 domain), suggesting a high selection pressure of S protein. Topologically, the S gene is more representative of the evolutionary relationship at the genome-wide level than are other genes. Structurally, the evolutionary pattern is highly S1 domain-related. The haplotype networks of the S gene showed that the strains are obviously clustered geographically in the lineages corresponding to genotypes GI and GII. The alignment analysis on representative strains of the main haplotypes revealed three distinguishable nucleic acid sites among those strains, suggesting a putative evolutionary mechanism in PEDV. These findings provide several new fundamental insights into the evolution of PEDV and guidance for developing effective prevention countermeasures against PEDV.

Codon Usage and Context Analysis of Genes Modulated during SARS-CoV-2 Infection and Dental Inflammation

The overexpression of SARS-CoV-2 primary receptors and co-receptors (ACE2, TMPRSS2, FURIN, and CD147) enhance the likeliness of SARS-CoV-2 infection. The genes for same receptors are overexpressed in the periodontal tissues of periodontitis patients. On the other hand, BMAL1 is recognized to play a crucial role in regulating pulmonary inflammation and enhancing susceptibility to viral infection. Silenced BMAL1 disrupts circadian transcriptional regulations, enhances vulnerability to SARS-CoV-2 infections, and may trigger the further production of TNF-α and other pro-inflammatory cytokines that propagate the cytokine storm and exacerbate periodontal inflammation. Therefore ACE2, TMPRSS2, FURIN, CD147, and BMAL1 are the crossroads between SARS-CoV-2 and Periodontitis genes. The enhanced expression of ACE2, TMPRSS2, FURIN, and CD147 and the diminished expression of BMAL1 may be a strategy to check both ailments simultaneously. In gene manipulation techniques, oligos are introduced, which contain all the necessary information to manipulate gene expression. The data are derived from the studies on genes' molecular patterns, including nucleotide composition, dinucleotide patterns, relative synonymous codon usage, codon usage bias, codon context, and rare and abundant codons. Such information may be used to manipulate the overexpression and underexpression of the genes at the time of SARS-CoV-2 infection and periodontitis to mitigate both ailments simultaneously; it can be explored to uncover possible future treatments.

Host adaptation of codon usage in SARS-CoV-2 from mammals indicates potential natural selection and viral fitness

SARS-CoV-2 infection, which is the cause of the COVID-19 pandemic, has expanded across various animal hosts, and the virus can be transmitted particularly efficiently in minks. It is still not clear how SARS-CoV-2 is selected and evolves in its hosts, or how mutations affect viral fitness. In this report, sequences of SARS-CoV-2 isolated from human and animal hosts were analyzed, and the binding energy and capacity of the spike protein to bind human ACE2 and the mink receptor were compared. Codon adaptation index (CAI) analysis indicated the optimization of viral codons in some animals such as bats and minks, and a neutrality plot demonstrated that natural selection had a greater influence on some SARS-CoV-2 sequences than mutational pressure. Molecular dynamics simulation results showed that the mutations Y453F and N501T in mink SARS-CoV-2 could enhance the binding of the viral spike to the mink receptor, indicating the involvement of these mutations in natural selection and viral fitness. Receptor binding analysis revealed that the mink SARS-CoV-2 spike interacted more strongly with the mink receptor than the human receptor. Tracking the variations and codon bias of SARS-CoV-2 is helpful for understanding the fitness of the virus in virus transmission, pathogenesis, and immune evasion.

Strategies and Patterns of Codon Bias in Molluscum Contagiosum Virus

Trends associated with codon usage in molluscum contagiosum virus (MCV) and factors governing the evolution of codon usage have not been investigated so far. In this study, attempts were made to decipher the codon usage trends and discover the major evolutionary forces that influence the patterns of codon usage in MCV with special reference to sub-types 1 and 2, MCV-1 and MCV-2, respectively. Three hypotheses were tested: (1) codon usage patterns of MCV-1 and MCV-2 are identical; (2) SCUB (synonymous codon usage bias) patterns of MCV-1 and MCV-2 slightly deviate from that of human host to avoid affecting the fitness of host; and (3) translational selection predominantly shapes the SCUB of MCV-1 and MCV-2. Various codon usage indices viz. relative codon usage value, effective number of codons and codon adaptation index were calculated to infer the nature of codon usage. Correspondence analysis and correlation analysis were performed to assess the relative contribution of silent base contents and significance of codon usage indices in defining bias in codon usage. Among the tested hypotheses, only the second and third hypotheses were accepted.

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.

Genes Common in Primary Immunodeficiencies and Cancer Display Overrepresentation of Codon CTG and Dominant Role of Selection Pressure in Shaping Codon Usage

Primary immunodeficiencies (PIDs) are disorders of the immune system that involve faulty cellular, humoral, or both cellular and humoral functions. PIDs are present at the crossroad between infections, immune dysregulation, and cancers. A panel encompassing 42 genes involved in both PIDs and cancer has been investigated for the genes’ compositional properties, codon usage patterns, various forces affecting codon choice, protein properties, and gene expression profiles. In the present study, the codon choice of genes was found to be dependent upon the richness of the nucleotide; the viz AT nucleotide rich genome preferred AT ending codons. The dinucleotide TpA adversely affected protein expression, while CpG did not. The CTG codon was the most overrepresented codon in 80.95% of genes. Analysis of various protein properties, including GRAVY, AROMA, isoelectric point, aliphatic index, hydrophobicity, instability index, and numbers of acidic, basic, and neutral amino acid residues revealed that the hydrophobicity index, instability index, and numbers of acidic and basic amino acid residues are the factors affecting gene expression. Based on neutrality analysis, parity analysis, ENc-GC3 analysis, and regression analysis of nucleotides present at the first and third positions of the codon, it was determined that selection pressure, mutation pressure, and compositional constraints all participated in shaping codon usage. The study will help determine the various evolutionary forces acting on genes common to both PIDs and cancer. Codon usage analysis might be helpful in the future to augment both diseases simultaneously. The research also indicates a peculiar pattern adapted by a set of genes involved in any disease.

Analysis of codon usage patterns of the chloroplast genome in Delphinium grandiflorum L. reveals a preference for AT-ending codons as a result of major selection constraints

Background

Codon usage bias analysis is a suitable strategy for identifying the principal evolutionary driving forces in different organisms. Delphinium grandiflorum L. is a perennial herb with high economic value and typical biological characteristics. Evolutionary analysis of D. grandiflorum can provide a rich resource of genetic information for developing hybridization resources of the genus Delphinium.

Methods

Synonymous codon usage (SCU) and related indices of 51 coding sequences from the D. grandiflorum chloroplast (cp) genome were calculated using Codon W, Cups of EMBOSS, SPSS and Microsoft Excel. Multivariate statistical analysis combined by principal component analysis (PCA), correspondence analysis (COA), PR2-plot mapping analysis and ENC plot analysis was then conducted to explore the factors affecting the usage of synonymous codons.

Results

The SCU bias of D. grandiflorum was weak and codons preferred A/T ending. A SCU imbalance between A/T and G/C at the third base position was revealed by PR2-plot mapping analysis. A total of eight codons were identified as the optimal codons. The PCA and COA results indicated that base composition (GC content, GC₃ content) and gene expression were important for SCU bias. A majority of genes were distributed below the expected curve from the ENC plot analysis and up the standard curve by neutrality plot analysis. Our results showed that with the exception of notable mutation pressure effects, the majority of genetic evolution in the D. grandiflorum cp genome might be driven by natural selection.

Discussions

Our results provide a theoretical foundation for elucidating the genetic architecture and mechanisms of D. grandiflorum, and contribute to enriching D. grandiflorum genetic resources.

Comprehensive Analysis of Synonymous Codon Usage Bias for Complete Genomes and E2 Gene of Atypical Porcine Pestivirus

Atypical porcine pestivirus (APPV) is an emerging novel pestivirus causing the congenital tremor (CT) in piglets. The worldwide distribution characteristic of APPV make it a threat to global swine health. E2 is the major envelope glycoprotein of APPV and the crucial target for vaccine development. Considering the genetic variability of APPV complete genomes and its E2 gene as well as gaps for codon analysis, a comprehensive analysis of codon usage patterns was performed. Relative synonymous codon usage (RSCU) and effective number of codon (ENC) analyses showed that a relatively instable change existed and a slight low codon usage bias (CUB) were displayed in APPV genomes. ENC-plot analysis and correlation analyses of nucleotide compositions and ENC showed that mutation pressure and natural selection both affected the codon usage bias of the APPV and natural selection had a more obvious influence for E2 gene compared with complete genomes. Principal component analysis (PCA) and correlation analyses confirmed the above results. Correlation analyses between Gravy and Aromaticity values and the codon bias showed that natural selection played an important role in shaping the synonymous codon bias. Furthermore, neutrality plot analysis showed that natural selection was the main force while mutation pressure was a minor force influencing the codon usage pattern of the APPV E2 gene and complete genomes. The results could illustrate the codon usage patterns of APPV genomes and provided valuable basic data for further fundamental research of evolution of APPV.

Biological Characteristics and Patterns of Codon Usage Evolution for the African Genotype Zika Virus

We investigated temporal trends of codon usage changes for different host species to determine their importance in Zika virus (ZIKV) evolution. Viral spillover resulting from the potential of codon adaptation to host genome was also assessed for the African genotype ZIKV in comparison to the Asian genotype. To improve our understanding on its zoonotic maintenance, we evaluated in vitro the biological properties of the African genotype ZIKV in vertebrate and mosquito cell lines. Analyses were performed in comparison to Yellow fever virus (YFV). Despite significantly lower codon adaptation index trends than YFV, ZIKV showed evident codon adaptation to vertebrate hosts, particularly for the green African monkey Chlorocebus aethiops. PCA and CAI analyses at the individual ZIKV gene level for both human and Aedes aegypti indicated a clear distinction between the two genotypes. African ZIKV isolates showed higher virulence in mosquito cells than in vertebrate cells. Their higher replication in mosquito cells than African YFV confirmed the role of mosquitoes in the natural maintenance of the African genotype ZIKV. An analysis of individual strain growth characteristics indicated that the widely used reference strain MR766 replicates poorly in comparison to African ZIKV isolates. The recombinant African Zika virus strain ArD128000*^E/NS5 may be a good model to include in studies on the mechanism of host tropism, as it cannot replicate in the tested vertebrate cell line.

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.

Codon Usage Bias Analysis of Bluetongue Virus Causing Livestock Infection

Bluetongue virus (BTV) is a double-stranded RNA virus with multiple segments and belongs to the genus Orbivirus within the family Reoviridae. BTV is spread to livestock through its dominant vector, biting midges of genus Culicoides. Although great progress has been made in genomic analyses, it is not fully understood how BTVs adapt to their hosts and evade the host's immune systems. In this study, we retrieved BTV genome sequences from the National Center for Biotechnology Information (NCBI) database and performed a comprehensive research to explore the codon usage patterns in 50 BTV strains. We used bioinformatic approaches to calculate the relative synonymous codon usage (RSCU), codon adaptation index (CAI), effective number of codons (ENC), and other indices. The results indicated that most of the overpreferred codons had A-endings, which revealed that mutational pressure was the major force shaping codon usage patterns in BTV. However, the influence of natural selection and geographical factors cannot be ignored on viral codon usage bias. Based on the RSCU values, we performed a comparative analysis between BTVs and their hosts, suggesting that BTVs were inclined to evolve their codon usage patterns that were comparable to those of their hosts. Such findings will be conducive to understanding the elements that contribute to viral evolution and adaptation to hosts.

Analysis of preferred codon usage in the coronavirus N genes and their implications for genome evolution and vaccine design

•

The nucleotide variations among the N genes of 13 different coronaviruses (CoVs) were interpreted.

•

Overall, 18 amino acids observed with varying preferred codons.

•

The effective number of codon values ranged from 40.43 to 53.85, revealing a slight codon bias.

•

A highly significant correlation between GC3s and ENc values was observed in porcine epidemic diarrhea CoV, followed by Middle East respiratory syndrome CoV.

The nucleocapsid (N) protein of a coronavirus plays a crucial role in virus assembly and in its RNA transcription. It is important to characterize a virus at the nucleotide level to discover the virus’s genomic sequence variations and similarities relative to other viruses that could have an impact on the functions of its genes and proteins. This entails a comprehensive and comparative analysis of the viral genomes of interest for preferred nucleotides, codon bias, nucleotide changes at the 3^rd position (NT3s), synonymous codon usage and relative synonymous codon usage. In this study, the variations in the N proteins among 13 different coronaviruses (CoVs) were analysed at the nucleotide and amino acid levels in an attempt to reveal how these viruses adapt to their hosts relative to their preferred codon usage in the N genes. The results revealed that, overall, eighteen amino acids had different preferred codons and eight of these were over-biased. The N genes had a higher AT% over GC% and the values of their effective number of codons ranged from 40.43 to 53.85, indicating a slight codon bias. Neutrality plots and correlation analyses showed a very high level of GC3s/GC correlation in porcine epidemic diarrhea CoV (pedCoV), followed by Middle East respiratory syndrome-CoV (MERS CoV), porcine delta CoV (dCoV), bat CoV (bCoV) and feline CoV (fCoV) with r values 0.81, 0.68, -0.47, 0.98 and 0.58, respectively. These data implied a high rate of evolution of the CoV genomes and a strong influence of mutation on evolutionary selection in the CoV N genes. This type of genetic analysis would be useful for evaluating a virus’s host adaptation, evolution and is thus of value to vaccine design strategies.

Codon Usage in the Iflaviridae Family Is Not Diverse Though the Family Members Are Isolated from Diverse Host Taxa

All iflavirus members belong to the unique genus, Iflavirus, of the family, Iflaviridae. The host taxa and sequence identities of these viruses are diverse. A codon usage bias, maintained by a balance between selection, mutation, and genetic drift, exists in a wide variety of organisms. We characterized the codon usage patterns of 44 iflavirus genomes that were isolated from the classes, Insecta, Arachnida, Mammalia, and Malacostraca. Iflaviruses lack a strong codon usage bias when they are evaluated using an effective number of codons. The odds ratios of the majority of dinucleotides are within the normal range. However, the dinucleotides at the 1st-2nd codon positions are more biased than those at the 2nd-3rd codon positions. Plots of effective numbers of codons, relative neutrality analysis, and PR2 bias analysis all indicate that selection pressure dominates mutations in shaping codon usage patterns in the family, Iflaviridae. When these viruses were grouped into their host taxa, we found that the indices, including the nucleotide composition, effective number of codons, relative synonymous codon usage, and the influencing factors behind the codon usage patterns, all show that there are non-significant differences between the six host-taxa-groups. Our results disagree with our assumption that diverse viruses should possess diverse codon usage patterns, suggesting that the nucleotide composition and codon usage in the family, Iflaviridae, are not host taxa-specific signatures.

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.

Comprehensive codon usage analysis of porcine deltacoronavirus

Porcine deltacoronavirus (PDCoV) is a newly identified coronavirus of pigs that was first reported in Hong Kong in 2012. Since then, many PDCoV isolates have been identified worldwide. In this study, we analyzed the codon usage pattern of the S gene using complete coding sequences and complete PDCoV genomes to gain a deeper understanding of their genetic relationships and evolutionary history. We found that during evolution three groups evolved with a relatively low codon usage bias (effective number of codons (ENC) of 52). The factors driving bias were complex. However, the primary element influencing the codon bias of PDCoVs was natural selection. Our results revealed that different natural environments may have a significant impact on the genetic characteristics of the strains. In the future, more epidemiological surveys are required to examine the factors that resulted in the emergence and outbreak of this virus.

Development and validation of a monoclonal antibody-based competitive ELISA for detection of antibodies against porcine epidemic diarrhoea virus (PEDV)

Porcine epidemic diarrhoea virus (PEDV), belongs to the genus Alphacoronavirus in the family Coronaviridae and causes severe diarrhoea, vomiting, dehydration and high mortality in seronegative newborn piglets. Thus, a precise and rapid diagnosis of PEDV infection is important for the application of control measures to limit viral dissemination. In this investigation, a monoclonal antibodies (MAbs)-based competitive enzyme-linked immunosorbent assay (ELISA) for detecting antibodies against PEDV was developed and validated. The diagnostic performance of the test was evaluated by receiver operating characteristic (ROC) analysis using a panel of 829 known sera collected from different commercial pig farms, with or without a history of PED presence and from experimentally infected pigs. The competitive ELISA showed excellent diagnostic performance and discriminatory power with high sensitivity (Se) and specificity (Sp) values (Se = 96.5%, 95% IC 94.1-98.1; Sp = 98.2%, 95% IC 96.3-99.3). Moreover, this competitive ELISA method has other properties, such as high feasibility of testing sera without pre-treatment and automatic and instrument-mediated revealing, that makes it appropriate for large-scale screenings of affected pig farms in endemic regions or for monitoring plans in PEDV-free areas.

Comprehensive Analysis of Codon Usage on Rabies Virus and Other Lyssaviruses

Rabies virus (RABV) and other lyssaviruses can cause rabies and rabies-like diseases, which are a persistent public health threat to humans and other mammals. Lyssaviruses exhibit distinct characteristics in terms of geographical distribution and host specificity, indicative of a long-standing diversification to adapt to the environment. However, the evolutionary diversity of lyssaviruses, in terms of codon usage, is still unclear. We found that RABV has the lowest codon usage bias among lyssaviruses strains, evidenced by its high mean effective number of codons (ENC) (53.84 ± 0.35). Moreover, natural selection is the driving force in shaping the codon usage pattern of these strains. In summary, our study sheds light on the codon usage patterns of lyssaviruses, which can aid in the development of control strategies and experimental research.

Analysis of the codon usage pattern in Middle East Respiratory Syndrome Coronavirus

Middle East Respiratory Syndrome Coronavirus (MERS-CoV), which first broken out in Jeddah in 2012, causes a severe acute respiratory illness with a high mortality rate. To better understand the molecular characteristics of isolated MERS-CoV genomes, we first analysed the codon usage pattern of the zoonotic MERS-CoV strains comprehensively to gain an insight into the mechanism of cross-species transmission. We found that MERS human/camel isolates showed a low codon usage bias. Both mutation and nature selection pressure have contributed to this low codon usage bias, with the former being the main determining factor. We also observed that gene function, evolution time and the different host species of the virus all contributed to the bias of MERS-CoV, to some extent. Additionally, the codon usage pattern of MERS-CoV isolates is different from other related Nidovirales viruses isolated from bats and hedgehogs. In the future, more epidemiological surveys are required to examine the factors that resulted in the emergence and outbreak of this virus.

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

Senecavirus A (SVA) infection was recently confirmed in pigs in Brazil, United States of America and Canada. To better understand the molecular characteristics of isolated SVA genomes, we first reported genome-wide comprehensive analyses of codon usage and various factors that have contribute to the molecular evolution in SVA. The effective number of codons (ENC) ranged from 54.51 to 55.54 with an average of 54.87 ± 0.285, which reveals a relatively stable nucleotide composition. We found that codon usage bias of the SVA was low. Mutational pressure acted as an increasingly dominant factor for the evolution of the virus compared with the natural selection. Notably, codon usage bias was also affected by the geographic distribution and isolated time. The first systemic analysis on the codon usage bias of the SVA provides important information for the understanding of the evolution of the SVA and has fundamental and theoretical benefits.

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.

Codon usage bias in the N gene of rabies virus

Since its emergence, rabies virus (RABV) has been a major worldwide concern especially in developing countries. The nucleoprotein (N) of RABV is highly conserved and key for genetic typing, thus a better understanding of the N gene evolutionary trajectory can assist the development of control measures. We found that the N gene of RABV has a low codon usage bias with a mean effective number of codons (ENC) value of 56.33 influenced by both mutation pressure and natural selection. However, neutrality analysis indicated that natural selection dominates over mutation pressure. Additionally, we found that dinucleotide bias partly contributed to RABV codon usage bias. On the other hand, based on the clades of phylogenetic tree, we found that the evolutionary rate of the Africa 2 clade was the highest with a mean value of 3.75×10^-3 substitutions per site per year. Above all, our results regarding N gene of RABV codon usage will serve future RABV evolution research.

Evolutionary and genetic analysis of the VP2 gene of canine parvovirus

BACKGROUND

Canine parvovirus (CPV) type 2 emerged in 1978 in the USA and quickly spread among dog populations all over the world with high morbidity. Although CPV is a DNA virus, its genomic substitution rate is similar to some RNA viruses. Therefore, it is important to trace the evolution of CPV to monitor the appearance of mutations that might affect vaccine effectiveness.

RESULTS

Our analysis shows that the VP2 genes of CPV isolated from 1979 to 2016 are divided into six groups: GI, GII, GIII, GIV, GV, and GVI. Amino acid mutation analysis revealed several undiscovered important mutation sites: F267Y, Y324I, and T440A. Of note, the evolutionary rate of the CPV VP2 gene from Asia and Europe decreased. Codon usage analysis showed that the VP2 gene of CPV exhibits high bias with an ENC ranging from 34.93 to 36.7. Furthermore, we demonstrate that natural selection plays a major role compared to mutation pressure driving CPV evolution.

CONCLUSIONS

There are few studies on the codon usage of CPV. Here, we comprehensively studied the genetic evolution, codon usage pattern, and evolutionary characterization of the VP2 gene of CPV. The novel findings revealing the evolutionary process of CPV will greatly serve future CPV research.

Genome-wide analysis of codon usage bias in Bovine Coronavirus

BACKGROUND

Bovine coronavirus (BCoV) belong to the genus Betacoronavirus of the family Coronaviridae. BCoV are widespread around the world and cause enteric or respiratory infections among cattle, leading to important economic losses to the beef and dairy industry worldwide. To study the relation of codon usage among viruses and their hosts is essential to understand host-pathogen interaction, evasion from host's immune system and evolution.

METHODS

We performed a comprehensive analysis of codon usage and composition of BCoV.

RESULTS

The global codon usage among BCoV strains is similar. Significant differences of codon preferences in BCoV genes in relation to codon usage of Bos taurus host genes were found. Most of the highly frequent codons are U-ending. G + C compositional constraint and dinucleotide composition also plays a role in the overall pattern of BCoV codon usage.

CONCLUSIONS

The results of these studies revealed that mutational bias is a leading force shaping codon usage in this virus. Additionally, relative dinucleotide frequencies, geographical distribution, and evolutionary processes also influenced the codon usage pattern.

Analysis of the codon usage of the ORF2 gene of feline calicivirus

Feline calicivirus (FCV) is a highly prevalent pathogen of the domestic cat that causes acute infections of the oral and upper respiratory tract. The E region of the ORF2 protein is responsible for the induction of virus-neutralizing antibodies, thus it is important to understand the codon usage of this gene. Here, analysed 90 coding sequences of ORF2 and show that it undergoes a low codon usage bias. In addition, although mutational bias is one of the factors shaping the codon usage bias of this gene, natural selection plays a more significant role. Our results reveal part of the mechanisms driving FCV evolution, which will lay foundation for the further research of FCV.

•

The codon usage of the ORF2 gene of FCV was first analysed in this study.

•

The ORF2 gene of FCV shows a low codon usage bias.

A detailed analysis of codon usage patterns and influencing factors in Zika virus

Recent outbreaks of Zika virus (ZIKV) in Africa, Latin America, Europe, and Southeast Asia have resulted in serious health concerns. To understand more about evolution and transmission of ZIKV, detailed codon usage analysis was performed for all available strains. A high effective number of codons (ENC) value indicated the presence of low codon usage bias in ZIKV. The effect of mutational pressure on codon usage bias was confirmed by significant correlations between nucleotide compositions at third codon positions and ENCs. Correlation analysis between Gravy values, Aroma values and nucleotide compositions at third codon positions also indicated some influence of natural selection. However, the low codon adaptation index (CAI) value of ZIKV with reference to human and mosquito indicated poor adaptation of ZIKV codon usage towards its hosts, signifying that natural selection has a weaker influence than mutational pressure. Additionally, relative dinucleotide frequencies, geographical distribution, and evolutionary processes also influenced the codon usage pattern to some extent.

Codon usage in Alphabaculovirus and Betabaculovirus hosted by the same insect species is weak, selection dominated and exhibits no more similar patterns than expected

Mutations shape synonymous codon usage bias in certain organism genomes, while selection shapes it in others. Lepidopteran-specific Alphabaculovirus and Betabaculovirus are two large genera in the family of Baculoviridae. In this study, we analyzed the codon usage patterns in 17 baculoviruses, including 10 alphabaculoviruses and 7 betabaculoviruses, which were isolated from seven insect species, and we characterized the codon usage patterns between Alphabaculovirus and Betabaculovirus. Our results show that all the baculoviruses possessed a general weak trend of codon bias. The differences of ENc (effective number of codons) values, nucleotide contents and the impacts of nucleotide content on ENc value within alpha-/betabaculovirus pairs were independent of whether the host species are the same or different. Furthermore, the majority of amino acid sequences adopted codons unequally in all viruses, but the numbers of common preferred codons between alpha- and betabaculoviruses hosted by the same insect species were not significantly different from the differences observed between alpha- and betabaculoviruses hosted by different insect species. In addition, the amino acids that adopt the same synonymous codon composition between alpha- and betabaculoviruses hosted by the same insect species were statistically as few as those between alpha- and betabaculoviruses hosted by different insect species. Correspondence analysis revealed that no major factors resulted in the codon bias in these baculoviruses, implying multiple minor influential factors exist. Neutrality plot analysis indicated that selection pressure dominated mutations in shaping the codon usage. However, the levels of selection pressure were not significantly different among viruses hosted by the same insect species. We expect that evolution would cause the alpha- and betabaculoviruses hosted by the same insect species to share more patterns, but this effect was not observed.

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.