Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

COVID-19 cumulative incidence, asymptomatic infections, and fatality in Long Island, NY, January–August 2020: A cohort of World Trade Center responders

  • Olga Morozova ,

    Roles Data curation, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

    olga.morozova@stonybrookmedicine.edu

    Affiliations Department of Family, Population and Preventive Medicine, Renaissance School of Medicine, Stony Brook University (SUNY), Stony Brook, New York, United States of America, Program in Public Health, Stony Brook University (SUNY), Stony Brook, New York, United States of America

  • Sean A. P. Clouston,

    Roles Conceptualization, Formal analysis, Methodology, Visualization, Writing – original draft, Writing – review & editing

    Affiliations Department of Family, Population and Preventive Medicine, Renaissance School of Medicine, Stony Brook University (SUNY), Stony Brook, New York, United States of America, Program in Public Health, Stony Brook University (SUNY), Stony Brook, New York, United States of America

  • Jennifer Valentine,

    Roles Data curation, Investigation, Validation, Writing – review & editing

    Affiliation World Trade Center Health and Wellness Program, Stony Brook University (SUNY), Stony Brook, New York, United States of America

  • Alexander Newman,

    Roles Data curation, Investigation, Validation, Writing – review & editing

    Affiliation World Trade Center Health and Wellness Program, Stony Brook University (SUNY), Stony Brook, New York, United States of America

  • Melissa Carr,

    Roles Data curation, Investigation, Project administration, Supervision, Writing – review & editing

    Affiliation World Trade Center Health and Wellness Program, Stony Brook University (SUNY), Stony Brook, New York, United States of America

  • Benjamin J. Luft

    Roles Conceptualization, Funding acquisition, Investigation, Resources, Supervision, Writing – review & editing

    Affiliations World Trade Center Health and Wellness Program, Stony Brook University (SUNY), Stony Brook, New York, United States of America, Department of Medicine, Renaissance School of Medicine, Stony Brook University (SUNY), Stony Brook, New York, United States of America

Abstract

Background

New York City and Long Island, NY were early foci of the COVID-19 epidemic in the US. The effects of COVID-19 on different sub-populations, and its key epidemiologic parameters remain unknown or highly uncertain. We investigated the epidemiology of COVID-19 from January to August of 2020 in an established academic monitoring cohort of N = 9,697 middle-aged World Trade Center responders residing in Long Island, NY.

Methods

A seroprevalence survey and a series of cross-sectional surveys were nested in a prospective cohort study. Measures included IgG antibody testing, SARS-CoV-2 polymerase chain reaction (PCR) testing, review of electronic medical records, and surveys of symptoms. Correlates of infection were analyzed with multivariable logistic regression.

Results

The cohort was predominantly men in their mid-fifties; 6,597 cohort members were successfully contacted (68%); 1,042 (11%) individuals participated in the seroprevalence survey; and 369 individuals (5.6% of 6,597 study participants) underwent PCR testing. The estimated standardized cumulative incidence was 21.9% (95%CI: 20.1–23.9%), the asymptomatic proportion was 16.4% (36/219; 95%CI: 11.8–22.0%), the case hospitalization ratio was 9.4% (36/385; 95%CI: 6.6–12.7%), the case fatality ratio was 1.8% (7/385; 95%CI: 0.7–3.7%), and the hospitalization fatality ratio was 8.3% (3/36; 95%CI: 1.8–22.5%). Confirmed SARS-CoV-2 infection was associated with younger age, race/ethnicity, and being currently employed.

Conclusions

The results of the present study suggest a high cumulative incidence of SARS-CoV-2 among WTC responders in the spring and summer of 2020 and contribute to narrowing the plausible range of the proportion of infections that exhibit no symptoms. An increased risk of infection among younger employed individuals is likely to reflect a higher probability of exposure to the virus, and the racial disparities in the infection risk warrant further investigation.

Introduction

Since late 2019, SARS-CoV-2, the causative agent of the respiratory disease called COVID-19, has caused a global pandemic, with New York City (NYC) being an early focus of the outbreak in the US [1]. Long Island, NY, where essential workers for NYC often reside, was also affected early in the epidemic, with large numbers of cases and deaths attributed to COVID-19 [2, 3]. The first confirmed cases of COVID-19 were detected on March 1st, 2020 in NYC and on March 5th, 2020 in Long Island. In response to an outbreak, starting on March 12th, the State of New York began mandatory work-from-home and social distancing orders to combat the epidemic. Several studies have estimated the seroprevalence in NYC or an average for New York State [36]. To our knowledge, only one study, through convenience sampling, has provided regional estimates that included Long Island [3].

