Quantification of occupational and community risk factors for SARS-CoV-2 seropositivity among healthcare workers in a large U.S. healthcare system
==================================================================================================================================================

* Julia M. Baker
* Kristin N. Nelson
* Elizabeth Overton
* Benjamin A. Lopman
* Timothy L. Lash
* Mark Photakis
* Jesse T. Jacob
* John Roback
* Scott K. Fridkin
* James P. Steinberg

## Abstract

**Background** Quantifying occupational risk factors for SARS-CoV-2 infection among healthcare workers can inform efforts to improve healthcare worker and patient safety and reduce transmission. This study aimed to quantify demographic, occupational, and community risk factors for SARS-CoV-2 seropositivity among healthcare workers in a large metropolitan healthcare system.

**Methods** We analyzed data from a cross-sectional survey conducted from April through June of 2020 linking risk factors for occupational and community exposure to COVID-19 with SARS-CoV-2 seropositivity. A multivariable logistic regression model was fit to quantify risk factors for infection. Participants were employees and medical staff members who elected to participate in SARS-CoV-2 serology testing offered to all healthcare workers as part of a quality initiative, and who completed a survey on exposure to COVID-19 and use of personal protective equipment. Exposures of interest included known demographic risk factors for COVID-19, residential zip code incidence of COVID-19, occupational exposure to PCR test-positive healthcare workers or patients, and use of personal protective equipment. The primary outcome of interest was SARS-CoV-2 seropositivity.

**Results** SARS-CoV-2 seropositivity was estimated to be 5.7% (95% CI: 5.2%-6.1%) among 10,275 healthcare workers. Community contact with a person known or suspected to have COVID-19 (aOR=1.9, 95% CI:1.4-2.5) and zip code level COVID-19 incidence (aOR: 1.4, 95% CI: 1.0-2.0) increased the odds of infection. Black individuals were at high risk (aOR=2.0, 95% CI:1.6-2.4). Overall, occupational risk factors accounted for 27% (95% CI: 25%-30%) of the risk among healthcare workers and included contact with a PCR test-positive healthcare worker (aOR=1.2, 95% CI:1.0-1.6).

**Conclusions** Community risk factors, including contact with a COVID-19 positive individual and residential COVID-19 incidence, are more strongly associated with SARS-CoV-2 seropositivity among healthcare workers than exposure in the workplace.

## Introduction

Healthcare workers (HCWs) are presumed to be at high risk for coronavirus disease 2019 (COVID-19) through occupational exposure to infected patients or coworkers. Studies have reported a wide range of seroprevalence of SARS-CoV-2, the virus that causes COVID-19, among HCWs. This variation has in part been attributed to differential risk of exposure related to COVID-19 incidence in the community.1 Indeed, recent studies have shown that a substantial number of infections among HCWs could not be traced to occupational exposures 2 and that community exposures were as or more strongly associated with infection.3 While previous studies have compared seroprevalence in HCWs with that of the general population, few have rigorously considered workplace risk factors alongside community risk factors among HCWs to estimate their relative contribution to overall infection risk.4,5

Accounting for the role of community risk, which may be large, is especially important because reports of occupational risk factors for SARS-CoV-2 infection among HCWs have been inconsistent.3-8 While some studies have shown weak associations between involvement in clinical care, care of COVID-19 patients, and exposure to coworkers with COVID-19, 6-8 others have shown that these are, in fact, risk factors for infection.4,5 Moreover, previous studies have not accounted for potential participation bias, though have cited it as a major limitation.1,2,4,7 Adjusting for participation bias while considering both workplace and community risk factors for infection can bring us closer to an accurate understanding of which workplace exposures confer the highest risk of infection for HCWs. This information can inform strategies to protect HCWs as the COVID-19 pandemic continues.

We aimed to quantify occupational, community, and demographic risk factors for SARS-CoV-2 seropositivity among HCWs in a large university-based healthcare system, adjusting for possible bias due to voluntary participation in testing. We estimate the proportion of risk attributable to occupational factors among HCWs.

