SARS-CoV-2 seroprevalence in Chattogram, Bangladesh before the Delta surge, March-June 2021 =========================================================================================== * Taufiqur Rahman Bhuiyan * Juan Dent Hulse * Sonia T. Hegde * Marjahan Akhtar * Md. Taufiqul Islam * Zahid Hasan Khan * Ishtiakul Islam Khan * Shakeel Ahmed * Md Mamunur Rashid * Rumana Rashid * Emily S. Gurley * Tahmina Shirin * Ashraful Islam Khan * Andrew S. Azman * Firdausi Qadri ## Abstract In a representative serosurvey conducted March–June 2021, 64.1% (95%CrI 60.0– 68.1%) of Sitakunda subdistrict (Bangladesh) had anti-SARS-CoV-2 IgG antibodies after adjusting for age, sex, household clustering and test performance. Before the surge of Delta, most of the population had been infected despite low incidence of virologically-confirmed COVID-19. Keywords * SARS-CoV-2 * Seroepidemiologic studies * Bangladesh * COVID-19 ## Text Since the start of the COVID-19 pandemic, Bangladesh has reported more than 1.55 million cases and 27,270 deaths (*1*). Reported mortality and incidence rates appear to be substantially lower than many other countries. Without population-based seroprevalence estimates, it is difficult to know whether these differences are due to limited surveillance and healthcare seeking or true differences in incidence resulting from interventions or differing biological responses to infection. In early March 2021, cases across Bangladesh began to rise at the same time the Delta variant was detected in neighboring India. Publicly available sequencing data (*2*) indicates the SARS-CoV-2 Delta variant was first detected in the Chattogram region of Bangladesh in mid-May 2021 and since July 1, 2021, 99% (n=98/99) of the submitted viral genomes have been the Delta variant, similar to national trends. ## The Study We conducted a representative serosurvey to understand the prevalence of anti-SARS-CoV-2 IgG antibodies in the Sitakunda subdistrict (Chattogram district) of Bangladesh, which has a significant urban-rural gradient starting from Bangladesh’s second largest city, Chattogram. The survey was conducted in two periods, from 27-March to 13-April and from 23-May to 13-June (after Eid al-Fitr), due to a national COVID-19-related lockdown. We used two-stage sampling based on digitized satellite imagery, where we first selected a 1km2 grid-cell proportional to the estimated number of households in each grid, then randomly selected a structure weighted by whether it was multi- or single-story. We attempted to enroll all individuals ≤1 years old in each household. We tested participant serum for total antibodies (IgG, IgM and IgA) against the S1-RBD of SARS-CoV-2 using the Wantai SARS-CoV-2 Ab ELISA kit (Beijing Wantai Biological Pharmacy Enterprise Co., China) following manufacturer instructions. We corrected seroprevalence estimates for imperfect test performance, household clustering and individual-level covariates (e.g., age) using a previously documented Bayesian modeling approach and post-stratified results to match the target population structure (*3*). This study was approved by the icddr,b Research and Ethics Review Committees and the Johns Hopkins Bloomberg School of Public Health Institutional Review Board. Additional methodological details, including details on sampling and enrollment, are in the supplement, and code and data to reproduce analyses are at [https://github.com/HopkinsIDD/sitakunda-sarscov2-round1](https://github.com/HopkinsIDD/sitakunda-sarscov2-round1). Given the limited data on the immunoassay’s performance in South Asia and more generally on its performance months after infection, we conducted a validation study to estimate test performance. We estimated the sensitivity and the specificity of the assay by testing samples from 214 healthy participants from a 2014 cholera vaccine study (pre-pandemic) and 81 samples from 52 symptomatic PCR-confirmed SARS-CoV-2 infections. None of the positive controls were hospitalized for COVID-19 nor vaccinated, and samples were collected 3 to 275 days post-symptom onset. We estimated a specificity of 99.1% (95%CI 96.7-99.9%, n=212/214) and sensitivity of 87.7% (95%CI 78.5-93.9, n=71/81) for detecting previous infection with little evidence of decreasing sensitivity with time since infection (**Supplementary Table 4**). We enrolled 580 households and 2,307 individuals who provided a blood sample. A majority of participants were female (54%, n=1235/2307), with a median age of 28 (IQR: 16–45) years. Most participants reported working at home (37%), going to school (29%), or conducting business outside of their home (20%) as their main occupation in the month prior to enrollment. Among all study participants, 22 individuals (0.95%) reported ever having been tested for COVID-19 including 3 with positive results (3/3 were seropositive). One hundred and twenty-five (5.4%) reported being vaccinated with at least one dose of the Covishield (ChAdOx1, n=117) or Pfizer (n=1) SARS-CoV-2 vaccines. As of 19-June-2021, 6 days after the end of the survey, 6.0% of the entire Chattogram district population is reported to have received at least one dose of any vaccine including 4.6% with two doses (*4*). There were 1,443 (63%) seropositive individuals. Nearly all (98%) who reported having been partially (n=47/49) or completely vaccinated (n=75/76) were seropositive. At least one individual was seropositive in 85% of enrolled households and, on average, 62% of participants in each household were seropositive. We estimate that 31% of the total variability in seropositivity in the community is attributed to variation in seropositivity between households (intraclass correlation coefficient: 0.31; 95%CI: 0.27–0.36). We found evidence of a gradient in seropositivity by population density, with those living in areas of higher population density being significantly more likely to be seropositive; 69% of individuals living in the most population dense areas were seropositive compared with 52% of individuals living in the least population dense areas (p<0.0001, Supplemental Table 1). We found similar results using alternative metrics related to urbanicity (Supplemental Table 1). Adjusting for age, sex, household clustering, and test performance we estimated the seroprevalence of SARS-CoV-2 in Sitakunda to have been 64.1% (95% CrI 60.0–68.1%) among all individuals and 63.4% (95% CrI 59.2-67.6%) when considering only unvaccinated individuals (Table 1, Supplemental Table 3). We estimated a 7% higher risk of being seropositive in men compared to women (95% CrI: 1-13%). Risk generally increased with age, with those less than 10 years old having ≥33% lower risk than other age groups (Table 1, Supplemental Table 4). We found a similar adjusted seroprevalence in the population recruited before the lockdown (63.1%, 95% CrI 56.2-69.8%; n=665) and those collected after the lockdown (65.3%, 95% CrI 60.6-69.9%; n=1,643). View this table: [Table 1.](http://medrxiv.org/content/early/2021/09/30/2021.07.16.21260611/T1) Table 1. Estimated seroprevalence of SARS-CoV-2 in Sitakunda Upazila adjusted for sex, age, household clustering, and test performance among all individuals, both vaccinated and unvaccinated. In the catchment area of this serosurvey, only one healthcare facility provided SARS-CoV-2 testing (Bangladesh Institute of Tropical and Infectious Diseases). Among the 2,400 individuals that were tested from April 2020 until May 31, 2021, 705 (29%) tested positive for SARS-CoV-2 by RT-PCR. Crudely extrapolating our serologic estimates, by multiplying the estimated population size by adjusted seroprevalence among unvaccinated, suggests that more than 200,000 infections occurred during the same period, a much higher detected case to infection ratio than has been documented in most settings across the world (*5-6*). However, though the infection to case ratio is high, a substantial portion of missed infections were likely symptomatic; among the 161 respondents who reported having any COVID-19 like symptoms, only 57% (n=92) said they actively sought healthcare for treatment and only 1% (n=15/1442) of those who were seropositive were ever tested. ## Conclusions These results illustrate that prior to the June 2021 surge in cases in Bangladesh fueled by the Delta-variant, the majority of the population in Sitakunda had already been infected despite a relatively low incidence of reported virologically confirmed SARS-CoV-2 infections. These data also provide reassuring evidence on the immunogenicity of the Covishield/ChAdOx1 vaccine in this population, where almost all self-reported vaccinees were seropositive. Key limitations to these results include the relatively small geographic area covered by the survey and that we only assessed circulating antibodies to a single SARS-CoV-2 epitope, which does not fully capture the immune profile of individuals. In Bangladesh, where cases captured by surveillance are limited by healthcare seeking, even in population dense settings, representative seroprevalence surveys can help us continue to track the evolution of this pandemic. In addition to providing important validation data on a widely used immunoassay, our results can help lay the foundation for understanding the role of variant strains on key epidemiologic parameters, including our understanding of re-infection, and help set expectations for SARS-CoV-2 control in the months to come in the study area and beyond. ## Supporting information Supplement [[supplements/260611_file09.docx]](pending:yes) ## Data Availability Data supporting these estimates are publicly available in the sitakunda-sarscov2-round1 repository in Github. [https://github.com/HopkinsIDD/sitakunda-sarscov2-round1](https://github.com/HopkinsIDD/sitakunda-sarscov2-round1) ## Conflict of Interest We declare no competing interests. ## Author Bio Dr. Bhuiyan is an Associate Scientist at the International Centre for Diarrhoeal Disease Research, Bangladesh. His research interests include immunology, vaccinology, and enteric infections. View this table: [Supplemental Table 1.](http://medrxiv.org/content/early/2021/09/30/2021.07.16.21260611/T2) Supplemental Table 1. Descriptive statistics for serosurvey participants (n=2,307) in Sitakunda Upazila by seropositivity. This table includes sociodemographic factors, measures of urbanicity, COVID-like symptoms, COVID testing and vaccination, and COVID-related behaviors. View this table: [Supplemental Table 2.](http://medrxiv.org/content/early/2021/09/30/2021.07.16.21260611/T3) Supplemental Table 2. Descriptive statistics for unvaccinated, serosurvey participants (n=2,181) by seropositivity in Sitakunda Upazila. This table includes sociodemographic factors, measures of urbanicity, COVID-like symptoms, COVID testing and vaccination, and COVID-related behaviors. View this table: [Supplemental Table 3.](http://medrxiv.org/content/early/2021/09/30/2021.07.16.21260611/T4) Supplemental Table 3. Estimated seroprevalence of SARS-CoV-2 in Sitakunda Upazila adjusted for sex, age, household clustering, and test performance among unvaccinated individuals. View this table: [Supplementary Table 4.](http://medrxiv.org/content/early/2021/09/30/2021.07.16.21260611/T5) Supplementary Table 4. The number of positive controls used to estimate the empirical sensitivity of the Wantai total Ab assay and SARS-CoV-2 positivity by time since symptom onset. **Supplemental Figure 1.** Map of the study population in the Sitakunda Upazila (in green) in the Chattogram District of Bangladesh. The 580 enrolled households sampled in the serosurvey by enrollment time (pre versus post lockdown) and the two healthcare facilities in Sitakunda (BITID and the Sitakunda UHC) are indicated on the right side of the figure. **Supplemental Figure 2.** Map of the sampled clusters and sampled dwellings within clusters in the Sitakunda Upazila by samples drawn at three different periods: first draw (March 2021), second draw (May 2021), third draw (June 2021). Three separate sample draws were conducted due to interruption from the nationally imposed lockdown and a large percentage of non-residential structures among housing structures sampled from the satellite imagery. (A) Clusters were sampled 41 times with 14 structures each during the first two draws and 12 structures during the last draw. (B) We oversampled the number of structures by 40% to account for nonresidential buildings for a total of 574 sampled structures for the first two draws and 492 for the last draw. (C) Households were enrolled across the entire subdistrict of Sitakunda during each enrollment period and by sample draw (households enrolled pre-lockdown were only drawn from the first sample). **Supplemental Figure 3.** Average seropositivity among all participants in each sampled cluster (1 km2) by (A) population density (B) the minutes required to travel one minute (shown with all cluster data points and excluding the cluster with an outlier value) and (C) average household distance to Chittagong Port, the center of the Chattogram City. Loess smoothing is shown in blue with 95% confidence intervals shown in grey. ![Figure1](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/09/30/2021.07.16.21260611/F1.medium.gif) [Figure1](http://medrxiv.org/content/early/2021/09/30/2021.07.16.21260611/F1) ![Figure2](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/09/30/2021.07.16.21260611/F2.medium.gif) [Figure2](http://medrxiv.org/content/early/2021/09/30/2021.07.16.21260611/F2) ![Figure3](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/09/30/2021.07.16.21260611/F3.medium.gif) [Figure3](http://medrxiv.org/content/early/2021/09/30/2021.07.16.21260611/F3) ![Figure4](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/09/30/2021.07.16.21260611/F4.medium.gif) [Figure4](http://medrxiv.org/content/early/2021/09/30/2021.07.16.21260611/F4) ![Figure5](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/09/30/2021.07.16.21260611/F5.medium.gif) [Figure5](http://medrxiv.org/content/early/2021/09/30/2021.07.16.21260611/F5) ![Figure6](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/09/30/2021.07.16.21260611/F6.medium.gif) [Figure6](http://medrxiv.org/content/early/2021/09/30/2021.07.16.21260611/F6) ![Figure7](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/09/30/2021.07.16.21260611/F7.medium.gif) [Figure7](http://medrxiv.org/content/early/2021/09/30/2021.07.16.21260611/F7) ## Acknowledgements This work was supported by Bill and Melinda Gates Foundation [grant number INV-021879]. The authors would like to thank the Ministry of Health and Family Welfare (MOHFW) of Bangladesh. The authors would also like to express their sincere thanks to the staff members of icddr,b for their dedicated work in the field and laboratory during this pandemic situation. icddr,b is thankful and supported by the Governments of Bangladesh, Canada, Sweden, and the UK. ## Footnotes * Updated to include more data from after lockdown. * Received July 16, 2021. * Revision received September 30, 2021. * Accepted September 30, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.Bangladesh: WHO Coronavirus Disease (COVID-19) Dashboard With Vaccination Data [Internet]. [cited 2021 Sep 17]. Available from: [https://covid19.who.int](https://covid19.who.int) 2. 2.GISAID Initiative [Internet]. 2021 [cited 2021 Sep 22]. Available from: [http://www.gisaid.org](http://www.gisaid.org) 3. 3.Stringhini S, Wisniak A, Piumatti G, Azman AS, Lauer SA, Baysson H, et al. Seroprevalence of anti-SARS-CoV-2 IgG antibodies in Geneva, Switzerland (SEROCoV-POP): a population-based study. The Lancet. 2020 Aug 1;396(10247):313–319. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s0140-6736(20)31304-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F30%2F2021.07.16.21260611.atom) 4. 4.Bangladesh Directorate General of Health Services: Ministry of Health and Family Welfare. COVID-19 [Internet]. [cited 2021 Jul 28]. Available from: [http://103.247.238.92/webportal/pages/covid19-vaccination-update.php](http://103.247.238.92/webportal/pages/covid19-vaccination-update.php) 5. 5.Sitakunda Upazila - Banglapedia [Internet]. [cited 2021 Jun 26]. Available from: [https://en.banglapedia.org/index.php/Sitakunda\_Upazila](https://en.banglapedia.org/index.php/Sitakunda_Upazila) 6. 6.Chen X, Chen Z, Azman AS, Deng X, Sun R, Zhao Z, et al. Serological evidence of human infection with SARS-CoV-2: a systematic review and meta-analysis. The Lancet Global Health. 2021 May 1;9(5): e598–e609.