Genetic and codon usage bias analyses of major capsid protein gene in Ranavirus

Analysis of codon usage bias of exonuclease genes in invertebrate iridescent viruses

Invertebrate iridescent viruses (IIVs) are double-stranded DNA viruses that belong to the Iridoviridae family. IIVs result diseases that vary in severity from subclinical to lethal in invertebrate hosts. Codon usage bias (CUB) analysis is a versatile method for comprehending the genetic and evolutionary aspects of species. In this study, we analyzed the CUB in 10 invertebrate iridescent viruses exonuclease genes by calculating and comparing the nucleotide contents, effective number of codons (ENC), codon adaptation index (CAI), relative synonymous codon usage (RSCU), and others. The results revealed that IIVs exonuclease genes are rich in A/T. The ENC analysis displayed a low codon usage bias in IIVs exonuclease genes. ENC-plot, neutrality plot, and parity rule 2 plot demonstrated that besides mutational pressure, other factors like natural selection, dinucleotide content, and aromaticity also contributed to CUB. The findings could enhance our understanding of the evolution of IIVs exonuclease genes.

Circular replication-associated protein-encoding single-stranded DNA virus with risk of spillover is widely prevalent in domestic animals in China

Circular replication-associated protein (Rep)-encoding single-stranded (CRESS) DNA viruses are highly diverse and have a broad range of hosts. In this study, we report the detection of Bo-Circo-like virus AH20-1 in the feces of diarrheal cattle. The virus has a circular genome of 3,912 nucleotides, three major putative open reading frames, and encodes a Rep gene of 310 amino acids. We found that the virus is closely related to the Bo-Circo-like virus CH strain, which belongs to the novel Kirkoviridae family. Furthermore, we conducted a nationwide surveillance program and found that the virus is prevalent in China (23.6%, 205/868), with the BCLa subtype being the predominant strain. Our findings suggest that the virus can infect sheep, highlighting the potential for cross-species transmission. Our pressure analysis indicates that the CRESS-DNA Kirkoviridae family has broad host adaptation, and that selection pressure played an important role in the evolution of its Rep genes. Our study underscores the need for continued epidemiological surveillance of this virus due to its widespread prevalence in our ruminant population and potential for cross-species transmission.

Sign of APOBEC editing, purifying selection, frameshift, and in-frame nonsense mutations in the microevolution of lumpy skin disease virus

The lumpy skin disease virus (LSDV), which mostly affects ruminants and causes huge-economic loss, was endemic in Africa, caused outbreaks in the Middle East, and was recently detected in Russia, Serbia, Greece, Bulgaria, Kazakhstan, China, Taiwan, Vietnam, Thailand, and India. However, the role of evolutionary drivers such as codon selection, negative/purifying selection, APOBEC editing, and genetic variations such as frameshift and in-frame nonsense mutations in the LSDVs, which cause outbreaks in cattle in various countries, are still largely unknown. In the present study, a frameshift mutation in LSDV035, LSDV019, LSDV134, and LSDV144 genes and in-frame non-sense mutations in LSDV026, LSDV086, LSDV087, LSDV114, LSDV130, LSDV131, LSDV145, LSDV154, LSDV155, LSDV057, and LSDV081 genes were revealed among different clusters. Based on the available complete genome sequences, the prototype wild-type cluster-1.2.1 virus has been found in other than Africa only in India, the wild-type cluster-1.2.2 virus found in Africa were spread outside Africa, and the recombinant viruses spreading only in Asia and Russia. Although LSD viruses circulating in different countries form a specific cluster, the viruses detected in each specific country are distinguished by frameshift and in-frame nonsense mutations. Furthermore, the present study has brought to light that the selection pressure for codons usage bias is mostly exerted by purifying selection, and this process is possibly caused by APOBEC editing. Overall, the present study sheds light on microevolutions in LSDV, expected to help in future studies towards disturbed ORFs, epidemiological diagnostics, attenuation/vaccine reverts, and predicting the evolutionary direction of LSDVs.