A rapidly growing body of research is investigating COVID-19 symptoms, clinical features, and infectivity [713]. Current epidemiologic work examining and modeling infection dynamics in the US population [14] relies on estimates of basic epidemiologic parameters largely derived from reported cases (e.g. [15, 16]). Population-level seroprevalence studies conducted early in the epidemic [4, 1720] have suggested low infection ascertainment rates, although some of these studies might have had selection bias. Several key epidemiologic parameters—particularly the proportions of cases that are asymptomatic and cases that were severe and required hospitalization—remain highly uncertain. A recent review has found estimates of the asymptomatic fraction ranging between 6% and 96% [21]. Given the high variability in the available estimates, concerns have been raised regarding the methods in many of these studies [22]. How the asymptomatic proportion varies depending on demographic characteristics and comorbid conditions remains unclear. Several studies have estimated the proportion of cases categorized as severe [15, 23], mostly among confirmed cases in China, and have provided results that may not be generalizable to other countries [24, 25]. Although studies have been conducted among residents of congregate settings [2628], healthcare workers [29], and patients with specific conditions [30, 31], the effects of COVID-19 on non-institutionalized and non-hospitalized populations remain largely unknown.

To fill these knowledge gaps, we used a unique opportunity to investigate COVID-19 epidemiology in an established academic monitoring cohort of World Trade Center (WTC) responders residing in Long Island, NY [32]. The primary goals of the present study were 1) to estimate the cumulative incidence of SARS-CoV-2 between the epidemic onset (mid-February) and August 2020, and 2) to estimate the asymptomatic proportion, case hospitalization ratio (CHR), case fatality ratio (CFR), and hospitalization fatality ratio (HFR) among the cohort of WTC responders. Our secondary goals included: 1) describing COVID-19 antibody response in a nested seroprevalence survey; 2) identifying correlates of infection, asymptomatic infection, hospitalization, and death in this cohort; and 3) describing the epidemic curve among symptomatic patients.

Materials and methods

The World Trade Center Health and Wellness Program (WTC-HWP) has been described previously [33]. In brief, the WTC-HWP monitors men and women who worked or volunteered at the WTC sites during and after the tragic events of 9/11/2001. Since July 2002, WTC responders have been eligible for prospective monitoring for WTC-related conditions. Soon after the state-wide lockdown was implemented on March 12th, 2020, all in-person visits were cancelled, and the WTC-HWP clinic began conducting telehealth visits. With this directive, research staff made efforts to shift research activities online, thus providing an opportunity to survey cohort members regarding their COVID-19 symptoms, testing, and status.

Sample: Seroprevalence survey

Living WTC responders participating in an ongoing monitoring program (N = 9,697) were eligible for this study. To determine cumulative incidence and the fraction of infections that did not exhibit symptoms, we conducted a seroprevalence survey by testing study participants for antibodies to SARS-CoV-2 by using the Abbott Architect IgG antibody assay (sensitivity = 93.4% 21 days after symptom onset and specificity = 100%) [34]. Testing was conducted among a consecutive sample of WTC-HWP cohort members who presented for routine annual monitoring visits and agreed to be tested. The testing dates ranged between March 25th and September 18th, 2020, and most tests were performed in June and July. The target sample size was N = 1,000; however, some patients received antibody testing at partner laboratories, thus resulting in a final sample size of N = 1,042 patients with available antibody testing results. We provide a comparison of demographic characteristics among the seroprevalence survey sample and the remaining N = 8,655 WTC-HWP cohort members.

Sample: Outreach surveys

To identify WTC-HWP cohort members affected by COVID-19, the entire cohort was asked, through multiple modalities including telehealth visits, in-person visits, via email and text messaging, and through direct calls, to provide information on their COVID-19 status at multiple time points between March 15th and August 31st, 2020. First, between March 17th and May 28th, a sample of WTC responders was seen during a comprehensive telehealth visit, at which time they were able to discuss COVID-19 symptoms (N = 1,053). Second, WTC responders were asked to respond to each of six waves of a brief needs assessment survey (N = 4,665). The survey was administered via e-mail and text messaging between March 16th and June 2nd. Third, a brief text-based survey was sent on July 9th to all WTC responders, who were asked whether they had been infected and diagnosed with COVID-19 (N = 3,417). Fourth, information was available for an additional 53 patients who neither responded to the outreach efforts listed above nor participated in the seroprevalence survey. Information about these patients’ COVID-19 history was primarily obtained during monitoring visits, during the calls made to schedule monitoring visits, and via the electronic medical records (EMR) notification system. Some patients responded multiple times to different outreach efforts, thus providing a final outreach analytic sample of N = 6,093 WTC responders (63% of eligible cohort members, hereafter the “main outreach sample”), of whom N = 770 cohort members also participated in the seroprevalence survey. The S1 Appendix provides details regarding the sample composition and overlap among the outreach modalities. To determine when symptoms began and whether individuals had been diagnosed or tested and when, we scheduled a standard follow-up call for all individuals who reported any symptoms [35] or diagnoses, or tested positive for antibodies. The follow-up calls were performed by trained research staff, certified medical assistants, trained case managers, or healthcare providers.

Noting that prior efforts had missed responders who ignored passive outreach efforts, potentially because they lacked symptoms, we conducted a brief telephone-based assessment in a sample of a sub-cohort of people who had not completed any outreach survey form. Among the 3,604 cohort members who had not responded to the main outreach surveys (“non-participants”), we selected a random sample of size N = 320, and between July 16th and September 18th, 2020, we successfully contacted N = 255 of them, thus resulting in an 80% response rate (hereafter the “confirmatory outreach sample”). No statistically significant differences were observed between these responders and the remaining N = 65 individuals in terms of age, sex, race, and employment status. Among the 255 patients in the confirmatory outreach sample, 23 participated in the seroprevalence survey.

