key: cord-0758618-xmgi40hg
authors: Dhar, Mahesh S; Marwal, Robin; Radhakrishnan, VS; Ponnusamy, Kalaiarasan; Jolly, Bani; Bhoyar, Rahul C.; Sardana, Viren; Naushin, Salwa; Rophina, Mercy; Mellan, Thomas A; Mishra, Swapnil; Whittaker, Charles; Fatihi, Saman; Datta, Meena; Singh, Priyanka; Sharma, Uma; Ujjainiya, Rajat; Bhatheja, Nitin; Divakar, Mohit Kumar; Singh, Manoj K; Imran, Mohamed; Senthivel, Vigneshwar; Maurya, Ranjeet; Jha, Neha; Mehta, Priyanka; Vivekanand, A; Sharma, Pooja; Arvinden, VR; Chaudhary, Urmila; Soni, Namita; Thukral, Lipi; Flaxman, Seth; Bhatt, Samir; Pandey, Rajesh; Dash, Debasis; Faruq, Mohammed; Lall, Hemlata; Gogia, Hema; Madan, Preeti; Kulkarni, Sanket; Chauhan, Himanshu; Sengupta, Shantanu; Kabra, Sandhya; Gupta, Ravindra K.; Singh, Sujeet K; Agrawal, Anurag; Rakshit, Partha
title: Genomic characterization and Epidemiology of an emerging SARS-CoV-2 variant in Delhi, India
date: 2021-11-19
journal: Science
DOI: 10.1126/science.abj9932
sha: 78fddad1e95f5eaeb02bce2f127a5299bb57b88b
doc_id: 758618
cord_uid: xmgi40hg

Delhi, the national capital of India, has experienced multiple SARS-CoV-2 outbreaks in 2020 and reached population seropositivity of over 50% by 2021. During April 2021, the city became overwhelmed by COVID-19 cases and fatalities, as a new variant B.1.617.2 (Delta) replaced B.1.1.7 (Alpha). A Bayesian model explains the growth advantage of Delta through a combination of increased transmissibility and reduced sensitivity to immune responses generated against earlier variants (median estimates; ×1.5-fold, 20% reduction). Seropositivity of an employee and family cohort increased from 42% to 87.5% between March and July 2021, with 27% reinfections, as judged by increased antibody concentration after a previous decline. The likely high transmissibility and partial evasion of immunity by the Delta variant contributed to an overwhelming surge in Delhi.

After escaping relatively unscathed during the first wave of the COVID-19 pandemic, India witnessed a ferocious second COVID-19 wave, starting in March 2021 and accounting for about half of global cases by the first week of May. SARS-CoV-2 had spread widely throughout India in the first wave, with the third national serosurvey in January 2021 finding that 21.4% of adults and 25.3% of 10-17 year old adolescents were seropositive (1) . Delhi, the national capital, was not included in the national serosurvey but had undergone multiple periods of high transmission in 2020 (Fig 1A) . In a district-wise stratified serosurvey conducted by the Delhi Government in January 2021, overall seropositivity was reported to be 56.1 % (95% CI, 55.5 -56.8%), ranging from 49.1% to 62.2% across eleven districts (2) . This was expected to confer some protection from future outbreaks.

Despite high seropositivity, Delhi was amongst the most affected cities during the second wave. The rise in new cases was exceptionally rapid in April, going from approximately 2000 to 20000 between March 31 and April 16. This was accompanied by a rapid rise in hospitalizations and ICU admissions (Fig 1B) . In this emergency situation with saturated bed occupancy by April 12, major private hospitals were declared by the state as full COVID care-only and senior medical students, including from alternative medicine branches, were pressed into service (3). Deaths rose proportionately ( Fig 1C) and the case-fatality ratio (CFR), estimated as the scaling factor between time-advanced cases and deaths (Fig 1D) , was stable (mean, SD; 1.9, 0.3%). Population spread of SARS-CoV-2 is underestimated by test positive cases alone (1, 2) . To better understand the degree of spread and the factors leading to the unexpectedly severe outbreak, we used all available data including testing, sequencing, serosurveys, and serially followed cohorts.

In the absence of finely resolved or serial data from national and state surveys, we focussed on data for Delhi participants of a national serosurvey of Council of Scientific and Industrial Research (CSIR, India) employees and their family members (Fig 2A, table S1 ). Samples were initially collected from the end of July to mid-September 2020 (Phase I). Subsequently 2 nd and 3 rd surveys were held in January/February (Phase II) and end of May to early July 2021 (Phase III), bracketing the time period of interest. The cohort details and serosurvey methodology have been previously published (4).

