key: cord-0718050-jb8lhsc1 authors: Shrestha, Nabin K; Marco Canosa, Francisco; Nowacki, Amy S; Procop, Gary W; Vogel, Sherilynn; Fraser, Thomas G; Erzurum, Serpil C; Terpeluk, Paul; Gordon, Steven M title: Distribution of Transmission Potential during Non-Severe COVID-19 Illness date: 2020-06-29 journal: Clin Infect Dis DOI: 10.1093/cid/ciaa886 sha: 5fbec5482df364bab6b604e29fca26b7f7d4f42d doc_id: 718050 cord_uid: jb8lhsc1 BACKGROUND: Patients recovering from coronavirus disease 2019 (COVID-19) often continue to test positive for the causative virus by PCR even after clinical recovery, thereby complicating return-to-work plans. The purpose of this study was to evaluate transmission potential of COVID-19 by examining viral load with respect to time. METHODS: Health care personnel (HCP) at Cleveland Clinic diagnosed with COVID-19, who recovered without needing hospitalization, were identified. Threshold cycle (Ct) for positive PCR tests were obtained and viral loads calculated. Association of viral load with days since symptom onset was examined in a multivariable regression model, which was reduced by stepwise backward selection to only keep variables significant at a level of 0.05. Viral loads by day since symptom onset were predicted using the model and transmission potential evaluated by examination of a viral load-time curve. RESULTS: Over six weeks, 230 HCP had 528 tests performed. Viral loads declined by orders of magnitude within a few days of symptom onset. The only variable significantly associated with viral load was time since onset of symptoms. Of the area under the curve (AUC) spanning symptom onset to 30 days, 96.9% lay within the first 7 days, and 99.7% within 10 days. Findings were very similar when validated using split-sample and 10-fold cross-validation. CONCLUSIONS: Among patients with non-severe COVID-19, viral loads in upper respiratory specimens peak by two or three days from symptom onset and decrease rapidly thereafter. The vast majority of the viral load-time AUC lies within 10 days of symptom onset. A large number of health care personnel (HCP) have contracted coronavirus disease 2019 [1] , the disease caused by the Severe Acute Respiratory Syndrome (SARS)associated Coronavirus-2 (SARS-CoV-2). When infected HCP recover, it becomes necessary to allow them back to work without placing their patients or coworkers at risk of acquiring the infection from them. In our institution, we adopted a policy which required at least 7 days since onset of symptoms, at least 3 days of improvement in respiratory symptoms and fever, and 2 negative SARS-CoV-2 PCR tests in a 24-hour period, an approach consistent with CDC guidance for returning essential personnel to work at the time. When HCP began to get repeat testing, however, we discovered that many continued to have positive PCR tests. Recognizing the reality of testing challenges in many parts of the country, the CDC updated its guidance for returning essential personnel to work, with an option of adopting either a testbased strategy requiring two negative tests 24 hours apart, or alternatively a symptom-based or time-based strategy that would allow return to work if personnel met specified clinical criteria including at least 3 days since improvement in symptoms and at least 10 days since onset of symptoms, or at least 10 days since a positive test in the absence of symptoms [2] . A non test-based strategy conveniently allows HCP back to work without having to contend with persistently positive test results. However, it would be more reassuring for healthcare organizations to embrace a non test-based strategy if data showed that such an approach is unlikely to pose a risk of infection to others from returning HCP. The risk of transmission of an infectious agent is not the same throughout the course of illness. Transmission risk at any instant would be expected to depend on the quantity of the infectious agent in relevant clinical sites at the time. Thus, for a respiratory viral infection, the risk of transmission at any instant depends on the viral load in respiratory sites at that instant. Hence, examination of viral load-time curves should provide good insight into transmission potential over different intervals during an episode of illness. A c c e p t e d M a n u s c r i p t The purpose of this study was to evaluate transmission potential by examining viral load with respect to time since onset of symptoms in HCP infected with SARS-CoV-2. A c c e p t e d M a n u s c r i p t This was a cohort study examining viral loads in HCP diagnosed with COVID-19. The HCP cohort was chosen for the study because of documentation of the date of onset of illness in templated notes of monitoring calls from Occupational Health. The study was approved by the Cleveland Clinic Institutional Review Board (IRB no. 20-1631). A waiver of informed consent and waiver of HIPAA authorization were approved to allow access to PHI by the research team, with the understanding that sharing or releasing identifiable data to anyone other than the study team was not permitted without additional IRB approval. The Cleveland Clinic began testing for SARS-CoV-2 by PCR on March 12, 2020, with a separate streamlined process to test HCP started at a similar time. All HCP tested between March 16 and April 20, 2020 were identified from the laboratory information system and screened for inclusion in the study. All those with at least one positive test were included. Those with severe illness (who required hospitalization during the course of their illness), and those for whom a symptom onset date could not be ascertained, were excluded. Introduction of the SARS-CoV-2 PCR test was done in concert with a standardized work flow for collection of specimens. Because of overwhelming demand, and consequently unacceptable delays, a separate queue was created for HCP. Specimen collection consisted of obtaining a nasopharyngeal swab by trained personnel at a single drive-through location and transporting it to the laboratory in universal/viral transport medium. A c c e p t e d M a n u s c r i p t To maintain consistency, a decision was made to test all HCP specimens on the same testing platform. The reverse-transcription (RT) -PCR test developed at the CDC was used, and internal validation studies were performed in-house before introduction of the test. Using commercially-available plasmids that contained the nucleocapsid (N) gene, the limit of detection was found to be 20 copies/µL for upper respiratory specimens and 2 copies/µL for lower respiratory specimens. The threshold cycle (Ct) for a real time PCR test is the cycle at which products of amplification first register at or above the detection threshold. The higher the load of the target in the specimen, the lower the Ct. The minimum detectable viral load (MDVL) is the smallest viral load that can be detected by the test. For this study, the Ct for a specimen containing the MDVL (CtMDVL) is 40 cycles, as any specimen registering a Ct higher than 40 is considered a negative test. Thus CtMDVL -Ct is zero for a specimen containing the MDVL, and is proportionally higher for any sample with a higher viral load. A c c e p t e d M a n u s c r i p t Each unit difference in CtMDVL -Ct represents a difference of one PCR cycle. Given that each PCR cycle results in a doubling of the amplified target, the quantity of the viral target as a multiple of the MDVL can be determined by the formula: 2^(CtMDVL -Ct). If the efficiency of the PCR reaction is less than 100%, calculated numbers will be proportionally smaller, but relative differences between any two viral loads will still be the same. Thus 2^(CtMDVL -Ct) is the actual viral load in the specimen, expressed as a multiple of the MDVL. Ct values for each PCR component (N1, N2, and N3) of the positive tests were obtained. The mean of the three Ct values was considered the Ct for the test. Relevant demographic information, comorbidity data, and date of onset of symptoms, were obtained from the electronic medical record. In addition to time since onset of illness, factors that were considered to possibly influence viral load were age, gender, race, body mass index (BMI), and comorbid conditions such as chronic lung disease, smoking status, immunosuppressed status, chronic heart disease, chronic kidney disease, and chronic liver disease. Chronic heart disease was defined as coronary artery disease or valvular heart disease with moderate to severe stenosis or regurgitation. Chronic lung disease was defined as chronic obstructive pulmonary disease, bronchiectasis, asthma, or chronic interstitial lung disease. Chronic kidney disease was defined as an estimated glomerular filtration rate of 60 mL/min or less, or maintenance hemodialysis or peritoneal dialysis. Chronic liver disease was defined as cirrhosis from any cause. Immunocompromised state was defined as being on biologic agent therapy or other immunosuppressive medications, having a hematological malignancy, or a history of hematological or solid organ transplant. M a n u s c r i p t Analyses were performed by N.K.S. and A.S.N. using R version 4.0.0 [3] . BMI data were missing for 15 study subjects (6.5%). Examination of missing values with margin plots from the R package VIM [4] , showed that the missing BMI values were missing completely at random. Values for these missing ones were imputed using a method of multivariate imputation by chained equations using the R package mice [5] . Scatter plots of viral load over time were investigated. The association of viral load with time since onset of symptoms was explored using restricted cubic spline regression with the R package splines [3] . Splines allow for the modeling of nonlinear relationships among variables. On 10-fold crossvalidation, a restricted cubic spline regression model with 4 degrees of freedom yielded the lowest error rate. This also produced a viral load-time curve that appeared to be biologically plausible (i.e. viral load increasing at the time of onset of symptoms). All numeric variables, diabetes mellitus, smoking status, and other categorical variables where the smaller group contained at least 10% of subjects were considered as covariates. Variables that did not remain significant at a level of 0.05 were dropped by backward elimination in a stepwise fashion. Viral loads on successive days since onset of symptoms were predicted using the regression model. Areas under the viral load-time curve were determined by integration of the function representing the relationship between the predicted viral load and days since onset of symptoms. Model performance was validated using both split-sample (80:20) and 10-fold crossvalidation. Graphics were created using the R base and ggplot2 packages [3, 6] . A c c e p t e d M a n u s c r i p t During the 6-week study period, 252 HCP tested positive for SARS-CoV-2. Twenty (8%) were excluded because they were hospitalized during the course of their illness and 2 (<1%) were excluded because the date of onset of symptoms could not be ascertained on review of their records. The remaining 230 were included in the study, and had a total of 528 tests performed. Of these patients, 106 (46%) had a single test done, 28 (12%) had two, 40 (17%) had three, 39 (17%) had four, 12 (5%) had five, and 5 (2%) had six tests done. Days from onset of symptoms to first test was distributed with a mean (SD) of 4.9 (4.9) days. Characteristics of the included HCPs are shown in Table 1 . Figure 1 ). Fluctuations in viral load over time seem to occur once viral loads are at lower levels, but these viral loads generally remain low. In a restricted cubic spline regression model the only variable that remained significantly associated with viral load was time since onset of symptoms. The performance of the final model is shown in Figure 2 . A c c e p t e d M a n u s c r i p t A c c e p t e d M a n u s c r i p t This study shows that in otherwise healthy subjects with non-severe illness from COVID-19, the SARS-CoV-2 viral load is very high within 2-3 days of onset of symptoms and falls rapidly by orders of magnitude within a few days. It should be noted that none of these HCP were treated with hydroxychloroquine or other COVID-19-related treatments. The premise of this study was that viral load in the nasopharynx should be a reasonable marker of viral loads at other locations in the respiratory tract. This is a reasonable assumption. Viral loads in saliva and endotracheal aspirates decreased proportionally over time in one study [7] ; and viral loads in tongue, nasal, and mid-turbinate swabs were linearly correlated with viral loads in nasopharyngeal swabs in another [8] . Modeling suggests that the contagiousness of respiratory viral infection is proportional to the quantity of microorganism at respiratory sites [9] . Such an association of viral load with transmission risk has been shown in other viral infections such as HIV and HHV6 [10, 11] . Estimates from our study imply that, in individuals with non-severe illness, transmission potential of COVID-19 is greatly diminished by 7-10 days since the onset of symptoms. The distribution of transmission potential found in this study is consistent with related findings in other studies. A rapid decline in viral load over a few days has been demonstrated in several studies [7, [12] [13] [14] [15] [16] . A positive PCR test does not necessarily correlate with presence of transmissible virus. In prior studies, live virus could not be isolated after 8 days since symptom onset even with positive PCR results [13, 17] , and virus could not be isolated from specimens with Ct values ≥24 [17] , supporting the hypothesis that transmission risk is low when viral load is low. A study that estimated infectiousness based on transmission events, not on viral loads, concluded that infectiousness declined quickly within 7 days [15] . A contact tracing assessment found no infections among those exposed to the index cases after more than five days since symptom onset [18] . The main strength of this study is that it demonstrates a way to quantitate transmission potential of COVID-19 over the course of illness. It is very difficult to determine the absolute risk of transmission A c c e p t e d M a n u s c r i p t of any viral infection at any specified viral load. However, knowledge of relative viral loads over time can provide a very good estimate of relative transmission potential over different intervals during the course of illness. Our study methodology using Ct values from PCR tests is simple enough that any healthcare organization can replicate the study internally using local data to help guide decisions in their own organizations. A limitation of our study is that specimens were not collected in all subjects at pre-specified time points. Once two negative test were obtained, the subject would not have undergone subsequent tests. Thus data points in the study could be expected to have been increasingly biased towards higher viral loads as time from onset of symptoms increased. Thus the further out, the more likely that the viral load would be an over-estimate of the true viral load, thereby explaining the upward tick at the tail end of the curve in figure 2 . However, this limitation has minimal effect on the overall conclusion that the majority of transmission potential lies early in the course of illness. The extremely high viral loads within two or three days since onset of symptoms suggest that viral loads may be almost as high in the immediate pre-symptomatic period, suggesting substantial transmission potential in the pre-symptomatic period. Since our study did not have actual viral loads in the pre-symptomatic period for any patient, concerns about pre-symptomatic transmission based on our study alone are hypothetical. However, pre-symptomatic transmission has been documented elsewhere [15, 19] . A study that examined 77 infector-infectee transmission pairs estimated that 44% of secondary cases were infected during the index cases' pre-symptomatic stage [15] . The potential for substantial transmission in the pre-symptomatic stage supports the use of masks in public during the pandemic in the current situation where there is active community transmission, the majority of the population remains susceptible, and there is no available vaccine to provide protection. Although this study examined HCP, findings can be expected to be similar in non-HCP with nonsevere illness, thereby allowing non-healthcare organizations to make decisions on appropriate timing for returning recovering employees to work based on this study. Findings from this study should not, however, be extrapolated to individuals who have severe illness, as it is possible such patients may have A c c e p t e d M a n u s c r i p t more prolonged shedding of virus [20] . Similarly these findings cannot be used to make decisions about when isolation can be discontinued in hospitalized patients. Testing to return personnel to work after COVID-19 poses challenges. PCR tests may continue to remain positive long after clinical recovery. Our study's finding that transmission potential is largely over by 10 days from onset of symptoms, supports a policy that would allow personnel with non-severe COVID-19 who have recovered clinically, to return to work after 10 days since onset of symptoms, without a need for laboratory testing. There will undoubtedly be outliers who shed virus for longer periods of time, but it is likely that the viral load in such situations is low, and it is not clear if this represents transmissible virus [13, 17] . A requirement of wearing masks for a period of time after returning to work would protect against such outliers. In our institution, our return-to-work policy for HCP infected with COVID-19 was updated to require at least 10 days since onset of symptoms, at least 3 days since clinical recovery, and a requirement to wear a mask for 2 weeks after return, without a need for testing to return to work. In conclusion, among patients with non-severe COVID-19, viral loads in upper respiratory specimens peak by two or three days from onset of symptoms and decrease greatly within a few days. The M a n u s c r i p t A c c e p t e d M a n u s c r i p t A c c e p t e d M a n u s c r i p t A c c e p t e d M a n u s c r i p t Figure 3 Characteristics of Health Care Personnel with COVID-19 -United States ?CDC_AA_refVal=https%3A%2F%2Fwww.cdc.gov%2Fcoronavirus%2F2019-ncov%2Fhealthcare-facilities%2Fhcp-return-work.html. Accessed 1 R: A Language and Environment for Statistical Computing Imputation with the {R} Package {VIM} Multivariate Imputation by Chained Equations in R Elegant Graphics for Data Analysis Temporal profiles of viral load in posterior oropharyngeal saliva samples and serum antibody responses during infection by SARS-CoV-2: an observational cohort study Swabs Collected by Patients or Health Care Workers for SARS-CoV-2 Testing Multiscale model within-host and betweenhost for viral infectious diseases Viral load and heterosexual transmission of human immunodeficiency virus type 1 Estimating the Risk of Human Herpesvirus 6 and Cytomegalovirus Transmission to Ugandan Infants from Viral Shedding in Saliva by Household Contacts Presymptomatic SARS-CoV-2 Infections and Transmission in a Skilled Nursing Facility Virological assessment of hospitalized patients with COVID-2019 Natural History of Asymptomatic SARS-CoV-2 Infection Temporal dynamics in viral shedding and transmissibility of COVID-19 SARS-CoV-2 viral load in upper respiratory specimens of infected patients Predicting infectious SARS-CoV-2 from diagnostic samples Contact Tracing Assessment of COVID-19 Transmission Dynamics in Taiwan and Risk at Different Exposure Periods Before and After Symptom Onset Presymptomatic Transmission of SARS-CoV-2 Clinical course and risk factors for mortality of adult inpatients with A c c e p t e d M a n u s c r i p t