Introduction

COVID-19 is caused by the novel SARS-CoV-2 virus [1]. SARS-CoV-2 infection can be clinically manifested in varying degrees of severity, from asymptomatic infection to serious illness, and sometimes death. The first local case of COVID-19 in Israel was detected in March 2020. Since then, Israel experienced two major waves of COVID-19 outbreaks. A total of 413,004 individuals underwent RT-PCR testing from February 21, 2020 until the end of May 2020, from which 17,078 cases were confirmed (a 4.1% positive rate). In May 2020, the influence of the commercial, social, and travel restrictions brought the number of detected cases to less than 15 per day. The second wave started in June 2020. By September 2020, a peak of over 8000 newly confirmed cases were detected daily. From June 1 to September 14, 2020, 1,626,064 individuals underwent RT-PCR testing, and 159,435 new cases were confirmed (a 9.8% positive rate) [2]. Up to September 14, 2020 6,983,391 were infected worldwide with a reported death rate of 4.14% [3].

Under the National Health Insurance law [4] all Israeli citizens must be registered in one of the country’s four Health Maintenance Organizations (HMOs). HMO community healthcare services are widespread throughout the country.

Serologic testing may serve as an internationally important tool for evaluating the genuine rate of all (those who were diagnosed and those who were not) infected individuals in the community regardless of their clinical presentation of the disease (symptomatic or asymptomatic). Individuals with or without symptoms who have not been tested by RT-PCR due to lack of awareness or willingness to be tested, as well as due to limited access to medical services or limited availability of tests, may lead to an underestimation of the extent of the COVID-19 burden in a country. Therefore, serology testing may better evaluate the actual magnitude of infection.

Worldwide, several sero-epidemiological studies were conducted [5,6,7,8,9] utilizing serological tests to reveal the true extent of COVID-19 infection in the general population. The sero-epidemiological surveys conducted thus far included a sample to population ratio of: 0.19% in England, [7] 0.1% in Spain, [5] 0.05% in US [9] and Sweden [8] and 0.026% in Brazil [6]. Compared to other studies conducted thus far, the current study covers much greater proportions of the target population, 0.58% sample to population ratio. This is the first nationwide seroprevalence study conducted in Israel. The vast and well-established data obtained in the current study might contribute to determination of testing strategies and policy decision making in other countries as well.

Methods

Study design and participants

This study was approved by the legal department of the Israeli Ministry of Health. The Israeli first nationwide seroepidemiologic survey began on June 28, 2020 and ended on September 14, 2020. In order to obtain an appropriate representation of the entire population, all municipalities in Israel were characterized according to the following parameters: ethnicity, religious affiliation, socioeconomic status, number of residents, and geographical location (district). The study population included insured individuals who arrived at the HMOs to undergo a blood test for any reason. Due to technical difficulties, the HMOs were only able to provide the number of performed tests and not the number of referred participants. The final sample included all participants that underwent the serology testing, excluding participants with missing data (Fig. 1). For further information, see Online Resource methods.

Fig. 1
figure 1

Flowchart of participants in the seroepidemiological nationwide survey for SARS–CoV–2

Full size image

A 3–5 ml blood sample was collected in a standard vacuum gel tube. Data on the participant's sex, age, and home address were retrieved from the medical file. Data on RT-PCR status, including date performed and result, were obtained from the Ministry of Health Registry.

Detection of SARS-CoV-2 antibodies

Two serologic tests were used, both using qualitative detection of IgG antibodies to SARS-CoV-2: SARS-CoV-2 IgG assay by Abbott® and LIAISON® SARS-CoV-2 S1/S2 IgG by DiaSorin. Specificity and sensitivity of the analytical methods were calculated at the Israeli Central Virology Lab and described in details elsewhere [10].

In order to achieve the highest specificity it was determined that the sample would be analyzed with the Abbott assay and, if found to be positive, would be validated with the DiaSorin assay. Only participants with a positive result in both tests were classified as being seropositive for SARS-Cov-2. The combined specificity after utilizing both tests was 99.99% (95%CI 99.8–100.0). In accordance, due to the double assay test, the combined sensitivity was reduced to 76.7% (95%CI 73.9–80.0).

Statistical analysis

We estimated seroprevalence as being the proportion of individuals who had positive serology results within the survey's population and assessed the effect of selected risk factors and determinants (age, sex, RT-PCR status, time period, and municipality characteristics) on the probability to receive a seropositive result. The serology test result was compared to the RT-PCR result in order to determine the RT-PCR capability as a predictor for receiving a seropositive result. Due to the changes in the epidemiological curve during the survey period, the analysis was divided into two-week periods when relevant. Since all of the municipalities included in the study were divided into strata according to their characteristics, we approached each municipal characteristic as an independent parameter in the univariate analysis.

