key: cord-0904524-8w1wspjd
authors: Knudtzen, Fredrikke Christie; Jensen, Thøger Gorm; Lindvig, Susan Olaf; Rasmussen, Line Dahlerup; Madsen, Lone Wulff; Hoegh, Silje Vermedal; Bek-Thomsen, Malene; Laursen, Christian B.; Nielsen, Stig Lønberg; Johansen, Isik Somuncu
title: SARS-CoV-2 viral load as a predictor for disease severity in outpatients and hospitalised patients with COVID-19: A prospective cohort study
date: 2021-10-12
journal: PLoS One
DOI: 10.1371/journal.pone.0258421
sha: c6c15a5cb52c30843e6f8cfe260c98ccac5b0375
doc_id: 904524
cord_uid: 8w1wspjd

INTRODUCTION: We aimed to examine if severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) polymerase chain reaction (PCR) cycle quantification (C(q)) value, as a surrogate for SARS-CoV-2 viral load, could predict hospitalisation and disease severity in adult patients with coronavirus disease 2019 (COVID-19). METHODS: We performed a prospective cohort study of adult patients with PCR positive SARS-CoV-2 airway samples including all out-patients registered at the Department of Infectious Diseases, Odense University Hospital (OUH) March 9-March 17 2020, and all hospitalised patients at OUH March 10-April 21 2020. To identify associations between C(q)-values and a) hospital admission and b) a severe outcome, logistic regression analyses were used to compute odds ratios (OR) and 95% Confidence Intervals (CI), adjusting for confounding factors (aOR). RESULTS: We included 87 non-hospitalised and 82 hospitalised patients. The median baseline C(q)-value was 25.5 (interquartile range 22.3–29.0). We found a significant association between increasing C(q)-value and hospital-admission in univariate analysis (OR 1.11, 95% CI 1.04–1.19). However, this was due to an association between time from symptom onset to testing and C(q)-values, and no association was found in the adjusted analysis (aOR 1.08, 95% CI 0.94–1.23). In hospitalised patients, a significant association between lower C(q)-values and higher risk of severe disease was found (aOR 0.89, 95% CI 0.81–0.98), independent of timing of testing. CONCLUSIONS: SARS-CoV-2 PCR C(q)-values in outpatients correlated with time after symptom onset, but was not a predictor of hospitalisation. However, in hospitalised patients lower C(q)-values were associated with higher risk of severe disease.

We aimed to examine if severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) polymerase chain reaction (PCR) cycle quantification (C q ) value, as a surrogate for SARS-CoV-2 viral load, could predict hospitalisation and disease severity in adult patients with coronavirus disease 2019 (COVID-19).

We performed a prospective cohort study of adult patients with PCR positive SARS-CoV-2 airway samples including all out-patients registered at the Department of Infectious Diseases, Odense University Hospital (OUH) March 9-March 17 2020, and all hospitalised patients at OUH March 10-April 21 2020. To identify associations between C q -values and a) hospital admission and b) a severe outcome, logistic regression analyses were used to compute odds ratios (OR) and 95% Confidence Intervals (CI), adjusting for confounding factors (aOR).

We included 87 non-hospitalised and 82 hospitalised patients. The median baseline C qvalue was 25.5 (interquartile range 22.3-29.0). We found a significant association between increasing C q -value and hospital-admission in univariate analysis (OR 1.11, 95% CI 1.04-1. 19 ). However, this was due to an association between time from symptom onset to testing and C q -values, and no association was found in the adjusted analysis (aOR 1.08, 95% CI 

As the novel Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) sweeps through the world, detection of viral RNA by polymerase chain reaction (PCR) has become the gold standard for diagnosing coronavirus disease 2019 (COVID-19) [1, 2] . Nasopharyngeal or oropharyngeal swabs make up the majority of tests, since most patients are unable to produce sputum despite higher sensitivity of the latter [3] .

