key: cord-0764107-2cai67qb
authors: Lieberman-Cribbin, Wil; Galanti, Marta; Shaman, Jeffrey
title: Socioeconomic Disparities in Severe Acute Respiratory Syndrome Coronavirus 2 Serological Testing and Positivity in New York City
date: 2021-10-17
journal: Open Forum Infect Dis
DOI: 10.1093/ofid/ofab534
sha: 30c6147e6d4475afa0bdf9b22746c69fd9c2baa8
doc_id: 764107
cord_uid: 2cai67qb

BACKGROUND: We characterized severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) antibody test prevalence and positive test prevalence across New York City (NYC) in order to investigate disparities in testing outcomes by race and socioeconomic status (SES). METHODS: Serologic data were downloaded from the NYC Coronavirus data repository (August 2020–December 2020). Area-level characteristics for NYC neighborhoods were downloaded from United States census data and a socioeconomic vulnerability index was created. Spatial generalized linear mixed models were performed to examine the association between SES and antibody testing and positivity. RESULTS: The proportion of Hispanic population (posterior median, 0.001 [95% credible interval, 0.0003–0.002]), healthcare workers (0.003 [0.0001–0.006]), essential workers (0.003 [0.001–0.005]), age ≥65 years (0.003 [0.00002–0.006]), and high SES (SES quartile 3 vs 1: 0.034 [0.003–0.062]) were positively associated with antibody tests per 100000 residents. The White proportion (–0.002 [–0.003 to –0.001]), SES index (quartile 3 vs 1, –0.068 [–0.115 to –0.017]; quartile 4 vs 1, –0.077 [–0.134 to –0.018]) and age ≥65 years (–0.005 [–0.009 to –0.002]) were inversely associated with positive test prevalence (%), whereas the Hispanic (0.004 [0.002–0.006]) and essential worker (0.008 [0.003–0.012]) proportions had positive coefficients. CONCLUSIONS: Disparities in serologic testing and seropositivity exist on SES and race/ethnicity across NYC, indicative of excess coronavirus disease burden in vulnerable and marginalized populations.

Following its first reported case of novel coronavirus disease 2019 (COVID-19) on 1 March 2020 [1] , New York City (NYC) quickly became the initial US epicenter of COVID-19. Antibody testing became widely available shortly thereafter in April 2020 and is used to confirm past infection with the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19 [2] . Antibody testing continues to be conducted among those with or without symptomatology [2, 3] . These test results are useful both at individual and population-level scales, and can be used to measure exposure to SARS-CoV-2 within different geographies and age and populations groups, gain information for contact tracing, and identify possible population sources for monoclonal antibody treatments [3] . However, the longevity of SARS-CoV-2 antibodies after infection remains uncertain, and more definitive quantification of adaptive immunity is needed to inform future testing, vaccination efforts, reinfection risk, and reopening plans.

Within NYC, there have been multiple efforts to quantify seroprevalence, although these have often been restricted to healthcare settings and occurred early in the pandemic. One study in the Northwell Health System studied immunoglobulin G (IgG) antibodies among 46 000 personnel, with testing performed between March and June 2020 [4] . They reported a seroprevalence of 13.7%, which was similar to results of SARS-CoV-2 IgG testing (14.0%) performed randomly across New York State during April 2020 [5] . However, higher seroprevalence estimates were reported in other healthcare populations. Among healthcare workers in the South Bronx, seroprevalence of IgG antibodies was 27% in May 2020 [6] . Subsequent estimates of the general NYC population found similar seroprevalence estimates. In a May-July 2020 sample of 45 367 NYC residents, the prevalence of IgG antibodies was 23.6% [7] , while a June-October 2020 study relying on both blood samples to test for IgG antibodies and self-reported antibody test results found a prevalence of 24.3% [8] .

Data from NYC hospital centers early in the pandemic revealed racial disparities in incidence and hospitalizations, with Black patients (1.3-fold) and Hispanic patients (1.5-fold) more likely to test positive, and Black patients 1.89-fold more likely to require hospitalization than White patients [9] . Hispanic patients were also found to be at increased risk for in-hospital mortality (1.84-fold increase) compared to Whites [10] . These findings corroborate both national-level racial variations in mortality [11] and hospitalizations [12] , as well as populationlevel data from NYC [13, 14] . To date, the age-adjusted case, hospitalization, and death rates in NYC are all highest among Black and Hispanic/Latino groups [14] . Moreover, a cross-sectional spatial analysis conducted in April 2020 reported that in areas with lower median income, a greater percentage of individuals who identify as non-White and/or Hispanic/Latino, a greater percentage of essential workers, and a greater percentage of healthcare workers had greater subway use during the early pandemic [15] , suggesting a greater potential for SARS-CoV-2 exposure. They also reported that the essential worker population was a main driver of subway use in lower socioeconomic status (SES) areas and communities of color, suggesting an interplay among multiple sociodemographic factors affecting COVID-19 spread.

