Research Article

Effects of IFIH1 rs1990760 variants on systemic inflammation and outcome in critically ill COVID-19 patients in an observational translational study

Unidad de Cuidados Intensivos Cardiológicos. Hospital Universitario Central de Asturias, Spain
Centro de Investigación Biomédica en Red (CIBER)-Enfermedades Respiratorias. Instituto de Salud Carlos III, Spain
Instituto de Investigación Sanitaria del Principado de Asturias, Spain
Instituto Universitario de Oncología del Principado de Asturias, Spain
Unidad de Cuidados Intensivos Polivalente. Hospital Universitario Central de Asturias, Spain
Red de Investigación Renal (REDINREN), Spain
Servicio de Inmunología. Hospital Universitario Central de Asturias, Spain
Area de Gestión Clínica de Pediatría. Hospital Universitario Central de Asturias, Spain
Department of Anesthesiology and Operative Intensive Care Medicine CCM/CVK, Charité - Universitätsmedizin Berlin, corporate member of Freie Universität Berlin, Humboldt Universität zu Berlin, and Berlin Institute of Health, Germany
Servicio de Genética Molecular. Hospital Universitario Central de Asturias, Spain
Servicio de Microbiología. Hospital Universitario Central de Asturias, Spain
Universidad de Oviedo, Spain
Servicio de Informática. Hospital Universitario Central de Asturias, Spain
Critical Care Center, Hospital Universitari Parc Taulí, Institut D'Investigació I Innovació Parc Taulí I3PT, Universitat Autònoma de Barcelona, Spain
Servicio de Bioquímica Clínica. Hospital Universitario Central de Asturias, Spain

Jan 21, 2022

https://doi.org/10.7554/eLife.73012

Open access
Copyright information

Version of Record

Accepted for publication after peer review and revision.

Download
Cite
Share
CommentOpen annotations (there are currently 0 annotations on this page).

Version of Record published: January 21, 2022 (This version)
Accepted: December 6, 2021
Received: August 12, 2021
Preprint posted: July 7, 2021 (Go to version)

1. Part of Collection
COVID-19: A Collection of Articles

Edited by Diane M Harper et al.
Further reading

Abstract
Editor's evaluation
eLife digest
Introduction
Materials and methods
Results
Discussion
Appendix 1
Data availability
References
Article and author information
Metrics

Abstract

Background:

Variants in IFIH1, a gene coding the cytoplasmatic RNA sensor MDA5, regulate the response to viral infections. We hypothesized that IFIH1 rs199076 variants would modulate host response and outcome 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 modeled using rs1990760 variants and dexamethasone as factors in this cohort and in-silico clinical trials.

Results:

