COVID-19 dynamics in an Ohio prison =================================== * Wasiur R. KhudaBukhsh * Sat Kartar Khalsa * Eben Kenah * Grzegorz A. Rempala * Joseph H. Tien ## Abstract **Importance** Incarcerated individuals are a vulnerable population for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Understanding SARS-CoV-2 dynamics in prisons is crucial for curbing transmission both within correctional facilities and in the surrounding community. **Objective** The purpose of this study was to identify transmission scenarios that could underly rapid, widespread SARS-CoV-2 infection among inmates in Marion Correctional Institution (MCI). **Design** Publicly available data reported by the Ohio Department of Rehabilitation and Corrections (ODRC) was analyzed using mathematical and statistical models. **Setting** We consider SARS-CoV-2 transmission dynamics among MCI inmates prior to and including April 16, 2020. **Participants** This study uses de-identified, publicly available SARS-CoV-2 test result data for MCI inmates. **Exposures** Inmates at MCI were considered exposed to potential infection with SARS-CoV-2. **Main outcome and measures** Results from mass testing conducted on April 16, 2020 were analyzed together with time of first reported SARS-CoV-2 infection among MCI inmates. **Results** Rapid, widespread infection of MCI inmates was reported, with nearly 80% of inmates infected within three weeks of first reported inmate case. These data are consistent with i) a basic reproduction number greater than 14, together with a single initially infected inmate, ii) an initial super-spreading event resulting in several hundred initially infected inmates, together with a basic reproduction number of approximately three, and iii) earlier undetected circulation of virus among inmates prior to April. **Conclusions and relevance** Mass testing data are consistent with extreme transmissibility, super-spreading events, or undetected circulation of virus among inmates. All scenarios consistent with these data attest to the vulnerabilities of prisoners to COVID-19. **Question** To identify transmission characteristics consistent with timing and extent of SARS-CoV-2 infection among inmates in Marion Correctional Institution. **Findings** Mathematical and statistical modeling finds three scenarios that are consistent with the observed widespread infection in Marion Correctional Institution: i) very high transmissibility corresponding to a basic reproduction number in the double digits, ii) an initial super-spreading event involving exposure of several hundred inmates, iii) undetected circulation of virus prior to the first documented case among inmates. **Meaning** High transmissibility, super-spreading events, and challenges with disease surveillance all attest to the vulnerabilities of prison populations to SARS-CoV-2. ## 1 Introduction Incarcerated individuals are a vulnerable population for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Factors contributing to SARS-CoV-2 transmission in prisons include shared housing, crowding, inability to social distance, challenges for hygiene, and more. Transmission can be extensive: as of November 6, 2020, various testing protocols across states and jurisdictions have identified 169,286 COVID-19 cases among prisoners and 38,438 cases among prison staff in the United States [1]. Outbreak sizes within facilities can be high: infections in more than 80% of prisoners at Marion Correctional Institution (MCI) in Ohio have been identified1, and similarly high levels of infection have been observed in correctional facilities under other jurisdictions [3]. The vulnerability of prisoners to COVID-19 and the potential for the penal system to serve as hot spots for transmission have been noted by many authors [4–8]. These outbreak sizes indicate widespread SARS-CoV-2 transmission. Perhaps even more remarkable is the apparent *intensity* of transmission within prisons. Here we analyze publicly available COVID-19 time series data from the Ohio Department of Rehabilitation and Corrections (ODRC). This analysis highlights the extreme potential for transmission in prisons, as well as the need for high quality, publicly available data on COVID-19 in correctional facilities. ## 2 Data and Methods ### Case data Mass RT-PCR testing at MCI of all inmates and partial testing of staff was conducted on April 16, 2020. Total number of inmates and number of inmates and staff at MCI testing positive for SARS-CoV-2 over time were obtained from public ODRC reports [2]. To account for lags in data entry, we accumulate the cases reported at MCI over April 16-23, 2020 as a single mass testing data point. ### Mathematical model We use a susceptible-exposed-infectious-recovered (SEIR) framework to model SARS-CoV-2 dynamics in MCI. Under the standard SEIR model, the proportions of individuals in the susceptible (*S**t*), exposed (*E**t*), infectious (*I**t*), and recovered (*R**t*) compartments (at time *t*) satisfy the following system of differential equations: ![Formula][1] where the positive parameters *β, α*, and *γ* denote the infection rate, the incubation rate, and the recovery rate. ### Statistical analysis We derive a likelihood function for observing *n* positives out of *N* incarcerated individuals on day *u* as follows. Using the Dynamic Survival Analysis approach of [9], we treat *S**t* as an improper survival function. Given observation up to time *T >* 0, the time *T**E* that an individual becomes infected and enters the *E* compartment follows the conditional probability density function ![Formula][2] where *τ**T* = 1 − *S**T*. The time *T**I* to becoming infectious has the conditional density ![Formula][3] and the recovery time *T**R* has the conditional density ![Formula][4] Note that the random variables *T**E*, *T**I* − *T**E*, and *T**R* − *T**E* are mutually independent, and that *T**I* − *T**E* and *T**R* − *T**E* are exponentially distributed with rates *α* and *γ*. Mass testing yields the number of individuals who test positive and the total number of tests administered on the day of mass testing. To use these data, let *T**N* denote the time when virus first becomes undetectable in an individual. We then describe the epidemic process by the pair of random variables (*T**E*, *T**N*). Let *ε* = *T**N* − *T**E*. Then, the conditional probability of an individual testing positive on the day of the mass testing (at time *u*) is given by ![Formula][5] We fix *ε* = 21 days, corresponding to detectable virus for three weeks following an individual becoming infectious [10, 11]. We set 1*/α* = 5.1*/* log(2) days (corresponding to a median incubation period of 5.1 days [12]), and assume a mean infectious period 1*/γ* of 5.6 days [13]. If *n* out of *N* individuals test positive on the day of the mass testing *u*, the likelihood function is given by ![Formula][6] with the probability of testing positive *p**u* as described in (2.5). ## 3 Results ### Reported outbreak time course According to ODRC reports [2], the first identified COVID-19 case at MCI was an infected staff member on March 29. Following this initial primary case, precautions such as cohorting and modified movement were enacted in order to restrict mixing and reduce transmission.2 The first COVID-19 case among inmates was identified on April 3. Mass RT-PCR testing of all inmates and partial testing of staff was conducted on April 16. By April 20, SARS-CoV-2 infection had been identified in 79% (1950/2453) of inmates and 35% (154/446) of staff.3 A time series of COVID-19 cases at MCI is shown in Figure 1. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/15/2021.01.14.21249782/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2021/01/15/2021.01.14.21249782/F1) Figure 1: Positive tests over time among (a) inmates and (b) staff at MCI, as reported in [2]. ### Basic reproduction number and initial exposure size To examine what values for the basic reproduction number *ℛ* are consistent with such rapid spread, we turn to the SEIR model (system (2.1)). In order for mass screening to identify 80% of the population as positive for COVID-19 on April 16, at least 80% of the population must have been infected by that date. Figure 2 shows the time needed to infect 80% of the population in the SEIR model as a function of the basic reproduction number (*ℛ*) and the initial number of exposed individuals (*E*). While *ℛ* values of two or larger are able to eventually infect 80% or more of the population, for modest *ℛ* values this can take on the order of months. Reproduction numbers larger than 14 are needed before outbreaks originating from a single exposed individual are able to generate infection levels in a three week span that are consistent with those reported for MCI. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/15/2021.01.14.21249782/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2021/01/15/2021.01.14.21249782/F2) Figure 2: Time to infect 80% of population in SEIR model with median incubation period of 5.1 days [12] and mean infectious period of 5.6 days. An alternative explanation is that the out-break involved more than one initially infected prisoner. From Figure 2, we see that for a fixed *ℛ* value, increasing *E* decreases the time needed to infect 80% of the population. However, an initial condition of *E* *>* 471 is needed for an outbreak with *ℛ* = 3.27 to infect 80% of the population within three weeks. ### Time of initial outbreak circulation A third possibility is that SARS-CoV-2 was circulating among prisoners prior to April 3. Figure 3 shows the log-likelihood (2.6) for observing the mass testing results in MCI according to the SEIR model (2.1) as a function of *ℛ* and the outbreak onset date, with *E* fixed at one. The outbreak onset date and *ℛ* are indistinguishable from the mass testing data alone, with the yellow “wishbone” in Figure 3 corresponding to outbreak onset and *ℛ* pairs that are equally likely given the data. Outbreak onsets in late March or later correspond to *ℛ* values larger than ten, while earlier outbreak onsets correspond to smaller *ℛ* values. Note that onset dates prior to March are required to give *ℛ* values of less than five. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/15/2021.01.14.21249782/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2021/01/15/2021.01.14.21249782/F3) Figure 3: *Left*. Log-likelihood in the *ℛ*-outbreak onset date plane, for mass testing data under a SEIR model with three week test-positive window following infectious onset. *Right. ℛ* values that maximize the likelihood as a function of outbreak onset date. This curve corresponds to following the right-most branch of the ‘wishbone’ in the likelihood plot. ## 4 Discussion The official reports from ODRC describe widespread infection of MCI inmates with SARS-CoV-2 within the span of three weeks. Our analysis indicates three (non-exclusive) possible explanations for this rapid spread: (i) values for the basic reproduction number that are far higher than the *ℛ* values between two and three that have been estimated for non-prison settings in the United States [13], (ii) initial exposure of a large number of infected prisoners, consistent with an extreme super-spreading event, and (iii) early, undetected circulation of SARS-CoV-2 among prisoners prior to April 3. All three possibilities speak to the vulnerabilities of the prison population for COVID-19. The permissive conditions for spread within correctional facilities, the challenges for disease surveillance and care in these settings, and the inextricable link between COVID-19 within correctional facilities and disease spread in the broader community, have been eloquently discussed by others [4–8]. None of these issues can be adequately addressed without collection and analysis of data regarding COVID-19 in prisons. Officially reported data such as those considered here can already be used for statistical and mathematical analyses [14, 15]. Further data collection and dissemination, for example distinguishing symptomatic versus mass testing and including both positive and negative test results, would facilitate additional analyses that would deepen our understanding of transmission. We urge health departments and corrections departments to undertake collection of accurate data and to make these data available for analysis. These efforts are imperative for preventing COVID-19 spread in prisons. ## Data Availability Code for this project is available at [https://github.com/wasiur/PrisonCOVID19Analysis](https://github.com/wasiur/PrisonCOVID19Analysis). ## Footnotes * 1 The July 21 ODRC report [2] reports that 2023 inmates had recovered from COVID-19. We take the May 5, 2020 listing of 2453 inmates at MCI as the denominator. * 2 As stated in the publicly available ODRC report from March 30, 2020: “Based on a staff member reporting a positive COVID test, MCI is operating under modified movement and the population is being separated by unit along with other precautionary measures. Every inmate at MCI is monitored daily and has their temperature taken along with a check for symptoms. Currently, there are no inmates symptomatic for COVID-19.” * 3 These numbers come directly from data on the ODRC website. There is a lag of a few days between when mass testing occurred (April 16) and when jumps in case counts are reported in the ODRC data (April 18 and April 19 for inmates; April 20 for staff). This lag may reflect delay in data entry. * Received January 14, 2021. * Revision received January 14, 2021. * Accepted January 15, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. [1].The Marshall Project. COVID tracking in prisons. [https://github.com/themarshallproject/COVID\_prison\_data](https://github.com/themarshallproject/COVID_prison_data), 2020. 2. [2].Ohio Department of Rehabilitation & Correction. COVID-19 inmate testing. [https://coronavirus.ohio.gov/static/reports/DRCCOVID-19Information.pdf](https://coronavirus.ohio.gov/static/reports/DRCCOVID-19Information.pdf), 2020. 3. [3]. Liesl M Hagan, Samantha P Williams, Anne C Spaulding, Robin L Toblin, Jessica Figlenski, Jeanne Ocampo, Tara Ross, Heidi Bauer, Justine Hutchinson, Kimberley D Lucas, Matthew Zahn, Chun Chiang, Timothy Collins, Alexis Burakoff, Juli Bettridge, Ginger Stringer, Ran- dolph Maul, Kristen Waters, Courtney Dewart, Jennifer Clayton, Sietske de Fijter, Radha Sadacharan, Linda Garcia, Naomi Lockett, Kirstin Short, Laxman Sunder, and Senad Handanagic. Mass testing for SARS-CoV-2 in 16 prisons and jails - six jurisdictions, United States, April-May 2020. MMWR. Morbidity and mortality weekly report, 69(33):1139–1143, 08 2020. 4. [4]. Matthew J. Akiyama, Anne C. Spaulding, and Josiah D. Rich. Flattening the curve for incarcerated populations — Covid-19 in jails and prisons. New England Journal of Medicine, 382(22):2075–2077, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMp2005687&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.14.21249782.atom) 5. [5]. Talha Burki. Prisons are “in no way equipped” to deal with COVID-19. Lancet (London, England), 395(10234):1411–1412, 05 2020. 6. [6]. Oluwadamilola T. Oladeru, Nguyen-Toan Tran, Tala Al-Rousan, Brie Williams, and Nickolas Zaller. A call to protect patients, correctional staff and healthcare professionals in jails and prisons during the COVID-19 pandemic. Health & Justice, 8(1):17, 2020. 7. [7]. R. Rubin. The challenge of preventing COVID-19 spread in correctional facilities. JAMA, 323(18):1760–1761. doi.10.1001/jama.2020.5427, 2020. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.14.21249782.atom) 8. [8]. Alysse G Wurcel, Emily Dauria, Nicholas Zaller, Ank Nijhawan, Curt Beckwith, Kathryn Nowotny, and Lauren Brinkley-Rubinstein. Spotlight on Jails: COVID-19 Mitigation Policies Needed Now. Clinical Infectious Diseases, 03 2020. ciaa346. 9. [9]. W. R. KhudaBukhsh, B. Choi, E. Kenah, and G. A. Rempala. Survival dynamical sys- tems: individual-level survival analysis from population-level epidemic models. Interface Focus, 10(1):20190048, 2020. 10. [10]. Xi He, Eric H. Y. Lau, Peng Wu, Xilong Deng, Jian Wang, Xinxin Hao, Yiu Chung Lau, Jessica Y. Wong, Yujuan Guan, Xinghua Tan, Xiaoneng Mo, Yanqing Chen, Baolin Liao, Weilie Chen, Fengyu Hu, Qing Zhang, Mingqiu Zhong, Yanrong Wu, Lingzhai Zhao, Fuchun Zhang, Benjamin J. Cowling, Fang Li, and Gabriel M. Leung. Temporal dynamics in viral shedding and transmissibility of covid-19. Nature Medicine, 26(5):672–675, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7326/M20-3012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.14.21249782.atom) 11. [11]. Ai Tang Xiao, Yi Xin Tong, and Sheng Zhang. Profile of RT-PCR for SARS-CoV-2: a preliminary study from 56 COVID-19 patients. Clinical Infectious Diseases, 04 2020. ciaa460. 12. [12]. S. A. Lauer, K. H. Grantz, Q. Bi, F. K. Jones, Q. Zheng, H. R. Meredith, A. S. Azman, N. G. Reich, and J. Lessler. The incubation period of coronavirus disease 2019 (COVID-19) from publicly reported confirmed cases: estimation and application. Annals of Internal Medicine, doi:10.7326/M20-0504, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7326/M20-0504&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32150748&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.14.21249782.atom) 13. [13]. S. M. Moghadas, A. Shoukat, M. C. Fitzpatrick, C. R. Wells, P. Sah, A. Pandey, J. D. Sachs, Z. Wang, L. A. Meyers, B. H. Singer, and A. P. Galvani. Projecting hospital utilization during the COVID-19 outbreaks in the United States. Proceedings of the National Academy of Sciences USA, [www.pnas.org/cgi/doi/10.1073/pnas.2004064117](http://www.pnas.org/cgi/doi/10.1073/pnas.2004064117), 2020. 14. [14]. M. C. Jimenez, T. L. Cowger, L. E. Simon, M. Behn, N. Cassarino, and M. T. Bassett. Epidemiology of COVID-19 among incarcerated individuals and staff in Massachusetts jails and prisons. JAMA Network Open, 3(8):e2018851. doi:10.1001/jamanetworkopen.2020.18851, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamanetworkopen.2020.18851&link_type=DOI) 15. [15]. Brendan Saloner, Kalind Parish, Julie A. Ward, Grace DiLaura, and Sharon Dolovich. COVID- 19 Cases and Deaths in Federal and State Prisons. JAMA, 324(6):602–603, 08 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2020.12528&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F15%2F2021.01.14.21249782.atom) [1]: /embed/graphic-1.gif [2]: /embed/graphic-2.gif [3]: /embed/graphic-3.gif [4]: /embed/graphic-4.gif [5]: /embed/graphic-5.gif [6]: /embed/graphic-6.gif