Seroprevalence of SARS-CoV-2 IgG antibodies in two regions of Estonia (KoroSero-EST-1)

Study design

The KoroSero-EST-1 was a cross-sectional seroepidemiological study conducted from May 8 to July 31, 2020 in two general practitioner (GP) practices in capital Tallinn in Harju county and in Saaremaa with total number of patients 13,260 and 7,525, respectively.

Participants were selected using stratified random sampling. Strata were formed by classifying patients from each GP into 10-year age groups, except the age group of 80 years or older due to their small number. From each stratum individuals were randomly identified by Estonian Health Insurance Fund with the aim to include at least 110 participants per age group from both GP practices to achieve desirable precision for the seroprevalence estimates.

Participants filled in a questionnaire based on World Health Organization recommendations [9] that included the presence of acute respiratory illness since March 1, 2020 and if present, its symptoms, known contact with PCR-confirmed COVID-19 cases, previously performed PCR test for SARS-CoV-2 and its result (Supplementary Information).

All participants and/or their legal guardians provided written informed consent. The study was approved by the Research Ethics Committee of the University of Tartu.

Sampling and antibody measurements

Blood samples (3.5 mL) were drawn and stored for 48 hours at +4°C until transported to laboratory, where sera were stored at -30° C until testing in SYNLAB Estonia Central Laboratory in Tallinn or at the research laboratories of the University of Tartu, Estonia.

First, all samples were tested by chemiluminescent microparticle immunoassay for detection of IgG against SARS-CoV-2 nucleoprotein (N) (Abbott Architect SARS-CoV-2 IgG with ARCHITECT i2000SR analyzer; Abbott Laboratories, USA) according to the manufacturer’s instructions. Samples close to cut-off value (signal/cut-off ratio 0.3 to 1.39) were additionally tested by chemiluminescence immunoassay for detection of IgG against S1 and S2 proteins (LIAISON® SARS-CoV-2 S1/S2 IgG test; DiaSorin S.p.A., Italy) to increase sensitivity without influencing specificity [10].

Next, sera positive for anti-N or anti-S1/S2 IgG or those from participants with previously positive PCR test and two age- and sex-matched negative controls or a random subset of them (Figure 1) were tested for anti-N IgG by a rapid Standard Q COVID-19 IgM/IgG Duo test (SD BIOSENSOR, Republic of Korea) according to the manufacturer’s instructions and by in-house luciferase immunoprecipitation system (LIPS) assays as described elsewhere [11].

• Download figure
• Open in new tab

Figure 1.

Flowchart of recruitment of study participants, testing and definition for being considered seropositive

Finally, an in-house neutralization assay was performed in duplicates for all PCR-positive samples and samples positive for anti-N or anti-S1/S2 IgG. Briefly, patients’ sera were two-fold serially diluted in Virus Growth Media (VGM – DMEM/0.2% BSA/Pen-Strep) starting from 1:4 to 1:4096 in a volume of 50 μl per well. The viral stock was diluted in VGM to contain 100 pfu (plaque forming units) per 50 μl. 100 plaque forming units of the virus stock per well (a local Estonian SARS-CoV-2 isolate from a COVID-19-positive patient’s nasal swab; from SYNLAB) was added to the wells containing sera dilutions, and the mixtures were incubated for 1 hour at 37°C. After this, 4×10⁴ Vero E6 cell suspension in 100 μl of VGM per well was added to the wells. Plates were incubated at 37°C and 5% CO2 for up to 5 days. Neutralization titer was calculated by microscopically evaluating cytopathic effect (CPE) in infected wells (round/floating cells compared to untouched monolayer in non-infected wells). The dilution of sera that did not show CPE was considered as neutralization titer.

Statistical analysis

The software program Stata 14.2 was used for calculating seroprevalence estimates and infection fatality rates and R (version 3.6.1) for other statistical analyses.

Total cumulative number of confirmed COVID-19 cases per 10,000 inhabitants was calculated as of two weeks before the midpoint of sample collection in GP practice. Daily numbers of confirmed COVID-19 cases (positive PCR tests for SARS-CoV-2 from nasopharyngeal swabs) in each age group in Harju county (assumed to reflect the spread of SARS-CoV-2 in Tallinn; 72.3% of Harju county inhabitants live in Tallinn) and Saaremaa were obtained from Health Board [2] and the number of inhabitants as of January 1, 2020 from Statistics Estonia.

Sample size calculation

