Molecular characterization & recombination analysis of complete enterovirus-88 isolated from acute flaccid paralysis cases in India

Background & objectives

Focus on non-polio enteroviruses (NPEVs) causing acute flaccid paralysis (AFP) due to myelitis has increased with the containment of the poliovirus. Enterovirus-B88 (EV-B88) has been associated with the AFP cases in Bangladesh, Ghana, South Africa, Thailand and India. In India, EV-B88 infection was linked to AFP a decade ago; however, to date, no complete genome has been made available. In this study, the complete genome sequence of EV-B88 was identified and reported from two different States (Bihar and Uttar Pradesh) in India using the next-generation sequencing technique.

Methods

Virus isolation was performed on the three AFP suspected cases as per the WHO-recommended protocol. Samples showing cytopathic effects in the human Rhabdocarcinoma were labelled as NPEVs. Next-generation sequencing was performed on these NPEVs to identify the aetiological agent. The contiguous sequences (contigs) generated were identified, and reference-based mapping was performed.

Results

EV-B88 sequences retrieved in our study were found to be 83 per cent similar to the EV-B88 isolate from Bangladesh in 2001 (strain: BAN01-10398; Accession number: AY843306.1). Recombination analyses of these samples demonstrate recombination events with sequences from echovirus-18 and echovirus-30.

Interpretation & conclusions

Recombination events in the EV-B serotypes are known, and this work reconfirms the same for EV-B88 isolates also. This study is a step in increasing the awareness about EV-B88 in India and emphasizes future studies to be conducted in the identification of other types of EV present in India.

Community transmission of SARS-CoV-2 with B.1.1.7 lineage in Mumbai, India

The B.1.1.7 (Alpha) variant has been detected in Mumbai, India during February 2021. Subsequently, we retrieved 43 sequences from specimens of 51 COVID-19 cases from Mumbai. The sequence analysis revealed that the cases were mainly affected with Alpha variant which suggests its role in community transmission of SARS-CoV-2 in Mumbai, India.

Delta variant SARS-CoV-2 infections in pediatric cases during the second wave in India

BACKGROUND

During October 2020, Delta variant was detected for the first time in India and rampantly spread across the globe. It also led to second wave of pandemic in India which affected millions of people. However, there is limited information pertaining to the SARS-CoV-2 strain infecting the children in India.

METHODS

Here, we assessed the SARS-CoV-2 lineages circulating in the pediatric population of India during the second wave of the pandemic. Clinical and demographic details linked with the nasopharyngeal/oropharyngeal swabs (NPS/OPS) collected from SARS-CoV-2 cases (n = 583) aged 0-18 year and tested positive by real-time RT-PCR were retrieved from March to June 2021.

RESULTS

Symptoms were reported among 37.2% of patients and 14.8% reported to be hospitalized. The E gene CT value had significant statistical difference at the point of sample collection when compared to that observed in the sequencing laboratory. Out of these 512 sequences 372 were VOCs, 51 were VOIs. Most common lineages observed were Delta, followed by Kappa, Alpha and B.1.36, seen in 65.82%, 9.96%, 6.83% and 4.68%, respectively in the study population.

CONCLUSION

Overall, it was observed that Delta strain was the leading cause of SARS-CoV-2 infection in Indian children during the second wave of the pandemic. We emphasize on the need of continuous genomic surveillance in SARS-CoV-2 infection even amongst children.

Zika a Vector Borne Disease Detected in Newer States of India Amidst the COVID-19 Pandemic

Background

During the second wave of the COVID-19 pandemic, outbreaks of Zika were reported from Kerala, Uttar Pradesh, and Maharashtra, India in 2021. The Dengue and Chikungunya negative samples were retrospectively screened to determine the presence of the Zika virus from different geographical regions of India.

Methods

During May to October 2021, the clinical samples of 1475 patients, across 13 states and a union territory of India were screened and re-tested for Dengue, Chikungunya and Zika by CDC Trioplex Real time RT-PCR. The Zika rRTPCR positive samples were further screened with anti-Zika IgM and Plaque Reduction Neutralization Test. Next generation sequencing was used for further molecular characterization.

Results

The positivity was observed for Zika (67), Dengue (121), and Chikungunya (10) amongst screened cases. The co-infections of Dengue/Chikungunya, Dengue/Zika, and Dengue/Chikungunya/Zika were also observed. All Zika cases were symptomatic with fever (84%) and rash (78%) as major presenting symptoms. Of them, four patients had respiratory distress, one presented with seizures, and one with suspected microcephaly at birth. The Asian Lineage of Zika and all four serotypes of Dengue were found in circulation.

Conclusion

Our study indicates the spread of the Zika virus to several states of India and an urgent need to strengthen its surveillance.

Booster dose of the inactivated COVID-19 vaccine BBV152 (Covaxin) enhances the neutralizing antibody response against Alpha, Beta, Delta and Omicron variants of concern

