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

  1. Laura Amado-Rodríguez
  2. Estefania Salgado del Riego
  3. Juan Gomez de Ona
  4. Inés López Alonso
  5. Helena Gil-Pena
  6. Cecilia López-Martínez
  7. Paula Martín-Vicente
  8. Antonio Lopez-Vazquez
  9. Adrian Gonzalez Lopez
  10. Elias Cuesta-Llavona
  11. Raquel Rodriguez-Garcia
  12. Jose Antonio Boga
  13. Marta Elena alvarez-Arguelles
  14. Juan Mayordomo-Colunga
  15. Jose Ramon Vidal-Castineira
  16. Irene Crespo
  17. Margarita Fernandez
  18. Loreto Criado
  19. Victoria Salvadores
  20. Francisco Jose Jimeno-Demuth
  21. Lluis Blanch
  22. Belen Prieto
  23. Alejandra Fernandez-Fernandez
  24. Carlos Lopez-Larrea
  25. Eliecer Coto
  26. Guillermo M Albaiceta  Is a corresponding author
  1. Unidad de Cuidados Intensivos Cardiológicos. Hospital Universitario Central de Asturias, Spain
  2. Centro de Investigación Biomédica en Red (CIBER)-Enfermedades Respiratorias. Instituto de Salud Carlos III, Spain
  3. Instituto de Investigación Sanitaria del Principado de Asturias, Spain
  4. Instituto Universitario de Oncología del Principado de Asturias, Spain
  5. Unidad de Cuidados Intensivos Polivalente. Hospital Universitario Central de Asturias, Spain
  6. Red de Investigación Renal (REDINREN), Spain
  7. Servicio de Inmunología. Hospital Universitario Central de Asturias, Spain
  8. Area de Gestión Clínica de Pediatría. Hospital Universitario Central de Asturias, Spain
  9. 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
  10. Servicio de Genética Molecular. Hospital Universitario Central de Asturias, Spain
  11. Servicio de Microbiología. Hospital Universitario Central de Asturias, Spain
  12. Universidad de Oviedo, Spain
  13. Servicio de Informática. Hospital Universitario Central de Asturias, Spain
  14. Critical Care Center, Hospital Universitari Parc Taulí, Institut D'Investigació I Innovació Parc Taulí I3PT, Universitat Autònoma de Barcelona, Spain
  15. Servicio de Bioquímica Clínica. Hospital Universitario Central de Asturias, Spain
7 figures, 1 table and 5 additional files

Figures

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 …

Figure 2 with 3 supplements
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 …

Figure 2—figure supplement 1
IFIH1 expression per rs1990760 variants.
Figure 2—figure supplement 2
Gene networks involved in the regulation of the inflammatory response identified by Ingenuity Pathway Analysis among the genes with differential expression in patients with rs1990760 CC/CT and TT variants.
Figure 2—figure supplement 3
In-silico predictions of IFIH1 downregulation over a gene network involved in the regulation of the inflammatory response.
Figure 3 with 1 supplement
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 …

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
Figure 3—figure supplement 1
Correlation between measured and estimated (from deconvolution analysis) lymphocyte percentages in peripheral blood.
Figure 4 with 1 supplement
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 …

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
Figure 4—figure supplement 1
Serum inflammatory mediators in patients with a CC or CT variant in rs1990760 polymorphisms.

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 …

Figure 5 with 1 supplement
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 …

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
Figure 5—figure supplement 1
Predicted effects of IFIH1 downregulation in absence (A) or in presence (B) of dexamethasone.

A number of molecules downregulated in response to lower IFIH1 expression are upregulated when dexamethasone is added (arrowheads).

Figure 6 with 1 supplement
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) …

Figure 6—figure supplement 1
Heatmaps with the differentially expressed genes in response to steroids in patients of each rs1990760 variant.

Transcriptomic signatures in each variant only share three genes.

Figure 7 with 1 supplement
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 …

Figure 7—figure supplement 1
Simulated survival curves corresponding to in-silico clinical trials testing dexamethasone in 6000 COVID-19 patients from populations with different allelic frequencies.

