Evolution of the SARS-CoV-2 pandemic in India.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a bat-derived betacoronavirus, that emerged around December 2019. In spite of the lesser genomic diversity of CoVs in general, a steady accumulation of mutations spread over its genome have been noted, resulting in the emergence of several clades and lineages. Majority of these mutations are random and non-functional changes; however a few variants of concern (VOC) and variants of interest (VOI) designated by the WHO since late 2020 have implications to diagnostics, pathogenicity and immune escape. This review discusses the various nomenclatures depicting the SARS-CoV-2 evolution, the designated VOCs and VOIs and the mutations characterizing these variants. The evolution of SARS-CoV-2 in India and the implications to vaccine efficacy and breakthrough infections is also addressed.

Perspectives for repurposing drugs for the coronavirus disease 2019.

The newly emerged 2019 novel coronavirus (CoV), named as severe acute respiratory syndrome CoV-2 (SARS-CoV-2), like SARS-CoV (now, SARS-CoV-1) and Middle East respiratory syndrome CoV (MERS-CoV), has been associated with high infection rates with over 36,405 deaths. In the absence of approved marketed drugs against coronaviruses, the treatment and management of this novel CoV disease (COVID-19) worldwide is a challenge. Drug repurposing that has emerged as an effective drug discovery approach from earlier approved drugs could reduce the time and cost compared to de novo drug discovery. Direct virus-targeted antiviral agents target specific nucleic acid or proteins of the virus while host-based antivirals target either the host innate immune responses or the cellular machineries that are crucial for viral infection. Both the approaches necessarily interfere with viral pathogenesis. Here we summarize the present status of both virus-based and host-based drug repurposing perspectives for coronaviruses in general and the SARS-CoV-2 in particular.

Full-genome sequences of the first two SARS-CoV-2 viruses from India.

Background & objectives: Since December 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has globally affected 195 countries. In India, suspected cases were screened for SARS-CoV-2 as per the advisory of the Ministry of Health and Family Welfare. The objective of this study was to characterize SARS-CoV-2 sequences from three identified positive cases as on February 29, 2020. Methods: Throat swab/nasal swab specimens for a total of 881 suspected cases were screened by E gene and confirmed by RdRp (1), RdRp (2) and N gene real-time reverse transcription-polymerase chain reactions and next-generation sequencing. Phylogenetic analysis, molecular characterization and prediction of B- and T-cell epitopes for Indian SARS-CoV-2 sequences were undertaken. Results: Three cases with a travel history from Wuhan, China, were confirmed positive for SARS-CoV-2. Almost complete (29,851 nucleotides) genomes of case 1, case 3 and a fragmented genome for case 2 were obtained. The sequences of Indian SARS-CoV-2 though not identical showed high (~99.98%) identity with Wuhan seafood market pneumonia virus (accession number: NC 045512). Phylogenetic analysis showed that the Indian sequences belonged to different clusters. Predicted linear B-cell epitopes were found to be concentrated in the S1 domain of spike protein, and a conformational epitope was identified in the receptor-binding domain. The predicted T-cell epitopes showed broad human leucocyte antigen allele coverage of A and B supertypes predominant in the Indian population. Interpretation & conclusions: The two SARS-CoV-2 sequences obtained from India represent two different introductions into the country. The genetic heterogeneity is as noted globally. The identified B- and T-cell epitopes may be considered suitable for future experiments towards the design of vaccines and diagnostics. Continuous monitoring and analysis of the sequences of new cases from India and the other affected countries would be vital to understand the genetic evolution and rates of substitution of the SARS-CoV-2.

Chikungunya phylogeography reveals persistent global transmissions of the Indian Ocean Lineage from India in association with mutational fitness.

Since the resurgence of chikungunya virus (CHIKV) in India in 2005, the Indian subcontinent sublineage of the Indian Ocean lineage (IOL) has continued transmission in India and also radiation from India causing additional outbreaks in surrounding countries. This study was undertaken for an in-depth understanding of the evolutionary dynamics of the IOL, the global transmission routes in the Indian context and possible association with mutational fitness. The whole genome sequencing of Indian isolates representing CHIKV outbreaks (2014-2018) from selected States of India was carried out, followed by phylogeography analysis of the IOL using the Bayesian Markov chain Monte Carlo method and selection pressure analysis. Phylogeography analysis of IOL strains revealed indigenous evolution in India at least at three time points, with specific mutations that conferred viral fitness in the Aedes vector species. Further dispersal of the strains from India was noted to neighbouring and distant countries with multiple exportations to Sri Lanka, Bangladesh and China. The study reveals India as an endemic reservoir for CHIKV and persistent global transmissions from India. Though natural selection does not appear to play a major role in establishment of the IOL, sustainable efforts towards vector control can help address the issues.