The neutralizing antibody responses were significantly elevated after 3rd dose of BBV152/Covaxin against the B.1 (19.11 fold) variant including Delta (16.51 fold), Beta (14.70 fold) and Omicron (18.53 fold) Variants of concern in comparison with two dose vaccine, providing assurance of a protective immune response of booster in recipients.

Immune responses against different variants of SARS-CoV-2 including Omicron following 6 months of administration of heterologous prime-boost COVID-19 vaccine

Comparative analysis at one- and six-months post-vaccination showed modest reduction in S1-RBD IgG antibody and NAb titers against B.1, Alpha, Beta and Delta variants in heterologous and homologous vaccine recipients groups. However, significant reduction in NAb titers against Omicron in vaccinees’ sera post-six months underlines need for cautious prospective follow-up.

SARS-CoV-2 Kappa Variant Shows Pathogenicity in a Syrian Hamster Model

Objectives: The emergence of SARS-CoV-2 lineage B.1.617 variants in India has been associated with a surge in the number of daily infections. We investigated the pathogenic potential of Kappa (B.1.617.1) variant in Syrian golden hamsters. Methods: Two groups of Syrian golden hamsters (18 each) were inoculated intranasally with SARS-CoV-2 isolates, B.1 (D614G) and Kappa variant, respectively. The animals were monitored daily for the clinical signs and body weight. Throat swab, nasal wash, and organ samples (lungs, nasal turbinate, trachea) were collected and screened using SARS-CoV-2-specific RT-qPCR. Histopathologic evaluation of the lung samples was performed. Results: The hamsters infected with the Kappa variant demonstrated increased body weight loss compared to the B.1 lineage isolate. The highest viral RNA load was observed in the nasal turbinate and lung specimens of animals infected with both variants. A significantly higher sgRNA load was observed in the nasal swabs (7 DPI), trachea (3 DPI), and lungs (3 DPI) of hamsters infected with the Kappa variant. Neutralizing antibody response generated in the B.1 lineage-infected hamster sera were comparable against both B.1 and Kappa variant in contrast to Kappa variant-infected hamsters, which showed lower titers against B.1 lineage isolate. Gross and microscopic evaluation of the lung specimens showed severe lung lesions in hamsters infected with Kappa variant compared to B.1. Conclusions: The study demonstrates pathogenicity of Kappa variant in hamsters evident with reduced body weight, high viral RNA load in lungs, and pronounced lung lesions. Both Kappa variant- and B.1-infected hamsters produced neutralizing antibodies against both variants studied.

Pathogenicity of SARS-CoV-2 Omicron (R346K) variant in Syrian hamsters and its cross-neutralization with different variants of concern

Background

SARS-CoV-2 Omicron variant is rampantly spreading across the globe. We assessed the pathogenicity and immune response generated by BA.1.1 sub-lineage of SARS-CoV-2 [Omicron (R346K) variant] in 5 to 6-week old Syrian hamsters and compared the observations with that of Delta variant infection.

Methods

Virus shedding, organ viral load, lung disease and immune response generated in hamsters were sequentially assessed.

Findings

The disease characteristics of the Omicron (R346K) variant were found to be similar to that of the Delta variant infection in hamsters like viral replication in the respiratory tract and interstitial pneumonia. The Omicron (R346K) infected hamsters demonstrated lesser body weight reduction and viral RNA load in the throat swab and nasal wash samples in comparison to the Delta variant infection. The viral load in the lungs and nasal turbinate samples and the lung disease severity of the Omicron (R346K) infected hamsters were found comparable with that of the Delta variant infected hamsters. Neutralizing antibody response against Omicron (R346K) variant was detected from day 5 and the cross-neutralization titre of the sera against other variants showed severe reduction ie., 7 fold reduction against Alpha and no titers against B.1, Beta and Delta.

Interpretation

This preliminary data shows that Omicron (R346K) variant infection can produce moderate to severe lung disease similar to that of the Delta variant and the neutralizing antibodies produced in response to Omicron (R346K) variant infection shows poor neutralizing ability against other co-circulating SARS-CoV-2 variants like Delta which necessitates caution as it may lead to increased cases of reinfection.

Funding

This study was supported by 10.13039/501100001411Indian Council of Medical Research as an intramural grant (COVID-19) to ICMR-National Institute of Virology, Pune.

Detection and isolation of SARS-CoV-2 Eta variant from the international travelers and local residents of India

International travel has been the major source for the rapid spread of new SARS‐CoV‐2 variants across the globe. During SARS‐CoV‐2 genomic surveillance, a total of 212 SARS‐CoV‐2 positive clinical specimens were sequenced using next‐generation sequencing. A complete SARS‐CoV‐2 genome could be retrieved from 90 clinical specimens. Of them, 14 sequences belonged to the Eta variant from clinical specimens of international travelers (n = 12) and local residents (n = 2) of India, and 76 belonged to other SARS‐CoV‐2 variants. Of all the Eta‐positive specimens, the virus isolates were obtained from the clinical specimens of six international travelers. Many variants of interest have been found to cause substantial community transmission or cluster infections. The detection of this variant with lethal E484K mutation across the globe and India necessitates persistent genomic surveillance of the SARS‐CoV‐2 variants, which would aid in taking preventive action.