The chi square was applied for the univariate analysis, when appropriate, and 95% CIs were calculated with the Clopper-Pearson interval. The seroprevalence-associated risk factors were examined with a multivariate generalized linear mixed model (MGLMM). The MGLMM was fitted with age, sex, time period, and RT-PCR status as fixed effects, and the municipal strata were taken as a random effect. The variance component was selected as the random effect covariance. Receiver operating characteristic (ROC) curves were created from which area under the ROC curve (AUC) was calculated as a metric of model performance.

In order to reflect the seroprevalence of the general population of Israel, we used sampling weights for sex, age, municipal strata and RT-PCR status. For the weights procedure, we weighted the age by sex distribution for each municipal stratum to that of the general population. Data regarding the general population was obtained from the Israeli Central Bureau of Statistics. The statistical analysis was performed with the IBM SPSS Statistics for Windows, Version 20.0. Armonk, NY: IBM Corp.

Results

During the recruitment period, 55,302 individuals consented to participate in the survey, and 54,357 of them were included in the final analysis (Fig. 1). The study population consisted of 58.6% males. The mean age of the entire cohort was 45.6 years (± standard deviation 19.83) (Online Resource Table 1). The overall survey seroprevalence was 3.8% (95%CI 3.7–4.0). After weighing for age, sex, and RT-PCR status, the overall seroprevalence was 4.6% (95%CI 4.4–4.8).

Table 1 demonstrates the seroprevalence according to sociodemographic characteristics. Males had a higher prevalence compared to females (4.9% vs. 3.1%, p < 0.001). The highest seroprevalence was among children (5.6%) and adolescents (7.8%). The highest seroprevalence within the adolescent group was among those between 10 and 13 years of age (8.6%, 95%CI 7.1–10.4), with a slightly lower rate among those who were 14–19 y of age (7.4%, 95%CI 6.5–8.4) (Online Resource Table 2).

Table 1 SARS–CoV–2 seroprevalence by population characteristics
Full size table

The seroprevalence decreased as the socioeconomic status increased, from 6.5% in the low socioeconomic municipalities to 1.6% in the high socioeconomic municipalities (p < 0.001). The prevalence among participants living in large municipalities was higher than that in medium and small municipalities.

The seroprevalence increased over time, along with the rise of the epidemic curve during the study period, from 2.7% in the first two weeks to 16% in the last two weeks (Online Resource Fig. 1). Comparison of the percentages of participants who tested positive both by RT-PCR and serology according to the time interval between the two tests revealed the highest seroprevalence rate between four to eight weeks after a positive RT-PCR test (83.8%) (Online Resource Table 3). Between 0 and 4 weeks after a positive RT-PCR result, only 63.6% of the samples had detectible IgG antibodies. Due to the difference in the prevalence according to the time interval between the two tests, it was decided that only RT-PCR tests that had been conducted at least 4 weeks prior to the serology test would be included in further analysis.

A comparison of the seroprevalence according to RT-PCR status revealed that the seroprevalence rate among the participants who underwent RT-PCR testing was higher compared to those who did not, regardless of their RT-PCR result (21.9 vs. 2.1%, respectively, p < 0.001). In total, 74% of the participants who had a positive RT-PCR result were also positive for the serologic test that had been carried out at least 4 weeks after the RT-PCR test (Table 2).

Table 2 SARS-CoV-2 seroprevalence by RT–PCR status
Full size table

We then stratified the data of the individuals who underwent both RT-PCR and serology tests into five periods of two weeks. The percentage of participants that underwent RT-PCR testing increased with time, from 7.6% at the first two weeks of the survey to 28.5% at the last two weeks. At the respective time intervals, the seroprevalence increased from 18 to 35% (p < 0.001). The percentage of people who had RT-PCR positive test results increased from approximately 18% at the first two weeks of the survey to approximately 37% at the last two weeks (Online Resource Table 4).

Table 3 presents the results from a multiple generalized linear mixed model demonstrating the odds ratio (OR) for a positive serology test result. Participants with previous RT-PCR results had an 11-fold higher probability of being serologically positive. Males had 1.5 higher odds (OR = 1.56, p < 0.001) to be seropositive in comparison to females. For adolescents aged 10-19 years, the odds to be seropositive increased by 3-fold in comparison to adults 70 years of age or more (OR = 3.3, p < 0.001). Participants who were tested during the last two weeks of the current survey had increased odds (Adj OR = 5.4, p < 0.001) to be seropositive in comparison to participants that underwent serology testing during the first two weeks of the survey. In a ROC analysis, the AUC was 0.87 (p < 0.001, 95%CI 0.86–0.88), demonstrating high compatibility of the model.

