A comprehensive analysis of genome composition and codon usage patterns of emerging coronaviruses

Optimization and Deoptimization of Codons in SARS-CoV-2 and Related Implications for Vaccine Development

The spread of coronavirus disease 2019 (COVID-19), caused by severe respiratory syndrome coronavirus 2 (SARS-CoV-2), has progressed into a global pandemic. To date, thousands of genetic variants have been identified among SARS-CoV-2 isolates collected from patients. Sequence analysis reveals that the codon adaptation index (CAI) values of viral sequences have decreased over time but with occasional fluctuations. Through evolution modeling, it is found that this phenomenon may result from the virus's mutation preference during transmission. Using dual-luciferase assays, it is further discovered that the deoptimization of codons in the viral sequence may weaken protein expression during virus evolution, indicating that codon usage may play an important role in virus fitness. Finally, given the importance of codon usage in protein expression and particularly for mRNA vaccines, it is designed several codon-optimized Omicron BA.2.12.1, BA.4/5, and XBB.1.5 spike mRNA vaccine candidates and experimentally validated their high levels of expression. This study highlights the importance of codon usage in virus evolution and provides guidelines for codon optimization in mRNA and DNA vaccine development.

Analysis of 3.5 million SARS-CoV-2 sequences reveals unique mutational trends with consistent nucleotide and codon frequencies

BACKGROUND

Since the onset of the SARS-CoV-2 pandemic, bioinformatic analyses have been performed to understand the nucleotide and synonymous codon usage features and mutational patterns of the virus. However, comparatively few have attempted to perform such analyses on a considerably large cohort of viral genomes while organizing the plethora of available sequence data for a month-by-month analysis to observe changes over time. Here, we aimed to perform sequence composition and mutation analysis of SARS-CoV-2, separating sequences by gene, clade, and timepoints, and contrast the mutational profile of SARS-CoV-2 to other comparable RNA viruses.

METHODS

Using a cleaned, filtered, and pre-aligned dataset of over 3.5 million sequences downloaded from the GISAID database, we computed nucleotide and codon usage statistics, including calculation of relative synonymous codon usage values. We then calculated codon adaptation index (CAI) changes and a nonsynonymous/synonymous mutation ratio (dN/dS) over time for our dataset. Finally, we compiled information on the types of mutations occurring for SARS-CoV-2 and other comparable RNA viruses, and generated heatmaps showing codon and nucleotide composition at high entropy positions along the Spike sequence.

RESULTS

We show that nucleotide and codon usage metrics remain relatively consistent over the 32-month span, though there are significant differences between clades within each gene at various timepoints. CAI and dN/dS values vary substantially between different timepoints and different genes, with Spike gene on average showing both the highest CAI and dN/dS values. Mutational analysis showed that SARS-CoV-2 Spike has a higher proportion of nonsynonymous mutations than analogous genes in other RNA viruses, with nonsynonymous mutations outnumbering synonymous ones by up to 20:1. However, at several specific positions, synonymous mutations were overwhelmingly predominant.

CONCLUSIONS

Our multifaceted analysis covering both the composition and mutation signature of SARS-CoV-2 gives valuable insight into the nucleotide frequency and codon usage heterogeneity of SARS-CoV-2 over time, and its unique mutational profile compared to other RNA viruses.

Componential usage patterns in dengue 4 viruses reveal their better evolutionary adaptation to humans

There have been at least four types of dengue outbreaks in the past few years. The evolutionary characteristics of dengue viruses have aroused great concerns. The evolutionary characteristics of dengue 4 viruses are studied in the present study based on their base usage patterns and codon usage patterns. The effective number of codons and relative synonymous codon usage (RSCU) values of four types of dengue viruses were counted or calculated. The Kullback–Leibler (K–L) divergences of relative synonymous codon usage from dengue viruses to humans and the Kullback–Leibler divergences of amino acid usage patterns from dengue viruses to humans were calculated to explore the adaptation levels of dengue viruses. The results suggested that: (1) codon adaptation in dengue 4 viruses occurred through an evolutionary process from 1956 to 2021, (2) overall relative synonymous codon usage values of dengue 4 viruses showed more similarities to humans than those of other subtypes of dengue viruses, and (3) the smaller Kullback–Leibler divergence of amino acid usage and relative synonymous codon usage from dengue viruses to humans indicated that the dengue 4 viruses adapted to human hosts better. All results indicated that both mutation pressure and natural selection pressure contributed to the codon usage pattern of dengue 4 viruses more obvious than to other subtypes of dengue viruses and that the dengue 4 viruses adapted to human hosts better than other types of dengue viruses during their evolutionary process.

Covariance predicts conserved protein residue interactions important for the emergence and continued evolution of SARS-CoV-2 as a human pathogen

SARS-CoV-2 is one of three recognized coronaviruses (CoVs) that have caused epidemics or pandemics in the 21st century and that likely emerged from animal reservoirs. Differences in nucleotide and protein sequence composition within related β-coronaviruses are often used to better understand CoV evolution, host adaptation, and their emergence as human pathogens. Here we report the comprehensive analysis of amino acid residue changes that have occurred in lineage B β-coronaviruses that show covariance with each other. This analysis revealed patterns of covariance within conserved viral proteins that potentially define conserved interactions within and between core proteins encoded by SARS-CoV-2 related β-coronaviruses. We identified not only individual pairs but also networks of amino acid residues that exhibited statistically high frequencies of covariance with each other using an independent pair model followed by a tandem model approach. Using 149 different CoV genomes that vary in their relatedness, we identified networks of unique combinations of alleles that can be incrementally traced genome by genome within different phylogenic lineages. Remarkably, covariant residues and their respective regions most abundantly represented are implicated in the emergence of SARS-CoV-2 and are also enriched in dominant SARS-CoV-2 variants.

Codon usage of host-specific P genotypes (VP4) in group A rotavirus

Background

Group A rotavirus (RVA) is a common causative agent of acute gastroenteritis in infants and young children worldwide. RVA P genotypes, determined by VP4 sequences, have been confirmed to infect humans and animals. However, their codon usage patterns that are essential to obtain insights into the viral evolution, host adaptability, and genetic characterization remained unclear, especially across animal hosts.

Results