Nipah Virus Outbreak in Kerala State, India Amidst of COVID-19 Pandemic

We report here a Nipah virus (NiV) outbreak in Kozhikode district of Kerala state, India, which had caused fatal encephalitis in a 12-year-old boy and the outbreak response, which led to the successful containment of the disease and the related investigations. Quantitative real-time reverse transcription (RT)-PCR, ELISA-based antibody detection, and whole genome sequencing (WGS) were performed to confirm the NiV infection. Contacts of the index case were traced and isolated based on risk categorization. Bats from the areas near the epicenter of the outbreak were sampled for throat swabs, rectal swabs, and blood samples for NiV screening by real-time RT-PCR and anti-NiV bat immunoglobulin G (IgG) ELISA. A plaque reduction neutralization test was performed for the detection of neutralizing antibodies. Nipah viral RNA could be detected from blood, bronchial wash, endotracheal (ET) secretion, and cerebrospinal fluid (CSF) and anti-NiV immunoglobulin M (IgM) antibodies from the serum sample of the index case. Rapid establishment of an onsite NiV diagnostic facility and contact tracing helped in quick containment of the outbreak. NiV sequences retrieved from the clinical specimen of the index case formed a sub-cluster with the earlier reported Nipah I genotype sequences from India with more than 95% similarity. Anti-NiV IgG positivity could be detected in 21% of Pteropus medius (P. medius) and 37.73% of Rousettus leschenaultia (R. leschenaultia). Neutralizing antibodies against NiV could be detected in P. medius. Stringent surveillance and awareness campaigns need to be implemented in the area to reduce human-bat interactions and minimize spillover events, which can lead to sporadic outbreaks of NiV.

Genome Sequencing Reveals a Mixed Picture of SARS-CoV-2 Variant of Concern Circulation in Eastern Uttar Pradesh, India

Uttar Pradesh is the densely populated state of India and is the sixth highest COVID-19 affected state with 22,904 deaths recorded on November 12, 2021. Whole-genome sequencing (WGS) is being used as a potential approach to investigate genomic evolution of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus. In this study, a total of 87 SARS-CoV-2 genomes-49 genomes from the first wave (March 2020 to February 2021) and 38 genomes from the second wave (March 2021 to July 2021) from Eastern Uttar Pradesh (E-UP) were sequenced and analyzed to understand its evolutionary pattern and variants against publicaly available sequences. The complete genome analysis of SARS-CoV-2 during the first wave in E-UP largely reported transmission of G, GR, and GH clades with specific mutations. In contrast, variants of concerns (VOCs) such as Delta (71.0%) followed by Delta AY.1 (21.05%) and Kappa (7.9%) lineages belong to G clade with prominent signature amino acids were introduced in the second wave. Signature substitution at positions S:L452R, S:P681R, and S:D614G were commonly detected in the Delta, Delta AY.1, and Kappa variants whereas S:T19R and S:T478K were confined to Delta and Delta AY.1 variants only. Vaccine breakthrough infections showed unique mutational changes at position S:D574Y in the case of the Delta variant, whereas position S:T95 was conserved among Kappa variants compared to the Wuhan isolate. During the transition from the first to second waves, a shift in the predominant clade from GH to G clade was observed. The identified spike protein mutations in the SARS-CoV-2 genome could be used as the potential target for vaccine and drug development to combat the effects of the COVID-19 disease.

Immunogenicity and safety of a heterologous prime-boost COVID-19 vaccine schedule: ChAdOx1 vaccine Covishield followed by BBV152 Covaxin

The evidence for effectiveness of heterologous priming of COVID-19 vaccine is very limited. Here, we studied eighteen participants who received heterologous vaccination regimen of AstraZeneca’s ChAdOx1-nCov-19 followed by inactivated whole virion BBV152. Heterologous group participant doesn’t report any adverse event following immunization and demonstrated high humoral and neutralizing antibody response.

Neutralization of Beta and Delta variant with sera of COVID-19 recovered cases and vaccinees of inactivated COVID-19 vaccine BBV152/Covaxin

SARS-CoV-2 variant of concern mainly Delta and Beta is global public health concern due to its impact on existing vaccines. Here, we assessed the neutralization of sera from COVID-19 recovered cases and BBV152 vaccinees against Beta and Delta variants. BBV152 vaccine found to confer significant protection against both the variants.

Identification of Phasi Charoen-Like Phasivirus in Field Collected Aedes aegypti from Karnataka State, India

