Evolutionary characterization of Tembusu virus infection through identification of codon usage patterns

Stem-Loop I of the Tembusu Virus 3'-Untranslated Region Is Responsible for Viral Host-Specific Adaptation and the Pathogenicity of the Virus in Mice

Tembusu virus (TMUV), an avian mosquito-borne flavivirus, was first identified from Culex tritaeniorhynchus in 1955. To validate the effects of the 3'-untranslated region (3'UTR) in viral host-specific adaptation, we generated a set of chimeric viruses using CQW1 (duck strain) and MM 1775 (mosquito strain) as backbones with heterogeneous 3'UTRs. Compared with rMM 1775, rMM-CQ3'UTR (recombinant MM 1775 virus carrying the 3'UTR of CQW1) exhibited enhanced proliferation in vitro, with peak titers increasing by 5-fold in duck embryonic fibroblast (DEF) cells or 12-fold in baby hamster kidney (BHK-21) cells; however, the neurovirulence of rMM-CQ3'UTR was attenuated in 14-day-old Kunming mice via intracranial injection, with slower weight loss, lower mortality, and reduced viral loads. In contrast, rCQ-MM3'UTR showed similar growth kinetics in vitro and neurovirulence in mice compared with those of rCQW1. Then, the Stem-loop I (SLI) structure, which showed the highest variation within the 3'UTR between CQW1 and MM 1775, was further chosen for making chimeric viruses. The peak titers of rMM-CQ3'UTRSLI displayed a 15- or 4-fold increase in vitro, and the neurovirulence in mice was attenuated, compared with that of rMM 1775; rCQ-MM3'UTRSLI displayed comparable multiplication ability in vitro but was significantly attenuated in mice, in contrast with rCQW1. In conclusion, we demonstrated that the TMUV SLI structure of the 3'UTR was responsible for viral host-specific adaptation of the mosquito-derived strain in DEF and BHK-21 cells and regulated viral pathogenicity in 14-day-old mice, providing a new understanding of the functions of TMUV 3'UTR in viral host switching and the pathogenicity changes in mice. IMPORTANCE Mosquito-borne flaviviruses (MBFVs) constitute a large number of mosquito-transmitted viruses. The 3'-untranslated region (3'UTR) of MBFV has been suggested to be relevant to viral host-specific adaptation. However, the evolutionary strategies for host-specific fitness among MBFV are different, and the virulence-related structures within the 3'UTR are largely unknown. Here, using Tembusu virus (TMUV), an avian MBFV as models, we observed that the duck-derived SLI of the 3'UTR significantly enhanced the proliferation ability of mosquito-derived TMUV in baby hamster kidney (BHK-21) and duck embryonic fibroblast (DEF) cells, suggesting that the SLI structure was crucial for viral host-specific adaptation of mosquito-derived TMUVs in mammalian and avian cells. In addition, all SLI mutant viruses exhibited reduced viral pathogenicity in mice, indicating that SLI structure was a key factor for the pathogenicity in mice. This study provides a new insight into the functions of the MBFV 3'UTR in viral host switching and pathogenicity changes in mice.

Review of antiviral peptides for use against zoonotic and selected non-zoonotic viruses

Viruses remain one of the leading causes of animal and human disease. Some animal viral infections spread sporadically to human populations, posing a serious health risk. Particularly the emerging viral zoonotic diseases such as the novel, zoonotic coronavirus represent an actual challenge for the scientific and medical community. Besides human health risks, some animal viral infections, although still not zoonotic, represent important economic loses to the livestock industry. Viral infections pose a genuine concern for which there has been an increasing interest for new antiviral molecules. Among these novel compounds, antiviral peptides have been proposed as promising therapeutic options, not only for the growing body of evidence showing hopeful results but also due to the many adverse effects of chemical-based drugs. Here we review the current progress, key targets and considerations for the development of antiviral peptides (AVPs). The review summarizes the state of the art of the AVPs tested in zoonotic (coronaviruses, Rift Valley fever viruses, Eastern Equine Encephalitis Virus, Dengue and Junín virus) and also non-zoonotic farm animal viruses (avian and cattle viruses). Their molecular target, amino acid sequence and mechanism of action are summarized and reviewed. Antiviral peptides are currently on the cutting edge since they have been reported to display anti-coronavirus activity. Particularly, the review will discuss the specific mode of action of AVPs that specifically inhibit the fusion of viral and host-cell membranes for SARS-CoV-2, showing in detail some important features of the fusion inhibiting peptides that target the spike protein of these risky viruses.

The spread of Tembusu virus in China from 2010 to 2019