We performed a comprehensive codon usage analysis of eight host-specific RVA P genotypes, including human RVA (P[4] and P[8]), porcine RVA (P[13] and P[23]), and zoonotic RVA (P[1], P[6], P[7] and P[19]), based on 233 VP4 complete coding sequences. Nucleotide composition, relative synonymous codon usage (RSCU), and effective number of codons (ENC) were calculated. Principal component analysis (PCA) based on RSCU values was used to explore the codon usage patterns of different RVA P genotypes. In addition, mutation pressure and natural selection were identified by using ENC-plot, parity rule 2 plot, and neutrality plot analyses. All VP4 sequences preferred using A/U nucleotides (A: 0.354-0.377, U: 0.267-0.314) than G/C nucleotides across genotypes. Similarly, majority of commonly used synonymous codons were likely to end with A/U nucleotides (A: 9/18-12/18, U: 6/18-9/18). In PCA, human, porcine, and zoonotic genotypes clustered separately in terms of RSCU values, indicating the host-specific codon usage patterns; however, porcine and zoonotic genotypes were partly overlapped. Human genotypes, P[4] and P[8], had stronger codon usage bias, as indicated by more over-represented codons and lower ENC, compared to porcine and zoonotic genotypes. Moreover, natural selection was determined to be a predominant driver in shaping the codon usage bias across the eight P genotypes. In addition, mutation pressure contributed to the codon usage bias of human genotypes.

Conclusions

Our study identified a strong codon usage bias of human RVA P genotypes attributable to both natural selection and mutation pressure, whereas similar codon usage bias between porcine and zoonotic genotypes predominantly attributable to natural selection. It further suggests possible cross-species transmission. Therefore, it warrants further surveillance of RVA P genotypes for early identification of zoonotic infection.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08730-2.

Natural selection plays a significant role in governing the codon usage bias in the novel SARS-CoV-2 variants of concern (VOC)

The ongoing prevailing COVID-19 pandemic caused by SARS-CoV-2 is becoming one of the major global health concerns worldwide. The SARS-CoV-2 genome encodes spike (S) glycoprotein that plays a very crucial role in viral entry into the host cell via binding of its receptor binding domain (RBD) to the host angiotensin converting enzyme 2 (ACE2) receptor. The continuously evolving SARS-CoV-2 genome results in more severe and transmissible variants characterized by the emergence of novel mutations called 'variants of concern' (VOC). The currently designated alpha, beta, gamma, delta and omicron VOC are the focus of this study due to their high transmissibility, increased virulence, and concerns for decreased effectiveness of the available vaccines. In VOC, the spike (S) gene and other non-structural protein mutations may affect the efficacies of the approved COVID-19 vaccines. To understand the diversity of SARS-CoV-2, several studies have been performed on a limited number of sequences. However, only a few studies have focused on codon usage bias (CUBs) pattern analysis of all the VOC strains. Therefore, to evaluate the evolutionary divergence of all VOC S-genes, we performed CUBs analysis on 300,354 sequences to understand the evolutionary relationship with its adaptation in different hosts, i.e., humans, bats, and pangolins. Base composition and RSCU analysis revealed the presence of 20 preferred AU-ended and 10 under-preferred GC-ended codons. In addition, CpG was found to be depleted, which may be attributable to the adaptive response by viruses to escape from the host defense process. Moreover, the ENC values revealed a higher bias in codon usage in the VOC S-gene. Further, the neutrality plot analysis demonstrated that S-genes analyzed in this study are under 83.93% influence of natural selection, suggesting its pivotal role in shaping the CUBs. The CUBs pattern of S-genes was found to be very similar among all the VOC strains. Interestingly, we observed that VOC strains followed a trend of antagonistic codon usage with respect to the human host. The identified CUBs divergence would help to understand the virus evolution and its host adaptation, thus help design novel vaccine strategies against the emerging VOC strains. To the best of our knowledge, this is the first report for identifying the evolution of CUBs pattern in all the currently identified VOC.

Variation in synonymous nucleotide composition among genomes of sarbecoviruses and consequences for the origin of COVID-19

The subgenus Sarbecovirus includes two human viruses, SARS-CoV and SARS-CoV-2, respectively responsible for the SARS epidemic and COVID-19 pandemic, as well as many bat viruses and two pangolin viruses. Here, the synonymous nucleotide composition (SNC) of Sarbecovirus genomes was analysed by examining third codon-positions, dinucleotides, and degenerate codons. The results show evidence for the eight following groups: (i) SARS-CoV related coronaviruses (SCoVrC including many bat viruses from China), (ii) SARS-CoV-2 related coronaviruses (SCoV2rC; including five bat viruses from Cambodia, Thailand and Yunnan), (iii) pangolin sarbecoviruses, (iv) three bat sarbecoviruses showing evidence of recombination between SCoVrC and SCoV2rC genomes, (v) two highly divergent bat sarbecoviruses from Yunnan, (vi) the bat sarbecovirus from Japan, (vii) the bat sarbecovirus from Bulgaria, and (viii) the bat sarbecovirus from Kenya. All these groups can be diagnosed by specific nucleotide compositional features except the one concerned by recombination between SCoVrC and SCoV2rC. In particular, SCoV2rC genomes have less cytosines and more uracils at third codon-positions than other sarbecoviruses, whereas the genomes of pangolin sarbecoviruses show more adenines at third codon-positions. I suggest that taxonomic differences in the imbalanced nucleotide pools available in host cells during viral replication can explain the eight groups of SNC here detected among Sarbecovirus genomes. A related effect due to hibernating bats and their latitudinal distribution is also discussed. I conclude that the two independent host switches from Rhinolophus bats to pangolins resulted in convergent mutational constraints and that SARS-CoV-2 emerged directly from a horseshoe bat sarbecovirus.

Analysis of SARS-CoV-2 synonymous codon usage evolution throughout the COVID-19 pandemic

SARS-CoV-2, the seventh coronavirus known to infect humans, can cause severe life-threatening respiratory pathologies. To better understand SARS-CoV-2 evolution, genome-wide analyses have been made, including the general characterization of its codons usage profile. Here we present a bioinformatic analysis of the evolution of SARS-CoV-2 codon usage over time using complete genomes collected since December 2019. Our results show that SARS-CoV-2 codon usage pattern is antagonistic to, and it is getting farther away from that of the human host. Further, a selection of deoptimized codons over time, which was accompanied by a decrease in both the codon adaptation index and the effective number of codons, was observed. All together, these findings suggest that SARS-CoV-2 could be evolving, at least from the perspective of the synonymous codon usage, to become less pathogenic.

Covariance predicts conserved protein residue interactions important to the emergence and continued evolution of SARS-CoV-2 as a human pathogen

