EXCLI J EXCLI Journal 1611-2156 Leibniz Research Centre for Working Environment and Human Factors 2021-3413 10.17179/excli2021-3413 Doc199 Original article All-cause mortality among patients treated with repurposed antivirals and antibiotics for COVID-19 in Mexico City: A real-world observational study Mancilla-Galindo Javier 1 2 García-Méndez Jorge Óscar 2 3 4 Márquez-Sánchez Jessica 5 Reyes-Casarrubias Rodrigo Estefano 2 3 Aguirre-Aguilar Eduardo 6 Rocha-González Héctor Isaac 7 Kammar-García Ashuin http://orcid.org/0000-0002-3875-0945 * 6 7 Unidad de Investigación UNAM-INC, Instituto Nacional de Cardiología Ignacio Chávez, Mexico City, Mexico Facultad de Medicina, Universidad Nacional Autónoma de México, Mexico City, Mexico Departamento de Posgrado y Educación Médica Continua, Instituto Nacional de Cancerología, Mexico City, Mexico Departamento de Infectología, Fundación Clínica Médica Sur, Mexico City, Mexico Departamento de Infectología, Instituto Nacional de Pediatría, Mexico City, Mexico Departamento de Atención Institucional Continua y Urgencias, Instituto Nacional de Ciencias Médicas y Nutrición Salvador Zubirán, Mexico City, Mexico Sección de Estudios de Posgrado e Investigación, Escuela Superior de Medicina, Instituto Politécnico Nacional, Mexico City, Mexico *To whom correspondence should be addressed: Ashuin Kammar-García, Departamento de Atención Institucional Continua y Urgencias, Instituto Nacional de Ciencias Médicas y Nutrición Salvador Zubirán, Mexico City, Mexico, Vasco de Quiroga 15, Tlalpan, Col. Belisario Domínguez Sección XVI, CP 14080 Ciudad de México, México; Tel: +52 5554870900, ext. 5010, E-mail: kammar_nutrition@hotmail.com 04 02 2021 2021 20 199 222 19 01 2021 01 02 2021 Copyright © 2021 Mancilla-Galindo et al. 2021