Background: A wide range of insect-specific viruses (ISVs) have been reported worldwide. There are no studies from India that have reported ISVs. The current study describes the identification of Phasi Charoen-like virus (PCLV) from Aedes aegypti mosquito-pools from six districts of Karnataka state, India. Materials and Methods: During the Chikungunya virus (CHIKV) outbreak in the Bangalore Urban district in 2019, using conventional PCR, it was found that both human and mosquito samples were positive for CHIKV. For retrieve the complete genome sequence, mosquito samples were subjected to next generation sequencing (NGS) analysis and PCLV was also found. During 2019, as part of a vector-borne disease surveillance, we received 50 mosquito pool samples from 6 districts of the state, all of them were subjected to NGS to identify PCLV. Results: The A. aegypti mosquito-pools samples were subjected to the NGS platform that led to identification of an ISV, PCLV. PCLV was identified in 26 A. aegypti mosquito-pools collected from 6 districts. We also found mixed infection of PCLV with the Dengue virus (DENV; genotypes 1 and 3) and CHIKV from five pools. The nucleotide identity for the L gene of Indian PCLV sequences ranged between 97.1% and 98.3% in comparison with the Thailand sequences. Conclusions: To the best of our knowledge, this is the first report of PCLV dual infection with DENV and CHIKV in India. The present study confirms the presence of PCLV in A. aegypti mosquitoes from Karnataka state. The study adds India in the global geographical distribution of PCLV.

Clinical Characterization and Genomic Analysis of Samples from COVID-19 Breakthrough Infections during the Second Wave among the Various States of India

From March to June 2021, India experienced a deadly second wave of COVID-19, with an increased number of post-vaccination breakthrough infections reported across the country. To understand the possible reason for these breakthroughs, we collected 677 clinical samples (throat swab/nasal swabs) of individuals from 17 states/Union Territories of the country who had received two doses (n = 592) and one dose (n = 85) of vaccines and tested positive for COVID-19. These cases were telephonically interviewed and clinical data were analyzed. A total of 511 SARS-CoV-2 genomes were recovered with genome coverage of higher than 98% from both groups. Analysis of both groups determined that 86.69% (n = 443) of them belonged to the Delta variant, along with Alpha, Kappa, Delta AY.1, and Delta AY.2. The Delta variant clustered into four distinct sub-lineages. Sub-lineage I had mutations in ORF1ab A1306S, P2046L, P2287S, V2930L, T3255I, T3446A, G5063S, P5401L, and A6319V, and in N G215C; Sub-lineage II had mutations in ORF1ab P309L, A3209V, V3718A, G5063S, P5401L, and ORF7a L116F; Sub-lineage III had mutations in ORF1ab A3209V, V3718A, T3750I, G5063S, and P5401L and in spike A222V; Sub-lineage IV had mutations in ORF1ab P309L, D2980N, and F3138S and spike K77T. This study indicates that majority of the breakthrough COVID-19 clinical cases were infected with the Delta variant, and only 9.8% cases required hospitalization, while fatality was observed in only 0.4% cases. This clearly suggests that the vaccination does provide reduction in hospital admission and mortality.

Clinical Characterization and Genomic Analysis of Samples from COVID-19 Breakthrough Infections during the Second Wave among the Various States of India

From March to June 2021, India experienced a deadly second wave of COVID-19, with an increased number of post-vaccination breakthrough infections reported across the country. To understand the possible reason for these breakthroughs, we collected 677 clinical samples (throat swab/nasal swabs) of individuals from 17 states/Union Territories of the country who had received two doses (n = 592) and one dose (n = 85) of vaccines and tested positive for COVID-19. These cases were telephonically interviewed and clinical data were analyzed. A total of 511 SARS-CoV-2 genomes were recovered with genome coverage of higher than 98% from both groups. Analysis of both groups determined that 86.69% (n = 443) of them belonged to the Delta variant, along with Alpha, Kappa, Delta AY.1, and Delta AY.2. The Delta variant clustered into four distinct sub-lineages. Sub-lineage I had mutations in ORF1ab A1306S, P2046L, P2287S, V2930L, T3255I, T3446A, G5063S, P5401L, and A6319V, and in N G215C; Sub-lineage II had mutations in ORF1ab P309L, A3209V, V3718A, G5063S, P5401L, and ORF7a L116F; Sub-lineage III had mutations in ORF1ab A3209V, V3718A, T3750I, G5063S, and P5401L and in spike A222V; Sub-lineage IV had mutations in ORF1ab P309L, D2980N, and F3138S and spike K77T. This study indicates that majority of the breakthrough COVID-19 clinical cases were infected with the Delta variant, and only 9.8% cases required hospitalization, while fatality was observed in only 0.4% cases. This clearly suggests that the vaccination does provide reduction in hospital admission and mortality.

Sequential determination of viral load, humoral responses and phylogenetic analysis in fatal and non-fatal cases of Crimean-Congo hemorrhagic fever patients from Gujarat, India, 2019

Background

Thirty-four CCHF cases (17 fatal; 17 survived) were confirmed from Gujarat state, India during the year 2019. We aimed to find out the viral load, antibody kinetics, cytokine profile and phylogenetic analysis between fatal and non- fatal cases.

Methods

Thirty four cases were included in this study. Blood and urine samples were collected from all the cases on the day of admission to hospital. Non-fatal cases were followed weekly for understanding the profile of viral kinetics, anti-CCHFV IgM and IgG antibodies. We also quantified the cytokines in both fatal and non-fatal cases. For epidemiological correlation, livestock were screened for anti-CCHF IgG antibodies and the tick pool specimens were tested by real time RT-PCR. Virus isolation was attempted on tick pools and human specimens and phylogenetic analysis performed on human and ticks complete genome sequences.