SARS-CoV-2 is one of three recognized coronaviruses (CoVs) that have caused epidemics or pandemics in the 21st century and that likely emerged from animal reservoirs. Differences in nucleotide and protein sequence composition within related β-coronaviruses are often used to better understand CoV evolution, host adaptation, and their emergence as human pathogens. Here we report the comprehensive analysis of amino acid residue changes that have occurred in lineage B β-coronaviruses that show covariance with each other. This analysis revealed patterns of covariance within conserved viral proteins that potentially define conserved interactions within and between core proteins encoded by SARS-CoV-2 related β-coranaviruses. We identified not only individual pairs but also networks of amino acid residues that exhibited statistically high frequencies of covariance with each other using an independent pair model followed by a tandem model approach. Using 149 different CoV genomes that vary in their relatedness, we identified networks of unique combinations of alleles that can be incrementally traced genome by genome within different phylogenic lineages. Remarkably, covariant residues and their respective regions most abundantly represented are implicated in the emergence of SARS-CoV-2 are also enriched in dominant SARS-CoV-2 variants.

A detailed comparative analysis of codon usage bias in Alongshan virus

Alongshan virus (ALSV) is an emerging tick-borne pathogen that infects humans, causing febrile disease. ALSV uses Ixodes Persulcatus ticks to infect humans with a wide range of signs, from asymptomatic to encephalitis-like syndrome. There is an increasing public health concern about the ALSV infection. To get insight into the impacts of viral relations with their hosts on viral ability, survival, and evasion from hosts immune systems remain unknown. The codon usage is a driving force in viral genome evolution; therefore, we enrolled 41 ALSV strains in codon usage analysis to elucidate the molecular evolutionary dynamics of ALSV. The results indicate that the overall codon usage among ALSV isolates is relatively similar and slightly biased. Base compositions for the cds were in order of G >A >C >U and in the third position of codons G3 >A3 >C3 >T3. The RSCU values revealed that the more frequently used codons were mostly GC ended. Different codon preferences in ALSV genes in relation to codon usage of H. sapiens and Ixodes Persulcatus genes were found. Neutrality plot was determined to reveal the superiority of natural selection over directional mutation pressure in causing CUB based on GC12 versus GC3 contents. The results of these studies suggest that the emergence of ALSV in China, Russia and Finland may also be reflected in ALSV codon usage. Altogether, the presence of both mutation pressure and natural selection effect in shaping the codon usage patterns of ALSV.

Bilirubin levels as an independent predictor of myocarditis in patients with COVID-19

Background

Myocardial damage worsens the clinical course and prognosis of coronavirus disease 2019 (COVID-19) patients. High total bilirubin levels have been associated with a poor prognosis in COVID-19. This study aimed to investigate the predictive value of the total bilirubin level, a marker of heme oxygenase-1 enzyme activity, in determining myocarditis in patients with COVID-19.

Results

A total of 190 patients diagnosed with COVID-19 were enrolled in the study. The patients were divided into two groups based on their troponin positivity. The study group (n = 95) consisted of patients with high troponin, and the control group (n = 95) consisted of patients without high troponin levels. The D-dimer (727 [572–995] vs. 591 [440–790], p = 0.001), C-reactive protein (CRP) (30.0 [10–48] vs. 10.3 [5.8–15.9], p < 0.001), and total bilirubin (9.5 [8.2–12.1] vs. 7.0 [5.3–8.0], p < 0.001) levels were significantly higher in the study group. In multivariate analysis, CRP (odds ratio [OR]: 1.103; 95% confidence interval [CI]: 1.060–1.148; p < 0.001) and total bilirubin (OR: 1.612; 95% CI: 1.330–1.954; p < 0.001) levels were independent predictors of myocarditis in COVID-19.

Conclusions

Total bilirubin levels can be used as an early predictor of myocarditis in COVID-19 and can contribute to therapy management.

Review: Development of SARS-CoV-2 immuno-enhanced COVID-19 vaccines with nano-platform

Vaccination is the most effective way to prevent coronavirus disease 2019 (COVID-19). Vaccine development approaches consist of viral vector vaccines, DNA vaccine, RNA vaccine, live attenuated virus, and recombinant proteins, which elicit a specific immune response. The use of nanoparticles displaying antigen is one of the alternative approaches to conventional vaccines. This is due to the fact that nano-based vaccines are stable, able to target, form images, and offer an opportunity to enhance the immune responses. The diameters of ultrafine nanoparticles are in the range of 1-100 nm. The application of nanotechnology on vaccine design provides precise fabrication of nanomaterials with desirable properties and ability to eliminate undesirable features. To be successful, nanomaterials must be uptaken into the cell, especially into the target and able to modulate cellular functions at the subcellular levels. The advantages of nano-based vaccines are the ability to protect a cargo such as RNA, DNA, protein, or synthesis substance and have enhanced stability in a broad range of pH, ambient temperatures, and humidity for long-term storage. Moreover, nano-based vaccines can be engineered to overcome biological barriers such as nonspecific distribution in order to elicit functions in antigen presenting cells. In this review, we will summarize on the developing COVID-19 vaccine strategies and how the nanotechnology can enhance antigen presentation and strong immunogenicity using advanced technology in nanocarrier to deliver antigens. The discussion about their safe, effective, and affordable vaccines to immunize against COVID-19 will be highlighted.

Variability in Codon Usage in Coronaviruses Is Mainly Driven by Mutational Bias and Selective Constraints on CpG Dinucleotide

The Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the third human-emerged virus of the 21st century from the Coronaviridae family, causing the ongoing coronavirus disease 2019 (COVID-19) pandemic. Due to the high zoonotic potential of coronaviruses, it is critical to unravel their evolutionary history of host species breadth, host-switch potential, adaptation and emergence, to identify viruses posing a pandemic risk in humans. We present here a comprehensive analysis of the composition and codon usage bias of the 82 Orthocoronavirinae members, infecting 47 different avian and mammalian hosts. Our results clearly establish that synonymous codon usage varies widely among viruses, is only weakly dependent on their primary host, and is dominated by mutational bias towards AU-enrichment and by CpG avoidance. Indeed, variation in GC3 explains around 34%, while variation in CpG frequency explains around 14% of total variation in codon usage bias. Further insight on the mutational equilibrium within Orthocoronavirinae revealed that most coronavirus genomes are close to their neutral equilibrium, the exception being the three recently infecting human coronaviruses, which lie further away from the mutational equilibrium than their endemic human coronavirus counterparts. Finally, our results suggest that, while replicating in humans, SARS-CoV-2 is slowly becoming AU-richer, likely until attaining a new mutational equilibrium.