In viral diseases, the PCR Quantification Cycle-value or Cycle threshold-value (C q or C t -value) can be used as a surrogate for viral load, with inverse correlation between the C q -value and viral load. The use of C q -value as a prognostic marker for disease severity in viral respiratory infections has been tested with varying results [4] [5] [6] . For coronaviruses, there is evidence from cohort studies supporting a correlation between C q -values in upper airway samples and disease severity for human coronavirus in children, and between upper airway sample viral loads and disease severity for SARS-CoV-1 and Middle East Respiratory Syndrome (MERS) coronavirus infections in adults [7] [8] [9] .

The most common symptoms of COVID-19 are fever, cough and dyspnea, and the course of disease can be complicated with Acute Respiratory Distress Syndrome (ARDS), respiratory failure and death [10] [11] [12] . There is limited data on whether viral load of SARS-CoV-2 correlate with disease severity. Two Chinese studies found that both in-and outpatients with COVID-19 had lower C q -values indicating higher viral loads early in their disease course [13, 14] . A German study including hospitalised patients diagnosed with COVID-19 found that viral loads were high in the initial oropharyngeal samples and declining in 1-2 weeks [15] . Hospitalised patients in China with severe disease were found to have higher initial viral loads and prolonged time to reach PCR-negativity compared with patients with mild disease [16, 17] .

With this study, we aimed to examine if baseline PCR C q -values can identify 1) SARS-CoV-2 positive patients at increased risk of hospitalisation, and 2) hospitalised COVID-19 patients at increased risk of severe disease.

We hypothesized that the initial PCR C q -values were lower among hospitalised patients, as a surrogate for higher viral loads, compared with non-hospitalised patients. We also hypothesized that due to a failure to reduce viral burden after the initial infection phase, lower PCR C q -values were related to severe disease in hospitalised patients.

Odense University Hospital (OUH) serves as a tertiary hospital for the Region of Southern Denmark (approximately 1.2 million inhabitants) as well as a secondary hospital for the island of Funen (approximately 0.5 million inhabitants) [18] . The Danish public healthcare system supplies free, tax-funded healthcare for all residents.

Initially, the Danish national COVID-19 strategy was based on containment, where individuals who met the case definition were tested for SARS-CoV-2. The strategy later changed to mitigation, where only patients with symptoms of COVID-19 requiring hospital admission were tested for SARS-CoV-2.

We used the unique 10-digit personal identification number assigned to all individuals in Denmark at birth or upon immigration to link the following two registries electronically with laboratory data: 

For the hospital cohort, demography-, clinical-, laboratory-, management-and outcome data were gathered through review of medical records [11] . For the outpatient cohort, all eligible patients were invited to participate in an online survey two months after symptom onset. By signing an electronic consent form, a survey could be filled-out and electronically retrieved into a database. The data included information on demography, disease exposure, clinical symptoms of COVID-19, days until recovery and remaining symptoms (see S1 Appendix). Data on PCR assays, type of airway samples (naso-and/or oropharyngeal swab, or sputum), PCR C q -values and sample dates were collected from the Department of Clinical Microbiology, OUH and the Department of Clinical Microbiology, Lillebaelt Hospital.

SARS-CoV-2 detection was established on three different analysis platforms-the fully automated high throughput Cobas 6800 (Roche), the commercially available kit RealStar1 SARS--CoV-2 RT-PCR kit 1.0 (Altona Diagnostics) and a laboratory developed real-time (RT)-PCR.

On Cobas 6800 a 650 μl respiratory sample (oropharyngeal swab sample or sputum) was applied onto the system and subsequently RNA extraction, reverse transcription, PCR analysis and detection were performed. SARS-CoV-2 detection on Cobas 6800 included an internal RNA control, primers and probes targeting the ORF1a/b non-structural region that is unique for SARS-CoV-2 (target 1) and a conserved region in the structural protein envelope E gene that is shared by the Sarbecovirus subgenus (target 2).