Results

CCHF cases were detected throughout year in 2019 with the peak in August. Out of 34 cases, eight secondary CCHF cases were reported. Cases were predominantly detected in males and in 19–45 years age group (55.88%). The persistence of viremia was observed till 76^th POD (post onset date) in one case whereas anti-CCHFV IgM and IgG was detected amongst these cases from the 2^nd and 20^th POD respectively. Positivity observed amongst livestock and tick pools were was 21.57% and 7.4% respectively. The cytokine analysis revealed a significant increase in the level of serum IL-6, IL-10 and IFN-γ during the acute phase of the infection, but interestingly IL-10 lowered to normal upon clearance of the virus in the clinically recovered case. Fatal cases had high viral RNA copy numbers. Bleeding from one or two mucosal sites was significantly associated with fatality (OR-16.47;p-0.0034 at 95% CI). We could do CCHF virus isolation from two cases. Phylogenetic analysis revealed circulation of re-assortment of Asian-West African genotypes in humans and ticks.

Conclusions

The persistence of CCHF viral RNA was detected till 76^th POD in one of the survivors. The circulation of a re-assortment Asian-West African genotype in a CCHF case is also reported first time from India.

Crimean Congo hemorrhagic fever is a zoonotic tick-borne viral hemorrhagic disease. This disease is reported from Europe, Mediterranean, north-western China, central Asia, Africa, and the Middle East. Several outbreaks of CCHF were reported from Gujarat and Rajasthan states, India from 2011 to 2019. In this study, we discuss the clinical, molecular, serological, and the cytokine data of 34 CCHF cases (17 fatal and 17 survived) which were detected from Gujarat state in the year 2019. A sequential weekly follow up of the CCHF survivors was performed to understand the viral kinetics and the antibody profile. Interestingly, the presence of persistence CCHF viral RNA was observed till 76^th POD in one of the survivors. To our knowledge, we are reporting this long term persistence of viremia for the first time. We also observed that the anti-CCHFV IgM detection in the serum samples starts as soon as 2^nd POD but anti-CCHFV IgG antibody could be detected in the majority of the cases only after the 28^th POD. The cytokine analysis revealed a significant increase in the level of serum IL-6, IL-10 and IFN-γ during the acute phase of the infection, but interestingly IL-10 lowered to normal upon clearance of the virus in the clinically recovered case. We did the phylogenetic analysis and concluded the circulation of the Asian-West African re-assortment genotype in humans, which has not been reported from India prior to this study.

An Epidemiological Analysis of SARS-CoV-2 Genomic Sequences from Different Regions of India

The number of Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) cases is increasing in India. This study looks upon the geographic distribution of the virus clades and variants circulating in different parts of India between January and August 2020. The NPS/OPS from representative positive cases from different states and union territories in India were collected every month through the VRDLs in the country and analyzed using next-generation sequencing. Epidemiological analysis of the 689 SARS-CoV-2 clinical samples revealed GH and GR to be the predominant clades circulating in different states in India. The northern part of India largely reported the ‘GH’ clade, whereas the southern part reported the ‘GR’, with a few exceptions. These sequences also revealed the presence of single independent mutations—E484Q and N440K—from Maharashtra (first observed in March 2020) and Southern Indian States (first observed in May 2020), respectively. Furthermore, this study indicates that the SARS-CoV-2 variant (VOC, VUI, variant of high consequence and double mutant) was not observed during the early phase of virus transmission (January–August). This increased number of variations observed within a short timeframe across the globe suggests virus evolution, which can be a step towards enhanced host adaptation.

Isolation and characterization of the new SARS-CoV-2 variant in travellers from the United Kingdom to India: VUI-202012/01 of the B.1.1.7 lineage

We have isolated the new severe acute respiratory syndrome coronavirus-2 variant of concern 202 012/01 from the positive coronavirus disease 2019 cases that travelled from the UK to India in the month of December 2020. This emphasizes the need for the strengthened surveillance system to limit the local transmission of this new variant.

Detection of Nipah virus in Pteropus medius in 2019 outbreak from Ernakulam district, Kerala, India

Background

In June 2019, Nipah virus (NiV) infection was detected in a 21-year-old male (index case) of Ernakulum, Kerala, India. This study was undertaken to determine if NiV was in circulation in Pteropus species (spp) in those areas where the index case had visit history in 1 month.

Methods

Specialized techniques were used to trap the Pteropus medius bats (random sampling) in the vicinity of the index case area. Throat and rectal swabs samples of 141 bats along with visceral organs of 92 bats were collected to detect the presence of NiV by real-time reverse transcriptase-polymerase chain reaction (qRTPCR). Serum samples of 52 bats were tested for anti-NiV Immunoglobulin (Ig) G antibodies by Enzyme-Linked Immunosorbent Assay (ELISA). The complete genome of NiV was sequenced by next-generation sequencing (NGS) from the tissues and swab samples of bats.