Full genome based sequence and structural characterization of an unusual group A rotavirus G12P[11] isolated from neonates in Pune, western India.

Studies conducted at neonatal intensive care units in Pune, western India, suggested early exposure to rotaviruses and predominance of unusual human-bovine-like G12P[11] strains. The whole genome sequencing and phylogenetic analyses of a naturally attenuated, culture adapted neonatal strain, (NIV-1740121) revealed multiple-gene reassortment events, containing ROTAVAC® vaccine strain, 116E-like VP4, VP6, NSP3, NSP5 genes, VP7 gene of G12 origin and VP3 gene of porcine ancestry in a human Wa-like backbone. Analysis of 3D structure modeling of the VP7 and VP4 proteins with respect to 116E suggested amino acid variations in the major neutralizing epitopes of VP7, contributed to a modified charge density. Visualization of receptor-glycan interaction structures of NIV-1740121 and 116E VP8* showed type I glycan binds with a similar conformation at the same active site as represented in the available crystal structure of G10P[11] VP8*. The study adds to the knowledge of age restricted tropism of P[11] strains in neonates.

Phylogeography of Kyasanur Forest Disease virus in India (1957-2017) reveals evolution and spread in the Western Ghats region.

The Kyasanur Forest Disease (KFD) has become a major public health problem in the State of Karnataka, India where the disease was first identified and in Tamil Nadu, Maharashtra, Kerala, and Goa covering the Western Ghats region of India. The incidence of positive cases and distribution of the Kyasanur Forest Disease virus (KFDV) in different geographical regions raises the need to understand the evolution and spatiotemporal transmission dynamics. Phylogeography analysis based on 48 whole genomes (46 from this study) and additionally 28 E-gene sequences of KFDV isolated from different regions spanning the period 1957-2017 was thus undertaken. The mean evolutionary rates based the E-gene was marginally higher than that based on the whole genomes. A subgroup of KFDV strains (2006-2017) differing from the early Karnataka strains (1957-1972) by ~2.76% in their whole genomes and representing spread to different geographical areas diverged around 1980. Dispersal from Karnataka to Goa and Maharashtra was indicated. Maharashtra represented a new source for transmission of KFDV since ~2013. Significant evidence of adaptive evolution at site 123 A/T located in the vicinity of the envelope protein dimer interface may have functional implications. The findings indicate the need to curtail the spread of KFDV by surveillance measures and improved vaccination strategies.

Whole-genome-based characterization of three human Rotavirus C strains isolated from gastroenteritis outbreaks in Western India and a provisional intra-genotypic lineage classification system.

The number of whole-genome sequences of human rotavirus C (RVC) strains available in public databases is recently increasing. Thus far from India only a single whole genome of human RVC of a sporadic case was available. In this study, nearly full-length genome sequencing and phylogenetic analyses of three RVC strains isolated from three different gastroenteritis outbreaks during 2010-2014 in Western India was carried out. Further, an intra-genotypic lineage classification system for human RVCs based on the nucleotide divergence cut-off values was proposed by using the algorithm of the Rotavirus Classification Working Group. Two lineages could be defined for all the genes except the VP7 gene and the M3 VP3 genotype. Provisional classification of the lineages indicated the absence of reassortment events in the genomic constellation of Indian strains, contrary to earlier reports. The comparatively higher variability of the NSP1, NSP3, NSP5 and M2 VP3 genotype, emphasizes their utility in lineage classification.

Molecular mechanism by which residues at position 481 and 546 of measles virus hemagglutinin protein define CD46 receptor binding using a molecular docking approach.

The hemagglutinin (H) protein of measles viruses (MeV) mediates binding to the cellular receptors, CD46,human signaling lymphocyte activation molecule and nectin-4. Vaccine strains primarily contain H-proteins possessing MeV-H: Y481 and can utilize CD46. Reports suggest that a single amino acid change in MeV-H at position 481 in wild type strains renders them inefficient in utilizing CD46. The in-depth molecular mechanism by which substitutions at 481 and another reported critical residue position 546 affects CD46 binding affinity however remains elusive. We used molecular docking studies of CD46 with MeV-H possessing Y481 N/D to understand the in-depth molecular mechanism involved. It was found that loss in either of the hydrogen bond (H-bond) contacts (MeV-H:481-CD46:65, MeV-H:546-CD46:63) in the central contact region prevented efficient CD46 binding. Y481 N could form the specific H-bond, while G546S H-bond could be formed only in conjunction with Y481, revealing the significance of these residues in determining CD46 receptor binding potential. Elucidating the underlying molecular mechanism of receptor usage by the MeV has implications to understanding cellular tropism, viral pathogenesis and therapy.