About 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 an MDA5 agonist showed an increase in IL6 expression after dexamethasone treatment. All patients with the TT variant not treated with steroids survived their ICU stay (hazard ratio [HR]: 2.49, 95% confidence interval [CI]: 1.29–4.79). Patients with a TT variant treated with dexamethasone showed an increased hospital mortality (HR: 2.19, 95% CI: 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.

Funding:

Centro de Investigación Biomédica en Red (CB17/06/00021), Instituto de Salud Carlos III (PI19/00184 and PI20/01360), and Fundació La Marató de TV3 (413/C/2021).

Editor's evaluation

Authors present an interesting approach of COVID-19 multiorgan failure and loss of homeostasis attributed to modified adaptive immune response due to presence/absence of MDA5 polymorphism(rs1990760). Data presented by authors are very interesting representing novel scientific work that adds new and potentially pivotal information to existing knowledge. Key data presented strongly support an individualized approach during COVID-19 pandemic that could eventually alter therapeutic algorithm regarding dexamethasone administration implemented for these patients.

https://doi.org/10.7554/eLife.73012.sa0

eLife digest

Patients with severe COVID-19 often need mechanical ventilation to help them breathe and other types of intensive care. The outcome for many of these patients depends on how their immune system reacts to the infection. If the inflammatory response triggered by the immune system is too strong, this can cause further harm to the patient.

One gene that plays an important role in inflammation is IFIH1 which encodes a protein that helps the body to recognize viruses. There are multiple versions of this gene which each produce a slightly different protein. It is possible that this variation impacts how the immune system responds to the virus that causes COVID-19.

To investigate, Amado-Rodríguez, Salgado del Riego et al. analyzed the IFIH1 gene in 227 patients admitted to an intensive care unit in Spain for severe COVID-19 between March and December 2020. They found that patients with a specific version of the gene called TT experienced less inflammation and were more likely to survive the infection.

Physicians typically treat patients with moderate to severe COVID-19 with corticosteroid drugs that reduce the inflammatory response. However, Amado-Rodríguez, Salgado del Riego et al. found that patients with the TT version of the IFIH1 gene were at greater risk of dying if they received corticosteroids.

The team then applied the distribution of IFIH1 variants among different ethnic ancestries to data from a previous clinical trial, and simulated the effects of corticosteroid treatment. This ‘mock’ clinical trial supported their findings from the patient-derived data, which were also validated by laboratory experiments on immune cells from individuals with the TT gene.

The work by Amado-Rodríguez, Salgado del Riego et al. suggests that while corticosteroids benefit some patients, they may cause harm to others. However, a real-world clinical trial is needed to determine whether patients with the TT version of the IFIH1 gene would do better without steroids.

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 (Grasselli et al., 2020). In these severe cases, systemic response to infection may be associated with multiorgan failure and death (Du et al., 2020). 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 (WHO Rapid Evidence Appraisal for COVID-19 Therapies (REACT) Working Group et al., 2020) and blockade of the IL-6 pathway (REMAP-CAP Investigators et al., 2021), 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 (Gorman et al., 2017). 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 pro-inflammatory cytokine levels without a major reduction in viral clearance (Domsgen et al., 2016).

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 (Maiti, 2020). 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.

Materials and methods

Study design

Request a detailed protocol

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. The study started after approval from the ethics committee and finished in December 2020, after the second pandemic wave.

All patients with confirmed SARS-CoV-2 infection (Escudero et al., 2020) and meeting the Kigali modification (Riviello et al., 2016) of ARDS criteria (to include patients without mechanical ventilation) from March 16, 2020 to December 10, 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. Patients were followed until hospital discharge, and clinical and analytical data collected. The main outcome was intensive care unit (ICU) discharge alive and spontaneously breathing. Secondary outcome was hospital discharge.

All biochemical analyses (cytokine measurements, RNAseq, estimation of cell populations…, see below) were performed by researchers blinded to genotype and outcome.

SARS-CoV-2 detection and quantification

Request a detailed protocol

The presence of SARS-CoV-2 was analyzed by detecting viral genome using a multiple quantitative retrotranscriptase (RT)-PCR. Nucleic acids were purified by MagNa Pure 96 System (Roche, Geneva, Switzerland) from the swabs transport medium. The extracts were subjected to an amplification reaction using TaqMan Fast 1-Step Master Mix (Life Technologies, Carlsbad, CA) supplemented with a mixture of primers (Thermo Fisher Scientific, Walthan, MA) and TaqMan MGB probes (Applied Biosystems, Foster City, CA) directed against ORF1ab and N genes (Supplementary file 1a). Amplifications and subsequent analysis were carried out using the Applied Biosystems 7500 Real-Time PCR System (Escudero et al., 2020). Amplification of viral genes with a Ct number lower than 35 was considered positive. Viral load was normalized by the number of cells and expressed as copies/1000 cells, as previously described (Gómez-Novo et al., 2018). Clearance of SARS-CoV-2 was evaluated by fitting viral load in tracheobronchial samples after its peak value over time using an exponential decay function and calculation of viral clearance half-life (λ_{Viral Clearance}).

Genotyping

Request a detailed protocol

DNA was extracted from total blood leukocytes (1 ml) with an automated equipment (Promega Maxwell). IFIH1 rs1990760 C/T polymorphism determined by real-time PCR with TaqMan genotyping master mix (Life Technologies) and TaqMan probes (Thermo Fisher Scientific, assay C_2780299_30) in an ABI-7500 device. The genotyping strategy was validated by Sanger sequencing of selected individuals from the three genotypes. Based on previous data showing that patients with a TT genotype have a different immune response compared to CC and CT (Wawrusiewicz-Kurylonek et al., 2020; Zhang et al., 2018), patients with these two variants were grouped and compared against the TT group.

Blood sampling

Request a detailed protocol

Two blood samples were taken in the first 24 hr after ICU admission. No patient had received steroids at the time of sampling. About 3 ml of blood was collected in Tempus tubes (Thermo Fisher Scientific) for RNA isolation and immediately stored at –80°C. Additional 5 ml was collected in a Vacutainer serum tube (BD Biosciences), and isolated by centrifugatio, and stored at –80°C until analysis.

RNAseq

View detailed protocol

After thawing, blood from Tempus tubes was diluted in phosphate-buffered saline (PBS) (1:3 v/v) and centrifuged at 3000 rpm for 30 min. Supernatant was discarded and the pellet resuspended in TRIzol (Sigma-Aldrich, Poole, UK) and precipitated overnight with isopropanol at –20°C. After centrifugation, RNA pellets were washed with 70% ethanol and resuspended in RNAse-free water. RNA quality was assessed using a TapeStation, and only samples with a RIN (RNA integrity number) above 8 were analyzed.

RNA sequencing was performed using Ion AmpliSeq Transcriptome Human Gene Expression Kit, in an Ion S5 GeneStudio sequencer (Ion Torrent). Briefly, 10 ng of total RNA was retrotranscripted and the obtained cDNA used for library synthesis using Ion AmpliSeq Transcriptome kits to amplify all the canonical human transcripts. After template preparation in an automated Ion Chef Instrument, semiconductor 540 chips were run in an Ion S5 GeneStudio sequencer. Torrent Suite software was used for base calling, alignment, and sequence quality controls. The generated FASTQ files (available at Gene Expression Omnibus, accession numbers GSE168400 and GSE 177025) were mapped against a reference transcriptome (obtained from http://refgenomes.databio.org) and transcripts counted using Salmon v1.4 software (Patro et al., 2017). Raw counts were compared between genotypes using the DESeq2 library (Love et al., 2014). The log₂-fold change between variants for each gene and the adjusted p-value (corrected using a false discovery rate of 0.05) were calculated and analyzed using Ingenuity Pathway Analysis (Qiagen, USA) to identify overrepresented gene sets and networks. Over the identified network, in-silico effects of IFIH1 downregulation and addition of exogenous dexamethasone were performed.

Peripheral blood cell populations

Request a detailed protocol

Circulating cell populations were estimated from gene expression using a previously validated deconvolution algorithm (Vallania et al., 2018). Using a reference expression matrix, proportions of 20 different cell lines were calculated. Only cell lines present in more than five samples were considered. As a quality check, we analyzed the correlation between estimated and measured lymphocyte percentages. The obtained correlation coefficient was 0.61 (Figure 3—figure supplement 1). It must be noted, however, that the deconvolution method estimates proportions over transcriptionally active cells, which may not be equivalent to the obtained cell count (as there may be inactive cells in the latter).

Inflammatory mediators

Request a detailed protocol

A panel of inflammatory mediators was studied in serum from patients not receiving steroids during their first ICU day. Serum concentrations of interferons (IFN)-β, -γ and -λ, tumor necrosis factor (TNF)-α, interleukins (IL)-1β and -6, and chemokines CXCL8, CXCL9, CXCL10, CXCL16, CCL2, CCL3, CCL4, and CCL7 were measured using a multiplexed assay (Luminex custom panel). Concentrations below the lower limit of detection for a given mediator were replaced with half of that limit.

Ex-vivo experiments

Request a detailed protocol

To study the potential interferences between IFIH1 rs1990760 variants and dexamethasone predicted by in-silico analyses, an ex-vivo experiment was designed. Blood samples from healthy volunteers (genotyped for IFIH1 rs1990760 variants using DNA obtained from buccal swabs) were collected in EDTA tubes and immediately processed. Peripheral blood mononuclear cells (PBMCs) were isolated via density-gradient centrifugation with Lymphoprep (Axis-Shield PoC AS, Oslo, Norway). Cells were washed with red blood cell lysis buffer (NH₄Cl 0.1 M, KHCO₃ 0.01 M, EDTA 0.1 mM in dH₂O, pH 7.33) and PBS before being resuspended in RPMI-1640 (Gibco, USA) + 10% fetal bovine serum (FBS). Cells were seeded in 12-well plates at a final concentration of 1.5×10⁶ PBMC/ml and cultured at 37°C and 5% CO₂ in presence of medium/FBS, medium/FBS plus a MDA5 ligand (1 µg/ml high-molecular weight poly-I:C/LyoVec, Invivogen, USA), or medium plus the MDA5 ligand and 10 µM dexamethasone (Kern Pharma, Spain). Although poly-I:C may bind to both RIG-I and MDA-5, long (high-molecular weight) poly-I:C binds specifically to MDA-5 (Kato et al., 2008). After 24 hr, cells were collected and homogenized in TRIzol for RNA extraction. 500 ng of total RNA was retrotranscribed into cDNA using an RT-PCR Kit (High-capacity cDNA rt Kit, Applied Biosystems, USA). Expression of STAT1, STAT3, FOXO3, IL6, and GAPDH was quantified using 5 ng of cDNA per well and in triplicate for each sample. SYBR Green Power up (Thermo Fisher Scientific) and 10 µM of the corresponding primers (Supplementary file 1b) were used in all the experiments. The relative expression of each gene was calculated as 2^{−ΔCT(gene of interest) – ΔCT(GAPDH)}.

Clinical trial simulations

Request a detailed protocol

Data from the RECOVERY trial were extracted and used to estimate risk ratios (RRs) for each rs1990760 variant, with and without steroids (see online Appendix 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 (HRs) for 28-day mortality were calculated from these simulations.

Statistical analysis

Request a detailed protocol

Data are shown as median (interquartile range). All data points and sample sizes represent individual patients/biological replicates. 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 model were fitted to a mixed-effects model including experimental group and genotype as covariables, and post hoc comparisons evaluated using Holm’s correction. No outliers were discarded.

Survival was analyzed using competing risks, Cox regression model, with ICU/hospital discharge alive and spontaneously breathing and death as competing risks, using the Aalen-Johansen estimator, as previously described (Barbaro et al., 2020). This competing-risks framework is needed as patients discharged alive have a low probability of death, so censoring at the time of discharge using a standard Kaplan-Meier approach would lead to biased observation, as the probability of death is different than in those still followed (i.e., informative censoring). Patients with an rs1990760 CC/CT variant not treated with steroids were considered the reference category in all the analyses.

All the analyses and plots were performed using the R 4.0.1 statistical environment (R Development Core Team, 2020) with the packages data.table (Dowle and Srinivasan, 2021), multcomp (Hothorn et al., 2008), survival (Therneau, 2020), MetaIntegrator (Haynes et al., 2017), and ggplot2 (Wickham, 2016, p. 2).

Results

Study cohort

About 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 PaO₂/FiO₂ ratio in patients with a TT variant (Table 1).

Figure 1

Download asset Open asset

Patient study flow chart.

From 250 patients admitted to the participant ICUs during the study period, 227 were included in the study and classified according to the *IFIH1* rs1990760 variants and prescription of steroid therapy. Peripheral blood samples for gene expression (N=41) and quantification of immune mediators (N=28) were taken in the first day after ICU admission in patients not receiving steroids. All patients were followed up until death or hospital discharge.

Table 1

Clinical characteristics of the study cohort.

COPD: Chronic Obstructive Pulmonary Disease. Values are shown as absolute count or median (interquartile range). PBW: Predicted body weight. NIV: Non-invasive ventilation. PEEP: Positive End-Expiratory Pressure. *p-values calculated for proportion over the number of intubated patients.

	rs1990760
	Overall (n=227)	TT (n=65)	CC/CT (n=162)	p-value
Demographics
Age (years)	67 (60–75)	68.5 (63–75.25)	66 (58–75)	0.123
Sex
Male	174 (77)	47 (72)	127 (78)	0.327
Female	53 (23)	18 (28)	35 (22)
Race
Black	3 (1.5)	1 (2)	2 (1.5)	0.527
White	207 (91)	62 (95)	146 (90)
Latino	15 (7)	2 (3)	13 (8)
Asian	1 (0.5)	0	1 (0.5)
Body mass index (kg/m²)	30 (27–33)	30 (26–33)	30 (27–32)	0.592
Days since symptom onset	8 (6–11)	8 (6–10)	8 (6–11)	0.873
Days from hospital admission to ICU admission	2 (1–4)	2 (1–3)	2 (1–4)	0.841
APACHE-II score	15 (12–18)	15.5 (13–18)	14 (11—18)	0.177
Comorbidities
Arterial hypertension	132 (58)	38 (59)	94 (58)	0.932
Diabetes	51 (23)	14 (22)	37 (23)	1
Chronic kidney disease	15 (7)	5 (8)	10 (6)	0.743
COPD	17 (8)	7 (11)	10 (6)	0.339
Cirrhosis	2 (1)	1 (2)	1 (0.5)	1
Neoplasms
No	215 (95)	62 (95)3 (5)0	154 (95)	0.669
Active	10 (4)	3 (5)	7 (3.5)
Past	2 (1)	0	2 (1.5)
Immunosuppresive drugs (incl. steroids)	5 (2)	3 (5)	2 (1.5)	0.142
Lung function at ICU admission
Ventilation at admission
Spontaneous /NIV	38 (17)	11 (17)	27 (17)	0.278
Controlled invasive ventilation	187 (82.5)	53 (82)	135 (83)
Pressure support ventilation	1 (0.5)	1 (2)	0
FiO₂	0.5 (0.4–0.6)	0.5 (0.4–0.7)	0.5 (0.4–0.6)	0.104
PaO₂/FiO₂ (mmHg)	204 (155–267)	185 (142–226)	218 (162–282)	0.007
PaCO₂ (mmHg)	43 (39–47)	44 (39–47)	43 (39–48)	0.924
Respiratory rate (min^–1)	18 (16–20)	18 (16–20)	18 (16–22)	0.929
Arterial pH	7.38 (7.32–7.41)	7.38 (7.33–7.42)	7.37 (7.32–7.41)	0.805
Tidal volume/PBW (ml/Kg)	7.6 (6.9–8.4)	8 (6.9–9.2)	7.5 (7–8.1)	0.109
Plateau pressure (cmH₂O)	24 (21–28)	24 (22–29)	24 (21–27)	0.221
PEEP (cmH₂O)	12 (10–14)	12 (10–14)	12 (10–14)	0.397
Driving pressure (cmH₂O)	12 (10–14)	12 (10–14)	12 (10–14)	0.496
Respiratory system compliance (ml/cmH₂O)	38 (32–50)	36 (31–44)	42 (32–50)	0.11
Laboratory results
Leukocytes (×10³/μl)	8.34 (6.05–11.38)	8.77 (5.96–11.52)	8.29 (6.08–11.18)	0.773
Lymphocytes (×10³/μl)	0.65 (0.48–0.93)	0.62 (0.47–0.94)	0.67 (0.49–0.93)	0.266
Serum creatinine (mg/dl)	0.83 (0.62–1.13)	0.87 (0.62–1.14)	0.82 (0.62–1.11)	0.859
Serum ferritin (ng/ml)	1107 (717–1667)	965 (634–1263)	1262 (776–2374)	0.011
D-dimer (ng/ml)	1053 (681–1967)	1085 (643–1805)	1036 (683–2126)	0.854
Additional treatments at ICU admission
Steroids	178 (78)	51 (79)	127 (78)	1
Vasoactive drugs	102 (45)	27 (42)	75 (46)	0.686
Invasive mechanical ventilation	210 (93)	59 (91)	151 (93)	0.724
Neuromuscular blocking agents	99 (44)	30 (46)	69 (43)	0.733*
Prone ventilation	119 (52)	32 (49)	87 (54)	0.657*
ECMO	5 (2)	2 (3)	3 (2)	0.544*

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 Figure 2—figure supplement 1). 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 2). 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 Figure 2—figure supplement 2). In-silico predictions suggest that IFIH1 downregulation (such as in patients with the rs1990760 TT variant) decreases the expression of pro-inflammatory pathways (Figure 2—figure supplement 3).