Sample size of at least 110 participants per stratum was needed to obtain binomial exact 95% CI with width of 5%, assuming seroprevalence was 1% or to estimate seroprevalence with precision ±10%, assuming seroprevalence was more than 10%. Assuming non-response of about 50% (up to 70% in children), 2341 individuals were sampled from both GP practice to whom invitations were sent until about 110 per stratum agreed.

Seroprevalence estimation

Seroprevalence was estimated as the proportion of participants considered seropositive, i.e. with positive result in at least two different antibody assay out of four assays performed (Figure 1), overall and in each age group, separately for the GP practices. To account for stratification, for each age group in a GP practice sampling weight was used that was further adjusted for response probability by post-stratification. Due to large sampling fractions, finite population correction factor was applied. Binomial exact 95% CI were calculated for seroprevalences in age groups.

Seroprevalence adjusted for test performance of anti-N IgG assay (Abbott Architect SARS-CoV-2 IgG with ARCHITECT i2000SR analyzer) was calculated as where r is sensitivity, s specificity and q denotes prevalence [12]. Nonparametric bootstrap confidence intervals were calculated for adjusted seroprevalence estimates by drawing 10,000 bootstrap samples at random from our dataset and binomial distribution of the assay validation dataset. Each bootstrap sample was adjusted for sensitivity of 92.7% and specificity of 99.9% [13] and 2.5th and 97.5th percentiles of the bootstrap distributions were reported. Infection fatality rate was calculated as number of reported death cases related to COVID-19 as of July 31, 2020 divided by the estimated total number of seropositive individuals. Confidence intervals were calculated assuming that number of deaths follow Poisson distribution.

Analysis of symptoms

The association between the presence of symptoms and seropositivity was analyzed by multiple correspondence analysis (MCA) using R package FactoMineR. Only dimensions that explained >1/Q, where Q is the number of variables used, of the total inertia were retained for interpretation. Association between seropositivity and symptoms was determined from dimensions with high contribution and quality (cos²>0.1) of seropositivity. Hierarchical clustering using Euclidean distance and Ward’s method followed by k-means consolidation was performed on all dimensions of MCA to partition the individuals into groups. The number of clusters was based on elbow of the barplot of the gain in within-inertia.

For comparing categorical variables, Fisher’s exact test was used.

Analysis of antibody concentrations/titer

Multiple linear regression was used to analyze the association between IgG and neutralizing (log-transformed and undetectable titers equaled to 0) antibodies concentrations/titer and the presence of acute respiratory illness or its symptoms or hospitalization adjusted for age, sex and time since onset of symptoms, shown to influence antibody concentrations [5, 14].

Multiple imputation was used to replace missing values of the dates of time of onset of symptoms with possible dates of infection in seropositive individuals without acute respiratory illness. Total of 1000 datasets were created where for each individual without acute respiratory illness the date was imputed from the dataset of confirmed COVID-19 cases [2] within the same county and in the inhabitants of the same age group detected until 14 days prior to drawing the individual’s blood sample. Coefficients of multiple linear regression models and their standard errors were calculated by Rubin’s rules [15].

Study population

In total 3608 subjects were invited to participate of whom 1960 (54.3%) agreed and provided blood samples (Figure 1); all except one filled in the questionnaire. Of the participants 809 (41.3%) were male, 182 (9.3%) had been in contact with confirmed COVID-19 case, 221 (11.3%) had had PCR test for SARS-CoV-2 and 157 (8.0%) had had acute respiratory illness (Table 1). During the study period, the spread of SARS-CoV-2 had largely ended in Saaremaa, while in Tallinn it was ongoing (Figure S1).

Laboratory testing

In total, 73 participants were seropositive. Of those all had detectable concentrations of anti-RBD by LIPS assay, 71 (97.3%) anti-N IgG by Abbott or anti-S1/S2 by DiaSorin assay, 68 (93.2%) anti-N IgG by LIPS assay, 60 (82.2%) neutralizing antibodies, 44 (60.3%) anti-N IgG by Biosensor assay (Figure 2). All 146 samples from seronegative age- and sex-matched controls were negative for anti-N IgG by Biosensor assay; of 22 samples tested by LIPS assays, one was positive for anti-N IgG.

• Download figure
• Open in new tab

Figure 2.

Upset plot of number of seropositive individuals with positive test results in different assays. Each bar shows number of individuals with positive results from the assays indicated by black circles beside the plot. Bars are colored according to number of individuals who had had acute respiratory illness (red) or not (green). Circles in the bottom show for each set of individuals negative (transparent) or positive (black) test result by the respective assay. Solid lines connect black circles.

