key: cord-0908629-8pjg8uyw
authors: Alamer, Ahmad; Alrashed, Ahmed A.; Alfaifi, Mashael; Alosaimi, Bandar; AlHassar, Fatimah; Almutairi, Malak; Howaidi, Jude; Almutairi, Wedad; Mohzari, Yahya; Sulaiman, Tarek; Al-jedai, Ahmed; Alajami, Hamdan N.; Alkharji, Fatima; Alsaeed, Ali; Alali, Alaa H.; Baredhwan, Abdullah A.; Abraham, Ivo; Almulhim, Abdulaziz S.
title: Effectiveness and safety of favipiravir compared to supportive care in moderately to critically ill COVID-19 patients: a retrospective study with propensity score matching sensitivity analysis
date: 2021-05-19
journal: Current medical research and opinion
DOI: 10.1080/03007995.2021.1920900
sha: 0b6e42f6a626e7cb01f5973c454bb0ebba5d0a8d
doc_id: 908629
cord_uid: 8pjg8uyw

INTRODUCTION: Favipiravir is a repurposed drug to treat coronavirus 2019 (COVID-19). Due to a lack of available real-world data, we assessed its effectiveness and safety in moderately to critically ill COVID-19 patients. METHODS: This retrospective study was conducted in two public/specialty hospitals in Saudi Arabia. We included patients ≥18 years) admitted April–August 2020 with confirmed SARS-CoV-2 diagnosed by real-time polymerase chain reaction (RT-PCR) from nasopharyngeal swab. Patients received either favipiravir (1800 mg or 1600 mg twice daily loading dose, followed by 800 mg or 600 mg twice daily) or supportive-care treatment. Patients were excluded if they were outside the study period, classified as having a mild form of the disease per WHO criteria, or had an incomplete patient file. Kaplan–Meier (KM) models were used to estimate median time to discharge. Discharge ratios, progression to mechanical ventilation, and mortality outcomes were estimated across the severity spectrum using Cox proportional-hazards models. As a sensitivity analysis, we performed propensity score-matching (PSM) analysis. RESULTS: Overall, median time to discharge was 10 days (95%CI = 9–10) in the favipiravir arm versus 15 days (95%CI = 14–16) in the supportive-care arm. The accelerated discharge benefit was seen across the COVID-19 spectrum of severity. The adjusted discharge ratio was 1.96 (95%CI = 1.56–2.46). Progression to mechanical ventilation was slower with favipiravir (HR(adj) = 0.10, 95%CI = 0.04–0.29). There was no significant effect on mortality (HR(adj) = 1.56, 95%CI = 0.73–3.36). There was a statistically non-significant trend toward worse outcomes in the critical category (HR(adj) = 2.80, 95%CI = 0.99–7.89). Age was an independent risk factor for mortality in mechanically ventilated patients. PSM analyses confirmed these findings. CONCLUSION: Favipiravir was associated with clinical benefits, including accelerated discharge rate and less progression to mechanical ventilation; however, no overall mortality benefits were seen across the severity spectrum.

The coronavirus disease 2019 pandemic is caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which belongs to a large family of singlestranded RNA (þssRNA) implicated in past epidemics over

The clinical course can be fatal in patients with comorbidities and in elderly populations. In just seven days, patients can deteriorate to acute respiratory distress syndrome (ARDS), sepsis, and/or multiple organ failure 4 .