T allele frequencies correspond to an African/American/Asian population (0.13), mechanically ventilated sample of the RECOVERY trial (0.47), overall sample of the RECOVERY trial (0.53), and a White …

Tables

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 …

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
 Male174 (77)47 (72)127 (78)0.327
 Female53 (23)18 (28)35 (22)
Race
 Black3 (1.5)1 (2)2 (1.5)0.527
 White207 (91)62 (95)146 (90)
 Latino15 (7)2 (3)13 (8)
 Asian1 (0.5)01 (0.5)
Body mass index (kg/m2)30 (27–33)30 (26–33)30 (27–32)0.592
Days since symptom onset8 (6–11)8 (6–10)8 (6–11)0.873
Days from hospital admission to ICU admission2 (1–4)2 (1–3)2 (1–4)0.841
APACHE-II score15 (12–18)15.5 (13–18)14 (11—18)0.177
Comorbidities
Arterial hypertension132 (58)38 (59)94 (58)0.932
Diabetes51 (23)14 (22)37 (23)1
Chronic kidney disease15 (7)5 (8)10 (6)0.743
COPD17 (8)7 (11)10 (6)0.339
Cirrhosis2 (1)1 (2)1 (0.5)1
Neoplasms
 No215 (95)62 (95)3 (5)0154 (95)0.669
 Active10 (4)3 (5)7 (3.5)
 Past2 (1)02 (1.5)
Immunosuppresive drugs (incl. steroids)5 (2)3 (5)2 (1.5)0.142
Lung function at ICU admission
Ventilation at admission
 Spontaneous /NIV38 (17)11 (17)27 (17)0.278
 Controlled invasive ventilation187 (82.5)53 (82)135 (83)
 Pressure support ventilation1 (0.5)1 (2)0
FiO20.5 (0.4–0.6)0.5 (0.4–0.7)0.5 (0.4–0.6)0.104
PaO2/FiO2 (mmHg)204 (155–267)185 (142–226)218 (162–282)0.007
PaCO2 (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 pH7.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 (cmH2O)24 (21–28)24 (22–29)24 (21–27)0.221
PEEP (cmH2O)12 (10–14)12 (10–14)12 (10–14)0.397
Driving pressure (cmH2O)12 (10–14)12 (10–14)12 (10–14)0.496
Respiratory system compliance (ml/cmH2O)38 (32–50)36 (31–44)42 (32–50)0.11
Laboratory results
Leukocytes (×103/μl)8.34 (6.05–11.38)8.77 (5.96–11.52)8.29 (6.08–11.18)0.773
Lymphocytes (×103/μ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
Steroids178 (78)51 (79)127 (78)1
Vasoactive drugs102 (45)27 (42)75 (46)0.686
Invasive mechanical ventilation210 (93)59 (91)151 (93)0.724
Neuromuscular blocking agents99 (44)30 (46)69 (43)0.733*
Prone ventilation119 (52)32 (49)87 (54)0.657*
ECMO5 (2)2 (3)3 (2)0.544*

Additional files

Supplementary file 1

Supplementary Tables.

1a. Primers and probes used to detect SARS-CoV-2. 1b: Primers used for qPCR of human genes. 1c: Clinical characteristics of patients according to IFIH1 rs1990760 genotypes and dexamethasone therapy. 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.

https://cdn.elifesciences.org/articles/73012/elife-73012-supp1-v1.docx
Supplementary file 2

Results of differential expression analysis between rs1990760 variants.

https://cdn.elifesciences.org/articles/73012/elife-73012-supp2-v1.csv
Supplementary file 3

Changes in gene expression after steroid treatment in rs1990760 variants.

https://cdn.elifesciences.org/articles/73012/elife-73012-supp3-v1.csv
Transparent reporting form
https://cdn.elifesciences.org/articles/73012/elife-73012-transrepform1-v1.docx
Reporting standard 1

STROBE check list for cohort studies.

https://cdn.elifesciences.org/articles/73012/elife-73012-repstand1-v1.doc

Download links