This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (http://creativecommons.org/licenses/by/4.0/) You are free to copy, distribute and transmit the work, provided the original author and source are credited.

This article is available from https://www.excli.de/vol20/excli2021-3413.pdf

The aim of our study was to evaluate all-cause mortality risk in patients with laboratory-confirmed COVID-19 in Mexico City treated with repurposed antivirals and antibiotics. This real-world retrospective cohort study contemplated 395,343 patients evaluated for suspected COVID-19 between February 24 and September 14, 2020 in 688 primary-to-tertiary medical units in Mexico City. Patients were included with a positive RT-PCR for SARS-CoV-2; those receiving unspecified antivirals, excluded; and groups of antivirals prescribed in < 30 patients, eliminated. Survival and mortality risks were determined for patients receiving antivirals, antibiotics, both, or none. We assessed the effect of early (<2 days) versus late (>2 days) use of antivirals on mortality in a sub-cohort of patients. Multivariable adjustment, propensity score matching, generalized estimating equations, and calculation of E-values were performed to limit confounding. 136,855 patients were analyzed; mean age 44.2 (SD:16.8) years; 51.3 % were men. 16.6 % received antivirals (3 %), antibiotics (10 %), or both (3.6 %). Antivirals studied were Oseltamivir (n=8414), Amantadine (n=319), Lopinavir-Ritonavir (n=100), Rimantadine (n=61), Zanamivir (n=39), and Acyclovir (n=36). Survival with antivirals (73.7 %, p<0.0001) and antibiotics (85.8 %, p<0.0001) was lower than no antiviral/antibiotic (93.6 %). After multivariable adjustment, increased risk of death occurred with antivirals (HR=1.72, 95 % CI: 1.61-1.84) in ambulatory (HR=4.7, 95 % CI: 3.94-5.62) and non-critical (HR=2.03, 95 % CI: 1.86-2.21) patients. Oseltamivir increased mortality risk in the general population (HR=1.72, 95 % CI: 1.61-1.84), ambulatory (HR=4.79, 95 % CI: 4.01-5.75), non-critical (HR=2.05, 95 % CI: 1.88-2.23), and pregnancy (HR=8.35, 95 % CI: 1.77-39.30); as well as hospitalized (HR=1.13, 95 % CI: 1.01-1.26) and critical patients (HR=1.22, 95 % CI: 1.05-1.43) after propensity score-matching. Early versus late oseltamivir did not modify the risk. Antibiotics were a risk factor in general population (HR=1.13, 95 % CI: 1.08-1.19) and pediatrics (HR=4.22, 95 % CI: 2.01-8.86), but a protective factor in hospitalized (HR=0.81, 95 % CI: 0.77-0.86) and critical patients (HR=0.67, 95 % CI: 0.63-0.72). No significant benefit for repurposed antivirals was observed; oseltamivir was associated with increased mortality. Antibiotics increased mortality risk in the general population but may increase survival in hospitalized and critical patients.

SARS-CoV-2 COVID-19 oseltamivir antibiotics pharmacoepidemiology
Introduction

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the etiologic agent of the coronavirus disease (COVID-19) pandemic, one of the most devastating infectious diseases of this century. Non-pharmacological interventions are the most effective means of limiting the impact of COVID-19 (Islam et al., 2020[13]). However, several countries have not been able to contain the disease (Sachs et al., 2020[32]).

One of the main strategies for finding ways to combat COVID-19 is drug repurposing since developing novel antivirals against SARS-CoV-2 may be protracted (Saul and Einav 2020[33]). Repurposing existing antivirals is attractive due to their relative safeness and potential anti-SARS-CoV-2 mechanisms (Senanayake, 2020[34]). Neuraminidase inhibitors (i.e. oseltamivir, zanamivir) and HIV protease inhibitors (i.e. lopinavir-ritonavir) have been hypothesized to inhibit SARS-CoV-2 proteases involved in the degradation of polyproteins that control viral replication (Muralidharan et al., 2020[26]). Adamantanes (i.e. amantadine, rimantadine) are thought to disrupt lysosomal trafficking, thereby impeding the release of SARS-CoV-2 ribonucleic acid (RNA) into the cell (Araújo et al., 2020[2]), and by inhibiting conductance of the envelope (E) protein (Aranda Abreu et al., 2020[1]). Acyclovir, a nucleotide analog antiviral, was found as a candidate drug for COVID-19 by potentially counteracting gene expression changes observed after SARS-CoV-2 infection (Li and Yang, 2020[16]).

Up to December 8, 2020 there were 413 registered studies to test antivirals for COVID-19, of which 400 were still active (NIHR, 2020[28]). Most trials on the World Health Organization (WHO) platform were for lopinavir/ritonavir (156), favipiravir (62), remdesivir (52), oseltamivir (19), and ribavirin (16) (WHO, 2020[39]). Other common antivirals are not tested for COVID-19, but may be widely used in the general population and in hospitals, as most antivirals are not recommended in practice guidelines due to lack of evidence (Bhimraj et al., 2020[4]; NIH, 2020[27]), others advise against most (Mexican Secretariat of Health, 2020[24]; WHO, 2020[38]), or recommend oseltamivir empirically during the influenza season (Mexican Secretariat of Health, 2020[23]) and when coinfection exists (Falavigna et al., 2020[8]).

The WHO recommends early empiric antibiotic therapy with continuous reassessment for de-escalation in patients with severe COVID-19, but not for mild-to-moderate disease (WHO, 2020[38]). However, questions have arisen on the effectiveness of antibiotics due to the viral origin of the disease, whilst indiscriminate use of antibiotics could aggravate antimicrobial resistance.

In this population-based study, we hypothesized that repurposed antivirals and antibiotics could be widely used in real-world settings. To determine their impact on mortality, we studied survival and all-cause mortality risk in patients with laboratory-confirmed COVID-19 in Mexico City receiving these drugs in both ambulatory and in-hospital settings.

MethodsStudy design

We conducted a real-world multicenter retrospective cohort study in patients who received medical attention for suspected COVID-19 in any of the 688 registered and accredited COVID-19 medical units in Mexico City, to evaluate all-cause mortality (main outcome) in those receiving antivirals, antibiotics, both, or none (exposition groups).

We considered 395,343 patients for eligibility evaluated for COVID-19 in 688 medical units (primary-to-tertiary care) between February 24, 2020 and September 14, 2020. All patients with a positive RT-PCR for SARS-CoV-2 were included to maximize the power and generalizability. Patients treated with unspecified antivirals were excluded. To perform reliable analyses, a cut-off value of 30 patients receiving the same antiviral was set and groups of antivirals with <30 patients were eliminated.

Source of data

We used the COVID-19 open dataset available on Mexico City Government's Open Data platform (Government of Mexico City, 2020[9]), collected and updated daily by the Secretariat of Health of Mexico City. Patients meeting criteria of suspected COVID-19 case have been included in this dataset starting on February 24, 2020 when the first suspected cases arrived in Mexico.

Criteria for suspected COVID-19 case in Mexico included having at least two of three signs/symptoms (cough, fever, or headache) plus at least one other (dyspnea, arthralgias, myalgias, sore throat, rhinorrhea, conjunctivitis, or chest pain) in the last 7 days. This operational definition was changed on August 24, 2020 to increase sensitivity (Directorate General of Epidemiology of Mexico, 2020[7]): at least one of four signs/symptoms (cough, fever, dyspnea, or headache), plus at least one other (myalgias, arthralgias, sore throat, chills, chest pain, rhinorrhea, anosmia, dysgeusia, or conjunctivitis) in the last 10 days.

For epidemiologic purposes, two strategies are outlined in the National COVID-19 Epidemiologic Surveillance Plan (Directorate General of Epidemiology of Mexico, 2020[6]): 1. testing of 10 % of ambulatory patients with mild symptoms of respiratory disease and 100 % of patients with respiratory distress at evaluation in monitoring units of viral respiratory disease (USMER, for its acronym in Spanish), and 2. testing 100 % of patients who meet diagnostic criteria of Severe Acute Respiratory Infection (defined as shortness of breath, temperature ≥38 °C, cough, and ≥1 of the following: chest pain, tachypnea, or acute respiratory distress syndrome) in non-USMER units.

Upon evaluating a patient suspected of having COVID-19, healthcare professionals are required to fill out a format (Supplementary Appendix) containing demographic, clinical, epidemiological, and treatment variables, later complemented with follow-up by accredited hospital epidemiologists (inpatients) and healthcare professionals in primary care units (ambulatory patients). For ambulatory patients, follow-up is performed daily for a minimum of 7 days and patients are considered recovered 14 days after the onset of symptoms if alive and not hospitalized. For hospitalized patients, follow-up is done daily until death or discharge; follow-up time for patients discharged from hospital is highly variable since no consensus or requirements by authorities exist but may extend from 14 days to 3-6 months after discharge. Duration of follow-up for each patient is not provided in the dataset and cannot be calculated.

For every medical unit there is only one responsible authority who ultimately uploads data into the Respiratory Diseases Epidemiologic Surveillance System and is accountable for accuracy. Results of diagnostic RT-PCR for SARS-CoV-2 are directly uploaded by the diagnostic facility; accreditation of diagnostic procedures by the Mexican Institute of Diagnostics and Epidemiological Reference is required to upload results. Reporting of all deaths of COVID-19 suspected or confirmed cases is obligatory and must be done in the first 48 hours after occurrence; in cases of deaths occurring in patients who had completed follow-up, registries are matched to death certificates and updated. There have been concerns that patients tested more than once may be duplicated. Since no variables that could lead to identification of patients are released, we searched for patients with identical demographic variables and only one registry was kept.

To determine whether prescription of oseltamivir occurred for cases of COVID-19 or influenza during the pandemic period, dates of prescription of oseltamivir and dates of admission of patients with laboratory-confirmed COVID-19 were obtained from the previously mentioned COVID-19 dataset and grouped according to epidemiological weeks. Data of patients who had a positive test for influenza in Mexico City for every epidemiological week were obtained from the Directorate General of Epidemiology's Weekly Surveillance Reports of Influenza (Directorate General of Epidemiology of Mexico, 2020[7]).

Management of variables

All categorical variables were classified as dummy variables (present/absent). Polytomous variables were created from frequencies of use of antivirals and antibiotics (no antiviral/antibiotic, antiviral only, antibiotic only, and antiviral plus antibiotic), type of antiviral with ≥30 patients, and the combination of every individual antiviral with antibiotics. These were considered as the exposition groups. Special populations for subgroup analyses were defined as: children and adolescents (<18 years), pregnancy, puerperium, and non-pregnant/puerperal adults (≥18 years). Further subgroups included ambulatory and hospitalized patients, as well as patients admitted to ICU and those requiring IMV. A variable of critical patients was built by grouping patients admitted to ICU and/or requiring IMV, whereas non-critical patients did not meet any of both.

Since it has been hypothesized that early use of antivirals for COVID-19 could diminish hospitalization rate (Benlloch et al., 2020[3]) and detain disease progression (Lipsitch et al., 2020[17]), thereby decreasing mortality, we distinguished early (≤2 days from symptom onset to initiation of antivirals) from late (>2 days) use of antivirals, and studied their relation to hospitalization rates and mortality; only patients who received antivirals before being evaluated in an accredited COVID-19 unit were included for this analysis since dates of initiation of antivirals are only collected for such patients. Therefore, we studied these patients as a different sub-cohort.

Occupations were grouped as follows: technical services (laborers), education (students and teachers), healthcare (dentists, nurses, diagnostic laboratorian, physicians, and other healthcare workers), agricultural activities (peasants), commerce (drivers, informal commerce, employees, and business people), unemployed, stay-at-home (stay-at-home parents and retired/pensioners), and other occupations (others, and other professions).

Variables for adjustment of models were: sex, age, indigenous self-identification, diabetes, COPD, immunosuppression, hypertension, HIV/AIDS, cardiovascular disease, obesity, CKD, smokers, unemployed, time from symptom onset to medical attention, fever, cough, sore throat, shortness of breath, irritability, diarrhea, chest pain, chills, headache, myalgias, arthralgias, abrupt deterioration, rhinorrhea, polypnea, vomit, abdominal pain, conjunctivitis, cyanosis, and sudden onset of symptoms.

Statistical analysis

Descriptive data were calculated and are provided as frequencies, percentages, mean with standard deviation (SD) or median with interquartile range (IQR). Qualitative comparisons were made with χ2 or Fisher´s exact test. Independent-samples t-test and ANOVA were used for quantitative comparisons. Survival was calculated for all treatment groups (antiviral only, antibiotic only, antiviral plus antibiotic, and no antiviral/antibiotic) and specific antivirals (acyclovir, amantadine, lopinavir-ritonavir, oseltamivir, rimantadine, and zanamivir) alone or combined with antibiotics; survival curves were created for the general population, ambulatory, hospitalized, non-critical, and critical patients. Survival between groups receiving distinct treatments were compared through the Log-Rank test against patients not receiving antivirals/antibiotics. Cox regression models were applied for the general population, ambulatory, hospitalized, non-critical, and critical patients to determine mortality risk in patients receiving any treatment compared with no antivirals/antibiotics (reference). Resulting hazard ratios (HR) were adjusted for demographic and clinical variables considered as risk factors in the univariate analysis for every group; all variables with p<0.1 were included in the final model using the Enter method.

To account for multicenter variability, adjusted risk was calculated through generalized estimating equations (GEE), setting the medical unit with the lowest case-fatality rate (CFR) and the highest number of patients for every subgroup as the reference value. Further subgroup survival analyses and multivariable Cox regression models were applied for special populations (children and adolescents, pregnancy, puerperium, and non-pregnant/puerperal adults), invasive mechanical ventilation (IMV), and intensive care unit (ICU). To quantify the minimal association strength of an unmeasured confounding factor that could reduce the risk conferred by exposures in our study, E-values were calculated for the point estimate and lower limit of the confidence interval.

To reduce potential confounding and selection bias, we applied a propensity score analysis. Propensity scores were calculated with a logistic regression model adjusted for: sex, age, signs and symptoms (fever, dyspnea, arthralgias, myalgias), and comorbidities (hypertension, diabetes, chronic kidney disease, cardiovascular disease, obesity, and immunosuppression). A 1:1 matching was performed through the nearest-neighbor algorithm matching. Density functions of treated patients against controls before and after matching were graphed to determine appropriate matching. Subsequently, multivariable Cox regression models were applied to determine mortality risk in patients treated with antivirals, antibiotics, and oseltamivir versus a control group.

A two-sided p value <0.05 was used to define statistical significance. Analyses and figures were created with SPSS software v.21 and GraphPad Prism v.8.0.1. Propensity score matching and analyses were done with R software v.3.6.2 and the 'MatchIt' package.

Results

No duplicated registries were found. After selection of eligible participants (Figure 1(Fig. 1)), 136,855 patients from all 688 medical units were analyzed. 97.83 % (n=133,887) were residents of the Mexico City Metropolitan Area, conformed by 17 municipalities of Mexico City (83.29 %, n=111,768), and 60 municipalities (16.71 %, n=22,119) of the State of Mexico. The remaining 2.17 % (n=2,968) sought medical attention from all other 30 states of the republic.

Of all patients, 10.0 % (n=13,743) received antibiotics only; 3.0 % (n=4,044), antivirals only; 3.6 % (n=4,925), antivirals plus antibiotics, and 83.4 % (n=114,143), none (Table 1(Tab. 1)). More symptomatic ambulatory patients received antivirals and antibiotics more frequently (Supplementary Table 1); hospitalized patients with more signs/ symptoms had greater use of antivirals, but less antibiotics (Supplementary Table 2).

Baseline and follow-up characteristics of survivors (91.47 %, n=136,855) and non-survivors (8.53 %, n=11,679) are shown in Supplementary Table 3. CFR in special populations were: 8.92 % (95 % CI: 8.76-9.07 %), for non-pregnant/puerperal adults; 1.72 % (95 % CI: 0.66-2.77), pregnancy; 0.97 % (95 % CI: 0.00-2.90), puerperium; and 0.69 % (95 % CI: 0.48-0.90), children and adolescents. Of all deaths, 92.7 % (95 % CI: 92.2-93.2) and 99.6 % (95 % CI: 99.5-99.7) occurred by day 28 and 56, respectively.

Patients treated only with antivirals had a lower survival rate than those not receiving antivirals/antibiotics in the general population (Figure 2a(Fig. 2)), ambulatory (Figure 2b(Fig. 2)), hospitalized (Figure 2c(Fig. 2)), non-critical (Figure 3a(Fig. 3)), critical (Figure 3c(Fig. 3)), IMV (Supplementary Table 4), ICU (Supplementary Table 5) and non-pregnant/puerperal adults (Supplementary Table 6); for children and adolescents (Supplementary Table 7) and pregnancy (Supplementary Table 8) differences in survival were not significant; there were not enough events for analysis in puerperal women. Increased survival with only antibiotics was observed in hospitalized, critical, and IMV, whereas decreased survival occurred in the general population, non-pregnant/puerperal adults, ambulatory, non-critical, ICU, and children and adolescents; there were no differences for pregnancy. Antivirals plus antibiotics resulted in decreased survival in the general population, ambulatory, non-critical, non-pregnant/puerperal adults, children and adolescents, pregnancy, and ICU; increased survival, in hospitalized; and no differences, in critical and IMV groups.

Decreased survival with oseltamivir was observed in the general population (Figure 2d(Fig. 2)), ambulatory (Figure 2e(Fig. 2)), non-critical (Figure 3b(Fig. 3)), ICU (Supplementary Table 5), non-pregnant/puerperal adults (Supplementary Table 6), children and adolescents (Supplementary Table 7), and pregnancy (Supplementary Table 8); no differences in survival occurred in hospitalized (Figure 2f(Fig. 2)), critical (Figure 3d(Fig. 3)), and IMV (Supplementary Table 4). Survival rates for amantadine, zanamivir, rimantadine, acyclovir, and lopinavir-ritonavir are shown in the same figures and tables as oseltamivir.

Unadjusted (Supplementary Table 9) and adjusted (Table 2(Tab. 2)) risk of death for the general population, ambulatory, hospitalized, non-critical and critical patients, as well as for other subgroups (Supplementary Tables 10-14) were calculated. E-values for statistically significant risk groups are provided in Supplementary Tables 15-16. After adjusting for center through GEE, we found no statistically significant variability in mortality risk for antivirals, antibiotics, or both in all groups. Oseltamivir presented variability in hospitalized and critical patients, with the largest increases in risk occurring in public hospitals. However, center was not a modifying risk factor after logistic regression analysis in hospitalized and critical patients.

After matching patients receiving antivirals (Supplementary Table 17), oseltamivir (Supplementary Table 18), or antibiotics (Supplementary Table 19) to controls, mortality risks were similar to those in the unmatched cohort (Table 3(Tab. 3)). However, contrary to the main analyses, antivirals in hospitalized patients, as well as oseltamivir in hospitalized and critical patients were a risk factor for death. Antibiotics were a protective factor in hospitalized and critical patients, but not in ambulatory and non-critical patients. Density functions before and after matching are shown in (Supplementary Figures 1-3).

Of all 8,969 patients receiving antivirals, 10 % (n=903) received antivirals before evaluation in COVID-19 accredited units; their baseline and follow-up characteristics are available in Supplementary Table 20. 25.2 % (n=227) were admitted to hospital. Median time from symptom onset to initiation of antivirals was 1 day (IQR:0-4) for both ambulatory and hospitalized patients; time from symptom onset to ambulatory care in accredited units was 5 days (IQR:3-8) and 6 days (IQR:4-9) for hospitalization. Time from initiation of antivirals to hospitalization was 3 days (IQR:0-6). Time-to-initiation of antivirals and time-to-hospitalization for specific antivirals are shown in Supplementary Figure 4.

Early (≤2 days) and late (>2 days) initiation of antivirals occurred in 64.2 % (n=580) and 35.8 % (n=323) patients, respectively. Overall survival in early (91.3 %) and late (88.9 %) groups was not different (p=0.2). Survival for early/late use of antivirals is shown in Supplementary Table 21. Oseltamivir was associated with increased risk of death in both early (HR=3.00, 95 % CI: 2.14-4.20) and late (HR=2.99, 95 % CI: 1.83-4.89) groups, as well as late use of lopinavir-ritonavir (HR=9.9, 95 % CI: 2.49-39.83); all other early/late antivirals did not reach statistical significance. There were no differences in hospitalization rates between early and late groups for every antiviral (Supplementary Figure 5).

Prescription of oseltamivir during the pandemic period was greater than confirmed influenza cases (Figure 4A(Fig. 4)) and followed a similar pattern to that of weekly new cases of patients with a positive SARS-CoV-2 test (Figure 4B(Fig. 4)).

Discussion

To the best of our knowledge, this is the first observational study evaluating amantadine, rimantadine, zanamivir, and acyclovir for COVID-19; no registered studies to evaluate these drugs exist (WHO, 2020[39]). Only one study has evaluated risk of death for oseltamivir (Liu et al., 2020[18]); lopinavir-ritonavir have been evaluated in clinical trials (Horby et al., 2020[12]).

We hypothesized that antivirals and antibiotics could be widely used in real-world settings as repurposed drugs for COVID-19. Therefore, we studied mortality in laboratory-confirmed COVID-19 patients treated with common antivirals and/or antibiotics in Mexico City. Most patients did not receive them (83.4 %), although a substantial proportion received antivirals alone (3.0 %) or combined with antibiotics (3.6 %) despite national guidelines advising against antivirals out of clinical trials (Mexican Secretariat of Health, 2020[24]). Patients receiving antivirals and antibiotics were overall more symptomatic, suggesting that florid clinical presentations and not evidence may be guiding decision to treat.

Our cohort included mostly middle-aged adults (44.2 [SD:16.8] years), of which non-survivors were younger and with a higher burden of hypertension, diabetes and obesity compared with other countries, findings consistent with previous studies in Mexican patients (Mancilla-Galindo et al., 2020[22], 2021[21]; Vera-Zertuche et al., 2021[37]). Rates of patients with cardiovascular diseases, COPD, asthma, and smokers were low, suggesting that underdiagnosis could be occurring. Only 20-30 % of patients had been vaccinated against influenza during the prior season.

Of patients treated with antivirals, ~10 % (n=903) started antivirals before seeking medical attention in accredited COVID-19 units. Hospitalization rates were not different when antivirals were used early (≤2 days) vs late (>2 days). In Mexico, antibiotics and most antivirals are sold under prescription. However, private pharmacy-associated clinics are a rapidly growing sector not included in our study where physicians tend to have lower experience, qualifications, compliance with regulations, and higher prescription rates (Pérez-Cuevas et al., 2014[29]; Díaz-Portillo et al., 2015[5]; López-Manning and García-Díaz, 2017[20]). Self-medication with amantadine could be occurring since it is a widely available over-the-counter medication. Importantly, the high prescription of oseltamivir in Mexico City during the pandemic period cannot be attributed to influenza outbreaks, since patients receiving oseltamivir largely exceed those with a positive test for influenza and follows a similar weekly pattern to that of newly diagnosed COVID-19 patients.

In one single-center study, oseltamivir was associated with decreased mortality risk in COVID-19-hospitalized patients (HR=0.21; 95 % CI: 0.10-0.43) (Liu et al., 2020[18]). Contrary to Liu et al., we found no benefit for oseltamivir in hospitalized patients (HR=1.07; 95 % CI: 0.99-1.15) which is consistent with studies of oseltamivir for SARS-CoV (HR=0.87; 95 % CI: 0.55-1.38) (Shi et al., 2020[35]). Combination of oseltamivir with antibiotics in hospitalized patients in our study resulted in decreased risk of death (HR=0.92; 95 % CI: 0.87-0.98), which could explain findings by Liu et al.[18], since most patients in their cohort (87.7 %) received antibiotics. Decreased mortality is likely driven by antibiotics since hospitalized patients in our study receiving only antibiotics had lower risk of dying (HR=0.81, 95 % CI: 0.77-0.86) than antibiotics plus oseltamivir.

After multivariable adjustment, oseltamivir was associated with increased mortality in the general population (HR=1.72, 95 % CI: 1.61-1.84), ambulatory (HR=4.79, 95 % CI: 4.01-5.75), non-critical (HR=2.05, 95 % CI: 1.88-2.23), and pregnant (HR=8.35, 95 % CI: 1.77-39.30) patients. Mortality risk was also higher in the cohort of 903 patients with both early (HR=3.00, 95 % CI: 2.14-4.20) and late (HR=2.99, 95 % CI: 1.83-4.89) use of oseltamivir. We performed propensity score matching analyses to further account for potential confounders and found that hospitalized (HR=1.13, 95 % CI: 1.01-1.26) and critical patients (HR:1.22, 95 % CI: 1.05-1.43) also had an increased risk of death when treated with oseltamivir. Further efforts to limit bias were adjustment for center through GEE and calculation of E-values which support a strong association between oseltamivir and increased mortality risk, since potential confounders not accounted for in our study should have substantially large hazard ratios and lower limit confidence intervals (i.e. 9.05 and 7.48, respectively, for ambulatory patients) to explain our findings.

The potentially inhibitory activity of proteases by oseltamivir (Muralidharan et al., 2020[26]) was found to be weak through molecular modeling, while inhibition of SARS-CoV-2 in vitro and reduction of symptoms in hospitalized patients failed (Tan et al., 2020[36]). Furthermore, oseltamivir nearly suppresses the production of pro-inflammatory cytokines by dendritic cells, polymorphonuclear leukocytes, and CD8+ T cells through inhibition of endogenous neuraminidase (sialidase) (Hama, 2016[10]), which could impair the immune response and limit the capacity to eliminate the infection. Use of oseltamivir for infections caused by viruses lacking a neuraminidase gene (i.e. respiratory syncytial virus) decrease viral clearance (Moore et al., 2007[25]), which could explain increased mortality in patients receiving oseltamivir. Additionally, sudden (cardiac arrest, respiratory suppression, hypothermia, and neuropsychiatric events) and delayed-onset (impaired renal function, hyperglycemic events, prolonged QTc interval, immunological impairment, and bleeding) adverse reactions could explain increased risk of death in patients with COVID-19 treated with oseltamivir (Hama, 2016[10]). Antiviral drug-related heart damage is a concern since some antivirals may be cardiotoxic, aggravating myocardial damage caused by SARS-CoV-2 (Long et al., 2020[19]). Renal and psychiatric adverse events in patients receiving oseltamivir were higher compared to placebo in a systematic review (Jefferson et al., 2014[14]).

In the RECOVERY study, there were no differences in mortality risk between hospitalized patients receiving lopinavir-ritonavir vs placebo (HR=1.03, 95 % CI: 0.91-1.17) (Horby et al., 2020[12]), which is consistent with our findings in hospitalized patients. Notably, ambulatory and late (>2 days) use of lopinavir-ritonavir were risk factors for death in our study.

Antibiotics were a risk factor for death in the general population after multivariable analysis but a protective factor in both ambulatory and hospitalized patients. Nonetheless, univariate models showed no overall effect of antibiotics in ambulatory patients; when adjusting only for demographic variables no effect persisted but were protective after adjusting only for clinical variables. This is explained by the fact that more symptomatic patients received antibiotics more often. Supporting this conclusion, no benefit was observed for antibiotics in non-critical patients in both matched and unmatched cohorts.

We observed a potential benefit for use of antibiotics in hospitalized, IMV, and critical patients. Increased survival in these groups could be due to prevention or treatment of concomitant bacterial infections, thereby supporting current WHO recommendations (WHO, 2020[38]). However, categorization of antibiotics as a single category in this dataset limits our study. Randomized clinical trials testing antibiotics in severe COVID-19 adult patients to prevent secondary infections would be important to corroborate this hypothesis.

For children and adolescents, antibiotics were a risk factor for death (HR=4.22, 95 % CI: 2.01-8.86). However, we did not differentiate ambulatory from hospitalized pediatric patients and current recommendations include using antibiotics in hospitalized patients with multisystem inflammatory syndrome (Harwood et al., 2021[11]). The lack of benefit from antivirals included in our study in pediatric patients supports current guidelines discouraging their use (Ye et al., 2020[40]).

The main limitation of our study is that we were not able to assess cointerventions being studied for COVID-19 since only data for antivirals and antibiotics were available. Steroids increase survival in patients requiring oxygen administration and decrease survival in patients without supplementary oxygen (RECOVERY Collaborative Group et al., 2020[30]). However, modifications of mortality risks by steroids in our population should be large according to E-values, reducing the likelihood that these associations could be due to steroids.

Another potential limitation is that Mexico has had a low diagnostic testing rate for SARS-CoV-2 (0.08 daily tests per 1,000 people) (Roser et al., 2020[31]). However, health authorities require 100 % of patients with severe disease to be tested. Since we only studied mortality, an outcome expected to occur in patients who progress to severe disease, our study feasibly included most events. Nonetheless, excess mortality rates suggest there could be an undercounting of deaths in Mexico City (Roser et al., 2020[31]). These patients could have refrained from seeking medical attention or received medical care in non-accredited COVID-19 units where mortality, quality of care, and use of antivirals/antibiotics could be different.

Also, the number of ICU beds in Mexico City was relatively low in March 2020 (6.0 per 100,000 population) compared to most European countries (5 to 33.9 per 100,000) in the pre-pandemic period; this capacity was expanded to 29.5 ICU beds per 100,000 by September 2020 (Karagiannidis et al., 2020[15]; Roser et al., 2020[31]). Mortality rates, especially in patients younger than 60 years, are lower under high availability of ICU beds (Karagiannidis et al., 2020[15]). Thus, mortality rates could have varied throughout our study period.

Although we were not able to determine duration of follow-up in our study, the mechanisms and resources used by epidemiologic authorities in Mexico are robust enough to guarantee adequate matching of patients who had completed follow-up with death certificates. Thus, our finding that 92.7 % (95 % CI: 92.2-93.2) and 99.6 % (95 % CI: 99.5-99.7) of deaths occurred by day 28 and 56, respectively, could be important for the interpretation and design of COVID-19 clinical trials assessing short-term mortality.

Conclusions

Repurposed antivirals (oseltamivir, zanamivir, amantadine, rimantadine, acyclovir, and lopinavir/ritonavir) for COVID-19 did not provide additional benefit, whereas oseltamivir was consistently associated with increased mortality risk. Antibiotics were associated with increased risk of death in the general population but lower mortality risk only in hospitalized and critical patients, thereby supporting a potential benefit for their use in patients with severe-to-critical disease; this hypothesis should be tested in clinical trials.

Acknowledgements

We thank Dr. César Sebastián Salinas Nájera, Head of Epidemiology in “Hospital Materno Infantil Magdalena Contreras” of Mexico City's Health Secretariat for his valuable clarifications of epidemiologic surveillance mechanisms in Mexico City. We want to acknowledge the Government of Mexico City for their effort to increase transparency and aid research by providing open datasets on their Open Data Platform; this sort of actions is most needed in low and middle-income countries to understand and solve the problems of our people. JMG would like to thank “Dirección General de Calidad y Educación en Salud” for supporting his participation in “Programa Nacional de Servicio Social en Investigación en Salud”.

Conflict of interest

All authors declare no competing interests.

Funding

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Data availability statement

The data that support the findings of this study are openly available on the Open Data Platform of Mexico City’s Government at https://datos.cdmx.gob.mx/dataset/base-covid-sinave, and the Directorate General of Epidemiology of Mexico City’s Weekly Epidemiological Surveillance Reports of Influenza 2020 at https://www.gob.mx/salud/documentos/informes-semanales-para-la-vigilancia-epidemiologica-de-influenza-2020.

Ethics statement

This is a retrospective analysis of an open-source dataset of patients evaluated for suspected COVID-19 in Mexico City. The Ministry of Health of Mexico City approved the collection and release of publicly available datasets.

Author contributions

Conception and design: JMG and AKG; Analysis and interpretation: AKG, JMG, JMS, JOGM, and HIRG; Drafting the manuscript, making tables and figures: All authors; Revising the manuscript for important intellectual content: AKG, JOGM, and HIRG. All authors approved the final draft of the manuscript.

Supplementary Material Supplementary appendix Supplementary information Aranda Abreu GE Hernández Aguilar ME Herrera Covarrubias D Rojas Durán F Amantadine as a drug to mitigate the effects of COVID-19 Med Hypotheses 2020 140 109755 Araújo R Aranda-Martínez JD Aranda-Abreu GE Amantadine treatment for people with COVID-19 Arch Med Res 2020 51 739–40 Benlloch J-M Cortés J-C Martínez-Rodríguez D Julián R-S Villanueva R-J Effect of the early use of antivirals on the COVID-19 pandemic. A computational network modeling approach Chaos, Solitons & Fractals 2020 140 110168 Bhimraj A Morgan RL Shumaker AH Lavergne V Baden L Cheng VC Infectious Diseases Society of America Guidelines on the Treatment and Management of Patients with COVID-19. IDSA Guidelines 2020 8 December 2020 Available from: https://www.idsociety.org/practice-guideline/covid-19-guideline-treatment-and-management/ Díaz-Portillo SP Idrovo ÁJ Dreser A Bonilla FR Matías-Juan B Wirtz VJ Consultorios adyacentes a farmacias privadas en México: Infraestructura y características del personal médico y su remuneración Clinics adjacent to private pharmacies in Mexico: infrastructure and characteristics of the physicians and their remuneration, in Spanish Salud Publica Mex 2015 57 320 328 Directorate General of Epidemiology of Mexico Standardized Guideline for Epidemiologic and Laboratory Surveillance of viral respiratory diseases, August 2020 8 December 2020 Mexico City Mexican Secreatariat of Health 1–72 Available from: https://www.gob.mx/salud/documentos/lineamiento-estandarizado-para-la-vigilancia-epidemiologica-y-por-laboratorio-de-la-enfermedad-respiratoria-viral Directorate General of Epidemiology of Mexico Weekly Epidemiological Surveillance Reports of Influenza 2020 2020 8 December 2020 Available from: https://www.gob.mx/salud/documentos/informes-semanales-para-la-vigilancia-epidemiologica-de-influenza-2020 Falavigna M Colpani V Stein C Azevedo LCP Bagattini AM de Brito GV Guidelines for the pharmacological treatment of COVID-19. The task force/consensus guideline of the Brazilian Association of Intensive Care Medicine, the Brazilian Society of Infectious Diseases and the Brazilian Society of Pulmonology and Tisiology Rev Bras Ter Intensiva 2020 32 166–96 Government of Mexico City Covid-19 National Epidemiologic Surveilance System (SINAVE) Mexico City. Open Data Platform 2020 14 September 2020 Available from: https://datos.cdmx.gob.mx/dataset/base-covid-sinave Hama R The mechanisms of delayed onset type adverse reactions to oseltamivir Infect Dis (Auckl) 2016 48 651–60 Harwood R Allin B Jones CE Whittaker E Ramnarayan P Ramanan AV A national consensus management pathway for paediatric inflammatory multisystem syndrome temporally associated with COVID-19 (PIMS-TS): Results of a national Delphi process Lancet Child Adolesc Health 2021 5 133–41 Horby PW Mafham M Bell JL Linsell L Staplin N Emberson J Lopinavir–ritonavir in patients admitted to hospital with COVID-19 (RECOVERY): a randomised, controlled, open-label, platform trial Lancet 2020 396 1345–52 Islam N Sharp SJ Chowell G Shabnam S Kawachi I Lacey B Physical distancing interventions and incidence of coronavirus disease 2019: Natural experiment in 149 countries BMJ 2020 370 m2743 Jefferson T Jones MA Doshi P Del Mar CB Hama R Thompson MJ Neuraminidase inhibitors for preventing and treating influenza in adults and children Cochrane Database Syst Rev 2014 2014 4 CD008965 Karagiannidis C Mostert C Hentschker C Voshaar T Malzahn J Schillinger G Case characteristics, resource use, and outcomes of 10 021 patients with COVID-19 admitted to 920 German hospitals: An observational study Lancet Respir Med 2020 8 853–62 Li Z Yang L Underlying mechanisms and candidate drugs for COVID-19 based on the connectivity map database Front Genet 2020 11 558557 Lipsitch M Perlman S Waldor MK Testing COVID-19 therapies to prevent progression of mild disease Lancet Infect Dis 2020 20 1367 Liu J Zhang S Wu Z Shang Y Dong X Li G Clinical outcomes of COVID-19 in Wuhan, China: A large cohort study Ann Intensive Care 2020 10 99 Long B Brady WJ Koyfman A Gottlieb M Cardiovascular complications in COVID-19 Am J Emerg Med 2020 38 1504–7 López-Manning M García-Díaz R Doctors adjacent to private pharmacies: The new ambulatory care provider for Mexican health care seekers Value Heal Reg Issues 2017 14 81–8 Mancilla-Galindo J Kammar-García A Martínez-Esteban A Meza-Comparán A-K Mancilla-Ramírez J Galindo-Sevilla N COVID-19 patients with increasing age experience differential time to initial medical care and severity of symptoms Cambridge Open Engag 06 Jan 2021 Version 1 10.33774/coe-2021-sjbcf Mancilla-Galindo J Vera-Zertuche JM Navarro-Cruz AR Segura-Badilla O Reyes-Velázquez G Tepepa-López FJ Development and validation of the patient history COVID-19 (PH-Covid19) scoring system: A multivariable prediction model of death in Mexican patients with COVID-19 Epidemiol Infect 2020 148 e286 Mexican Secretariat of Health Diagnostic algorithm for acute respiratory disease in the context of the COVID19-influenza syndemic. Guidance documents for healthcare professionals 2020 8 December 2020 20 Available from: https://coronavirus.gob.mx/wp-content/uploads/2020/07/Algoritmo_Dx_Personal_Salud_08072020.pdf Mexican Secretariat of Health Recommendations for treatment of SARS-CoV-2 infection, the etiologic agent of COVID-19. Vol. July, Guidance documents for healthcare professionals 2020 8 December 2020 7 Available from: https://coronavirus.gob.mx/wp-content/uploads/2020/07/Recomendaciones_para_tratamiento_SARS-CoV2.pdf Moore ML Chi MH Zhou W Goleniewska K O’Neal JF Higginbotham JN Cutting edge: Oseltamivir decreases T cell GM1 expression and inhibits clearance of respiratory syncytial virus: Potential role of endogenous sialidase in antiviral immunity J Immunol 2007 178 2651–4 Muralidharan N Sakthivel R Velmurugan D Gromiha MM Computational studies of drug repurposing and synergism of lopinavir, oseltamivir and ritonavir binding with SARS-CoV-2 protease against COVID-19 J Biomol Struct Dyn 2020 epub ahead of print 1 6 10.1080/07391102.2020.1752802 NIH, National Institutes of Health COVID-19 Treatment Guidelines Panel. Coronavirus Disease 2019 (COVID-19) Treatment Guidelines 2020 7 October 2020 Available from: https://www.covid19treatmentguidelines.nih.gov/ NIHR, National Institute for Health Research NIHR Innovation Observatory. COVID-19 therapeutics 2020 8 December 2020 Available from: http://www.io.nihr.ac.uk/report/covid-19-therapeutics/ Pérez-Cuevas R Doubova SV Wirtz VJ Servan-Mori E Dreser A Hernández-Ávila M Effects of the expansion of doctors’ offices adjacent to private pharmacies in Mexico: Secondary data analysis of a national survey BMJ Open 2014 4 e004669 RECOVERY Collaborative Group, Horby P, Lim WS, Emberson JR, Mafham M, Bell JL, Dexamethasone in hospitalized patients with Covid-19 - preliminary report N Engl J Med 2020 NEJMoa2021436 Roser M Ritchie H Ortiz-Ospina E Hasell J Coronavirus pandemic (COVID-19) 2020 20 September 2020 Available from: https://ourworldindata.org/coronavirus Sachs JD Abdool Karim S Aknin L Allen J Brosbøl K Cuevas Barron G Lancet COVID-19 Commission Statement on the occasion of the 75th session of the UN General Assembly Lancet 2020 396 1102–24 Saul S Einav S Old drugs for a new virus: Repurposed approaches for combating COVID-19 ACS Infect Dis 2020 6 2304–18 Senanayake SL Drug repurposing strategies for COVID-19 Future Drug Discov 2020 2 fdd 2020 10.4155/fdd-2020-0010 Shi Q Zhou Q Wang X Liao J Yu Y Wang Z Potential effectiveness and safety of antiviral agents in children with coronavirus disease 2019: A rapid review and meta-analysis Ann Transl Med 2020 8 10 624 Tan Q Duan L Ma Y Wu F Huang Q Mao K Is oseltamivir suitable for fighting against COVID-19: In silico assessment, in vitro and retrospective study Bioorg Chem 2020 104 104257 Vera-Zertuche JM Mancilla-Galindo J Tlalpa-Prisco M Aguilar-Alonso P Aguirre-García MM Segura-Badilla O Obesity is a strong independent risk factor for short-term mortality and adverse outcomes in Mexican patients with COVID-19: A national observational study Cambridge Open Engag 06 Jan 2021 Version 1 10.33774/coe-2021-qf4r6 WHO, World Health Organization Clinical management of COVID-19 - interim guidance, 27 May 2020 2020 8 December 2020 Geneva WHO Available from: https://apps.who.int/iris/handle/10665/332196 WHO, World Health Organization International Clinical Trials Registry Platform (ICTRP) 2020 8 December 2020 Geneva WHO Available from: https://www.who.int/ictrp/en/ Ye Z Rochwerg B Wang Y Adhikari NK Murthy S Lamontagne F Treatment of patients with nonsevere and severe coronavirus disease 2019: An evidence-based guideline Can Med Assoc J 2020 192 E536–45 Baseline characteristics of patients with laboratory confirmed COVID-19 who were treated with or without antivirals/antibiotics, in 688 accredited COVID-19 medical units in Mexico City Adjusted mortality risk in laboratory-confirmed COVID-19 patients receiving antivirals, antibiotics, both, or none in 688 accredited COVID-19 medical units in Mexico City Adjusted mortality risk in laboratory-confirmed COVID-19 patients receiving antivirals, oseltamivir, or antibiotics after propensity score matching Flow diagram of patients assessed for eligibility Survival of patients (general population, ambulatory, and hospitalized) treated with antivirals and/or antibiotics.

Survival curves are shown according to treatment modality in the general population (a), ambulatory (b), and hospitalized (c) patients. Survival in patients receiving specific antivirals, antibiotics, , or none in the general population (d), ambulatory (e), and hospitalization (f) settings.

Survival of patients (non-critical and critical) treated with antivirals and/or antibiotics.

Survival curves of non-critical patients according to modality of treatment (a) and those receiving specific antivirals, antibiotics, or none (b). Survival curves of critical patients according to modality of treatment (c) and those receiving specific antivirals, antibiotics, or none (d).

Weekly prescription of oseltamivir and newly diagnosed patients with influenza (A) and COVID-19 (B) during the COVID-19 pandemic in Mexico City.

*Data of prescription of oseltamivir for patients tested only for influenza was not collected.