ABSTRACT
Sero-monitoring provides context to the epidemiology of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections and changes in population immunity following vaccine introduction. Here, we describe results of a cross-sectional hospital-based study of anti-spike seroprevalence in New York City (NYC) from February 2020 to July 2022, and a follow-up period from August 2023 to October 2023. Samples from 55,092 individuals, spanning five epidemiological waves were analyzed. Prevalence ratios (PR) were obtained using Poisson regression. Anti-spike antibody levels increased gradually over the first two waves, with a sharp increase during the 3rd wave coinciding with SARS-CoV-2 vaccination in NYC resulting in seroprevalence levels >90% by July 2022. Our data provide insights into the dynamic changes in immunity occurring in a large and diverse metropolitan community faced with a new viral pathogen and reflects the patterns of antibody responses as the pandemic transitions into an endemic stage.
Introduction
Coronavirus disease 2019 (COVID-19) has severely impacted human health worldwide. Within the United States, New York City emerged as an early epicenter of infection and cases, with the first reported case on February 29, 20201. However, public reporting of cases over time was inconsistent due to limited testing capacity early in the pandemic1. Moreover, testing was differentially available to groups depending on their socioeconomic status2, impacting understanding of the course of the epidemic in different populations.
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, has rapidly evolved to increase its transmission potential and to evade immunity derived from prior infection and vaccination, contributing to several waves of infections3. The emergence of viral variants capable of overcoming pre-established immune responses is mainly attributed to mutations in the spike surface glycoprotein. During the first wave of infections, the wild type (WT) virus acquired the D614G mutation that increased virus infectivity and transmissibility4, 5. These strains prevailed for several months until the second wave of infections caused by the Alpha, Beta and Gamma variants occurred in different parts of the world. A 3rd wave of infections was caused by the Delta variant which showed higher transmissibility and enhanced severe disease. A significant reduction of viral neutralization by sera from convalescent or vaccinated individuals was detected6, 7. In December 2021, the emergence of the Omicron BA.1 strain, which carried over 30 spike protein mutations, resulted in a 4th wave of infections. Since then, distinct Omicron sub-lineages have circulated worldwide, including BA.2, BA.5, BQ.1, BQ.1.1, BF.7, and more recently XBB.1, XBB.1.5, BA.2.86, EG.5.1, JN.1, and derived variants.
Throughout the pandemic, monitoring the level of immunity in the population has been crucial for public health. Sero-monitoring has guided safety measures and vaccination policies worldwide to reduce viral infection, transmission, and severe disease8. Laboratories can measure binding antibodies to the spike (S) protein, which indicates past infection or vaccination, or antibodies to the nucleoprotein (NP), which indicates a past infection. The level of anti-SARS-CoV-2 antibodies and T cell responses correlate with an individual’s capacity to combat the virus and with protection against infection and disease9, 10. Importantly, binding antibody levels correlate with the number of prior exposures, and with the level of neutralizing antibodies11, 12. Hence, sero-surveillance for SARS-CoV-2 provides an estimate of the level of immunity in the population and reflects the patterns of the pandemic. Moreover, existing case reporting systems for COVID-19 are limited by changing testing availability and behaviors over time.
Here, we report the results of a prospective, hospital-based cross-sectional study of anti-S and anti-NP sero-monitoring from the beginning of February 2020 through July 2022, and follow-up measurements in August-October 2023. We describe the longitudinal seroprevalence in two different groups: a) an ‘urgent care’ group that was enriched for acute cases of COVID-19 early during the pandemic, and b) a ‘routine care’ group, resembling the general population, as previously described8. Stratification by zip code allowed us to study the geographical distribution of the seroprevalence and antibody titers in residents of NYC throughout the pandemic.
Methods
Study participants and human samples
Sample collection and testing started on February 9th, 2020 and ended on July 18th, 2022. Follow up testing started on August 28st 2023 and ended on October 2nd 2023. Collection was performed on a weekly basis (n=608/week on average) and in a blinded manner. Plasma was collected from ethylenediaminetetraacetic acid (EDTA)-treated blood specimens that remained from standard-of-care testing at the Mount Sinai Hospital (MSH) Blood Bank. Samples were sorted to collection week, visit location and practice including: OB/GYN, labor and deliveries; oncology; surgery; cardiology; emergency department; and other related hospital admissions. Up to 152 specimens at each location in a defined collection week were randomly selected.
During sample aliquoting, specimens were captured into a coded database with randomly assigned identifiers and verified patient categories. Additional specimens obtained from the same patient within the same week were not considered for analysis. Samples collected 7 days apart or more from the first specimen collected were considered as independent samples of the population. Typically, 67-314 specimens were selected weekly from the ‘routine care’ group (OB/GYN, labor and deliveries, oncology, surgery, and cardiology) and 67 to 314 samples from the ‘urgent care’ group consisting of specimens from patient visits to the emergency department and urgent care. Some individuals had a PCR test to diagnose SARS-CoV-2 infection. The analyses of samples collected during the week ending on February 9th 2020 to July 5th 2020 were reported previously8. We incorporate all samples with additional demographics data in the current analysis in order to describe seroprevalence and antibody titers since the beginning of the pandemic. We stratify our analyses by ‘urgent care’ vs ‘routine care’ groups to capture the fact that access to non-urgent medical care was limited at the beginning of the pandemic13.
Specimen metadata included information on the individuals’ age, sex, race, insurance status, and whether they had received a COVID-19 vaccination prior to specimen collection. Due to the low number of children in our study, we limit our main analysis to adults, and divided them into 18-44 years, 45-64 years, and 65+ years age groups. Information on race was based on self-reported categories (e.g., Asian and Pacific Islander, African American, White, Other, and unknown). Information on ethnicity (e.g., Hispanic or non-Hispanic) was not consistently available and, thus, not included in our analysis. Medicaid and Medicare were collapsed into a public insurance category. All data used in this study was anonymized and recoded through the use of an honest broker following local regulations.
Epidemiological waves
We categorized samples by time into 5 epidemiological waves, that correspond to successive peaks in COVID-19 incidence in NYC: Wave 1 (February 9th to August 30th, 2020; ancestral SARS-CoV-2 and D614G mutant circulating), Wave 2 (August 31st, 2020, to June 20th, 2021; Iota (B.1.529) and Alpha (B.1.1.7) circulating), Wave 3 (June 21st to October 31st, 2021; Delta (B.1.617.2) circulating), Wave 4 (November 1st, 2021, to March 6th, 2022; Delta (B.1.617.2) and Omicron BA.1 (B.1.1.529.1) circulating), and Wave 5 (March 7th to July 18th, 2022; Omicron BA.2 (B.1.1.529.2) and Omicron BA.5 (B.1.1.529.5) circulating). Epidemiological waves and circulating variants are depicted in Figures 1A and 1B respectively. The follow up period (August 28st 2023 - October 2nd 2023), was not associated to a specific wave timeframe, and was included to assess the maintenance of seroprevalence and durability of spike binding antibody titers.
Geographical maps
The neighborhood tabulation area (NTA) map created by the Department of City Planning and sourced from NYC OpenData was used as the base geometry for our data analyses. For each NTA, the average of the common logarithm of spike titer values was calculated and plotted as a shaded map using the ggplot2 and sf packages in R. For plots showing progression over time, slices were taken from each data range and plotted as separate facets.
Statistical analysis
We describe the distribution of results graphically and through proportions, stratified by ‘urgent care’ vs ‘routine care’ groups. Visualizations of antibody prevalence and titers were generated using Matplotlib (3.7.2), Seaborn (0.11.2), Pandas (1.5.3), Numpy (1.26.0), and Scipy (1.11.3) packages within python 3.11.5. Geographical charts, and live data explorer were generated using R version 4.3.2, with the packages sf (1.0.14), DT (0.31), zoo (1.8.12), Hmisc (5.1.1), shiny (1.8.0), scico (1.5.0), ggpubr (0.6.0), bsplus (0.1.4), forcats (1.0.0), ggplot2 (3.4.4), ggbreak (0.1.2), patchwork (1.1.3), lubridate (1.9.3), hrbrthemes (0.8.0), data.table (1.14.10), shinyWidgets (0.8.0), RColorBrewer (1.1.3), shinydashboard (0.7.2), and shinycssloaders (1.0.0).
Based on a priori considerations, we include age, sex, race, and insurance status in our multivariable models: age representing differential humoral response to infection and increased risk of severe disease14, 15 and therefore a marker of potentially enhanced social distancing behaviors16, sex linked to health seeking behaviors17 as well as gynecological appointments, and race and insurance status measures of poverty, health seeking behaviors, and exposure to SARS-CoV-2 through a variety of mechanisms18.
Differences in seropositivity by demographic group over time were assessed through multivariable, Poisson regression models that included robust variance estimates19. These models output prevalence ratios (PR) and 95% confidence intervals (CI), which were interpreted relative to a reference group. These models estimated the direct effect of age, sex, race, and insurance status on seropositivity. These models are initially stratified by epidemiological wave. To estimate changes in seropositivity by demographic group across wave, we also specified a larger model, where an interaction term was placed between wave and every demographic variable. The p-value from this interaction term indicates if there has been a significant shift in differences across wave, for a given demographic variable.
We conduct several sensitivity analyses. In a separate analysis, we also include vaccination, but limit the timing to data collected after January 11, 2021, when the vaccine started becoming available to certain age groups in the general population20. And in a third set of models, we exclude individuals 65+ based on their ability to access Medicare. The regression models were implanted in SAS version 9.4 (SAS Institute, Cary, NC, USA). We used an alpha level of 0.05 to assess significance.
Ethical approval
The study protocol HS 20-00308 was reviewed by the Mount Sinai Health System Institutional Review Board, Icahn School of Medicine at Mount Sinai, and it was exempt from human research as defined by regulations of the Department of Health and Human Services (45 CFR 46. 104).
Data availability
All data are available under ImmPort accession #xxx.
Results
Overall, there are 55,092 individuals sampled between February 9, 2020 and July 18, 2022: 21,075 from urgent care group and 34,017 from the routine care group. Of these 882 were children <18, and 54,210 were adults. The demographic distribution among adults is shown in Table 1. In the routine care groups, a plurality of individuals were in the 18-44 year age group (16,499, 49%), a majority female (23,155, 69%), a plurality self-identified as White (14,971, 45%), and about half had private insurance (16,386, 49%). The distribution of these variables for children is shown in Supplementary Table 1. COVID-19 vaccination began in January 2021 in the general population and varied across a number of groups, including race. Uptake was faster in Asian and Pacific Islanders and White New Yorkers compared to African Americans (e.g., at wave 3, uptake in these three groups was 63%, 55%, and 50%, respectively) (Supplementary Table 5).
We measured spike antibody prevalence and titers over five epidemiological waves based on average case counts (7-day) of SARS-CoV-2 infections in New York City, reported by the NYC Department of Health and Mental Hygiene (NYC DOHMH) (Figure 1A). Viral variant prevalence through this period continuously changed, from the initial circulating of ancestral variants to the appearance of variants of the Omicron lineage (Figure 1B). Overall, among adults in the routine care group, spike protein seropositivity increased from an average of 13% in wave 1, to 36%, 76%, 88%, and 93% in the following four waves (Figure 1, C-D; Supplementary Table 2), consistent with seroprevalence estimates reported by the NYC Department of Health and Mental Hygiene (DOHMH)21. A sharp increase in antibody prevalence was detected during waves 2 and 3, starting on February 2021. This increase was corresponded temporally with expanded availability of SARS-CoV-2 vaccinations, the emergence of the Iota, Alpha, and Delta variants, and changes in public health measures. Importantly, antibody prevalence was maintained during the follow up period from August 28st 2023 to October 2nd 2023 (Figure 1, C-D). Vaccinated individuals showed significantly higher antibody prevalence than non-vaccinated as early as January 2021, with the largest differences just prior to and during wave 3 (Figure 2, G-H) owing to incomplete vaccination in the earlier periods and unvaccinated individuals experiencing high infection rates during wave 3. By wave 5, spike protein seropositivity was uniformly high (≥90% in all categories). Similar trends were observed in populations stratified by gender and age, with no significant differences among intra-group strata (Figure 2, A-D). Differences by race were observed, for example, 14% of African Americans and 18% of those within the Other racial category were seropositive in wave 1, versus 7% of Asians and Pacific Islanders (Figure 2, E-F).
The multivariable Poisson regression model of spike protein seropositivity reveals several trends within demographic group and across time (Table 2). At wave 1, seropositivity was higher in younger age groups compared to older (1.66 times higher in those 18-44 vs 65+, 95% CI: 1.41, 1.96). By wave 5, this had largely attenuated (interaction p-value <0.0001). Larger disparities by race were also observed at the beginning of the pandemic. Compared to White New Yorkers, Black (prevalence ratios (PR): 1.24, 95% CI: 0.94, 1.63) and Other New Yorkers (PR: 1.64, 95% CI: 1.44, 1.88) had higher seropositivity, and Asians and Pacific Islanders had relatively low seropositivity (PR: 0.67, 95% CI: 0.51, 0.87). During waves 3 and 4, some of these associations reversed, with lower seropositivity in Black and Other groups compared to White New Yorkers. Around a similar or earlier time (waves 2 and 3), vaccination coverage was relatively low among Black compared to White New Yorkers (Supplementary Table 5). An analysis of spike protein seropositivity across racial ethnic groups stratified by vaccination status reveals significant differences among the unvaccinated (Supplementary Table 6), with seropositivity relatively high in waves 1 through 2 and then relatively low in waves 3 and 4 in Black New Yorkers. However, among the vaccinated, seropositivity remained equivalent across racial groups. Indeed, in the models with vaccination status that cover periods once the vaccine had been rolled out (Supplementary Table 3), this factor emerges as a significant contributor to seropositivity, with a larger impact in wave 2, where seropositivity was 1.95 times higher among those vaccinated compared to unvaccinated (95% CI: 1.87, 2.03), while it was only 1.04 times higher in wave 5 (95% CI: 1.02, 1.07) (interaction p-value <0.0001).
Anti-spike antibody titers in seropositive individuals were induced at moderate to high levels during the first waves of infections, with modest antibody decay over time (Figure 1, E-F; Supplementary Figure 1). Amid the second wave, titers increased gradually reaching the highest levels during wave 5. No significant differences were detected when titers were stratified by gender, age or race (Figure 3, A-F). As expected, the introduction of vaccination during wave 2 resulted in increased antibody prevalence and titers at the time of sample collection, while individuals without a record of vaccination displayed significantly lower antibody prevalence through waves 2, 3 and 4 (Figures 2 and 3, G-H). These differences narrowed over time, likely due to multiple infections in unvaccinated individuals. Moreover, differences in antibody prevalence and titers in vaccinated vs non-vaccinated individuals were preserved when stratification was based on sex, age, and race (Supplementary Figures 2-7).
For select weeks (July 6th-July 14th 2020, February 15th-March 1st 2021, June 1st – June 14th 2021, August 16th – August 30th 2021, May 23rd – May 30th 2022, August 28th – Oct 2nd 2023) we also measured seropositivity and titers for nucleoprotein (NP). Although initially in 2020 NP seropositivity was similar to the spike seropositivity, this changed in 2021 as NP seropositivity fell behind spike positivity by a relatively large margin (Figure 4). This was influenced not only by vaccination (which only induces anti-spike antibodies) but also by waning of NP antibodies over time and the low induction of NP antibodies during breakthrough infections. Our estimates during mid 2021 (June 1st – June 14st, 57% (95% CI = 53.1–60.9) to late 2021 (August 16th – August 30th, 54% (95% CI = 50.5–57.5)) are higher than the reported nationwide NP seroprevalence (December 2021 33.5% (95% CI = 33.1–34.0))22. These differences may be due to discrepancies in the time points of antibody measurements. Antibodies induced or boosted during the Delta wave may have waned over the subsequent months, due to the relatively fast decay rates of NP antibodies induced by infection11. Alternatively, differences in the sensitivity of the methods used or higher COVID-19 incidence in NYC compared to nationwide levels23 could have played a role.
To analyze the geographical distribution of antibody prevalence and titers through the five boroughs of NYC (Manhattan, Brooklyn, Queens, The Bronx, and Staten Island), we used neighborhood tabulation areas (NTAs) to build maps of the city divided into 216 NTAs (Figure 5). We detected higher antibody prevalence in areas located outside Manhattan during all waves of infections, that were more evident during waves 1-3, when antibody prevalence was lower (Figure 5, A-C). Importantly, seroprevalence reached high levels (above 80%) almost homogenously throughout the five boroughs by wave 5 (Figure 5, E). Likewise, higher titers were detected in boroughs other than Manhattan (Figure 6), with areas reaching the highest titers located in Queens (Figure 6, E).
Discussion
Serosurveillance systems provide means to describe population patterns in infection and vaccination over time. Overall, in a hospital-based study in New York City, we previously described a sharp increase in anti-spike antibody prevalence early during 2020 due to SARS-CoV-2 infections8. By May 2021, seroprevalence reached levels above 70% due to both infections, and vaccination programs. Here we report on the continued sero-monitoring efforts to capture antibody prevalence and titers over the course of the pandemic. We show that antibody prevalence gradually increased over time, reaching levels above 90% by January 2022. Seroprevalence and titers were maintained over time, with similar trends observed during a follow up period in August-October, 2023, likely due to constant antigenic exposures by vaccine boosters and breakthrough infections. Moreover, we analyzed the geographical distribution of antibody prevalence and titers through the five boroughs of NYC. Finally, we analyzed NP titers – indicative of immunity derived by infection – at discreet timepoints during the pandemic.
Our seroprevalence estimates for the most part overlap with estimates from the NYC DOHMH. However, there are differences during the early phase of the pandemic and later during waves 4 and 5. Specifically, the NYC DOHMH reports a large proportion (>50%) of seropositive individuals in wave 1 who were tested for SARS-CoV_2 by PCR. Individuals seeking testing were probably more likely to be infected and potentially already antibody positive given the long initial incubation time1. This is similar to what we observed in the urgent care group that was biased towards SARS-CoV-2 positives in the first wave. In more recent waves, seropositivity is lower according to the NYC DOHMH than in our serosurveillance study. This could reflect differences in the assays used or differences in the study populations.
The COVID-19 pandemic intensified racial and economic disparities in health24. Our data collected in the early waves, indicate that seropositivity was relatively high in adults self-identifying as Black and Others, pointing to higher infection rates in essential workers lacking the possibility to socially distance or due socioeconomic constraints leading to more crowded living conditions25, 26. Our serosurveillance system revealed other sociodemographic trends by age. These reflect epidemiological risk patterns and vaccination policies. For example, early in the pandemic, seropositivity was greater in younger adults, in accordance with risk messaging directed towards older adults. As vaccination became available and was prioritized for older ages27, seropositivity became relatively higher in older age groups.
The spike antibody response induced towards the emerging Omicron sub-lineages is mostly cross-reactive with ancestral strains28, hence the seroprevalence estimates presented here are likely to be preserved towards the spike of Omicron sub-variants. Although the neutralization capacity of sera from individuals exposed to ancestral SARS-CoV-2 antigens through infection or vaccination is severely reduced against Omicron strains29, non-neutralizing antibodies may protect through alternative mechanisms mediated by Fc-dependent antibody effector functions30 and the presence of antibodies likely also indicates the presence of T-cell responses which largely cross-react between different variants of concern 31. Dissecting the contribution of non-neutralizing cross-reactive antibodies upon exposure to emerging variants of the Omicron lineage will be important to have a broader picture of an individual’s immunity.
Our analysis has some limitations. Our catchment area was tied to individuals visiting the Mount Sinai Health System, and we may have limited scope in generalizing outside the NYC metropolitan area. Residual, de-identified samples are collected, which impedes following individuals longitudinally. Representation of samples by different zones was limited, which limited the geographical distribution analysis.
Conclusions
New York City was an early epicenter of the COVID-19 pandemic. Leveraging a hospital-based serosurveillance system, we detect several patterns in line with the epidemiology of SARS-CoV-2 and the roll-out of the COVID-19 vaccines. Notably, by early 2022, virtually the entire population had some evidence of COVID-19 seropositivity. Importantly our data clearly show increasing overall seropositivity in the general population as SARS-CoV-2 vaccines became available. The overall sustained high level of implied immunity in the total population also is noteworthy. Finally, we believe our data set provides a great resource for additional modelling studies and could help to predict immunity and seroprevalence scenarios during future pandemics.
Conflict of interest
The Icahn School of Medicine at Mount Sinai has filed patent applications relating to SARS-CoV-2 serological assays and NDV-based SARS-CoV-2 vaccines which list Florian Krammer as co-inventor. Dr. Simon is listed on the SARS-CoV-2 serological assays patent. Mount Sinai has spun out a company, Kantaro, to market serological tests for SARSCoV-2. Dr. Krammer has consulted for Merck, Seqirus, CureVac and Pfizer in the past, and is currently consulting for Pfizer, 3rd Rock Ventures, GSK and Avimex and he is a co-founder and scientific advisory board member of CastleVax. The Krammer laboratory has been collaborating with Pfizer on animal models for SARS-CoV-2.
Author contributions
F.K., V.S., H.v.B., A.G. and E.M.S. conceived and designed the study. H.v.B. wrote and maintained the Institutional Review Board protocol. J.T., A.R., D.B., G.S., S.M., M.F., T.Y., and L.S. performed the serological assays. J.M.C. and F.K. supervised serological assays. D.B., S.M, and B.M. collected, organized and aliquoted plasma samples. G.S. produced and purified recombinant proteins. J.M.C., J.T., B.M., and A.L.W. collected and analyzed the data. J.M.C., B.M., and D.F. generated figures. J.M.C. and A.L.W. wrote original manuscript draft. All authors edited and approved the manuscript.
Acknowledgements
This effort was supported by the Serological Sciences Network (SeroNet) in part with Federal funds from the National Cancer Institute, National Institutes of Health, under Contract No. 75N91019D00024, Task Order No. 75N91021F00001. The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does mention of trade names, commercial products or organizations imply endorsement by the U.S. Government. This work was also partially funded by the Centers of Excellence for Influenza Research and Surveillance (CEIRS, contract # HHSN272201400008C), the Centers of Excellence for Influenza Research and Response (CEIRR, contract # 75N93021C00014), by the Collaborative Influenza Vaccine Innovation Centers (CIVICs contract # 75N93019C00051) and by institutional funds.
Footnotes
↵# PARIS study group: Fatima Amanat, Guha Asthagiri Arunkumar, Christina Capuano, Jordan Ehrenhaus, Shelcie Fabre, Matthew M Hernandez, Kaijun Jiang, Brian Lerman, Meagan McMahon, Daniel Stadlbauer, Jessica Tan, Catherine Teo, Kathryn Twyman,