Infection was determined by anti-nucleocapsid assay, which is not affected by immunization with spike protein based vaccines. The presence of neutralizing antibodies (NAb) to wild-type SARS-CoV-2 spike protein was estimated by a surrogate viral neutralization test (sVNT, Genscript). Previous results from the full cohort have been comparable to government serosurveys, but Delhi cohort values have been slightly lower. This may be due to an overrepresentation of members with the ability to reduce exposure and avoid public transport (Fig 2A) . Within these limitations, the Delhi cohort showed a rise in seropositivity from 14.7% in Phase I (95% CI, 12.6-17.0%) to 42.1% in Phase II (95% CI, 39.0-45.2%). About one-third of NAb +ve subjects at Phase I became NAb -ve by Phase II, with most showing declining inhibition on sVNT assays ( Fig 2B) (5) . Phase III seropositivity increased to 87.5% (95% CI, 85.0-89.7%) amongst unvaccinated subjects. New infections between March and June 2021 are thus likely to have vastly exceeded known cases. Amongst ninety one previously infected subjects with serial measurements at three time points including Phase III (T3), 25 (27.5%, 95% CI 18.4-37.5%) had a pattern of declining antibody concentration between T1 and T2, followed by a sharp rise at T3, indicative of reinfection ( Fig 2C) . Confirmation of reinfection by either RT-PCR (n=8) or symptomatic illness (n=2) was available for ten of the twenty-five subjects. No severe illness or hospitalization was reported in reinfections.

Time periods of increased transmission were associated with declining RT-PCR cycle threshold (Ct) values (Fig 3A, fig. S1 ), attributable to a higher proportion of recently infected individuals with high viral loads (6) . However, the Ct decline was far greater in April 2021 (dCt, SE; -4.06, 0.27, p<0.001) than seen previously (dCt Nov 2020 ~ -1.5). Comparing April 2021 with November 2020, high viral load samples (Ct<20) doubled in clinical samples (p<0.001), and nearly doubled in campus surveillance testing data, where most positives were from recently infected individuals (15% (n=297) vs 9% (n=358), p=0.02).

Genome analysis trends, in representative samples drawn from the general population over the same period, showed seeding and expansion of B.1.1.7 (Alpha), B.1.617.1 (Kappa) and B.1.617.2 (Delta) lineages in 2021, with Delta becoming the dominant lineage in Delhi during April ( Fig 3B) . The proportion of Delta variant was strongly correlated to the rise in cases and healthcare stress ( fig. S2 ). Overall, the genomic and epidemiological data were most consistent with the hypothesis that a new variant with higher infectivity, Delta, was driving the unexpected overwhelming surge in Delhi. Recent in vitro data supports the possibility of a higher replication rate for Delta, thereby explaining potentially higher viral loads in RT-PCR data and greater transmissibility (7).

We further investigated the sequence of seeding and spread of Alpha, Kappa, and Delta variants of SARS-CoV-2. Phylogenetic analysis showed common origins between Alpha variants in Delhi and Punjab, and between Kappa or Delta variants in Delhi and Maharashtra, where Kappa and Delta were first sequenced (Fig 3C, fig. S3 ). There was a recurring pattern of initial smaller outbreaks with Alpha variant, followed by larger outbreaks coinciding with Alpha to Delta transition across all neighbouring states ( fig. S4 ). The substantial relative growth advantage of Delta was explored in terms of transmissibility and/or immune escape. The rise of Delta, but not other lineages, was temporally coincident S5 ). While overall vaccination levels were only about 5% in Delhi, most healthcare workers had received one or two doses of ChAdOx1-nCov19 (Astra-Zeneca / Serum Institute, India) or BBV152 (Bharat Biotech, India) (8, 9) . We sequenced twenty-four breakthrough infections starting at least one week after the first dose, collected between 22 March and 28 April at NCDC. The ratio of Delta to non-Delta lineages was 850:1211 from 20 March to 30 April.

In contrast, Delta to other lineage ratio was 13:3 in sixteen breakthroughs post one dose and 7:1 in eight breakthroughs post the 2 nd dose of vaccine. While the small sample size and lack of a formal control group preclude definitive analysis, estimated higher odds for Delta in vaccination breakthroughs (OR, 7.1; 95% CI, 2.4-20.9) corroborate other reports of reduced vaccine effectiveness against Delta (10).

