Insights from a Pan India Sero-Epidemiological survey (Phenome-India Cohort) for SARS-CoV2
Abstract
To understand the spread of SARS-CoV2, in August and September 2020, the Council of Scientific and Industrial Research (India) conducted a serosurvey across its constituent laboratories and centers across India. Of 10,427 volunteers, 1058 (10.14%) tested positive for SARS-CoV2 anti-nucleocapsid (anti-NC) antibodies, 95% of which had surrogate neutralization activity. Three-fourth of these recalled no symptoms. Repeat serology tests at 3 (n = 607) and 6 (n = 175) months showed stable anti-NC antibodies but declining neutralization activity. Local seropositivity was higher in densely populated cities and was inversely correlated with a 30-day change in regional test positivity rates (TPRs). Regional seropositivity above 10% was associated with declining TPR. Personal factors associated with higher odds of seropositivity were high-exposure work (odds ratio, 95% confidence interval, p value: 2.23, 1.92–2.59, <0.0001), use of public transport (1.79, 1.43–2.24, <0.0001), not smoking (1.52, 1.16–1.99, 0.0257), non-vegetarian diet (1.67, 1.41–1.99, <0.0001), and B blood group (1.36, 1.15–1.61, 0.001).
Introduction
The World Health Organization declared SARS-CoV-2 infection as a pandemic on March 11, 2020 (WHO, 2020). Within 2 weeks, India announced a lockdown strategy that severely influenced the growth of the pandemic, which was initially very focal in the large cities, gathering pace and spreading to smaller cities and towns as the nation unlocked for societal and economic considerations.
Early literature pointed towards asymptomatic transmission of SARS-CoV-2 and raised the need for extended testing (Nishiura et al., 2020; Bai et al., 2020). While RT-PCR was an undisputed choice for establishing a positive infection, serosurveillance revealed that many more were probably getting infected without manifesting symptoms (Brown and Walensky, 2020; Oran and Topol, 2020). Initial estimates of asymptomatic infection rate from the West were around 40–45% (Oran and Topol, 2020).
In India, the first case of Covid was reported on January 30, 2020 (Andrews et al., 2020). Serological surveys have confirmed that spread beyond the Indian megacities was minimal in early May–June, with less than 1% seropositivity outside the designated containment zones, suggesting that the lockdown had been effective in limiting the spread (Murhekar et al., 2020a). This was not without human and economic costs. By the end of June, migrant workforces caught in the cities during the lockdown were sent to their rural homes, which may have led to subsequent rapid, multifocal rise in cases in July 2020. By October 2020, total new cases began to decline with further outbreaks limited to a few geographies. Thus, the period of August–September 2020 is an important transition point. The present study is one of two studies at a national scale that was designed to assess the spread of infection. A national serosurvey in 70 districts of India, conducted by the Indian Council of Medical Research (ICMR), had a reported crude positivity rate of about 10% (Murhekar et al., 2020b). Since only a few districts of India account for the majority of urban areas, the ICMR survey is not representative of Indian cities (Murhekar et al., 2020a), where a combination of higher population density and indoor lives has led to greater spread of disease than rural areas. Existing city-wise serosurveys are variable in target population as well as methodology and cannot be easily compared to determine the relative course of the pandemic (Ray et al., 2020; Satpati et al., 2020; Babu et al., 2020; Sharma et al., 2020). The current study was launched by the Council of Scientific and Industrial Research (CSIR) in its more than 40 constituent laboratories and centers spread all over the country, representing a wide range of ethnicities, geosocial habitats, and occupational exposures in the form of a longitudinal cohort (Phenome-India Cohort). Though limited by not being a population-denominated study, the cohort includes permanent staff, families of staff members, students, and temporary employees proving support services such as security, sanitation, housekeeping, etc. This is a diverse microcosm of India encompassing multiple socioeconomic groups and has the advantage of permitting deeper assessments such as questionnaire surveys and periodic reassessment of humoral antibody response in those found to be seropositive. Our data is thus the first that permits valid comparisons between many important urban regions of India. Here, we report results from phase 1 of this study covering the critical period of August–September 2020.
Results
Seropositivity varied widely across India
In 10,427 subjects from over 17 states and 2 union territories, the average seropositivity was 10.14% (95% confidence interval [CI] 9.6–10.7), but varied widely across locations (Figure 1). We found that 95% of the seropositive subjects also had significant neutralizing activity, suggesting at least partial immunity (Figure 1—source data 1).
Seropositivity, population density, and trajectory of new infections
As expected from the known outward spread of infection from large Indian cities, seropositivity was greater in regions with higher population density (Figure 2A). Changes in test positivity rate (TPR) are a relatively robust marker of the local level of transmission and are preferred when absolute number of tests or test rates are variable, as was the case here. Lab-wise seropositivity was correlated with the regional change in TPR. By this measure, regional transmission of SARS-CoV2 was inversely correlated to local seropositivity (Figure 2B). Seropositivity of 10% or more was associated with reductions in TPR, which may mean declining transmission (Figure 2—source data 1).
Survey-based correlates of seropositivity
Out of 861 seropositive subjects who also provided data on symptomatology, 647 subjects (75.3%) did not recall any of the nine symptoms asked for (two of these did not provide gender data) (Supplementary file 1). Among the minority of subjects with symptoms, the most reported symptom constellation was those of a mild flu-like disease with fever (~50%) as the most frequent symptom. Loss of smell or taste was uncommonly reported (~25% of symptomatic subjects) (Supplementary file 2).
We further examined associations of other available variables with seropositivity to explore potential factors that modulate risk of infection in India. Apart from gender and age, distribution of the other variables recorded in CSIR-cohort (prevalence of smoking, diet, physiological parameters like ABO blood group type) was similar to the national averages and the sample can be considered representative (Patidar and Dhiman, 2021; Mohan et al., 2018; Government of India, 2014). The univariate associations are shown in Figure 3, separately for each gender. Due to gender imbalance and possible confounding between various parameters, significance of associations was further tested in a balanced iterative logistic regression (Figure 3—figure supplement 1). The strongest gender-independent associations were with occupation and mode of transport. Outsourced staff performing support services such as security, housekeeping, etc., and subjects using public transport were more likely to be seropositive. In males, smoking and vegetarian diet was associated with lower seropositivity.
Blood group type was reported for 7496 subjects. Blood group distribution of subjects in our study was similar to national reference based on a recent systematic review (Patidar and Dhiman, 2021). Seropositivity was significantly different between different groups, being highest for blood group type AB (10.19%) followed by B (9.94%), O (7.09%), and A (6.97%). Blood group O was found to be associated with a lower seropositivity rate, with an odds ratio (OR) of 0.76 (95% CI 0.64–0.91, p=0.018) vs. non-O blood group types, while B appeared to be at high risk with an OR of 1.36 (95% CI 1.15–1.61, p=0∙001). While blood group A had an OR of 0.78, the association was not found to be significant (p=0.10) and a similar observation was made with blood group AB (p=0.35), it had an OR of 1.27 (Supplementary file 3). Rh factor was not found to have a significant association with seropositivity (p=0.35).
Stability of humoral response to SARS-CoV-2
Of 607 subjects whose samples were collected again at 3 months, anti-nucleocapsid (anti-NC) antibody levels were similar or higher for most, with 17 (2.8%) becoming seronegative (Figure 4A). In contrast, 34 subjects (5.6%) did not have neutralizing activity based on a surrogate measure (>20% inhibition of receptor-spike protein binding; Figure 4B). Of 175 subjects whose samples were collected again at 5–6 months, 8 (4.6%) became seronegative (Figure 4C). In contrast, neutralizing activity was not present in 31 (17.7%) subjects (Figure 4D, Figure 4—source data 1).
Discussion
This study, which recruited subjects from 24 cities in India, provides an important and timely snapshot of the spread of SARS-CoV2 pandemic across India shortly before the peak of new cases. It confirms that by September 2020 a large pool of recovered Indians with at least partial immunity existed. Between our study and the other national serosurvey at the same time, more than a hundred million Indians were likely to belong to this category. The subsequent nationwide decline of new cases starting in October 2020 can be well understood through these observations. There is some evidence of declining transmission in high-seropositivity regions within September, based on falling TPRs, but due to changes in tests and expansion of testing, this can only be indirectly inferred. As shown in our study, the fraction of such recovered subjects with resistance to reinfection was more than double among people performing high-contact jobs and using public transport. Thus, it is not surprising that in combination with a strong emphasis on masking and distancing, new cases started declining soon after this serosurvey. As of January 2021, despite the onset of winter and new year festivities, India has not seen major outbreaks. It is important to note that the crude seropositivity rate reported by us needs adjustments for demographics, fraction of infected subjects who may not develop antibodies, test characteristics, among others, to be a true population seropositivity rate. Here, we intentionally avoid such adjustments since it provides a false sense of precision with too many unknowns and unmeasurables. We focus on the more meaningful variation of the crude rate, its clinical correlates, and public health implications. We are confident from the data that a very large pool of recovered and immune subjects existed by September 2020, as stated, and expect that downward adjustments for national demographics will be counterbalanced by upward adjustments for almost 20% of RT-PCR-proven asymptomatic infections who develop transient antibody responses (Koopmans and Haagmans, 2020).
Apart from the seropositivity rate, our data also reveals important associations between demographic, physiological, lifestyle-related, and occupational attributes with susceptibility to infection. The workforce in our cohort comprised adult population, and no major difference was observed in seropositivity among different age groups or those with comorbidities. Males were found to be more susceptible, in agreement with other published reports (Ortolan et al., 2020). However, there were fewer females in our study and many of the occupational responsibilities with higher chances of exposure, like that of security personnel, were skewed towards males. On iteratively ran regression models, we did not find gender to be a predictor for seropositivity.
ABO blood group type has been shown to be associated with SARS-CoV-2 infection, but the results are variable in different studies. Most studies found that the O group is associated with a lower risk of infection or severity and blood group A was reported to be high risk in some studies (Zhao et al., 2020; Wu et al., 2020; Göker et al., 2020; Ellinghaus et al., 2020). In a meta-analysis authored by Golinelli et al., 2020, they observed positive association with A blood group, while blood group O was to be associated with lesser positivity using a random effects model. Another study from India observed blood group O to be associated with less mortality while blood group B with higher mortality when they analyzed the national data available (Padhi et al., 2020). A complex molecular interaction is said to play a significant role, and the molecular pathways for the same need to be elucidated for the effects observed especially with blood group O, which was also confirmed by our study.
The associations with diet and smoking are intriguing but preliminary. It has been proposed that a fiber-rich diet may play an important role in COVID-19 through anti-inflammatory properties by modification of gut microbiota (Conte and Toraldo, 2020). A recent review has highlighted the role of trace elements, nutraceuticals, and probiotics in COVID-19 (Jayawardena et al., 2020). The negative association with smoking has been reported elsewhere, but not shown to be causal (Fontanet et al., 2020; Makoto Miyara et al., 2020; Rossato et al., 2020; Petrilli et al., 2020; Farsalinos et al., 2020). Further exploration is necessary before reaching any conclusions, especially since seropositivity is an imperfect marker of infection risk and may equally well be explained by altered antibody response and dynamics.
The stability of antibody response is poorly understood, especially in India (Figueiredo-Campos et al., 2020; Seow et al., 2020; Deshpande et al., 2020). Our data shows that while anti-NC antibodies provide long-lasting evidence of viral exposure or infection, about 20% of seropositive individuals lack meaningful neutralization activity after 5–6 months. Using more stringent measures (more than 30% inhibition of surrogate receptor-spike protein binding), the loss of neutralization may be even higher. We speculate that this may be related to recurrence of outbreaks in March 2021, after the peak in September 2020.
The aggregate seropositivity of 10.14% in our multicentric study suggests that India had a large pool of recovered immune subjects by September 2020, especially among its high-contact workers and people using public transport, leading to a decline in new infections. However, the duration of such immunity may not be sufficient to prevent future outbreaks, even in highly affected regions.
Materials and methods
Study design, sampling, and data collection
Request a detailed protocolThe longitudinal cohort study was approved by the Institutional Human Ethics Committee of CSIR-IGIB vide approval CSIR-IGIB/IHEC/2019–20. The participation was voluntary, and participants had to fill an online informed consent with a consent to publish and share findings and deidentified data. Online statistical tool was utilized to calculate the minimum sample size for estimating seropositivity of about 5% with 10% precision (0.005) with 95% confidence to be 7300 (Dhand and Khatkar, 2014). In this study, we enrolled >10,000 subjects. 10,427 adult subjects working in the CSIR laboratories and their family members enrolled for the study based on voluntary participation. Informed consent was obtained from all the participants, and the samples were collected maintaining all recommended precautions. Blood samples (6 ml) were collected in EDTA vials from each participant, and antibodies to SARS-CoV-2 NC antigen were measured using an Electro-chemiluminescence Immunoassay (ECLIA)-Elecsys Anti-SARS-CoV-2 kit (Roche Diagnostics) as per the manufacturer’s protocol. This approved assay is considered a method of choice when a single test is to be deployed (Krüttgen et al., 2021). A Cut-off index COI >1 was considered seropositive. Positive samples were further tested for neutralizing antibody (NAB) response directed against the spike protein using GENScript cPass SARS-CoV-2 Neutralization Antibody Detection Kit (GenScript, USA), according to the manufacturer’s protocol. This is a blocking ELISA used for qualitative detection of total neutralizing antibodies against SARS-CoV-2 virus in plasma. A value of 20% or above was considered to have neutralizing ability.
All the participants were requested to fill an online questionnaire, which included information on date of birth, gender, blood group, type of occupation, history of diabetes, hypertension, cardiovascular disease, liver and kidney disease, diet preferences, mode of travel, contact history, and hospital visits. These forms were then downloaded in MS Excel data format and merged with registration forms filled at the time of sample collection based on unique IDs.
Data and statistical analysis
Request a detailed protocolRegion-wise and total seropositivity was calculated from the fraction of samples positive for antibodies to SARS-CoV-2 NC antigen. Data regarding RT-PCR/rapid antigen testing and positive cases was gathered from http://www.covid19India.org. Change in TPR, a robust parameter for estimating the level of infection transmission when the level of testing is variable, was calculated as per the following equation:
*DOC = date of collection.
IGIB, New Delhi, and NAL, Bengaluru, were removed from this change in TPR analysis for the sample collection was spread over 2–3 weeks in these labs. Questionnaire-based variables were assessed for response types and blank fields, that is, responses that were not provided by the participants of the survey. Based on multiple response types for each variable, categories were made to assign the response to either of the categories. For visualization, ggpubr (v0.4.0), ggrepel (v0.8.2), and ggplot2 (v3.3.2) packages were used in R. No data imputation was carried out. Chi-square test was performed to evaluate variables that had a significant association with outcome of being tested positive (p<0.05) along with OR with 95% CI. An adjusted p value was obtained through Bonferroni correction method for multiple comparison testing. Following the chi-square test, an iterative logistic regression was carried out on a balanced dataset. Variance inflation factor (VIF) was separately evaluated to assess multicollinearity. Statistical analysis and model development were carried out with visualization in R programming environment version 3.6.1, MS Excel 2016, and OriginPro V2021; faraway (v1.0.7) package was utilized for estimation of VIF.
Role of the funding source
Request a detailed protocolThe sponsor of this study had no role in the study design, data collection, data analysis, data interpretation, or writing of the report. The corresponding authors had full access to all the data in the study and had final responsibility for the decision to submit for publication.
Data availability
Source data files have been provided for Figures 1, 2,3 and 4. For Figure 2B, data is available in source data file and the district wise data listed in the source file on the number of confirmed cases and tests was taken from https://www.covid19india.org/. From this site, the data has to be manually extracted for desired time points by manually hovering over the number of cases and tests done graph when a district is selected after selection of the state. Code of Data Analysis and Model Definition Files- is available on Github at https://github.com/rintukutum/pan-india-csir-sero-epi (copy archived at https://archive.softwareheritage.org/swh:1:rev:7b789acbdb09c6842ab208c0ea84b52e88be43d7).
References
-
First confirmed case of COVID-19 infection in India: a case reportIndian Journal of Medical Research 151:490–492.https://doi.org/10.4103/ijmr.IJMR_2131_20
-
Targeting the gut–lung microbiota axis by means of a high-fibre diet and probiotics may have anti-inflammatory effects in COVID-19 infectionTherapeutic Advances in Respiratory Disease 14:175346662093717.https://doi.org/10.1177/1753466620937170
-
Neutralizing antibody responses to SARS-CoV-2 in COVID-19 patientsIndian Journal of Medical Research 152:82–87.https://doi.org/10.4103/ijmr.IJMR_2382_20
-
WebsiteStatulator: an online statistical calculatorSample Size Calculator for Estimating a Single Proportion. Accessed June 17, 2020.
-
Genomewide association study of severe Covid-19 with respiratory failureThe New England Journal of Medicine 383:1522–1534.https://doi.org/10.1056/NEJMoa2020283
-
Seroprevalence of anti-SARS-CoV-2 antibodies in COVID-19 patients and healthy volunteers up to 6 months post disease onsetEuropean Journal of Immunology 50:2025–2040.https://doi.org/10.1002/eji.202048970
-
The effects of blood group types on the risk of COVID-19 infection and its clinical outcomeTurkish Journal of Medical Sciences 50:679–683.https://doi.org/10.3906/sag-2005-395
-
Enhancing immunity in viral infections, with special emphasis on COVID-19: a reviewDiabetes & Metabolic Syndrome: Clinical Research & Reviews 14:367–382.https://doi.org/10.1016/j.dsx.2020.04.015
-
Determination of SARS-CoV-2 antibodies with assays from diasorin, Roche and IDvetJournal of Virological Methods 287:113978.https://doi.org/10.1016/j.jviromet.2020.113978
-
ReportLow Rate of Daily Active Tobacco Smoking in Patients with Symptomatic COVID-19Qeios.
-
Assessment of tobacco consumption and control in IndiaIndian Journal of Clinical Medicine 9:117991611875928.https://doi.org/10.1177/1179916118759289
-
Prevalence of SARS-CoV-2 infection in India: findings from the national serosurvey, May-June 2020Indian Journal of Medical Research 152:48–60.https://doi.org/10.4103/ijmr.IJMR_3290_20
-
Estimation of the asymptomatic ratio of novel coronavirus infections (COVID-19)International Journal of Infectious Diseases 94:154–155.https://doi.org/10.1016/j.ijid.2020.03.020
-
Prevalence of asymptomatic SARS-CoV-2 infection : a narrative reviewAnnals of Internal Medicine 173:362–367.https://doi.org/10.7326/M20-3012
-
Does gender influence clinical expression and disease outcomes in COVID-19? A systematic review and meta-analysisInternational Journal of Infectious Diseases 99:496–504.https://doi.org/10.1016/j.ijid.2020.07.076
-
ABO blood group system is associated with COVID-19 mortality: an epidemiological investigation in the indian populationTransfusion Clinique Et Biologique 27:253–258.https://doi.org/10.1016/j.tracli.2020.08.009
-
Distribution of ABO and rh (D) Blood groups in India: a systematic reviewISBT Science Series 16:37–48.https://doi.org/10.1111/voxs.12576
-
Current smoking is not associated with COVID-19European Respiratory Journal 55:2001290.https://doi.org/10.1183/13993003.01290-2020
Article and author information
Author details
Funding
CSIR (MLP-2007)
- Shantanu Sengupta
CSIR (MLP-2002)
- Satyartha Prakash
Department of Science and Technology, India (GAP-0192)
- Nitin Bhatheja
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
SSG would like to acknowledge CSIR grant MLP 2007 for this work. SN, AKB, and RajatU would like to thank CSIR for their fellowships. SPrakash would like to acknowledge CSIR grant MLP-2002 (CSIR-IGIB) for the fellowship. NB would like to address CSIR grant GAP-0192 (CSIR-IGIB) for this work. We would also like to thank the directors of all the CSIR Institutes for facilitating the study. We also thank Pushpesh Ranjan, Jitendra, Neeraj Kumar, Abhijeet, and Rajkumar from AMPRI. V Santosh Kumar from CCMB. Dr. Chandra Prakash Pandey from CDRI. Vipul Sharma, Akansh Agarwal, Hansraj Choudhary, Vijay Chatterjee, Narendra Meena, Ved Prakash, Alok Mishra, Navin Singhal, Ankit Shukla, and Sudeep Rathore from CEERI. Avilash S, Rani C, and Naveen Shashidhar Kumbar from CFTRI. Swachchha Majumdar from CGCRI. Dr. Dyaneshwar Umrao Bawankule, Dr. Debabrata Chanda, Pankaj Shukla, Sanjay Singh, Dr. Dayanandan Mani, Ravi Kumar, Pankaj Yadav, and Parmanand Kumar from CIMAP. D C Sharma, Dr. Neelam J Gupta, A K Jain, and Sudhansu Bhagat from CRRI. Pankaj Pandey, Rajesh, Dr. Mohammed Faruq, and Ajay Pratap Singh from IGIB. Yogita Singh and Karvan Kaushal from CSIO. Jaykumar Patel, Shrikant Khandare, and Dr. Kannan Srinivasan from CSMCRI. Dr. Prakash L Ganapathi, G Bhatt, Shashikala U, and Shashidhar KN from NAL. Dr. Vidyadhar Mudkavi, Ravichandran C, and Sunil Babu MG from 4PI. Dr. Robin Singh, Mahesh S, Mohit Kumar Swarnakar, and Dr. Pankaj Kulurkar from IHBT. Saikat Chakrabarti and Sandip Paul from IICB. Siva Ranjith and B Vijay Kumar from IICT. Sajad Ahmed from IIIM. Rene Christina, Neha Bansal, and Ayan Banerjee from IIP. SK Goyal from NEERI-Delhi. Dr. PR Meganathan and Dr. Shaik Basha from NEERI-Hyderabad. Antara Sharma from NEIST. Dr. Shuchismita Benzwal and Chaitanya Dinesh from NGRI. Shana S Nair from NIIST. N Anandavalli and P Vasudevan from SERC. Vibha Malhotra Swaney from TKDL. Rashmi Arya and Prafulla Malwadkar from URDIP. Ajeet Singh and Dr. RK Sinha from HRDC.
Ethics
Human subjects: The longitudinal cohort study was approved by Institutional Human Ethics Committee of CSIR-IGIB vide approval CSIR-IGIB/IHEC/2019-20. The participation was voluntary, and participants had to fill an online informed consent with a consent to publish and share findings and deidentified data.
Copyright
© 2021, Naushin et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 6,873
- views
-
- 747
- downloads
-
- 21
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Medicine
- Microbiology and Infectious Disease
- Epidemiology and Global Health
- Immunology and Inflammation
eLife has published articles on a wide range of infectious diseases, including COVID-19, influenza, tuberculosis, HIV/AIDS, malaria and typhoid fever.
-
- Cell Biology
- Medicine
The IncRNA Malat1 was initially believed to be dispensable for physiology due to the lack of observable phenotypes in Malat1 knockout (KO) mice. However, our study challenges this conclusion. We found that both Malat1 KO and conditional KO mice in the osteoblast lineage exhibit significant osteoporosis. Mechanistically, Malat1 acts as an intrinsic regulator in osteoblasts to promote osteogenesis. Interestingly, Malat1 does not directly affect osteoclastogenesis but inhibits osteoclastogenesis in a non-autonomous manner in vivo via integrating crosstalk between multiple cell types, including osteoblasts, osteoclasts, and chondrocytes. Our findings substantiate the existence of a novel remodeling network in which Malat1 serves as a central regulator by binding to β-catenin and functioning through the β-catenin-OPG/Jagged1 pathway in osteoblasts and chondrocytes. In pathological conditions, Malat1 significantly promotes bone regeneration in fracture healing. Bone homeostasis and regeneration are crucial to well-being. Our discoveries establish a previous unrecognized paradigm model of Malat1 function in the skeletal system, providing novel mechanistic insights into how a lncRNA integrates cellular crosstalk and molecular networks to fine tune tissue homeostasis, remodeling and repair.