Seroprevalence

Overall seroprevalence was 1.5% (95% CI 0.9-2.5) in Tallinn and 6.3% (95% CI 5.0-7.9) in Saaremaa (Table 1). There were no difference in seroprevalence between males and females, age groups or household sizes. Seroprevalence was higher in the participants with contact with confirmed case or acute respiratory illness compared with those without. The seroprevalence estimates adjusted for test performance were similar to the estimates not adjusted for test performance, both overall (1.5% (95% CI 0.6-2.4) in Tallinn, 6.7% (95% CI 5-8.4) in Saaremaa) and in age groups (Table S1).

Number of deaths due to COVID-19 was 7 in Tallinn and 26 in Saaremaa and the respective infection fatality rates were 0.1% (95% CI 0.0-0.2) and 1.3% (95% CI 0.4-2.1).

Seroprevalence vs confirmed cases

The cumulative number of confirmed COVID-19 cases in Tallinn (as of June 10, 2020) was 11.3 and in Saaremaa (as of June 15, 2020) 167.2 per 10,000 inhabitants. Seroprevalence of COVID-19 in Tallinn was 13.3 (range 8.0-22.1 based on 95% CI of seroprevalence estimate) and in Saaremaa 3.8 (3.0-4.7) times greater than the cumulative incidence of COVID-19 reported in the respective regions.

Of individuals with contact with confirmed case 43.0% (95/221) and of those with acute respiratory illness 33.1% (52/57) had PCR test performed. Considering participants without PCR test performed, in those with contact with confirmed case the proportion of seropositive individuals was 13.5% (17/126) and in those without contact 1.5 % (29/1612) (p<0.001); in those with acute respiratory illness the proportion of seropositive individuals was 6.7% (7/105) and in those without illness 2.4% (39/1633) (p=0.018).

Symptoms

Of seropositive individuals 19.2% (14/73) and of seronegative individuals 7.6% (143/1886) self-reported having had acute respiratory illness (p=0.001) with median (interquartile range) of 6 (5-7) and 5 (3-7) symptoms, respectively. The frequency of different symptoms was similar in seropositive and seronegative individuals with acute respiratory illness, but seropositive individuals required more often hospitalization due to the acute respiratory illness compared with seronegative individuals (Table 2).

In MCA of all individuals with acute respiratory illness, the first two dimensions discriminated poorly seropositive and seronegative subjects (Figure S2). Seropositivity had high quality on dimensions three and four and was associated with the absence of cough and the absence of sore throat (Figure S2). In MCA of 102 participants aged less than 50 years, of whom 5 were seropositive, it was not possible to associate seropositivity with any symptom by the first 6 dimensions due to low quality and contribution of seropositivity on dimensions 1-4 or symptoms on 5-6. In MCA of 55 participants aged 50 years or more, of whom 9 were seropositive, according to the first two dimensions of MCA seropositivity was associated with the absence of cough and runny nose and the presence of diarrhea and fever (Figure 3) and majority (7/9) of seropositive individuals were grouped into a single k-means cluster containing no seronegative individuals (Figure S3).

• Download figure
• Open in new tab

Figure 3.

Biplot of the first two dimensions of multiple correspondence analysis on association between the presence (filled large dots) or absence (empty large dots) of symptoms and seropositivity (crosses) in all individuals aged 50 years or more with acute respiratory illness. The variables included in the analysis are colored according to the sum of cos² of the variable on dimensions one and two. Grey and black small dots represent seronegative and seropositive individuals, respectively. Practices of gravity of seropositive and seronegative individuals are shown by squares surrounded by 95% confidence ellipses.

Antibody concentrations/titer

According to multiple linear regression, anti-N IgG and anti-RBD IgG concentrations by LIPS assay and anti-N IgG by Abbott assay were significantly higher in seropositive individuals with fever (Table 3). Anti-RBD and anti-N IgG by LIPS assay were significantly higher in hospitalized patients, anti-N IgG by LIPS and/or Abbott assay in patients with difficulty breathing, shortness of breath or chest pain and anti-N IgG by Abbott assay in patients with diarrhea. There were no statistically significant associations between titer of neutralizing antibodies and characteristics of acute respiratory illness.

Seroprevalence