Codon usage bias analysis of the gene encoding NAD+-dependent DNA ligase protein of Invertebrate iridescent virus 6

The genome of Invertebrate iridescent virus 6 (IIV6) contains a sequence that shows similarity to eubacterial NAD+-dependent DNA ligases. The 615-amino acid open reading frame (ORF 205R) consists of several domains, including an N-terminal domain Ia, followed by an adenylation domain, an OB-fold domain, a helix-hairpin-helix (HhH) domain, and a BRCT domain. Notably, the zinc finger domain, typically present in NAD+-dependent DNA ligases, is absent in ORF 205R. Since the protein encoded by ORF 205R (IIV6 DNA ligase gene) is involved in critical functions such as DNA replication, modification, and repair, it is crucial to comprehend the codon usage associated with this gene. In this paper, the codon usage bias (CUB) in DNA ligase gene of IIV6 and 11 reference iridoviruses was analyzed by comparing the nucleotide contents, relative synonymous codon usage (RSCU), effective number of codons (ENC), codon adaptation index (CAI), relative abundance of dinucleotides and other indices. Both the base content and the RCSU analysis indicated that the A- and T-ending codons were mostly favored in the DNA ligase gene of IIV6. The ENC value of 35.64 implied a high CUB in the IIV6 DNA ligase gene. The ENC plot, neutrality plot, parity rule 2 plot, correspondence analysis revealed that mutation pressure and natural selection had an impact on the CUB of the IIVs DNA ligase genes. Additionally, the analysis of codon adaptation index demonstrated that the IIV6 DNA ligase gene is strongly adapted to its host. These findings will improve our comprehension of the CUB of IIV6 DNA ligase and reference genes, which may provide the required information for a fundamental evolutionary analysis of these genes.

Genetic Diversity and Characterization of Circular Replication (Rep)-Encoding Single-Stranded (CRESS) DNA Viruses

The CRESS-DNA viruses are the ubiquitous virus detected in almost all eukaryotic life trees and play an essential role in the maintaining ecosystem of the globe. Still, their genetic diversity is not fully understood. Here, we bring to light the genetic diversity of replication (Rep) and capsid (Cap) proteins of CRESS-DNA viruses. We divided the Rep protein of the CRESS-DNA virus into 10 clusters using CLANS and phylogenetic analyses. Also, most of the Rep protein in Rep cluster 1 (R1) and R2 (Circoviridae, Smacoviridae, Nanoviridae, and CRESSV1-5) contain the Viral_Rep superfamily and P-loop_NTPase superfamily domains, while the Rep protein of viruses in other clusters has no such characterized functional domain. The Circoviridae, Nanoviridae, and CRESSV1-3 viruses contain two domains, such as Viral_Rep and P-loop_NTPase; the CRESSV4 and CRESSV5 viruses have only the Viral_Rep domain; most of the sequences in the pCRESS-related group have only P-loop_NTPase; and Smacoviridae do not have these two domains. Further, we divided the Cap protein of the CRESS-DNA virus into 20 clusters using CLANS and phylogenetic analyses. The Rep and Cap proteins of Circoviridae and Smacoviridae are grouped into a specific cluster. Cap protein of CRESS-DNA viruses grouped with one cluster and Rep protein with another cluster. Further, our study reveals that selection pressure plays a significant role in the evolution of CRESS-DNA viruses' Rep and Cap genes rather than mutational pressure. We hope this study will help determine the genetic diversity of CRESS-DNA viruses as more sequences are discovered in the future. IMPORTANCE The genetic diversity of CRESS-DNA viruses is not fully understood. CRESS-DNA viruses are classified as CRESSV1 to CRESSV6 using only Rep protein. This study revealed that the Rep protein of the CRESS-DNA viruses is classified as CRESSV1 to CRESSV6 groups and the new Smacoviridae-related, CRESSV2-related, pCRESS-related, Circoviridae-related, and 1 to 4 outgroups, according to the Viral_Rep and P-loop_NTPase domain organization, CLANS, and phylogenetic analysis. Furthermore, for the first time in this study, the Cap protein of CRESS-DNA viruses was classified into 20 distinct clusters by CLANS and phylogenetic analysis. Through this classification, the genetic diversity of CRESS-DNA viruses clarifies the possibility of recombinations in Cap and Rep proteins. Finally, it has been shown that selection pressure plays a significant role in the evolution and genetic diversity of Cap and Rep proteins. This study explains the genetic diversity of CRESS-DNA viruses and hopes that it will help classify future detected viruses.