Results

One rectal swab sample and three bats visceral organs were found positive for the NiV. Interestingly, 20.68% (12/58) of Pteropus were positive for anti-NiV IgG antibodies. NiV sequences of 18,172; 17,200 and 15,100 nucleotide bps could be retrieved from three Pteropus bats.

Conclusion

A distinct cluster of NiV sequences, with significant net-evolutionary nucleotide divergence, was obtained, suggesting the circulation of new genotype (I-India) in South India. NiV Positivity in Pteropus spp. of bats revealed that NiV is circulating in many districts of Kerala state, and active surveillance of NiV should be immediately set up to know the hotspot area for NiV infection.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-021-05865-7.

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.

Transcriptome & viral growth analysis of SARS-CoV-2-infected Vero CCL-81 cells

Background & objectives:

The genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), belonging to the family Coronaviridae, encodes for structural, non-structural, and accessory proteins, which are required for replication of the virus. These proteins are encoded by different genes present on the SARS-CoV-2 genome. The expression pattern of these genes in the host cells needs to be assessed. This study was undertaken to understand the transcription pattern of the SARS-CoV-2 genes in the Vero CCL-81 cells during the course of infection.

Methods:

Vero CCL-81 cells were infected with the SARS-CoV-2 virus inoculum having a 0.1 multiplicity of infection. The supernatants and cell pellets were harvested after centrifugation at different time points, post-infection. The 50% tissue culture infective dose (TCID₅₀) and cycle threshold (C_t) values of the E and the RdRp-2 genes were calculated. Next-generation sequencing of the harvested sample was carried out to observe the expression pattern of the virus by mapping to the SARS-CoV-2 Wuhan HU-1 reference sequence. The expressions were in terms of the reads per kilobase million (RPKM) values.

Results:

In the inital six hours post-infection, the copy numbers of E and RdRp-2 genes were approximately constant, which raised 10 log-fold and continued to increase till the 12 h post-infection (hpi). The TCID₅₀ was observed in the supernatant after 7 hpi, indicating the release of the viral progeny. ORF8 and ORF7a, along with the nucleocapsid transcript, were found to express at higher levels.

Interpretation & conclusions:

This study was a step towards understanding the growth kinetics of the SARS-CoV-2 replication cycle. The findings indicated that ORF8 and ORF7b gene transcripts were expressed in higher amounts indicating their essential role in viral replication. Future studies need to be conducted to explore their role in the SARS-CoV-2 replication.

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.

Detection of coronaviruses in Pteropus & Rousettus species of bats from different States of India

Background & objectives:

Bats are considered to be the natural reservoir for many viruses, of which some are potential human pathogens. In India, an association of Pteropus medius bats with the Nipah virus was reported in the past. It is suspected that the recently emerged severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) also has its association with bats. To assess the presence of CoVs in bats, we performed identification and characterization of bat CoV (BtCoV) in P. medius and Rousettus species from representative States in India, collected during 2018 and 2019.

Methods:

Representative rectal swab (RS) and throat swab specimens of Pteropus and Rousettus spp. bats were screened for CoVs using a pan-CoV reverse transcription-polymerase chain reaction (RT-PCR) targeting the RNA-dependent RNA polymerase (RdRp) gene. A single-step RT-PCR was performed on the RNA extracted from the bat specimens. Next-generation sequencing (NGS) was performed on a few representative bat specimens that were tested positive. Phylogenetic analysis was carried out on the partial sequences of RdRp gene sequences retrieved from both the bat species and complete viral genomes recovered from Rousettus spp.

Results:

Bat samples from the seven States were screened, and the RS specimens of eight Rousettus spp. and 21 Pteropus spp. were found positive for CoV RdRp gene. Among these, by Sanger sequencing, partial RdRp sequences could be retrieved from three Rousettus and eight Pteropus bat specimens. Phylogenetic analysis of the partial RdRp region demonstrated distinct subclustering of the BtCoV sequences retrieved from these Rousettus and Pteropus spp. bats. NGS led to the recovery of four sequences covering approximately 94.3 per cent of the whole genome of the BtCoVs from Rousettus bats. Three BtCoV sequences had 93.69 per cent identity to CoV BtRt-BetaCoV/GX2018. The fourth BtCoV sequence was 96.8 per cent identical to BtCoV HKU9-1.

Interpretation & conclusions:

This study was a step towards understanding the CoV circulation in Indian bats. Detection of potentially pathogenic CoVs in Indian bats stresses the need for enhanced screening for novel viruses in them. One Health approach with collaborative activities by the animal health and human health sectors in these surveillance activities shall be of use to public health. This would help in the development of diagnostic assays for novel viruses with outbreak potential and be useful in disease interventions. Proactive surveillance remains crucial for identifying the emerging novel viruses with epidemic potential and measures for risk mitigation.

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.

Kinetics of viral RNA, immunoglobulin-M & G antibodies in Kyasanur forest disease

Background & objectives:

Kyasanur forest disease (KFD) is an infectious disease discovered in Karnataka State of India in 1957; since then, the State has been known to be enzootic for KFD. In the last few years, its presence was observed in the adjoining five States of the Western Ghats of India. The present study was conducted to understand the kinetics of viral RNA, immunoglobulin M (IgM) and IgG antibody in KFD-infected humans for developing a diagnostic algorithm for KFD.

Methods:

A prospective follow up study was performed among KFD patients in Sindhudurg district of Maharashtra State, India. A total of 1046 suspected patients were tested, and 72 KFD patients were enrolled and followed for 17 months (January 2016 to May 2017). Serum samples of KFD patients were screened for viral RNA, and IgM and IgG antibodies.

Results:

KFD viral positivity was observed from 1^st to 18^th post-onset day (POD). Positivity of anti-KFD virus (KFDV) IgM antibodies was detected from 4^th till 122^nd POD and anti-KFDV IgG antibodies detected from 5^th till 474^th POD. A prediction probability was determined from statistical analysis using the generalized additive model in R-software to support the laboratory findings regarding viral kinetics.

Interpretation & conclusions:

This study demonstrated the presence of KFD viral RNA till 18^th POD, IgM antibodies till 122^nd POD and IgG till the last sample collected. Based on our study an algorithm was recommended for accurate laboratory diagnosis of KFDV infection. A sample collected between 1 and 3 POD can be tested using KFDV real-time reverse transcriptase polymerase chain reaction (RT-PCR); between 4 and 24 POD, the combination of real-time RT-PCR and anti-KFDV IgM enzyme-linked immunosorbent assay (ELISA) tests can be used; between POD 25 and 132, anti-KFDV IgM and IgG ELISA are recommended.

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.

Isolation and phylogenomic analysis of buffalopox virus from human and buffaloes in India

Three genome sequences of Buffalopox virus (BPVX) were retrieved from a human and two buffaloes scab samples. Phylogenomic analysis of the BPXV indicates that it shares a most recent common ancestor with Lister and closely related vaccine strains when compared to potential wild-type VACV strains (like Horsepox virus).

Identification and phylogenetic analysis of herpes simplex virus-1 from clinical isolates in India

Human herpes simplex virus (HSV)-1 infection is acquired in childhood and persists throughout a person’s lifetime. Here, we present two cases from India; one showing symptoms of postpartum haemorrhage with disseminated intravascular coagulation, and the second one showing signs of acute encephalitis syndrome. The aetiological agent in both cases was identified as HSV-1 using the PCR method. The next-generation sequencing method retrieved ~97 % of the viral genome from the isolates of the clinical samples. The phylogenetic analysis of the retrieved genomes revealed that they belong to clade II of HSV-1. This study identifies a few sequence variations in the glycoprotein region of HSV-1 during two different clinical manifestations. There are a couple of papers that analyse variations in the glycoprotein region of clinical samples. Further, this study also highlights the importance of considering HSV-1 during differential diagnosis when analysing the nosocomial infection.

Characterization of a strain of quaranfil virus isolated from soft ticks in India. Is quaranfil virus an unrecognized cause of disease in human and animals?"

The soft ticks collected during a field survey in Karnataka state, India, in 1983, yielded a novel virus isolate, which caused mortality in an infant mouse upon inoculation. Attempts at characterizing the virus using the conventional methods were unsuccessful, which prompted us to study it by Next-Generation Sequencing (NGS). This virus isolate was obtained from the viral repository of National Institute of Virology, and an initial virus stock was prepared as a mouse brain homogenate. The virus stock showed cytopathic effects in different cell-lines and was used in NGS. Based on the complete genome sequence, obtained using de novo and reference mapping approach, the virus isolate was identified as a Quaranfil virus (QRFV) belonging to the family Orthomyxoviridae, genus Quaranjavirus. The genome size of the virus is 11,427 nucleotides which consist of 6 segments encoding six proteins. Homology analysis suggested this isolate as similar to QRFV of Afghanistan. In silico analysis showed the HA protein secondary structure to be a class III penetrance similar to Thogotovirus. QRFV was first isolated in 1953 from ticks [Cairo, Egypt] and subsequently reported from other geographical areas. This is the first report describing the presence of QRFV from India. This discovery emphasizes the need for investigating mild febrile illness cases with influenza-like symptoms, particularly in the area of high risk for tick bites.

Equine Encephalosis Virus in India, 2008

A virus isolated from a sick horse from India in 2008 was confirmed by next-generation sequencing analysis to be equine encephalosis virus (EEV). EEV in India is concerning because several species of Culicoides midge, which play a major role in EEV natural maintenance and transmission, are present in this country.

Complete genome sequencing of Kaisodi virus isolated from ticks in India belonging to Phlebovirus genus, family Phenuiviridae