Relative synonymous codon usage of ORF1ab in SARS-CoV-2 and SARS-CoV

Background

COVID-19, as a novel coronavirus disease caused by new coronavirus SARS-CoV-2, spreads all over the world, and brings harm to human in many countries. Humans suffered a lot from both SARS-CoV-2 now and by SARS-CoV in the year 2003. It is important to understand the differences and the relationships between these two types of viruses.

Objective

To compare relative synonymous codon usage of ORF1ab gene in SARS-CoV-2 and SARS-CoV, relative synonymous codon usage of their genomes are studied in this paper from the bioinformatics perspective.

Methods

The ORF1ab gene, which is an important non-structural polyprotein coding gene and now used for nucleic acid detection markers in many measurement method, in both SARS-CoV-2 (30 strains) and SARS-CoV (20 strains) are considered to be the research object in the present paper. The relative synonymous codon usage values of the ORF1ab gene are calculated to characterize the differences and the evolutionary characteristics among 50 strains.

Results

There is a significant difference between SARS-CoV and SARS-CoV-2 when the relative synonymous codon usage value of ORF1ab genes is concerned. The results suggest that codon usage pattern of SARS-CoV is more similar to human than that of the SARS-CoV-2, and that the inner difference in SARS-CoV-2 strains is larger than that of SARS-CoV, which denote the larger diversity exits in the SARS-CoV-2 virus.

Conclusion

These results show that the relative synonymous codon usage values in the coronavirus could be used for further research on their evolutionary phenomenon.

Edging on Mutational Bias, Induced Natural Selection From Host and Natural Reservoirs Predominates Codon Usage Evolution in Hantaan Virus

The molecular evolutionary dynamics that shape hantaviruses’ evolution are poorly understood even now, besides the contribution of virus-host interaction to their evolution remains an open question. Our study aimed to investigate these two aspects in Hantaan virus (HTNV)—the prototype of hantaviruses and an emerging zoonotic pathogen that infects humans, causing hemorrhagic fever with renal syndrome (HFRS): endemic in Far East Russia, China, and South Korea—via a comprehensive, phylogenetic-dependent codon usage analysis. We found that host- and natural reservoir-induced natural selection is the primary determinant of its biased codon choices, exceeding the mutational bias effect. The phylogenetic analysis of HTNV strains resulted in three distinct clades: South Korean, Russian, and Chinese. An effective number of codon (ENC) analysis showed a slightly biased codon usage in HTNV genomes. Nucleotide composition and RSCU analyses revealed a significant bias toward A/U nucleotides and A/U-ended codons, indicating the potential influence of mutational bias on the codon usage patterns of HTNV. Via ENC-plot, Parity Rule 2 (PR2), and neutrality plot analyses, we would conclude the presence of both mutation pressure and natural selection effect in shaping the codon usage patterns of HTNV; however, natural selection is the dominant factor influencing its codon usage bias. Codon adaptation index (CAI), Relative codon deoptimization index (RCDI), and Similarity Index (SiD) analyses uncovered the intense selection pressure from the host (Human) and natural reservoirs (Striped field mouse and Chinese white-bellied rat) in shaping HTNV biased codon choices. Our study clearly revealed the evolutionary processes in HTNV and the role of virus-host interaction in its evolution. Moreover, it opens the door for a more comprehensive codon usage analysis for all hantaviruses species to determine their molecular evolutionary dynamics and adaptability to several hosts and environments. We believe that our research will help in a better and deep understanding of HTNV evolution that will serve its future basic research and aid live attenuated vaccines design.

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.

Genome-Wide Analysis of Codon Usage Patterns of SARS-CoV-2 Virus Reveals Global Heterogeneity of COVID-19

The ongoing outbreak of coronavirus disease COVID-19 is significantly implicated by global heterogeneity in the genome organization of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The causative agents of global heterogeneity in the whole genome of SARS-CoV-2 are not well characterized due to the lack of comparative study of a large enough sample size from around the globe to reduce the standard deviation to the acceptable margin of error. To better understand the SARS-CoV-2 genome architecture, we have performed a comprehensive analysis of codon usage bias of sixty (60) strains to get a snapshot of its global heterogeneity. Our study shows a relatively low codon usage bias in the SARS-CoV-2 viral genome globally, with nearly all the over-preferred codons' A.U. ended. We concluded that the SARS-CoV-2 genome is primarily shaped by mutation pressure; however, marginal selection pressure cannot be overlooked. Within the A/U rich virus genomes of SARS-CoV-2, the standard deviation in G.C. (42.91% ± 5.84%) and the GC3 value (30.14% ± 6.93%) points towards global heterogeneity of the virus. Several SARS-CoV-2 viral strains were originated from different viral lineages at the exact geographic location also supports this fact. Taking all together, these findings suggest that the general root ancestry of the global genomes are different with different genome's level adaptation to host. This research may provide new insights into the codon patterns, host adaptation, and global heterogeneity of SARS-CoV-2.

Genetic and evolutionary analysis of emerging HoBi-like pestivirus

HoBi-like pestivirus, an emerging species within the Pestivirus genus, is an important pathogen associated with a variety of clinical manifestations of ruminants, especially cattle. HoBi-like pestiviruses were identified in several countries and from different hosts, and raised concerns with regard to their acute and persistent infections, which is implicated in economic losses for cattle farmers. However, the transmission path, codon usage bias, and host adaptation of the virus has not been studied. Hence, we performed the analysis the spatio-temporal transmission based on the available 5'-UTR sequences of HoBi-like pestivirus, and then conducted codon analysis of the complete coding sequence of the virus. The results show the virus appeared in 1952 (95% HPD: 1905-1985) and may have originated in India. In addition, Italy is the hub for the spread of the virus. Moreover, six potential recombination events and two complex recombination events were discovered. Analysis of codon usage patterns revealed that the effective number of codon (ENC) values with an average of 50.85, and the codon usage bias is greatly affected by natural selection, which is different from the previous BVDV-1, 2. Finally, codon adaptation index (CAI) analysis shows that pigs may be the potential origin species of the HoBi-like pestivirus. These findings will contribute to more effective control of the spread of the virus, extend the knowledge about the genetic and evolutionary features of HoBi-like viruses and provide some information for vaccine research.