Evolution of monkeypox virus from 2017 to 2022: In the light of point mutations

Monkeypox virus (MPXV) causing multi-country outbreak-2022 is related to viruses caused outbreak-2017-2018 in West Africa. Still not fully understood which proteins of the MPXV discovered in Nigeria in 2017 have mutated through different lineages to the extent that it could cause a multi-country outbreak in 2022; similarly, codon usage bias, host adaptation indices, and the role of selection or mutation pressure in the mutated genes are also not fully studied. Here we report that according to the available sequence data this monkeypox virus acquires point mutations in multiple proteins in each period, and these point mutations accumulate and become a virus that can root outbreak-2022. Viruses exported from Nigeria to Singapore, Israel, and the United Kingdom in 2018-2019 were developed as evolutionary ancestors to B.1 viruses (MPXVs causing multi-country outbreak-2022) through MPXV/United States/2021/MD virus. Although these exported viruses have different amino acid mutations in different proteins, amino acid mutations in 10 proteins are common among them. The MPXV-United Kingdom-P2 virus evolved with only mutations in these 10 proteins and further evolved into MPXV/United States/2021/MD with amino acid mutations in 26 (including amino acid mutations in 10 proteins of the MPXV-United States-P2) proteins. It is noteworthy that specific amino acid mutations in these 22/26 (presence in MPXV/United States/2021/MD) proteins are present in B.1 viruses. Further, analysis of Relative Synonymous Codon Usage (RSCU), Synonymous Codon Usage Fraction (SCUF), and Effective Number of Codons (ENc) revealed codon usage bias in genes that exhibited nucleotide mutations in lineage B.1. Also, host adaptation indices analyzes such as Codon Adaptation Index (CAI), Expected-CAI (eCAI), Relative Codon Deoptimization Index (RCDI) and Expected value for the RCDI (eRCDI) analyzes reveal that the genes that demonstrated nucleotide mutations in lineage B.1 are favorable for human adaptation. Similarly, ENc-GC3s plot, Neutrality plot, and Parity Rule 2 (PR2)-bias plot analyzes suggest a major role of selection pressure than mutation pressure in the evolution of genes displaying nucleotide mutations in lineage B.1. Overall, from 2017 to 2022, MPXV's mutation and spread suggests that this virus continues to evolve through point mutation in the genes according to the available sequence data.

SARS-CoV-2 gained a novel spike protein S1-N-Terminal Domain (S1-NTD)