To better characterise how the properties of Delta might differ from other SARS-CoV-2 lineages previously circulating in the city, we used a Bayesian model of SARS-CoV-2 transmission and mortality that simultaneously models the dynamics of two-categories of virus ("B.1.617.2" and "non-B.1.617.2") (11), whilst also explicitly incorporating natural waning of immunity derived from prior infection, with the duration of immunity consistent with the results of recent longitudinal cohort studies (12, 13) . Details of the model and input data are in Supplementary Methods. Briefly, the model is fitted to COVID-19 mortality data, genomic sequence data presented here and from GISAID (with Pangolin lineage classification) (14, 15, 16) and serological data presented alongside an additional longitudinal serosurvey carried out in the city in the period July-December 2020 (17) . Substantial uncertainties remain as to the date of introduction of B.1.617.2 into Delhi and the degree of COVID-19 death under-ascertainment. We, therefore, explored a range of different scenarios in which we varied under-ascertainment (10%, 33%, 50%, and 66%) and introduction dates (15 January 2021, 31 January 2021, 14 February 2021 and 28 February 2021).

Using this framework for an introduction date of 14 February 2021 and mortality under ascertainment of 50%, our results shown in Fig. 4 indicate that the Delta variant is ×1.3-1.7 (50% bCI) fold more transmissible than earlier/co-circulating SARS-CoV-2 lineages in Delhi, including the Alpha variant. Importantly, the model also indicates that the Delta variant can partially evade immunity elicited by prior infection, with prior infection providing only 50 to 90% (50% bCI) of the protection against infection with Delta variant that it provides against previous lineages. We note there is an inherent trade-off between transmissibility and immune escape, and that the worst-case scenario of both very high transmissibility and immune escape is rejected a posteriori by the data. ; Tables S2-3 ). The results are valid for a population where the majority of immunity arose from prior infection (rather than vaccination), which is true for Delhi. Based on median estimates of the model (Fig 4A) and high transmissibility of the background Alpha lineage (18) , Delta should potentially be at least twice as transmissible as the wildtype lineage.

We note that after the massive Delta outbreak, new cases in Delhi and other North Indian states have stayed extremely low, with TPR in Delhi below 0.1% as of September 2021. This fits the serological picture presented in Fig 2, with a very high fraction of the population being recently infected and with good immunity to Delta. However, Delta outbreaks have continued in parts of India and elsewhere in the world, despite moderately high seropositivity or vaccination levels that were previously considered to be adequate (7, 19, 20) . We conclude that the Delta variant is capable of initiating fast-rising outbreaks in populations with immune responses to prior variants, resulting in reinfections and vaccination breakthroughs. Public health strategies may need to be revised to account for variants with heightened transmissibility and immune escape.

Refer to Web version on PubMed Central for supplementary material.

Suggestions by Prof Narendra K Arora, Prof G. Padmanabhan and member of the Scientific and Clinical Advisory Group of INSACOG are acknowledged. Coordinators of CSIR labs and centers are acknowledged for CSIR-Cohort. We gratefully acknowledge the authors from the originating laboratories and the submitting laboratories, who generated and shared via GISAID genetic sequence data on which this research is based. A full list acknowledging GISAID contributors for the data used in this study can be found in data file S3. 

All data used in the analysis is accessible through supplementary materials as data files or repository links (https://github.com/banijolly/ncov-Delhi-Epidemiology, DOI:10.5281/ zenodo.5521073) including accession numbers for all genomic data that has been deposited at GISAID. Code and data for replication of the Bayesian inference are available in the following repository: https://github.com/ImperialCollegeLondon/Delta_Variant_Delhi (DOI: 10.5281/zenodo.5524442)

This work is licensed under a Creative Commons Attribution 4.0 International (CC BY 4.0) license, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. To view a copy of this license, visit https:// creativecommons.org/licenses/by/4.0/. This license does not apply to figures/photos/artwork or other content included in the article that is credited to a third party; obtain authorization from the rights holder before using such material.

Delhi experienced an overwhelming surge of COVID-19 cases and fatalities peaking in May 2021 as the highly transmissible and immune evasive Delta variant replaced the Alpha variant. A) Seropositivity in CSIR cohort, sub-divided by nature of employment and use of public transport, is plotted for different time-periods (Phase I to Phase III, proportion +/-95% CI). Details are in table S1. B) Variability and temporal decline in neutralization capacity estimated by sVNT assay between Phase I and II (n=52). C) Serial antibody concentration measurements in initially seropositive subjects (n=91). Pattern suggestive of reinfections is shown (decline followed by rise, n=25). D) shows remaining data (n=66), with four indeterminate reinfection cases indicated by arrowheads. Values were inferred from a two category Bayesian transmission model fitted to mortality, serosurvey and genomic data from Delhi, India. A. Joint posterior distribution, with isoclines corresponding to the 90% and 50% enclosures of posterior density of the Delta variant immune escape and transmissibility increase relative to non-Delta categories. Immune escape has a median of 20% with (10%-50%) bCI50%, and transmissibility increase has a median of ×1.5 with (1.3-1.7) bCI50%. B. Delta fraction over time inferred by the model. Black dots represent genome sampling data points, with exact binomial confidence intervals. C. Serosurvey data (black dots) and inferred cumulative incidence for Delta and non-Delta variant categories. D. Mortality data (black dots) and inferred deaths assuming 50% under reporting. Other under-ascertainment scenarios are presented in the Supplementary Information.