Base Composition and Host Adaptation of the SARS-CoV-2: Insight From the Codon Usage Perspective

The novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been spreading rapidly all over the world and has raised grave concern globally. The present research aims to conduct a robust base compositional analysis of SARS-CoV-2 to reveal adaptive intricacies to the human host. Multivariate statistical analysis revealed a complex interplay of various factors including compositional constraint, natural selection, length of viral coding sequences, hydropathicity, and aromaticity of the viral gene products that are operational to codon usage patterns, with compositional bias being the most crucial determinant. UpG and CpA dinucleotides were found to be highly preferred whereas, CpG dinucleotide was mostly avoided in SARS-CoV-2, a pattern consistent with the human host. Strict avoidance of the CpG dinucleotide might be attributed to a strategy for evading a human immune response. A lower degree of adaptation of SARS-CoV-2 to the human host, compared to Middle East respiratory syndrome (MERS) coronavirus and SARS-CoV, might be indicative of its milder clinical severity and progression contrasted to SARS and MERS. Similar patterns of enhanced adaptation between viral isolates from intermediate and human hosts, contrasted with those isolated from the natural bat reservoir, signifies an indispensable role of the intermediate host in transmission dynamics and spillover events of the virus to human populations. The information regarding avoided codon pairs in SARS-CoV-2, as conferred by the present analysis, promises to be useful for the design of vaccines employing codon pair deoptimization based synthetic attenuated virus engineering.

Predicting the animal hosts of coronaviruses from compositional biases of spike protein and whole genome sequences through machine learning

The COVID-19 pandemic has demonstrated the serious potential for novel zoonotic coronaviruses to emerge and cause major outbreaks. The immediate animal origin of the causative virus, SARS-CoV-2, remains unknown, a notoriously challenging task for emerging disease investigations. Coevolution with hosts leads to specific evolutionary signatures within viral genomes that can inform likely animal origins. We obtained a set of 650 spike protein and 511 whole genome nucleotide sequences from 222 and 185 viruses belonging to the family Coronaviridae, respectively. We then trained random forest models independently on genome composition biases of spike protein and whole genome sequences, including dinucleotide and codon usage biases in order to predict animal host (of nine possible categories, including human). In hold-one-out cross-validation, predictive accuracy on unseen coronaviruses consistently reached ~73%, indicating evolutionary signal in spike proteins to be just as informative as whole genome sequences. However, different composition biases were informative in each case. Applying optimised random forest models to classify human sequences of MERS-CoV and SARS-CoV revealed evolutionary signatures consistent with their recognised intermediate hosts (camelids, carnivores), while human sequences of SARS-CoV-2 were predicted as having bat hosts (suborder Yinpterochiroptera), supporting bats as the suspected origins of the current pandemic. In addition to phylogeny, variation in genome composition can act as an informative approach to predict emerging virus traits as soon as sequences are available. More widely, this work demonstrates the potential in combining genetic resources with machine learning algorithms to address long-standing challenges in emerging infectious diseases.

Molecular features similarities between SARS-CoV-2, SARS, MERS and key human genes could favour the viral infections and trigger collateral effects

In December 2019, rising pneumonia cases caused by a novel β-coronavirus (SARS-CoV-2) occurred in Wuhan, China, which has rapidly spread worldwide, causing thousands of deaths. The WHO declared the SARS-CoV-2 outbreak as a public health emergency of international concern, since then several scientists are dedicated to its study. It has been observed that many human viruses have codon usage biases that match highly expressed proteins in the tissues they infect and depend on the host cell machinery for the replication and co-evolution. In this work, we analysed 91 molecular features and codon usage patterns for 339 viral genes and 463 human genes that consisted of 677,873 codon positions. Hereby, we selected the highly expressed genes from human lung tissue to perform computational studies that permit to compare their molecular features with those of SARS, SARS-CoV-2 and MERS genes. The integrated analysis of all the features revealed that certain viral genes and overexpressed human genes have similar codon usage patterns. The main pattern was the A/T bias that together with other features could propitiate the viral infection, enhanced by a host dependant specialization of the translation machinery of only some of the overexpressed genes. The envelope protein E, the membrane glycoprotein M and ORF7 could be further benefited. This could be the key for a facilitated translation and viral replication conducting to different comorbidities depending on the genetic variability of population due to the host translation machinery. This is the first codon usage approach that reveals which human genes could be potentially deregulated due to the codon usage similarities between the host and the viral genes when the virus is already inside the human cells of the lung tissues. Our work leaded to the identification of additional highly expressed human genes which are not the usual suspects but might play a role in the viral infection and settle the basis for further research in the field of human genetics associated with new viral infections. To identify the genes that could be deregulated under a viral infection is important to predict the collateral effects and determine which individuals would be more susceptible based on their genetic features and comorbidities associated.

A comparative genomics-based study of positive strand RNA viruses emphasizing on SARS-CoV-2 utilizing dinucleotide signature, codon usage and codon context analyses

The novel corona virus disease or COVID-19 caused by a positive strand RNA virus (PRV) called SARS-CoV-2 is plaguing the entire planet as we conduct this study. In this study a multifaceted analysis was carried out employing dinucleotide signature, codon usage and codon context to compare and unravel the genomic as well as genic characteristics of the SARS-CoV-2 isolates and how they compare to other PRVs which represents some of the most pathogenic human viruses. The main emphasis of this study was to comprehend the codon biology of the SARS-CoV-2 in the backdrop of the other PRVs like Poliovirus, Japanese encephalitis virus, Hepatitis C virus, Norovirus, Rubella virus, Semliki Forest virus, Zika virus, Dengue virus, Human rhinoviruses and the Betacoronaviruses since codon usage pattern along with the nucleotide composition prevalent within the viral genome helps to understand the biology and evolution of viruses. Our results suggest discrete genomic dinucleotide signature within the PRVs. Some of the genes from the different SARS-CoV-2 isolates were also found to demonstrate heterogeneity in terms of their dinucleotide signature. The SARS-CoV-2 isolates also demonstrated a codon context trend characteristically dissimilar to the other PRVs. The findings of this study are expected to contribute to the developing global knowledge base in countering COVID-19.

The Long-Term Evolutionary History of Gradual Reduction of CpG Dinucleotides in the SARS-CoV-2 Lineage

