Does adaptation to vertebrate codon usage relate to flavivirus emergence potential?

Analysis of synonymous codon usage bias of Lassa virus

Lassa virus genome consists of two single-stranded, negative-sense RNA segments that lie in the genus Arenavirus. The disease associated with the Lassa virus is distributed all over the world, with approximately 3,000,000-5,000,000 infections diagnosed annually in West Africa. It shows high health risks to the human being. Previous research used the evolutionary time scale and adaptive evolution to describe the Lassa virus population pattern. However, it is still unclear how the Lassa virus takes advantage of synonymous codons. In this study, we analyzed the codon usage bias in 162 Lassa virus strains 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 disclosed that LASV strains are rich in A/T. The average ENC value indicated a low codon usage bias in LASVs. The ENC-plot, neutrality plot and parity rule 2 plot demonstrated that, besides mutational pressure, other factors like natural selection also contributed to codon usage bias. This study is significant because it described the pattern of codon usage in the genomes of the Lassa viruses and provided the information needed for a fundamental evolutionary study of them.

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.

Comprehensive analysis of codon usage pattern in Withania somnifera and its associated pathogens: Meloidogyne incognita and Alternaria alternata

Meloidogyne incognita (Root-knot nematode) and Alternaria alternata (fungus) were among the dominant parasites of the medicinal plant Withania somnifera. Despite the fatal nature of their infection, a comprehensive study to explore their evolution and adaptation is lacking. The present study elucidates evolutionary and codon usage bias analysis of W. somnifera (host plant), M. incognita (root-knot nematode) and A. alternata (fungal parasite). The results of the present study revealed a weak codon usage bias prevalent in all the three organisms. Based on the nucleotide analysis, genome of W. somnifera and M. incognita was found to be A-T biased while A. alternata had GC biased genome. We found high similarity of CUB pattern between host and its nematode pathogen as compared to the fungal pathogen. Inclusively, both the evolutionary forces influenced the CUB in host and its associated pathogens. However, neutrality plot indicated the pervasiveness of natural selection on CUB of the host and its pathogens. Correspondence analysis revealed the dominant effect of mutation on CUB of W. somnifera and M. incognita while natural selection was the main force affecting CUB of A. alternata. Taken together the present study would provide some prolific insight into the role of codon usage bias in the adaptability of pathogens to the host’s environment for establishing parasitic relationship.

The Evolution of Severe Acute Respiratory Syndrome Coronavirus-2 during Pandemic and Adaptation to the Host

Severe Acute Respiratory Syndrome Coronavirus-2 is a zoonotic virus with a possible origin in bats and potential transmission to humans through an intermediate host. When zoonotic viruses jump to a new host, they undergo both mutational and natural selective pressures that result in non-synonymous and synonymous adaptive changes, necessary for efficient replication and rapid spread of diseases in new host species. The nucleotide composition and codon usage pattern of SARS-CoV-2 indicate the presence of a highly conserved, gene-specific codon usage bias. The codon usage pattern of SARS-CoV-2 is mostly antagonistic to human and bat codon usage. SARS-CoV-2 codon usage bias is mainly shaped by the natural selection, while mutational pressure plays a minor role. The time-series analysis of SARS-CoV-2 genome indicates that the virus is slowly evolving. Virus isolates from later stages of the outbreak have more biased codon usage and nucleotide composition than virus isolates from early stages of the outbreak.

Mechanisms Underlying Host Range Variation in Flavivirus: From Empirical Knowledge to Predictive Models

