Seroprevalence of COVID-19 virus infection in Guilan province, Iran =================================================================== * Maryam Shakiba * Seyed Saeed Hashemi Nazari * Fardin Mehrabian * Seyed Mahmoud Rezvani * Zahra Ghasempour * Abtin Heidarzadeh ## Abstract **Background** The extent of infection by coronavirus disease 2019 has not been well known. In this study we aimed to determine seropositivity of COVID-19 virus infection in population of a highly affected area in north of Iran. **Methods** In a population-based cluster random sampling design through phone call invitation, a total of 196 household including 551 subjects agreed to participate in this study. Each participant were taken 50ml blood sample at health care center. Rapid test kits were used to detect antibody against COVID-19. Crude, population-weight adjusted and test performance adjusted prevalence of antibody seropositivity to SARS-CoV-2 were reported. **Results** The prevalence of antibody seropositivity was 22% (95%CI: 19-26%). The population weight adjusted estimate was 21% (95%CI: 14-29%) and test performance adjusted prevalence was 33% (95%CI: 28-39%). Based on these estimates the range of infected people in this province would be between 518000 and 777000. **Conclusion** The population seropositivity prevalence of COVID-19 virus infection indicated that the asymptomatic infection is much higher than the number of confirmed cases of COVID-19. This estimate can be used to better detect infection fatality rate and decide for public policy guidelines. Keywords * Asymptomatic Infections * COVID-19 * Seroprevalence ## Introduction Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) that was first reported in China has now spread to all countries throughout the world. The disease that subsequently called as coronavirus disease 2019 (COVID-19) has an extensive spectrum of manifestation including asymptomatic infection, mild disease of upper respiratory system, severe viral pneumonia with acute respiratory syndrome, and death (1). Asymptomatic or subclinical infections are one of the major challenges of COVID-19 with public health concerns as it is presumed that they can spread the infection and remained undiscovered in the community (2-4). According to the report on “Diamond Princess”, 51% of all confirmed cases including 10 crew and 308 passengers were asymptomatic(5). The prevalence of asymptomatic infection in the community is unclear yet, however, this is important to be identified for at least two reasons. First, the unknown infection fatality rate that is an important public health measures of COVID-19, would properly be estimated by adequate knowledge of proportion of all persons infected with novel coronavirus (6, 7). Current fatality rate estimates are based on confirmed cases multiplied by a factor representing asymptomatic cases (8, 9). Moreover, for achieving the so-called temporary herd immunity, the minimum level of population immunity to halt the spread of infection in the community would be 1-(1/R) (10). Where, R is the number of secondary cases generated by an infected individual in an ongoing epidemic and estimated to be 2 for Iran (11). Therefore, the information about the percentage of previously infected and hence immune people can help us to identify herd immunity, project epidemic and decide about public policy guidelines. According to World health organization (WHO), surveillance of antibody seropositivity including immunoglobulin G (IgG) and immunoglobulin M (IgM) against COVID-19 in a population can provide inferences about the extent and cumulative incidence of infection in the population (12). Guilan province was the second large city located at north of Iran that had several earliest number of COVID-19 soon after epidemic initiation in china. The epidemic curve in this province now has been stepped down making it an appropriate population-level location to test for the presence of active and past infections through seroprevalence survey. Therefore in this study we aimed to provide population-based seropositivity estimate of COVID-19 virus infection based on WHO protocol. ## Methods ### Study type and study population This was a cross-sectional population-based study to investigate seroprevalence of COVID-19 virus infection in Guilan province, north of Iran in April 2020. The study had been approved by IRB of Guilan University of medical sciences, Rasht, Iran. Individuals permanently resided in Guilan province irrespective of age were eligible for recruitment into the study. All persons living in the household including children invited to participate in the study. Household were defined as at least two people living in the same place. Residence of Institutional living centers such as nursing homes, prisons, and boarding schools and those subjects who refused to give informed consent, were under active treatment of COVID-19, or had any contraindication to venipuncture were excluded from the study. With a prior prevalence of 50% for coronavirus infection with a 5% precision and design effect of 1.24, a total of 525 subjects were considered for sample size. Guilan province with a population of 2354848 is located in north of Iran and consisted of 16 counties. The study conducted in three high incidence counties (Rasht, Anzali, Lahijan) and two low incidence counties (Astara, Roudbar) according to the protocol of World health organization. A multistage cluster random sampling approach was used to select the participants. The number of clusters in each county was assigned proportional to the population size of that county. Therefore, 15 clusters were considered for Rasht which is the capital city of province, and 5 cluster for each of the remaining 4 counties. Clusters were identified at random. Then in each cluster, the household were selected using random digit dialing of mobile telephone number of head of the household that was previously registered in a community-based surveillance of coronavirus infection by health care centers. ### Data collection A brief explanation of the research objectives followed by formal invitation of all family members to the study were given to the head of the household upon the phone call. For acknowledgment and to increase rate of participation a package of incentives including ethanol alcohol, face masks and kids incentives were considered for households. All invitation schedules were set to space enough to follow social distancing among households. At the day of attendance, after providing concise information about the risk and benefits accompanied by participation in this study, an informed consent was taken from the head of the household. An interviewer completed an electronic questionnaire containing participants’ demographic and exposure history information. Then, sample collectors in personal protective equipment drew 50 µL of capillary blood into an EDTA-coated microtainer. Tubes were labeled with participant ID. VivaDiag COVID-19 IgM/IgG from VivaChek was used for COVID-19-specific serological assay. According to manufacturer’s instruction, 10□µL of serum or whole blood sample was added into the sample port followed by adding 2 to 3 drops (70-100□µL) of dilution buffer. Test kits were read after about 15 minutes. According to the manufacturer, test performance characteristics showed a sensitivity of 97.1% in 11-14 days from infection and 81.25% in 4-10 days from infection. Test performance also verified on a sample of 30 PCR negative volunteer and 30 hospitalized patients. The sensitivity was 63.3%for both IgM and IgG (13). The specifity was reported to be 100% in both documents. ### Statistical analysis Demographic characteristics described as frequency and percent. The crude age-specific population estimate of COVID-19 positive test were estimate as proportion of total sample size. Then, in order to account for 3-stage cluster sampling, survey design with probability weights equal to the inverse probability of selection of counties, clusters in each county and individuals in each cluster and also finite population correction for each stage was used. The cumulative incidence estimates were further adjusted for Rapid test sensitivity and specifity. We used Blaker’s methods to estimate standard error for cumulative incidence which account for the uncertainty in the sensitivity and specificity. Blaker’s method provide smaller exact two-sided confidence interval estimates (14). For test performance adjustment we applied the results of previously published sensitivity and specifity as 63% and 100%, repectively. The raw and the survey-weight adjusted estimates were performed in Stata version 13 (StataCorp LP, College Station, TX). Adjustment for test performance characteristics was performed using epiR package in R software. ## Results From 632 household contact, 196 household including 551 subjects (49% male) were participated in this study. The major reasons for non-participation were concerning about getting infection in health facility center (17%), busy schedule (20%), lack of assurance to system (2%), no response (31%), and other non-specified reasons (28%). Table 1 provided demographic characteristic of unadjusted sample and population-weight adjusted of the sample and Guilan province estimates. The sample distribution was not meaningfully deviated from that of the province except for living place (32% in sample lived in village, 43% in the province). View this table: [Table 1.](http://medrxiv.org/content/early/2020/05/01/2020.04.26.20079244/T1) Table 1. Demographic characteristics of sample relative to Guilan province estimates From total, 11 subjects could not be tested and 12 had invalid test results (i.e. test kit showed no control band). From remaining 528, 117 subjects were positive for either IgM or IgG, crude prevalence rate of 22% (95% confidence interval (CI): 19-26%). The population-weight adjusted prevalence estimates was 21% (95%CI: 14-29%). Table 2 presents the unadjusted, population weight adjusted prevalence estimates, and improved estimation by adjusting for test performance results. The adjusted test performance results yielded prevalence rate as 33% (95%CI: 28-39%). View this table: [Table 2.](http://medrxiv.org/content/early/2020/05/01/2020.04.26.20079244/T2) Table 2. Seropositivity prevalence estimates in Guilan province View this table: [Table 3.](http://medrxiv.org/content/early/2020/05/01/2020.04.26.20079244/T3) Table 3. Seropositivity prevalence estimates according to the study variables Age groups-specific prevalence rate showed that the highest prevalence of seropositivity was among subjects older than 60 years and children less than 5 years had the lowest seropositivity. There was no significant difference in seropositivity prevalence estimates in terms of sex, age group, place of residence, job and county of residential place. From total 528 participants who provided blood sample, 366 subjects (69%) had none of the symptoms attributed to COVID-19 including general, respiratory, or gastrointestinal symptoms. Among them, 18% (65 subjects) had positive antibody test results which indicate the percent of quite asymptomatic cases in our sample. ## Discussion In this study, the seroprevalence estimate of SARS-COV-2 antibodies after adjusting for population and test performance characteristics was between 22 and 33%. The current results is much higher than previously seroprevalence estimate in California, USA which was between 2.49% and 4.16% (15). Our estimate is also higher than other preliminary findings in Italy (16) with 10%, Germany (17) with 14%, Colorado (18) point to 2% seropositivity, and New York city around 14% (19). In this study test performance characteristic had the highest influence on reported estimates as test kit had lower sensitivity resulting in about 40% of false negative cases that might be missed by the test. This study was among the first population-based seroprevalence study according to WHO protocol which incorporated cluster random sampling approach among community households of a highly affected area. In comparison to the estimates of western countries and even preliminary results of other provinces in our country, the high seroprevalence in our province may be explained by large amount of economic relationship with China in a free trade zone that is located at north of the Guilan province. The zone has been placed in Anzali that has the highest seropositivity prevalence in Guilan province. Along with other seroprevalence surveys reported globally, the results of our study showed that the number of infection is much higher than the number of confirmed case reported based on PCR results. The community seropositivity prevalence may also represent the percentage of asymptomatic infection by novel coronavirus. In our study, the percentage of seropositivity in non PCR-confirmed or symptomatic COVID cases was 18% which may represent the rate of asymptomatic infection. The asymptomatic proportion of novel coronavirus infection in our study is also higher than previously reported SARS and MERS coronavirus epidemic which has reported between 0-7.5% (20-24). Our seroprevalence estimate can be used for approximating infection fatality rate rather than case fatality rate. Based on the estimates of current study, the range of infected people in this province would be between 518000 and 777000. As of April 23 (the completion date of survey) 617 people have died of confirmed COVID-19 in Guilan province. This number corresponds to an infection fatality rate between 0.08-0.12% that is much lower than currently reported estimates of case fatality rate for COVID-19. The range of infection fatality rate in California was 0.12-0.2% (15). Previous estimate of case fatality rate using lag time for fatality in china was between 0.25-3.0% (25) far higher than currently estimated infection fatality rate. In our study, the household affirmative response rate was 31% and it may induce some sampling bias favoring those with prior COVID-like symptoms seeking antibody test. However, in our study only 11 subjects (2%) had a history of COVID confirmed diagnosis and 125 subjects (23%) had a combination of fever, dry cough, ache and tiredness as the major symptoms of COVID-19. Otherwise, bias favoring individuals in good health who are capable of participating in the test center may result in an underestimation of true prevalence. ## Conclusion This study revealed that the prevalence of antibody seropositivity in a highly affected area was between 21 and 33%. The resulting number of infection in the population is much higher than the number of confirmed cases currently reported in Guilan province. These estimates have implications for future public policy guidelines and showed that the true infection fatality rate would be much lower than current estimates. ## Data Availability Data are not publicly available. ## Acknowledgment This research has been sponsored by Iran’s Ministry of health and performed under a research grant from deputy of research at Guilan University of medical sciences. The authors are grateful to all health centers’ employees and personnel for cooperating and conducting the survey as well as all participants that took part in this research. * Received April 26, 2020. * Revision received April 26, 2020. * Accepted May 1, 2020. * © 2020, 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.Zhu N, Zhang D, Wang W, Li X, Yang B, Song J, et al. A Novel Coronavirus from Patients with Pneumonia in China, 2019. The New England journal of medicine. 2020;382(8):727–33. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2001017&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F01%2F2020.04.26.20079244.atom) 2. 2.Bai Y, Yao L, Wei T, Tian F, Jin DY, Chen L, et al. Presumed Asymptomatic Carrier Transmission of COVID-19. Jama. 2020. 3. 3.Wang M, Wu Q, Xu W, Qiao B, Wang J, Zheng H, et al. Clinical diagnosis of 8274 samples with 2019-novel coronavirus in Wuhan. medRxiv. 2020:2020.02.12.20022327. 4. 4.Gao WJ, Li LM. [Advances on presymptomatic or asymptomatic carrier transmission of COVID-19]. Zhonghua liu xing bing xue za zhi = Zhonghua liuxingbingxue zazhi. 2020;41():485–8. 5. 5.National Institute of Infectious Diseases (NIID). Field Briefing: Diamond Princess COVID-19 Cases. Tokyo: NIID2020 [cited 2020 April 15]. Available from: [https://www.niid.go.jp/niid/en/2019-ncov-e/9417-covid-dp-fe-02.html](https://www.niid.go.jp/niid/en/2019-ncov-e/9417-covid-dp-fe-02.html). 6. 6.Flahault A. Has China faced only a herald wave of SARS-CoV-2? Lancet (London, England). 2020;395(10228):947. 7. 7.Russell TW, Hellewell J, Jarvis CI, van Zandvoort K, Abbott S, Ratnayake R, et al. Estimating the infection and case fatality ratio for coronavirus disease (COVID-19) using age-adjusted data from the outbreak on the Diamond Princess cruise ship, February 2020. Euro surveillance : bulletin Europeen sur les maladies transmissibles= European communicable disease bulletin. 2020;25(12). 8. 8.Onder G, Rezza G, Brusaferro S. Case-Fatality Rate and Characteristics of Patients Dying in Relation to COVID-19 in Italy. Jama. 2020. 9. 9.Spychalski P, Błażyńska-Spychalska A, Kobiela J. Estimating case fatality rates of COVID-19. The Lancet Infectious Diseases. 2020. 10. 10.Anderson R, May R. Infectious Diseases of Humans: Dynamics and Control Oxford, UK: Oxford Science Publications; 1992. 11. 11.Lytras T, Panagiotakopoulos G, Tsiodras S. Estimating the ascertainment rate of SARS-CoV-2 infection in Wuhan, China: implications for management of the global outbreak. medRxiv. 2020:2020.03.24.20042218. 12. 12.World Heath Organization. Population-based age-stratified seroepidemiological investigation protocol for COVID-19 virus infection 2020 [cited 2020 April 01]. Available from: [https://www.who.int/publications-detail/population-based-age-stratified-seroepidemiological-investigation-protocol-for-covid-19-virus-infection](https://www.who.int/publications-detail/population-based-age-stratified-seroepidemiological-investigation-protocol-for-covid-19-virus-infection). 13. 13.Cassaniti I, Novazzi F, Giardina F, Salinaro F, Sachs M, Perlini S, et al. Performance of VivaDiag COVID-19 IgM/IgG Rapid Test is inadequate for diagnosis of COVID-19 in acute patients referring to emergency room department. Journal of medical virology.n/a(n/a). 14. 14.Blaker H. Confidence curves and improved exact confidence intervals for discrete distributions. Canadian Journal of Statistics. 2000;28:783–98. 15. 15.Bendavid E, Mulaney B, Sood N, Shah S, Ling E, Bromley-Dulfano R, et al. COVID-19 Antibody Seroprevalence in Santa Clara County, California. medRxiv. 2020:2020.04.14.20062463. 16. 16.What the first data of Robbio reveal, the only Italian country to take immunity tests on all citizens. 2020 [cited 2020 April 23]. Available from: [https://it.businessinsider.com/esclusiva-cosa-rivelano-i-primi-test-di-robbio-primo-paese-italiano-afare-i-test-sullimmunita-a-tutti-i-cittadini/](https://it.businessinsider.com/esclusiva-cosa-rivelano-i-primi-test-di-robbio-primo-paese-italiano-afare-i-test-sullimmunita-a-tutti-i-cittadini/). 17. 17.Blood tests show 14% of people are now immune to covid-19 in one town in Germany 2020 [cited 2020 April 23]. Available from: [https://www.technologyreview.com/2020/04/09/999015/blood-tests-show-15-of-people-are-nowimmune-to-covid-19-in-one-town-in-germany/](https://www.technologyreview.com/2020/04/09/999015/blood-tests-show-15-of-people-are-nowimmune-to-covid-19-in-one-town-in-germany/). 18. 18.San Miguel County suspends, then recommits to its universal coronavirus testing experiment over 3-day span 2020 [cited 2020 April 23]. Available from: [https://coloradosun.com/2020/04/09/san-miguel-county-telluride-universal-coronavirus-testing/](https://coloradosun.com/2020/04/09/san-miguel-county-telluride-universal-coronavirus-testing/). 19. 19.New York antibody study estimates 13.9% of residents have had the coronavirus, Gov. Cuomo says 2020 [cited 2020 April 24]. Available from: [https://www.cnbc.com/2020/04/23/new-york-antibody-study-estimates-13point9percent-of-residents-have-had-the-coronavirus-cuomo-says.html](https://www.cnbc.com/2020/04/23/new-york-antibody-study-estimates-13point9percent-of-residents-have-had-the-coronavirus-cuomo-says.html). 20. 20.Leung GM, Lim WW, Ho LM, Lam TH, Ghani AC, Donnelly CA, et al. Seroprevalence of IgG antibodies to SARS-coronavirus in asymptomatic or subclinical population groups. Epidemiology and infection. 2006;134(2):211–21. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16490123&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F01%2F2020.04.26.20079244.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000236138600001&link_type=ISI) 21. 21.Shi YL, Li LH, Chu XW, Chen Q, Wang YL, Ma QJ, et al. Seroepidemiological Study on SARS-CoV IgG Antibodies of Different Populations from Several Areas. International journal of biomedical science : IJBS. 2005;1(1):72–5. 22. 22.Yu WC, Tsang TH, Tong WL, Ng TK, Lim W, Yeung HC, et al. Prevalence of subclinical infection by the SARS coronavirus among general practitioners in Hong Kong. Scandinavian journal of infectious diseases. 2004;36(4):287–90. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/00365540410019525&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15198186&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F01%2F2020.04.26.20079244.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000221162100008&link_type=ISI) 23. 23.Degnah AA, Al-amri SS, Hassan AM, Almasoud AS, Mousa M, Almahboub SA, et al. Seroprevalence of MERS-CoV in healthy adults in western Saudi Arabia, 2011–2016. Journal of Infection and Public Health. 2020. 24. 24.Wilder-Smith A, Teleman MD, Heng BH, Earnest A, Ling AE, Leo YS. Asymptomatic SARS coronavirus infection among healthcare workers, Singapore. Emerging infectious diseases. 2005;11(7):1142–5. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid1107.041165&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16022801&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F01%2F2020.04.26.20079244.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000230106600032&link_type=ISI) 25. 25.Wilson N, Kvalsvig A, Barnard LT, Baker MG. Case-Fatality Risk Estimates for COVID-19 Calculated by Using a Lag Time for Fatality. Emerging infectious diseases. 2020;26(6).