Based on the overall lower amino acid sequence identity (Supplementary Table 1) and on our previous work with SARS-CoV-2 10, spike protein subunit 1 (S1) was selected as the antigen to set up an enzyme immunoassay (EIA) for the detection of HCoV S binding IgG antibodies. S1 subunits were produced as mFc-fusion (S1-mFc) proteins in HEK-293F cells, purified (Supplementary Fig. 1) and used in EIA.
Seropositivity and HCoV spike-specific IgG antibody levels in children of 1 to 3 years of age
Serum specimens collected in 2009-2013 from 140 healthy children at 1, 2, and 3 years of age were analyzed by EIA for antibodies against four seasonal HCoV (229E, HKU1, NL63, OC43), MERS and SARS-CoV-2 S1 proteins. The absorbance values were converted to EIA units enabling reliable identification of seropositive samples and to compare antibody levels against different coronavirus species (Fig. 1). An increase in the geometric mean antibody levels between 1- and 2-year samples was significant for all four seasonal HCoVs (p <0.0001) while the change between 2 and 3 years was significant only for NL63 and OC43 (p <0.0001 and p =0.0004, respectively). One participant had low levels of MERS and SARS-CoV-2 S1-binding IgG antibodies at 2 and 3 years of age.
The rate of antibody positive children increased by age for the four seasonal HCoVs and the antibody levels remained elevated in most of the seropositive children (Table 1). In all age groups the seropositivity was the highest for OC43 with a seropositivity of 31% in 1-year-olds and 81% in 3-year-olds. At 3 years of age the seropositivity rate for 229E (37%) was lower than for the other seasonal HCoVs (59% for HKU1, 76% for NL63, and 81% for OC43). It is noteworthy that the cumulative seropositivity for HKU1 declined between 2 and 3 years of age, while in 229E, NL63 and OC43 the annual and cumulative seropositivity rates increased almost at the same rate (Table 1). A decline in anti-HKU1 S1 antibodies was observed in previously seropositive children whose sequential samples showed a decrease in geometric mean IgG antibody levels (GMALs) from 2 years to 3 years (44 to 33 EIA units). For 229E, NL63, and OC43 the mean IgG antibody levels in seropositive children remained relatively stable during the 3-year follow-up.
A subset of seropositive 1-year-old or 2-years-old children showed a diagnostic increase in HCoV S1-specific antibody levels in subsequent samples at 2 or 3 years of age, respectively, indicating a likely reinfection or re-exposure (below referred to as reinfection). Altogether 31% (27/88) of the children seropositive for any HCoV at 1 year showed a diagnostic/significant increase (more than 20 EIA units between sequential samples) in anti-S1 HCoV antibodies by the age of 2 years. A similar rate of likely reinfection for any HCoV (34%) was observed between 2 and 3 years for children who were seropositive at the age of 2 years.
Correlation of anti-HCoV S1 and anti-HCoV N antibody levels
HCoV infections have been shown to induce antibodies for both HCoV S and N proteins 9,10. We have previously analyzed anti-HCoV N IgG antibody levels for the same set of serum samples from children (n=420)18 allowing us to use the previous data to analyze the correlation of N and S1 protein-specific antibody responses. We calculated the correlation coefficients for S1 and N EIA data of each HCoV (Fig. 2). The correlation coefficients of S1 and N IgG antibody levels for 229E and NL63 were relatively high (r=0.66 and r=0.65, respectively, p <0.0001). The rates of correlation for anti-HKU1 S1 and N, and anti-OC43 S1 and N IgG antibodies were somewhat lower, yet significant (r=0.44 and r=0.55, respectively, p <0.0001). However, the negative samples may contribute considerably to the higher rates of correlation since the removal of samples that were negative in both assays weakened the correlation of all assay pairs (Supplementary Fig. 2).
Homologous and cross-reactive S1 antibodies in immunized rabbits and guinea pigs
To validate further the antigens chosen for EIA, we immunized rabbits and guinea pigs with the HCoV S1 without mFc-fusion (except for OC43 for which S1 with mFc was used due to unsuccessful production of OC43 S1 without mFc, Supplementary Fig. 1). Antigens without mFc-fusion were used in immunization to avoid formation of antibodies against the mFc domain. The specificity of immune sera was examined by immunofluorescence (IF) assay for homologous and cross-reactive antibodies using 229E, OC43, and NL63 virus-infected cells (Huh7 or LLC-MK2). Virus-infected cells were detected with high IF titers (1:1000 to 1:5000) by the corresponding anti-S1 sera of immunized rabbits and guinea pigs (Table 2, Supplementary Fig. 3). In addition to homologous reactivity of sera, cross-reactivity was also observed. Cells infected with 229E were recognized by anti-HKU1 S1 rabbit serum (titer 50) and by both rabbit and guinea pig anti-NL63 S1 sera (titer 200; Supplementary Table 1, Supplementary Fig. 4). OC43-infected cells were recognized by rabbit and guinea pig anti-HKU1 S1 sera (titer 200, Supplementary Fig. 5) and NL63-infected cells were recognized by rabbit anti-OC43 S1-mFc sera (titer 200, Supplementary Fig. 6).
COVID-19 vaccination has no effect on HCoV S1 IgG antibody levels in adults
Results from HCoV S1 immunized animals suggested potential weak cross-reactivity between HCoV-specific antibodies. To investigate whether COVID-19 vaccination induces the production of HCoV S1 cross-reactive antibodies, we analyzed sera of COVID-19 vaccinees with HCoV S1 protein-specific EIAs. Sera were collected from 113 HCWs (25–65 years old) before vaccination with BNT162b2 (0D; August 2020 to January 2021), three weeks after the second BNT162b2 dose (2D; January to March 2021), and three weeks after the third vaccination dose with BNT162b2 or mRNA-1273 (3D; October 2021 to January 2022). Geometric mean IgG antibody levels (GMALs) for the seasonal HCoV S1 proteins decreased from 0D to 2D but increased from 2D to 3D (44 EIA units at 0D vs. 41 EIA units at 2D vs. 44 EIA units at 3D for 229E, both p < 0.0001; 35 vs. 33 vs. 36 for HKU1, p = 0.0007 and p = 0.0008; 50 vs. 46 vs. 50 for NL63, p < 0.0001; 43 vs. 41 vs. 46, for OC43 p < 0.0001; Fig. 3). GMALs for MERS S1 remained practically negative while two doses of BNT162b2 increased the SARS-CoV-2 S1 GMALs from negative to high levels (GMAL of 1.6 EIA units at 0D vs. 105 at 2D, p < 0.0001). The third mRNA vaccine dose (BNT162b2 or mRNA-1273) further boosted the SARS-CoV-2 antibody levels (GMAL of 105 EIA units at 2D vs. 113 at 3D, p < 0.0001).
Prevalence of HCoV infections in 2020 – 2021
To investigate how the circulation of different HCoVs associate with the changes in HCoV antibody levels in the study population, we analyzed the monthly number and prevalence of HCoV PCR-positive samples in hospitalized patients in Southwest Finland health district in 2020 – 2021 (Fig. 4). The multiplex RT-qPCR with Allplex Respiratory Panel 3 (Seegene Inc.), which detected 229E, NL63, and OC43 (but not HKU1) RNA, showed the circulation of these viruses in early 2020 (January – March) with a peak of 12% of positive HCoV samples in March 2020 (7 229E-positive, 15 NL63-positive, and 19 OC43-positive out of 349 tested samples). This was followed by 12 months (April 2020 to March 2021) of very low number of positive samples (only 4 NL63-positive samples out of 2356 tested samples). Starting from April 2021 the circulation of OC43 was observed monthly until the end of the study period. NL63 was detected from May to July 2021 but no circulation was seen later in the year. 229E was detected sporadically after May 2021 and occasionally between October and December 2021. Detection of SARS-CoV-2 with laboratory developed RT-qPCR showed higher number of positive cases starting from autumn 2020 with peaks in March and December 2021.
Congruence of seasonal HCoV occurrence and S1 IgG antibody levels
The absence of serologically observable seasonal HCoV reinfections (defined as >20 EIA unit increase in antibody levels of sequential samples) between 0D and 2D (Fig. 3a-d) correlated with the low number of seasonal HCoV-positive samples in Autumn 2020 to March 2021 in Southwest Finland health district. The increases in seasonal HCoV S1 IgG antibody levels between the second and third vaccine doses (8 months apart) on the other hand suggested potential circulation of HCoVs and reinfection amongst HCWs (Fig. 5). The PCR-data on the rates of 229E, NL63, and OC43 detections from May 2021 onwards matched well with the serological changes among some of the HCWs on the same period. Diagnostic rises in antibody levels from 2D to 3D indicated a reinfection rate of 5.3% (6/113) for 229E, 6.2% (7/113) for HKU1, 3.5% (4/113) for NL63, and 14.2% (16/113) for OC43. Interestingly, each participant (n=7) showing an increase in HKU1 S1 antibodies had also an increase in OC43 S1 antibodies (Fig. 5).
Correlation of HCoV S1-binding antibody levels between HCoVs
To evaluate the presence of homologous and potentially cross-reactive anti-S1 antibodies, we used the data from 339 serum samples of COVID-19 vaccinated HCWs to analyze the correlation of seasonal HCoV S1-binding IgG antibodies. The highest rates of correlation were observed for anti-229E and anti-NL63 S1 antibodies (r=0.524, p<0.0001), and anti-HKU1 and anti-OC43 S1 antibodies (r=0.631, p<0.0001) but also other assay pairs showed moderate correlation coefficients (r>0.334, p<0.0001 for other anti-HCoV S1 antibody level pairs) (Fig. 6).
In children (420 serum samples, Supplementary Fig. 7), a moderate correlation was observed for anti-HKU1 S1 and anti-OC43 S1 binding antibodies (r=0.589, p<0.0001) and for anti-229E S1 and anti-NL63 S1 binding antibodies (r=0.512, p<0.0001) while the correlation coefficient values of other anti-HCoV S1 pairs was lower (r<0.26). Pairwise comparison of S1 amino acid sequences (Supplementary Table 1) showed that the pairs of 229E and NL63, as well as HKU1 and OC43 shared more identical amino acids (50% and 59%, respectively) in comparison to other sequence pairs (10-22%), which may contribute to sequence identity-related immunological correlation.
The correlation coefficient values for anti-SARS-CoV-2 S1 antibodies and other anti-HCoV S1 antibodies after COVID-19-vaccination were low at 2D and 3D time points (Supplementary Fig. 8). Despite the low correlation coefficient values, likely due to mainly seronegative specimens, the correlation of anti-SARS-CoV-2 S1 antibodies was statistically significant (p<0.05) with anti-OC43 antibodies at 2D (r=0.3033, p=0.0011), and with anti-MERS antibodies both at 2D and 3D (r=0.2462, p=0.0086 and r=0.2098, p=0.0257, respectively).