744
Views
0
CrossRef citations to date
0
Altmetric
Original Article

SARS-CoV-2 RT-qPCR Ct values in saliva and nasopharyngeal swab samples for disease severity prediction

ORCID Icon, ORCID Icon, ORCID Icon, ORCID Icon, ORCID Icon, ORCID Icon, ORCID Icon, ORCID Icon, ORCID Icon, ORCID Icon, ORCID Icon, ORCID Icon & ORCID Icon show all
Article: 2213106 | Received 02 Sep 2022, Accepted 08 May 2023, Published online: 17 May 2023

ABSTRACT

Background

Comparison of clinical value of RT-qPCR-based SARS-CoV-2 tests performed on saliva samples (SSs) and nasopharyngeal swab samples (NPSs) for prediction of the COVID-19 disease severity.

Methods

Three paired SSs and NPSs collected every 3 days from 100 hospitalised COVID-19 patients during 2020 Jul-2021 Jan were tested by RT-qPCR for the original SARS-CoV-2 virus and compared to 150 healthy controls. Cases were divided into mild+moderate (Cohort I, N = 47) and severe disease (Cohort II, N = 53) cohorts and compared.

Results

SARS-CoV-2 was detected in 65% (91/140) vs. 53% (82/156) of NPSs and 49% (68/139) vs. 48% (75/157) of SSs collected from Cohort I and II, respectively, resulting in the total respective detection rates of 58% (173/296) vs. 48% (143/296) (P = 0.017). Ct values of SSs were lower than those of NPSs (mean Ct = 28.01 vs. 30.07, P = 0.002). Although Ct values of the first SSs were significantly lower in Cohort I than in Cohort II (P = 0.04), it became negative earlier (mean 11.7 vs. 14.8 days, P = 0.005). Multivariate Cox proportional hazards regression analysis showed that Ct value ≤30 from SSs was the independent predictor for severe COVID-19 (HR = 10.06, 95% CI: 1.84–55.14, P = 0.008).

Conclusion

Salivary RT-qPCR testing is suitable for SARS-CoV-2 infection control, while simple measurement of Ct values can assist in prediction of COVID-19 severity.

Introduction

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a novel coronavirus that emerged in late 2019 and spread across the globe, with the number of newly confirmed cases growing to more than 760 million [Citation1]. Nasopharyngeal specimen (NPS) quantitative reverse transcription PCR (RT-qPCR) is still considered the gold-standard laboratory technique for the etiological diagnosis of SARS-CoV-2 infection, having the highest sensitivity of all upper respiratory tract specimens [Citation2,Citation3]. However, collecting NPSs is technically challenging, and there is a risk of droplet transmission of the infection to healthcare professionals [Citation4]. In addition, the number of false-negative results depends on the quality of the sample, sensitivity and specificity of the used reagents and on the experience and training of the person performing the procedure as well [Citation5]. Multiple studies support non-invasive specimen collection from saliva as a suitable alternative to NPSs [Citation6–11]. Saliva specimens are already used in some countries for pandemic control as a new gold-standard sample type [Citation12]. However, the clinical application of the cycle threshold (Ct) values of saliva and NPS samples is understudied, and definitive data to support the prediction of disease severity is lacking [Citation13]. Quantitative changes of viral SARS-CoV-2 RNA at different stages of disease could provide clinically significant information regarding the severity of the disease [Citation14]. Hence, viral load changes might serve as a guideline for therapy prescription and impact clinical outcomes of the disease [Citation15].

In this study, we compared RT-qPCR measures of paired and serial SARS-CoV-2 nasopharyngeal and saliva swab samples to assess the prognostic test value and on predicting the COVID-19 disease severity.

Materials and methods

Participant cohorts

NPSs and corresponding SSs were collected from 100 COVID-19 patients treated at the Center of Infectious Diseases of Vilnius University Hospital Santaros Klinikos (VUHSK, Vilnius, Lithuania) from 2020 July to 2021 January for the original SARS-CoV-2 testing. For the control group, NPSs and corresponding SSs were collected from 150 persons visiting VUHSK for routine procedures (drug prescriptions, vaccines, etc.). Criteria for inclusion of the COVID-19 patient’s group were as follows: persons over the age of 18; a diagnosis of COVID-19 has been established and are being treated in an inpatient facility. Criteria for inclusion of the control group in the study were as follows: persons over the age of 18; not diagnosed with COVID-19. Exclusion criteria for COVID-19 patients and control groups were as follows: pregnant women; persons under the age of 18; persons with mental disabilities; persons suffering from oncological diseases.