The clue behind the SARS-CoV-2 origin is still a matter of debate. Here, we report that SARS-CoV-2 has gained a novel spike protein S1-N-terminal domain (S1-NTD). In our CLuster ANalysis of Sequences (CLANS) analysis, SARS-CoV/SARS-CoV-2 S1-NTDs displayed a close relationship with OC43 and HKU1. However, in the complete and S1-NTD-free spike protein, SARS-CoV/SARS-CoV-2 revealed closeness with MERS-CoV. Further, we have divided the S1-NTD of SARS-CoV-2 related viruses into three distinct types (Type-I to III S1-NTD) and the S1-NTD of viruses associated with SARS-CoVs into another three classes (Type-A to C S1-NTD) using CLANS and phylogenetic analyses. In particular, the results of our study indicate that SARS-CoV-2, RaTG13, and BANAL-20-52 viruses carry Type-I-S1-NTD and other SARS-CoV-2-related-bat viruses have Type-II and III. In addition, it was revealed that the Pangolin-GX and Pangolin-Guangdong lineages inherited Type-I-like and Type-II-like S1-NTD, respectively. Then our CLANS study shows the potential for evolution of Type-I and Type-III S1-NTD from SARS-CoV-related viruses Type-A and Type-B S1-NTDs, respectively. Furthermore, our analysis clarifies the possibility that Type-II S1-NTDs may have evolved from Type-A-S1-NTD of SARS-CoV-related viruses through Type-I S1-NTDs. We also observed that BANAL-20-103, BANAL-20-236, and Pangolin-Guangdong-lineage viruses containing Type-II-like S1-NTD are very close to SARS-CoV-2 in spike genetic areas other than S1-NTD. Possibly, it suggests that the common ancestor spike gene of SARS-CoV-2/RaTG13/BANAL-20-52-like virus may have evolved by recombining the Pangolin-Guangdong/BANAL-20-103/BANAL-20-236-like spike gene to Pangolin-GX-like Type-I-like-S1-NTD in the unsampled bat or undiscovered intermediate host or possibly pangolin. These may then have evolved into SARS-CoV-2, RaTG13, and BANAL-20-52 virus spike genes by host jump mediated evolution. The potential function of the novel Type-I-S1-NTD and other types of S1-NTDs needs to be studied further to understand better its importance in the ongoing COVID-19 outbreak and for future pandemic preparedness.

Epidemiology and Evolution of Emerging Porcine Circovirus-like Viruses in Pigs with Hemorrhagic Dysentery and Diarrhea Symptoms in Central China from 2018 to 2021

Porcine circovirus-like virus (PCLV) is a type of circular Rep-encoding single-stranded DNA virus and may be associated with the development of diarrheal symptoms in pigs. In this study, we retrospectively analyzed three years of past cases in Anhui, China, and reported a case of hemorrhagic enteritis and death in a pregnant sow possibly caused by PCLV. In addition, we analyzed the evolutionary characteristics of PCLV and found that mutation, recombination and selective pressure all played an important role in the evolution of PCLV. We identified N15D and T17S as well as L56T, T58R, K59Q, M62R, L75I and R190K mutations in two different branches, and we noted recombination events in the Rep of a group of Chinese strains. Analysis of selection pressure revealed that PCLV gained more positive selection, indicating that the virus is in a continuous evolutionary state. The PR2 plot, ENC-plot and neutrality analysis showed a greater role of natural selection than that of mutational pressure in the formation of codon usage patterns. This study is the first to identify PCLV in sows with hemorrhagic dysentery and death, and it provides new epidemiological information on PCLV infection in pigs in China.

Evolutionary Patterns of Codon Usage in Major Lineages of Porcine Reproductive and Respiratory Syndrome Virus in China

Porcine reproductive and respiratory syndrome virus (PRRSV) is economically important and characterized by its extensive variation. The codon usage patterns and their influence on viral evolution and host adaptation among different PRRSV strains remain largely unknown. Here, the codon usage of ORF5 genes from lineages 1, 3, 5, and 8, and MLV strains of type 2 PRRSV in China was analyzed. A compositional property analysis of ORF5 genes revealed that nucleotide C is most frequently used at the third position of codons, accompanied by rich GC3s. The effective number of codon (ENC) and codon pair bias (CPB) values indicate that all ORF5 genes have low codon bias and the differences in CPB scores among four lineages are almost not significant. When compared with host codon usage patterns, lineage 1 strains show higher CAI and SiD values, with a high similarity to pig, which might relate to its predominant epidemic propensity in the field. The CAI, RCDI, and SiD values of ORF5 genes from different passages of MLV JXA1R indicate no relation between attenuation and CPB or codon adaptation decrease during serial passage on non-host cells. These findings provide a novel way of understanding the PRRSV's evolution, related to viral survival, host adaptation, and virulence.