The seroprevalence estimates found in our study are within the range of 0.4-10.9% reported in other population-based seroepidemiological studies [4-8] and far from the threshold of 60 to 70% required for herd immunity [16]. Considering the number of cumulative cases per 10,000 inhabitants of 167.2, we would have expected larger seroprevalence in Saaremaa, as seroprevalence estimates were between 4.6-10.9% in Spain, United Kingdom and Geneva [4, 6, 7] with smaller number of cumulative cases per 10,000 (44.7-103.2) and thereby larger infection-to-case ratios 10.6-14.3 [1, 4]. However, these seroepidemiological studies were conducted shortly after the peak of the spread when the daily incidence curve had flattened at low levels, but the spread was still ongoing [1, 4, 6, 7], suggesting the presence of undocumented cases and thereby underestimation of cumulative cases per 10,000 inhabitants. Similar pattern of ongoing spread during sample collection for seroepidemiological study and infection-to-case ratio 11.3 was in Tallinn. In contrast, in Saaremaa the spread of SARS-CoV-2 had been successfully contained and only few cases were detected during the study period, a situation similar to Reykjavik and Luxembourg, where the cumulative cases per 10,000 inhabitants were 50.4 and 53.3 and seroprevalence estimates 0.4% and 2.0%, respectively [5, 8]. Such discrepancies between infection-to-case ratios underline once again the importance of seroepidemiological studies in estimating the actual prevalence of infection. It is believed that incidence of SARS-CoV-2 infection was severely underreported in spring not only in Estonia but elsewhere [17]. This is further supported by the age distribution of confirmed cases with predominance of older age groups (46.1% and 68.4% aged 50 years or more in Tallinn and Saaremaa, respectively [2]) among reported cases while age distribution in our study was even in our seroepidemiological study.

Symptoms

Of individuals with acute respiratory illness, two thirds had not been tested by PCR, but they were significantly more likely to be seropositive, showing underdetection of cases among them. Detection of cases by individual symptoms has low specificity [18, 19], particularly during the season of respiratory infections [18], but this could be increased by a combination of them [20]. Small number of individuals with acute respiratory illness impeded detection of pattern in our study. Although the proportion of asymptomatic cases among PCR-positive individuals has been up to 87.9% [21], we may have missed some symptomatic cases due to focusing only on acute respiratory illness. However, the proportion of COVID-19 cases presenting without any respiratory symptoms, such as solely by fever or gastrointestinal symptoms, is low, up to 21% [22]. Furthermore, it is notable that while the proportion of fever, the most specific symptom of COVID-19 [18], in cases detected by PCR is larger than that of cough [18, 23], in seropositive individuals fever is mostly less common or occurs with rate comparable to cough [8, 24]. This may partly corroborate that the proportion of symptoms in seropositive cases attributable causes other than COVID-19 may be notable, up to 31.5% [25]. Additionally, this may suggest underdetection of cases among individuals with respiratory infections presenting without fever, supported by smaller proportion of people with cough being tested compared to those with fever [26], warranting focusing particularly on respiratory infections, as in our study. In adults aged 50 years or more, we found that the symptoms associated with seropositivity were the absence of runny nose and cough and the presence fever and diarrhea, similarly to a previous study [27], but the lack of information on anosmia or ageusia that are most strongly associated with COVID-19 [24, 26] may have compromised our analysis. Nevertheless, our results suggest that focusing only on respiratory symptoms may delay accurate diagnosis of SARS-CoV-2 infection in older people.

Antibody concentrations/titer

In line with several other studies, we found that hospitalized patients or patients with symptoms suggesting more severe disease, such as fever, or shortness of breath had higher anti-N IgG concentrations [5, 14, 28]. However, these factors were not associated with titer of neutralizing antibodies or were less often associated with anti-RBD IgG that contain large proportion of neutralizing antibodies [29], in contrast with other studies [14, 30]. While these studies are mostly conducted on samples collected within two months after development of symptoms, long-term response of neutralizing antibodies is scarcely described [28]. Notably, in contrast to anti-RBD and anti-N IgG, neutralizing antibodies decline to undetectable titer two months after the onset of symptoms in 20-25%of patients, as in our study, despite development of detectable peak in nearly all of them [28]. Moreover, titers of neutralizing antibodies decline faster compared with anti-RBD and anti-N IgG [14, 28]. Therefore, no associations found for neutralizing antibodies could be attributable to such faster decline and possible attainment of a plateau concentration where differences are less clear, as suggested [28], but further studies on longevity of neutralizing antibodies are warranted.