This was a prospective, observational study. The primary outcome was the ability of salivary SARS-CoV-2 RNA RT-qPCR testing to predict the severity of COVID-19. The primary predictive variable was salivary SARS-CoV-2 RNA RT-qPCR test Ct values, dichotomized into Ct ≤ 30 and Ct > 30. The secondary outcomes were suitability of salivary SARS-CoV-2 RNA RT-qPCR test to detect severe COVID-19, and time till SARS-CoV-2 virus clearance, defined as duration of positive salivary or nasopharyngeal SARS-CoV-2 RNA RT-qPCR test.

The bioethics approval was obtained by the Regional Bioethics Committee (No. 2020/7-N4˗1245˗725), and the study was conducted following the Declaration of Helsinki. The informed consent forms were gained from all the participants before the beginning of the study. The SARS-CoV-2 infection was confirmed by RT-qPCR from NPS before the hospitalization at VUHSK. According to the disease severity, patients were dichotomized into two cohorts: mild and moderate disease cases (N = 47, Cohort I) and severe disease cases (N = 53, Cohort II). The patients were divided into groups dependent on the disease severity according to the National Institutes of Health guidelines [Citation16], which briefly define the severe cases as individuals having an oxygen saturation measured by pulse oximetry on room air at sea level (SpO2) <94%, a ratio of arterial partial pressure of oxygen to fraction of inspired oxygen (PaO2/FiO2) <300 mm Hg, and a respiratory rate >30 breaths/min, or lung infiltrates >50%; moderate – SpO2≥94% and evidence of lower respiratory disease during clinical assessment or imaging; mild – any of the various signs and symptoms of COVID-19 (e.g. fever, cough, sore throat, malaise, headache, muscle pain, nausea, vomiting, diarrhoea, loss of taste and smell) but do not have shortness of breath, dyspnoea, or abnormal chest imaging.

Sample collection

