Molnupiravir and risk of post-acute sequelae of covid-19: cohort study
BMJ 2023; 381 doi: https://doi.org/10.1136/bmj-2022-074572 (Published 25 April 2023) Cite this as: BMJ 2023;381:e074572- Yan Xie, senior clinical epidemiologist12,
- Taeyoung Choi, biostatistician12,
- Ziyad Al-Aly, director12345
- 1Clinical Epidemiology Center, Research and Development Service, VA Saint Louis Health Care System, Saint Louis, Missouri, USA
- 2Veterans Research and Education Foundation of Saint Louis, Saint Louis, Missouri, USA
- 3Department of Medicine, Washington University School of Medicine, Saint Louis, Missouri, USA
- 4Nephrology Section, Medicine Service, VA Saint Louis Health Care System, Saint Louis, Missouri, USA
- 5Institute for Public Health, Washington University in Saint Louis, Saint Louis, Missouri, USA
- Correspondence to: Z Al-Aly VA Saint Louis Health Care System, Saint Louis, MO 63106, USA zalaly{at}gmail.com (or @zalaly on Twitter)
- Accepted 14 March 2023
Abstract
Objective To examine whether treatment with the antiviral agent molnupiravir during the first five days of SARS-CoV-2 infection is associated with reduced risk of post-acute adverse health outcomes.
Design Cohort study.
Setting US Department of Veterans Affairs.
Participants 229 286 participants who tested positive for SARS-CoV-2 between 5 January 2022 and 15 January 2023, had at least one risk factor for progression to severe covid-19, and survived the first 30 days after testing positive were enrolled. 11 472 participants received a prescription for molnupiravir within five days of the positive test result and 217 814 received no covid-19 antiviral or antibody treatment (no treatment group).
Main outcome measures Risks of post-acute sequelae of SARS-CoV-2 (PASC, defined based on a prespecified set of 13 post-acute sequelae), post-acute death, post-acute hospital admission, and each individual post-acute sequela between the molnupiravir group and no treatment group were examined after application of inverse probability weighting to balance the treatment and no treatment groups. Post-acute outcomes were ascertained from 30 days after the first SARS-CoV-2 positive test result until end of follow-up. Risks on the relative scale (relative risk or hazard ratio) and absolute scale (absolute risk reduction at 180 days) were estimated.
Results Compared with no treatment, molnupiravir use within five days of a positive SARS-CoV-2 test result was associated with reduced risk of PASC (relative risk 0.86 (95% confidence interval 0.83 to 0.89); absolute risk reduction at 180 days 2.97% (95% confidence interval 2.31% to 3.60%)), post-acute death (hazard ratio 0.62 (0.52 to 0.74); 0.87% (0.62% to 1.13%)), and post-acute hospital admission (0.86 (0.80 to 0.93); 1.32% (0.72% to 1.92%)). Molnupiravir was associated with reduced risk of eight of the 13 post-acute sequelae: dysrhythmia, pulmonary embolism, deep vein thrombosis, fatigue and malaise, liver disease, acute kidney injury, muscle pain, and neurocognitive impairment. Molnupiravir was also associated with reduced risk of PASC in people who had not received a covid-19 vaccine, had received at one or two vaccine doses, and had received a booster dose, and in people with primary SARS-CoV-2 infection and reinfection.
Conclusions In people with SARS-CoV-2 infection and at least one risk factor for progression to severe covid-19, compared with no treatment, molnupiravir use within five days of infection was associated with reduced risk of PASC in people who had not received a covid-19 vaccine, had received one or two vaccine doses, and had received a booster dose, and in those with primary SARS-CoV-2 infection and reinfection. Among people at high risk of progression to severe covid-19, molnupiravir use within five days of SARS-CoV-2 infection may be a viable approach to reduce the risk of PASC.
Introduction
Long covid represents the constellation of post-acute and chronic sequelae of SARS-CoV-2 infection.123 At least 65 million individuals globally are estimated to have long covid.4 Despite progress in characterizing long covid over the past three years, however, drug interventions for its prevention or treatment are not yet available.4
Nearly all the proposed mechanistic hypotheses to explain the myriad manifestations of long covid involve SARS-CoV-2 as the initiating agent.5 We previously reported that use of the antiviral nirmatrelvir (in combination with ritonavir) within five days of SARS-CoV-2 infection is associated with reduced risk of long covid.6 Collectively, these observations suggest that treatment of SARS-CoV-2 with antivirals during the early, acute phase may reduce the risk of post-acute adverse health outcomes.
Molnupiravir is an antiviral agent that impairs SARS-CoV-2 replication.7 Molnupiravir use within five days of testing positive for SARS-CoV-2 has been shown to reduce the risk of hospital admission or death in people at risk of progression to severe acute covid-19.89 In December 2021, the US Food and Drug Administration issued an emergency use authorization for molnupiravir in the treatment of people with acute SARS-CoV-2 infection at high risk of progression to severe covid-19 for whom alternative covid-19 treatment options were not accessible or were clinically inappropriate. Subsequently, millions of people in the United States have been treated with molnupiravir. The question of whether molnupiravir use during the acute phase of SARS-CoV-2 infection reduces the risk of long covid has been highlighted as a major knowledge gap that should be urgently addressed.10 Addressing this question will help inform treatment approaches for SARS-CoV-2 infections and guide the development and testing of additional drugs to prevent and treat long covid.
We used the electronic health record records of the US Department of Veterans Affairs to identify patients with a positive test result for SARS-CoV-2 between 5 January 2022 and 15 January 2023 and at least one risk factor for progression to severe covid-19. Weighting was then used to balance the personal and health characteristics of those who received a prescription for molnupiravir within five days of testing positive for SARS-CoV-2 (n=11 472) and those who did not receive outpatient covid-19 antiviral or antibody treatment within the first 30 days of infection (n=217 814). We evaluated whether treatment with molnupiravir was associated with reduced risk of post-acute sequelae of SARS-CoV-2 (PASC, based on 13 prespecified post-acute sequelae), post-acute death, post-acute hospital admission, and each individual post-acute sequela, where post-acute outcomes were ascertained 30 days after a first positive SARS-CoV-2 test result until end of follow-up.
Methods
Setting
Study data were collected from 1293 Veterans Affairs healthcare facilities (171 Veterans Affairs medical centers and 1112 outpatient sites) across the US. Veterans Affairs operates the largest integrated healthcare system in the US and provides healthcare services to veterans of the US armed forces. Healthcare provided includes preventive and health maintenance, outpatient care, inpatient hospital care, prescriptions, mental health care, home healthcare, primary care, specialty care, geriatric and extended care, medical equipment, and prosthetics.
Data sources
We extracted data from the US Department of Veterans Affairs healthcare databases. The databases are updated daily and include information collected during patients’ routine healthcare encounters. Data domains used were outpatient and inpatient diagnoses, outpatient and inpatient pharmacy records, and laboratory test results. Vaccination status was collected from the VA COVID-19 Shared Data Resource. Information about care outside of Veterans Affairs was collected from non-Veterans Affairs care program integrity tools data domain.11 The area deprivation index, a composite measure of income, education, employment, or housing, was used as a summary measure of contextual disadvantage at participants’ residential locations.12
Cohort
Supplementary figures S1 and S2 show the flowchart and timeline plots of cohort construction, respectively. We enrolled 339 719 participants with a positive SARS-CoV-2 test result between 5 January 2022 and 15 January 2023, with the date of a first positive result considered as T0. In accordance with the emergency use authorization issued by the FDA, which recommends initiation of molnupiravir within five days of SARS-CoV-2 infection, 13 007 of those who tested positive received a prescription for molnupiravir within five days of T0 and were selected for the molnupiravir group. We specified the inclusion and exclusion criteria in accordance with the MOVe-OUT (Efficacy and Safety of Molnupiravir in Non-Hospitalized Adult Participants With COVID-19) trial and the emergency use authorization issued by the FDA.813 Overall, we selected 12 684 participants with at least one of the following risk factors for progression to severe acute covid-19: age >60 years, body mass index (BMI) >30, cancer, cardiovascular disease, chronic kidney disease, chronic lung disease, diabetes, and immune dysfunction. We excluded participants with contraindications to treatment: end stage kidney disease or an estimated glomerular filtration rate (eGFR) <30 mL/min/1.73m2 (remaining n=11 699). To examine the effect of molnupiravir on post-acute events, we further excluded those who used other outpatient covid-19 antiviral or antibody treatments within 30 days after T0, leaving 11 570 participants, and those who died within 30 days after the testing positive, resulting in a cohort comprising 11 472 participants.
We constructed the no treatment group using the same inclusion and exclusion criteria as for the molnupiravir group. The no treatment control group were selected from the 326 712 participants with a positive SARS-CoV-2 test result between 5 January 2022 and 15 January 2023 and who did not receive a prescription for molnupiravir within five days of testing positive. The specified inclusion and exclusion criteria enabled the selection of participants eligible for molnupiravir in accordance with the MOVE-OUT trial and the emergency use authorization issued by the FDA.813 Overall, we selected 292 408 of the participants with at least one risk factor for progression to severe acute covid-19 (see supplementary figure S1). We also excluded participants with a contraindication to molnupiravir—those with end stage kidney disease or an eGFR <30 mL/min/1.73m2 (remaining n=281 876). We then further excluded those who used outpatient covid-19 antiviral or antibody treatments within 30 days after T0, leaving 230 714 participants, and those who died within 30 days after testing positive, resulting in a cohort comprising 217 814 participants. The cohort was followed until 15 February 2023.
Outcomes
The outcomes measured included post-acute death and post-acute hospital admission and a composite of post-acute death or hospital admission. We also examined the risk of developing 13 prespecified post-acute sequelae of SARS-CoV-2: incident ischemic heart disease, dysrhythmia, deep vein thrombosis, pulmonary embolism, fatigue and malaise, liver disease, acute kidney injury, muscle pain, diabetes, neurocognitive impairment, dysautonomia, shortness of breath, and cough. Outcomes were ascertained as follows: deaths from vital status data, hospital admissions from inpatient encounter data, and individual sequelae from inpatient and outpatient ICD-10 (international classification of diseases, 10th revision) diagnosis codes, pharmacy records, and laboratory test results.12141516171819202122232425262728293031 Incident outcomes were assessed in those with no history of the related outcome within one year before testing positive.
We built the PASC score by assigning weights to each sequela (using the methods of the Global Burden of Disease Long COVID Collaborators32); this approach accounts for the number of sequelae experienced by an individual participant and the influence of each sequela (compared with other diseases) on health.333435 The PASC score was then constructed for each participant as the sum of weights of all the incident sequelae in that participant (see https://github.com/yxie618/Molnupiravir_PASC for information on the weights and PASC score computation). To provide generalizable estimates for the burden of PASC on the absolute scale, we standardized the rate distribution of the PASC score to the rate distribution of the outcome, which was defined as having at least one sequela. All post-acute outcomes were ascertained 30 days after the first positive SARS-CoV-2 test result between 5 January 2022 and 15 January 2023 until 15 February 2023.
Covariates
Based on literature review and previous knowledge, we selected baseline covariates that might influence the treatment assignment and the occurrence or assessment of outcomes.127363738 Supplementary figure S3 shows the directed acyclic graph, which guided the selection of covariates. The ascertainment period for covariates was within three years before the baseline positive SARS-CoV-2 test result unless specified otherwise, and we used covariate values before and closest to the test date. Personal and health factors for age, race (white, black, and other), sex (men or women), area of deprivation index, state of residency, BMI, smoking status (never, former, and current), history of SARS-CoV-2 infection, and use of long term care were selected. We also controlled for clinical factors—use of steroids, eGFR, systolic and diastolic blood pressure, cancer, chronic lung disease, dementia, diabetes, hyperlipidemia, and immune dysfunction. Healthcare utilization variables were history of influenza vaccination, number of outpatient and inpatient encounters, number of laboratory tests, number of outpatient prescriptions, and number of Medicare related outpatient and inpatient encounters. We additionally specified pandemic related characteristics for week of the positive SARS-CoV-2 test result, hospital bed capacity, and hospital bed occupancy at the participants’ healthcare facility within the week of the positive test result. Hospital bed capacity was computed as the sum of number of beds occupied and number of beds available based on information reported from each healthcare facility. Hospital bed occupancy was computed as number of beds occupied divided by hospital bed capacity.36 Because the FDA’s emergency use authorization specified that molnupiravir may be considered when alternative options for covid-19 treatment (ie, nirmatrelvir-ritonavir) are not accessible or clinically appropriate, medicines that would have drug-drug interaction with and might preclude use of nirmatrelvir-ritonavir were classified at four levels (require an alternative covid-19 treatment to nirmatrelvir-ritonavir, require temporary suspension of concomitant drug, require adjustment of concomitant drug dosing, and require monitoring for adverse effects) and were adjusted for in the models. Within the study cohort, 8.0% of the eGFR, 9.1% of the BMI, and 3.9% of the blood pressure measurements were missing. We imputed missing data using the predictive mean matching method based on multivariate imputation by chained equations.39 To account for potential non-linear relationships, we transformed continuous variables into restricted cubic spline functions.
Statistical analysis
Baseline characteristics of the overall cohort and the molnupiravir and no treatment groups were reported, when appropriate, as mean and standard deviation or frequency and percentage. Balance of baseline characteristics between the two study groups was evaluated using absolute standardized differences, with a value <0.1 considered evidence of good balance.
To evaluate the risk of outcomes in the two study groups, we used inverse probability weighting to balance baseline characteristics between the groups. Multivariate logistic regression was applied to estimate the probability of receiving molnupiravir (the propensity score) given the covariates. The inverse probability weights to estimate treatment effect within treated participants was then constructed as 1 for those in the molnupiravir group and as propensity score/(1−propensity score) for those in the no treatment group. Although we prespecified our analytic approach to truncate weights that were larger than 10 to reduce the influence of extreme weights, no weights were that large and hence no truncation was required. To estimate the effect of molnupiravir on PASC score, we applied the inverse probability weights to a zero inflated Poisson regression, where the follow-up time was set to be the offset of the model. From the Poisson regression we estimated the relative risk (incident rate ratio) and absolute risk reduction at 180 days. The inverse probability weights were also applied to Cox survival models to estimate the effect of molnupiravir on risk of post-acute death, post-acute hospital admission, post-acute death or hospital admission, and individual post-acute sequela. Hazard ratio and absolute risk reduction were estimated from the models, where risk reduction was estimated as the difference in survival probability between the molnupiravir group and no treatment group at 180 days.
The effect of molnupiravir on risk of PASC was also examined in those who had received no covid-19 vaccine, received one or two vaccine doses, and received at least one booster dose. We also examined the effect in those with primary SARS-CoV2 infection and reinfection. We then further examined the effect within prespecified subgroups of age (≤65 and >65 years), race (white and black), sex (women and men), smoking status (never, former, and current), and presence or absence of risk factors for progression to severe acute covid-19: cancer, cardiovascular disease, kidney disease, chronic lung disease, diabetes, and immune dysfunction. Interaction analyses were conducted to examine whether subgroup status (eg, being in different age subgroups) modified the association between molnupiravir and PASC. Inverse probability weighting and outcome model were constructed independently within each subgroup examined. The effect of molnupiravir on risk of PASC across ages was also examined, where age was treated as a cubic spline function.
We conducted multiple sensitivity analyses for the outcome of PASC and the composite outcome of post-acute death or hospital admission (when applicable) to examine the robustness of findings by: balancing baseline characteristics through overlap weighting, instead of inverse probability weighting used in the primary approach; additionally adjusting for covariates in the weighted outcome model using a doubly robust approach, compared with the primary approach, which only applied inverse probability weighting to the outcome model; utilizing additional 100 covariates selected from a high dimensional variable selection algorithm to construct the propensity score, compared with the primary approach which used only prespecified covariates; using a negative binomial model to estimate the effect, compared with the primary approach, which used a Poisson model; enrolling those who died during the acute phase of infection (within 30 days of a SARS-CoV-2 positive test result) and applying inverse probability of censoring weight to account for the death, compared with the primary approach, which removed those who died during the acute phase from the analyses; examining effects within those who received treatment on the same day as the positive SARS-CoV-2 test result and, separately, within those who received treatment 2-5 days after testing positive, compared with the primary approach, which estimated the effect in people who received treatment within five days of testing positive; examining the effect in those who were not admitted to hospital during the acute phase of infection, compared with the primary approach, which included participants who were and were not admitted to hospital during the acute phase of infection; additionally adjusting for hospital admission, intensive care unit admission, and ventilator use during the acute phase of infection, compared with the primary approach, which used covariate information at T0; additionally adjusting for healthcare utilization, including outpatient encounters, hospital admissions, laboratory tests, and medicines during follow-up in the outcome model, compared with the primary approach, which used covariate information at T0; defining outcomes based on events that occurred 90 days after infection, compared with the primary approach, which defined outcomes based on events that occurred 30 days after infection; defining outcomes based on events that occurred any time after the infection, compared with the primary approach, which defined outcomes based on events that occurred 30 days after infection; and assigning the same weight to each individual sequela, compared with the primary approach, in which each individual sequela was weighted based on its health burden.
We considered a relative risk or hazard ratio value (and 95% confidence interval) that did not cross 1 and an absolute risk reduction value (and 95% confidence interval) that did not cross 0 to be statistically significant. Data management and analyses were performed with SAS Enterprise Guide, version 8.2 (SAS Institute, Cary, NC). Data visualizations were performed in R 4.2.2 (R Foundation for Statistical Computing, Vienna, Austria).
Patient and public involvement
No patients or members of the public were directly involved in developing the research question, study design, and implementation. We did, however share drafts of the manuscript with members of the public to obtain feedback, including informing the selection of outcomes, which were incorporated in the revisions.
Results
The overall cohort comprised 229 286 participants with a positive SARS-CoV-2 test result between 5 January 2022 and 15 January 2023 and at least one risk factor for progression to severe covid-19. Of these participants, 11 472 received a prescription for molnupiravir within five days after testing positive and 217 814 received no covid-19 outpatient treatment within the first 30 days after testing positive. Supplementary table S1 presents the baseline personal and health characteristics of the molnupiravir and no treatment groups before inverse probability weighting. Table 1 presents the characteristics of the study groups after inverse probability weighting. Values for the absolute standardized mean differences of personal and health characteristics between the molnupiravir and no treatment groups after weighting were <0.1, indicating good balance between the groups (table 1).
Risk of PASC, post-acute death, post-acute hospital admission, and post-acute sequalae
Molnupiravir use was associated with a reduced risk of PASC (relative risk 0.86 (95% confidence interval 0.83 to 0.89) and absolute risk reduction at 180 days 2.97% (95% confidence interval 2.31% to 3.60%)) (fig 1, fig 2, table 2, and supplementary table S2) compared with no treatment. The corresponding event rates of PASC at 180 days were 18.58% (17.94% to 19.24%) in the molnupiravir group and 21.55% (21.04% to 22.07%) in the no treatment group.
Compared with no treatment, molnupiravir use was associated with reduced risk of post-acute death (hazard ratio 0.62 (95% confidence interval 0.52 to 0.74) and absolute risk reduction 0.87% (0.62% to 1.13%)), post-acute hospital admission (0.86 (0.80 to 0.93) and 1.32% (0.72% to 1.92%)), and composite outcome of post-acute death or hospital admission (0.82 (0.77 to 0.88) and 2.00% (1.37% to 2.62%) (fig 1, fig 2, and table 2).
Molnupiravir use was associated with reduced risk of eight of the 13 individual post-acute sequela examined compared with no treatment—dysrhythmia, pulmonary embolism, deep vein thrombosis, fatigue and malaise, liver disease, acute kidney injury, muscle pain, and neurocognitive impairment. We did not observe statistically significant associations between molnupiravir and ischemic heart disease, diabetes, dysautonomia, shortness of breath, and cough (fig 3 and supplementary table S3).
Risk of PASC in subgroups
Compared with no treatment, molnupiravir use was associated with reduced risk of PASC in those who had not received a covid-19 vaccine, received one or two vaccine doses, and received a booster dose. Molnupiravir use was also associated with reduced risk of PASC in those with primary SARS-CoV-2 infection and reinfection (fig 4 and supplementary table S4).
In subgroup analyses based on age (≤65 and >65 years), race (white and black), sex (men and women), smoking status (never, former, and current), and disease status (cancer, cardiovascular disease, chronic kidney disease, chronic lung disease, diabetes, and immune dysfunction), molnupiravir use was associated with reduced risk of PASC compared with no treatment (fig 5, supplementary figure S4, supplementary table S4). No statistically significant interaction was observed in any of the subgroups.
Sensitivity analyses
Multiple sensitivity analyses were conducted to assess the robustness of our findings. First, we balanced baseline characteristics through overlap weighting, instead of inverse probability weighting used in the primary approach. Second, we additionally adjusted for covariates in the weighted outcome model using a doubly robust methodology, instead of inverse probability weighting to the outcome model used in the primary approach. Third, we utilized additional 100 covariates selected using a high dimensional variable selection algorithm to construct the propensity score, instead of using only prespecified covariates as in the primary approach. Fourth, we used a negative binomial model to estimate the effect, instead of a Poisson model as used in the primary approach. Fifth, we enrolled those who died during the acute phase of infection (within 30 days of a positive SARS-CoV-2 test result) and applied inverse probability of censoring weight to account for the death, instead of excluding those who died during the acute phase from analyses as was used in the primary approach. Sixth, we examined effect within those who received treatment on the same day as a positive SARS-CoV-2 test result, and separately within those who received treatment 2-5 days after the date of testing positive, compared with the primary approach, which estimated the effect in people who received treatment within five days of testing positive. Seventh, we examined the effect in those who were not admitted to hospital during the acute phase of infection, compared with the primary approach, which included participants who were and were not admitted to hospital during the acute phase. Eighth, we additionally adjusted for hospital admission, intensive care unit admission, and ventilator use during the acute phase of infection, instead of using covariate information at T0 as was applied in the primary approach. Ninth, we additionally adjusted for healthcare utilization during follow-up (outpatient encounters, hospital admissions, laboratory tests, and drugs) in the outcome model, compared with the primary approach, in which we used covariate information at T0. Tenth, we defined outcomes based on events that occurred 90 days after infection, compared with the primary approach, which defined outcomes based on events that occurred 30 days after infection Eleventh, we defined outcomes based on events that occurred any time after the infection, instead of defining outcomes based on events that occurred 30 days after infection as in the primary approach. Finally, we assigned the same weight to each individual sequela, compared with the primary approach, in which we weighted each individual sequela based on its health burden. We obtained consistent results in all sensitivity analyses, suggesting that, compared with no treatment, molnupiravir use was associated with reduced risk of PASC and post-acute death or hospital admission (see supplementary table S5).
Discussion
Among 229 286 people included in this study with SARS-CoV-2 infection and at least one risk factor for progression to severe covid-19, use of the antiviral molnupiravir within five days of testing positive for SARS-CoV-2 compared with no treatment was associated with reduced risk of PASC, post-acute death, and post-acute hospital admission. Molnupiravir use was associated with reduced risk of PASC in those who had not received a covid-19 vaccine, received one or two vaccine doses, and received a booster dose, and in people with primary SARS-CoV-2 infection and reinfection. Molnupiravir use was also associated with reduced risk of PASC in various subgroups based on age, race, sex, smoking status, cancer, cardiovascular disease, chronic kidney disease, chronic lung disease, diabetes, and immune dysfunction. Compared with no treatment, molnupiravir use was also associated with reduced risk of eight out of the 13 prespecified sequelae—dysrhythmia, pulmonary embolism, deep vein thrombosis, fatigue and malaise, liver disease, acute kidney injury, muscle pain, and neurocognitive impairment. Overall, the findings suggest that among people at high risk of progression to severe covid-19, molnupiravir use within five days of a positive SARS-CoV-2 test result was associated with reduced risk of post-acute adverse health outcomes.
Findings in relation to other studies
Our study findings provide evidence that use of molnupiravir within the first five days of a positive SARS-CoV-2 test result is associated with reduced risk of PASC. Molnupiravir use was associated with 14% less risk of PASC, 38% lower risk of post-acute death, and 14% lower risk of post-acute hospital admission; the magnitude of risk reduction on the absolute scale was also substantial—a 2.97%, 0.87%, and 1.32% lower absolute risk of PASC, post-acute death, and post-acute hospital admission for every 100 treated people between 30 and 180 days of infection. These findings align with previous observations that use of combined antivirals nirmatrelvir-ritonavir is associated with reduced risk of PASC, post-acute death, and post-acute hospital admission.6 Together, these findings suggest that use of antivirals in the acute phase of infection may not only help reduce the risk of progression to severe covid-19 but may be a viable strategy to reduce the risk of post-acute adverse health outcomes, including PASC. The clinical decision to use antivirals during acute SARS-CoV-2 infection should consider the effectiveness of antivirals in reducing the burden of death and disease in both the acute and the post-acute phases of covid-19.
We specified our enrollment criteria to follow the FDA’s emergency use authorization for molnupiravir—which is generally used in people with at least one risk factor for progression to severe covid-19 and for whom the alternative treatment nirmatrelvir is not accessible or clinically appropriate. Evidence points to the possibility that although molnupiravir might be effective in the acute phase of SARS-CoV-2 in adults at high risk of progression to severe covid-19, these beneficial effects may not translate to effectiveness in the acute phase among populations with relatively lower risks.8940 Whether the reduced risk of PASC (and the other post-acute health outcomes) observed in this report (in people with at least one risk factor for progression to severe covid-19) extends to other populations (eg, those without risk factors for progression to severe disease) is not known. Furthermore, whether higher dose or more prolonged duration of treatment results in higher risk reduction is also not clear. Also, for the millions of people already experiencing long covid, the question of whether initiation of antiviral treatment in the post-acute phase of covid-19 is effective in treating long covid demands urgent attention.
Our results show risk reduction in eight out of 13 prespecified post-acute sequelae. It is possible that these sequelae may be driven by various mechanistic pathways with different sensitivity to antiviral treatment during the acute disease phase; it is also possible that the putative effect of molnupiravir on these sequelae is variable and that a higher level of power and precision is required for the observation of statistically significant results.
The mechanisms underpinning the association between molnupiravir and PASC are not entirely clear. Because PASC is a consequence of SARS-CoV-2 infection, and because the incidence of PASC—although evident in those with mild infection—is generally associated with severity of acute infection, it is reasonable to hypothesize that impairing viral replication and reducing severity of acute infection by early use of antivirals may lead to reduced risk of PASC.69 This hypothesis is consistent with the putative mechanisms that are being proposed to explain the occurrence of PASC, including viral persistence, immune dysfunction, occult organ injury during the acute phase that may manifest in the post-acute phase, and microbiome dysbiosis. All these putative mechanisms of long covid originate with SARS-CoV-2 as the initiating agent, suggesting the hypothesis that reduced viral replication (through use of antivirals) may ameliorate these manifestations.
Strengths and limitations of this study
This study has several strengths. We conducted this study using real world data from the Veterans Affairs, which operates the largest nationally integrated healthcare system in the US. We comprehensively leveraged the breadth and depth of Veterans Affairs data sources to incorporate information in multiple domains, such as personal characteristics, diagnoses, laboratory tests, drugs, vital signs, and healthcare utilization as well as contextual factors from the previous three years before the start of follow-up. For this cohort study, we enrolled a larger number of participants than are typically included in randomized trials and evaluated a prespecified set of 13 post-acute sequelae informed by our extensive research on long covid.1241319232728363841 We conceptualized the PASC score by assigning weights to each individual sequela not only to reflect its occurrence but also to reflect the number of sequelae experienced by an individual, and the influence of each sequela (compared with other diseases) on health. We assessed robustness of our approach in multiple sensitivity analyses, and results were consistent.
This study has several limitations. Our cohort comprised veterans, who are predominately white men, which may limit generalizability of study findings. We conducted an observational cohort study using data from the Veterans Affairs electronic record system; although we carefully designed our enrollment criteria to reduce differences between the study groups and subsequently used a large set of covariates from multiple data domains to balance these groups, biases, such as misclassification and residual confounding, may still exist. We relied on Veterans Affairs pharmacy records to define use of molnupiravir and did not capture drug use outside the Veterans Affairs system. If a large proportion of participants in the no treatment group received molnupiravir or other covid-19 treatment outside of the Veterans Affairs healthcare system, and if association in those participants was directionally similar to that reported here, this could have biased the results toward the null. Our analyses were restricted to a set of 13 sequelae and did not comprehensibly include all possible post-acute sequelae of SARS-CoV-2 infection. Finally, the pandemic remains highly dynamic, SARS-CoV-2 continues to mutate and new variants and subvariants are constantly emerging, and immunity levels in populations are also changing as a function of infection, reinfection, and vaccine uptake. Several factors may lead to the emergence of antiviral resistance; PASC itself may change both in terms of clinical features and in magnitude of risk. As a result, the real world effectiveness of molnupiravir on post-acute outcomes may also change dynamically over time.
Conclusions
Among people with SARS-CoV-2 infection and at least one risk factor for progression to severe covid-19, molnupiravir use during the first five days of infection was associated with reduced risk of PASC, regardless of an individual’s vaccination status and history of infection. These findings suggest that treatment with molnupiravir in the acute phase of SARS-CoV-2 infection may be a viable strategy to reduce the risk of post-acute adverse health outcomes.
What is already known on this topic
Among people with SARS-CoV-2 infection, molnupiravir use within five days of symptom onset has been shown to reduce the risk of hospital admission or death in people at risk of progression to severe covid-19
What this study adds
In people with SARS-CoV-2 infection and at least one risk factor for progression to severe covid-19, compared with no treatment, molnupiravir use during the first five days after a positive SARS-CoV-2 test result was associated with reduced risk of post-acute sequelae of SARS-CoV-2
The reduced risk in the molnupiravir group was evident in those who had not received a covid-19 vaccine, received at least one vaccine dose, and had received a booster dose, and in those with primary SARS-CoV-2 infection and reinfection
Ethics statements
Ethical approval
This study was approved, and waiver of consent granted, by the institutional review board of the VA Saint Louis Health Care System (protocol No 1606333).
Data availability statement
The data that support the findings of this study are available from the US Department of Veterans Affairs (VA). VA data are made freely available to researchers behind the VA firewall with an approved VA study protocol (see https://www.virec.research.va.gov or contact the VA Information Resource Center at VIReC@va.gov for more information).
Acknowledgments
This study used data from the Veterans Affairs COVID-19 Shared Data Resource.
Footnotes
Contributors: ZAA and YX contributed to the development of the study concept and design. ZAA, YX, and TC contributed to data analysis and interpretation of results. ZAA and YX drafted the manuscript. ZAA and YX contributed to critical revision of the manuscript. ZAA provided administrative, technical, and material support, as well as supervision and mentorship. Each author contributed important intellectual content during manuscript drafting or revision and accepts accountability for the overall work by ensuring that questions pertaining to the accuracy or integrity of any portion of the work are appropriately investigated and resolved. All authors approved the final version of the report. ZAA is the guarantor. The contents do not represent the views of the US Department of Veterans Affairs or the US government. The corresponding author attests that all the listed authors meet the authorship criteria and that no others meeting the criteria have been omitted.
Funding: This research was funded by the US Department of Veterans Affairs (ZAA). The funders had no role in considering the study design, or in the collection, analysis, interpretation of data, writing of the report, or decision to submit the article for publication.
Competing interests: All authors have completed the ICMJE uniform disclosure form at www.icmje.org/coi_disclosure.pdf and declare: support from the US Department of Veterans Affairs for the submitted work. ZAA reports receiving consultation fees from Gilead Sciences and receipt of funding (unrelated to this work) from Tonix pharmaceuticals. ZAA and YX reports consulting (uncompensated) for Pfizer; no other relationships or activities that could appear to have influenced the submitted work.
The lead author (ZAA) affirms that the manuscript is an honest, accurate, and transparent account of the study being reported; that no important aspects of the study have been omitted; and that any discrepancies from the study as planned have been explained.
Dissemination to participants and related patient and public communities: The authors have no specific plans to disseminate the findings to study participants. However, study results will be disseminated via a press release with mainstream media and will also be promoted on Twitter (via the corresponding senior author’s account @zalaly). We will also share the findings with public health agencies and our contacts in the US government (the US covid-19 task force), British government, and the Canadian government (the Office of the Chief Science Advisor of Canada). The results will be shared with advocacy groups and will be presented in scientific meetings and at academic institutions.
Provenance and peer review: Not commissioned; externally peer reviewed.
This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.