To date, there have been few studies investigating racial disparities in seroprevalence [7, 16] , and those that do exist are not NYC-specific. Among a national sample of dialysis patients, non-Hispanic Black and Hispanic neighborhoods were more likely to be seropositive than patients from non-Hispanic White areas [17] . This finding was also observed in regional areas of the US in both healthcare settings and the general population. Among healthcare workers in the Midwest, Black, Asian, and mixed-race workers were more likely to be seropositive than White workers [18] , and a seroprevalence study conducted in Baton Rouge, Louisiana, found that seroprevalence was highest (7.5%) among Black participants [19] .

Seroprevalence testing has emerged as a central means of characterizing the burden of SARS-CoV-2 infection and for determining how infections differ by SES and area-level characteristics. To investigate the presence of potential racial and socioeconomic disparities with respect to seroprevalence in NYC, we characterize the spatial and temporal distribution of antibody test prevalence and antibody positivity prevalence and investigate the relationship between community-level characteristics and testing outcomes.

The weekly number of SARS-CoV-2 antibody tests performed and the number of positive test results were downloaded from the NYC Coronavirus (COVID-19) data repository hosted by the NYC Department of Health and Mental Hygiene (DOHMH) from August 2020 through 10 December 2020 [20] . Data were downloaded at the modified zip (postal) code tabulation area (MODZCTA) level, which are areas that approximate an individual's zip code. As used here, ZCTA refers to a person's residence, and not the location where testing occurred. Data in the repository include people of all ages and include any type of serologic test reported to NYC DOHMH. The date of first data release in the repository by MODZCTA was 20 August 2020, which includes data from 5 April-20 August 2020. Cumulative data after 20 August 2020 were downloaded from the NYC Coronavirus data repository weekly.

An index of socioeconomic vulnerability was created at ZCTA level from 2018 American Community Survey 5-year estimates [21] . This index has been utilized previously [22] to describe area-level SES with respect to SARS-CoV-2 testing and includes data on median household income during the past 12 months, median gross rent, percentage living on income <150% of the poverty line, education, percentage working class, percentage unemployed, and percentage living with >1 occupant per room. The education measure was calculated based from the population aged ≥25 years and is a weighted combination of the population percentage holding a high school degree, only a high school degree, and more than a high school degree, with a greater value indicating higher educational attainment.

The variables used in this analysis differ from the original set [22] through incorporation of the occupants per room measure while removing the median home value to make the index more relevant to NYC and SARS-CoV-2. These variables were then combined using principal component analysis and the first eigenvector was used to create an SES index score for each ZCTA. ZCTAs were classified into quartiles based on this SES score, with a score of 1 corresponding to the lowest SES quartile (ie, lowest resourced) and a score of 4 representing the highest SES quartile (highest resourced).

In addition to the variables included in the SES index, the racial (White alone proportion), Hispanic proportion, essential worker proportion, healthcare essential worker proportion, population with insurance coverage, and population >65 years of age for each ZCTA were obtained from 2019 American Community Survey 5-year estimates. The percentage of nonhealthcare essential workers in a ZCTA was defined by the fraction of the civilian employed population ≥16 years employed in agriculture, forestry, fishing and hunting, mining, construction, manufacturing, wholesale trade, retail trade, transportation and warehousing, and utilities, whereas the percentage of healthcare essential workers was defined by healthcare practitioners and healthcare technical occupations, as described previously [15] .

Two-week periods of SARS-CoV-2 antibody test prevalence and positivity were compared across SES index quartiles with the spatial distribution of test prevalence and positivity mapped according to ZCTA of residence. Cluster analysis was performed using Anselin Local Moran I over the entire study period to determine clusters of high-test prevalence and positive prevalence. The presence of spatial autocorrelation among outcome variables was investigated with the Moran I test. To account for the spatial structure of our data, a spatial weights matrix was created where each ZCTA polygon was treated as a unique feature. A queen contiguity structure was utilized, where each touching ZCTA was treated as a neighbor, with neighbors considered more closely related than nonneighbors.