Tembusu virus (TMUV) is a positive-sense RNA virus that is associated with severe reduction in egg production and even death in ducks. TMUV infection shows high incidence and is a threat to the global duck industry. However, the possible origin, genotype, and codon usage bias of TMUV are not very clear. Here, we addressed these questions by analyzing the available genomic sequences from China. The results showed that the ancestor of avian TMUV was most likely a mosquito TMUV. Moreover, three TMUV clades were identified by three different phylogenetic analysis methods. The TMUV genome exhibits a stronger mutation pressure than natural selection pressure. Our findings provide important insights that reveal the ongoing TMUV spread in China and can aid in future prevention and control.

Codon usage bias in the H gene of canine distemper virus

Canine distemper virus (CDV), a non-segmented single negative-stranded RNA (ssRNA), is the etiological agent of canine distemper. Canine distemper is a highly contagious and lethal viral disease in domestic dogs and wild carnivores. Study of the evolution of CDV presents an essential key to improve the vaccine efficacy. In this study, a total of 328 full-length CDV hemagglutinin (H) gene sequences were subjected to phylogenetic, amino acid mutations, and codon usage analysis. In accordance with previous study, CDV genotypes consisted of fifteen lineages. The unique amino acid substitution sites in each CDV lineages have been identified for the first time, including America-1 (Q330H), America-2 (I585S), Asia-1 (A359V), Asia-2 (H61R), Asia-3 (P108Q), Asia-4 (K213T), India-1/Asia-5(S497P), Arctic (S20L), Africa-1(N489S), Colombian (V41I), EWL (I44V), Europe (D560E), Europe-1/South America-1(K161Q), South America-2 (R580Q), and East African (S214A). Codon usage analysis indicated that H gene exhibited low codon usage bias and further neutrality plot analysis demonstrated that natural selection played a dominated role in driving CPV evolution. The effective number of codons (ENC) plots show that all the different sequences are below the standard curve, indicating that mutational pressure is not the only factor affecting CUB but other forces, including natural selection. The neutrality analysis showed that the slope of the regression line was 0.1501, indicating natural selection dominates directional mutation pressure in driving the codon usage pattern. In addition, nucleotide composition, relative synonymous codon usage value, dinucleotide content, and geographical distribution have been proven to influence the codon usage bias of the CDV H gene. The novel findings enhanced the understanding of CDV evolution.

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.

Avian Flavivirus Infection of Monocytes/Macrophages by Extensive Subversion of Host Antiviral Innate Immune Responses

Avian Tembusu virus (TMUV) is a newly emerging avian pathogenic flavivirus in China and Southeast Asia with features of rapid spread, an expanding host range, and cross-species transmission. The mechanisms of its infection and pathogenesis remain largely unclear. Here, we investigated the tropism of this arbovirus in peripheral blood mononuclear cells of specific-pathogen-free (SPF) ducks and SPF chickens and identified monocytes/macrophages as the key targets of TMUV infection. In vivo studies in SPF ducks and SPF chickens with monocyte/macrophage clearance demonstrated that the infection of monocytes/macrophages was crucial for viral replication, transmission, and pathogenesis. Further genome-wide transcriptome analyses of TMUV-infected chicken macrophages revealed that host antiviral innate immune barriers were the major targets of TMUV in macrophages. Despite the activation of major pattern recognition receptor signaling, the inductions of alpha interferon (IFN-α) and IFN-β were blocked by TMUV infection on transcription and translation levels, respectively. Meanwhile, TMUV inhibited host redox responses by repressing the transcription of genes encoding NADPH oxidase subunits and promoting Nrf2-mediated antioxidant responses. The recovery of either of the above-mentioned innate immune barriers was sufficient to suppress TMUV infection. Collectively, we identify an essential step of TMUV infection and reveal extensive subversion of host antiviral innate immune responses.IMPORTANCE Mosquito-borne flaviviruses include a group of pathogenic viruses that cause serious diseases in humans and animals, including dengue, West Nile, and Japanese encephalitis viruses. These flaviviruses are zoonotic and use animals, including birds, as amplifying and reservoir hosts. Avian Tembusu virus (TMUV) is an emerging mosquito-borne flavivirus that is pathogenic for many avian species and can infect cells derived from mammals and humans in vitro Although not currently pathogenic for primates, the infection of duck industry workers and the potential risk of TMUV infection in immunocompromised individuals have been highlighted. Thus, the prevention of TMUV in flocks is important for both avian and mammalian health. Our study reveals the escape of TMUV from the first line of the host defense system in the arthropod-borne transmission route of arboviruses, possibly helping to extend our understanding of flavivirus infection in birds and refine the design of anti-TMUV therapeutics.

Analysis of Synonymous Codon Usage Bias in Flaviviridae Virus

Background