Table 3 A multivariate generalized linear mixed modela for a positive serology test result
Full size table

In order to examine the effect of a given municipality’s socioeconomic status as a risk factor for a positive serology test result, we stratified the data by socioeconomic status and used a multivariable logistic analysis. Neither age nor sex differ when a high socioeconomic status was compared to medium and low socioeconomic statuses (Online Resource Table 5).

Figure 2 demonstrates the ratio between the percentage of positive serology test results and positive RT-PCR test results at the corresponding time point (at least 4 weeks prior to the serology test date). For each time period, the cumulative seroprevalence was calculated for all participants that underwent the serological test. In addition, the incidence of SARS-CoV-2 in the general population, as was verified by RT-PCR testing, was calculated for each time period. After weighting for age, sex, and district of the survey's population to the respective characteristics of the entire population that underwent RT-PCR testing, a positive correlation was found between the incidence of SARS-CoV-2 and seroprevalence (r = 0.94, p = 0.017). The ratio between the seroprevalence and the incidence was 15.7 (3.3 vs. 0.2%, respectively) during the first two weeks of the survey. That ratio decreased over time to a ratio of 4.6 during the last two weeks of the survey (4.6 vs. 1.0%, respectively)

Fig. 2
figure 2

Seroprevalencea among survey participants and positive RT–PCR results in the general population by two–week intervalsb. aAdjusted for age, sex, and district. bFor each time period the RT–PCR test date is at least 4 weeks prior to the serology test date

Full size image

Discussion

This nationwide Israeli survey is one of the largest SARS-CoV-2 seroprevalence studies conducted thus far throughout the world. Since we predicted a low seroprevalence rate, we favored greater specificity over sensitivity. This goal was achieved by combining two different test kits. Our results indicated a seroprevalence of 3.8% (weighted 4.6%). Since the majority of participants (89.7%) underwent the serology test up to August 1, 2020 (the study closure was September 14, 2020), the results mainly reflect the end of the first outbreak until the middle of the second outbreak in Israel. Our findings are important for establishing a baseline level of seroprevalence for future investigations.

While most studies reported no difference in seroprevalence between sexes [5, 6, 11, 12] our finding suggested a 1.6-fold (95%CI 1.4–1.7) higher prevalence in males compared to females. In agreement with our results, Rosenberg et al. [13] and Iversen et al. [14] also reported higher seroprevalence in males. In contrast, other studies reported a higher seroprevalence in females vs males, [15, 16] so the issue remains controversial.

In agreement with other publications, [6, 11, 15] our results indicated differences in prevalence between age groups. We found the highest prevalence rates, both adjusted and unadjusted, in the group aged 10–19 years (7.8%). These findings may be explained by the higher numbers of social interactions among adolescence in comparison to other age groups. Children and adolescents are less likely to undergo routine blood tests, therefore, their representation was relatively low in the study (10.8%) compared to their proportion in the Israeli general population (36%). However, the survey included 5864 children and adolescents, which is sufficient to draw a representative and reliable conclusion. In agreement with other studies, [5, 7, 8, 11] the lowest seroprevalence was found in participants over the age of 60 years. Since the elderly comprise the major risk group for COVID-19 complications, the low prevalence may result from their preference to minimize social interactions. Alternatively, the elderly immune system might be suppressed, whereupon this age group might represent lower or even undetectable levels of antibodies [17].

A municipality’s socioeconomic status was highly associated with seropositive results, with a four-fold higher prevalence among low compared to high status. Others have also reported high seroprevalence within populations with low socioeconomic characteristics, such as lower levels of education, [18] lower income, and big household size. [6, 7] Greater exposure to the SARS-CoV-2 virus may be related to higher population density, crowded households, and lack of awareness, all of which result in greater difficulty in observing social distancing. In Israel, ultraorthodox municipalities are characterized by a relatively low socioeconomic status, high population density, and a very unique community structure that makes social distancing difficult to follow, all of which may explain the high seroprevalence detected in ultraorthodox municipalities.

Recent studies have indicated that SARS-CoV-2 IgG antibodies titer continues to rise for three to four weeks after symptom onset [19,20,21,22,23]. In our study, the highest seroprevalence was detected four to eight weeks after a positive RT-PCR test result. In accordance, only when a RT-PCR test took place at least four weeks prior to a serology test was it considered as having been performed.