NPSs were collected by following manufacturers’ protocols and hospital’s procedures. At least 2 mL of saliva samples from 50 patients were collected into sterile tubes, whereas 2 mL saliva samples from another 50 patients and 150 controls were taken by using a Saliva Collection Kit (#SCK, Diagnolita, Vilnius, Lithuania) and mixed with 2 mL stabilizing media immediately after the collection as described in recent publication [Citation11]. All the samples were stored at −20°C until use.

When possible, two subsequent SS and NPS collections from the same hospitalized patients were performed 3 days apart. Altogether, 296 infectious NPSs (I 100, II 97, III 99) and 296 SSs (I 99, II 98, III 99) were taken, in total assembling a collection of 892 samples (592 infectious, 300 healthy).

RNA extraction

Total RNA was extracted by using Saliva NA Purification Kit (#SPK, Diagnolita) and following the manufacturer’s protocol. Briefly, the purification solution, consisting of stabilization media, binding solution, magnetic beads, reducing reagent, and proteinase K, was mixed with 300 µL of nasopharyngeal swab medium or 250 µL of unstabilized saliva (490 µL of the saliva-stabilizer mix for stabilized saliva samples) and incubated in thermoshaker for 10 min at 65°C temperature and 900 rpm [Citation11]. After two washing steps with PB1 and PB2 buffers (Diagnolita), the residual ethanol was eliminated by heating. Purified nucleic acids were eluted into 50 µL elution buffer. Eluted nucleic acids were immediately used for downstream applications or stored at −20 until use.

cDNA synthesis and quantitative reverse transcription PCR (RT-qPCR)

Five microlitre of total RNA was used for one-step cDNA synthesis and RT-qPCR reaction in a total volume of 20 µL consisting of 20× Primer-Probe Mix for SARS-CoV-2 Detection (#PPM, Diagnolita) [Citation11], 4× TaqMan Fast Virus 1-Step Master Mix (Thermo Fisher Scientific (TFS), Waltham, MA, USA), and RNase-free water. For cDNA synthesis, reaction mixes were incubated at 50°C for 5 min, and the subsequent RT-qPCR runs consisted of enzyme activation at 95°C for 20 s and 40 cycles of 95°C for 3 s and 60°C for 30 s. RT-qPCR reactions were run on the QuantStudio 5 Real-Time PCR System by using QuantStudio Design & Analysis Software v1.4.1 (both from Applied Biosystems, TFS).

Statistical analysis

Continuous variables are presented as means and standard deviation (SD) or as median and interquartile range (IQR). Categorical variables are reported as numbers and percentages (%). Differences in characteristics among the illness severity groups were tested using Student‘s t-test for normally distributed data and using Mann–Whitney U-test for non-normally distributed data. Paired Student’s t-test was used to compare characteristics of two related groups. For categorical variables, the Pearson‘s chi-squared test or Fisher’s exact test was used. Pearson linear or Spearman rank (for ordinal variables and as a more robust measure for interval variables with highly influential points) correlation coefficients were used to find the associations between clinical, laboratory characteristics and SARS-CoV-2 RT-qPCR Cycle threshold (Ct) values. FDR correction was applied for the analyses that involved multiple comparisons to adjust the resulting P-values. Cox regression analysis was used to determine which factors had a significant impact for prognosis of severe COVID-19 infection. The Kaplan–Meier method (log-rank testing) was used to predict the positivity for SARS-CoV-2. A P-value ≤ 0.05 was defined as statistically significant. Data were analysed using the Statistical Package for Social Science software version 21.0 for Windows (IBM Corp., Armonk, NY, USA). Graphs were depicted by using GraphPad Prism v8.0.1 (GraphPad Software, Inc., San Diego, CA, USA).

Results

Demographics and clinical profiles

Demographic and clinical characteristics of patients according to the disease severity (Cohort I vs. Cohort II) are presented in . The mean age was 56.7 ± 13.2 years (range: 25–90 years), 52.0% were male. Cohort II patients were significantly older (mean: 61.6 vs. 51.2 years, P < 0.001), had higher National Early Warning (NEW) score on admission (median: 3 vs. 1, P < 0.001) and bigger Charlson comorbidity index (median: 2 vs. 1, P < 0.001). Treatment with dexamethasone (20 vs. 8 patients, P = 0.020); dexamethasone and remdesivir combination (28 vs. 12 patients, P=0.043) or antibacterial treatment (47 vs. 22 patients, P < 0.001) was applied more often for Cohort II patients as well. Three Cohort II patients were treated with a high flow nasal cannula, one patient needed intubation and invasive mechanical ventilation and three patients died. More patients in Cohort I were included in study with a disease duration ≤10 days, while more patients in Cohort II were included in the study with a disease duration >10 days (P = 0.001). The median duration from symptoms onset to tests was significantly longer in Cohort II group compared to Cohort I (first sample median: 11 vs. 8 days, P = 0.003).

Table 1. The clinical characteristics, laboratory tests, treatment, and outcomes of patients with COVID-19 infection depending on the severity of the disease.

Table 2. The results of nasopharyngeal swab and saliva sample tests for SARS-CoV-2 infection depending on the severity of the disease.

Laboratory tests profiles

On admission, Cohort II patients had a higher concentration of interleukin-6 (IL-6, median: 43.5 vs. 16.9 pg/ml, P = 0.001), ferritin (median: 652 vs. 270.2 ng/ml, P = 0.010), C-reactive protein (median: 102.5 vs. 25.7 mg/dl, P < 0.001), D-dimer (median: 465 vs. 285 ng/ml, P = 0.001), lactate dehydrogenase (median: 342 vs. 238 IU/l, P < 0.001) compared to Cohort I patients. There were significant changes in leukocyte counts between these groups with a bigger leukocytosis (median count: 7.2 vs. 4.9 per mm3, P = 0.003), neutrophilia (median count: 5.2 vs. 3.4 per mm3, P < 0.001) and lymphopenia less than 1.1 per mm3 (64.2% vs. 44.7%, P = 0.040) in Cohort II patients. The neutrophil to lymphocyte ratios were significantly higher in these patients too (median: 6.9 vs. 2.6, P < 0.001).

Nasopharyngeal swabs and saliva samples profiles

In total, SARS-CoV-2 was detected in 53% (316/592) of the samples. Out of the 296 tested NPSs and 296 SSs, significantly more SARS-CoV-2 positive samples were in NPS group − 173 (58%) vs. 143 (48%), respectively; P = 0.017 (). Analysing according to the cohorts, the same difference was found in Cohort I (P = 0.011) and was especially evident on the first RT-qPCR test (P = 0.010) (). For the comparison, NPSs and SSs collected from 150 healthy persons were tested and SARS-CoV-2 was detected in 2 NPS and none of SS samples, reaching the specificity of 99% and 100%, respectively.

Figure 1. Detection rates of SARS-CoV-2 in saliva and nasopharyngeal samples collected from COVID-19 patients and controls (a) and in patients dichotomised into mild/moderate (Cohort I, b) and severe cases (Cohort II, c).

Figure 1. Detection rates of SARS-CoV-2 in saliva and nasopharyngeal samples collected from COVID-19 patients and controls (a) and in patients dichotomised into mild/moderate (Cohort I, b) and severe cases (Cohort II, c).

RT-qPCR test Ct values of SSs were significantly lower (reflecting higher viral load) than in NPS samples (mean 28.01 vs. 30.07, P = 0.002, ), and the significant differences between SSs and NPSs were retained when the samples were dichotomised into Cohort I and Cohort II (P =0.022 and P = 0.035, respectively, ), or the Ist test was analysed only (P = 0.05, ). No significant differences in the mean Ct values were found between the cohorts using SS or NPT samples (P > 0.05, ).

Figure 2. Comparison of Cycle threshold (Ct) values from positive saliva specimens and nasopharyngeal swabs between all samples altogether (a), all samples dichotomised into Cohort I vs. Cohort II, (b) and the first test only (c).

Figure 2. Comparison of Cycle threshold (Ct) values from positive saliva specimens and nasopharyngeal swabs between all samples altogether (a), all samples dichotomised into Cohort I vs. Cohort II, (b) and the first test only (c).

The correlation analysis of paired NPS and SS test results for SARS-CoV-2 repeated three times every 3 days was performed. A statistically significant moderately strong correlation on the first (r = 0.54, P < 0.001, ) and weak correlation on the second (r = 0.38, P = 0.029, ) tests between the Ct values of SSs and NPSs was observed; however, there was no significant correlation on the third test (). Significant correlations were found between all SS and NPT samples and between Cohort I and Cohort II samples analysed separately as well ().

Figure 3. Correlation analysis of cycle threshold (Ct) values from positive saliva specimens and nasopharyngeal swabs: corresponding first (a), second (b), and third (c) sample tests, altogether (d) and according to cohorts (e and f).

Figure 3. Correlation analysis of cycle threshold (Ct) values from positive saliva specimens and nasopharyngeal swabs: corresponding first (a), second (b), and third (c) sample tests, altogether (d) and according to cohorts (e and f).

Virus dynamics over time

Analysis of virus detection dynamics over time confirmed that virus loads reflected by lower Ct values in both, NPSs and SSs were higher in the Ist samples and gradually decreased during time in the second and third tests (P < 0.05; ). Separate cohort analysis showed the same tendencies (). Also, the correlation analysis between Ct values from the first SSs and NPSs and days since the onset of symptoms was done (). There was a significant and moderately strong correlation between SSs SARS-CoV-2 RT-qPCR Ct and duration of the disease (r = 0.47, P < 0.001), as well as in the case of NPSs (r = 0.41, P < 0.001). In Cohort I, a stronger correlation was noticed for SS samples than for NPTs (r = 0.55, P = 0.002 vs. r = 0.38, P = 0.017); however, in Cohort II, vice versa association was found (r = 0.34, P = 0.037 vs. r = 0.43, P = 0.009). SARS-CoV-2 RT-qPCR test became negative earlier in SSs in Cohort I as compared to Cohort II (mean 11.7 vs. 14.8 days, P = 0.005), while there was no such difference in NPS tests.

Figure 4. Analysis of serial saliva and nasopharyngeal samples: comparison of cycle threshold (Ct) values (a–c) and correlation of Ct values from saliva, nasopharyngeal and days since onset of symptoms (d–f). r, Pearson correlation.

Figure 4. Analysis of serial saliva and nasopharyngeal samples: comparison of cycle threshold (Ct) values (a–c) and correlation of Ct values from saliva, nasopharyngeal and days since onset of symptoms (d–f). r, Pearson correlation.

SARS-CoV-2 Ct values association with disease severity markers

SSs (N = 61) and NPSs (N = 73) first tests’ SARS-CoV-2 RT-qPCR Ct values correlations with clinical and laboratory characteristics are presented in . There was a strong positive correlation between the number of days since the onset of symptoms and Ct values in both saliva and nasopharyngeal analyses (P < 0.001). Also, Ct values of SSs correlated with NEW score, ferritin, CRP, D-dimer concentrations, WBC, neutrophil and lymphocyte counts, whereas NPSs Ct values correlated with CRP, D-dimer and LDH concentrations, and WBC and neutrophil counts.

Table 3. The correlations between demographic, clinical, laboratory characteristics, and SARS-CoV-2 RT-qPCR Ct values in saliva samples and nasopharyngeal swabs.

Salivary SARS-CoV-2 Ct value as an independent predictor of severe COVID-19

Univariate and multivariate Cox proportional hazards regression analyses for prognosis of severe COVID-19 infection are given in . The results of univariate analyses showed that age, NEW score, Charlson comorbidity index, IL-6, ferritin, and CRP concentration, WBC and neutrophil count, NLR, LDH on admission, and Ct values ≤30 from SSs were associated with severe COVID-19 (). Ct values of NPS were not a prognostic factor for disease severity. In multivariate analysis, only Ct values from SSs retained significant prognostic value and therefore were the independent predictors for severe COVID-19 ().

Table 4. Univariate and multivariate Cox proportional hazards regression analyses for prognosis of severe COVID-19 infection.

In the Kaplan–Meier analysis, no significant difference was found in positivity for SARS-CoV-2 duration between SSs and NPSs when the entire cohort of 100 was analysed (median: 15.0 days vs. 17.0 days, P = 0.076) (). The median duration of SSs positivity in Cohort II was significantly longer than in Cohort I by the log-rank test (16.0 days [CI, 15.2–16.7 days] vs. 14.0 days [CI, 12.0–15.9 days]; P = 0.041) (). However, there was no difference in positivity duration in different COVID-19 disease severity groups using NPSs (16.0 days [CI, 14.5–17.5] vs. 17.0 days [CI, 14.1–19.8]; P = 0.155) ().

Figure 5. Kaplan–Meier analysis of positivity for SARS-CoV-2 prediction from saliva samples and nasopharyngeal swabs when the entire cohort (a) and mild to moderate illness (Cohort I) vs. severe disease (Cohort II) cases are analysed in saliva samples (b) and nasopharyngeal swab samples (c) separately.

Figure 5. Kaplan–Meier analysis of positivity for SARS-CoV-2 prediction from saliva samples and nasopharyngeal swabs when the entire cohort (a) and mild to moderate illness (Cohort I) vs. severe disease (Cohort II) cases are analysed in saliva samples (b) and nasopharyngeal swab samples (c) separately.

Discussion

Quantitative reverse transcription polymerase chain reaction (RT-qPCR) is considered a gold standard confirmatory test for coronavirus disease 2019 (COVID-19). However, the relationship between SARS-CoV-2 Ct values which are directly related to the viral load and the severity of the disease is not fully understood and there is still a lack of tools for timely disease severity prediction. The most commonly used NPSs are uncomfortable and require trained healthcare staff to perform. There is also a risk of generating aerosol, which increases the viral transmission risk in healthcare workers [Citation4]. Saliva samples (SSs) may be obtained safely and are compatible with self-collection. A number of studies with SSs revealed their suitability for the SARS-CoV-2 RNA detection [Citation6–12,Citation17–21]. The performed meta-analyses suggest that SARS-CoV-2 detection using SSs is of similar sensitivity or slightly less sensitive than NPSs [Citation10,Citation22]. Our study supported that SSs are a suitable source for COVID-19 diagnosis and even have some advantages over NPSs.

One of our main findings is that SS Ct values, unlike NPS, are associated with the severity of COVID-19 disease. Although conflicting, there are data that shows that NPSs Ct values are not associated with the severity of the disease [Citation23–26]. We found no differences in NPSs Ct values in mild or severe disease in our study also. However, SSs SARS-CoV-2 RT-qPCR test became negative earlier in mild/moderate disease compared to severe disease (11.7 vs. 14.8 days), while there was no such difference in NPS tests. Even more, we found that Ct ≤ 30 from SSs were the independent predictors for severe COVID-19, while the virus load in NPSs was not a prognostic factor for disease severity. According to our data, we may suppose that in the case of mild/moderate disease, SARS-CoV-2 RT-qPCR test in SSs becomes negative earlier than in NPSs; however, in the case of severe COVID-19 infection, the SS RT-qPCR Ct values stay lower for a longer time. There are not many studies analysing SS SARS-CoV-2 association with the disease severity. Aydin et al. found that the saliva viral load in the early stages of COVID-19 infection may have a prognostic value in predicting disease progression in patients over 45 years of age [Citation20]. Similar results were presented in a study where the mild disease pattern showed the high initial viral load in sputum and saliva after deep cough specimens in the first week with a peak in the viral load during the second week, followed by a decrease, and the severe disease pattern showed high levels even in the third and fourth weeks [Citation27]. More studies are needed, but it seems that longer-lasting low Ct values in saliva could help to suspect a severe form of COVID-19 disease. Therefore, patients with low Ct values of simple saliva SARS-CoV-2 testing could be monitored more carefully, enabling timely hospitalization and administration of appropriate treatment measures by healthcare professionals.

Other interesting finding was that viral load reflected by Ct values in SSs was significantly higher than in NPSs, regardless of the severity of the disease. Low Ct value designates an elevated concentration of genetic material reflecting a higher viral load, typically correlated with a high infection risk [Citation28]. The high value of Ct specifies a low concentration of viral genetic material with less infectivity risk; however, higher values may be seen in a later convalescent stage [Citation29]. Both SSs and NPSs RT-qPCR Ct correlated clearly with duration of the disease in our study, with lower values in the first tests and higher in the third tests. Studies show [Citation20,Citation30,Citation31] that the ACE2 receptor-rich minor salivary gland ductal epithelium and oral mucosal epithelium are the early targets for SARS-CoV-2 infection, and this fact could explain the high viral load in saliva at the early phases of infection. Thus, we may conclude that saliva samples are reliable materials for screening and detection of SARS-CoV-2.

The relationship between SARS-CoV-2 virus load determined by RT-qPCR and biochemical or haematological markers of COVID-19 disease is being extensively investigated; however, the data are ambiguous. Some studies [Citation19,Citation32–35] suggest associations between lower Ct values on admission samples and disease severity markers, including elevated levels of IL-6, LDH, CRP, lymphopenia and an elevated neutrophil/lymphocyte ratio. Interestingly, our data showed the opposite phenomenon – the positive correlation of higher LDH, D-dimer, CRP, leukocyte and neutrophil levels with higher Ct values in NPSs and ferritin, CRP, leukocytes and neutrophil levels in SSs, and obvious correlation with the number of days since the onset of the disease for both SSs and NPSs. However, there was a positive correlation in SSs (but not NPSs) between lymphocyte count and Ct values, with lymphopenia observed more often in the cases of higher viral load in saliva. Chua et al. in the paediatric COVID-19 study also found that children with lymphopenia were more likely to have a higher SSs viral load, but not NPSs viral load, which could support the interaction between COVID-19 infection, especially via the salivary glands, and the evasion of the immune system [Citation36]. The discrepancies in associations between Ct values and biochemical or haematological markers between studies should be interpreted with caution because there is a wide heterogeneity in fluid samples collected at different phases of the disease. Our data showed that despite the differences in inflammatory and haematological markers, the mean Ct values in either SSs or NPSs were similar in severe disease or mild/moderate patient groups. However, there were more patients in the severe disease group included in our study at a later time and more patients in the mild/moderate disease group included in our study earlier. We may suppose that heightened inflammatory markers in severe patient group were influenced by the later stage of the disease, when hyperactive inflammatory response was involved in the disease process more than the virus by itself.

Despite the relevance of the data, there are some limitations of our study. The protocol of the study did not provide the inclusion of all subjects on a specific day of illness, so that the dynamics of Ct, the association with disease severity or laboratory parameters could be more accurately assessed. Consequently, despite the different clinical course, the differences in Ct values between the two cohorts could have been influenced by the time difference in symptom onset and sample collection. The inequalities between the positive corresponding NPT and SS samples (positive NPT, but negative SS and vice versa) might have also affected the final result. Another drawback to look out for is that RT-qPCR results cannot distinguish between viable and non-viable virus and reflect the total amount of viral RNA. Moreover, the SARS-CoV-2 quantification was assessed only in Ct values, not in viral load, which could have also limited the evaluation of active disease.

Conclusions

Our study supports evidence that salivary SARS-CoV-2 RNA RT-qPCR testing is a valuable diagnostic tool for COVID-19 showing that cycle threshold values ≤30 are associated with severe COVID-19. Moreover, long-lasting positivity in salivary SARS-CoV-2 RNA RT-qPCR testing provides a longer window for the severe COVID-19 diagnosis, enabling healthcare professionals to select appropriate treatment options for patients.

Author contributions

KS and BZ performed data analysis and drafted the manuscript. BZ, RS, KZ, EM, TS, MP, and JG participated in sample collection and handled the clinical data. KS, BZ, RZ, and ML performed statistical analysis. JG, AP, LJ, and SJ carried out the study design. All authors read and approved the final version of the manuscript.

Acknowledgments

The authors would like to acknowledge Jonas Belevicius, Vaida Kurmauskaite, and Ieva Asmenaviciute for laboratory testing.

Disclosure statement

Mantvydas Lopeta and Julius Gagilas are employees at JSC Diagnolita. JSC Diagnolita has developed the kits used in this research. The authors declare no other conflicts of interest.

Additional information

Funding

No funding was received for conducting this study.

References

  • Weekly epidemiological update on COVID-19 - 6 Apr 2023 [Internet]. [updated 2023 Apr 7]. Available from: https://www.who.int/publications/m/item/weekly-epidemiological-update-on-covid-19—6-april-2023
  • Lippi G, Simundic AM, Plebani M. Potential preanalytical and analytical vulnerabilities in the laboratory diagnosis of coronavirus disease 2019 (COVID-19). Clin Chem Lab Med. 2020 Jun 25;58(7):1070–12. DOI:10.1515/cclm-2020-0285.
  • Khiabani K, Amirzade-Iranaq MH. Are saliva and deep throat sputum as reliable as common respiratory specimens for SARS-CoV-2 detection? A systematic review and meta-analysis. Am J Infect Control. 2021 Sep;49(9):1165–1176. doi:10.1016/j.ajic.2021.03.008.
  • Harding H, Broom A, Broom J. Aerosol-generating procedures and infective risk to healthcare workers from SARS-CoV-2: the limits of the evidence. J Hosp Infect. 2020 Aug;105(4):717–725. doi:10.1016/j.jhin.2020.05.037.
  • Karligkiotis A, Arosio A, Castelnuovo P. How to obtain a nasopharyngeal Swab Specimen. N Engl J Med. 2020 Jul 16;383(3):e14. DOI:10.1056/NEJMc2015949.
  • Bruxvoort K, Tenggardjaja CF, Slezak J, et al. Variation in SARS-CoV-2 molecular test sensitivity by specimen types in a large sample of emergency department patients. Am J Emerg Med. 2021 Dec 1;50:381–387. DOI:10.1016/j.ajem.2021.08.034.
  • LeGoff J, Kernéis S, Elie C, et al. Evaluation of a saliva molecular point of care for the detection of SARS-CoV-2 in ambulatory care. Sci Rep. 2021 Oct 26;11(1):21126. DOI:10.1038/s41598-021-00560-8.
  • Uddin MKM, Shirin T, Hossain ME, et al. Diagnostic performance of self-collected saliva versus nasopharyngeal swab for the molecular detection of SARS-CoV-2 in the clinical setting. Microbiol Spectr. 0(0): e00468-21. 10.1128/Spectrum.00468-21.
  • Williams E, Bond K, Zhang B, et al. Saliva as a noninvasive specimen for detection of SARS-CoV-2. J Clin Microbiol. 2020 Jul 23;58(8): e00776-20. DOI: 10.1128/JCM.00776-20.
  • Czumbel LM, Kiss S, Farkas N, et al. Saliva as a candidate for COVID-19 diagnostic testing: a meta-analysis. Front Med. 2020 Aug 4;7:465. DOI:10.3389/fmed.2020.00465.
  • Palikša S, Lopeta M, Belevičius J, et al. Saliva testing is a robust non-invasive method for SARS-CoV-2 RNA detection. Infect Drug Resist. 2021;14:2943–2951.
  • Tan SH, Allicock O, Armstrong-Hough M, et al. Saliva as a gold-standard sample for SARS-CoV-2 detection. Lancet Respir Med. 2021;9:562–564.
  • IDSA and AMP joint statement on the use of SARS-CoV-2 PCR cycle threshold (Ct) values for clinical decision-making [Internet]. [updated 2022 May 10]. Available from: https://www.idsociety.org/globalassets/idsa/public-health/covid-19/idsa-amp-statement.pdf.
  • Yoon JG, Yoon J, Song JY, et al. Clinical significance of a high SARS-CoV-2 viral load in the saliva. J Korean Med Sci. 2020 May 20;35(20):e195. DOI:10.3346/jkms.2020.35.e195.
  • Shenoy S. SARS-CoV-2 (COVID-19), viral load and clinical outcomes; lessons learned one year into the pandemic: a systematic review. World J Crit Care Med. 2021 Jul 9;10(4):132–150. DOI:10.5492/wjccm.v10.i4.132.
  • COVID-19 treatment guidelines [Internet]. [updated 2023 March 31]. Available from: https://files.covid19treatmentguidelines.nih.gov/guidelines/covid19treatmentguidelines.pdf
  • Leung ECM, Chow VCY, Lee MKP, et al. Deep throat saliva as an alternative diagnostic specimen type for the detection of SARS-CoV-2. J Med Virol. 2021 Jan;93(1):533–536. doi:10.1002/jmv.26258.
  • Wong SCY, Tse H, Siu HK, et al. Posterior oropharyngeal saliva for the detection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Clin Infect Dis off Publ Infect Dis Soc Am. 2020 Jun 21 ciaa797;71:2939–2946. 10.1093/cid/ciaa797.
  • Azzi L, Carcano G, Gianfagna F, et al. Saliva is a reliable tool to detect SARS-CoV-2. J Infect. 2020 Jul;81(1):e45–50. doi:10.1016/j.jinf.2020.04.005.
  • Aydin S, Benk IG, Geckil AA. May viral load detected in saliva in the early stages of infection be a prognostic indicator in COVID-19 patients? J Virol Methods. 2021 Aug;294:114198. DOI:10.1016/j.jviromet.2021.114198.
  • Nagura-Ikeda M, Imai K, Tabata S, et al. Clinical evaluation of self-collected saliva by quantitative reverse transcription-PCR (RT-qPCR), direct RT-qPCR, reverse transcription-loop-mediated isothermal amplification, and a rapid antigen test to diagnose COVID-19. J Clin Microbiol. 2020 Aug 24;58(9). 10.1128/JCM.01438-20 e01438-20.
  • Ibrahimi N, Delaunay-Moisan A, Hill C, et al. Screening for SARS-CoV-2 by RT-PCR: saliva or nasopharyngeal swab? Rapid review and meta-analysis. PLoS ONE. 2021;16(6):e0253007.
  • Salto-Alejandre S, Berastegui-Cabrera J, Camacho-Martínez P, et al. SARS-CoV-2 viral load in nasopharyngeal swabs is not an independent predictor of unfavorable outcome. Sci Rep. 2021 Jun 21;11(1):12931. DOI:10.1038/s41598-021-92400-y.
  • Maltezou HC, Raftopoulos V, Vorou R, et al. Association between upper respiratory tract viral load, comorbidities, disease severity, and outcome of patients with SARS-CoV-2 infection. J Infect Dis. 2021 Apr 8;223(7):1132–1138. DOI:10.1093/infdis/jiaa804.
  • Argyropoulos KV, Serrano A, Hu J, et al. Association of initial viral load in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) patients with outcome and symptoms. Am J Pathol. 2020 Sep;190(9):1881–1887. doi:10.1016/j.ajpath.2020.07.001.
  • Lee S, Kim T, Lee E, et al. Clinical course and molecular viral shedding among asymptomatic and symptomatic patients with SARS-CoV-2 infection in a community treatment center in the Republic of Korea. JAMA Intern Med. 2020 Nov 1;180(11):1447–1452. DOI:10.1001/jamainternmed.2020.3862.
  • Zheng S, Fan J, Yu F, et al. Viral load dynamics and disease severity in patients infected with SARS-CoV-2 in Zhejiang province, China, January-March 2020: retrospective cohort study. BMJ. 2020 Apr 21;369:m1443.
  • Rao SN, Manissero D, Steele VR, et al. A systematic review of the clinical utility of cycle threshold values in the context of COVID-19. Infect Dis Ther. 2020 Sep;9(3):573–586. doi:10.1007/s40121-020-00324-3.
  • Sarkar B, Sinha RN, Sarkar K. Initial viral load of a COVID-19-infected case indicated by its cycle threshold value of polymerase Chain Reaction could be used as a predictor of its transmissibility - an experience from Gujarat, India. Indian J Community Med off Publ Indian Assoc Prev Soc Med. 2020 Sep;45(3):278–282. doi:10.4103/ijcm.IJCM_593_20.
  • Liu L, Wei Q, Alvarez X, et al. Epithelial cells lining salivary gland ducts are early target cells of severe acute respiratory syndrome coronavirus infection in the upper respiratory tracts of rhesus macaques. J Virol. 2011 Apr;85(8):4025–4030. doi:10.1128/JVI.02292-10.
  • Xu H, Zhong L, Deng J, et al. High expression of ACE2 receptor of 2019-nCov on the epithelial cells of oral mucosa. Int J Oral Sci. 2020 Feb 24;12(1):1–5. DOI:10.1038/s41368-020-0074-x.
  • Rodríguez-Serrano DA, Roy-Vallejo E, Zurita Cruz ND, et al. Detection of SARS-CoV-2 RNA in serum is associated with increased mortality risk in hospitalized COVID-19 patients. Sci Rep. 2021 Jun 23;11(1):13134. DOI:10.1038/s41598-021-92497-1.
  • Chen X, Zhao B, Qu Y, et al. Detectable serum severe acute respiratory syndrome coronavirus 2 viral load (RNAemia) is closely correlated with drastically elevated interleukin 6 level in critically Ill patients with coronavirus disease 2019. Clin Infect Dis off Publ Infect Dis Soc Am. 2020 Nov 5;71(8):1937–1942. DOI:10.1093/cid/ciaa449.
  • Xia XY, Wu J, Liu HL, et al. Epidemiological and initial clinical characteristics of patients with family aggregation of COVID-19. J Clin Virol off Publ Pan Am Soc Clin Virol. 2020 Jun;127:104360. DOI:10.1016/j.jcv.2020.104360.
  • Liu Y, Liao W, Wan L, et al. Correlation between relative nasopharyngeal virus RNA load and lymphocyte count disease severity in patients with COVID-19. Viral Immunol. 2021 Jun;34(5):330–335. doi:10.1089/vim.2020.0062.
  • Chua GT, Wong JSC, KKW T, et al. Saliva viral load better correlates with clinical and immunological profiles in children with coronavirus disease 2019. Emerg Microbes Infect. 2021 Dec;10(1):235–241. doi:10.1080/22221751.2021.1878937.