Flaviviridae viruses are single-stranded, positive-sense RNA viruses, which threat human constantly mediated by mosquitoes, ticks, and sandflies. Considering the recent increase in the prevalence of the family virus and its risk potential, we investigated the codon usage pattern to understand its evolutionary processes and provide some useful data to develop the medications for most of Flaviviridae viruses.

Results

The overall extent of codon usage bias in 65 Flaviviridae viruses is low with the average value of GC contents being 50.5% and the highest value being 55.9%; the lowest value is 40.2%. ENC values of Flaviviridae virus genes vary from 48.75 to 57.83 with a mean value of 55.56. U- and A-ended codons are preferred in the Flaviviridae virus. Correlation analysis shows that the positive correlation between ENC value and GC content at the third nucleotide positions was significant in this family virus. The result of analysis of ENC, neutrality plot analysis, and correlation analysis revealed that codon usage bias of all the viruses was affected mainly by natural selection. Meanwhile, according to correspondence analysis (CoA) based on RSCU and phylogenetic analysis, the Flaviviridae viruses mainly are made up of two groups, Group I (Yellow fever virus, Apoi virus, Tembusu virus, Dengue virus 1, and others) and Group II (West Nile virus lineage 2, Japanese encephalitis virus, Usutu virus, Kedougou virus, and others).

Conclusions

All in, the bias of codon usage pattern is affected not only by compositional constraints but also by natural selection. Phylogenetic analysis also illustrates that codon usage bias of virus can serve as an effective means of evolutionary classification in Flaviviridae virus.

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.

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.

An updated review of avian-origin Tembusu virus: a newly emerging avian Flavivirus

Tembusu virus (TMUV, genus Flavivirus, family Flaviviridae) was first isolated in 1955 from Culex tritaeniorhynchus mosquitoes in Kuala Lumpur, Malaysia. In April 2010, duck TMUV was first identified as the causative agent of egg-drop syndrome, characterized by a substantial decrease in egg laying and depression, growth retardation and neurological signs or death in infected egg-laying and breeder ducks, in the People's Republic of China. Since 2010, duck TMUV has spread to most of the duck-producing regions in China, including many of the coastal provinces, neighbouring regions and certain Southeast Asia areas (i.e. Thailand and Malaysia). This review describes the current understanding of the genome characteristics, host range, transmission, epidemiology, phylogenetic and immune evasion of avian-origin TMUV and the innate immune response of the host.

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.

Significant inhibition of Tembusu virus envelope and NS5 gene using an adenovirus-mediated short hairpin RNA delivery system

Tembusu virus (TMUV) is a mosquito-borne flavivirus, which was first isolated in the tropics during the 1970s. Recently, a disease characterized by ovarian haemorrhage and neurological symptoms was observed in ducks in China, which threatens poultry production. However, there is no suitable vaccination strategy or effective antiviral drugs to combat TMUV infections. Consequently, there is an urgent need to develop a new anti-TMUV therapy. In this study, we report an efficient short hairpin RNA (shRNA) delivery strategy for the inhibition of TMUV production using an adenovirus vector system. Using specifically designed shRNAs based on the E and NS5 protein genes of TMUV, the vector-expressed viral genes, TMUV RNA replication and infectious virus production were downregulated at different levels in Vero cells, where the shRNA (NS52) was highly effective in inhibiting TMUV. Using the human adenovirus type 5 shRNA delivery system, the recombinant adenovirus (rAd-NS52) inhibited TMUV multiplication with high efficiency. Furthermore, the significant dose-dependent inhibition of viral RNA copies induced by rAd-NS52 was found in TMUV-infected cells, which could last for at least 96h post infection. Our results indicated that the adenovirus-mediated delivery of shRNAs could play an active role in future TMUV antiviral therapeutics.

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.

Antigen distribution of TMUV and GPV are coincident with the expression profiles of CD8α-positive cells and goose IFNγ

Both Tembusu virus (TMUV) and goose parvovirus (GPV) are causative agents of goose disease. However, the host immune response of the goose against these two different categories of virus has not been well documented. Here, we compared the clinical symptoms and pathological characteristics, antigen distribution and intensity, and expression of immune-related genes in TMUV- and GPV- infected goose. The immunohistochemistry analysis demonstrated that GPV was primarily located in the liver, lung, small intestine, and rectum, while TMUV was situated in the liver, brain, spleen, and small intestine. The induction of IFNγ and proinflammatory cytokines is highly associated with the distribution profiles of antigen and CD8α+ molecules. The effector function of CD8 T cells may be accomplished by the secretion of IFNγ together with high expression of proinflammatory cytokines such as IL1 and IL6. Remarkably, significant increases in the transcription of immune genes were observed after infection, which suggested that both GPV and TMUV can effectively induce immune response in goose PMBCs. This study will provide fundamental information for goose molecular immunology in defending against pandemic viruses.