Remdesivir (RDV) was synthesized in 2017 to fight the Ebola outbreak and has shown some in vivo and in vitro activity against SARS-CoV and Middle East respiratory syndrome coronavirus (MERS-CoV) 5 . The Adaptive COVID-19 Treatment Trial (ACTT-1) study group published a preliminary clinical trial report suggesting several clinical benefits in COVID-19 patients. In this double-blind, randomized, placebo-controlled trial (RCT), intravenous (IV) RDV was shown to reduce time to recovery by 11 days (95% confidence interval [CI], [9] [10] [11] [12] compared with 15 days (95%CI, [13] [14] [15] [16] [17] [18] [19] . In the subgroup analysis, recovery-rate ratio was significant in patients on oxygen therapy at 1.47 (95%CI 1. 17-1.84) . No effect on mortality (14 or 28 days) was seen, with a hazard ratio of 0.70; 95%CI, 0.47-1.04 6 . The US Food and Drug Administration (FDA) recently approved the use of RDV in COVID-19 patients 7 .

Favipiravir (FPV), an RNA-dependent RNA polymerase (RdRP) inhibitor with a broad spectrum of activity, is an oral drug synthesized in 2002 and approved in 2004 to treat influenza in Japan 8 . Proposed dosage, though still being formally investigated for coronaviruses, ranges between 1600 and 2000 mg BID (twice daily) followed by a maintenance dosage 300-1800 mg BID 9 . Data on FPV's efficacy and safety are very limited. A recent open-label RCT by Cai et al. found that a 14-day FPV regimen (dosage 1600 mg BID followed by 600 mg) plus interferon-alfa (aerosol inhalation 5 million units) cleared SARS-CoV-2 infection more rapidly than treatment with lopinavir-ritonavir (400 mg/100 mg) with interferon-alfa combination therapy. This was evidenced by significant improvements in chest imaging compared with the control arm (91.43% versus 62.22%; p ¼ .004). In this study, fewer side effects were reported in the FPV arm than in the control arm (11.43% versus 55.56%; p < .001) 10 . Another open-label, prospective RCT compared FPV (600 mg) to umifenovir (another antiinfluenza treatment approved in Russia and China) and found that the clinical-recovery rate at day 7 did not differ significantly in the two arms (p ¼ .139). However, FPV has led to shorter latencies of relief for both pyrexia and cough (p < .001) 11 .

Regarding safety, a systematic review by Pilkington et al. showed that FPV is well-tolerated except for elevation in uric acid levels. Per the systematic review, teratogenicity and QTc prolongation have not yet been adequately studied 12 . With its recent availability in Saudi Arabia, the Ministry of Health (MOH) included FPV as an option for treating moderate to critically ill COVID-19 patients with the following dosage: 1800 mg BID followed by 800 mg BID 13 . Due to the urgent nature of the pandemic and the relative infeasibility of a prospective RCT in the short term, we designed a retrospective trial to evaluate the effectiveness and safety of FPV in moderate-to-severe/critically ill COVID-19 patients when compared to supportive care.

This was a retrospective study conducted in two cluster hospitals in Riyadh: King Fahad Medical City (KFMC; 1200 beds) and Prince Mohammed Bin Abdulaziz Hospital (PMAH; 500 beds), both of which are nationally designated hospitals for COVID-19 and to which complicated patients from all over Saudi Arabia are transferred. The study was approved by the Institutional Review Boards at KFMC and PMAH (IRB 20-477E, July 2020). Informed consent was waived as this was considered an exempt study. In our report, we adhered to the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) Statement checklist 14 .

A list of deceased or discharged patients who had received FPV or supportive care and had been admitted between April and August 2020 to KFMC or PMAH for COVID-19 management was obtained. Through a random-selection process, we screened patients for eligibility. In random selection, each patient's record has an equal opportunity to be selected for coding to minimize sampling bias 15 . The inclusion criteria were as follows: age !18 years, having received at least two doses of FPV (1800 mg or 1600 mg BID loading dose, followed by 800 mg or 600 mg BID or supportive care with a confirmed diagnosis of SARS-CoV-2 by real-time polymerase chain reaction (RT-PCR) from nasopharyngeal swab. Patients were excluded if they were outside the study period, classified as having a mild form of the disease per WHO criteria, or had an incomplete patient file.

Data were manually extracted from electronic health records (EHRs) from a master list and entered into the Research Electronic Data Capture (REDCap) system in a de-identified manner 16 . Data quality during the collection process was ensured to assess consistency and resolve any discrepancies.

Disease severity was defined according to WHO classification, which was adopted by the Saudi MOH protocol 17 . Moderate illness was defined as the development of non-severe pneumonia not requiring supplemental oxygen. Severe illness was defined as fever plus !1 of the following symptoms: respiratory rate !30/min, dyspnea, respiratory distress, SpO2: 93% on room air, PaO 2 /FiO 2 ratio <300 or lung infiltrate >50% of lung field within 24-48 hr. Critical illness was defined by one or more of the following presentations: ARDS, septic shock, altered consciousness, multi-organ failure. The criteria for clinical discharge used in our analyses were those prescribed by protocol by the Saudi MOH for symptomatic and confirmed COVID-19 cases. The initial MOH definition for clinical discharge was that a patient must have passed three days beyond recovery (defined as fever resolution without antipyretic medications and resolution of respiratory symptoms [cough or shortness of breath]), followed by two negative respiratory samples !24 h apart 13 . Later, an updated definition by the WHO (May 2020) was adopted by the Saudi MOH and became the standard protocol 18 , in which clinical discharge can be after fever resolution without antipyretic medications and resolution of respiratory symptoms for at least 3 days, in addition to at least 10 days having passed since symptoms first appeared 13, 18 . Upon review of the discharge dates for patients in our study, we realized that all patients were discharged after 1 June, and therefore that the updated WHO-based definition applied to all discharges.

The primary endpoint was the time to discharge, which served as a surrogate for clinical recovery. Other endpoints were mechanical ventilation progression and mortality. For the mortality outcome, patients on mechanical ventilation at baseline were analyzed as a subgroup. Safety endpoints included adverse events as defined by the Common Terminology Criteria for Adverse Events (CTCAE) 19 .

Assuming that 4% of the subjects in the control group and 3% in the FPV group would experience the outcome of death (based on a relative risk of 0.8 between the supportive group and FPV group) and after applying continuity correction, the study was estimated to require a sample size of 234 for each group (i.e. a total sample size of 468, assuming equal groups) to achieve a power of 80% and a level of significance of 5%. In order to declare the test drug to be superior to the supportive care, a À10% margin of superiority was set (assuming that a smaller proportion is desirable).

Using R Core Team (2020) software (R Foundation for Statistical Computing, Version 4.0.1, Vienna, Austria), continuous variables in the baseline characteristics were presented as means with standard deviations (±SD) and medians with interquartile ranges (IQRs). Student's t-test or Mann-Whitney U test was used for arm comparison. Categorical data were reported as frequencies and percentages and analyzed using either the Chi-square test for nxm tables or Fisher's exact test for 2 Â 2-table group comparisons. For modeling purposes, we quantified the correlation between the WHO severity classification and baseline oxygen therapy using Cram er's V test. We performed the survival analysis using the "survival package" in R 20 . We fitted Kaplan-Meier (KM) models to estimate median time to discharge for live patients and each WHO severity category. Death was considered as a competing risk with discharge events, and therefore, we analyzed the outcome on the live-patient dataset. Semi-parametric Cox proportional hazard models were fitted with clinically and statistically significant covariates to obtain hazard ratios. We used a rule of 10 events for each covariate considered in the Cox models. Proportional hazard assumption was checked using Schoenfeld residuals plots, and nonlinearity for continuous variables was checked using the Martingale residual. The possibility of time-varying coefficients was rigorously tested using the "timereg" package 21 . There were few missing data for few pertinent variables such as body mass index (BMI), renal function, Acute Physiology and Chronic Health Evaluation II (APACHE II). We determined that these variables were missing at random (MAR); therefore, we performed multivariate imputation with the chained equations (MICE) technique using the "mice" package that included the Nelson-Aalen estimator for hazards in the imputation model, as recommended by the literature 22, 23 . We obtained 50 imputations (5 imputed datasets with 10 iterations) and then pooled estimates across datasets. The subgroup analyses for the primary outcomes were performed on the baseline characteristics. In the case of low events, crude odds ratios were estimated with 0.5 correction in the case of zero events in accordance with Deeks and Higgins's recommendations 24 . All statistical inferences were drawn with 95% confidence intervals with p < .05.

The PSM was conducted as a confirmatory analysis. We used the "MatchThem" package on the imputed datasets (50 imputations, 5 datasets with 10 iterations). The use of this package in R was recently illustrated by Pishgar et al. 25 . The matching-procedure specifications followed the nearestneighbor matching within datasets approach (using Rubin's rules) and 1:1 ratio with a caliper of 0.1 and no replacement. Adequate propensity scores distribution for the matched datasets was evaluated using the standardized mean difference (SMD) 26 . SMD values <0.1 indicated adequate balance and fruitful matching. We produced love plots for the SMD distribution ( Figure S1 (A-C) and Table S1 in the supplementary material detailing SMD values for the covariates used) and a mirror-diagram distribution to evaluate propensity scores distribution ( Figure S2 (A-C)). Cox proportional-hazard models were then fitted to the matched datasets to obtain hazards ratios. The variance was estimated using the "Survey" package 27 .

A total of 661 patients were identified and through random selection we screened 480 records of patients admitted to the two hospitals (KFMC and PMAH) between April and August 2020. Reasons for exclusion are illustrated in detail in Figure 1 . The FPV arm had 234 patients and the supportivecare arm had 223 patients, the latter slightly below the target but sufficient for adequate statistical power. Table 1 presents baseline characteristics for each arm. There were statistical differences for renal function, hypertension, criticalillness category, non-invasive and invasive oxygen therapy, and several of the medications received during the hospital stay, showing a relative imbalance between both treatment arms. The majority of patients (86.7%) received the 1600 mg CURRENT MEDICAL RESEARCH AND OPINION 3 BID loading dose with median duration (IQR) of 7 days (6-8 days). The median time (IQR) to start FPV therapy was 2 days (1-3 days) from admission. A moderate association was found between the baseline WHO category and baseline oxygen therapy with a Cram er's V of 0.57.

Live patients on FPV had a median time to discharge of 10 days (95%CI 9-10) compared to 15 days (95%CI 14-16) in the supportive-care arm. The Kaplan-Meier (K-M) curve for discharge events is shown in Figure 2 (A). The adjusted discharge ratio was 1.96 (95%CI 1.56-2.46, p < .001) favoring the FPV arm. Progression to mechanical ventilation was lower in the FPV arm compared to the supportive-care arm with adjusted HR (HR adj ) ¼ 0.10 (95%CI 0.04-0.29, p < .001). Clinical outcomes are detailed in Table 2 , and Supplemental materials ( Figure S3 ) shows K-M curves for specific categories. FPV was not associated with a decrease in hospital mortality, with adjusted HR adj ¼1.56 (95%CI 0.73-3.36, p ¼ .251). The FPV arm was associated with higher mortality in the critical-illness category, with HR ¼ 2.92 (95%CI 1.04-8.20, p ¼ .041); however, when adjusted for age, HR was 2.80 (95%CI 0.99-7.89, p ¼ .051). Table 3 shows detailed results for the mortality outcome. The survival K-M curve displayed in Figure 2 (B) and Supplemental materials ( Figure S4 ) shows K-M curves for specific categories for mortality outcome.

The subgroup analysis for the discharge outcome (live dataset) ( Figure 3 ) shows greater benefit in male compared to female patients (HR ¼ 1. 

The distribution balance for the covariates is presented in Table S3 in the Supplementary section. Further illustration of the balance of covariates is presented in Figure S1 (Love plots), and propensity scores distribution is shown in Figure  S2 (mirror diagram). As shown in 

There were no statistically significant differences in adverse events between the two study arms, except for hyperuricemia 12 (5.1%) and diarrhea 13 (5.6%) (p < .001). See Table 5 for detailed adverse-event results.

To date, the US FDA has approved only RDV to treat COVID-19, despite the WHO's large SOLIDARITY clinical trial, which found RDV to be ineffective in reducing hospital mortality 

Progression to mechanical ventilation, n (%) Cox proportional hazard model with the following variables: age, sex, BMI (body mass index), ethnicity, baseline oxygen therapy, hypertension, CKD-EPI (Chronic Kidney Disease Epidemiology Collaboration), diabetes, intravenous steroid, dyslipidemia, established cardiovascular disease, community-acquired pneumonia regimen (azithromycin and ceftriaxone) and angiotensin-converting enzyme inhibitor/receptor blockers (ACEIs/ARBs).

In the case of critical cases, due to sample size the adjusted model included only age as a covariate. The missing values were determined to be missing at random (namely CKD-EPI and BMI) and were estimated using a multivariate imputation for the chained equations (MICE) technique that included the Nelson-Aalen estimator for hazards in the imputation model. c In the case of low events in particular categories; crude odds per Deeks and Higgins's recommendations 24 . The adjusted Cox hazard model included age, hypertension, and non-invasive baseline oxygen.

Abbreviations. CI, confidence interval; FPV, favipiravir; SC, supportive care. and progression to mechanical ventilation 7, 28 . Recently, the FDA granted approval for a new monoclonal antibody therapy for COVID-19 29 . Data supporting COVID-19 treatment modalities beyond antivirals, such as tocilizumab and convalescent plasma, remain conflicting 30, 31 . However, dexamethasone was found to be associated with lower mortality in specific COVID-19 populations 32 . Novel antiviral drugs with activity against COVID-19 are under investigation, and molecular-docking studies of existing antiviral agents (e.g. protease inhibitors and FPV) have demonstrated that modified antiviral analogues have high drug scores and similarities to parent drugs with improved drug properties 33 . Treatment guidelines published by the National Institutes of Health (NIH) do not include FPV as a treatment option in patients with confirmed COVID-19 infection, possibly due to its lack of availability in the US 34 . By contrast, other institutions and/or country-specific guidelines list FPV as one of the treatment options for patients with confirmed mild-to-severe COVID-19 13, 35 . This study examined whether FPV improved COVID-19 patients' clinical outcomes in terms of time to discharge and mortality compared to supportive care. To the best of our knowledge, the present study is the first of its kind to compare FPV to supportive care in patients with COVID-19 across the spectrum of severity definition in a real-world setting. The overall analysis showed that FPV treatment for a median duration of 7 days resulted in a faster discharge rate in the overall analysis (median 10 days versus 15 and discharge ratio: 1.96). Although this statistical significance was not observed, a trend was seen following the adjustment of confounder variables in the moderate subgroup (median 9 days versus 11; discharge ratio: 1.48, 95%CI 0.97-2.27 (Table 2 ; Supplementary Figure S3(B) ). Interestingly, the discharge ratio remained statistically significant in both severe and critical subgroups (Table 2; Figure S3 ). Chen and colleagues conducted a multicenter, open-label randomized trial that evaluated clinical recovery in the FPV arm compared to umifenovir and found no statistical difference in either the severe or critical subgroups 11 . By contrast, Cai et al. conducted an open-label randomized trial to compare FPV in combination with interferon-alpha to lopinavir/ritonavir in combination with interferon-alpha. The FPV plus interferonalpha combination was associated with rapid clearance of SARS-CoV-2, evidenced by improvements in chest imaging and rapid viral clearance compared with the control arm 10 . Our results align with those of Cai et al. where FPV was associated with a higher discharge ratio. The benefit of FPV was also evident in regard to reduced progression to mechanical ventilation in the overall analysis and the severe COVID-19 category ( Figure S3(E) ; Table 2 ). The lack of statistical significance in the moderate group could be explained by the small sample size.

Although not statistically significant, the association with mortality, mainly in the critically ill group, requires further studyconsidering also that no previous studies have explored the possible effect of FPV on mortality when compared to supportive care or active control 10, 11 . The subgroup analysis suggested that mechanically ventilated patients Hazard ratios >1 indicate an association of harm in the favipiravir arm; <1 is associated with benefit and 1 has no association. The significance of this association is interpreted by considering the confidence interval.

The unadjusted model is a simple Cox model with arm as a covariate. Odds ratios was estimated in the case of low zero events for the moderate category (calculated according to Altman, 1991).

b Cox proportional hazard model for baseline oxygen therapy for the overall analysis. For the overall model, we adjusted for age, Acute Physiology and Chronic Health Evaluation (APACHE II) score, and steroid use. The severe category was adjusted for age and APACHE II score only. The critical illness category was adjusted for age only. Age was an independent risk factor for mortality in the overall analysis of severe and critical illness.

Missing data computed using the multivariate imputation by chained equations (MICE) technique that included the Nelson-Aalen estimator for hazards in the imputation model. Crude odds ratios were estimated in the case of low zero events for moderate category, per Deeks and Higgins's recommendations 24 .

Abbreviations. FPV, favipiravir; SC, supportive care. were at higher risk of dying in the FPV group (Figure 4 ). This result prompted us to investigate the mortality outcome in this subset of patients. The Cox regression model indicated that age was the only independent risk factor associated with this possible increase in mortality in mechanically ventilated patients. While reassuring, given the retrospective and non-randomized nature of our study, we implemented the PSM procedure to achieve balance in the data and to compensate for the non-randomized study design. This analysis confirmed the benefits seen in the discharge ratio; however, it presented a possible mortality signal, as the 95%CI of the HR was completely outside the (admittedly arbitrary) superiority margin of 0.90, thus signaling a possible association of FPV treatment with mortality. There could be a number of explanations for the mortality results. Although FPV is well-distributed to the lungs, it is a mild inhibitor of RdRP compared to RDV (a closely related medication in the same class); therefore, its efficacy in the severe/critical subgroups may be limited 36 . Yet RDV efficacy was also examined recently by Beigel and colleagues' RCT in Table 4 . Clinical outcomes (discharge, progression to mechanical ventilation, and mortality) using propensity score matching procedure. COVID-19 patients. Similar to our findings, RDV was associated with shorter time to recovery but did not result in any mortality benefits 6 . We believe that the current sample size is sufficient; however, we urge scientific communities to supplement our study with additional data and, possibly, a larger sample size. Rattanaumpawan and colleagues recently published a retrospective study to determine factors that could predict poor outcomes in FPV recipients. Older age and lower loading dose were associated with poor outcomes in the logistic regression multivariable model 37 . No difference in regard to age was found in our subgroup analysis (Figure 4) . Also, since this is a real-world study, it revealed the prescribing pattern in our institutions with a majority of the FPV patients (87%) receiving 1600 mg BID loading dose rather than 1800 mg BID. The initial Saudi MOH treatment protocol (or guidelines) for FPV recommended a loading dose of 1600 mg BID, followed by 600 mg BID; however, the updated protocol recommends a loading dose of 1800 mg BID followed by 800 mg BID. Whether higher doses are associated with better outcomes in COVID-19 patients remains undetermined. Evidence from a recently published pharmacokinetic study indicates that when FPV was dosed at 1600 mg loading dose followed by 600 mg BID, the trough concentration in most samples was lower the limit of quantification (1 mg/mL) and half the maximal effective concentration (9.7 mg/mL) in critically ill patients compared to healthy volunteers 38 . Due to the limited number of patients receiving 1800 mg BID, comparing doses was not possible. Despite this limitation, our study confirmed that FPV was associated with a faster discharge rate and less progression to mechanical ventilation compared to supportive care. In addition, the subgroup analysis provided an insight into which group of patients benefited from FPV therapy. Similar to RDV trials, patients on oxygen therapy benefited more from FPV therapy compared to nooxygen therapy or mechanical ventilation 6 . This benefit was also independent of steroid use since the Cox proportional hazard models included steroid use as a confounder. Moreover, as we included these covariates (steroid and baseline oxygen therapy) as matching covariates, we eliminated the differences between the two groups in the propensity scores sensitivity analysis.

It can be argued that patients started on FPV treatment were discharged sooner regardless of whether they met the strict clinical criteria for recovery, e.g. chest X-ray clearance or viral clearance. For this study, we were limited by the clinical discharge definition as described in the Saudi MOH guidelines. In the study by Cai et al., viral clearance was the discharge criterion 10 . This was defined as "the presence of two consecutive negative results with quantitative polymerase chain reaction (qPCR) detection over an interval of 24 h." In that study, the median time (IQR) for viral clearance was 4 days (95% CI 2.5-9 days) for FPV treatment versus 10 days (95% CI 9-10) in our study. This could be due to several factors. First, our patients had more severe cases of COVID-19 than those in the study by Cai et al., which excluded patients with severe COVID-19 disease. Second, in that study, FPV was combined with interferon (IFN)-alfa, whereas most of our patients were treated with FPV as the main therapy for COVID-19. Lastly, in cases of a lack of availability of PCR testing, providers had to wait 3 days before a decision could be made, and one may argue that this prolonged discharge rates. However, our understanding of the virus evolved during the pandemic, and the WHO's latest recommendations (updated May 2020) added the criterion that patients could be discharged 3 days after symptom resolution in addition to 10 days having elapsed since initial presentation of symptoms 18 . This is independent of the PCR results, as later evidence showed that some patients may have prolonged shedding of the virus, thus delaying the discharge decision. The delay in discharge can be consequential as it affects access to the healthcare system and the individual's wellbeing. Moreover, repeated PCR can be challenging due to the unavailability of testing capacity in many parts of the world. We believe that the defined discharge criteria in our case provided an adequate surrogate marker for clinical recovery because patients in Saudi Arabia must meet the discharge criteria based on the Saudi MOH discharge protocol (similar to the WHO's updated recommendation). Moreover, after a review of our data for discharge dates we found out that all the included patients were discharged after June 2020. Therefore, we believe that many providers may have decided based on the second updated criterion as it provides several advantages to the healthcare system as a whole.

Arguably the cellular damage and the fluid accumulation within the lungs induced by the cytokine storm may not be mitigated by FPV when started late during the course of therapy (i.e. >48 h of symptoms onset). In our study, the median time to initiation was 2 days (IQR: 1-3 day) from hospital admission. Given the retrospective nature of our study, lack of evidence from prospective studies, the evolving nature of the disease during the first wave of the pandemic, and inconsistency in the prescriptive pattern, initiating FPV within 48 of symptoms in all patients may not be feasible. In addition, the use of antibiotics differed between the two arms. This may have been due to differences in the change of patients' clinical courses. We have merely described the pattern of usage, as antibiotics are not indicated for COVID-19 disease. We have not attempted to include these variables in the propensity score matching procedure as they can be considered as post-treatment covariates 39 . Moreover, the inclusion of variables as multivariable Cox models may present an unnecessary statistical challenge because it could easily violate the basic assumption of hazards proportionality 40 .

Due to the retrospective nature of our study, the pandemic, and the need for a rapid response with effective treatment options, clinicians initially utilized FPV without checking uric acid levels. They began checking levels as their knowledge of this medication expanded with more data being introduced in the literature concerning the safety of this antiviral agent. Furthermore, we have not collected some baseline data, such as C-reactive protein (CRP), procalcitonin, D-dimer, aspartate transaminase (AST), alanine transaminase (ALT), and lymphocyte count, that may have affected treatment allocation.

Another limitation was the imbalance between the two arms in regards to variables in the baseline characteristics. To address this limitation, we followed two analytical approaches (multivariable Cox models and PSM procedures) to account for clinically important confounders. Lastly, missing data were addressed with the powerful MICE procedure and could have affected the discharge outcome only. However, the outcomes in the complete case analysis and our MICE estimation with missing data showed a negligible numerical difference that did not materially change the outcome.

This independent study suggests that FPV treatment is welltolerated and is associated with faster discharge compared to supportive care across the COVID-19 severity spectrum. However, FPV treatment is not protective against mortality, and age was the only independent risk factor for mortality in mechanically ventilated patients. To confirm or refute our findings, a large RCT is needed.

No funding or sponsorship was received for this study or for the publication of this article.

The authors have no relevant affiliations or financial involvement with any organization or entity with a financial interest in or financial conflict with the subject matter or materials discussed in the manuscript. This includes employment, consultancies, honoraria, stock ownership or options, expert testimony, grants or patents received or pending, or royalties. Peer reviewers on this manuscript have no relevant financial or other relationships to disclose.

COVID-19, SARS and MERS: are they closely related? Clinical microbiology and infection: the official publication of the European Society of Clinical Microbiology and Infectious Diseases

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

Olfactory and gustatory function impairment in COVID-19 patients: Italian objective multicenter-study

Clinical features of patients infected with 2019 novel coronavirus in Wuhan

The journey of remdesivir: from Ebola to COVID-19

Remdesivir for the treatment of Covid-19 -preliminary report. Reply

FDA approves first treatment for COVID-19

T-705), a novel viral RNA polymerase inhibitor

Experimental treatment with favipiravir for ebola virus disease (the JIKI Trial): a historically controlled, single-arm proof-of-concept trial in Guinea

Experimental treatment with favipiravir for COVID-19: an open-label control study

Favipiravir versus arbidol for COVID-19: a randomized clinical trial

A review of the safety of favipiravir -a potential treatment in the COVID-19 pandemic?

Saudi MoH Protocol for Patients Suspected of/Confirmed with COVID-19

The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) Statement: guidelines for reporting observational studies

The retrospective chart review: important methodological considerations

The REDCap consortium: building an international community of software platform partners

Coronavirus Disease COVID-19 Guidelines, v1.3

World Health Organization. Criteria for releasing COVID-19 patients from isolation

Department of Health and Human Services & National Institutes of Health. Common Terminology Criteria for Adverse Events (CTCAE) Version 5.0

A package for survival analysis in R. 2020

Analyzing competing risk data using the R timereg package

Multiple imputation for missing data in epidemiological and clinical research: potential and pitfalls

Multivariate imputation by chained equations in R

Cochrane handbook for systematic reviews of interventions

MatchThem: matching and weighting after multiple imputation

A comparison of two methods of estimating propensity scores after multiple imputation

Analysis of complex survey samples

Repurposed antiviral drugs for COVID-19 -interim WHO Solidarity trial results

COVID-19) update: FDA authorizes monoclonal antibodies for treatment of COVID-19

Treatment of COVID-19 with convalescent plasma: lessons from past coronavirus outbreaks

Time to reassess Tocilizumab's role in COVID-19 pneumonia

Dexamethasone in hospitalized patients with COVID-19

Combination of QSAR, molecular docking, molecular dynamic simulation and MM-PBSA: analogues of lopinavir and favipiravir as potential drug candidates against COVID-19

COVID-19) treatment guidelines

Brigham and women's. Brigham and women's hospital inpatient COVID 19 infectious diseases treatment guidelines

Tissue distributions of antiviral drugs affect their capabilities of reducing viral loads in COVID-19 treatment

Realworld experience with favipiravir for treatment of COVID-19 in Thailand: results from a multi-center observational study

Pharmacokinetics of favipiravir in critically ill patients with COVID-19

Propensity scores in cardiovascular research

Survival analysis part II: multivariate data analysis-an introduction to concepts and methods

The authors would like to thank the staff at the Research Center at King Fahad Medical City, Riyadh, for their valuable technical support provided for the manuscript. This publication was supported by the Deanship of Scientific Research at Prince Sattam bin Abdulaziz University. The authors would also like to express their gratitude to the staff at the Ministry of Health in Saudi Arabia for their support.

All data generated or analyzed for this study are included in this published article (and its Supplementary Information files) .

The study was approved by the Institutional Review Boards at KFMC and PMAH (IRB 20-477E, July 2020). Informed consent was waived as this study was considered exempt.