Selection of avian influenza A (H9N2) virus with reduced susceptibility to neuraminidase inhibitors oseltamivir and zanamivir.

Identification of amino-acid substitutions in the neuraminidase (NA) of low-pathogenic avian influenza (AI) H9N2 viruses is important to study the susceptibility to NA inhibitors (NAI). To identify mutations under NAI selective pressure, the virus was serially passaged with increasing levels of either oseltamivir or zanamivir in ovo, and the growth of the viruses in the presence and absence of NAI's compared. Mutations R292 K in the presence of oseltamivir and E119D in presence of zanamivir were observed within passage one and two respectively. The R292 K mutation reduced oseltamivir susceptibility significantly (2,523-fold) and moderately reduced susceptibility to zanamivir. The E119D mutation significantly reduced susceptibility to zanamivir (415-fold) and remained susceptible to oseltamivir. Genetic stability of the mutations assessed by serial passages of the mutant viruses in eggs without drug pressure resulted in the loss of these mutations, making the virus susceptible to both the drugs. Molecular modeling and dynamics simulations revealed that the R292 K mutation disrupted oseltamivir binding similar to other group 2 NAs, while a different mechanism was noted for zanamivir binding for both R292 K and E119D mutations. The study highlights the need for regular susceptibility screening of circulating AI viruses.

Global spatiotemporal transmission dynamics of measles virus clade D genotypes in the context of the measles elimination goal 2020 in India.

Measles viruses (MeV) circulating in India mainly belong to genoypes D4 and D8 of clade D. In the context of measles elimination goal 2020 in India, molecular clock and phylogeography studies would help to identify the timescales of evolution and track the transmission pathways of MeV. We used nucleoprotein gene sequences (n = 756) from GenBank, representing 86 countries (1973-2016), to study the spatiotemporal transmission dynamics of clade D. Genotype D4 was introduced into India around 1991 and genotype D8 around 1994. Recent transmissions of the D4 genotype of measles virus (MeV) were noted from India to the United States of America and East Asia region while D8 genotype importations from North America were noted in recent years.

Genetic characterization of chikungunya viruses isolated during the 2015-2017 outbreaks in different states of India, based on their E1 and E2 genes.

During 2015-2017, chikungunya virus (CHIKV) showed a resurgence in several parts of India with Karnataka, Maharashtra and New Delhi accounting for a majority of the cases. E2-E1 gene based characterization revealed Indian subcontinent sublineage strains possessing Aedes aegypti mosquito-adaptive mutations E1: K211E and E2:V264A, with the 211 site positively selected. Novel mutational sites E1: K16E/Q, E1: K132Q/T, E1: S355T, E2: C19R and E2:S185Y could be associated with epitopes or virulence determining domains. The study examines the role of host, vector and viral factors and fills gaps in our molecular epidemiology data for these regions which are known to possess a dynamic population.

Computational analysis of the effect of polymerase acidic (PA) gene mutation F35L in the 2009 pandemic influenza A (H1N1) virus on binding aspects of mononucleotides in the endonuclease domain.

An F35L mutation in the N-terminal domain of the polymerase acidic protein (PA-Nter), which contains the active site of the endonuclease, has been reported to result in higher polymerase activity in mouse-adapted strains of the 2009 pandemic influenza A H1N1 virus. We modeled wild and mutant complexes of uridine 5'-monophosphate (UMP) as the endonuclease substrate and performed molecular dynamics simulations. The results demonstrated that the F35L mutation could result in a changed orientation of a helix containing active site residues and improve the ligand affinity in the mutant strain. This study suggests a molecular mechanism of enhanced polymerase activity.

A molecular modelling approach to understand the effect of co-evolutionary mutations (V344M, I354L) identified in the PB2 subunit of influenza A 2009 pandemic H1N1 virus on m7GTP ligand binding.

