G+C content differs in conserved and variable amino acid residues of flaviviruses and other evolutionary groups

Comparative analysis of human coronaviruses focusing on nucleotide variability and synonymous codon usage patterns

The prevailing COVID-19 pandemic has drawn the attention of the scientific community to study the evolutionary origin of Severe Acute Respiratory Syndrome Corona Virus 2 (SARS-CoV-2). This study is a comprehensive quantitative analysis of the protein-coding sequences of seven human coronaviruses (HCoVs) to decipher the nucleotide sequence variability and codon usage patterns. It is essential to understand the survival ability of the viruses, their adaptation to hosts, and their evolution. The current analysis revealed a high abundance of the relative dinucleotide (odds ratio), GC and CT pairs in the first and last two codon positions, respectively, as well as a low abundance of the CG pair in the last two positions of the codon, which might be related to the evolution of the viruses. A remarkable level of variability of GC content in the third position of the codon among the seven coronaviruses was observed. Codons with high RSCU values are primarily from the aliphatic and hydroxyl amino acid groups, and codons with low RSCU values belong to the aliphatic, cyclic, positively charged, and sulfur-containing amino acid groups. In order to elucidate the evolutionary processes of the seven coronaviruses, a phylogenetic tree (dendrogram) was constructed based on the RSCU scores of the codons. The severe and mild categories CoVs were positioned in different clades. A comparative phylogenetic study with other coronaviruses depicted that SARS-CoV-2 is close to the CoV isolated from pangolins (Manis javanica, Pangolin-CoV) and cats (Felis catus, SARS(r)-CoV). Further analysis of the effective number of codon (ENC) usage bias showed a relatively higher bias for SARS-CoV and MERS-CoV compared to SARS-CoV-2. The ENC plot against GC3 suggested that the mutational bias might have a role in determining the codon usage variation among candidate viruses. A codon adaptability study on a few human host parasites (from different kingdoms), including CoVs, showed a diverse adaptability pattern. SARS-CoV-2 and SARS-CoV exhibit relatively lower but similar codon adaptability compared to MERS-CoV.

Flavivirus integrations in Aedes aegypti are limited and highly conserved across samples from different geographic regions unlike integrations in Aedes albopictus

Mosquitoes of the genus Aedes are the main vectors of many viruses, e.g. dengue and Zika, which affect millions of people each year and for which there are limited treatment options. Understanding how Aedes mosquitoes tolerate high viral loads may lead to better disease control strategies. Elucidating endogenous viral elements (EVEs) within vector genomes may give exploitable biological insights. Previous studies have reported the presence of a large number of EVEs in Aedes genomes. Here we investigated if flavivirus EVEs are conserved across populations and different Aedes species by using ~ 500 whole genome sequence libraries from Aedes aegypti and Aedes albopictus, sourced from colonies and field mosquitoes across continents. We found that nearly all flavivirus EVEs in the Ae. aegypti reference genome originate from four separate putative viral integration events, and that they are highly conserved across geographically diverse samples. By contrast, flavivirus EVEs in the Ae. albopictus reference genome originate from up to nine distinct integration events and show low levels of conservation, even within samples from narrow geographical ranges. Our analysis suggests that flaviviruses integrated as long sequences and were subsequently fragmented and shuffled by transposable elements. Given that EVEs of Ae. aegypti and Ae. albopictus belong to different phylogenetic clades and have very differing levels of conservation, they may have different evolutionary origins and potentially different functional roles.

Molecular Basis for Arbovirus Transmission by Aedes aegypti Mosquitoes

Mosquitoes are considered the most important vectors for the transmission of pathogens to humans. Aedes aegypti is a unique species, not only by its highly anthropophilic and peridomestic habits but also because it can transmit an important variety of pathogenic viruses. Examples are dengue, yellow fever, chikungunya, Zika, and Mayaro viruses. After ingesting viremic blood, a wide range of mechanisms are activated in the mosquito to counteract viral infection. Nevertheless, these arboviruses possess strategies to overcome barriers in the mosquito and eventually reach the salivary glands to continue the transmission cycle. However, the infection and eventual transmission of arbovirus depends on multiple factors. The current review focuses in detail on the anatomic, physiological, and molecular characteristics of the mosquito A. aegypti that participate in response to a viral infection. In the past decades, the awareness of the importance of this mosquito as a disease vector and its impact on human health was largely recognized. We need to improve our comprehension of molecular mechanisms that determine the outcome of successful virus replication or control of infection for each arbovirus in the vector; this could lead to the design of effective control strategies in the future.

Computational screening of medicinal plant phytochemicals to discover potent pan-serotype inhibitors against dengue virus

Emergence of Dengue as one of the deadliest viral diseases prompts the need for development of effective therapeutic agents. Dengue virus (DV) exists in four different serotypes and infection caused by one serotype predisposes its host to another DV serotype heterotypic re-infection. We undertook virtual ligand screening (VLS) to filter compounds against DV that may inhibit inclusively all of its serotypes. Conserved non-structural DV protein targets such as NS1, NS3/NS2B and NS5, which play crucial role in viral replication, infection cycle and host interaction, were selected for screening of vital antiviral drug leads. A dataset of plant based natural antiviral derivatives was developed. Molecular docking was performed to estimate the spatial affinity of target compounds for the active sites of DV’s NS1, NS3/NS2B and NS5 proteins. The drug likeliness of the screened compounds was followed by ADMET analysis whereas the binding behaviors were further elucidated through molecular dynamics (MD) simulation experiments. VLS screened three potential compounds including Canthin-6-one 9-O-beta-glucopyranoside, Kushenol W and Kushenol K which exhibited optimal binding with all the three conserved DV proteins. This study brings forth novel scaffolds against DV serotypes to serve as lead molecules for further optimization and drug development against all DV serotypes with equal effect against multiple disease causing DV proteins. We therefore anticipate that the insights given in the current study could be regarded valuable towards exploration and development of a broad-spectrum natural anti-dengue therapy.

Codon usage variation of Zika virus: The potential roles of NS2B and NS4A in its global pandemic

A comprehensive demonstration of Zika virus (ZIKV) molecular evolution is essential for understanding its adaptation and expansion in its recent pandemics. Despite several studies on mutations and codon usage in ZIKVs, the variations in codon usage patterns across individual genes and their biological implication remains unclear. Here, we performed a gene-by-gene comparison of the codon usage variation in ZIKVs of the African and Asian lineages. We found that besides the evidence of positive selection (Ka/Ks >1) in the Asian lineage of the ZIKV genome, codon usage patterns were gene-specific and codon usage variation of ZIKV genes, was possibly constrained by their individual functional features, such as transmembrane domains, or antigenicity. In particular, the NS2B and NS4A genes showed distinct codon usage patterns, clearly separating them from the clusters of other genes in the correspondence analysis (CA). In the Asian lineage, the NS2B and NS4A genes showed the highest codon usage bias (ENC values: 51.01 ± 0.72 and 48.89 ± 0.99 respectively), and were subjected to the highest translation selection (ENC_obs/ENC_exp ratio: 0.847 ± 0.0297 and 0.828 ± 0.0233 respectively) in comparison to the African lineages of ZIKV. The CpG frequency of the NS2B showed a gradual ascending trend in the Asian ZIKV lineages, while in NS4A it was constrained along with the expansion of the Asian lineage. Furthermore, between the African and Asian lineages, differentiated and specific over-represented codons were more prominent in the NS2B and NS4A. Together, our study implies that ZIKVs are in the process of evolutionary fine tuning their codon as seen in the recent pandemics, and NS2B and NS4A could have played a potential role in the molecular evolution of the Asian lineage and their establishment.