Simple Summary

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) caused the coronavirus disease 2019 (COVID-19), a pandemic that infected over 81 million people worldwide. This has led the scientific community to characterize the genome of this virus, including its nucleotide composition. Investigation of the dinucleotide frequency revealed that the proportion of CG dinucleotides (CpG) is highly reduced in the viral genomes. Since CpG dinucleotides is the target site for the host antiviral zinc finger protein, it has been suggested that the reduction in the proportion of CpG is the viral response to escape from the host defense machinery. In the present study, we investigated the time of origin of reduction in the CpG content. Whole genome analyses based on all representative viral genomes of the group Betacoronavirus revealed that the CpG content in the lineage of SARS-CoV-2 has been progressively declining over the past 1213 years. The depletion of CpG was found to occur at neutral—as well as selectively constrained—positions of the viral genomes.

Abstract

Recent studies suggested that the fraction of CG dinucleotides (CpG) is severely reduced in the genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The CpG deficiency was predicted to be the adaptive response of the virus to evade degradation of the viral RNA by the antiviral zinc finger protein that specifically binds to CpG nucleotides. By comparing all representative genomes belonging to the genus Betacoronavirus, this study examined the potential time of origin of CpG depletion. The results of this investigation revealed a highly significant correlation between the proportions of CpG nucleotide (CpG content) of the betacoronavirus species and their times of divergence from SARS-CoV-2. Species that are distantly related to SARS-CoV-2 had much higher CpG contents than that of SARS-CoV-2. Conversely, closely related species had low CpG contents that are similar to or slightly higher than that of SARS-CoV-2. These results suggest a systematic and continuous reduction in the CpG content in the SARS-CoV-2 lineage that might have started since the Sarbecovirus + Hibecovirus clade separated from Nobecovirus, which was estimated to be 1213 years ago. This depletion was not found to be mediated by the GC contents of the genomes. Our results also showed that the depletion of CpG occurred at neutral positions of the genome as well as those under selection. The latter is evident from the progressive reduction in the proportion of arginine amino acid (coded by CpG dinucleotides) in the SARS-CoV-2 lineage over time. The results of this study suggest that shedding CpG nucleotides from their genome is a continuing process in this viral lineage, potentially to escape from their host defense mechanisms.

Potential therapeutic approaches of microRNAs for COVID-19: Challenges and opportunities

The coronavirus disease 2019 (COVID-19) emerges as current outbreak cause by Novel Severe Acute Respiratory Syndrome Corona Virus-2 (SARS-CoV-2). This infection affects respiratory system and provides uncontrolled systemic inflammatory response as cytokine storm. The main concern about SARS-CoV-2 pandemic is high viral pathogenicity with no specific drugs. MicroRNAs (miRs) as small non-coding RNAs (21–25 ​nt) regulate gene expression. The SARS-CoV-2 encoded-miRs affect human genes that involved in transcription, translation, apoptosis, immune response and inflammation. Also, they alter self-gene regulation and hijacked host miRs that provide protective environment to maintain its latency. On the other hand, Host miRs play critical role in viral gene expression to restrict infection. Over expression/inhibition of miRs might result in cell cycle irregularity, impaired immune response or cancer. In this manner, exact role of each miR should be specified. Mimic encoded-miRs like antagomirs showed successful result in phases of clinical trial prevent from negative effects of viral encoded-miRs. Products of mimic miRs are inexpensive corresponds to synthesis of primer; they are short and nanoscale in size. Although SARS-CoV-2 genome is undergoing evaluation, detection of exact molecular pathogenesis open up opportunities to for vaccine development. Salivaomics can evaluate SARS-CoV-2 genome, transcriptome, proteome and biomarkers like miRs in oral related and cancer disease. In this review, we studied the challenge and opportunities of miRs in therapeutic approach for SARS-CoV-2 infection, then overviewed the role of miRs in saliva droplet during SARS-CoV-2 infection and related cancer.

A Path toward SARS-CoV-2 Attenuation: Metabolic Pressure on CTP Synthesis Rules the Virus Evolution

In the context of the COVID-19 pandemic, we describe here the singular metabolic background that constrains enveloped RNA viruses to evolve toward likely attenuation in the long term, possibly after a step of increased pathogenicity. Cytidine triphosphate (CTP) is at the crossroad of the processes allowing SARS-CoV-2 to multiply, because CTP is in demand for four essential metabolic steps. It is a building block of the virus genome, it is required for synthesis of the cytosine-based liponucleotide precursors of the viral envelope, it is a critical building block of the host transfer RNAs synthesis and it is required for synthesis of dolichol-phosphate, a precursor of viral protein glycosylation. The CCA 3'-end of all the transfer RNAs required to translate the RNA genome and further transcripts into the proteins used to build active virus copies is not coded in the human genome. It must be synthesized de novo from CTP and ATP. Furthermore, intermediary metabolism is built on compulsory steps of synthesis and salvage of cytosine-based metabolites via uridine triphosphate that keep limiting CTP availability. As a consequence, accidental replication errors tend to replace cytosine by uracil in the genome, unless recombination events allow the sequence to return to its ancestral sequences. We document some of the consequences of this situation in the function of viral proteins. This unique metabolic setup allowed us to highlight and provide a raison d'être to viperin, an enzyme of innate antiviral immunity, which synthesizes 3'-deoxy-3',4'-didehydro-CTP as an extremely efficient antiviral nucleotide.

Natural Selection Plays an Important Role in Shaping the Codon Usage of Structural Genes of the Viruses Belonging to the Coronaviridae Family

Viruses belonging to the Coronaviridae family have a single-stranded positive-sense RNA with a poly-A tail. The genome has a length of ~29.9 kbps, which encodes for genes that are essential for cell survival and replication. Different evolutionary constraints constantly influence the codon usage bias (CUB) of different genes. A virus optimizes its codon usage to fit the host environment on which it savors. This study is a comprehensive analysis of the CUB for the different genes encoded by viruses of the Coronaviridae family. Different methods including relative synonymous codon usage (RSCU), an Effective number of codons (ENc), parity plot 2, and Neutrality plot, were adopted to analyze the factors responsible for the genetic evolution of the Coronaviridae family. Base composition and RSCU analyses demonstrated the presence of A-ended and U-ended codons being preferred in the 3rd codon position and are suggestive of mutational selection. The lesser ENc value for the spike ‘S’ gene suggests a higher bias in the codon usage of this gene compared to the other structural genes. Parity plot 2 and neutrality plot analyses demonstrate the role and the extent of mutational and natural selection towards the codon usage pattern. It was observed that the structural genes of the Coronaviridae family analyzed in this study were at the least under 84% influence of natural selection, implying a major role of natural selection in shaping the codon usage.