Preventing and controlling epidemics caused by vector-borne viruses are particularly challenging due to their diverse pool of hosts and highly adaptive nature. Many vector-borne viruses belong to the Flavivirus genus, whose members vary greatly in host range and specificity. Members of the Flavivirus genus can be categorized to four main groups: insect-specific viruses that are maintained solely in arthropod populations, mosquito-borne viruses and tick-borne viruses that are transmitted to vertebrate hosts by mosquitoes or ticks via blood feeding, and those with no-known vector. The mosquito-borne group encompasses the yellow fever, dengue, and West Nile viruses, all of which are globally spread and cause severe morbidity in humans. The Flavivirus genus is genetically diverse, and its members are subject to different host-specific and vector-specific selective constraints, which do not always align. Thus, understanding the underlying genetic differences that led to the diversity in host range within this genus is an important aspect in deciphering the mechanisms that drive host compatibility and can aid in the constant arms-race against viral threats. Here, we review the phylogenetic relationships between members of the genus, their infection bottlenecks, and phenotypic and genomic differences. We further discuss methods that utilize these differences for prediction of host shifts in flaviviruses and can contribute to viral surveillance efforts.

Guapiaçu virus, a new insect-specific flavivirus isolated from two species of Aedes mosquitoes from Brazil

Classical insect-flaviviruses (cISFVs) and dual host-related insect-specific flavivirus (dISFV) are within the major group of insect-specific flavivirus. Remarkably dISFV are evolutionarily related to some of the pathogenic flavivirus, such as Zika and dengue viruses. The Evolutionary relatedness of dISFV to flavivirus allowed us to investigate the evolutionary principle of host adaptation. Additionally, dISFV can be used for the development of flavivirus vaccines and to explore underlying principles of mammalian pathogenicity. Here we describe the genetic characterization of a novel putative dISFV, termed Guapiaçu virus (GUAPV). Distinct strains of GUAPV were isolated from pools of Aedes terrens and Aedes scapularis mosquitoes. Additionally, we also detected viral GUAPV RNA in a plasma sample of an individual febrile from the Amazon region (North of Brazil). Although GUAPV did not replicate in tested mammalian cells, 3'UTR secondary structures duplication and codon usage index were similar to pathogenic flavivirus.

Mycobacterium lepromatosis genome exhibits unusually high CpG dinucleotide content and selection is key force in shaping codon usage

Mycobacterium lepromatosis was identified as a causative agent for leprosy in the year 2008 in the United States and later more cases were identified in Canada, Singapore, Brazil, and Myanmar. It is known to cause diffuse lepromatosis leprosy among humans. Since it is invasive, the mortality rates are higher in comparison to the M. leprae. At genomic level, there exists 90.9% similarity between M. lepromatosis and M. leprae. Codon usage analysis based on analyses of 228 coding sequences (CDSs) of M. lepromatosis, revealed that the genome is GC rich. Among the total 16 dinucleotides, CpG dinucleotide possesses the highest dinucleotide frequency in M. lepromatosis, that is strikingly an unobvious observation since higher CpG is associated with higher proinflammatory cytokine production and NF-κB activation that eventually leads to high pathogenicity. To evade immune response, CpG content is generally less in pathogens. The unusually high CpG content can be explained by the fact that the nucleotide composition of M. lepromatosis is CG rich. Various forces interplay to shape codon usage pattern of any organism including selection; mutation, nucleotide composition as well as GC biased gene conversion. To understand the interplay between various forces; neutrality, parity, Nc-GC3 (Effective number of codons-GC content at 3rd position of the codon), aromaticity (AROMO) and the general average hydropathicity score (GRAVY) analyses have been carried out. The analyses revealed that selection force is the major contributory force. Along with the selection; mutation, nucleotide composition as well as GC biased gene conversion also play role in shaping codon usage bias in M. lepromatosis. This is the first report on the codon usage in M. lepromatosis.

Genetic Adaptations, Biases, and Evolutionary Analysis of Canine Distemper Virus Asia-4 Lineage in a Fatal Outbreak of Wild-Caught Civets in Thailand

