Abstract
Rationale Outcomes in patients with severe SARS-CoV-2 infection (COVID-19) are conditioned by virus clearance and regulation of inflammation. Variants in IFIH1, a gene coding the cytoplasmatic RNA sensor MDA5, regulate the response to viral infections.
Objective To characterize the impact of IFIH1 rs199076 variants on host response and outcomes after severe COVID-19.
Methods Patients admitted to an intensive care unit (ICU) with confirmed COVID-19 were prospectively studied and rs1990760 variants determined. Peripheral blood gene expression, cell populations and immune mediators were measured. Peripheral blood mononuclear cells from healthy volunteers were exposed to an MDA5 agonist and dexamethasone ex-vivo, and changes in gene expression assessed. ICU discharge and hospital death were modelled using rs1990760 variants and dexamethasone as factors in this cohort and in-silico clinical trials.
Measurements and Main Results 227 patients were studied. Patients with the IFIH1 rs1990760 TT variant showed a lower expression of inflammation-related pathways, an anti-inflammatory cell profile and lower concentrations of pro-inflammatory mediators. Cells with TT variant exposed to a MDA5 agonist showed an increase in IL6 expression after dexamethasone treatment. All patients with the TT variant not treated with steroids (N=14) survived their ICU stay (HR 2.49, 95% confidence interval 1.29–4.79). Patients with a TT variant treated with dexamethasone (N=50) showed an increased hospital mortality (HR 2.19, 95% confidence interval 1.01–4.87) and serum IL-6. In-silico clinical trials supported these findings.
Conclusions COVID-19 patients with the IFIH1 rs1990760 TT variant show an attenuated inflammatory response and better outcomes. Dexamethasone may reverse this anti-inflammatory phenotype.
Introduction
The spectrum of disease after infection by SARS-CoV-2 (COVID-19) may range from mild respiratory symptoms to a severe form of lung injury fulfilling acute respiratory distress syndrome (ARDS) criteria (1). In these severe cases, systemic response to infection may be associated to multiorgan failure and death (2). Therefore, outcomes in critically ill COVID-19 patients are related not only to viral clearance, but also to preservation of homeostasis.
The mechanisms responsible for the development of severe forms of COVID-19 have not been fully elucidated, but non-adaptative inflammatory responses play a central role. The only strategies that have decreased mortality in this population, steroids (3) and blockade of the IL-6 pathway (4), aim to limit this exacerbated immune response to prevent organ dysfunction.
The human gene IFIH1, located in the reverse strand of chromosome 2, encodes MDA5, a helicase that acts as a cytoplasmatic virus receptor. After binding to a viral RNA strand, MDA5 interacts with a mitochondrial adapter (MAVS, mitochondrial antiviral signaling protein), triggering the transcription of type-1 interferon genes and ultimately the systemic inflammatory response. In humans, the rs1990760 polymorphism encodes a variant of the IFIH1 gene (NM_022168: c.2836G>A [p.Ala946Thr]) that has been related to different susceptibility to viral infections and autoimmune disorders (5). By regulation of IFN-dependent pathways, IFIH1 participates in a feedback loop that ultimately modulates viral clearance and host inflammatory responses. In experimental models of Coxsackie virus infection, TT variants in rs1990760 result in lower proinflammatory cytokine levels without a major reduction in viral clearance (6).
Although the role of this variant in Coronavirus infections has not been explored, it has been proposed that the rs1990760 TT variant could confer resistance to SARS-CoV-2 infection and that differences in allelic frequencies could explain the epidemiological features of the pandemic in different countries (7). We hypothesized that, once infection is established, the inflammatory response in severe COVID-19 patients could be conditioned by IFIH1 variants. To test this hypothesis, we prospectively followed a cohort of critically ill patients with confirmed infection by SARS-CoV-2, in which peripheral blood gene expression, cell populations, concentrations of immune mediators and clinical outcomes were studied and related to rs1990760 variants.
Methods
A detailed description of the methods is provided in the online supplement.
Study design
This single-center prospective, observational study was approved (ref. 2020/188) by the Clinical Research Ethics Committee of Principado de Asturias (Spain). Informed consent was obtained from each participant or next of kin. Given the exploratory nature of the study objective and the absence of previous data, no formal sample size calculations were performed.
All patients with confirmed SARS-CoV-2 infection (8) and meeting the Kigali modification (9) of ARDS criteria (to include patients without mechanical ventilation) from 16-Mar-2020 to 10-Dec-2020 were included in the study. Exclusion criteria were age<18, any condition that could explain the respiratory failure other than COVID-19, do-not-resuscitate orders or terminal status, or refusal to participate. The presence of SARS-CoV-2 in respiratory samples was detected by reverse transcriptase qPCR (Supplementary table E1) and clearance half-life (λViral Clearance) calculated (10).
Patients were followed until hospital discharge, and clinical and analytical data collected. The main outcome was ICU discharge alive and spontaneously breathing. Secondary outcome was hospital discharge.
Biochemical measurements
IFIH1 rs1990760 C/T polymorphism was determined by PCR genotyping in DNA extracted from peripheral blood. Blood samples were taken in days 1 (first 24 hours after ICU admission, before steroid treatment) and 4 of ICU stay. RNA was extracted and sequenced. The generated sequences were mapped and counted using Salmon v1.4 software (11). Raw counts were compared between groups using DESeq2 (12). Genes with differential expression (adjusted p-value<0.05, corrected using a false discovery rate of 0.05) were analyzed using Ingenuity Pathway Analysis (Qiagen, USA) to identify overrepresented gene sets and networks. Over the identified networks, in-silico effects of IFIH1 down-regulation and addition of exogenous dexamethasone were tested. Circulating cell populations were estimated from gene expression using a previously validated deconvolution algorithm (13). The correlation between estimated and measured lymphocyte percentages was 0.61 (Supplementary Figure E1).
In a second sample, a panel of serum inflammatory mediators (interferons (IFN)-β, -γ and -λ, tumor necrosis factor (TNF)-α, interleukins (IL)-1β and -6, and chemokines CXCL8, CXCL9, CXCL10, CXCL16, CCL2, CCL3, CCL4 and CCL7) was measured using a multiplexed assay. Additional serum IL-6 measurements were done during clinical follow-up.
Ex-vivo experiments
Peripheral blood mononuclear cells (PBMC) from healthy volunteers, genotyped for IFIH1 rs1990760, were isolated and cultured in presence of medium, medium plus a MDA5 ligand (high-molecular weight poly-I:C), or medium plus the MDA5 ligand and dexamethasone (Kern Pharma, Spain). After 24h cells were collected and homogenized, RNA extracted and expression of STAT1, STAT3, FOXO3, IL6 and GAPDH quantified (Supplementary Table E2).
Clinical trial simulations
Data from the RECOVERY trial were extracted and used to estimate risk ratios (RR) for each rs1990760 variant, with and without steroids (see online supplement for details). With these data, a survival model was developed and used to simulate scenarios with different allelic frequencies and baseline risks of death. Hazard ratios (HR) for 28-day mortality were calculated from these simulations.
Statistical analysis
Data are shown as median (interquartile range). Comparisons between IFIH1 variants were done using Wilcoxon or ANOVA tests, and p-values corrected using a false discovery rate of 0.05. Results from the ex-vivo experiments were fitted to a mixed-effects model including donor as a fixed effect and experimental group and genotype as random effects. Survival was analyzed using competing risks models considering ICU/hospital discharge alive and spontaneously breathing and ICU/hospital death as different terminal events. IFIH1 variants and steroid treatment were added in these models as covariables with an interaction, thus resulting in 4 study groups. Patients with rs1990760 CC/CT variants not treated with steroids were considered the reference category in all the analyses. All the analyses and plots were performed using R 4.0.1 (14).
Results
Study cohort
250 patients were admitted in ICU due to suspected or confirmed COVID-19. Among these, 227 were included in the study. The study flow chart and reasons for exclusion are shown in Figure 1. Basic demographic and clinical data are shown in table 1. IFIH1 rs1990760 variants in this population met Hardy-Weinberg conditions (53 CC, 110 CT, 64 TT, chi-square p=0.19). There were no differences in comorbidities or clinical data at admission between genotypes other than a lower PaO2/FiO2 ratio in patients with a TT variant (Table 1).
IFIH1 rs1190760 variants and inflammatory response to SARS-CoV-2
Gene expression in peripheral blood during the first day of ICU admission was profiled in 42 patients who did not receive steroid therapy at that time (11, 19 and 12 with rs1990760 TT, CT and CC genotypes respectively). Expression of IFIH1 was significantly lower in patients with the TT genotype (Figure 2A and Supplementary Figure E2). Comparison of peripheral blood gene expression between patients with TT and CT/CC variants yielded significant differences in 179 genes (Figure 2B-C, Supplementary File 1). Visual inspection of the heatmap reveals that differences are quantitative rather than qualitative. Ingenuity pathway analysis revealed several gene networks involved in the regulation of the inflammatory response among the differentially expressed genes (Figure 2D and Supplementary Figure E3). In-silico predictions suggests that IFIH1 downregulation (such as in patients with the rs1990760 TT variant) decreases expression of proinflammatory pathways (Supplementary Figure E4).
Then we assessed the immunological consequences of these differences in gene expression and cell populations. Patients with the TT genotype have lower inferred percentages of classic CD14+ monocytes and circulating plasma cells, and an increase in hematopoietic precursors, myeloid dendritic cells, M2 macrophages and CD56dim NK cells (Figure 3).
Serum cytokines were measured at ICU admission in 28 patients (8, 10 and 10 with TT, CT and CC genotypes respectively). There were no differences in any of the measured interferons. Patients with the TT genotype showed lower levels of proinflammatory mediators including IL-6, CXCL10, CXCL16, and CCL7 (Figure 4).
Ex-vivo effects of dexamethasone in IFIH1 variants
During the study period, dexamethasone was added to COVID-19 treatment based on results from published trials (3). An in-silico analysis focused on the interactions between IFIH1 and steroids suggested that dexamethasone could disrupt some of the effects caused by IFIH1 downregulation (Supplementary Figure E5). Specifically, dexamethasone may change expression of several IFIH1-dependent genes, including STAT1, STAT3 or FOXO3 among others.
To test these predictions, an ex-vivo experiment using peripheral blood mononuclear cells from healthy volunteers with different IFIH1 rs1990760 variants (n=5, 6 and 7 for CC, CT and TT variants respectively) was performed. Cells were collected and exposed to poly I:C, to mimic SARS-Cov-2 infection, and dexamethasone. Compared to exposure to poly I:C alone, dexamethasone had no effect on STAT1 (Figure 5A) or STAT3 (Figure 5B) expression. However, the steroid increased the expression of FOXO3 (Figure 5C) and IL6 (Figure 5D) only in cells with the TT variant. These results suggest that dexamethasone may alter the inflammatory response triggered by MDA5 activation in those patients with the TT variant of the gene.
Impact of IFIH1 variants on outcomes in critically ill COVID-19 patients
Median follow-up was 25 days (interquartile range 16 – 39 days) and 27 days in survivors (interquartile range 17 – 44 days). There were no differences in ICU length of stay among groups (Supplementary table E2). ICU survival and hospital mortality were modelled using IFIH rs1990760 genotype and steroid treatment as interacting covariables. There were no significant differences in the main clinical characteristics among the resulting groups (Supplementary Table E2). Regarding ICU survival (Figure 6A and Supplementary Table E3), 27 out of 35 patients with CC/CT alleles who did not receive steroids were discharged alive and spontaneously breathing from the ICU (HR 1, used as reference). In patients with this variant, dexamethasone was not related to a better outcome (98 out of 128 patients discharged, HR 1.20 [0.78 – 1.38], p=0.41). All patients with the TT allele (n=14) who did not receive steroids were discharged alive (HR 2.49 [1.29 – 4.79], p=0.012). Steroid treatment in patients with the TT variant was related to the loss of this benefit (32 out of 50 patients discharged alive, HR=1.03 [0.62 – 1.72], p=0.91).
Hospital mortality is shown in Figure 6B and Supplementary Table E3. There were 9 hospital deaths in patients with the CT/CC genotypes not treated with steroids (out of a total of 35 patients), and 27 out of 128 hospital deaths when treated with steroids, resulting in a HR of 1.11 [0.52 – 2.37] (p=0.80). All patients with the TT allele who were not treated with steroids survived after their hospital stay (HR 1.23 × 10−7, confidence intervals and p value cannot be computed due to the absence of events). Patients with the TT allele and treated with steroids showed a worse outcome (19 deaths in 50 patients, HR 2.19 [1.01 – 4.87], p=0.05).
To investigate the causes responsible for the differences in mortality among groups, we first quantified viral load at diagnosis (Figure 6C) and clearance after its peak value (Figure 6D). There were no differences in these parameters related to IFIH1 variants or treatment groups. We also compared serum IL-6 levels at ICU admission and one week later. Compared to patients with CC/CT variants not treated with steroids, patients with the same genotype but receiving dexamethasone and patients with the TT variant not receiving steroids showed lower levels of IL-6 at admission (Figure 6E). After one week (Figure 6F), serum IL-6 levels remained lower in patients receiving dexamethasone with the CC/CT variant and in those with the TT variant and not treated with this drug. However, IL-6 levels were higher in those with the TT variant and treated with steroids, in line with our ex-vivo findings.
To explore the mechanisms behind these differences, peripheral blood gene expression at ICU day 4 was compared between patients treated with or without steroids for each rs1990760 variant. IFIH1 expression decreased in patients with rs1990760 CC/CT variants treated with steroids, but not in those with a TT variant (Figure 6G). In opposite, steroids increased FOXO3 expression only in patients with a TT variant (Figure 6H), resembling the results from the ex-vivo experiments. IL6 gene raw counts at day 4 were below 5 in all the patients, and thus not compared. When the whole transcriptomes were compared, steroids significantly changed the expression of 58 genes in patients with a CC/CT variant and 23 in patients with a TT variant. Overall, changes in gene expression were qualitatively different between variants, with only three genes in common (Figure 6I and Supplementary Figure E6).
In-silico clinical trials
These findings may support the hypothesis that COVID-19 populations with a low proportion of the rs1990760 TT variant would show a better response to dexamethasone. To explore this prediction, data from the RECOVERY trial (15) were examined. The T-allele frequency in populations from a Black/Asian ancestry is 0.13 (16), so only 1.7% has a TT variant. In these populations, RR was 0.7 (0.51-0.95), whereas in patients from a White ancestry (with a T allele frequency of 0.61), RR increased to 0.9 (0.8-1.2).
Using these data, RRs related to each rs1990760 variant, with and without steroids, were extracted and a survival model developed (see online supplement for details). A simulation including 500 patients from each rs1990760 variant and treatment arm (placebo or dexamethasone) showed a significant reduction in mortality with steroids in patients with a CC/CT variant, but not in those with a TT variant (Figure 7A). In patients not-receiving steroids, HRs related to a TT variant, compared to CC/CT variants, were independent of the allelic frequency (Figure 7B). Simulations of steroid therapy with different allelic frequencies showed that HRs in each specific variant remained constant whereas the effect of steroids in the overall population was dependent on the allele distribution (Figure 7C and Supplementary Figure E7).
Discussion
Our results show that critically ill patients with the TT variant in the IFIH1 rs1990760 polymorphism have an attenuated inflammatory response to severe SARS-CoV-2 infection, leading to a decreased mortality. In this selected population, treatment with steroids has no immunomodulatory effects and could be related to worse outcomes. These results confirm the impact of the host response on patients’ outcomes and suggest that patient geno/phenotypes should be taken into account to prescribe steroids in this setting.
Host response is a major determinant of outcome in critically ill patients, including those with COVID-19. Several genetic variants involved in the inflammatory response have been related to SARS-CoV-2 infection and its severity (17–19). Activation of immune responses causes local and systemic inflammation, aimed to block viral replication. However, exacerbation of these responses can cause organ damage even after virus clearance. Our results recapitulate previous findings in COVID-19 describing the release of proinflammatory mediators, NK cell exhaustion, monocyte dysregulation and emergency hematopoiesis (20–22). Some of our findings in patients with the rs1990760 TT variant, including lower levels of circulating proinflammatory molecules and a shift towards anti-inflammatory cell populations (hematopoietic precursors, M2 macrophages or CD56dim NK cells) have been linked to better outcomes in other observational studies (23). Collectively, our findings raise the hypothesis that the TT variant could be related to a lower IFIH1 expression and an attenuated proinflammatory response. In line with this observation, a TT variant has been associated to lower proinflammatory cytokine levels in patients with lupus (24, 25) and in experimental viral infections (6).
Development of tolerance to viral diseases can be considered a major evolutionary adaptative response (26), and inhibition of pro-inflammatory responses has been proposed to improve the outcome of severe COVID-19 patients. Notably, pangolins, an intermediate host of coronaviruses (27), lack functional IFIH1. It has been suggested that this deficiency reduces the inflammatory response to coronavirus infections (28). IL-6 blockade (4) or steroids (3) are the only treatments that have improved the outcome of critically ill COVID-19 patients. However, mortality rates are still around 27-32% (3, 4), so new therapeutic approaches are warranted.
According to our results, the impact of the treatment with steroids along a given population may depend on intrinsic individual characteristics, allowing a personalized approach based on a specific genomic biomarker. In our study sample, patients with the IFIH1 rs1990760 TT variant constitute a population with a better prognosis, in whom treatment with dexamethasone may be reconsidered as it was associated to higher mortality rates. Indeed, as shown by the in-silico trials, any population enriched for patients with CC/CT variants in rs1990760 will show higher mortality rates and a better response to steroids, due to the low proportion of TT variants. Dexamethasone in patients receiving mechanical ventilation at inclusion in the RECOVERY trial had a HR of 0.64, compared to a HR of 0.83 in the whole population. Of note, the proportion of Black/Asian patients (with a low frequency of a T allele) in the mechanically ventilated group was 29%, but only 18% in the whole study population. To date, no other genomic markers of personalized therapies in critically ill patients have been identified. The absence of differences in viral clearance and the late increase in IL-6 in these patients suggest that this worse outcome is more related to the disruption of ongoing immunoregulatory mechanisms than to antiviral responses. Our ex-vivo experiments do not fully elucidate the molecular mechanisms behind the interaction between TT variants and steroids, as the effects of FOXO3 upregulation and IL6 expression may be variable (29), but clearly illustrate that this specific combination may disrupt the ongoing self-regulation of inflammation. Moreover, steroid-induced changes in gene expression are qualitatively different in patients of each variant.
Our findings point towards the therapeutic potential of MDA5 modulation in COVID-19, either induced by steroids or targeted by other drugs. Amelioration of MDA5-dependent RNA sensing could avoid an exacerbated inflammatory response without impairing viral clearance. In this sense, it has been described that the interferon response triggered by MDA5 is unable to control viral replication (30). However, genetic ablation of IFIH1 results in increased viral loads and decreased cytokine production (31), so this approach must be viewed with caution.
It is unclear if these findings can be translated to other viral diseases. Steroids may decrease mortality in an unselected ARDS population (32), but the role of underlying genotypes has not been addressed. An experimental model of Coxsackie virus infection revealed lower IFIH1 and CXCL10 expression in cells with the TT variant (6), with no relevant differences in viral clearance. However, total absence of MDA5 results in impaired clearance of West Nile virus (33).
Our work has several limitations. First, the results must be validated in an independent cohort. Our in-silico simulations reinforce the external validity of our findings, but a pharmacogenomic analysis of patients included in clinical trials is warranted for confirmation. Second, steroid treatment was not randomized, so we cannot discard other underlying factors responsible for the observed differences. Although there could be confounding by indication, there were no baseline differences among groups suggesting a higher severity in steroid-treated patients. Moreover, there are significant differences between genotypes irrespective of the treatment received. In addition, in-silico and ex-vivo experiments are congruent with the observed clinical results, supporting the differential effects of steroid therapy in IFIH1 rs1990760 variants. Third, the favourable outcome of patients with a TT allele without steroids is based on a small sample size. As steroids are now the standard of care for severe COVID-19, this sample size cannot be increased outside a hypothetical clinical trial focused of personalized steroid therapy according to IFIH1 rs1990760 variants. However, similar reduced numbers have served to identify other variants in the immune response (34), and the related finding of increased mortality in this genotype after steroid therapy compared to all other variants is supported by a larger sample. In addition, our sample is representative of European population with a limited racial diversity that may have influenced our results, as the in-silico analyses further suggest. Finally, other studies (35, 36) have focused on the genetic variants linked to an increased risk of severe COVID-19, compared to non-infected populations. However, no genetic markers have been associated to mortality in cohorts of infected patients. According to our data, rs1990760 is linked to the outcome, but no inferences on susceptibility to severe COVID-19 can be extracted.
In summary, we have identified a genetic variant of IFIH1 that results in decreased expression of this gene and an ameliorated inflammatory response after severe SARS-CoV-2 infection. Patients with the rs1990760 TT genotype show a good outcome. However, this adaptative response was not observed in patients with a TT variant and treated with steroids. These findings suggest that the systemic response to severe COVID-19 is regulated by genetic factors that modulate the response to the infection and the prescribed therapy and, ultimately, may impact the outcome.
Data Availability
RNA-seq data has been deposited in Gene Expression Omnibus (accession numbers GSE168400 and GSE177025). Code for in-silico clinical trials is available at https://github.com/Crit-Lab/IFIH1_simulation. De-identified clinical data will be shared with researchers with any proposal after review and approval by the local Ethics Committee and a signed data access agreement. R code used for statistical analyses is also available on request. All requests should be sent to the corresponding author (gma{at}crit-lab.org).
https://github.com/Crit-Lab/IFIH1_simulation
https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE168400
https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE177025
Data sharing
RNA-seq data has been deposited in Gene Expression Omnibus (accession numbers GSE168400 and GSE177025). Code for in-silico clinical trials is available at https://github.com/Crit-Lab/IFIH1_simulation. De-identified clinical data will be shared with researchers with any proposal after review and approval by the local Ethics Committee and a signed data access agreement. R code used for statistical analyses is also available on request. All requests should be sent to the corresponding author (gma{at}crit-lab.org).
Acknowledgements
The authors want to thank all the personnel at the participating ICUs and laboratories for their support during the development of the study. In addition, they thank Salvador Villalgordo, Silvia Viñas, Salvador Balboa, Rodrigo Albillos, Cecilia del Busto, Emilio García-Prieto, Jose Antonio Gonzalo, Diego Parra, Lisardo Iglesias, José Alonso, Sérida Dominguez and Alfredo González their additional efforts for sample collection.
Footnotes
Funding Supported by Centro de Investigación Biomédica en Red (CB17/06/00021), Instituto de Salud Carlos III (grant numbers PI19/00184 and PI20/01360) and Fundació La Marató de TV3 (grant 413/C/2021). Instituto Universitario de Oncología is supported by Fundación Liberbank. RRG is the recipient of a grant from Instituto de Salud Carlos III (CM20/00083).
Declaration of interests None of the authors report a conflict of interest.