An unknown virus was repeatedly isolated from hard tick (Haemaphysalis spinigera) during a proactive arbovirus survey in ticks conducted in 1957, in India. The virus remained uncharacterized for a long time. The passages of this virus in different vertebrate and invertebrate cells along with human and monkey-derived cell culture showed no cytopathic effect. It was identified later to be a member of Kaisodi group among Phlebovirus genus in the family Phenuiviridae (Order: Bunyavirales) by serological methods. Due to its genomic diversity, sequencing of this virus was a challenge for a while. In this study, we were able to sequence the complete genome of this virus isolate using next-generation sequencing (NGS) platform. The unknown virus was identified to be Kaisodi virus (KASDV) using NGS analysis. De novo genome assembly derived three genomic segments for the KASDV which encode for RNA-dependent RNA polymerase, glycoprotein precursor, and nucleoprotein. Functional as well as conserved domains for Kaisodi serogroup viruses were predicted and compared to a known representative of the genus Phlebovirus. The phylogenetic tree revealed its closeness to Silverwater virus, of Kaisodi serogroup with nucleotide (69%, 62%, and 61%) and amino acid (52%, 51%, and 62%) identity for L, M, and S segment, respectively. The study demonstrates the presence of a conserved motif (⁷²TRGNK⁷⁶) around the RNA binding motif region in tick-borne phleboviruses. The intergenic region encompassing the S segment of Kaisodi serogroup was GC-rich whereas the other Phlebovirus had AT-rich genome. KASDV has the largest intergenic region and larger loops, suggesting stem-loops formed due to larger loops as a possible factor for instability and cause of transcription termination. This paper also describes the real-time RT-PCR and RT-PCR assays developed and used for the detection of KASDV RNA in ticks from Karnataka, Kerala and Maharashtra State, India. The KASDV positivity observed in the recently collected tick pools indicates that the KASDV, isolated from Karnataka state in 1957, is also circulating in the adjoining Kerala state. On the basis of the current study, it should be possible to develop diagnostic assays which would facilitate an in-depth field survey exploring the veterinary and medical significance of KASDV.

Clinico-epidemiological investigation on Varicella Zoster Virus indicates multiple clade circulation in Maharashtra state, India

BACKGROUND

Varicella Zoster Virus (VZV) is consistently in circulation and shows an increase in disease burden during the spring season. Due to a wide range of clinical presentation from a vesicular rash to bleeding or neurological complications, it makes the clinical diagnosis difficult. The present study aims to understand whether the same strain of virus is responsible for the increase in the seasonal outbreaks occurring in different parts of the country with reference to the samples from Maharashtra, Rajasthan and Gujarat states of India.

MATERIALS AND METHODS

This study reports the clinico-epidemiological and laboratory findings of suspected Varicella cases. To understand the circulating clade few representative real-time Polymerase Chain Reaction (PCR) positive were analyzed by conventional PCR and partial Open Reading Frame (ORF) 22, partial ORF 38 and partial ORF 54 were sequenced to identify single nucleotide polymorphisms responsible for clade determination. Further partial glycoprotein B gene was sequenced, and a phylogenetic tree was generated.

RESULTS

A total of 50 cases from Maharashtra (Mumbai district) and referred clinical samples of Rajasthan (Barmer district; n = 12) and Gujarat States (Gandhi Nagar, Surat districts; n = 17) were tested for the presence of VZV. Vesicular rash with fever was a common clinical presentation with 82% cases having contact history with VZV positive cases, suggesting higher secondary attack rate. The vesicular fluid of all 50 cases from Mumbai revealed the presence of VZV by real-time PCR. Urine, serum and throat swab samples showed positivity by real-time PCR. Healthcare provider's samples from Rajasthan showed 36.4% [4/11] positivity. Clinical samples from Gujarat had positivity of 41.2% [7/17].

CONCLUSIONS

This study analyses the clade based circulation of VZV in three states in India and suggests different clades circulating in Maharashtra state. Health education amongst the general population is suggested to reduce the secondary cases by early diagnosis, effective isolation policies and vaccination to reduce the burden of disease.

Identification and characterization of novel mosquito-borne (Kammavanpettai virus) and tick-borne (Wad Medani) reoviruses isolated in India

In 1954, a virus named Wad Medani virus (WMV) was isolated from Hyalomma marginatum ticks from Maharashtra State, India. In 1963, another virus was isolated from Sturnia pagodarum birds in Tamil Nadu, India, and named Kammavanpettai virus (KVPTV) based on the site of its isolation. Originally these virus isolates could not be identified with conventional methods. Here we describe next-generation sequencing studies leading to the determination of their complete genome sequences, and identification of both virus isolates as orbiviruses (family Reoviridae). Sequencing data showed that KVPTV has an AT-rich genome, whereas the genome of WMV is GC-rich. The size of the KVPTV genome is 18 234 nucleotides encoding proteins ranging 238-1290 amino acids (aa) in length. Similarly, the size of the WMV genome is 16 941 nucleotides encoding proteins ranging 214-1305 amino acids in length. Phylogenetic analysis of the VP1 gene, along with the capsid genes VP5 and VP7, revealed that KVPTV is likely a novel mosquito-borne virus and WMV is a tick-borne orbivirus. This study focuses on the phylogenetic comparison of these newly identified orbiviruses with mosquito-, tick- and Culicoides-borne orbiviruses isolated in India and other countries.