Spatial generalized linear mixed models accounting for random effects with a Besag-York-Mollie conditional autoregressive prior were employed to examine the association between SES and antibody testing and positivity. This was performed using the CARBayes package in R [23] , which uses Markov chain Monte Carlo simulations to generate model estimates. The model was applied, in turn, to all combinations of predictor variables (SES index, White proportion, Hispanic proportion, healthcare worker proportion, nonessential worker proportion, aged ≥65 years, insurance coverage); 45 000 samples were generated for each model, with a burn-in of 20 000 samples. In the model predicting antibody tests, the population of the ZCTA was used as an offset, whereas the total number of tests was used as an offset for predicting antibody-positive test results. Goodness-of-fit was assessed using the Watanabe-Akaike information criterion (wAIC).

The candidate models were then subject to multimodel inference in order to identify key, common explanatory features. To conduct the multimodel inference, candidate model weights were calculated according to the following formula:

and minwAIC is the wAIC value of the model with the lowest wAIC among all candidate models. The smallest subset of models with a total weight w j ≥ 0.95 was used for model averaging. The relative importance, or the sum of wAIC weights across models that contain a given predictor variable, as well as the probability this predictor was included in models, were also calculated. All analyses were performed in R version 1.4.1106.

Across all SES quartiles, antibody tests per 100 000 residents in each ZCTA decreased during the study period, whereas test positivity increased (Figure 1 There are 177 MODZCTAs in NYC; the spatial distribution of antibody tests performed (Figure 2A ) and positive test prevalence ( Figure 2B ) are displayed for each 2-week period. The highest rates of testing per 100 000 people in the ZCTA were during 20 August-3 September 2020 and 3 September-17 September 2020, when most testing was performed in Brooklyn, central Queens, and the Bronx. By 26 November-10 December 2020, testing decreased across NYC as a whole, but remained highest in southern Brooklyn, central Queens, and the Bronx. Table 1) . The individual best fitting model ( Table 2) Variable importance was similar across model averages for both testing and positivity outcomes ( Table 3 ). In the models of antibody testing for which cumulative wAIC weight ≥95%, insurance coverage and White proportion appeared in the greatest number of models (59 and 58, respectively) and had the highest weights sum (0.52). In the models of positivity for which cumulative wAIC weight ≥95%, Hispanic proportion appeared in the most models (62) and had the highest weights sum (0.67); White proportion (58 models, 0.62) and essential worker proportion (59 models, 0.59) were also important predictors.

The findings of this study show that the number of antibody tests performed and positive antibody test results for SARS-CoV-2 varied across NYC by ZCTA, and that SES, proportion of essential workers, the proportion of population aged ≥65 years, and the proportion of the population identifying as Hispanic were important predictors for testing and positivity outcomes. Specifically, areas with higher proportions of Hispanic residents and essential workers were more likely to have increased SARS-CoV-2 testing and seropositivity. Although areas with increased proportions of residents aged ≥65 years experienced more testing, these same areas experienced fewer positive test results. There was a notable trend in SES in which higher SES areas had a lower rate of antibody positivity, despite the absence of a statistically significant relationship between SES and testing.