Measures

The results of the seroprevalence survey were stratified by demographic characteristics, the history of symptoms, and the history and results of polymerase chain reaction (PCR) testing. In addition, we examined the time between the antibody testing and the onset of symptoms or PCR diagnosis, and analyzed its relationship to the antibody testing results and titers among those who tested positive. Antibody titers were available for a subset of patients who tested positive. We compared these results among patients with and without history of symptoms.

On the basis of multiple data sources, we classified all individuals as one of the following: confirmed COVID-19 positive, confirmed COVID-19 negative, symptomatic / PCR-negative / not tested for antibodies, symptomatic untested, and without known history of symptoms or testing. A case was defined as confirmed COVID-19 positive if a positive IgG antibody test or an EMR- or self-reported positive PCR test was present. A case was defined as confirmed COVID-19 negative if it was antibody-negative and was not classified as confirmed positive. Symptomatic patients with negative PCR test results and no antibody testing were classified as a separate category, because the PCR tests might have been false negative or performed outside the detection window. We show the epidemiologic curve of symptom onset dates stratified by confirmed or unconfirmed status and excluding symptomatic cases confirmed to be negative.

The COVID-19 positive cases were classified as asymptomatic, mildly symptomatic (with recovery without hospitalization), or severe (with hospitalization and/or death), on the basis of self-reported symptomatology and the EMR. Disease progression characteristics included CHR, CFR, and HFR.

Estimation of epidemiologic parameters

We provide estimates of the following epidemiologic parameters of interest: cumulative incidence of SARS-CoV-2 through August 31, 2020, asymptomatic proportion, CHR, CFR, and HFR.

Cumulative incidence was estimated as the proportion of seroprevalence survey participants who tested positive for IgG antibodies. Given the differences in some demographic characteristics between the seroprevalence survey sample and the remaining cohort members, we used logistic regression to perform standardization of cumulative incidence according to the demographic characteristics that significantly differed between groups [36]. In addition, we computed the cumulative incidence by including survey participants who tested negative for IgG antibodies, but who had a record of a positive PCR test.

The asymptomatic proportion was estimated among the seroprevalence survey participants who tested positive for IgG antibodies. To corroborate this result, we additionally estimated the asymptomatic proportion among the confirmed positive COVID-19 cases identified via the outreach survey, adjusting for potential self-selection bias based on the stratified results estimated in the main and the confirmatory outreach samples (estimation approach details are provided in the S1 Appendix).

Beyond the COVID-19 cases identified via the seroprevalence survey, the outreach survey efforts identified another 166 confirmed positive COVID-19 cases, thus increasing the sample size used to estimate CHR, CFR, and HFR, as well as to identify correlates of confirmed infection. CHR was estimated as a proportion of all confirmed positive COVID-19 cases aggregated from seroprevalence and outreach surveys with hospitalization for COVID-19 for at least 1 day. CFR was estimated as a proportion of all confirmed positive COVID-19 cases resulting in death, and with COVID-19 listed as either a primary or a contributing cause of death. HFR was estimated in a similar manner among patients hospitalized for COVID-19.

Correlates of infection

We examined the following characteristics as potential correlates of confirmed infection: Demographics including age in years, sex, race/ethnicity, and current employment status; medical factors including history of hypertension, diabetes, all-cause cancer, obstructive airway disease, upper respiratory disease, morbid obesity, depression, post-traumatic stress disorder (PTSD), and anxiety; and pharmaceutical factors including prescriptions of ibuprofen, angiotensin-converting enzyme (ACE) inhibitors, statins, and steroids. Correlates of confirmed infection were examined in the subsample of confirmed positive and confirmed negative individuals pooled from both seroprevalence and outreach surveys. Low counts of asymptomatic infections, hospitalizations, and deaths among the confirmed cases prohibited exploration of correlates of these outcomes.

Statistical analyses

Descriptive analyses used percentages and means alongside standard deviations. Sample means were compared with independent sample t-test, and proportions were compared with Pearson’s chi-squared test (p-values calculated using bootstrap with 1,000,000 iterations [37]). The variance and 95% confidence intervals (95% CI) for proportions were calculated with finite population correction if a given sample exceeded 5% of the target population [38]. The variance for the estimates of selection bias-corrected proportions was estimated conservatively by assuming independence between the main and confirmatory outreach samples. Correlates of confirmed infection were analyzed with multivariable logistic regression. Analysis was performed in the R statistical computing environment [39].

Ethics

