ABSTRACT
Background Epidemiologic risk factors for incident SARS-CoV-2 infection as determined via prospective cohort studies greatly augment and complement information from case-based surveillance and cross-sectional seroprevalence surveys.
Methods We estimated the incidence of SARS-CoV-2 infection and risk factors in a well-characterized, national prospective cohort of 6,738 U.S. adults, enrolled March-August 2020, a subset of whom (n=4,510) underwent repeat serologic testing between May 2020 and January 2021. We examined the crude associations of sociodemographic factors, epidemiologic risk factors, and county-level community transmission with the incidence of seroconversion. In multivariable Poisson models we examined the association of social distancing and a composite score of several epidemiologic risk factors with the rate of seroconversion.
Findings Among the 4,510 individuals with at least one serologic test, 323 (7.3%, 95% confidence interval [CI] 6.5%-8.1%) seroconverted by January 2021. Among 3,422 participants seronegative in May-September 2020 and tested during November 2020-January 2021, we observed 161 seroconversions over 1,646 person-years of follow-up (incidence rate of 9.8 per 100 person-years [95%CI 8.3-11.5]). In adjusted models, participants who reported always or sometimes social distancing with people they knew (IRRalways vs. never 0.43, 95%CI 0.21-1.0; IRRsometimes vs. never 0.47, 95%CI 0.22-1.2) and people they did not know (IRRalways vs. never 0.64, 95%CI 0.39-1.1; IRRsometimes vs. never 0.60, 95%CI 0.38-0.97) had lower rates of seroconversion. The rate of seroconversion increased across tertiles of the composite score of epidemiologic risk (IRRmedium vs. low 1.5, 95%CI 0.92-2.4; IRRhigh vs. low 3.0, 95%CI 2.0-4.6). Among the 161 observed seroconversions, 28% reported no symptoms of COVID-like illness (i.e., were asymptomatic), and 27% reported a positive SARS-CoV-2 diagnostic test. Ultimately, only 29% reported isolating and 19% were asked about contacts.
Interpretation Modifiable epidemiologic risk factors and poor reach of public health strategies drove SARS-CoV-2 transmission across the U.S during May 2020-January 2021.
Funding U.S. National Institutes of Allergy and Infectious Diseases (NIAID).
INTRODUCTION
A major challenge of controlling community transmission of SARS-CoV-2 is that the virus’ infectious period allows for onward spread without, or prior to, diagnosis of infection[1,2], including by fully-vaccinated individuals.[3] One national study in the United States (U.S.) prior to the vaccine era estimated that there were 5 undiagnosed infections for every diagnosed case.[2]
While SARS-CoV-2 is understood to be transmitted from person-to-person via airborne and droplet spread, to date, the incidence of SARS-CoV-2 infection and risk factors for incident infection have not been adequately characterized by routine case-based surveillance of SARS-CoV-2 diagnoses or by cross-sectional seroprevalence studies.[4–6] It is critical for prospective studies to investigate COVID-19’s evolving epidemiology and risk factors for SARS-CoV-2 acquisition in communities, including the uptake and impact of non-pharmaceutical interventions (NPIs)[7], and the reach of public health strategies aimed at controlling community transmission, including testing, quarantine, isolation, contact tracing, and vaccination.
Globally, few community-based prospective epidemiologic studies of SARS-CoV-2 incidence and risk factors have been undertaken. One recent global systematic review of observational studies of SARS-CoV-2 that employed serologic or polymerase chain reaction (PCR) testing found 18 prospective studies.[8] Most were focused on healthcare workers or other occupational groups, individuals in congregate settings, evacuees, or cruise ship patrons; none were community-based (i.e., focused on risk factors in communities vs other higher risk populations/settings).[8] Since then, there have been publications from a national community-based prospective cohort in the U.K. [9–11], but to date, there have been no prospective community-based studies with biomarkers in the U.S. Such studies are needed to help inform aspects of implementation of the public health response and policies, both to the current pandemic and future ones.
In late March 2020, we launched the prospective Communities, Households and SARS-CoV-2 Epidemiology (CHASING) COVID Cohort.[12] We describe the incidence of SARS-CoV-2 infection and risk factors for SARS-CoV-2 seroconversion during May 2020-January 2021, as well as the reach and uptake of public health strategies aimed at controlling community spread among those who seroconverted.
METHODS
Recruitment
We used internet-based strategies[13,14] to recruit a geographically and socio-demographically diverse cohort of adults into longitudinal follow-up with at-home specimen collection. To be eligible for inclusion in the cohort, individuals had to: 1) reside in the U.S. or a U.S. territory; 2) be >18 years old; 3) provide a valid email address for follow-up; and 4) demonstrate early engagement in study activities (provision of a baseline specimen or completion of >1 recruitment/enrollment visit). Details of the study design and recruitment procedures are described elsewhere.[15] The full cohort includes participants from all 50 U.S. states, the District of Columbia, Puerto Rico, and Guam (Supplemental Figure S1). Of the 6,738 participants in the full cohort, 4,510 (67%) had at least one serologic test and comprised the study population for this analysis (Table S1).
Data collection
Cohort recruitment and enrollment visits were completed between March 28-August 21, 2020, during which multiple rounds of interviews took place. Demographic and COVID-19 related risk factors were collected at baseline. From three follow-up interviews between August and November 2020, we obtained repeated measurements of COVID-19 symptoms, laboratory testing (PCR or serologic), hospitalizations, use of NPIs such as mask use and social distancing, public health strategies such as quarantine, isolation, and contact tracing encounters, and other time-varying factors.
During May-September 2020 (Period 1) and November 2020-January 2021 (Period 2), participants were invited to complete serologic testing using an at-home self-collected dried blood spot (DBS) specimen collection kit. DBS cards were sent from and returned to the study laboratory (Molecular Testing Laboratories [MTL], Vancouver, Washington) via the U.S. Postal Service using a self-addressed, stamped envelope containing a biohazard bag™.
All DBS specimens were tested by the study laboratory for total antibodies (Total Ab) using the Bio-Rad Platelia test for IgA, IgM, and IgG which targets the SARS-CoV-2 nucleocapsid protein (manufacturer sensitivity 98.0%, specificity 99.3%).[16] Other studies have independently validated this assay and found average sensitivity and specificity of 91.7% and 98.8%, respectively.[17–19] This assay was also validated for use with DBS by the study laboratory, which found 100% sensitivity and 100% specificity (MTL, personal communication).
Outcomes
Cumulative incidence of SARS-CoV-2 infection
Among participants who underwent serologic testing, we estimated the serology-based cumulative incidence of SARS-CoV-2 as the proportion of individuals with a positive Total Ab test in either of the two time periods (i.e., number of participants ever positive) divided by the total number of persons with one or more Total Ab tests (i.e., number of participants ever tested). We adjusted cumulative incidence estimates for laboratory test error, assuming a sensitivity of 91.7% and a specificity of 98.8%[17–19] using the following formula[20]:
Observed SARS-CoV-2 seroconversion
Among those individuals with two Total Ab tests, an observed seroconversion was defined as a negative Total Ab test in Period 1 followed by a positive Total Ab test in Period 2. We estimated person-years of follow-up using the collection dates for each specimen in Periods 1 and 2. When the collection date was missing, we used the date the laboratory received the sample. The seroconversion date was assigned as the midpoint between the first and second specimen collection dates for person-time calculations.
Exposures
Individual-level COVID-19 risk factors
We examined epidemiologic risk factors for SARS-CoV-2 individually and as part of a composite score. We considered the following risk behaviors and factors reported by participants prior to specimen collection: Household factors (household crowding [4 or more people living in a household], having a child in the household, and having a confirmed COVID-19 case in a household member prior to participant testing positive), spending time in public places (attending mass gatherings such as protests, spending time in an indoor restaurant or bar, spending time at an outdoor restaurant or bar, visiting places of worship, or visiting public parks or pools), mask use indoors (for grocery shopping, visiting non-household members, at work, and in salons or gyms), mask use outdoors, gathering in groups with more than 10 people (indoors only, outdoors only, and indoors and outdoors), travel during the pandemic (recent air travel and public transit use), and individual-level factors that may increase the risk of severe COVID-19 (substance use, binge drinking, and medically diagnosed comorbidities).
Global social distancing assessment
While social distancing in specific scenarios is addressed in some of the above individual risk factors (such as spending time in various public places), we were specifically interested in the association between social distancing in general and incident SARS-CoV-2 infection. We asked two global questions on social distancing, which were “In the past month, how often have you practiced social distancing with: a) people you know and b) people you do not know,” with possible response options of Always, Sometimes, or Never. These assessments were not included in the calculation of the composite risk score.
Composite score of COVID-19 risk factors
We computed a composite COVID-19 risk score as many of the above COVID-19 risk factors may be highly correlated. We applied Least Absolute Shrinkage Selection Operator (LASSO) regression to select the set of risk factors which best predicted seroconversion.[21] Scores were assigned to each participant based on whether they engaged in the risk factors selected by the LASSO regression model and were normalized between 0 and 100. High scores indicate engagement in high-risk activities. Details are in the Statistical Appendix. The composite score was divided into tertiles (low, medium, high) for analysis.
Statistical analysis
Cumulative incidence estimates were stratified by baseline characteristics and epidemiologic risk factors. Crude and adjusted Incidence Rate Ratios (IRRs) of seroincidence and associated 95% confidence intervals (CIs) were estimated using Poisson regression. We examined crude seroconversion rates by sociodemographic characteristics and each risk factor. Finally, we separately modeled three exposure variables: 1) social distancing with “people you know” (Yes/No); 2) social distancing with “people you don’t know” (Yes/No); and 3) the composite COVID-19 risk score(high/medium/low). Two multivariable models were constructed for each exposure variable, including models that adjusted for age, gender, race/ethnicity and comorbidities (Model 1); and a model that further controlled for changes in community-level COVID-19 transmission (Model 2). See Statistical Appendix for details. All data were cleaned and analyzed in R (version 4.0.3) and SAS (V9.4).
Ethical Approval
The study protocol was approved by the Institutional Review Board at the City University of New York (CUNY).
RESULTS
Sample characteristics
The characteristics of the study sample are shown in Table 1. A total of 4,232 persons underwent serologic testing in Period 1, and 3,883 were tested in Period 2 (Table S1). Of the 4,510 participants who tested at least once, 3,605 (80%) tested at both time points (Table 1). Differences between those participants testing in Period 1 and Period 2 were negligible (Table S1). The median time between specimen collection dates for participants providing specimens for both serologic tests was 190 days (IQR 152-201) (Figure S2).
Cumulative incidence of SARS-CoV-2 as of January 31, 2021
We observed 323 seropositives among the 4,510 participants who tested at least once during follow-up, for an overall crude and adjusted serology-based estimates of cumulative incidence of 7.3% (95% confidence interval [CI]: 6.5%-8.1%) and 6.7% (95% CI 5.9%-7.6%), respectively (Table S2).
SARS-CoV-2 seroincidence, May 2020-January 2021
There were 3,422 participants who were seronegative in Period 1 with a subsequent serologic test in Period 2 who were followed prospectively for the outcome of SARS-CoV-2 seroconversion. There were 161 observed seroconversions over 1,646 person years of follow-up, for an overall incidence rate of 9.8 per 100 person-years (95% CI 8.3-11.5) (Table 2). The rate of incident SARS-CoV-2 infection was lower for females compared to males (IRR=0.69, 95% CI:0.50, 0.94), and higher for Hispanics (IRR=2.09, 95% CI 1.41-3.05) and non-Hispanic Blacks (IRR=1.69, 95% CI 0.96-2.82) compared with non-Hispanic Whites. Essential workers had higher incidence than non-essential workers (IRR=1.65, 95% CI 1.10-2.26). Incidence rates were higher among those in rural versus urban areas (IRR=1.29, 95% CI 0.89-1.29), and among those in the Midwest (IRR=1.59, 95% CI 0.98-2.56), the South (IRR=1.67, 95% CI 1.08-2.59), and the West (IRR=1.32, 95% CI 0.83-2.11) compared to the Northeast.
Table 3 shows the seroincidence and crude incidence rate ratios by epidemiologic risk factors that were present prior to or between serologic tests. There was higher incidence among those who dined indoors at restaurants or bars (IRR=1.93, 95% CI 1.39-2.70); those who visited a place of worship (IRR=1.92, 95% CI 1.26-2.84); those who wore a mask only sometimes while grocery shopping (IRR=10.57, 95% CI 4.00-30.51); those who visited indoors with people not in their own household while sometimes wearing a mask (IRR=1.94; 95% CI 1.37-3.31) or while never wearing a mask (IRR=2.62; 95% CI 1.50-4.70); those working indoors at a place of employment while never wearing a mask (IRR=2.50, 95% CI 0.98-5.26); those who wore masks only sometimes while attending a salon or gym (IRR=3.23, 95% CI 1.90-5.23); and those who reported traveling by air during August-November 2020 (IRR=1.52, 95% CI 1.05-2.17).
Poisson models of SARS-CoV-2 seroconversion for three different exposures of interest, May 2020-January 2021
Global social distancing assessment
In crude analyses, participants who reported that they always or sometimes engaged in social distancing with people they know had a statistically significantly lower seroincidence (IRRalways vs never=0.30, 95% CI 0.15-0.72; IRRsometimes vs never=0.35, 95% CI 0.16-0.86) compared with those who said they never social distanced with people they know (Table 4). In multivariable analyses that controlled for sociodemographics and comorbidities (Model 1, aIRRalways vs never=0.37, 95% CI: 0.18-0.89), and additionally for community-level transmission (Model 2, aIRRalways vs never=0.43, 95% CI: 0.21-1.04), participants who reported always social distancing with those they know (versus never) had a significantly lower seroincidence, although the 95% confidence intervals were wider. Participants reporting that they socially distanced sometimes with people they knew also had lower incidence of seroconversion compared with those never socially distancing in both models (Table 4).
Participants who reported social distancing always or sometimes (vs. never) with people they don’t know had lower seroincidence rates in both crude analyses and adjusted models, however the association was weaker than it was for social distancing with people they did know, with IRRs and aIRRs ranging from 0.53-0.64.
Clinical and public health outcomes among persons with SARS-CoV-2 seroconversion during May 2020-January 2021
Among the 161 individuals who seroconverted during May 2020-January 2021, only 27 (16.8%) were aware that they had a prior SARS-CoV-2 infection (Table 5). A substantial proportion (28.0%) recalled no symptoms of COVID-like illness (i.e., were asymptomatic cases). In terms of public health outcomes, 60.3% said that they were ever tested for SARS-CoV-2 outside the study (Table 5), but only half of them (26.7% of total) reported ever having a positive SARS-CoV-2 test. Only 29.2% said that they had ever isolated themselves from people outside their household because of their infection, and, among those who did not live alone, even fewer (17.4% overall) said they ever isolated themselves from others within their household. In terms of contact tracing, only 19.3% of all seroconverters were asked about contacts following diagnosis and only 11.8% had been informed by a contact tracer that they may have had contact with someone confirmed to have SARS-CoV-2. Only 5.0% of all seroconverters were told by a contact tracer to stay home for a period of time because they had COVID-19. Findings were similar among all 323 seropositive participants (data not shown).
DISCUSSION
We report findings from a large community-based prospective epidemiologic study of SARS-CoV-2 incidence and risk factors in the U.S. during May 2020-January 2021. Using serologic tests, we longitudinally characterized the incidence of SARS-CoV-2 infection in relation to a range of modifiable risk factors. We found that social distancing with those whom participants knew as well as those whom they did not know was highly protective against SARS-CoV-2 infection, even after controlling for other risk factors and measured confounders. Finally, public health strategies such as quarantining, testing, isolation, and contact tracing appear to have had inadequate coverage and adoption during the infectious periods of those infected with SARS-CoV-2, limiting their effectiveness at reducing SARS-CoV-2 transmission in the community. Taken together, our study findings document some of the principal reasons why the U.S. has continued to experience sustained community transmission, hospitalizations, and deaths from COVID-19 in many areas, including into the vaccine and Delta variant eras.
The protective association observed for two global measures of social distancing persisted even after controlling for other differences. This protective association was strongest for participants reporting that they always or sometimes socially distanced from people whom they knew, but was also present among those who always or sometimes socially distanced from those they did not know. Our findings suggest a need for more effective and consistent messaging around social distancing.
We observed substantially increased risk, in dose response fashion, from a number of other key epidemiologic risk factors and exposures reflected in a composite risk score. Among participants in the top tertile of the risk score the risk of seroincidence was 3-fold higher, accounting for 55% of the observed seroconversions. Reducing multiple risk factors (e.g., through policies on masking, mass gatherings, indoor dining/bars, social distancing, air travel) would likely substantially reduce community transmission.
Our findings suggest that elevated risk among essential workers, observed early in the U.S. pandemic, persisted into the second phase of the pandemic. Essential workers risk exposure to SARS-CoV-2 not only in their workplaces, but also in their communities and as part of their commutes to and from work, especially when using public transportation. The increased burden of risk of SARS-CoV-2 infection in essential workers is shared with their household members, among whom transmission is very efficient.[28] National, state and local workplace safety measures, mandates, and policies that protect essential workers have the added benefit of protecting household members and other close contacts of workers.
Finally, detailed examination of 161 individuals with SARS-CoV-2 seroconversion showed major gaps in the reach of public health interventions aimed at reducing the spread of SARS-CoV-2. Most who seroconverted (73.3%) did not report a prior positive PCR test, and a substantial proportion were asymptomatic. Moreover, few people who seroconverted in our study reported being reached by contact tracers (11.8%). These results highlight the barriers to successful implementation of isolation, contact tracing, and quarantine. Now that rapid home tests are easily available, frequent proactive testing at home can be a more effective way to capture asymptomatic and presymptomatic cases early and prevent onward transmission.
Our study, as well as data from routine case surveillance, highlight that the drivers of racial/ethnic disparities in SARS-CoV-2 risk have not been addressed in the U.S. The ultimate drivers of these disparities need to be targeted by governments, health departments and researchers, and used to course-correct the public health response.[29] Structural factors, such as household crowding, the need to go to work to avoid income loss, and inequitable access to SARS-CoV-2 testing[30], create and perpetuate a disparate burden of SARS-CoV-2 exposure and incidence.[31] To date, no targeted strategies or policies have been deployed that aim to protect those who cannot afford missing work, including, but not only, essential workers. Public health leaders and policy makers should anticipate and proactively design pandemic response implementation strategies, with performance metrics related to inequalities[32,33], that account for and counteract the prevailing structural forces, including structural racism, that create, maintain, or exacerbate inequities in safety, health, and well-being during and after a public health crisis.[32–35]
Our study has limitations worth noting. The observed cumulative incidence in our cohort may be lower than the true cumulative incidence in our cohort because of waning of SARS-CoV-2 antibodies.[36] Recent studies suggest waning of antibodies to both nucleocapsid and spike proteins.[5] Combined with the timing of specimen collection relative to infection for many participants in our cohort (median of 190 days)[12], this could mean that we have underestimated the true cumulative incidence. Next, estimated associations between SARS-CoV-2 risk factors and incidence are subject to confounding. The crude associations we presented may also vary by setting, with interpretation for some associations further hampered by small sample sizes. Some risk behaviors may have been underreported, due to social desirability, which would bias observed associations toward the null. Finally, our study period for the current analysis pre-dated the vaccine era and the emergence of the highly transmissible and possibly more virulent SARS-CoV-2 Delta variant. We could not, therefore, examine risk factors for infection among vaccinated persons. However, given recent data showing decreasing vaccine effectiveness against the Delta variant for the outcome of infection[37,38], as well as major outbreaks and high viral loads among fully vaccinated persons[3,39], it is likely that many of our findings related to transmission risk factors also apply to vaccinated persons, as they remain at risk for breakthrough infection and onward transmission of SARS-CoV-2 when engaging in some of the same risk factors. A similar study performed in the Delta variant era might result in stronger risk factor and weaker protective factor associations with seroconversion, as well as a different estimate of proportion asymptomatic.
Conclusion
Modifiable risk factors and poor reach of public health strategies continue to drive transmission of SARS-CoV-2 across the U.S. While continuing to increase vaccine access and coverage, it remains critical for public health agencies to simultaneously reduce risk factors and address structural factors that contribute to high incidence and persistent inequities. Future research will include monitoring SARS-CoV-2 outcomes, including infections in the vaccine and Delta variant eras, at least through December 2021.
Data Availability
The data that support the findings of this study are available on request from the corresponding author, [DN]. The data are not yet publicly available, but we are preparing to post a deidentified, HIPAA compliant, public use version of our baseline and follow-up data on GitHub.
CONTRIBUTORS
DN, MSR, SGK, MMR, conceptualized the study. MSR, MC and DN performed statistical analyses. DN and MSR wrote the first draft of the paper. DN, MSR, MC, SGK, and AP contributed to interpreting the data, DN, RZ, MSR, MC, SGK, WY, AB, CM, SK, AM, MMR, DAW, AP, LW, and CG contributed to the writing and revising of the manuscript. SGK, WY, AB, CM, SK, and DN contributed to data collection, cleaning and management. DN, SGK, MMR and CG contributed to obtaining funding for the research.
FUNDING
Funding for this project is provided by The National Institute of Allergy and Infectious Diseases (NIAID), award number 3UH3AI133675-04S1 (MPIs: D Nash and C Grov), the CUNY Institute for Implementation Science in Population Health (cunyisph.org) and the COVID-19 Grant Program of the CUNY Graduate School of Public Health and Health Policy, and National Institute of Child Health and Human Development grant P2C HD050924 (Carolina Population Center).
STATISTICAL APPENDIX
LASSO
The study sample of 3,422 participants with 2 serosamples was split randomly into equally sized training and test data sets. A grouped LASSO regression was fit using training data and 10-fold cross validation was used to obtain the minimum value of lambda (the tuning parameter). The variables included in the LASSO model were household crowding, having children in the household, confirmed COVID-19 case in a household member, attending mass gatherings, indoor dining in a restaurant/bar, outdoor dining in a restaurant/bar, visiting places of worship, visiting public park or pool, gathering with groups of 10 or more indoors and/or outdoors, mask use while grocery shopping, mask use while visiting non-household members, mask use indoors in gym/salons, mask use indoors at work, mask use outdoors, using public transit, recent air travel, alcohol use,and substance use. Using the minimum lambda value, a grouped LASSO model was run on the entire dataset (training + test) to obtain factors that best predicted seroconversion. The LASSO model selected having a confirmed COVID-19 case in a household member, indoor dining in a bar/restaurant, gathering with groups of 10 or more outdoors, and mask use indoors in salons/gyms as the most predictive of seroconversion in our cohort. A logistic regression model was fitted using these variables selected by LASSO regression with seroconversion (Y/N) as the outcome. Coefficients of the variables were multiplied by 10 and rounded up to the nearest integer to create a score associated with that variable/risk factor. If a participant engaged in a given risk factor, they were assigned the score associated with that risk factor. Scores were then summed across risk factors, normalized between 0 and 100, and a total composite risk score was created for each participant.
County-level COVID-19 community transmission
We used population-based, county-level death rates, lagged by 23 days [22,23], as a proxy for community transmission. We assumed that data on the number of deaths for a given day represented community transmission that was occurring 23 days earlier, specifically 5 days from infection to symptom onset (reflecting the average incubation period) [24]; 5 days from symptom onset to severe disease/hospitalization [25]; and 13 days from hospitalization to death. [26]
We obtained the cumulative number of COVID-19 deaths per 100,000 persons for each county, using data from the New York Times Github website (from 01/21/2020 to 01/07/2021)[27] on the days samples were collected in Period 1 and Period 2 for each participant. We calculated the rate of change in cumulative COVID-19 deaths per 100,000 at three time points for each participant, including between: 1) January 2020 (approximate start of U.S. pandemic) and the date of specimen collection in Period 1; 2) the start of the pandemic and the date of specimen collection in Period 2; and 3) the difference between the specimen collection dates in Period 1 and Period 2. The rate of change in COVID-19 deaths per 100,000 were divided into tertiles of low, medium and high community transmission.
Supplemental Material
ACKNOWLEDGEMENTS
The authors wish to thank the participants of the CHASING COVID Cohort Study. We are grateful to you for your contributions to the advancement of science around the SARS-CoV-2 pandemic. We thank Prof. Patrick Sullivan and MTL for local validation work on the serologic assays for use with DBS that greatly benefited our study. We are also grateful to MTL Labs for processing specimen collection kits and serologic testing of our cohort’s specimens.
Footnotes
CONFLICTS OF INTEREST: None declared