Phylogenomics Analysis of SARS-CoV2 Genomes Reveals Distinct Selection Pressure on Different Viral Strains

We are witnessing a tremendous outbreak of a novel coronavirus (SARS-CoV2) across the globe. Upon exposure to different population and changing environment, the viral strain might experience different mutational bias that leads to genetic diversity among the viral population. Also, the diversification can be influenced by distinct selection pressure on different viral genomes. We have carried out a comparative genomic analysis of 82 SARS-CoV2 genomes. We have evaluated their evolutionary divergence, substitution pattern, and rates. Viral genomes under distinct selection pressure have been identified. Sites that experience strong selection pressure also have been identified. Our result shows that the translational preference of a few codons is strongly correlated with the mutational bias imposed by genome compositional constraint and influenced by natural selection. Few genomes are evolving with a higher mutational rate with a distinct signature of nucleotide substitution in comparison to others. Four viral strains are under the effect of purifying selection, while nine SARS-CoV2 genomes are under strong positive selection bias. Site analysis indicates a strong positive selection pressure on two codon positions at 3606th and 8439th positions. Our study elucidates adaptation of few SARS-CoV2 viral strain during the outbreak shaping by natural selection and genomic compositional constraints.

Clinical Characteristics, Diagnosis, and Treatment of Major Coronavirus Outbreaks

Human coronavirus infections have been known to cause mild respiratory illness. It changed in the last two decades as three global outbreaks by coronaviruses led to significant mortality and morbidity. SARS CoV-1 led to the first epidemic of the twenty first century due to coronavirus. SARS COV-1 infection had a broad array of symptoms with respiratory and gastrointestinal as most frequent. The last known case was reported in 2004. Middle East respiratory syndrome coronavirus (MERS-CoV) led to the second outbreak in 2012, and case fatality was much higher than SARS. MERS-CoV has a wide array of clinical presentations from mild, moderate to severe, and some patients end up with acute respiratory distress syndrome (ARDS). The third and recent outbreak by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) started in December 2019, which lead to a global pandemic. Patients with SARS-CoV2 infection can be asymptomatic or have a range of symptoms with fever, cough, and shortness of breath being most common. Reverse transcriptase-Polymerase chain reaction (RT-PCR) is a diagnostic test of choice for SARS CoV-1, MERS-CoV, and SARS CoV-2 infections. This review aims to discuss epidemiological, clinical features, diagnosis, and management of human coronaviruses with a focus on SARS CoV-1, MERS-CoV, and SARS CoV-2.

Deciphering the co-adaptation of codon usage between respiratory coronaviruses and their human host uncovers candidate therapeutics for COVID-19

Coronavirus disease 2019 (COVID-19) has caused thousands of deaths worldwide and has become an urgent public health concern. The extraordinary interhuman transmission of this disease has urged scientists to examine the various facets of its pathogenic agent, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Herein, based on publicly available genomic data, we analyzed the codon usage co-adaptation profiles of SARS-CoV-2 and other respiratory coronaviruses (CoVs) with their human host, identified CoV-responsive human genes and their functional roles on the basis of both the relative synonymous codon usage (RSCU)-based correlation of viral genes with human genes and differential gene expression analysis, and predicted potential drugs for COVID-19 treatment based on these genes. The relatively high codon adaptation index (CAI) values (>0.70) signposted the gene expressivity efficiency of CoVs in human. The ENc-GC3 plot indicated that SARS-CoV-2 genome was under strict selection pressure while SARS-CoV and MERS-CoV were under selection and mutational pressures. The RSCU-based correlation analysis indicated that the viral genomes shared similar codons with a panoply of human genes. The merging of RSCU-based correlation data and SARS-CoV-2-responsive differentially expressed genes allowed the identification of human genes potentially affected by SARS-CoV-2 infection. Functional enrichment analysis indicated that these genes were enriched in biological processes and pathways related to host response to viral infection and immune response. Using the drug-gene interaction database, we screened a list of drugs that could target these genes as potential COVID-19 therapeutics. Our findings not only will contribute in vaccine development but also provide a useful set of drugs that could guide practitioners in strategical monitoring of COVID-19. We recommend practitioners to scrupulously screen this list of predicted drugs in order to authenticate those qualified for treating COVID-19 symptoms.

Evolution of Transmissible Gastroenteritis Virus (TGEV): A Codon Usage Perspective

Transmissible gastroenteritis virus (TGEV) is a coronavirus associated with diarrhea and high mortality in piglets. To gain insight into the evolution and adaptation of TGEV, a comprehensive analysis of phylogeny and codon usage bias was performed. The phylogenetic analyses of maximum likelihood and Bayesian inference displayed two distinct genotypes: genotypes I and II, and genotype I was classified into subtypes Ia and Ib. The compositional properties revealed that the coding sequence contained a higher number of A/U nucleotides than G/C nucleotides, and that the synonymous codon third position was A/U-enriched. The principal component analysis based on the values of relative synonymous codon usage (RSCU) showed the genotype-specific codon usage patterns. The effective number of codons (ENC) indicated moderate codon usage bias in the TGEV genome. Dinucleotide analysis showed that CpA and UpG were over-represented and CpG was under-represented in the coding sequence of the TGEV genome. The analyses of Parity Rule 2 plot, ENC-plot, and neutrality plot displayed that natural selection was the dominant evolutionary driving force in shaping codon usage preference in genotypes Ia and II. In addition, natural selection played a major role, while mutation pressure had a minor role in driving the codon usage bias in genotype Ib. The codon adaptation index (CAI), relative codon deoptimization index (RCDI), and similarity index (SiD) analyses suggested that genotype I might be more adaptive to pigs than genotype II. Current findings contribute to understanding the evolution and adaptation of TGEV.

Characterization of codon usage pattern in SARS-CoV-2