SARS-CoV-2 seroprevalence among the general population and healthcare workers in India

Govt of National Capital Territory of Delhi (Twitter Feed of Health Minister

Declaration of major private hospitals as COVID care only and engagement of 4th and 5th year MBBS students in govt hospitals GNCTD dt 12

Insights from a Pan India Sero-Epidemiological survey (Phenome-India Cohort) for SARS-CoV2. eLife

Immunological memory to SARS-CoV-2 assessed for up to 8 months after infection

Estimating epidemiologic dynamics from cross-sectional viral load distributions

SARS-CoV-2 B.1.617.2 Delta variant replication and immune evasion

Safety and immunogenicity of the ChAdOx1 nCoV-19 vaccine against SARS-CoV-2: a preliminary report of a phase 1/2, single-blind, randomised controlled trial

Efficacy, safety, and lot to lot immunogenicity of an inactivated SARS-CoV-2 vaccine (BBV152): a, double-blind, randomised

Effectiveness of Covid-19 Vaccines against the B.1.617.2 (Delta) Variant

Estimating the effects of non-pharmaceutical interventions on COVID-19 in europe

Assessment of protection against reinfection with SARS-CoV-2 among 4 million PCR-tested individuals in Denmark in 2020: a population-level observational study

SARS-CoV-2 infection rates of antibody-positive compared with antibody-negative health-care workers in England: a large, multicentre, prospective cohort study (SIREN). The Lancet

Global initiative on sharing all influenza data -from vision to reality

A dynamic nomenclature proposal for SARS-CoV-2 lineages to assist genomic epidemiology

Assignment of Epidemiological Lineages in an Emerging Pandemic Using the Pangolin Tool. Virus Evolution

SARS-CoV-2 Seroprevalence in 12 Cities of India from

Estimated transmissibility and impact of SARS-CoV-2 lineage B.1.1.7 in England

Reduced sensitivity of SARS-CoV-2 variant Delta to antibody neutralization

Changing composition of SARS-CoV-2 lineages and rise of Delta variant in England

A new coronavirus associated with human respiratory disease in China

Nextstrain: real-time tracking of pathogen evolution

TreeTime: Maximum-likelihood phylodynamic analysis

IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era. Molecular biology and evolution

Using ggtree to Visualize Data on Tree-Like Structures. Current protocols in bioinformatics

Estimating the effects of non-pharmaceutical interventions on COVID-19 in Europe

Genomics and epidemiology of the P.1 SARS-CoV-2 lineage in Manaus

State-level tracking of COVID-19 in the United States

On the Integral Equation of Renewal Theory. The Annals of Mathematical Statistics

On Age-Dependent Binary Branching Processes

COVID-19 infection fatality ratio: Estimates from seroprevalence

Estimating COVID-19 infection fatality rate in Mumbai during 2020

A comparison of five epidemiological models for transmission of SARS CoV-2 in India

Reconstructing the COVID-19 epidemic in Delhi, India: Infection attack rate and reporting of deaths

Prevalence of SARS-CoV-2 infection in India: Findings from the national serosurvey

Estimating under-reporting of COVID-19 cases in Indian states: an approach using a delay-adjusted case fatality ratio

Underreporting COVID-19: the curious case of the Indian subcontinent

Estimating the effectiveness of routine asymptomatic PCR testing at different frequencies for the detection of SARS-CoV-2 infections

SARS-CoV-2 seroprevalence in 12 cities of India from

Quantifying antibody kinetics and RNA detection during early-phase SARS CoV-2 infection by time since symptom onset. eLife

Modeling of epidemics using hierarchical Bayesian models

Semi-Mechanistic Bayesian Modeling of COVID-19 with Renewal Processes

Genomic characterization and Epidemiology of an emerging SARS-CoV-2 variant in Delhi, India