This analysis corroborates findings identifying racial disparities in COVID-19 cases and deaths, both nationally [11, [24] [25] [26] and within NYC [4, 6, 27] , underscoring the vulnerability of marginalized, working-class, and lower SES populations. While seroprevalence studies have been performed in NYC using patient or hospital-level data [7, 9, 28] , NYCwide seroprevalence studies are needed to gain a picture of the larger NYC population [7, 8] . Previous serosurveys conducted across NYC [7, 8] reported greatest seroprevalence among Black and Hispanic respondents, people from highpoverty areas, and respondents employed in healthcare or essential worker industries. However, these results were based on a smaller number of surveyed residents. Our results are consistent with what has been previously reported, and are more comprehensive, utilizing city-wide data from a longer time period. The interplay between SES, essential worker population, and Hispanic composition in New York City is highlighted by examining the distribution of these variables in light of seroprevalence outcomes (Supplementary Figure 1) . Regions of high positivity such as southern Brooklyn have a low Hispanic composition, higher rates of essential workers, and average SES index scores (Supplementary Figure 2) . This is in contrast to high-positivity areas in the Bronx, which have low SES scores and high Hispanic composition, or those in central Queens, which also have majority Hispanic populations and average SES index scores. However, these 3 regions, in particular, all have a high proportion of essential worker populations, suggesting that this variable is relevant for these 3 communities. The diversity of NYC is indicative of the importance of multiple predictor variables that determine the distribution of antibody testing and positivity in NYC.

One limitation of this analysis is that the dataset used only partially captures antibody testing in NYC, as a number of These data also reflect changes in the accessibility and attitudes toward diagnostic and serologic testing over time. While these tests are currently free and easily available in NYC, this was not the case early in the pandemic, and in spring 2020 diagnostic testing was not encouraged in those with mild/moderate symptoms [24] . This analysis uses American Community Survey data from 2018 and 2019, which may not be reflective of community characteristics during 2020. As this was an ecological study, there are also likely other individual-level factors that may influence seeking SARS-CoV-2 antibody testing that are not accounted for here. It is also possible that populations vulnerable to COVID-19 disease are unable to seek testing due to work restrictions or potential loss of income, or are less likely to seek antibody testing as opposed to diagnostic testing. However, the association between higher SES and reduced positivity found in this study corroborates evidence of increased burden of COVID-19 in vulnerable populations that have been previously documented. This research does not study immune memory over time and does not intend to shed light on what fraction of the population previously exposed to SARS-CoV-2 is protected from reinfection. In trying to understand the risk of reinfection, there needs to be more comprehensive measurement of long-term immunity in the NYC population. However, identifying disparities in the prior burden of SARS-CoV-2 infection can be used to support future vaccination efforts in areas at greatest risk for further SARS-CoV-2 infection and adverse outcomes.

We found disparities in the burden of cumulative SARS-CoV-2 infection within NYC through the end of 2020, with lower SES communities and essential worker populations at greater risk of seropositivity. Future studies should not only use serological testing in NYC to estimate the extent of COVID-19 disease burden but also investigate why these disparities, such as access to testing, exist. Analyses of seroprevalence at the population level in tandem with research on the longevity of immune memory after infection can also indicate areas at risk of potential reinfection and areas to focus vaccination efforts.

Supplementary materials are available at Open Forum Infectious Diseases online. Consisting of data provided by the authors to benefit the reader, the posted materials are not copyedited and are the sole responsibility of the authors, so questions or comments should be addressed to the corresponding author.

Patient consent statement. This study does not include factors necessitating patient consent as individual-level data were not used.

Disclaimer. The funders had no role in study design, data collection and analysis, preparation of the manuscript, or decision to submit the manuscript for publication.

Financial support. This work was funded by the National Institute of Environmental Health Sciences (grant number T32 ES007322); the National Institute of Allergy and Infectious Diseases (grant number AI163023); and the Morris-Singer Foundation.

Potential conflicts of interest. J. S. and Columbia University disclose ownership of SK Analytics and J. S. discloses personal fees from BNI (Business Network International). All other authors report no potential conflicts of interest.

World Health Organization. Coronavirus disease 2019 (COVID-19) situation report-92

Testing for COVID-19

The promise and peril of antibody testing for COVID-19

Northwell Health COVID-19 Research Consortium. Prevalence of SARS-CoV-2 antibodies in health care personnel in the New York city area

Cumulative incidence and diagnosis of SARS-CoV-2 infection in New York

SARS-CoV-2 seroprevalence among health care workers in a New York City hospital: a cross-sectional analysis during the COVID-19 pandemic

SARS-CoV-2 following the largest initial epidemic wave in the United States: findings from

Prevalence of SARS-CoV-2 antibodies in New York City adults

Assessment of racial/ethnic disparities in hospitalization and mortality in patients with COVID-19 in New York city

Racial disparities in COVID-19 hospitalization and in-hospital mortality at the height of the New York City pandemic

Variation in racial/ethnic disparities in COVID-19 mortality by age in the United States: a cross-sectional study

Characteristics associated with racial/ethnic disparities in COVID-19 outcomes in an academic health care system

Unequal burdens: assessing the determinants of elevated COVID-19 case and death rates in New York City's racial/ethnic minority neighbourhoods

COVID-19 data

Socioeconomic disparities in subway use and COVID-19 outcomes

SARS-CoV-2 seroprevalence among parturient women in Philadelphia

Prevalence of SARS-CoV-2 antibodies in a large nationwide sample of patients on dialysis in the USA: a cross-sectional study

Disparities in seroprevalence of SARS-CoV-2 immunoglobulin antibodies in a large midwestern health care system

Racial and workplace disparities in seroprevalence of SARS-CoV-2

New York City Department of Health. NYC coronavirus disease 2019 (COVID-19) data

American Community Survey 5-year estimates

Socioeconomic status and breast cancer incidence in California for different race/ethnic groups

CARBayes: an R package for Bayesian Spatial Modeling with conditional autoregressive priors

Patterns of COVID-19 testing and mortality by race and ethnicity among United States veterans: a nationwide cohort study

COVID-19 and racial/ethnic disparities

Racial and ethnic health disparities related to COVID-19

Disparities in COVID-19 testing and positivity in New York City

Repeated cross-sectional sero-monitoring of SARS-CoV-2 in New York City

All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Conflicts that the editors consider relevant to the content of the manuscript have been disclosed.