## Methods

We analyzed data from Emory Healthcare, which includes 11 hospitals, 250 provider locations, and approximately 25,000 employees and medical staff members based in the Atlanta, Georgia, metropolitan area.

A voluntary serological survey was conducted among employees and medical staff members from April 19 through June 26, 2020 to inform process improvement activities. At the time of testing, HCWs completed a survey describing use of personal protective equipment (PPE) and possible exposure to COVID-19 inside and outside the workplace. These data were combined with employee demographic data. To assess community exposure to COVID-19, we used data from the Georgia Department of Public Health to calculate COVID-19 cumulative incidence by zip code and week. As a measure of exposure in their community, we assigned each participant the COVID-19 cumulative incidence in their zip code of residence two weeks prior to their test date to account for the lag from the time of infection to seroconversion (eFigure 1 in Supplement). The serologic test used to analyze participant samples was developed at the Emory Medical Laboratory and measures IgG antibody to the receptor binding domain of the SARS-CoV-2 spike protein.9

We fit a logistic regression model to estimate adjusted odds ratios (aORs) between potential risk factors and SARS-CoV-2 seropositivity. Model predictors (eTable 1 in Supplement) included demographics (age group; race; ethnicity), community exposure (cumulative incidence of COVID-19 by residential zip code two weeks before testing date; contact with confirmed/suspected COVID-19 cases outside the workplace), and occupational factors (workplace role and location; contact with COVID-positive patients and staff; use of PPE). Workplace location was self-reported and described the location(s) where a HCW spent the most time during the COVID-19 pandemic. Since over one-third (36.4%) of HCWs reported spending equal time in more than one location, the workplace location variable was structured as a hierarchy. We ranked locations based on the anticipated risk of encounters with COVID-19 positive patients (see eTable 1 in Supplement for details) with high encounter risk locations, such as the emergency department and COVID-19 focused units, highest and no patient contact and working from home ranked lowest. HCWs were categorized into the highest-ranking location where they reported working.

To account for potential selective participation, we used inverse probability of participation weighting to assess whether our results were sensitive to differences between those tested and all employees by age, race, and sex. We used the joint distribution of the age, race, and sex of all employees to weight individuals in the regression analysis such that demographic groups who were over-represented among survey participants compared to all employees were down-weighted, and vice versa. We calculated aORs using the weighted data and 95% confidence intervals (CIs) accounting for the weights.