The outbreak of coronavirus disease 2019 (COVID-19) due to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has posed significant threats to international health. The genetic traits as well as evolutionary processes in this novel coronavirus are not fully characterized, and their roles in viral pathogenesis are yet largely unknown. To get a better picture of the codon architecture of this newly emerging coronavirus, in this study we perform bioinformatic analysis, based on publicly available nucleotide sequences of SARS-CoV-2 along with those of other members of human coronaviruses as well as non-human coronaviruses in different hosts, to take a snapshot of the genome-wide codon usage pattern of SARS-CoV-2 and uncover that all over-represented codons end with A/U and this newly emerging coronavirus has a relatively low codon usage bias, which is shaped by both mutation pressure and natural selection. Additionally, there is slight variation in the codon usage pattern among the SARS-CoV-2 isolates from different geo-locations. Furthermore, the overall codon usage pattern of SARS-CoV-2 is generally similar to that of its phylogenetic relatives among non-human betacoronaviruses such as RaTG13. Taken together, we comprehensively analyze the characteristics of codon usage pattern in SARS-CoV-2 via bioinformatic approaches. The information from this research may not only be helpful to get new insights into the evolution of SARS-CoV-2, but also have potential value for developing coronavirus vaccines.

Analysis of codon usage of severe acute respiratory syndrome corona virus 2 (SARS-CoV-2) and its adaptability in dog

•

Codon analysis reveal natural selection and other undefined factors dominates the overall codon usage bias in SARS-CoV-2 rather than mutational pressure.

•

The host adaptation potential of SARS-CoV-2 is more in human as compared to dog.

Severe acute respiratory syndrome corona virus 2 (SARS-CoV-2) is recognized as one of the life-threatening viruses causing the most destructive pandemic in this century. The genesis of this virus is still unknown. To elucidate its molecular evolution and regulation of gene expression, the knowledge of codon usage is a pre-requisite. In this study, an attempt was made to document the genome-wide codon usage profile and the various factors influencing the codon usage patterns of SARS-CoV-2 in human and dog. The SARS-CoV-2 genome showed relative abundance of A and U nucleotides and relative synonymous codon usage analysis revealed that the preferred synonymous codons mostly end with A/U. The analysis of ENc-GC3s, Neutrality and Parity rule 2 plots indicated that natural selection and other undefined factors dominate the overall codon usage bias in SARS-CoV-2 whereas the impact of mutation pressure is comparatively minor. The codon adaptation index and relative codon deoptimization index of SARS-CoV-2 deciphered that human is more favoured host for adaptation compared to dog. These results enhance our understanding of the factors involved in evolution of the novel human SARS-CoV-2 and its adaptability in dog.

Natural Products as Potential Leads Against Coronaviruses: Could They be Encouraging Structural Models Against SARS-CoV-2?

New coronavirus referred to SARS-CoV-2 has caused a worldwide pandemic (COVID-19) declared by WHO. Coronavirus disease 2019 (COVID-19) is an infectious disease with severe acute respiratory syndrome caused by coronavirus-2 (SARS-CoV-2). SARS-CoV-2 is akin to SARS-CoV, which was the causative agent of severe acute respiratory syndrome (SARS) in 2002 as well as to that of Middle East respiratory syndrome (MERS) in 2012. SARS-CoV-2 has been revealed to belong to Coronaviridiae family as a member of β-coronaviruses. It has a positive-sense single-stranded RNA with the largest RNA genome. Since its genomic sequence has a notable similarity to that of SARS-CoV, antiviral drugs used to treat SARS and MERS are now being also applied for COVID-19 treatment. In order to combat SARS-CoV-2, many drug and vaccine development studies at experimental and clinical levels are currently conducted worldwide. In this sense, medicinal plants and the pure natural molecules isolated from plants have been reported to exhibit significant inhibitory antiviral activity against SARS-CoV and other types of coronaviruses. In the present review, plant extracts and natural molecules with the mentioned activity are discussed in order to give inspiration to researchers to take these molecules into consideration against SARS-CoV-2.

Implications of SARS-CoV-2 Mutations for Genomic RNA Structure and Host microRNA Targeting

The SARS-CoV-2 virus is a recently-emerged zoonotic pathogen already well adapted to transmission and replication in humans. Although the mutation rate is limited, recently introduced mutations in SARS-CoV-2 have the potential to alter viral fitness. In addition to amino acid changes, mutations could affect RNA secondary structure critical to viral life cycle, or interfere with sequences targeted by host miRNAs. We have analysed subsets of genomes from SARS-CoV-2 isolates from around the globe and show that several mutations introduce changes in Watson-Crick pairing, with resultant changes in predicted secondary structure. Filtering to targets matching miRNAs expressed in SARS-CoV-2-permissive host cells, we identified ten separate target sequences in the SARS-CoV-2 genome; three of these targets have been lost through conserved mutations. A genomic site targeted by the highly abundant miR-197-5p, overexpressed in patients with cardiovascular disease, is lost by a conserved mutation. Our results are compatible with a model that SARS-CoV-2 replication within the human host is constrained by host miRNA defences. The impact of these and further mutations on secondary structures, miRNA targets or potential splice sites offers a new context in which to view future SARS-CoV-2 evolution, and a potential platform for engineering conditional attenuation to vaccine development, as well as providing a better understanding of viral tropism and pathogenesis.

Natural compounds from Clerodendrum spp. as possible therapeutic candidates against SARS-CoV-2: An in silico investigation

The COVID-19 pandemic caused by the novel coronavirus SARS-CoV-2 has rattled global public health, with researchers struggling to find specific therapeutic solutions. In this context, the present study employed an in silico approach to assess the inhibitory potential of the phytochemicals obtained from GC-MS analysis of twelve Clerodendrum species against the imperative spike protein, main protease enzyme M^pro and RNA-dependent RNA polymerase (RdRp) of SARS-CoV-2. An extensive molecular docking investigation of the phytocompounds at the active binding pockets of the viral proteins revealed promising inhibitory potential of the phytochemicals taraxerol, friedelin and stigmasterol. Decent physicochemical attributes of the compounds in accordance with Lipinski’s rule of five and Veber’s rule further established them as potential therapeutic candidates against SARS-CoV-2. Molecular mechanics-generalized Born surface area (MM-GBSA) binding free energy estimation revealed that taraxerol was the most promising candidate displaying the highest binding efficacy with all the concerned SARS-CoV-2 proteins included in the present analysis. Our observations were supported by robust molecular dynamics simulations of the complexes of the viral proteins with taraxerol for a timescale of 40 nanoseconds. It was striking to note that taraxerol exhibited better binding energy scores with the concerned viral proteins than the drugs that are specifically targeted against them. The present results promise to provide new avenues to further evaluate the potential of the phytocompound taraxerol in vitro and in vivo towards its successful deployment as a SARS-CoV-2 inhibitor and combat the catastrophic COVID-19.

Communicated by Ramaswamy H. Sarma