The Stony Brook University Ethics Review Board approved this study (#604113). The participants provided signed written informed consent.

Results

Demographics

The WTC-HWP cohort (N = 9,697) comprised mostly men in their mid-fifties. Table 1 shows the demographic characteristics of the entire cohort, as well as a comparison between the seroprevalence survey sample (N = 1,042) and the remaining cohort (N = 8,655). WTC responders in the seroprevalence survey sample were slightly older and more likely to be retired than the rest of the cohort.

thumbnail
Table 1. Demographic characteristics of the World Trade Center responder cohort (N = 9,697) disaggregated by participation in the seroprevalence survey (N = 1,042 participants).

https://doi.org/10.1371/journal.pone.0254713.t001

S1 Table in the S1 Appendix shows the demographic characteristics of the main and confirmatory outreach samples. In brief, the main outreach sample was slightly older and more likely to be retired or unemployed.

Seroprevalence survey results

Table 2 shows the results of the seroprevalence survey, stratified by demographic characteristics, history of COVID-19 symptoms, and PCR testing. Overall, 21% (219/1042; 95%CI: 18.6–23.6%) of the seroprevalence survey participants tested positive for IgG antibodies. The proportion testing positive was significantly higher among younger working individuals and among black and Hispanic WTC responders than non-Hispanic white responders. Among the 760 survey participants who did not have COVID-19 symptoms, 4.7% (36/760; 95%CI: 3.3–6.5%) had a positive antibody test. A total of 166 patients, 163 of whom were symptomatic, received both PCR and antibody tests, and most results were concordant. Six of 130 (4.6%) participants who had a positive PCR test had a negative antibody test, in agreement with the antibody test sensitivity (the time between the PCR and antibody tests varied between 23 and 81 days). All of them were symptomatic and recovered; one was hospitalized. Fourteen of 36 (38.9%) patients who tested negative with PCR had a positive antibody test, all of whom were mildly symptomatic. Three asymptomatic cases had concordant positive PCR and antibody testing results.

thumbnail
Table 2. Proportion of the seroprevalence survey participants (N = 1,042) testing positive for IgG antibodies to SARS-CoV-2.

https://doi.org/10.1371/journal.pone.0254713.t002

Among 219 participants with a positive antibody test (IgG antibody test result > 1.0 AU/ml), titers information was available for 144 individuals: 119 symptomatic and 25 asymptomatic. The mean IgG antibody titers corresponded to 83 AU/ml (SD = 49) among symptomatic cases and 65 AU/ml (SD = 53) among asymptomatic cases (p = 0.099). In the available sample, we did not find any association between antibody positivity or titers and the time to antibody testing after symptom onset or PCR diagnosis.

Estimates of epidemiologic parameters

Cumulative incidence.

The proportion of seroprevalence survey participants who tested positive for IgG antibodies was 21.0% (219/1042; 95%CI: 18.6–23.6%). The cumulative incidence standardized by age (categorical variable), race, and employment status (binary variable) was similar: 21.3% (95%CI: 19.5–23.2%). As part of the outreach survey, we identified 6 seroprevalence survey participants who tested negative for IgG antibodies but had a history of a positive PCR test, in agreement with the antibody test sensitivity estimates. This adjustment increased the crude cumulative incidence to 21.6% (225/1042; 95%CI: 19.1–24.2%) and the standardized cumulative incidence to 21.9% (95%CI: 20.1–23.9%).

Asymptomatic proportion.

The proportion of seroprevalence survey participants who did not exhibit any COVID-19 symptoms among those who tested positive for IgG antibodies was 16.4% (36/219; 95%CI: 11.8–22.0%). Among the confirmed COVID-19 cases identified among the outreach survey participants, the respective proportion adjusted for selection bias was very similar, but the confidence interval was wider because of the smaller confirmatory sample size: 16.8% (95%CI: 6.8–26.9%).

Combined seroprevalence and outreach surveys identified a total of 385 confirmed positive COVID-19 cases: 219 were identified via the seroprevalence survey, and an additional 166 cases were identified via outreach survey efforts. Because some of the outreach survey participants were also sampled in the seroprevalence study, the total number of confirmed COVID-19 cases in the outreach survey samples was 353. S2 Table in the S1 Appendix shows SARS-CoV-2 testing and COVID-19 history information for the main and confirmatory outreach samples. No significant differences were found between samples, except for a higher asymptomatic proportion in the confirmatory sample.

CHR.

The proportion of all confirmed COVID-19 cases hospitalized for COVID-19 was 9.4% (36/385; 95%CI: 6.6–12.7%). Among a subset of COVID-19 cases identified via the seroprevalence survey, the respective proportion was lower, at 5.5% (12/219; 95%CI: 2.9–9.4%). This estimate can be interpreted as an approximation of the infection hospitalization ratio (IHR), which is naturally lower than the CHR.

CFR and HFR.

A total of 10 deaths occurred in the merged study sample during the observation period. Of these, 8 deaths were participants with confirmed COVID-19, of which 7 deaths were associated with COVID-19. One patient had confirmed COVID-19 and recovered, as confirmed by a negative PCR test, but died soon after recovery from terminal cancer. Of the 7 patients with COVID-related deaths, 3 were hospitalized, and 4 died rapidly around the time of diagnosis without hospitalization. The CFR among all confirmed COVID-19 cases was 1.8% (7/385; 95%CI: 0.7–3.7%), and the HFR was 8.3% (3/36; 95%CI: 1.8–22.5%).

Dates of symptom onset

The symptom onset dates histogram shown in Fig 1 suggests that most cases occurred during the initial outbreak period. Although the first case in Long Island was officially reported on March 5th, these results suggest that some infections might have occurred as early as January and that several cases were present in late February. Similarly to the regional case count trajectory, the curve demonstrated an exponential increase in March and a rapid decrease starting approximately 2 weeks after the shutdowns were initiated. The curve followed a similar pattern among confirmed (orange) and unconfirmed (gray) symptomatic cases during the spring and summer. During the winter, when testing was not widely available, most symptomatic cases were classified as unconfirmed. Some of them might have been associated with influenza or other respiratory infections.

thumbnail
Fig 1. Dates of symptom onset according to COVID-19 status among the World Trade Center responder cohort in Long Island, NY between November 2019 and August 2020.

Confirmed positive cases are shown in orange (N = 273), and unconfirmed cases are shown in gray (N = 129). Symptomatic cases confirmed to be negative for SARS-CoV-2 are not shown.

https://doi.org/10.1371/journal.pone.0254713.g001

Correlates of infection

Table 3 shows results of multivariable logistic regression to identify correlates of confirmed infection. Confirmed positive COVID-19 was significantly associated with younger age, being black or Hispanic vs. non-Hispanic white, being currently employed, and being prescribed ibuprofen. Being multiracial/other vs non-Hispanic white was negatively associated with confirmed COVID-19.

thumbnail
Table 3. Multivariable logistic regression: Correlates of confirmed SARS-CoV-2 infection among World Trade Center responders with known COVID-19 status (event N/total N: 385/1,202)a.

https://doi.org/10.1371/journal.pone.0254713.t003

Discussion

In this study, we conducted an investigation of COVID-19 epidemiology in a cohort of middle-aged men and women who responded to 9/11/2001 WTC events and resided in Long Island, NY. We used an array of data collection and analytic tools intended to minimize the bias and provide estimates of key epidemiologic characteristics in this population.

We estimated that approximately 22% of WTC responders were infected with SARS-CoV-2 from January to August of 2020. Our estimate of cumulative incidence is consistent with that previously reported in Long Island as of March 29, 2020, which was based on a convenience sample and indicated evidence of prior infection among 13.2% of study participants [3].

Our estimate of cumulative incidence is based on a seroprevalence survey among a consecutive sample of WTC responders presenting for health monitoring visits, corrected on the basis of available PCR testing results in the sample, and standardized by age, race/ethnicity, and employment status. Whereas significant differences in these demographic characteristics were observed between the seroprevalence survey sample and the rest of the cohort, standardization resulted in a very small change in the cumulative incidence estimate compared to the crude sample proportion, thus suggesting that the bias of the seroprevalence survey sample is likely to be minimal. However, unmeasured confounding might still be present. One potential source of confounding may be associated with the study participants’ knowledge of prior infection, which might have been an additional motivation to present for routine health monitoring visit rather than postponing or canceling, in which case the cumulative incidence would have been overestimated. Potential sources that could lead to the underestimation of cumulative incidence include lower than reported test sensitivity and the potential for rapid degradation of antibodies (in fewer than 4 months) among asymptomatic individuals, both of which are unlikely. Similarly to other studies (e.g., [40]), we found lower antibody levels in asymptomatic than symptomatic infections, although the difference was not statistically significant.

Although we found that older age was associated with a lower probability of infection, this phenomenon is likely to be a behavioral rather than biologic. Given that many WTC responders are essential workers, including police officers, exposure to infection in this cohort was probably comparable to, or somewhat higher than, that in the general population. Yet, we also found a higher PCR-based case ascertainment proportion (57%) among seroprevalence survey participants than what had been estimated in the general population [4], which may potentially be associated with regular testing requirements at work and better access to health care among many cohort members, as well as better testing availability in late spring and summer of 2020 than in early spring.

We estimated that approximately 16% of the SARS-CoV-2 infections in this population were asymptomatic. This estimate is consistent with the 17.9% estimated from an outbreak on the Diamond Princess cruise ship, in which the population was comparable in age to the WTC cohort [41]. This estimate can be interpreted as a lower bound for the true asymptomatic proportion in this population, owing to the potential unmeasured confounders discussed above, thus helping to narrow the range of currently available estimates—a topic of extensive scientific debate [22, 42].

Our estimates of CHR, CFR, and HFR are consistent with previously reported values, adjusting for age when available [15, 23, 25, 4347]. This finding was somewhat unexpected, given that the population predominately comprised men with at least one comorbid condition and a high prevalence of pre-existing respiratory conditions; however uncertainty regarding CFR and HFR is high. Many positive cases in this study were confirmed with serology rather than PCR alone, thus potentially explaining lower CHR and CFR among confirmed cases [48]. At the same time, whereas 7 COVID-19 related deaths were reported, only 3 of those who died had previously been hospitalized, thus suggesting that one possible source for under-reporting of COVID-related mortality may have emerged if hospitals were considered central clearing houses for COVID-related deaths, particularly during the first wave of the pandemic [49]. Recognizing that COVID-19 might have caused many deaths before the implementation of the present study, we examined all death events in the parent cohort between March and August of 2020. The crude COVID-19 related mortality rate was slightly higher among the survey participants than the rest of the cohort, thereby suggesting that the study sample was unlikely to have been affected by survivor bias.

In the available sample, confirmed SARS-CoV-2 infection was correlated with demographic characteristics, including younger age, black or Hispanic race / ethnicity, and being currently employed rather than unemployed or retired. These associations are likely to be proxy measures for greater exposure probability. In our sample, none of the pre-existing conditions were associated with SARS-CoV-2 infection. We found that ibuprofen prescription was correlated with SARS-CoV-2 infection. Ibuprofen usage has been previously reported to be associated with greater COVID-19 severity [50], although the interpretation of this association remains unclear [51]. Further replication analyses are needed to determine whether these results are generalizable to other populations.

This is the first study of its kind in WTC responders, and one of the first studies to report results from a large non-hospital cohort globally. While this study avoids many limitations inherent to hospital-based cross-sectional or cohort studies, several limitations should still be noted. Although we made every effort to reach cohort members, approximately one-third of those contacted did not respond to requests for information. We addressed this issue by recruiting a confirmatory sample to correct for potential self-selection bias. A seroprevalence survey, which provided a basis for population-level estimates, was conducted over several months. Despite the similarity between the crude and standardized estimates of cumulative incidence, our results might still have unmeasured confounding. At the same time, antibody testing has been demonstrated to be more reliable than PCR testing in identifying infections. If residual confounding was present, it was likely to have been in the direction of overestimating the cumulative incidence and underestimating the asymptomatic proportion. Additional outreach efforts used to identify confirmed COVID-19 cases allowed for higher precision in the estimation of disease progression characteristics and assessment of the correlates of SARS-CoV-2 infection among WTC responders.

Conclusions

In this study, we sought to investigate the extent to which COVID-19 affected middle-aged individuals in a hard-hit part of New York State through the end of summer 2020. We found a cumulative incidence of approximately 22% and a possibly earlier community transmission than has been generally accepted in Long Island, NY. Higher infection risk was associated with employment and younger age–variables likely to indicate higher exposure probability. We also found that black and Hispanic individuals had a higher risk of contracting the virus than non-Hispanic white individuals. None of the pre-existing conditions that we investigated were associated with increased infection risk.

Supporting information

S1 Appendix. Detailed description of the sample, estimation methodology, descriptive statistics of the outreach surveys, additional regression analysis results, and survey instruments.

https://doi.org/10.1371/journal.pone.0254713.s001

(PDF)

Acknowledgments

We thank Dr. Evelyn J. Bromet for helpful comments and discussion.

References

  1. 1. Goyal P, Choi JJ, Pinheiro LC, Schenck EJ, Chen R, Jabri A, et al. Clinical Characteristics of Covid-19 in New York City. The New England journal of medicine. 2020;382(24):2372–4. pmid:32302078
  2. 2. Dong E, Du H, Gardner L. An interactive web-based dashboard to track COVID-19 in real time. The Lancet infectious diseases. 2020;20(5):533–4. pmid:32087114
  3. 3. Rosenberg ES, Tesoriero JM, Rosenthal EM, Chung R, Barranco MA, Styer LM, et al. Cumulative incidence and diagnosis of SARS-CoV-2 infection in New York. Annals of epidemiology. 2020;48:23–9 e4. pmid:32648546
  4. 4. Havers FP, Reed C, Lim T, Montgomery JM, Klena JD, Hall AJ, et al. Seroprevalence of Antibodies to SARS-CoV-2 in 10 Sites in the United States, March 23-May 12, 2020. JAMA Intern Med. 2020. pmid:32692365
  5. 5. Stadlbauer D, Tan J, Jiang K, Hernandez M, Fabre S, Amanat F, et al. Seroconversion of a city: Longitudinal monitoring of SARS-CoV-2 seroprevalence in New York City. medRxiv. 2020:2020.06.28.20142190.
  6. 6. Rigatti SJ, Stout R. SARS-CoV-2 Antibody Prevalence and Association with Routine Laboratory Values in a Life Insurance Applicant Population. medRxiv. 2020:2020.09.09.20191296.
  7. 7. Siordia JA Jr. Epidemiology and clinical features of COVID-19: A review of current literature. J Clin Virol. 2020;127:104357. pmid:32305884
  8. 8. Borges do Nascimento IJ, Cacic N, Abdulazeem HM, von Groote TC, Jayarajah U, Weerasekara I, et al. Novel Coronavirus Infection (COVID-19) in Humans: A Scoping Review and Meta-Analysis. J Clin Med. 2020;9(4). pmid:32235486
  9. 9. Park M, Cook AR, Lim JT, Sun Y, Dickens BL. A Systematic Review of COVID-19 Epidemiology Based on Current Evidence. J Clin Med. 2020;9(4). pmid:32244365
  10. 10. Bar-On YM, Flamholz A, Phillips R, Milo R. SARS-CoV-2 (COVID-19) by the numbers. Elife. 2020;9. pmid:32228860
  11. 11. Wolfel R, Corman VM, Guggemos W, Seilmaier M, Zange S, Muller MA, et al. Virological assessment of hospitalized patients with COVID-2019. Nature. 2020;581(7809):465–9. pmid:32235945
  12. 12. Han R, Huang L, Jiang H, Dong J, Peng H, Zhang D. Early Clinical and CT Manifestations of Coronavirus Disease 2019 (COVID-19) Pneumonia. AJR Am J Roentgenol. 2020;215(2):338–43. pmid:32181672
  13. 13. Kanberg N, Ashton NJ, Andersson LM, Yilmaz A, Lindh M, Nilsson S, et al. Neurochemical evidence of astrocytic and neuronal injury commonly found in COVID-19. Neurology. 2020;95(12):e1754–e9. pmid:32546655
  14. 14. Centers for Disease Control and Prevention. CDC COVID Data Tracker: https://covid.cdc.gov/covid-data-tracker [accessed: December 2, 2020].
  15. 15. Verity R, Okell LC, Dorigatti I, Winskill P, Whittaker C, Imai N, et al. Estimates of the severity of coronavirus disease 2019: a model-based analysis. The Lancet infectious diseases. 2020;20(6):669–77. pmid:32240634
  16. 16. Lauer SA, Grantz KH, Bi Q, Jones FK, Zheng Q, Meredith HR, et al. The Incubation Period of Coronavirus Disease 2019 (COVID-19) From Publicly Reported Confirmed Cases: Estimation and Application. Annals of internal medicine. 2020;172(9):577–82. pmid:32150748
  17. 17. Sood N, Simon P, Ebner P, Eichner D, Reynolds J, Bendavid E, et al. Seroprevalence of SARS-CoV-2-Specific Antibodies Among Adults in Los Angeles County, California, on April 10–11, 2020. JAMA: the journal of the American Medical Association. 2020;323(23):2425–7. pmid:32421144
  18. 18. Menachemi N, Yiannoutsos CT, Dixon BE, Duszynski TJ, Fadel WF, Wools-Kaloustian KK, et al. Population Point Prevalence of SARS-CoV-2 Infection Based on a Statewide Random Sample—Indiana, April 25–29, 2020. MMWR Morb Mortal Wkly Rep. 2020;69(29):960–4. pmid:32701938
  19. 19. Biggs HM, Harris JB, Breakwell L, Dahlgren FS, Abedi GR, Szablewski CM, et al. Estimated Community Seroprevalence of SARS-CoV-2 Antibodies—Two Georgia Counties, April 28-May 3, 2020. MMWR Morb Mortal Wkly Rep. 2020;69(29):965–70. pmid:32701941
  20. 20. Bendavid E, Mulaney B, Sood N, Shah S, Ling E, Bromley-Dulfano R, et al. COVID-19 Antibody Seroprevalence in Santa Clara County, California. medRxiv. 2020:2020.04.14.20062463.
  21. 21. Oran DP, Topol EJ. Prevalence of Asymptomatic SARS-CoV-2 Infection: A Narrative Review. Annals of internal medicine. 2020;173(5):362–7. pmid:32491919
  22. 22. Meyerowitz EA, Richterman A, Bogoch II, Low N, Cevik M. Towards an accurate and systematic characterisation of persistently asymptomatic infection with SARS-CoV-2. The Lancet infectious diseases. 2020. pmid:33301725
  23. 23. Bi Q, Wu Y, Mei S, Ye C, Zou X, Zhang Z, et al. Epidemiology and transmission of COVID-19 in 391 cases and 1286 of their close contacts in Shenzhen, China: a retrospective cohort study. The Lancet infectious diseases. 2020;20(8):911–9. pmid:32353347
  24. 24. Onder G, Rezza G, Brusaferro S. Case-Fatality Rate and Characteristics of Patients Dying in Relation to COVID-19 in Italy. JAMA: the journal of the American Medical Association. 2020;323(18):1775–6. pmid:32203977
  25. 25. Lane EA, Barrett DJ, Casey M, McAloon CG, Collins AB, Hunt K, et al. Country differences in hospitalisation, length of stay, admission to Intensive Care Units, and mortality due to SARS-CoV-2 infection at the end of the first wave in Europe: a rapid review of available literature. medRxiv. 2020:2020.05.12.20099473.
  26. 26. Baggett TP, Keyes H, Sporn N, Gaeta JM. Prevalence of SARS-CoV-2 Infection in Residents of a Large Homeless Shelter in Boston. JAMA: the journal of the American Medical Association. 2020;323(21):2191–2. pmid:32338732
  27. 27. Arons MM, Hatfield KM, Reddy SC, Kimball A, James A, Jacobs JR, et al. Presymptomatic SARS-CoV-2 Infections and Transmission in a Skilled Nursing Facility. The New England journal of medicine. 2020;382(22):2081–90. pmid:32329971
  28. 28. Payne DC, Smith-Jeffcoat SE, Nowak G, Chukwuma U, Geibe JR, Hawkins RJ, et al. SARS-CoV-2 Infections and Serologic Responses from a Sample of U.S. Navy Service Members—USS Theodore Roosevelt, April 2020. MMWR Morb Mortal Wkly Rep. 2020;69(23):714–21. pmid:32525850
  29. 29. Gomez-Ochoa SA, Franco OH, Rojas LZ, Raguindin PF, Roa-Diaz ZM, Wyssmann BM, et al. COVID-19 in Healthcare Workers: A Living Systematic Review and Meta-analysis of Prevalence, Risk Factors, Clinical Characteristics, and Outcomes. Am J Epidemiol. 2020.
  30. 30. Anand S, Montez-Rath M, Han J, Bozeman J, Kerschmann R, Beyer P, et al. Prevalence of SARS-CoV-2 antibodies in a large nationwide sample of patients on dialysis in the USA: a cross-sectional study. Lancet. 2020.
  31. 31. Robilotti EV, Babady NE, Mead PA, Rolling T, Perez-Johnston R, Bernardes M, et al. Determinants of COVID-19 disease severity in patients with cancer. Nat Med. 2020;26(8):1218–23. pmid:32581323
  32. 32. Clouston SAP, Kuan P, Kotov R, Mukherjee S, Thompson-Carino P, Bromet EJ, et al. Risk factors for incident prostate cancer in a cohort of world trade center responders. BMC Psychiatry. 2019;19(1):389. pmid:31822278
  33. 33. Clouston SAP, Diminich ED, Kotov R, Pietrzak RH, Richards M, Spiro A 3rd, et al. Incidence of mild cognitive impairment in World Trade Center responders: Long-term consequences of re-experiencing the events on 9/11/2001. Alzheimers Dement (Amst). 2019;11:628–36. pmid:31517025
  34. 34. Baraniuk C. Covid-19 antibody tests: a briefing. BMJ. 2020;369:m2284. pmid:32580937
  35. 35. Centers for Disease Control and Prevention. Symptoms of Coronavirus: https://www.cdc.gov/coronavirus/2019-ncov/symptoms-testing/symptoms.html [accessed February 14, 2021].
  36. 36. Roalfe AK, Holder RL, Wilson S. Standardisation of rates using logistic regression: a comparison with the direct method. BMC health services research. 2008;8(1):1–7. pmid:19113996
  37. 37. Hope AC. A simplified Monte Carlo significance test procedure. Journal of the Royal Statistical Society: Series B (Methodological). 1968;30(3):582–98.
  38. 38. Cochran WG. Sampling techniques: John Wiley & Sons; 2007.
  39. 39. R Development Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing; 2013.
  40. 40. Long QX, Liu BZ, Deng HJ, Wu GC, Deng K, Chen YK, et al. Antibody responses to SARS-CoV-2 in patients with COVID-19. Nat Med. 2020;26(6):845–8. pmid:32350462
  41. 41. Mizumoto K, Kagaya K, Zarebski A, Chowell G. Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID-19) cases on board the Diamond Princess cruise ship, Yokohama, Japan, 2020. Euro Surveill. 2020;25(10).
  42. 42. Various authors. Comments to “Oran DP, Topol EJ. Prevalence of Asymptomatic SARS-CoV-2 Infection: A Narrative Review. Annals of internal medicine 2020;173(5):362–67” https://www.acpjournals.org/doi/full/10.7326/M20-3012 [accessed 7-October-2020].
  43. 43. Wu JT, Leung K, Bushman M, Kishore N, Niehus R, de Salazar PM, et al. Addendum: Estimating clinical severity of COVID-19 from the transmission dynamics in Wuhan, China. Nat Med. 2020;26(7):1149–50. pmid:32661399
  44. 44. Levin AT, Hanage WP, Owusu-Boaitey N, Cochran KB, Walsh SP, Meyerowitz-Katz G. Assessing the age specificity of infection fatality rates for COVID-19: Meta-analysis & public policy implications. National Bureau of Economic Research, 2020 0898–2937. pmid:33289900
  45. 45. Rinaldi G, Paradisi M. An empirical estimate of the infection fatality rate of COVID-19 from the first Italian outbreak. medRxiv. 2020:2020.04.18.20070912.
  46. 46. Ward H, Atchison CJ, Whitaker M, Ainslie KEC, Elliott J, Okell LC, et al. Antibody prevalence for SARS-CoV-2 in England following first peak of the pandemic: REACT2 study in 100,000 adults. medRxiv. 2020:2020.08.12.20173690.
  47. 47. Grewelle R, De Leo G. Estimating the Global Infection Fatality Rate of COVID-19. medRxiv. 2020:2020.05.11.20098780.
  48. 48. Liang L-L, Tseng C-H, Ho HJ, Wu C-Y. Covid-19 mortality is negatively associated with test number and government effectiveness. Scientific reports. 2020;10(1):1–7.
  49. 49. Woolf SH, Chapman DA, Sabo RT, Weinberger DM, Hill L, Taylor DDH. Excess Deaths From COVID-19 and Other Causes, March-July 2020. JAMA: the journal of the American Medical Association. 2020;324(15):1562–4. pmid:33044483
  50. 50. Day M. Covid-19: ibuprofen should not be used for managing symptoms, say doctors and scientists. British Medical Journal Publishing Group; 2020.
  51. 51. Sodhi M, Etminan M. Safety of Ibuprofen in Patients With COVID-19: Causal or Confounded? Chest. 2020;158(1):55–6. pmid:32243944