RNA used for the RealStar1 SARS-CoV-2 RT-PCR test (Altona Diagnostics) and for the laboratory developed test was either extracted from: 1) 500 μl respiratory sample material (oropharyngeal swab sample or sputum) using MagNA Pure 96 (Roche) with the extraction kit DNA and viral NA large volume kit (Roche) using the protocol Pathogen Universal, or 2) 300 μl respiratory sample material (naso-and oropharyngeal swab sample or sputum) using the Maxwell1 16 Viral Total Nucleic Acid Purification Kit (Promega) following the manufacturer's protocol. RealStar1 SARS-CoV-2 RT-PCR kit 1.0 included three PCR analyses for the qualitative detection of and differentiation of Sarbecovirus subgenus (E gene) and SARS-CoV-2 specific RNA (S gene) in addition to an internal control. The kit was used according to the manufacturer's instructions with 30 μl reaction volume, and the 1-step RT-PCR was performed using Lightcycler 480 (Roche) or Stratagene Mx3005P (Agilent) in 96 well formats.

The laboratory developed real-time PCR E gene assay used for SARS-CoV-2 detection has been described previously [19] . This assay targeted a conserved sequence in the E gene region that is shared by the Sarbecovirus subgenus group ( For all assays, the PCR C q -value cut-off for a negative test was set at 40 cycles. As for choice of baseline test assay when more than one test was available, in assays with both a target for pan-Sarbecovirus and Coronavirus SARS-CoV-2 (Cobas 6800 and RealStar1 SARS-CoV-2 RT-PCR test), the C q -value for the pan-Sarbecovirus was chosen if available, if not the C q -value for Coronavirus SARS-CoV-2 target was used. If one sample was tested with multiple assays, the assays were prioritized in the following order; 1. Cobas 6800 (Roche), 2. The laboratory developed real-time PCR, and 3. RealStar SARS-CoV-2 RT-PCR test (Altona Diagnostics). The order was chosen by an experienced molecular biologist and a senior clinical microbiologist. For both naso/oropharyngeal swabs and sputum samples, the baseline PCRsample for each patient was set to the first registered test for that patient (= day 0). If there were multiple tests for one patient within the first 3 days (Day 0, 1 and 2), the sample with the lowest C q -value within this period was chosen.

We conducted a retrospective case-control study consisting of two different comparisons of sub-groups: 1) a case-control study with the hospital cohort as cases and the outpatient cohort as controls, and 2) a case-control study of our hospital cohort where the hospitalised patients with severe disease defined as ARDS, admittance to the Intensive Care Unit (ICU) and/or death during admission were included as cases, and the hospitalised patients with moderate disease (not fulfilling the definition of severe disease) were controls. The criteria for ARDS and grading of severity of ARDS were based on current international guidelines [20, 21] .

PCR C q -values were used to estimate predictors for 1) hospital admission, and 2) ARDS, ICU admission and/or death.

For baseline variables, descriptive statistics were reported as numbers and percentages for categorical variables and medians with interquartile ranges (IQR) for continuous variables. Chisquared test or Fisher's exact test were used to compare categorical variables between groups, student's t-test and Wilcoxon Mann-Whitney test were used for parametric and non-parametric continuous variables, respectively.

We plotted the C q -values according to days since symptom onset and examined a possible association using linear regression. To identify whether the C q -value could predict 1) hospital admission and 2) ARDS, ICU admission and/or death, we used logistic regression to compute odds ratios (OR) and 95% confidence intervals (CI). Analyses were adjusted for potential confounding variables, which based on the current knowledge on COVID-19 was predetermined to be age, sex, comorbidities, Body Mass Index (BMI) and days from symptom debut to baseline PCR-sample (1) . To reduce the risk of over-fitting, we only included confounders considered most important (sex, age) in the final multiple regression model (2).

Data on all patients were registered in a REDCap database hosted by Open Patient data Explorative Network (OPEN) [22] . STATA version 15 (StataCorp LP, Texas) was used for data processing and analyses.

This study was registered as a quality development project at OUH, approved by the Danish Data Protection Agency (j.nr. 20/16169 and 20/20759) and the Danish Patient Safety Authority (Sagsnr. 31-1521-344). All data were handled in accordance with The General Data Protection Regulation (GDPR), the Danish Act on Data Protection, the Danish Act on Research Ethics Review of Health Research Projects and the Danish Health Act. The study adheres to the STROBE guidelines for observational studies. All patients gave informed consent for study participation prior to inclusion.