Lastly, the fraction of infection attributable to occupational risk factors was estimated using a model-based approach to calculating the population attributable fraction for the combined set of all occupational risk factors.10 The 95% CIs for the attributable fraction was estimated by applying a normal distribution to the predicted probability of seropositivity for each survey participant and generating 1,000 bootstrap samples from this distribution. From these samples, we calculated the resulting distribution of the attributable fraction to calculate a 95% CI around the median attributable fraction. Analyses were conducted in R version 4.0.2 using the *survey, boot* and *zipcode* packages. The code for this analysis is available at [https://github.com/lopmanlab/emory-hcw-serosurvey](https://github.com/lopmanlab/emory-hcw-serosurvey).

This study was approved by the Emory University Institutional Review Board.

## Results

Among 10,275 participating HCWs (35% of the workforce), the SARS-CoV-2 seroprevalence was 5.7% (95%CI: 5.2%-6.1%). A total of 665 participants were missing zip code data, leaving 9,610 participants available for further analysis.

The majority (71%) of participants were tested in May 2020 (eTable 2 in Supplement). Over three-quarters of participants were female (77.6%) with female and male participants having similar seropositivity (5.7% and 5.9%, respectively). Black HCW were notably underrepresented in the serosurvey, comprising 30% of those who volunteered for antibody testing and 49% of the workforce. Prior confirmed COVID-19 infection was reported for 133 participants (1.4%). Forty-four percent (245/555) of seropositive HCWs reported no fever or COVID-19-like symptoms since February 1, 2020.

Contact with a person known or suspected to have COVID-19 outside the workplace (aOR=1.9, 95% CI:1.4-2.5) and race were most strongly associated with seropositivity (Table 1). Black (aOR=2.0, 95%CI:1.6-2.4) and multiracial (aOR:1.7, 95%CI:0.8-3.4) HCWs had higher odds of infection than white HCWs. Higher residential zip code incidence of COVID-19 infection (aOR=1.4, 95%CI:1.0-2.0) was also associated with seropositivity. Age under 30 years was associated with seropositivity (aOR 1.3, 95% CI 0.9, 1.8).

View this table:
[Table 1.](http://medrxiv.org/content/early/2020/11/03/2020.10.30.20222877/T1)

Table 1. Association of demographic, community, and occupational risk factors with seropositivity in unadjusted and adjusted analyses (n=9,610).

In the workplace, participants who reported close contact with a COVID-positive HCW had increased odds of seropositivity (aOR=1.2, 95% CI:1.0-1.6). Although HCWs who reported caring for patients with COVID-19 did not have increased odds of seropositivity (aOR=0.9, 95% CI: 0.7, 1.3), working in clinical locations, such as inpatient non-COVID-19 focused areas (aOR=1.3, 95% CI: 0.8, 2.1), the emergency department (aOR=1.4, 95% CI: 0.9, 2.4) or COVID-focused units (aOR=1.5, 95% CI: 0.9, 2.4), was associated with higher odds of seropositivity. Respiratory therapists (aOR=0.9, 95% CI: 0.3, 2.4) and those who work in the operating room or procedure areas (aOR=1.0, 95% CI: 0.6, 1.7) were not at increased risk of being seropositive. Differences in seropositivity based on HCW role were generally present in the unadjusted but not the adjusted models. Accounting for participation bias did not result in appreciably different estimated associations (Table 1).

Overall, 27% (95% CI: 25%-30%) of risk could be attributed to the occupational risk factors included in our analysis.

## Discussion

In this large serological testing effort of U.S. HCWs in a university-based healthcare system, we found an overall SARS-CoV-2 seroprevalence of 5.7% (95% CI: 5.2%-6.1%) following the initial surge of the epidemic. This rate is similar to the overall 6.0% seroprevalence reported from 13 U.S. academic medical centers and is generally consistent with other estimates among HCWs in the early stages of the pandemic. 1,2,7,11 The percent of seropositive HCWs who reported no past COVID-19-like illness was also similar to past estimates of asymptomatic SARS-CoV-2 infections.1

We found that community and demographic factors--contact with a confirmed or suspected COVID-positive case, Black race, and age under 30 years--were more important predictors of seropositivity than occupational factors. Notably, racial disparities, now well-documented in the general population, 12-14 persist in HCWs 1,15 after accounting for other risk factors, underscoring the fundamental racial inequities that have become the hallmark of the COVID-19 pandemic. We partially adjusted for community risk by including zip code-level COVID-19 incidence in our model, but we are unable to account for more proximal factors likely to have contributed to higher risk of infection among Black HCWs, including higher likelihood of exposure at home or use of public transportation. Ultimately, these risk factors are tied to entrenched, systemic social processes that underlie many individual and population health disparities. 16,17

We found few strong risk factors for infection in the workplace. Those who reported caring for patients with COVID-19 infection and those working in procedure areas where aerosol generating procedures are routine were not more likely to be seropositive, supporting the efficacy of PPE practices for known COVID-19 patients. Working in clinical areas was associated with increased odds of seropositivity, although the aORs were imprecise. We were not able to determine if risk associated with workplace location was from patient exposure, including unsuspected COVID-19 infected patients or from co-workers. Risk from contact with HCWs later found to be COVID-19-positive could reflect transmission during pre-symptomatic or asymptomatic infection prior to universal use of masks or transmission in settings where masks were not worn. In our hospital system, contact tracing identified staff eating together in break rooms as a risk factor for transmission.

While there appeared to be a difference in seropositivity based on workplace role in the unadjusted model, aORs were similar for all workplace roles in the adjusted model, underscoring the importance of considering demographic factors in assessing seroprevalence risk. Other studies that have not controlled for demographic or other risk factors outside of the workplace have reported different seroprevalence rates depending on job role 4,5; these results need to be interpreted with caution. Overall, we found that 27% (95% CI: 25%-30%) of risk could be attributed to the set of occupational risk factors we considered (Table 1). While the majority of risk may originate from the community, workplace risk cannot be ignored. In this cross-sectional analysis including time periods before and after universal masking, we were not able to assess the impact of infection prevention policies in reducing infection risk among HCWs. Going forward, studies should investigate the role of specific exposures contributing to infection risk, including risk from HCW to HCW exposure, and the efficacy of interventions to prevent transmission. It is critical to know if the current level of interventions, including screening for asymptomatic viral shedding, PPE practices and efforts to prevent HCW to HCW transmission, have substantially reduced or eliminated the workplace risk identified in this and other studies.

There are several limitations to this analysis. First, serological testing was voluntary, which may introduce bias if groups more likely to participate were also more (or less) likely to be seropositive. We partially adjusted for participation using demographic characteristics of HCWs overall. This adjustment at least partially accounts for poor representation of Black HCWs, in whom seroprevalence was higher than the population average, among those who volunteered for testing. However, other factors related to infection risk may have influenced participation. Second, a test with imperfect specificity in a population where seroprevalence is low will likely result in some false positives which may bias aORs towards the null. Third, we could not account for rapidly evolving infection prevention practices early in the pandemic and social behavior inside or outside the workplace. Lastly, the large influx of COVID-19 patients caused major shifts in care delivery and personnel deployment. Many HCWs worked in multiple locations and even different roles.

In conclusion, using a model incorporating demographic, community, and occupational risk factors for infection, we quantified community and occupational risk of SARS-CoV-2 seropositivity in HCWs. We found that most of this risk is due to community exposure; ongoing efforts to keep the healthcare workforce safe should emphasize risk mitigation in and outside the workplace. After adjusting for a number of community and occupational risk factors, race remains a critical marker of infection risk. Future seroprevalence studies of HCWs need to account for these community and demographic factors.

## Supporting information

Supplemental Material [[supplements/222877_file02.pdf]](pending:yes)

## Data Availability

The data that support the findings of this study are not publicly available due to participant privacy. Requests for data should be submitted to the authors and require approval from Emory Healthcare.

## Funding statement

Kristin N. Nelson and Benjamin A. Lopman received funding support from the Emory COVID-19 Response Collaborative. Julia M. Baker received funding support from the National Institute of Allergy and Infectious Diseases (T32AI074492).

## Acknowledgements

The authors are grateful to the Emory Healthcare employees and medical staff who volunteered for serology testing and completed the accompanying survey. We thank Laura Edison and Melissa Tobin-D’Angelo from the Georgia Department of Public Health for providing zip code level data on COVID-19 incidence and for their review of the manuscript. The authors have no conflicts of interest to disclose.

*   Received October 30, 2020.
*   Revision received October 30, 2020.
*   Accepted November 3, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/)

## References

1.  1.Self WH TM, Stubblefield WB, et al. Seroprevalence of SARS-CoV-2 Among Frontline Health Care Personnel in a Multistate Hospital Network — 13 Academic Medical Centers, April–June 2020. 31 August 2020 2020.
    
    

2.  2.Kluytmans-van den Bergh MFQ, Buiting AGM, Pas SD, et al. Prevalence and Clinical Presentation of Health Care Workers With Symptoms of Coronavirus Disease 2019 in 2 Dutch Hospitals During an Early Phase of the Pandemic. JAMA Network Open. 2020;3(5):e209673–e209673.
    
    

3.  3.Ran L, Chen X, Wang Y, Wu W, Zhang L, Tan X. Risk Factors of Healthcare Workers With Coronavirus Disease 2019: A Retrospective Cohort Study in a Designated Hospital of Wuhan in China. Clinical Infectious Diseases. 2020.
    
    

4.  4.Moscola J, Sembajwe G, Jarrett M, et al. Prevalence of SARS-CoV-2 Antibodies in Health Care Personnel in the New York City Area. JAMA. 2020;324(9):893–895.
    
    

5.  5.Iversen K, Bundgaard H, Hasselbalch RB, et al. Risk of COVID-19 in health-care workers in Denmark: an observational cohort study. The Lancet Infectious Diseases.
    
    

6.  6.Lai X, Wang M, Qin C, et al. Coronavirus Disease 2019 (COVID-2019) Infection Among Health Care Workers and Implications for Prevention Measures in a Tertiary Hospital in Wuhan, China. JAMA Network Open. 2020;3(5):e209666–e209666.
    
    

7.  7.Steensels D, Oris E, Coninx L, et al. Hospital-Wide SARS-CoV-2 Antibody Screening in 3056 Staff in a Tertiary Center in Belgium. JAMA. 2020;324(2):195–197.
    
    

8.  8.Leeds JS, Raviprakash V, Jacques T, Scanlon N, Cundall J, Leeds CM. Risk factors for detection of SARS-CoV-2 in healthcare workers during April 2020 in a UK hospital testing programme. EClinicalMedicine.
    
    

9.  9.Suthar MS, Zimmerman MG, Kauffman RC, et al. Rapid Generation of Neutralizing Antibody Responses in COVID-19 Patients. Cell Reports Medicine. 2020;1(3).
    
    

10. 10.Benichou J. A review of adjusted estimators of attributable risk. Stat Methods Med Res. 2001;10(3):195–216.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/096228020101000303&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11446148&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F03%2F2020.10.30.20222877.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000169505300003&link_type=ISI) 