Figure 2 with 3 supplements see all

Download asset Open asset

Differences in gene expression according to *IFIH1* rs1990760 genotypes.

(A) Expression of IFIH1 in peripheral blood in patients with CC/CT (n=31) and TT (n=11) variants (p-value calculated using a Wilcoxon test). (B) Heatmap showing expression of the 160 genes with significant differences between variants. (C) Volcano plot showing the distribution of the magnitude of the differences in gene expression (Log₂ Fold change) and their statistical significance. Inflammation-related genes with differential expression and included in the network shown in panel (D) are labeled. (D) Inflammation-related gene network identified using Ingenuity Pathway Analysis on the RNAseq data. Points represent individual patient data. In boxplots, bold line represents the median, lower and upper hinges correspond to the first and third quartiles (25th and 75th percentiles) and upper and lower whiskers extend from the hinge to the largest or smallest value no further than 1.5 times the interquartile range.

Then we assessed the immunological consequences of these differences in gene expression and cell populations. There were no statistically significant differences in inferred percentages of classic CD14+ monocytes, circulating plasma cells, M2 macrophages, and CD56dim NK cells, between genotypes. Patients with the TT genotype showed an increase in hematopoietic precursors and myeloid dendritic cells (Figure 3).

Figure 3 with 1 supplement see all

Download asset Open asset

Differences in circulating cell populations according to *IFIH1* rs1990760 variants (31 and 11 samples from patients with a CC/CT or TT variant, respectively).

Proportions of each cell line were estimated by deconvolution of RNAseq data. P-values were calculated using a Wilcoxon test and adjusted using the Benjamini-Hochberg method for a false discovery rate of 5%. Points represent individual patient data. In boxplots, bold line represents the median, lower and upper hinges correspond to the first and third quartiles (25th and 75th percentiles) and upper and lower whiskers extend from the hinge to the largest or smallest value no further than 1.5 times the interquartile range.

Figure 3—source data 1 Raw data used in Figure 3, showing different peripheral cell populations obtained after deconvolution of RNAseq data.: https://cdn.elifesciences.org/articles/73012/elife-73012-fig3-data1-v1.csv
Download elife-73012-fig3-data1-v1.csv

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, which were below the limit of detection in a large proportion of patients. Patients with the TT genotype showed lower levels of pro-inflammatory mediators, including IL-6, CXCL10, CXCL16, and CCL7 (Figure 4). There were no significant differences in any of the measured cytokines between patients with CC and CT variants (Figure 4—figure supplement 1).

Figure 4 with 1 supplement see all

Download asset Open asset

Differences in serum immune mediators according to *IFIH1* rs1990760 variants (19 and 8 samples from patients with a CC/CT or TT variant, respectively).

P-values were calculated using a Wilcoxon test and adjusted using the Benjamini-Hochberg method for a false discovery rate of 5%. Points represent individual patient data. In boxplots, bold line represents the median, lower and upper hinges correspond to the first and third quartiles (25th and 75th percentiles) and upper and lower whiskers extend from the hinge to the largest or smallest value no further than 1.5 times the interquartile range.

Figure 4—source data 1 Raw data used in Figure 4, showing serum concentration of inflammatory mediators.: https://cdn.elifesciences.org/articles/73012/elife-73012-fig4-data1-v1.csv
Download elife-73012-fig4-data1-v1.csv

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 (WHO Rapid Evidence Appraisal for COVID-19 Therapies (REACT) Working Group et al., 2020). 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 (Figure 5—figure supplement 1). 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 PBMCs 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.

Figure 5 with 1 supplement see all

Download asset Open asset

Expression of transcription factors *STAT1* (A, p=0.012, no significant differences between variants in post hoc tests), *STAT3* (B, p=0.443) and *FOXO3* (C, p=0.002), and cytokine *IL6* (D, p<0.001) induced by MDA5 activation using high molecular weight poly-I:C (a viral RNA analog) and dexamethasone in peripheral mononuclear blood cells.

Results were fitted using a mixed-effects linear model including cell donor as a fixed effect and experimental group and genotype as random effects. Pairwise p-values lower than 0.05 (using Holm’s correction) are shown. Points represent individual patient data. In boxplots, bold line represents the median, lower and upper hinges correspond to the first and third quartiles (the 25th and 75th percentiles) and upper and lower whiskers extend from the hinge to the largest or smallest value no further than 1.5 times the interquartile range.

Figure 5—source data 1 Raw data used in Figure 5, showing gene expression in an ex-vivo model of MDA5 stimulation.: https://cdn.elifesciences.org/articles/73012/elife-73012-fig5-data1-v1.csv
Download elife-73012-fig5-data1-v1.csv

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 file 1c). ICU survival and hospital mortality were modeled 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 file 1c). Regarding ICU survival (Figure 6A and Supplementary file 1c), 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).

Figure 6 with 1 supplement see all

Download asset Open asset

Outcomes according to *IFIH1* rs1990760 variants and treatment with dexamethasone.

(A) Cumulative incidence of ICU discharge alive and spontaneously breathing for each group. (B) Cumulative incidence of hospital death for each group. Outcomes were modeled as competing events. (C) Viral load at diagnosis for each group (p=0.888). (D) Viral clearance evaluated as half-life of an exponential decay function calculated from the maximal viral load in each patient (p=0.056, no significant differences in post hoc tests). (**E, F**) Serum IL-6 concentrations during the first day (E, p<0.001) and 7 days (F, p=0.014) after ICU admission for each group. Differences among groups in panels (**C–F**) were evaluated using an analysis of the variance. Pairwise p-values lower than 0.05 (using Holm’s correction) are shown. (**G, H**) Changes in *IFIH1* (G) and *FOXO3* (H) expression in day 4 of ICU stay for each group (n=6 and 7 for CC/CT variants without and with dexamethasone, respectively, n = 4 and 4 for TT variant without and with dexamethasone, respectively). (I) Genes with dexamethasone-induced changes in gene expression in ICU day four in each rs1990760 variant. Points represent individual patient data. In boxplots, bold line represents the median, lower and upper hinges correspond to the first and third quartiles (25th and 75th percentiles) and upper and lower whiskers extend from the hinge to the largest or smallest value no further than 1.5 times the interquartile range.