Canine morbillivirus (CDV) is a serious pathogen that can cause fatal systemic disease in a wide range of domestic and wildlife carnivores. Outbreaks of CDV in wildlife species lead to questions regarding the dispersal of the CDV origin. In the present study, we identified a fatal CDV outbreak in caged wild-caught civets in Thailand. Full-length genetic analysis revealed that CDV from the Asia-4 lineage served as the likely causative agent, which was supported by the viral localization in tissues. Evolutionary analysis based on the CDV hemagglutinin (H) gene revealed that the present civet CDV has co-evolved with CDV strains in dogs in Thailand since about 2014. The codon usage pattern of the CDV H gene revealed that the CDV genome has a selective bias of an A/U-ended codon preference. Furthermore, the codon usage pattern of the CDV Asia-4 strain from potential hosts revealed that the usage pattern was related more to the codon usage of civets than of dogs. This finding may indicate the possibility that the discovered CDV had initially adapted its virulence to infect civets. Therefore, the CDV Asia-4 strain might pose a potential risk to civets. Further epidemiological, evolutionary, and codon usage pattern analyses of other CDV-susceptible hosts are required.

Molecular Characterization of Hamster-Adapted Yellow Fever Virus

We previously reported two hamster models for viscerotropic yellow fever virus (YFV) infection: one using a YFV strain (Jiménez), isolated from a fatal human case in Panama in 1974, and the other using the prototype YFV strain (Asibi). Asibi hamster passage 7 (P7) was associated with accumulation of seven amino acid substitutions, including five in the envelope protein. In this study we report the genome sequences of the hamster Jiménez P0 and P10 viruses in which we identified only two amino acid substitutions during passage, one each in the nonstructural proteins NS3 and NS5, indicating a role for the nonstructural proteins in increased YFV viscerotropism in the Jiménez hamster model. Thus, there are multiple molecular mechanisms involved in viscerotropism of YFV in the hamster model. Neither Asibi P7 nor Jiménez P10 viruses were viscerotropic in mice or guinea pigs. Thus, the hamster viscerotropic phenotype did not translate to other laboratory rodent species.

Molecular epidemiological study on Infectious Pancreatic Necrosis Virus isolates from aquafarms in Scotland over three decades

In order to obtain an insight into genomic changes and associated evolution and adaptation of Infectious Pancreatic Necrosis Virus (IPNV), the complete coding genomes of 57 IPNV isolates collected from Scottish aquafarms from 1982 to 2014 were sequenced and analysed. Phylogenetic analysis of the sequenced IPNV strains showed separate clustering of genogroups I, II, III and V. IPNV isolates with genetic reassortment of segment A/B of genogroup III/II were determined. About 59 % of the IPNV isolates belonged to the persistent type and 32 % to the low-virulent type, and only one highly pathogenic strain (1.79 %) was identified. Codon adaptation index calculations indicated that the IPNV major capsid protein VP2 has adapted to its salmonid host. Under-representation of CpG dinucleotides in the IPNV genome to minimize detection by the innate immunity receptors, and observed positive selection in the virulence determination sites of VP2 embedded in the variable region of the main antigenic region, suggest an immune escape mechanism driving virulence evolution. The prevalence of mostly persistent genotypes, together with the assumption of adaptation and immune escape, indicates that IPNV is evolving with the host.

Full-Genome Characterization and Genetic Evolution of West African Isolates of Bagaza Virus

Bagaza virus is a mosquito-borne flavivirus, first isolated in 1966 in Central African Republic. It has currently been identified in mosquito pools collected in the field in West and Central Africa. Emergence in wild birds in Europe and serological evidence in encephalitis patients in India raise questions on its genetic evolution and the diversity of isolates circulating in Africa. To better understand genetic diversity and evolution of Bagaza virus, we describe the full-genome characterization of 11 West African isolates, sampled from 1988 to 2014. Parameters such as genetic distances, N-glycosylation patterns, recombination events, selective pressures, and its codon adaptation to human genes are assessed. Our study is noteworthy for the observation of N-glycosylation and recombination in Bagaza virus and provides insight into its Indian origin from the 13th century. Interestingly, evidence of Bagaza virus codon adaptation to human house-keeping genes is also observed to be higher than those of other flaviviruses well known in human infections. Genetic variations on genome of West African Bagaza virus could play an important role in generating diversity and may promote Bagaza virus adaptation to other vertebrates and become an important threat in human health.