A total of 169 patients were included in the final cohort; 87 from the outpatient cohort (Fig 1) .

The baseline patient characteristics of the two cohorts are shown in Table 1 .

The hospital cohort was significantly older (median age 63 years (IQR 55-74) vs 46 years (IQR 36-54), p<0.001), had higher weights (median BMI 26.5 (IQR 23.7-30.1) vs 24.6 (IQR 23.1-27.2), p = 0.003) and had a significantly higher proportion of all comorbidities except pulmonary diseases compared with the outpatient cohort.

Compared with the hospital cohort, the outpatients had a significantly higher degree of known exposure to COVID-19 (Table 1 ). In the outpatient group, 65 patients (77.4%) had travelled to a COVID-19 hot-spot in the 14 days prior to symptom onset. Of these, 58 patients (66.7%) had been on skiing holidays in the Tyrol region of Austria. COVID-19 symptoms in the two cohorts are illustrated in Fig 2. Compared with the outpatient cohort, hospitalised patients more often had fever, cough, dyspnoea and gastrointestinal symptoms but less often rhinitis/throat pain and loss of smell/ taste.

The median baseline SARS-CoV-2 PCR C q -value for the entire study population was 25.5 (IQR 22.3-29.0). The outpatients had a significantly lower median baseline SARS-CoV-2 PCR C q -value (24.6, IQR 21.8-27.5) compared with the hospitalised patients (median C q -value 26.9, IQR 23.6-31.3), p = 0.001 (Fig 3A) .

We found a statistically significant association between an increasing baseline C q -value and higher risk of admission to hospital (OR 1.11, 95%CI 1.04-1.19, p = 0.002) when using unadjusted logistic regression (see S1 Table) . However, this was mainly due to a strong association between time from symptom onset and C q -value (coefficient 0.26, 95%CI 0.15-0.38, p<0.001), as the patients in the outpatient cohort were tested significantly earlier in their course of disease compared with the hospital cohort (median 3 days (IQR 2-4) vs. 8 days (IQR 5-11), p<0.001). When adjusting for this difference in timing of testing, we no longer found a significant association between C q -values and admission (OR 1.00, 95%CI 0.91-1.09, p = 0.97), irrespective of further adjustment for confounding factors (OR 1.08, 95%CI 0.94-1.24 p = 0.27).

We found no significant difference in median baseline C q -values between naso-and/or-oropharyngeal swabs (143 patients; 13 naso-and-oropharyngeal and 130 oropharyngeal) and observed a significantly shorter time from symptom onset to first PCR sample in patients tested with naso-and/or-oropharyngeal swabs compared with patients tested with sputum samples (median 3 days (IQR 2-7) vs 8 days (IQR 6-11), p<0.001).

A total of 31 of the 82 patients (38.0%) in the hospital cohort developed severe COVID-19 disease. Patients with moderate and severe disease did not differ with regards to sex, age, BMI, comorbidities, tobacco or alcohol consumption (see Table 2 ). Patients with severe disease had significantly lower baseline C q -values compared with patients with moderate disease (median 24.8 (IQR 21.0-28.8) vs 28.1 (IQR 24.3-33.2), p = 0.01). We found a statistically significant association between lower C q -values and higher risk of severe disease (OR 0.89, 95%CI 0.81-0.98, p = 0.018). This association was independent of timing of the test in relation to symptom onset as well as presence of confounding factors including type of airway sample.

For patients with moderate disease, we found a direct linear association between the C qvalue and time of baseline test (Fig 3B) . In contrast, we observed that patients with severe disease had a low baseline PCR C q -value irrespective of time of testing. However, the regression coefficient between these two curves did not differ statistically (coef.-0.59 95%CI -1.20-0.02, p = 0.056) conferring to no significant interaction between C q -value and time of test. 

SARS-CoV-2 viral load in COVID-19

Median symptom duration in the out-patient cohort was 11 days (IQR 5-16) when excluding fatigue and loss of taste/smell, which persisted two months after onset of COVID-19 disease in 15 (17.2%) and 27 patients (31.0%), respectively. A SARS-CoV-2 PCR test was repeated in 17 patients after a median of 8 days (IQR 6-9); all C q -values increased, and 6 patients were PCR negative (Fig 4A) .

