Zeng J, Xu H, Luo S, Zhou X, et. al., DOI https://doi.org/10.2147/IDR.S515224
Background:
Understanding the factors influencing viral clearance in hospitalized COVID-19 patients, including vaccination status and antiviral therapy, is critical for optimizing clinical management.
Methods: 1,424 hospitalized COVID-19 patients retrospectively included from four tertiary hospitals in Hainan Province between March and December 2022. Viral clearance was defined as the interval from hospital admission to the first of two consecutive RT-PCR tests with Ct values ≥ 35. Clinical data, vaccination history, and antiviral treatment were collected. A generalized linear mixed model and Robust regression were used to assess viral clearance dynamics and their predictors.
Results: Delayed viral clearance was independently associated with advanced age (p < 0.001), male sex (p = 0.006), hypertension (p < 0.001), coronary heart disease (p = 0.004), ICU admission (p < 0.001), and mechanical ventilation (p < 0.001).Patients receiving ≥ 2 inactivated vaccine doses had significantly higher baseline Ct values (median 29.75 vs 28.75, p = 0.014), shorter time to viral negativity (6.3 vs 7.4 days, p < 0.001), and reduced hospital stay (11.2 vs 12.7 days, p < 0.001). Among these, patients vaccinated ≥ 360 days prior had shortest negative conversion time (5.6 days) and shortest hospitalization (10.3 days).Antiviral therapy with Nirmatrelvir-ritonavir (N/R) accelerated viral clearance more effectively than Azvudine (2.29 vs 1.82 Ct/day, p = 0.045) and no antiviral treatment (1.88 Ct/day, p = 0.041), Although NAT-treated patients achieved viral negativity more rapidly (6.2 days, p = 0.013), N/R demonstrated superior clearance rate. Hospital stays were shorter with N/R than Azvudine (12.1 vs 13.5 days, p = 0.015).
Conclusion: Viral clearance dynamics in hospitalized COVID-19 patients are influenced by age, comorbidities, vaccination, and antiviral treatment. Administration of ≥ 2 inactivated vaccine doses—especially ≥ 360 days apart—and early N/R therapy may accelerate viral clearance and reduce hospital stay.
Introduction
SARS-CoV-2 persists globally, albeit with reduced transmission intensity.1 Despite widespread vaccination and antiviral use having alleviated the overall burden of COVID-19,2 a detailed understanding of viral clearance dynamics remains critical. Viral load not only correlates with disease severity,3 transmissibility, and treatment response,4–6 but also serves as a proposed surrogate marker for therapeutic efficacy. Accordingly, elucidating the kinetics of viral decline is essential for informing clinical decision-making.
Certain populations are disproportionately affected by severe COVID-19. Elderly individuals are more likely to develop acute respiratory distress syndrome (ARDS) and multiple organ dysfunction syndrome (MODS), while male sex is consistently associated with worse outcomes.7–10 Children, by contrast, exhibit more robust immune responses and milder disease courses.8,11 Comorbidities such as hypertension, coronary heart disease (CHD), and diabetes further increase severity risks.12 However, how these factors influence viral clearance dynamics remains poorly understood—representing a crucial gap in tailoring patient-specific therapeutic strategies for those underlying conditions.
Vaccination continues to serve as the cornerstone of COVID-19 prevention. Although mRNA vaccines demonstrate high efficacy,13 uncertainties persist regarding the durability of protection and optimal dosing intervals for inactivated vaccines. Addressing these knowledge gaps is crucial for maintaining herd immunity, particularly in settings where inactivated vaccines are the primary immunization strategy.
Antiviral therapies provide essential adjunctive benefits, especially among high-risk individuals with elevated viral loads. Existing evidence supports their role in reducing hospitalization rates and mortality.14,15 Nonetheless, the comparative effects of individual antiviral agents on viral clearance dynamics remain incompletely characterized—hindering the development of precision-based antiviral strategies.
Bridging these critical knowledge gaps requires comprehensive real-world data that integrate host factors, immunization timelines, and antiviral regimens. Large-scale observational studies are urgently needed to unravel these complex interactions and inform precision-guided interventions aimed at enhancing viral clearance, shortening hospitalization, and mitigating severe disease outcomes in at-risk cohorts.
Patients and Methods
Study Design and Setting
This retrospective multicenter cohort study included 1,862 patients hospitalized with confirmed COVID-19 between March and December 2022 at four tertiary hospitals in Hainan Province, China: Hainan General Hospital, Haikou People’s Hospital, Hainan Western Central Hospital, and Sanya Central Hospital. Standardized admission protocols were uniformly implemented in accordance with the Diagnosis and Treatment Protocol for COVID-19 (Trial Version 9, China).
Eligibility Criteria
Eligible patients were aged ≥14 years, had laboratory-confirmed SARS-CoV-2 infection (via RT-PCR), classic respiratory symptoms (eg, cough, fever, or dyspnea); A positive SARS-CoV-2 nucleic acid test (RT-PCR); and Radiological evidence of pulmonary involvement on chest X-ray or ground-glass opacities on CT, and had complete medical and virological records including serial Ct values. Patients were excluded if they had lacked RT-PCR data, incomplete vaccination or treatment history, or comorbidities such as malignancy, chronic kidney disease, chronic lung disease, or immunosuppressive conditions. Following exclusions, 1,424 patients were included in the final analytic cohort.
Data Collection and Variables
Clinical and demographic data were extracted from standardized electronic health records. Variables included age, sex, comorbidities (hypertension, diabetes, coronary heart disease), ICU admission, mechanical ventilation, COVID-19 vaccination status (Unvaccinated, Partially vaccinated 1 dose, Fully vaccinated 2 doses, Boosted more than 3 doses, interval to infection), and antiviral treatment regimen (Nirmatrelvir/Ritonavir, Azvudine, or no therapy). Patients were stratified accordingly for comparative analysis [Figure 1].
Outcome Definitions
The primary outcome was time to viral clearance, defined as the number of days from hospital admission to the first of two consecutive negative RT-PCR tests (Ct ≥35 for ORF1ab and N genes, ≥24 hours apart). Viral clearance rate (Ct/day) was calculated from serial Ct values. The secondary outcome was hospitalization duration.
Viral Screen
Nasopharyngeal swabs were collected every 2–3 days. SARS-CoV-2 RNA was quantified using CFDA-approved qRT-PCR assays (GeneoDX, Shanghai; Daan Gene, Guangzhou), targeting the ORF1ab and N genes. Viral load was reported as log10 RNA copies/mL and Ct values. Discharge required clinical improvement, ≥3 days of normal temperature, and two negative tests ≥24 hours apart.
Ethics Approval and STROBE List
The study was approved by the Institutional Review Board of Hainan General Hospital (Approval No: Med-Eth-Re [2025]152). As this was a retrospective cohort study based on anonymized medical records, the requirement for informed consent was waived by the IRB. All patient data were de-identified prior to analysis to ensure confidentiality. The study was conducted in accordance with the principles outlined in the Declaration of Helsinki and relevant local regulations governing medical research. This study is reported in accordance with the “Strengthening the Reporting of Observational Studies in Epidemiology” (STROBE) statement for cohort studies in Supplementary Table 1.
Statistical Analysis
All Statistical analyses were conducted using R software (version 4.2.3), employing the tidyverse, lme4, robustbase, and ggplot2 packages. Continuous variables were summarized as means ± standard deviations or medians with interquartile ranges, depending on distribution assessed via the Shapiro–Wilk test. Categorical variables were reported as counts and percentages.
Group comparisons were performed using chi-square or Fisher’s exact tests for categorical variables, and either independent-sample t-tests or Mann–Whitney U-tests for continuous variables, based on normality.
To examine temporal changes in viral load, a generalized linear mixed model (GLMM) was used to assess Ct value trajectories over time.
For predictors of viral clearance rate, univariate analyses were first conducted, followed by multivariate mixed-effects regression models adjusting for key confounders: age, sex, comorbidities, vaccination status, and antiviral therapy.
Additionally, robust linear regression was applied to assess the association of clinical variables with hospitalization duration and time to viral negativity, accounting for potential outliers and heteroscedasticity.
All statistical tests were two-tailed, with statistical significance defined as p < 0.05. Missing data were handled via complete case analysis, as the proportion of missingness was <5%.
Results
Detailed patient characteristics of 1,424 hospitalized COVID-19 patients included are provided in Table 1.
Male and Age
The male (n = 601, 42.2%) had significant slower viral clearance than female patients (P=0.006) as showed in trend linear [Figure 2A] and violin plots [Figure 2B], with males demonstrating a slower viral clearance rates compared to females. As shown in Figure 3A and B, higher age groups were associated with progressively slower viral clearance rates. The age group of 0–18 years old (n = 296, 20.8%) demonstrated significantly faster viral clearance compared to other age groups (0–18 years vs 18–35 years, P<0.001; 0–18 years vs 35–65 years, P<0.001; 0–18 years vs >65 years, P<0.001). Additionally, the age group of 18–35 years old (n = 297, 20.9%) showed faster viral clearance compared to the aged group of 35–65 years old (n = 508, 35.7%) (18–35 years vs 35–65 years, P<0.001). In contrast, the age group (>65 years) (n = 323, 22.7%) exhibited the slowest viral clearance rates in comparison (18–35 years vs >65 years, P<0.001; 35–65 years vs >65 years, P<0.001).
Comorbidities of Hypertension, CHD and Diabetes
Patients with hypertension (n = 205,14.4%, P<0.001) [Figure 4A and B] or CHD (n = 55,3.9%, P=0.004) [Figure 4C and D] exhibited significantly reduced viral clearance rates compared to those without hypertension and CHD. Conversely, diabetes (n = 88, 6.2%) was not significantly associated with Ct dynamics (P = 0.333) in our analysis [Figure 4E and F].
Vaccination
All COVID-19 vaccines administered in China are inactivated vaccines, Figure 5 illustrates the dynamic changes in Ct values, time to viral negativity, and duration of hospitalization across vaccination groups (0–1 dose vs 2–3 doses). There was no statistically significant difference in Ct values between patients who received a single dose of the COVID-19 vaccine and those who were unvaccinated. However, patients who received 2–3 doses of the COVID-19 vaccine (n = 749,52.6%) had higher initial Ct values (Ct median =29.75) compared to those vaccinated less than 2 doses (n = 675,47.4%) (Ct median =28.75, P=0.014) in trend linear [Figure 5A] and violin plots [Figure 5B], but Patients with 2–3 vaccine doses showed no significant difference in viral clearance rate (P = 0.809), Moreover, vaccination 2–3 doses showed rapid in negative conversion (6.3 vs 7.4 days, P<0.001) [Figure 5C and Table 2] and decrease hospital stays (11.2 vs 12.7 days,P<0.001) [Figure 5D and Table 2] as compared to vaccination less than 2 doses.
In patients who received 2–3 vaccine doses, the day last vaccination to hospital day divided according to 0–180 days,180–360 days, above 360 days and 0–1dose as reference [Figure 6A and B], as shows in Table 3, those whose last vaccination occurred more than 360 days before hospitalization (n = 152 of 749, 20.3%) had highest initial Ct values (Ct median =30.67) compared to those vaccinated within 360 days prior to hospitalization (0–180 days,Ct median =29.59;180–360 days,Ct median =29.64). This trend was consistent across comparisons, with significantly higher baseline Ct values observed in the >360 days group compared to the 180–360 days group (P = 0.013). Specifically, patients vaccinated over 360 days ago achieved negative conversion time in a mean of 5.6 days, significantly earlier than patients vaccinated less than 2 doses, who became viral negative in 7.4 days (P < 0.001) [Figure 6C]. Similarly, the mean hospital stay was shorter in patients vaccinated more than 360 days ago (10.3 days) compared to patients vaccinated less than 2 doses (12.7 days) (P < 0.001) [Figure 6D].
Viral Clearance
The study population included 87 patients in the N/R group (6.1%, mean [SD] age, 49.3 [19.1] years), 332 in the Azvudine group (23.3%, mean [SD] age, 51.8 [18.1] years), and 1005 in the NAT group (70.6%, mean [SD] age, 40.4 [24.0] years).
Despite lower initial cycle threshold (Ct) values in the N/R (25.92) and Azvudine (28.62) groups compared to NAT (31.16) (N/R vs Azvudine, P = 0.007; N/R vs NAT, P<0.001) [Figure 7A], N/R demonstrated superior therapeutic effects. N/R was associated with faster viral clearance rates (2.29 vs 1.82 increase in Ct value per day for Azvudine, P = 0.045) [Figure 7B and Table 4], a shorter time to viral negativity (7.3 vs 8.5 days for Azvudine, P = 0.013) [Figure 7C and Table 4], and reduced hospital stays (12.1 vs 13.5 days for Azvudine, P = 0.015) [Figure 7D and Table 4].
Compared to NAT, the N/R group also exhibited faster viral clearance (2.29 vs 1.88 increase in Ct value per day, P = 0.041) [Table 4], but NAT patients achieved viral negativity more quickly (6.2 vs 7.3 days, P = 0.013) [Figure 7C and Table 4]. No significant difference was observed in hospital stay between the N/R and NAT groups (12.1 vs 11.4 days, P = 0.27) [Figure 7D and Table 4].
In robust regression analysis (Figure 8 and Table 5), N/R showed a significant association with shorter negative conversion time (1.09 days [95% CI: 0.29–1.88], P = 0.007), compared to Azvudine (2.10 days [95% CI: 1.56–2.64], P < 0.001). Vaccination with two or more doses significantly reduced time to viral clearance (–1.33 days [95% CI: –1.73 to –0.93], P < 0.001), with the greatest effect observed in individuals vaccinated more than 360 days prior (–1.95 days [95% CI: –2.54 to –1.36], P < 0.001). Age showed a positive correlation with delayed clearance (0.02 days per year [95% CI: 0.01–0.03], P < 0.001), while gender and comorbidities were not statistically significant.
For hospitalization stay (Figure 9 and Table 6), Azvudine was associated with prolonged hospital stays (2.12 days [95% CI: 1.49–2.75], P < 0.001), while N/R was not statistically significant (0.64 days [95% CI: –0.41 to 1.68], P = 0.2309). Receiving two or more vaccine doses significantly reduced hospital days (–1.55 [95% CI: –2.05 to –1.05], P < 0.001), most notably for those vaccinated over 360 days ago (–2.35 [95% CI: –3.12 to –1.58], P < 0.001). Age remained a minor but significant factor (0.02 per year [95% CI: 0.01–0.03], P = 0.005).
ICU Admission and Mechanical Ventilation
The viral clearance rate in patients admitted to the ICU (n = 12[0.8%]) was significantly lower (P < 0.001) compared to non-ICU patients. ICU patients had significantly lower initial Ct values (23.3 vs 29.4), indicating higher viral load at baseline [Figure 10A]. Furthermore, patients requiring mechanical ventilation (n = 15[1.1%]) had a median Ct value of 25.7 on their first nucleic acid test, and their viral clearance rate was also significantly lower (P < 0.001) compared to those who did not require mechanical ventilation, with a median Ct value of 29.4 [Figure 10B].
Discussion
This multicenter retrospective cohort study is among the first to evaluate real-world SARS-CoV-2 viral clearance kinetics in hospitalized patients during the Omicron BA.5.1.3 wave in China, focusing on the impact of inactivated vaccination schedules and antiviral therapy. Three key findings emerged: (1) ≥2 doses of inactivated vaccines spaced ≥360 days apart were associated with higher baseline Ct values, faster viral negativity, and shorter hospitalization; (2) Delayed viral clearance was independently associated with advanced age, male sex, comorbidities, and ICU admission; (3) Early Nirmatrelvir/Ritonavir use enhanced viral clearance more effectively than Azvudine or no therapy. Patients vaccinated ≥360 days prior showed the best outcomes, supporting durable immune memory and the need for frequent boosters.
Inactivated vaccines, as the primary strategy in China, provided clear protection by reducing viral loads and hospital stays. Patients vaccinated over 360 days prior to hospitalization showed the highest Ct values and fastest viral clearance, implying prolonged immune efficacy. This novel observation implies that the immunity induced by whole-virus inactivated vaccines may have greater durability, supporting a less frequent booster schedule for populations.16–20
Our results complement earlier reports showing that mRNA vaccines reduce severity over short intervals,21,22 while inactivated vaccines maintain efficacy up to 20 weeks.23–25 We extend this understanding by demonstrating beneficial effects at >360 days, providing rare long-term immunogenicity data in real-world hospitalized settings. Immunologically, inactivated vaccines present the full viral proteome, stimulating broader epitope recognition and supporting T follicular helper cell–mediated memory.26,27 This may explain enhanced protection seen with longer intervals, and highlights the value of diversified antigenic stimulation in durable immunity.
Consistent with previous evidence, ICU patients had lower baseline Ct values and slower viral clearance, reaffirming Ct as a prognostic marker of severity and transmissibility.28–30 Patients requiring mechanical ventilation showed similar trends, reinforcing the need for early viral suppression strategies in high-risk groups.
Among antivirals, Nirmatrelvir/Ritonavir (N/R) was more effective than Azvudine in improving viral clearance metrics, especially in patients with higher initial viral loads. Although NAT-treated patients reached viral negativity faster in some comparisons, the rate of Ct rise was highest in the N/R group, indicating stronger antiviral activity.31–33 These findings support early targeted antiviral use, particularly in elderly or comorbid individuals.34
Comorbidities such as CHD, hypertension, and advanced age delayed viral clearance, likely via mechanisms of immune senescence and persistent inflammation.35,36 Upregulated ACE2 expression in cardiovascular tissues may facilitate enhanced viral replication, prolonging viral shedding.37–43
In contrast, pediatric and young adult patients had faster viral clearance, consistent with lower ACE2 expression and more robust adaptive responses.41,44–47 Age-related reductions in interferon-γ, IL-2, and memory B cell quality further explain impaired clearance in the elderly.48,49
This study has several limitations. Data were collected during China’s zero-COVID policy, potentially limiting generalizability. Residual confounding may persist due to unmeasured variables such as healthcare access and diagnostic delays. Selection bias is possible due to unequal antiviral availability. Additionally, patients with immunosuppressive comorbidities were excluded, narrowing external validity. Lastly, the absence of post-vaccine antibody titers precludes direct immunological correlation.
In conclusion, this study demonstrates that viral clearance among hospitalized COVID-19 patients is strongly shaped by vaccination interval, antiviral strategy, and host risk profile. Receipt of ≥2 doses of inactivated vaccines spaced ≥360 days apart was associated with lower viral load and improved clinical outcomes. Moreover, early N/R treatment enhanced viral suppression in high-risk groups. These findings support a precision-medicine strategy integrating vaccination schedules and antiviral timing to optimize outcomes.