11. 11.Barrett ES, Horton DB, Roy J, et al. Prevalence of SARS-CoV-2 infection in previously undiagnosed health care workers at the onset of the U.S. COVID-19 epidemic. medRxiv. 2020:2020.2004.2020.20072470.
    
    

12. 12.Moore JT RJ, Rose CE, Fuld J, Parise M, Kang GJ, Driscoll AK, Norris T, Wilson N, Rainisch G, Valverde E, Beresovsky V, Brune CA, Oussayef NL, Rose DA, Adams LE, Awel S, Villanueva J, Meaney-Delman D, Honein M. Disparities in Incidence of COVID- 19 Among Underrepresented Racial/Ethnic Groups in Counties Identified as Hotspots During June 5–18, 2020 — 22 States, February–June 2020. 2020.
    
    

13. 13.Adhikari S, Pantaleo NP, Feldman JM, Ogedegbe O, Thorpe L, Troxel AB. Assessment of Community-Level Disparities in Coronavirus Disease 2019 (COVID-19) Infections and Deaths in Large US Metropolitan Areas. JAMA Network Open. 2020;3(7):e2016938– e2016938.
    
    

14. 14.Webb Hooper M, Nápoles AM, Pérez-Stable EJ. COVID-19 and Racial/Ethnic Disparities. JAMA. 2020;323(24):2466–2467.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2020.8598&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32391864&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F03%2F2020.10.30.20222877.atom) 

15. 15.Nguyen LH, Drew DA, Graham MS, et al. Risk of COVID-19 among front-line health-care workers and the general community: a prospective cohort study. The Lancet Public Health.
    
    

16. 16.Arah OA. On the relationship between individual and population health. Med Health Care Philos. 2009;12(3):235–244.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11019-008-9173-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19107577&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F11%2F03%2F2020.10.30.20222877.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000272270300002&link_type=ISI) 

17. 17.Phelan JC, Link BG. Is Racism a Fundamental Cause of Inequalities in Health? Annual Review of Sociology. 2015;41(1):311–330.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev-soc-073014-112305&link_type=DOI)