In the hospital cohort, the median time from hospital admittance to either discharge (n = 78) or death (n = 4) was 7.5 days (IQR 3-11). Multiple PCR-samples were available in 33 patients with moderate disease and 18 patients with severe disease and showed a more complex pattern compared with the out-patient cohort. We observed a less linear increase in C q -values, longer PCR positivity and several patients with subsequently decreasing C q -values (Fig 4B+4C) .

To our knowledge, this prospective study is the first to compare SARS-CoV-2 PCR C q -values between non-hospitalised and hospitalised patients. Our most important findings were the strong linear association between C q -values and time of testing after symptom onset, the correlation between lower C q -values and increased disease severity in hospitalised patients and the lack of association between C q -values and risk of hospitalisation. Our results of a linear association between C q -values and timing of the test after symptom onset are in line with available data that suggest higher SARS-CoV-2 viral loads in airway samples at symptom presentation followed by a gradual decrease [13] [14] [15] 17] . In this way, the novel Coronavirus differs from SARS-CoV-1, where viral loads were found to increase in airway samples until day 12-14 after symptom onset before decreasing [23, 24] .

In hospitalised patients, we found that a lower C q -value was associated with a significantly higher risk of severe disease irrespective of time of sampling and confounding factors. These findings are in line with the initial Chinese studies by Zheng and Liu, where patients with clinically severe disease had lower C q -values and were PCR positive longer than patients with mild disease [16, 17] . Due to the limited size of our population, a specific PCR C q -cutoff-value for patients in high risk of severe disease could not be estimated. Other studies are needed to explore this further in order to use it in a prediction model.

We could not confirm our hypothesis of an association between lower baseline C q -values and higher risk of hospital admission when adjusting for timing of the test and confounding factors. To our knowledge, there are no available studies that have investigated this possible correlation.

Two systematic literature reviews regarding the use of PCR C q -values in SARS-CoV-2 have been published since we undertook our study [25, 26] . In accordance with our findings, both studies report evidence of increasing C q -values in respiratory samples over time, and an association between C q -values and disease severity in hospitalized patients. However, the evidence was not conclusive and more data is needed in this area.

Of symptoms of COVID-19, we found significantly more patients in the hospitalized cohort with cough, dyspnea, fever and gastrointestinal symptoms. On the other hand, significantly more non-hospitalized patients suffered from rhinitis/throat pain and change in taste and/or smell. Other studies have shown conflicting results regarding change in taste/smell and severity of disease [27, 28] .

The main strength of this study is the well-described cohort with near-complete data of high quality for all patients as well as electronically retrieved C q -values. Furthermore, patients in both cohorts have been tested for SARS-CoV-2 based on the national standardized guidelines.

Our study has some limitations. Available data show varying inter-test agreement between different SARS-CoV-2 assays, especially in samples with high C q -values [29] [30] [31] [32] . In our study, three different PCR assays were used. This may have affected the reproducibility of the results. The different assays used reflect time and availability of assays during the pandemic; in the beginning most patients were tested using the in-house Flow, which was later replaced by the Cobas 6800. We also included results from both naso-and/or-oropharyngeal swabs and sputum samples, the latter only from hospitalised patients. Though we did not find any significant difference between sputum and oropharyngeal baseline C q -values, this could be explained by the time of sampling, as the sputum samples were generally tested later in the patients' disease course. However, when including type of airway sample in our regression model for hospitalised patients, it did not alter the results. All airway-samples used in our study were sampled after clinical indication, and not as part of a research project. The airway swabs have therefore been sampled by different medical personnel. As this is an operator-dependent procedure, this lack of standardization may have affected the results. Whereas data on the hospital cohort was based on hospital files, data from the outpatient cohort was based on questionnaires filled out approximately two months after onset of disease. Therefore, recall bias cannot be excluded. Due to the small size of the two cohorts, we cannot exclude a risk of type 2 errors.

Finally, despite omission of confounding variables deemed not statistically significant, we cannot exclude some degree of over-fitting of the multivariate regression analyses. More research in this area is needed, and larger cohorts would be able to confirm our findings with greater certainty.