The cap binding domain of the polymerase basic 2 (PB2) subunit of influenza polymerases plays a critical role in mediating the 'cap-snatching' mechanism by binding the 5' cap of host pre-mRNAs during viral mRNA transcription. Monitoring variations in the PB2 protein is thus vital for evaluating the pathogenic potential of the virus. Based on selection pressure analysis of PB2 gene sequences of the pandemic H1N1 (pH1N1) viruses of the period 2009-2014, we identified a site, 344V/M, in the vicinity of the cap binding pocket showing evidence of adaptive evolution and another co-evolving residue, 354I/L, in close vicinity. Modelling of the three-dimensional structure of the pH1N1 PB2 cap binding domain, docking of the pre-mRNA cap analogue m7GTP and molecular dynamics simulation studies of the docked complexes performed for four PB2 variants observed showed that the complex possessing V344M with I354L possessed better ligand binding affinity due to additional hydrogen bond contacts between m7GTP and the key residues His432 and Arg355 that was attributed to a displacement of the 424 loop and a flip of the side chain of Arg355, respectively. The co-evolutionary mutations identified (V344M, I354L) were found to be established in the PB2 gene of the pH1N1 viral population over the period 2010-2014. The study demonstrates the molecular basis for the enhanced m7GTP ligand binding affinity with the 344M-354L synergistic combination in PB2. Furthermore, the insight gained into understanding the molecular mechanism of cap binding in pH1N1 viruses may be useful for designing novel drugs targeting the PB2 cap binding domain.

Genomic characterization of coxsackievirus type B3 strains associated with acute flaccid paralysis in south-western India.

Acute flaccid paralysis (AFP) associated with coxsackievirus type B3 (CV-B3) of the species Enterovirus B is an emerging concern worldwide. Although CV-B3-associated AFP in India has been demonstrated previously, the genomic characterization of these strains is unreported. Here, CV-B3 strains detected on the basis of the partial VP1 gene in 10 AFP cases and five asymptomatic contacts identified from different regions of south-western India during 2009-2010 through the Polio Surveillance Project were considered for complete genome sequencing and characterization. Phylogenetic analysis of complete VP1 gene sequences of global CV-B3 strains classified Indian CV-B3 strains into genogroup GVI, along with strains from Uzbekistan and Bangladesh, and into a new genogroup, GVII. Genomic divergence between genogroups of the study strains was 14.4 % with significantly lower divergence (1.8 %) within GVI (n = 12) than that within GVII (8.5 %) (n = 3). The strains from both AFP cases and asymptomatic contacts, identified mainly in coastal Karnataka and Kerala, belonged to the dominant genogroup GVI, while the GVII strains were recovered from AFP cases in north interior Karnataka. All study strains carried inter-genotypic recombination with the structural region similar to reference CV-B3 strains, and 5' non-coding regions and non-structural regions closer to other enterovirus B types. Domain II structures of 5' non-coding regions, described to modulate virus replication, were predicted to have varied structural folds in the two genogroups and were attributed to differing recombination patterns. The results indicate two distinct genomic compositions of CV-B3 strains circulating in India and suggest the need for concurrent analysis of viral and host factors to further understand the varied manifestations of their infections.

Molecular dynamics simulation of the effects of single (S221P) and double (S221P and K216E) mutations in the hemagglutinin protein of influenza A H5N1 virus: a study on host receptor specificity.

Avian influenza viruses of subtype H5N1 circulating in animals continue to pose threats to human health. The binding preference of the viral surface protein hemagglutinin (HA) to sialosaccharides of receptors is an important area for understanding mutations in the receptor binding site that could be the cause for avian-to-human transmission. In the present work, we studied the effect of two receptor binding site mutations, S221P singly and in combination with another mutation K216E in the HA protein of influenza A H5N1 viruses. Docking of sialic acid ligands corresponding to both avian and human receptors and molecular dynamics simulations of the complexes for wild and mutant strains of H5N1 viruses were carried out. The H5N1 strain possessing the S221P mutation indicated decreased binding to α2,3-linked sialic acids (avian receptor, SAα2,3Gal) when compared to the binding of the wild-type strain that did not possess the HA-221 mutation. The binding to α2,6-linked sialic acids (human receptor, SAα2,6Gal) was found to be comparable, indicating that the mutant strain shows limited dual receptor specificity. On the other hand, the S221P mutation in synergism with the K216E mutation in the binding site, resulted in increased binding affinity for SAα2,6Gal when compared to SAα2,3Gal, indicative of enhanced binding to human receptors. The in-depth study of the molecular interactions in the docked complexes could explain how co-occurring mutations in the HA viral protein can aid in providing fitness advantage to the virus, in the context of host receptor specificity in emerging variants of H5N1 influenza viruses.

Phylogeographic analysis of Japanese encephalitis virus in India (1956-2012).

Japanese encephalitis virus (JEV) isolates from India phylogenetically belong to two genotypes, III and I. We used envelope gene sequences from GenBank, representing different states of India and other countries, to study the spatiotemporal transmission histories of these two JEV genotypes separately. Genotype III was found to have been successively introduced in the 1930s, 1950s and 1960s, followed by genotype I twice around 2003-2006. Changes in JEV disease patterns in India over the last five decades could thus be attributed to multiple introductions of JEV strains from neighboring Asian countries along with increased transmission potential due to altered ecological settings.