Hospital mortality is shown in Figure 6B and Supplementary file 1c. There were nine 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 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 1 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 1 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 (Supplementary file 3). Overall, changes in gene expression were qualitatively different between variants, with only three genes in common (Figure 6I and Figure 6—figure supplement 1).

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 (RECOVERY Collaborative Group et al., 2021) were examined. The T-allele frequency in populations from a Black/Asian ancestry is 0.13 (data available at https://grasp.nhlbi.nih.gov/Covid19GWASResults.aspx), 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 Figure 7—figure supplement 1).

Figure 7 with 1 supplement see all

Download asset Open asset

In-silico clinical trials.

(A) Mortality curves modeled using data from the RECOVERY clinical trial, for each rs1990760 variant and treatment allocation, assuming a sample size of 500 patients per group. (B) Effect of rs1990760 TT variants (expressed as hazard ratio, HR) in simulated clinical trials including 6000 patients from populations with different minor allele frequencies, assigned to standard care. (C) Effect of dexamethasone therapy according to rs1990760 variants and different minor allele frequencies in simulated clinical trials including 6000 patients. Minor allele frequencies were estimated according to reported race in several subsets of patients from the RECOVERY trial: European population (0.61), overall population (0.53), mechanically ventilated patients (0.47), and African American/Asian population (0.13). An intermediate value of 0.30 was added for illustrative purposes. Points and error bars represent the estimated HR and the corresponding 95% confidence interval.

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 (Bovijn et al., 2020; Gómez, 2021; Zhang et al., 2020). 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 pro-inflammatory mediators, NK cell exhaustion, monocyte dysregulation, and emergency hematopoiesis (Schulte-Schrepping et al., 2020; Wen et al., 2020; Wilk et al., 2020). Some of our findings in patients with the rs1990760 TT variant, including lower levels of circulating pro-inflammatory molecules and a shift toward anti-inflammatory cell populations (hematopoietic precursors, M2 macrophages, or CD56^dim NK cells) have been linked to better outcomes in other observational studies (Maucourant et al., 2020).

Recognition of viral RNA by cytosolic receptors leading to the ultimate induction of expression of types I and II IFNs as part of the activation of innate antiviral signaling cascades, may not only relay on MDA5 receptor but also on other retinoid acid-inducible gene (RIG)-I-like receptors, or on the endosomal Toll-like receptor 3 (TLR3) (Takeuchi and Akira, 2010). Interestingly, no differences were found in serum interferons levels. This may be explained by several factors. Most patients included in our study had serum IFNs levels below the inferior detection threshold. This result is similar to the one reported by Galani et al., 2021, where IFN-lambda and type I IFN production were both diminished and delayed in moderate-to-severe COVID-19 patients followed up during hospitalization. In our study, IFNs levels were only measured on day 1 of ICU admission, precluding any conclusions regarding time-dependent changes on these mediators’ levels during, for instance, subacute or long-term infection.

More importantly, transcriptomic analysis did not reveal significantly different expression levels for IFIH1 gene after correction for multiple comparisons, although single-gene comparison showed a downregulated expression in patients with TT genotype. Although expression levels may be similar, full activation of an antiviral response triggering specific signaling intermediates and transcription factors, may vary between genotypes due to differences in MDA5 molecular architecture and function. The effect of rs1990760 on MDA5 function is related to modifications on the molecular architecture (Berke and Modis, 2012; Wu et al., 2013) and has been reported by several previous studies (Looney et al., 2015; Nejentsev et al., 2009).

Collectively, our findings raise the hypothesis that the TT variant could be related to an attenuated pro-inflammatory response. In line with this observation, a TT variant has been associated with lower pro-inflammatory cytokine levels in patients with lupus (Robinson et al., 2011; Zhang et al., 2018, p. 1) and in experimental viral infections (Domsgen et al., 2016).

Development of tolerance to viral diseases can be considered a major evolutionary adaptative response (Schneider and Ayres, 2008), 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 (Xiao et al., 2020), lack functional IFIH1. It has been suggested that this deficiency reduces the inflammatory response to coronavirus infections (Fischer et al., 2020). IL-6 blockade (REMAP-CAP Investigators et al., 2021) or steroids (WHO Rapid Evidence Appraisal for COVID-19 Therapies (REACT) Working Group et al., 2020) are the only treatments that have improved the outcome of critically ill COVID-19 patients. However, mortality rates are still around 27–32% (REMAP-CAP Investigators et al., 2021; WHO Rapid Evidence Appraisal for COVID-19 Therapies (REACT) Working Group et al., 2020), 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 an HR of 0.64, compared to an 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 (Joseph et al., 2016), 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 toward 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 (Rebendenne et al., 2021). However, genetic ablation of IFIH1 results in increased viral loads and decreased cytokine production (Yin et al., 2021), 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 (Villar et al., 2020), 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 (Domsgen et al., 2016), with no relevant differences in viral clearance. However, total absence of MDA5 results in impaired clearance of West Nile virus (Errett et al., 2013).

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 favorable 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 (Bastard, 2020), 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 Pairo-Castineira et al., 2020; Ellinghaus et al., 2020 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 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.

Genotype	Patients	Standard care	Dexamethasone
	4689	849/3139	401/1550
CC/CT	2954	a/1978	b/976
TT	1735	(849 a)/1161	(401-b)/574

Appendix 1

Hazard ratio estimation and in-silico clinical trials

Estimation of hazard ratios from the RECOVERY trial

Our findings raise the hypothesis that steroid therapy will have a larger effect in populations with a low proportion of individuals with a TT genotype. Populations with black/Asian ancestry have an allelic frequency of the T allele of 0.13 (COVID-19 GWAS Results. https://grasp.nhlbi.nih.gov/Covid19GWASResults.aspx). Therefore, the frequency of a TT genotype is 0.017, implying that the majority of the population has a CC/CT variant. Using data from the RECOVERY trial, we assume that the effects of steroids in this population correspond to the effect on patients with a CC/CT variant. In this population, steroids significantly decreased mortality, with a RR of 0.7 (0.51–0.95). Therefore, we assumed that treatment with steroids in patients with these non-TT variants decreases mortality with a RR of 0.7.

White populations have an allelic frequency of the T allele in rs1990760 of 0.61 (COVID-19 GWAS Results. https://grasp.nhlbi.nih.gov/Covid19GWASResults.aspx). Therefore, the proportion of patients with a TT variant is 0.37. GWAS data released by the UK biobank report a coefficient for hospital death in patients with a T allele (compared to positive survivors) of 0.03 (COVID-19 GWAS Results. https://grasp.nhlbi.nih.gov/downloads/COVID19GWAS/05092021/top_UKBB_death_ALLpstv_050921.txt.gz), so the RR for a TT variant (irrespective of the therapy) is e^2·0.03=1.06. Translating these findings to white patients included in the RECOVERY trial yields the following mortality table (assuming that randomization was independent of the rs1990760 variant):

where a and b are the number of deaths in patients with a CC/CT variant assigned to standard care or dexamethasone, respectively.

Considering that the previously calculated RR of steroids in CC/CT population is 0.7, we can estimate:

\frac{b}{976} = 0.7 \frac{a}{1978}

b = 0.3454 a

And assuming a global RR of a TT genotype of 1.06 (irrespective of therapy), we can estimate:

\frac{\frac{(849 - a) + (401 - b)}{1735}}{\frac{a + b}{2954}} = 1.06

Solving these two equations yields the following values:

a = 573

b = 198

Therefore, distribution of mortality according to therapy and rs1990760 variants was estimated as:

Genotype	Patients	Standard care	Dexamethasone
	4689	849/3139	401/1550
CC/CT	2954	573/1978	198/976
TT	1735	276/1161	203/574

From these values, the RR of a TT genotype compared to CC/CT variants in patients not receiving steroids is:

{R R}_{(T T v s \frac{C C}{C T}) N o s t e r o i d s} = \frac{\frac{276}{1161}}{\frac{573}{1978}} = \frac{0.238}{0.290} = 0.821

And the RR related to steroid therapy compared to standard care in patients with a TT variant is:

{R R}_{(D e x a m e t h a s o n e v s s t d . c a r e) T T v a r i a n t} = \frac{\frac{203}{574}}{\frac{276}{1161}} = \frac{0.354}{0.238} = 1.487

In-silico clinical trials

The previously estimated RRs were used to simulate mortality curves and perform in-silico clinical trials. Mortality was modeled using an asymptotic regression curve with a cumulative distribution function:

F (x) = a - a e^{- c t}

where a is the upper asymptote, c is a curvature parameter, e is the base of natural logarithms and t is time. Assuming that final mortality (this is, the upper asymptote) is 10% higher than the observed mortality at day 28 (as in the survival curves of our cohort), we obtain that the curvature parameter is 0.08564. This parameter is constant for all curves, so the differences among groups depend only on the upper asymptote. This model is appropriate to fit acute conditions in which most of the deaths occur in the short term, but a substantial proportion of patients survive the disease. The code for these simulations can be found at https://github.com/Crit-Lab/IFIH1_simulation (copy archived at swh:1:rev:343c92a4ce58b831df7b289bf40b539db5298d0a, Crit-Lab, 2022).

Different clinical trials including patients randomized into four groups depending on the rs1990760 variant (CC/CT or TT) and prescription of dexamethasone were simulated. Randomization to dexamethasone or placebo was 1:1. Minor allele frequencies from 0.13 to 0.61 (corresponding to populational distributions) and two baseline mortality rates (by adjusting the mortality of patients with a CC/CT variant not treated with steroids to 25.7% and 41.4%, corresponding to mortality rates of the overall RECOVERY trial and those receiving mechanical ventilation respectively) were tested. Each set of conditions was repeated 1000 times and the average RR was calculated. All the simulations were done using the Mediana package for R (available at https://CRAN.R-project.org/package=Mediana).

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 (copy archived at swh:1:rev:343c92a4ce58b831df7b289bf40b539db5298d0a). Source data files for figures 3, 4 and 5 have been provided. Patient-derived data related to this article are not publicly available out of respect for the privacy of the patients involved. 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. All requests should be sent to the corresponding author (gma@crit-lab.org).

The following data sets were generated

(2021) NCBI Gene Expression Omnibus
ID GSE168400. Peripheral blood gene expression according to IFIH1 rs1990760 variants in critically-ill COVID-19 patients.

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE168400
(2021) NCBI Gene Expression Omnibus
ID GSE177025. Effects of steroids on peripheral blood gene expression in COVID-19 patients according to rs1990760 variants.

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE177025

References

1. Barbaro RP
2. MacLaren G
3. Boonstra PS
4. Iwashyna TJ
5. Slutsky AS
6. Fan E
7. Bartlett RH
8. Tonna JE
9. Hyslop R
10. Fanning JJ
11. Rycus PT
12. Hyer SJ
13. Anders MM
14. Agerstrand CL
15. Hryniewicz K
16. Diaz R
17. Lorusso R
18. Combes A
19. Brodie D
20. Extracorporeal Life Support Organization
(2020) Extracorporeal membrane oxygenation support in COVID-19: an international cohort study of the Extracorporeal Life Support Organization registry
Lancet 396:1071–1078.
https://doi.org/10.1016/S0140-6736(20)32008-0
- PubMed
- Google Scholar
1. Bastard P
(2020) Autoantibodies against type I IFNs in patients with life-threatening COVID-19
Science 370:eabd4585.
https://doi.org/10.1126/science.abd4585
- Google Scholar
1. Berke IC
2. Modis Y
(2012) MDA5 cooperatively forms dimers and ATP-sensitive filaments upon binding double-stranded RNA
The EMBO Journal 31:1714–1726.
https://doi.org/10.1038/emboj.2012.19
- PubMed
- Google Scholar
(2020) Genetic variants mimicking therapeutic inhibition of IL-6 receptor signaling and risk of COVID-19
The Lancet. Rheumatology 2:e658–e659.
https://doi.org/10.1016/S2665-9913(20)30345-3
- PubMed
- Google Scholar
Software
1. Crit-Lab
(2022) IFIH1_simulation, version swh:1:rev:343c92a4ce58b831df7b289bf40b539db5298d0a
Software Heritage.

https://archive.softwareheritage.org/swh:1:dir:191b32f9b0a371a6c0610fe7688b6ae38d19544e;origin=https://github.com/Crit-Lab/IFIH1_simulation;visit=swh:1:snp:2946993d2d56beebd48df44b93b4447114a19070;anchor=swh:1:rev:343c92a4ce58b831df7b289bf40b539db5298d0a
(2016) An IFIH1 gene polymorphism associated with risk for autoimmunity regulates canonical antiviral defence pathways in Coxsackievirus infected human pancreatic islets
Scientific Reports 6:39378.
https://doi.org/10.1038/srep39378
- PubMed
- Google Scholar
Software
1. Dowle M
2. Srinivasan A
(2021) data.table: Extension of `data.frame’
Data.Frame.

https://cran.r-project.org
1. Du Y
2. Tu L
3. Zhu P
4. Mu M
5. Wang R
6. Yang P
7. Wang X
8. Hu C
9. Ping R
10. Hu P
11. Li T
12. Cao F
13. Chang C
14. Hu Q
15. Jin Y
16. Xu G
(2020) Clinical Features of 85 Fatal Cases of COVID-19 from Wuhan A Retrospective Observational Study
American Journal of Respiratory and Critical Care Medicine 201:1372–1379.
https://doi.org/10.1164/rccm.202003-0543OC
- PubMed
- Google Scholar
1. Ellinghaus D
2. Degenhardt F
3. Bujanda L
4. Buti M
5. Albillos A
6. Invernizzi P
7. Fernández J
8. Prati D
9. Baselli G
10. Asselta R
11. Grimsrud MM
12. Milani C
13. Aziz F
14. Kässens J
15. May S
16. Wendorff M
17. Wienbrandt L
18. Uellendahl-Werth F
19. Zheng T
20. Yi X
21. de Pablo R
22. Chercoles AG
23. Palom A
24. Garcia-Fernandez AE
25. Rodriguez-Frias F
26. Zanella A
27. Bandera A
28. Protti A
29. Aghemo A
30. Lleo A
31. Biondi A
32. Caballero-Garralda A
33. Gori A
34. Tanck A
35. Carreras Nolla A
36. Latiano A
37. Fracanzani AL
38. Peschuck A
39. Julià A
40. Pesenti A
41. Voza A
42. Jiménez D
43. Mateos B
44. Nafria Jimenez B
45. Quereda C
46. Paccapelo C
47. Gassner C
48. Angelini C
49. Cea C
50. Solier A
51. Pestaña D
52. Muñiz-Diaz E
53. Sandoval E
54. Paraboschi EM
55. Navas E
56. García Sánchez F
57. Ceriotti F
58. Martinelli-Boneschi F
59. Peyvandi F
60. Blasi F
61. Téllez L
62. Blanco-Grau A
63. Hemmrich-Stanisak G
64. Grasselli G
65. Costantino G
66. Cardamone G
67. Foti G
68. Aneli S
69. Kurihara H
70. ElAbd H
71. My I
72. Galván-Femenia I
73. Martín J
74. Erdmann J
75. Ferrusquía-Acosta J
76. Garcia-Etxebarria K
77. Izquierdo-Sanchez L
78. Bettini LR
79. Sumoy L
80. Terranova L
81. Moreira L
82. Santoro L
83. Scudeller L
84. Mesonero F
85. Roade L
86. Rühlemann MC
87. Schaefer M
88. Carrabba M
89. Riveiro-Barciela M
90. Figuera Basso ME
91. Valsecchi MG
92. Hernandez-Tejero M
93. Acosta-Herrera M
94. Baldini M
95. Cazzaniga M
96. Schulzky M
97. Cecconi M
98. Wittig M
99. Ciccarelli M
100. Rodríguez-Gandía M
101. Bocciolone M
102. Miozzo M
103. Montano N
104. Braun N
105. Sacchi N
106. Martínez N
107. Özer O
108. Palmieri O
109. Faverio P
110. Preatoni P
111. Bonfanti P
112. Omodei P
113. Tentorio P
114. Castro P
115. Rodrigues PM
116. Blandino Ortiz A
117. de Cid R
118. Ferrer R
119. Gualtierotti R
120. Nieto R
121. Goerg S
122. Badalamenti S
123. Marsal S
124. Matullo G
125. Pelusi S
126. Juzenas S
127. Aliberti S
128. Monzani V
129. Moreno V
130. Wesse T
131. Lenz TL
132. Pumarola T
133. Rimoldi V
134. Bosari S
135. Albrecht W
136. Peter W
137. D’ M
138. Duga S
139. Banales JM
140. Hov JR
141. Folseraas T
142. Valenti L
143. Franke A
144. Karlsen TH
145. Severe Covid-19 GWAS Group
(2020) Genomewide Association Study of Severe Covid-19 with Respiratory Failure
The New England Journal of Medicine 383:1522–1534.
https://doi.org/10.1056/NEJMoa2020283
- PubMed
- Google Scholar
1. Errett JS
2. Suthar MS
3. McMillan A
4. Diamond MS
5. Gale M
(2013) The essential, nonredundant roles of RIG-I and MDA5 in detecting and controlling West Nile virus infection
Journal of Virology 87:11416–11425.
https://doi.org/10.1128/JVI.01488-13
- PubMed
- Google Scholar
(2020) SARS-CoV-2 analysis on environmental surfaces collected in an intensive care unit: keeping Ernest Shackleton’s spirit
Intensive Care Medicine Experimental 8:68.
https://doi.org/10.1186/s40635-020-00349-5
- PubMed
- Google Scholar
(2020) Pangolins Lack IFIH1/MDA5, a Cytoplasmic RNA Sensor That Initiates Innate Immune Defense Upon Coronavirus Infection
Frontiers in Immunology 11:939.
https://doi.org/10.3389/fimmu.2020.00939
- PubMed
- Google Scholar
1. Galani IE
2. Rovina N
3. Lampropoulou V
4. Triantafyllia V
5. Manioudaki M
6. Pavlos E
7. Koukaki E
8. Fragkou PC
9. Panou V
10. Rapti V
11. Koltsida O
12. Mentis A
13. Koulouris N
14. Tsiodras S
15. Koutsoukou A
16. Andreakos E
(2021) Untuned antiviral immunity in COVID-19 revealed by temporal type I/III interferon patterns and flu comparison
Nature Immunology 22:32–40.
https://doi.org/10.1038/s41590-020-00840-x
- PubMed
- Google Scholar
1. Gómez J
(2021) The Interferon-induced transmembrane protein 3 gene (IFITM3) rs12252 C variant is associated with COVID-19
Cytokine 137:155354.
https://doi.org/10.1016/j.cyto.2020.155354
- Google Scholar
(2018) Human respiratory syncytial virus load normalized by cell quantification as predictor of acute respiratory tract infection
Journal of Medical Virology 90:861–866.
https://doi.org/10.1002/jmv.25020
- PubMed
- Google Scholar
1. Gorman JA
2. Hundhausen C
3. Errett JS
4. Stone AE
5. Allenspach EJ
6. Ge Y
7. Arkatkar T
8. Clough C
9. Dai X
10. Khim S
11. Pestal K
12. Liggitt D
13. Cerosaletti K
14. Stetson DB
15. James RG
16. Oukka M
17. Concannon P
18. Gale M
19. Buckner JH
20. Rawlings DJ
(2017) The A946T variant of the RNA sensor IFIH1 mediates an interferon program that limits viral infection but increases the risk for autoimmunity
Nature Immunology 18:744–752.
https://doi.org/10.1038/ni.3766
- PubMed
- Google Scholar
1. Grasselli G
2. Zangrillo A
3. Zanella A
4. Antonelli M
5. Cabrini L
6. Castelli A
7. Cereda D
8. Coluccello A
9. Foti G
10. Fumagalli R
11. Iotti G
12. Latronico N
13. Lorini L
14. Merler S
15. Natalini G
16. Piatti A
17. Ranieri MV
18. Scandroglio AM
19. Storti E
20. Cecconi M
21. Pesenti A
(2020) Baseline Characteristics and Outcomes of 1591 Patients Infected With SARS-CoV-2 Admitted to ICUs of the Lombardy Region, Italy
JAMA 323:1574–1581.
https://doi.org/10.1001/jama.2020.5394
- Google Scholar
Conference
1. Haynes WA
2. Vallania F
3. Liu C
4. Bongen E
5. Tomczak A
6. Andres-Terrè M
7. Lofgren S
8. Tam A
9. Deisseroth CA
10. Li MD
11. Sweeney TE
12. Khatri P
(2017) Pacific Symposium on Biocomputing
EMPOWERING MULTI-COHORT GENE EXPRESSION ANALYSIS TO INCREASE REPRODUCIBILITY. pp. 144–153.
https://doi.org/10.1142/9789813207813_0015
- Google Scholar
(2008) Simultaneous inference in general parametric models
Biometrical Journal. Biometrische Zeitschrift 50:346–363.
https://doi.org/10.1002/bimj.200810425
- PubMed
- Google Scholar
1. Joseph J
2. Ametepe ES
3. Haribabu N
4. Agbayani G
5. Krishnan L
6. Blais A
7. Sad S
(2016) Inhibition of ROS and upregulation of inflammatory cytokines by FoxO3a promotes survival against Salmonella typhimurium
Nature Communications 7:12748.
https://doi.org/10.1038/ncomms12748
- PubMed
- Google Scholar
1. Kato H
2. Takeuchi O
3. Mikamo-Satoh E
4. Hirai R
5. Kawai T
6. Matsushita K
7. Hiiragi A
8. Dermody TS
9. Fujita T
10. Akira S
(2008) Length-dependent recognition of double-stranded ribonucleic acids by retinoic acid-inducible gene-I and melanoma differentiation-associated gene 5
The Journal of Experimental Medicine 205:1601–1610.
https://doi.org/10.1084/jem.20080091
- PubMed
- Google Scholar
(2015) Effects of type 1 diabetes-associated IFIH1 polymorphisms on MDA5 function and expression
Current Diabetes Reports 15:96.
https://doi.org/10.1007/s11892-015-0656-8
- PubMed
- Google Scholar
1. Love MI
2. Huber W
3. Anders S
(2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2
Genome Biology 15:550.
https://doi.org/10.1186/s13059-014-0550-8
- PubMed
- Google Scholar
1. Maiti AK
(2020) The African-American population with a low allele frequency of SNP rs1990760 (T allele) in IFIH1 predicts less IFN-beta expression and potential vulnerability to COVID-19 infection
Immunogenetics 72:387–391.
https://doi.org/10.1007/s00251-020-01174-6
- PubMed
- Google Scholar
1. Maucourant C
2. Filipovic I
3. Ponzetta A
4. Aleman S
5. Cornillet M
6. Hertwig L
7. Strunz B
8. Lentini A
9. Reinius B
10. Brownlie D
11. Cuapio A
12. Ask EH
13. Hull RM
14. Haroun-Izquierdo A
15. Schaffer M
16. Klingström J
17. Folkesson E
18. Buggert M
19. Sandberg JK
20. Eriksson LI
21. Rooyackers O
22. Ljunggren HG
23. Malmberg KJ
24. Michaëlsson J
25. Marquardt N
26. Hammer Q
27. Strålin K
28. Björkström NK
(2020) Natural killer cell immunotypes related to COVID-19 disease severity
Science Immunology 5:eabd6832.
https://doi.org/10.1126/sciimmunol.abd6832
- Google Scholar
1. Nejentsev S
2. Walker N
3. Riches D
4. Egholm M
5. Todd JA
(2009) Rare variants of IFIH1, a gene implicated in antiviral responses, protect against type 1 diabetes
Science 324:387–389.
https://doi.org/10.1126/science.1167728
- PubMed
- Google Scholar
1. Pairo-Castineira E
2. Clohisey S
3. Klaric L
4. Bretherick AD
5. Rawlik K
6. Pasko D
7. Walker S
8. Parkinson N
9. Fourman MH
10. Russell CD
11. Furniss J
12. Richmond A
13. Gountouna E
14. Wrobel N
15. Harrison D
16. Wang B
17. Wu Y
18. Meynert A
19. Griffiths F
20. Oosthuyzen W
21. Kousathanas A
22. Moutsianas L
23. Yang Z
24. Zhai R
25. Zheng C
26. Grimes G
27. Beale R
28. Millar J
29. Shih B
30. Keating S
31. Zechner M
32. Haley C
33. Porteous DJ
34. Hayward C
35. Yang J
36. Knight J
37. Summers C
38. Shankar-Hari M
39. Klenerman P
40. Turtle L
41. Ho A
42. Moore SC
43. Hinds C
44. Horby P
45. Nichol A
46. Maslove D
47. Ling L
48. McAuley D
49. Montgomery H
50. Walsh T
51. Pereira A
52. Renieri A
(2020) Genetic mechanisms of critical illness in Covid-19
Nature 591:92–98.
https://doi.org/10.1038/s41586-020-03065-y
- Google Scholar
1. Patro R
2. Duggal G
3. Love MI
4. Irizarry RA
5. Kingsford C
(2017) Salmon provides fast and bias-aware quantification of transcript expression
Nature Methods 14:417–419.
https://doi.org/10.1038/nmeth.4197
- PubMed
- Google Scholar
Software
1. R Development Core Team
(2020) R: A language and environment for statistical computing
R Foundation for Statistical Computing, Vienna, Austria.

http://www.R-project.org
(2021) SARS-CoV-2 triggers an MDA-5-dependent interferon response which is unable to control replication in lung epithelial cells
Journal of Virology 20:JVI.02415-20.
https://doi.org/10.1128/JVI.02415-20
- PubMed
- Google Scholar
1. RECOVERY Collaborative Group
2. Horby P
3. Lim WS
4. Emberson JR
5. Mafham M
6. Bell JL
7. Linsell L
8. Staplin N
9. Brightling C
10. Ustianowski A
11. Elmahi E
12. Prudon B
13. Green C
14. Felton T
15. Chadwick D
16. Rege K
17. Fegan C
18. Chappell LC
19. Faust SN
20. Jaki T
21. Jeffery K
22. Montgomery A
23. Rowan K
24. Juszczak E
25. Baillie JK
26. Haynes R
27. Landray MJ
(2021) Dexamethasone in Hospitalized Patients with Covid-19
The New England Journal of Medicine 384:693–704.
https://doi.org/10.1056/NEJMoa2021436
- PubMed
- Google Scholar
1. REMAP-CAP Investigators
2. Gordon AC
3. Mouncey PR
4. Al-Beidh F
5. Rowan KM
6. Nichol AD
7. Arabi YM
8. Annane D
9. Beane A
10. van Bentum-Puijk W
11. Berry LR
12. Bhimani Z
13. Bonten MJM
14. Bradbury CA
15. Brunkhorst FM
16. Buzgau A
17. Cheng AC
18. Detry MA
19. Duffy EJ
20. Estcourt LJ
21. Fitzgerald M
22. Goossens H
23. Haniffa R
24. Higgins AM
25. Hills TE
26. Horvat CM
27. Lamontagne F
28. Lawler PR
29. Leavis HL
30. Linstrum KM
31. Litton E
32. Lorenzi E
33. Marshall JC
34. Mayr FB
35. McAuley DF
36. McGlothlin A
37. McGuinness SP
38. McVerry BJ
39. Montgomery SK
40. Morpeth SC
41. Murthy S
42. Orr K
43. Parke RL
44. Parker JC
45. Patanwala AE
46. Pettilä V
47. Rademaker E
48. Santos MS
49. Saunders CT
50. Seymour CW
51. Shankar-Hari M
52. Sligl WI
53. Turgeon AF
54. Turner AM
55. van de Veerdonk FL
56. Zarychanski R
57. Green C
58. Lewis RJ
59. Angus DC
60. McArthur CJ
61. Berry S
62. Webb SA
63. Derde LPG
(2021) Interleukin-6 Receptor Antagonists in Critically Ill Patients with Covid-19
The New England Journal of Medicine 384:1491–1502.
https://doi.org/10.1056/NEJMoa2100433
- PubMed
- Google Scholar
(2016) Hospital Incidence and Outcomes of the Acute Respiratory Distress Syndrome Using the Kigali Modification of the Berlin Definition
American Journal of Respiratory and Critical Care Medicine 193:52–59.
https://doi.org/10.1164/rccm.201503-0584OC
- PubMed
- Google Scholar
1. Robinson T
2. Kariuki SN
3. Franek BS
4. Kumabe M
5. Kumar AA
6. Badaracco M
7. Mikolaitis RA
8. Guerrero G
9. Utset TO
10. Drevlow BE
11. Zaacks LS
12. Grober JS
13. Cohen LM
14. Kirou KA
15. Crow MK
16. Jolly M
17. Niewold TB
(2011) Autoimmune disease risk variant of IFIH1 is associated with increased sensitivity to IFN-α and serologic autoimmunity in lupus patients
Journal of Immunology 187:1298–1303.
https://doi.org/10.4049/jimmunol.1100857
- PubMed
- Google Scholar
1. Schneider DS
2. Ayres JS
(2008) Two ways to survive infection: what resistance and tolerance can teach us about treating infectious diseases
Nature Reviews. Immunology 8:889–895.
https://doi.org/10.1038/nri2432
- PubMed
- Google Scholar
1. Schulte-Schrepping J
2. Reusch N
3. Paclik D
4. Baßler K
5. Schlickeiser S
6. Zhang B
7. Krämer B
8. Krammer T
9. Brumhard S
10. Bonaguro L
11. De Domenico E
12. Wendisch D
13. Grasshoff M
14. Kapellos TS
15. Beckstette M
16. Pecht T
17. Saglam A
18. Dietrich O
19. Mei HE
20. Schulz AR
21. Conrad C
22. Kunkel D
23. Vafadarnejad E
24. Xu C-J
25. Horne A
26. Herbert M
27. Drews A
28. Thibeault C
29. Pfeiffer M
30. Hippenstiel S
31. Hocke A
32. Müller-Redetzky H
33. Heim K-M
34. Machleidt F
35. Uhrig A
36. Bosquillon de Jarcy L
37. Jürgens L
38. Stegemann M
39. Glösenkamp CR
40. Volk H-D
41. Goffinet C
42. Landthaler M
43. Wyler E
44. Georg P
45. Schneider M
46. Dang-Heine C
47. Neuwinger N
48. Kappert K
49. Tauber R
50. Corman V
51. Raabe J
52. Kaiser KM
53. Vinh MT
54. Rieke G
55. Meisel C
56. Ulas T
57. Becker M
58. Geffers R
59. Witzenrath M
60. Drosten C
61. Suttorp N
62. von Kalle C
63. Kurth F
64. Händler K
65. Schultze JL
66. Aschenbrenner AC
67. Li Y
68. Nattermann J
69. Sawitzki B
70. Saliba A-E
71. Sander LE
72. Deutsche COVID-19 OMICS Initiative (DeCOI)
(2020) Severe COVID-19 Is Marked by a Dysregulated Myeloid Cell Compartment
Cell 182:1419–1440.
https://doi.org/10.1016/j.cell.2020.08.001
- PubMed
- Google Scholar
1. Takeuchi O
2. Akira S
(2010) Pattern recognition receptors and inflammation
Cell 140:805–820.
https://doi.org/10.1016/j.cell.2010.01.022
- PubMed
- Google Scholar
Software
1. Therneau TM
(2020) A Package for Survival Analysis in R
Survival.

https://cran.r-project.org
1. Vallania F
2. Tam A
3. Lofgren S
4. Schaffert S
5. Azad TD
6. Bongen E
7. Haynes W
8. Alsup M
9. Alonso M
10. Davis M
11. Engleman E
12. Khatri P
(2018) Leveraging heterogeneity across multiple datasets increases cell-mixture deconvolution accuracy and reduces biological and technical biases
Nature Communications 9:4735.
https://doi.org/10.1038/s41467-018-07242-6
- PubMed
- Google Scholar
(2020) Dexamethasone treatment for the acute respiratory distress syndrome: a multicentre, randomised controlled trial
The Lancet. Respiratory Medicine 8:267–276.
https://doi.org/10.1016/S2213-2600(19)30417-5
- PubMed
- Google Scholar
(2020) The interferon-induced helicase C domain-containing protein 1 gene variant (rs1990760) as an autoimmune-based pathology susceptibility factor
Immunobiology 225:151864.
https://doi.org/10.1016/j.imbio.2019.10.013
- PubMed
- Google Scholar
1. Wen W
2. Su W
3. Tang H
4. Le W
5. Zhang X
6. Zheng Y
7. Liu X
8. Xie L
9. Li J
10. Ye J
11. Dong L
12. Cui X
13. Miao Y
14. Wang D
15. Dong J
16. Xiao C
17. Chen W
18. Wang H
(2020) Immune cell profiling of COVID-19 patients in the recovery stage by single-cell sequencing
Cell Discovery 6:31.
https://doi.org/10.1038/s41421-020-0168-9
- PubMed
- Google Scholar
1. WHO Rapid Evidence Appraisal for COVID-19 Therapies (REACT) Working Group
2. Sterne JAC
3. Murthy S
4. Diaz JV
5. Slutsky AS
6. Villar J
7. Angus DC
8. Annane D
9. Azevedo LCP
10. Berwanger O
11. Cavalcanti AB
12. Dequin PF
13. Du B
14. Emberson J
15. Fisher D
16. Giraudeau B
17. Gordon AC
18. Granholm A
19. Green C
20. Haynes R
21. Heming N
22. Higgins JPT
23. Horby P
24. Jüni P
25. Landray MJ
26. Le Gouge A
27. Leclerc M
28. Lim WS
29. Machado FR
30. McArthur C
31. Meziani F
32. Møller MH
33. Perner A
34. Petersen MW
35. Savovic J
36. Tomazini B
37. Veiga VC
38. Webb S
39. Marshall JC
(2020) Association Between Administration of Systemic Corticosteroids and Mortality Among Critically Ill Patients With COVID-19: A Meta-analysis
JAMA 324:1330–1341.
https://doi.org/10.1001/jama.2020.17023
- PubMed
- Google Scholar
Book
1. Wickham H
(2016) Ggplot2: Elegant Graphics for Data Analysis
Cham: Springer-Verlag.
https://doi.org/10.1007/978-3-319-24277-4
- Google Scholar
1. Wilk AJ
2. Rustagi A
3. Zhao NQ
4. Roque J
5. Martínez-Colón GJ
6. McKechnie JL
7. Ivison GT
8. Ranganath T
9. Vergara R
10. Hollis T
11. Simpson LJ
12. Grant P
13. Subramanian A
14. Rogers AJ
15. Blish CA
(2020) A single-cell atlas of the peripheral immune response in patients with severe COVID-19
Nature Medicine 26:1070–1076.
https://doi.org/10.1038/s41591-020-0944-y
- PubMed
- Google Scholar
1. Wu B
2. Peisley A
3. Richards C
4. Yao H
5. Zeng X
6. Lin C
7. Chu F
8. Walz T
9. Hur S
(2013) Structural basis for dsRNA recognition, filament formation, and antiviral signal activation by MDA5
Cell 152:276–289.
https://doi.org/10.1016/j.cell.2012.11.048
- PubMed
- Google Scholar
1. Xiao K
2. Zhai J
3. Feng Y
4. Zhou N
5. Zhang X
6. Zou J-J
7. Li N
8. Guo Y
9. Li X
10. Shen X
11. Zhang Z
12. Shu F
13. Huang W
14. Li Y
15. Zhang Z
16. Chen R-A
17. Wu Y-J
18. Peng S-M
19. Huang M
20. Xie W-J
21. Cai Q-H
22. Hou F-H
23. Chen W
24. Xiao L
25. Shen Y
(2020) Isolation of SARS-CoV-2-related coronavirus from Malayan pangolins
Nature 583:286–289.
https://doi.org/10.1038/s41586-020-2313-x
- PubMed
- Google Scholar
1. Yin X
2. Riva L
3. Pu Y
4. Martin-Sancho L
5. Kanamune J
6. Yamamoto Y
7. Sakai K
8. Gotoh S
9. Miorin L
10. De Jesus PD
11. Yang CC
12. Herbert KM
13. Yoh S
14. Hultquist JF
15. García-Sastre A
16. Chanda SK
(2021) MDA5 Governs the Innate Immune Response to SARS-CoV-2 in Lung Epithelial Cells
Cell Reports 34:108628.
https://doi.org/10.1016/j.celrep.2020.108628
- PubMed
- Google Scholar
1. Zhang J
2. Liu X
3. Meng Y
4. Wu H
5. Wu Y
6. Yang B
7. Wang L
(2018) Autoimmune disease associated IFIH1 single nucleotide polymorphism related with IL-18 serum levels in Chinese systemic lupus erythematosus patients
Scientific Reports 8:9442.
https://doi.org/10.1038/s41598-018-27782-7
- PubMed
- Google Scholar
1. Zhang Q
2. Bastard P
3. Liu Z
4. Le Pen J
5. Moncada-Velez M
6. Chen J
7. Ogishi M
8. Sabli IKD
9. Hodeib S
10. Korol C
11. Rosain J
12. Bilguvar K
13. Ye J
14. Bolze A
15. Bigio B
16. Yang R
17. Arias AA
18. Zhou Q
19. Zhang Y
20. Onodi F
21. Korniotis S
22. Karpf L
23. Philippot Q
24. Chbihi M
25. Bonnet-Madin L
26. Dorgham K
27. Smith N
28. Schneider WM
29. Razooky BS
30. Hoffmann HH
31. Michailidis E
32. Moens L
33. Han JE
34. Lorenzo L
35. Bizien L
36. Meade P
37. Neehus AL
38. Ugurbil AC
39. Corneau A
40. Kerner G
41. Zhang P
42. Rapaport F
43. Seeleuthner Y
44. Manry J
45. Masson C
46. Schmitt Y
47. Schlüter A
48. Le Voyer T
49. Khan T
50. Li J
51. Fellay J
52. Roussel L
53. Shahrooei M
54. Alosaimi MF
55. Mansouri D
56. Al-Saud H
57. Al-Mulla F
58. Almourfi F
59. Al-Muhsen SZ
60. Alsohime F
61. Al Turki S
62. Hasanato R
63. van de Beek D
64. Biondi A
65. Bettini LR
66. D’Angio’ M
67. Bonfanti P
68. Imberti L
69. Sottini A
70. Paghera S
71. Quiros-Roldan E
72. Rossi C
73. Oler AJ
74. Tompkins MF
75. Alba C
76. Vandernoot I
77. Goffard JC
78. Smits G
79. Migeotte I
80. Haerynck F
81. Soler-Palacin P
82. Martin-Nalda A
83. Colobran R
84. Morange PE
85. Keles S
86. Çölkesen F
87. Ozcelik T
88. Yasar KK
89. Senoglu S
90. Karabela ŞN
91. Rodríguez-Gallego C
92. Novelli G
93. Hraiech S
94. Tandjaoui-Lambiotte Y
95. Duval X
96. Laouénan C
(2020) Inborn errors of type I IFN immunity in patients with life-threatening COVID-19
Science 370:eabd4570.
https://doi.org/10.1126/science.abd4570
- Google Scholar

Article and author information

Author details

Laura Amado-Rodríguez
1. Unidad de Cuidados Intensivos Cardiológicos. Hospital Universitario Central de Asturias, Oviedo, Spain
2. Centro de Investigación Biomédica en Red (CIBER)-Enfermedades Respiratorias. Instituto de Salud Carlos III, Madrid, Spain
3. Instituto de Investigación Sanitaria del Principado de Asturias, Oviedo, Spain
4. Instituto Universitario de Oncología del Principado de Asturias, Oviedo, Spain
Contribution
Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Supervision, Validation, Visualization, Writing – original draft, Writing – review and editing

Contributed equally with
Estefania Salgado del Riego

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-8793-0213
Estefania Salgado del Riego

Unidad de Cuidados Intensivos Polivalente. Hospital Universitario Central de Asturias, Oviedo, Spain

Contribution
Investigation, Resources, Writing – review and editing

Contributed equally with
Laura Amado-Rodríguez

Competing interests
No competing interests declared
Juan Gomez de Ona
1. Instituto de Investigación Sanitaria del Principado de Asturias, Oviedo, Spain
2. Red de Investigación Renal (REDINREN), Madrid, Spain
3. Servicio de Inmunología. Hospital Universitario Central de Asturias, Oviedo, Spain
Contribution
Data curation, Investigation

Competing interests
No competing interests declared
Inés López Alonso
1. Centro de Investigación Biomédica en Red (CIBER)-Enfermedades Respiratorias. Instituto de Salud Carlos III, Madrid, Spain
2. Instituto de Investigación Sanitaria del Principado de Asturias, Oviedo, Spain
3. Instituto Universitario de Oncología del Principado de Asturias, Oviedo, Spain
Contribution
Conceptualization, Data curation, Formal analysis, Investigation, Supervision, Validation, Writing – review and editing

Competing interests
No competing interests declared
Helena Gil-Pena
1. Instituto de Investigación Sanitaria del Principado de Asturias, Oviedo, Spain
2. Area de Gestión Clínica de Pediatría. Hospital Universitario Central de Asturias, Oviedo, Spain
Contribution
Data curation, Investigation, Supervision, Validation

Competing interests
No competing interests declared
Cecilia López-Martínez
1. Centro de Investigación Biomédica en Red (CIBER)-Enfermedades Respiratorias. Instituto de Salud Carlos III, Madrid, Spain
2. Instituto de Investigación Sanitaria del Principado de Asturias, Oviedo, Spain
3. Instituto Universitario de Oncología del Principado de Asturias, Oviedo, Spain
Contribution
Data curation, Formal analysis, Investigation, Methodology, Software, Visualization, Writing – review and editing

Competing interests
No competing interests declared
Paula Martín-Vicente
1. Centro de Investigación Biomédica en Red (CIBER)-Enfermedades Respiratorias. Instituto de Salud Carlos III, Madrid, Spain
2. Instituto de Investigación Sanitaria del Principado de Asturias, Oviedo, Spain
3. Instituto Universitario de Oncología del Principado de Asturias, Oviedo, Spain
Contribution
Data curation, Formal analysis, Investigation, Writing – review and editing

Competing interests
No competing interests declared
Antonio Lopez-Vazquez
1. Instituto de Investigación Sanitaria del Principado de Asturias, Oviedo, Spain
2. Red de Investigación Renal (REDINREN), Madrid, Spain
3. Servicio de Inmunología. Hospital Universitario Central de Asturias, Oviedo, Spain
Contribution
Investigation, Writing – review and editing

Competing interests
No competing interests declared
Adrian Gonzalez Lopez
1. Centro de Investigación Biomédica en Red (CIBER)-Enfermedades Respiratorias. Instituto de Salud Carlos III, Madrid, Spain
2. Department of Anesthesiology and Operative Intensive Care Medicine CCM/CVK, Charité - Universitätsmedizin Berlin, corporate member of Freie Universität Berlin, Humboldt Universität zu Berlin, and Berlin Institute of Health, Berlin, Germany
Contribution
Data curation, Formal analysis, Investigation, Software, Validation, Visualization, Writing – review and editing

Competing interests
No competing interests declared
Elias Cuesta-Llavona
1. Instituto de Investigación Sanitaria del Principado de Asturias, Oviedo, Spain
2. Servicio de Genética Molecular. Hospital Universitario Central de Asturias, Oviedo, Spain
Contribution
Data curation, Investigation, Writing – review and editing

Competing interests
No competing interests declared
Raquel Rodriguez-Garcia
1. Unidad de Cuidados Intensivos Cardiológicos. Hospital Universitario Central de Asturias, Oviedo, Spain
2. Centro de Investigación Biomédica en Red (CIBER)-Enfermedades Respiratorias. Instituto de Salud Carlos III, Madrid, Spain
3. Instituto de Investigación Sanitaria del Principado de Asturias, Oviedo, Spain
Contribution
Investigation, Writing – review and editing

Competing interests
No competing interests declared
Jose Antonio Boga
1. Instituto de Investigación Sanitaria del Principado de Asturias, Oviedo, Spain
2. Servicio de Microbiología. Hospital Universitario Central de Asturias, Oviedo, Spain
Contribution
Investigation, Methodology, Writing – review and editing

Competing interests
No competing interests declared
Marta Elena alvarez-Arguelles
1. Instituto de Investigación Sanitaria del Principado de Asturias, Oviedo, Spain
2. Servicio de Microbiología. Hospital Universitario Central de Asturias, Oviedo, Spain
Contribution
Investigation, Writing – review and editing

Competing interests
No competing interests declared
Juan Mayordomo-Colunga
1. Centro de Investigación Biomédica en Red (CIBER)-Enfermedades Respiratorias. Instituto de Salud Carlos III, Madrid, Spain
2. Instituto de Investigación Sanitaria del Principado de Asturias, Oviedo, Spain
3. Area de Gestión Clínica de Pediatría. Hospital Universitario Central de Asturias, Oviedo, Spain
4. Servicio de Genética Molecular. Hospital Universitario Central de Asturias, Oviedo, Spain
Contribution
Investigation, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-0997-4410
Jose Ramon Vidal-Castineira

Instituto de Investigación Sanitaria del Principado de Asturias, Oviedo, Spain

Contribution
Investigation, Writing – review and editing

Competing interests
No competing interests declared
Irene Crespo
1. Centro de Investigación Biomédica en Red (CIBER)-Enfermedades Respiratorias. Instituto de Salud Carlos III, Madrid, Spain
2. Instituto de Investigación Sanitaria del Principado de Asturias, Oviedo, Spain
3. Universidad de Oviedo, Oviedo, Spain
Contribution
Investigation, Validation, Writing – review and editing

Competing interests
No competing interests declared
Margarita Fernandez

Instituto Universitario de Oncología del Principado de Asturias, Oviedo, Spain

Contribution
Investigation, Writing – review and editing

Competing interests
No competing interests declared
Loreto Criado

Universidad de Oviedo, Oviedo, Spain

Contribution
Investigation, Writing – review and editing

Competing interests
No competing interests declared
Victoria Salvadores

Universidad de Oviedo, Oviedo, Spain

Contribution
Investigation, Writing – review and editing

Competing interests
No competing interests declared
Francisco Jose Jimeno-Demuth

Servicio de Informática. Hospital Universitario Central de Asturias, Oviedo, Spain

Contribution
Formal analysis, Investigation, Software, Writing – review and editing

Competing interests
No competing interests declared
Lluis Blanch
1. Centro de Investigación Biomédica en Red (CIBER)-Enfermedades Respiratorias. Instituto de Salud Carlos III, Madrid, Spain
2. Critical Care Center, Hospital Universitari Parc Taulí, Institut D'Investigació I Innovació Parc Taulí I3PT, Universitat Autònoma de Barcelona, Sabadell, Spain
Contribution
Methodology, Supervision, Validation, Writing – review and editing

Competing interests
No competing interests declared
Belen Prieto
1. Universidad de Oviedo, Oviedo, Spain
2. Servicio de Bioquímica Clínica. Hospital Universitario Central de Asturias, Oviedo, Spain
Contribution
Investigation, Writing – review and editing

Competing interests
No competing interests declared
Alejandra Fernandez-Fernandez

Servicio de Bioquímica Clínica. Hospital Universitario Central de Asturias, Oviedo, Spain

Contribution
Investigation, Writing – review and editing

Competing interests
No competing interests declared
Carlos Lopez-Larrea
1. Instituto de Investigación Sanitaria del Principado de Asturias, Oviedo, Spain
2. Red de Investigación Renal (REDINREN), Madrid, Spain
3. Area de Gestión Clínica de Pediatría. Hospital Universitario Central de Asturias, Oviedo, Spain
4. Universidad de Oviedo, Oviedo, Spain
Contribution
Conceptualization, Funding acquisition, Investigation, Methodology, Supervision, Validation, Writing – review and editing

Competing interests
No competing interests declared
Eliecer Coto
1. Instituto de Investigación Sanitaria del Principado de Asturias, Oviedo, Spain
2. Red de Investigación Renal (REDINREN), Madrid, Spain
3. Servicio de Genética Molecular. Hospital Universitario Central de Asturias, Oviedo, Spain
4. Universidad de Oviedo, Oviedo, Spain
Contribution
Conceptualization, Data curation, Funding acquisition, Investigation, Methodology, Supervision, Validation, Writing – review and editing

Competing interests
No competing interests declared
Guillermo M Albaiceta
1. Unidad de Cuidados Intensivos Cardiológicos. Hospital Universitario Central de Asturias, Oviedo, Spain
2. Centro de Investigación Biomédica en Red (CIBER)-Enfermedades Respiratorias. Instituto de Salud Carlos III, Madrid, Spain
3. Instituto de Investigación Sanitaria del Principado de Asturias, Oviedo, Spain
4. Instituto Universitario de Oncología del Principado de Asturias, Oviedo, Spain
5. Universidad de Oviedo, Oviedo, Spain
Contribution
Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review and editing

For correspondence
gma@crit-lab.org

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-9276-3253

Funding

Instituto de Salud Carlos III (PI19/00184)

Carlos Lopez-Larrea

Instituto de Salud Carlos III (PI20/01360)

Guillermo M Albaiceta

Centro de investigación biomédica en red (CIBER)-Enfermedades respiratorias (CB17/06/00021)

Guillermo M Albaiceta

Fundació la Marató de TV3 (grant 413/C/2021)

Laura Amado-Rodríguez

Instituto de Salud Carlos III (CM20/00083)

Raquel Rodriguez-Garcia

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

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, the authors 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.

Ethics

Human subjects: 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.

Version history

Preprint posted: July 7, 2021 (view preprint)
Received: August 12, 2021
Accepted: December 6, 2021
Version of Record published: January 21, 2022 (version 1)

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

1,237

views
173

downloads
16

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Laura Amado-Rodríguez
Estefania Salgado del Riego
Juan Gomez de Ona
Inés López Alonso
Helena Gil-Pena
Cecilia López-Martínez
Paula Martín-Vicente
Antonio Lopez-Vazquez
Adrian Gonzalez Lopez
Elias Cuesta-Llavona
Raquel Rodriguez-Garcia
Jose Antonio Boga
Marta Elena alvarez-Arguelles
Juan Mayordomo-Colunga
Jose Ramon Vidal-Castineira
Irene Crespo
Margarita Fernandez
Loreto Criado
Victoria Salvadores
Francisco Jose Jimeno-Demuth
Lluis Blanch
Belen Prieto
Alejandra Fernandez-Fernandez
Carlos Lopez-Larrea
Eliecer Coto
Guillermo M Albaiceta

(2022)

Effects of IFIH1 rs1990760 variants on systemic inflammation and outcome in critically ill COVID-19 patients in an observational translational study

eLife 11:e73012.

https://doi.org/10.7554/eLife.73012

Categories and tags

Research organism

Human

Share this article

Cite this article

Patient study flow chart.

Clinical characteristics of the study cohort.

Differences in gene expression according to IFIH1 rs1990760 genotypes.

Differences in circulating cell populations according to IFIH1 rs1990760 variants (31 and 11 samples from patients with a CC/CT or TT variant, respectively).

Figure 3—source data 1

Differences in serum immune mediators according to IFIH1 rs1990760 variants (19 and 8 samples from patients with a CC/CT or TT variant, respectively).

Figure 4—source data 1

Figure 5—source data 1

Outcomes according to IFIH1 rs1990760 variants and treatment with dexamethasone.

In-silico clinical trials.

Author details

Laura Amado-Rodríguez

Contribution

Contributed equally with

Competing interests

Estefania Salgado del Riego

Contribution

Contributed equally with

Competing interests

Juan Gomez de Ona

Contribution

Competing interests

Inés López Alonso

Contribution

Competing interests

Helena Gil-Pena

Contribution

Competing interests

Cecilia López-Martínez

Contribution

Competing interests

Paula Martín-Vicente

Contribution

Competing interests

Antonio Lopez-Vazquez

Contribution

Competing interests

Adrian Gonzalez Lopez

Contribution

Competing interests

Elias Cuesta-Llavona

Contribution

Competing interests

Raquel Rodriguez-Garcia

Contribution

Competing interests

Jose Antonio Boga

Contribution

Competing interests

Marta Elena alvarez-Arguelles

Contribution

Competing interests

Juan Mayordomo-Colunga

Contribution

Competing interests

Jose Ramon Vidal-Castineira

Contribution

Competing interests

Irene Crespo

Contribution

Competing interests

Margarita Fernandez

Contribution

Competing interests

Loreto Criado

Contribution

Competing interests

Victoria Salvadores

Contribution

Competing interests

Francisco Jose Jimeno-Demuth

Contribution

Competing interests

Lluis Blanch

Contribution

Competing interests

Belen Prieto

Contribution