There are still questions that need to be enlightened regarding why some patients get severe COVID-19 disease and others do not. Our findings suggest that clinicians cannot use the baseline C q -value in outpatients to predict risk of hospitalisation later in their disease course. However, treating physicians should be vigilant of admitted patients with initial low C q -values in their airway samples. When interpreting C q -values, time of symptom onset should be considered, and patients with continuously low C q -values should be closely monitored.

In conclusion, SARS-CoV-2 PCR C q -values correlated with time after onset of symptoms. Early in the disease course C q -values were low as a sign of high viral loads. We did not find C qvalues to be a predictor for hospitalisation. However, in hospitalised patients lower C q -values were found to be predictive of more severe disease. 

A new coronavirus associated with human respiratory disease in China

WHO. Laboratory testing for 2019 novel coronavirus (2019-CoV) in suspected human cases 2020

Clinical infectious diseases: an official publication of the Infectious Diseases Society of America

Influence of viral load in the outcome of hospitalized patients with influenza virus infection. European journal of clinical microbiology & infectious diseases: official publication of the European Society of Clinical Microbiology

Factors associated with real-time RT-PCR cycle threshold values among medically attended influenza episodes

Lower respiratory tract infection in the community: associations between viral aetiology and illness course

Journal of clinical virology: the official publication of the Pan American Society for Clinical Virology

Viral loads in clinical specimens and SARS manifestations

Viral Load Kinetics of MERS Coronavirus Infection. The New England journal of medicine

Predictors of infection, symptoms development, and mortality in people with SARS-CoV-2 living in retirement nursing homes

Low mortality of hospitalised patients with COVID-19 in a tertiary Danish hospital setting. International journal of infectious diseases: IJID: official publication of the International Society for Infectious Diseases

Clinical features, laboratory findings and predictors of death in hospitalized patients with COVID-19 in Sardinia, Italy

SARS-CoV-2 Viral Load in Upper Respiratory Specimens of Infected Patients. The New England journal of medicine

Clinical features and dynamics of viral load in imported and non-imported patients with COVID-19

Virological assessment of hospitalized patients with COVID-2019

Viral load dynamics and disease severity in patients infected with SARS-CoV-2 in Zhejiang province, China

Viral dynamics in mild and severe cases of COVID-19. The Lancet Infectious diseases

Statistikbanken Befolkning og valg Folketal 2020

Detection of 2019 novel coronavirus (2019-nCoV) by real-time RT-PCR. Euro surveillance: bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin

The Berlin definition of ARDS: an expanded rationale, justification, and supplementary material

Guidelines on the management of acute respiratory distress syndrome

Mandell, Douglas, and Bennett's Principles and Practice of Infectious Diseases

Viral shedding patterns of coronavirus in patients with probable severe acute respiratory syndrome

Laboratory diagnosis of SARS. Emerging infectious diseases

SARS-CoV-2 detection, viral load and infectivity over the course of an infection

A Systematic Review of the Clinical Utility of Cycle Threshold Values in the Context of COVID-19

Objective evaluation of anosmia and ageusia in COVID-19 patients: Single-center experience on 72 cases

Systemic inflammatory markers and psychophysical olfactory scores in coronavirus disease 2019 patients: is there any correlation?

Comparison of Cepheid Xpert Xpress and Abbott ID Now to Roche cobas for the Rapid Detection of SARS-CoV-2

Comparison of SARS-CoV-2 detection from nasopharyngeal swab samples by the Roche cobas 6800 SARS-CoV-2 test and a laboratory-developed real-time RT-PCR test

Comparison of seven commercial RT-PCR diagnostic kits for COVID-19

Comparison of Commercially Available and Laboratory-Developed Assays for In Vitro Detection of SARS-CoV-2 in Clinical Laboratories

For the laboratory-developed real time PCR, the primer and probe sequences targeting the internal control virus were kindly provided by Dr. Kurt Handberg, Department of Clinical Microbiology, Aarhus University Hospital. We thank Benedicte Christie Knudtzen for help with graphic design of Figs 1 and 2.