Participants who had undergone RT-PCR testing had an over 10-fold higher likelihood to be seropositive, regardless of their RT-PCR result. During the entire survey, the criterion for undergoing RT-PCR testing in Israel was suspicion for COVID-19, which might explain the higher seropositive rate among the RT-PCR-tested participants. Participants with negative RT-PCR results had a three-fold higher likelihood to be seropositive compared to participants who did not undergo RT-PCR testing. The same trend was observed by others [5, 7]. Ward et al. [7] reported similar results, indicating that individuals with a suspected case of SARS-CoV-2 infection (both by the doctor or by self-report) had a 8.5-fold higher probability to be seropositive than individuals with no suspicion.

In the current cohort, 6.2% of the individuals with negative RT-PCR test results were seropositive. This is more likely to be explained by the timing of the negative RT-PCR test rather than receiving false negative RT-PCR test result. Only 74% of the people with previous positive RT-PCR test results were also seropositive. This may be explained by the reduced combined sensitivity of the test kits (Abbott and DiaSorin), by the reduction in the antibodies titer over time, or by individuals that did not develop antibodies despite of being exposed to SARS-CoV-2 virus [10]. In agreement, several studies reported that asymptomatic carriers are less likely to develop antibodies in comparison to symptomatic patients [23,24,25]. Another explanation might be the reported rapid decay and short half-life of SARS-CoV-2 antibodies [26, 27].

Israel experienced two COVID-19 outbreaks between the beginning of March 2020 and the end of September 2020. The second wave was prolonged in comparison to the first one and characterized by higher numbers of detected cases. Therefore, both the incidence and seroprevalence rates increased over time.

In the current study, seroprevalence indicated 4.5- to 15.7-fold higher magnitude of infection than that identified by RT-PCR testing. Bendavid et al. [28] reported an estimated 22- to 95-fold higher seroprevalence compared to incidence, and Havers et al. [29] estimated a ratio of 6- to 24-fold seroprevalence over incidence. The relatively low seroprevalence-to-incidence ratio detected in the current study may be due to the high availability and accessibility of RT-PCR testing in Israel [3]. At the early stages of the outbreak, the ratio between the two parameters was the highest and it decreased over time. This trend might be explained by the greater availability and accessibility of RT-PCR testing to the Israeli public over time, resulting in a greater number of identified cases. It can be assumed that affording more individuals the possibility of undergoing routine RT-PCR tests would minimize the rate of undetected cases, which would otherwise be revealed only by serology.

Similarly to the sample collection design applied in other seroepidemiologic studyies, [9, 30] the population in the current survey is comprised of people undergoing blood tests for any reason. Therefore, despite the large sample size tested, this study might not reliably represent the entire population. However, there is no evidence indicating a higher susceptibility for COVID-19 among individuals undergoing blood tests in HMOs.

Serologic testing is important for evaluating both identified and unidentified infection. It may also be used as an intervention measure, as had been suggested by others [31]. In Israel, serologic test was recognized as a sufficient tool to define a person as having recovered from SARS-CoV-2 infection [32]. Recovery status exempts a seropositive person from 14 days of isolation in cases of close contact with a confirmed SARS-CoV-2 carrier, or after entering the country from abroad. Furthermore, serological diagnosis may support the deployment of employees with a positive serology test to ensure the maintenance of a stable and functional economy. Seroprevalence surveys among healthcare workers and other essential subpopulations are crucial for maintaining a functioning economy in times of COVID-19 outbreaks and may contribute both on personal and national levels.

On December 20, 2020 a wide vaccination program was launched in Israel, aiming at vaccinating the entire population against SARS-CoV-2 within a few months. The number of vaccines already purchased by the Israeli government along with the relatively small population size of the country may enable Israel to be the first country worldwide to achieve herd immunity against COVID-19.

While vaccinations against SARS-CoV-2 are becoming more available worldwide, the use of serology as a tool to manage the pandemic may shift from being a diagnostic tool to a vaccination prioritizing tool. For instance, individuals with a seronegative result may be prioritized over individuals with a seropositive result to achieve a more efficient use of vaccinations and gain "herd immunity" earlier.

In conclusion, our results indicate that the population of Israel is still far from being protected against SARS-Cov-2 by "herd immunity". Additional nationwide surveys are warranted to evaluate the effect of further outbreaks on the seroprevalence in Israel. Consecutive periodical surveys will make it possible to monitor SARS-CoV-2 infection in Israel over time. The findings of this study provide evidence-based data for public health decision-making, not only at the national level but also worldwide. The seroprevalence-to-incidence ratio emphasizes the benefits of serology testing, mainly where RT-PCR tests are less available or accessible to the population.