key: cord-0877829-a9fiwwb0
authors: Girona-Alarcon, Mònica; Argüello, Guillermo; Esteve-Sole, Ana; Bobillo-Perez, Sara; Paolo-Burgos, Xavier; Bonet-Carne, Elisenda; Mensa-Vilaró, Anna; Codina, Anna; Hernández-Garcia, María; Jou, Cristina; Alsina, Laia; Jordan, Iolanda
title: Low levels of CIITA and high levels of SOCS1 predict COVID-19 disease severity in children and adults
date: 2021-12-09
journal: iScience
DOI: 10.1016/j.isci.2021.103595
sha: f462e4ee6bc68bed44f23c1ea354b19977bcfe2f
doc_id: 877829
cord_uid: a9fiwwb0

It is unclear why COVID-19 ranges from asymptomatic to severe. When SARS-CoV-2 is detected, interferon (IFN) response is activated. When it is insufficient or delayed, it might lead to overproduction of cytokines and severe COVID-19. The aim was to compare cytokine and IFN patterns in children and adults with differing severity with SARS-CoV-2. It was a prospective, observational study, including 84 patients. Patients with moderate/severe disease had higher cytokines’ values than patients with mild disease (p<0.001). Two IFN genes were selected to build a decision tree for severity classification: SOCS1 (representative of the rest of the IFN genes) and CIITA (inverse correlation). Low values of CIITA and high values of SOCS1 indicated severe disease. This method correctly classified 33/38(86.8%) of children and 27/34 (79.4%) of adults. To conclude, patients with severe disease had an elevated cytokine pattern, which correlated with the IFN response, with low CIITA and high SOCS1 values.

It is unclear why COVID-19 ranges from asymptomatic to severe. When SARS-CoV-2 is detected, interferon (IFN) response is activated. When it is insufficient or delayed, it might lead to overproduction of cytokines and severe COVID-19. The aim was to compare cytokine and IFN patterns in children and adults with differing severity with SARS-CoV-2.

It was a prospective, observational study, including 84 patients. Patients with moderate/severe disease had higher cytokines' values than patients with mild disease (p<0.001).

Two IFN genes were selected to build a decision tree for severity classification: SOCS1

(representative of the rest of the IFN genes) and CIITA (inverse correlation). Low values of CIITA and high values of SOCS1 indicated severe disease. This method correctly classified

In December 2019, a severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was isolated. It has caused a pandemic that has been declared a significant threat to international health by the World Health Organisation (WHO). To date, 193,657,725 cases have been reported, with 4,154,660 deaths(WHO, 2021) .

In adults, the spectrum of clinical manifestations ranges from asymptomatic, or mild symptoms, to severe pneumonia with acute respiratory syndrome (ARDS) and death.

Concretely, the mortality rate in hospitalised patients is around 20-25% (Richardson et al., 2020) . Secondarily to the severe outcomes of the disease, many researchers are struggling to detect risk factors or biomarkers to detect vulnerable patients in advance. On the one hand, some risk factors such as being male, older age, hypertension, obesity, and diabetes have been reported to be associated with more severe disease (Gao et al., 2021; Richardson et al., 2020) .

For this reason, some risk scores that include these clinical factors have been developed, in order to help to stratify the patients as soon as possible (Bello-Chavolla et al., 2020) . On the other hand, immune factors have also been described as involved with a critical disease course, such as the ones related to the regulation of interferon (IFN). For instance, the presence of autoantibodies against IFN entails blocking IFN activity and, therefore, its capacity to control SARS-CoV-2 infection . Furthermore, genetic inborn errors of IFN-related genes might be also detrimental for the adequate COVID-19 disease control .

In children, the disease seems to be milder, with symptoms similar to other viral respiratory infections, with a low mortality rate of around 0-0.2% (Ludvigsson, 2020; Sankar et al., 2020) .

There are some reasons that have been proposed to explain the milder disease in children.

First, the higher levels of cross-neutralising antibodies, lower levels of ACE-2 receptors in nasal J o u r n a l P r e -p r o o f epithelium, and immature B and T cells, might be responsible for a less specific but more immediate response. In addition, greater regulatory T-cell responses, and lower IL-6 and TNFalpha production may limit the inflammatory response (Wong et al., 2020) .

However, in April 2020, an unexpected severe manifestation of COVID-19, named multisystem inflammatory syndrome (MIS-C), was reported in children. The diagnostic criteria are fever and elevated inflammatory biomarkers, evidence of COVID-19, as well as two of the following: rash, hypotension or shock, myocardial dysfunction, coagulopathy, and acute gastrointestinal problems(World Health Organization, 2020). It remains unclear why some children evolve to MIS-C while others manifest only a mild form of the disease (Ahmed et al., 2020) . MIS-C is the main reason for requiring PICU in children with COVID-19 (Tagarro et al., 2021) , and understanding the immune dysregulation underlying this disease will be fundamental to find other approaches for it.

Both in adults and children, long COVID is a concerning consequence still not well understood.

It is defined by the persistence of symptoms beyond 4 weeks from the onset of the disease, and some mechanisms have been purposed, such as virus-specific pathophysiological changes, immunological alterations and inflammatory damage, but it is still unclear (Asadi-Pooya et al., 2021; Naeije and Caravita, 2021) .

Researchers are struggling to detect clinical or analytical biomarkers that could be of use in determining the risk of an individual patient to evolve to a more severe disease course, including MIS-C. In line with this, some genes have been reported to be linked to a worse prognosis Zhang et al., 2020) , such as TLR7, whose loss-of-function variants have been associated with impaired IFN response ( Van der Made et al., 2020) , and other genes related to the binding of the ACE-2 cell surface receptor and also to entry at the initial stages, which could determine susceptibility to SARS-CoV-2 (Anastassopoulou et al., J o u r n a l P r e -p r o o f 2020). It thus seems that IFN is pivotal in controlling SARS-CoV-2 infection and avoiding a worse clinical course (Hadjadj et al., 2020) .

Indeed, when a viral infection is detected by the innate immune sensors, IFN-I response is activated, and IFN-I response genes (IRG) are secreted by infected cells. They have three functions: i) to limit the spread of the infectious agents, ii) to modulate innate immune responses, and iii) to activate the adaptive immune system (Ivashkiv and Donlin, 2014) . IFN-I binds to the IFN receptor (IFNAR), and it activates an antiviral program of interferonstimulated genes (ISGs) that can interfere with viral replication (Park and Iwasaki, 2020) .

SARS-CoV-2 leads to a lower antiviral transcriptional response, with low IFN-I and IFN-III levels and elevated cytokine expression. Therefore, it seems that patients with severe COVID-19 might have a low or delayed IFN response, which could lead to high cytokine release, ending up in a hyperinflammatory state (Park and Iwasaki, 2020) . That is to say, the infection has been proven to cause an immune dysregulation with overproduction of pro-inflammatory cytokines (IL-1, IL-6, IL-2R, IL-10, and TNF-alpha) in severe cases Ong et al., 2020; Pedersen and Ho, 2020) . What's more, in children with MIS-C, which is the severe manifestation of COVID-19, the cytokine response seems to be greater than in adults (Esteve-Sole et al., 2021) .

We hypothesised that the IFN response and the cytokines released during acute SARS-CoV-2 infection differed depending on the severity of the disease. Moreover, they could be different in children and in adults. Therefore, we hypothesized that we would be able to stratify the risk of patients using these analytical values.

The main aim of the study, then, was to analyse the IFN response and pro-inflammatory cytokines in children and in adults with acute SARS-CoV-2. A secondary objective was to compare the immune response depending on the severity of COVID-19 disease using a machine learning model. Table S3 ).

In adults, patients with severe disease showed higher values for all the studied cytokines than did patients with asymptomatic/mild disease ( Figure 2 , Table S3 ).

Regarding the differences between the cytokine profile in adult severe and children moderate/severe patients, the latter with MIS-C had higher values of all the cytokines, except for IL-6 and MCP-3, which showed no significant differences between the two populations ( Figure 2 , Table S3 ).

The combination of IL-6 and lymphocyte counts correlates with disease severity in children and adults IL-6 is the only cytokine that we are currently able to measure in a general routine blood test;

As it behaves similarly to other cytokines, it was chosen to build a severity decision tree. The severity decision tree included IL-6 circulating levels and lymphocyte counts. Severe patients (adults and children) had low values of lymphocytes at admission and also high IL-6 values. In contrast, mild patients (adults and children) had higher lymphocyte count and lower IL-6 values ( Figure 3 ).

The IFN transcriptional profile was analysed for 72 patients. There were positive correlations for almost all the IFN genes, indicating that when any of the IFN genes was high the others also tended to be high. Only CIITA showed a negative correlation with all the other variables ( Figure   1 ).

J o u r n a l P r e -p r o o f No significant differences were found in child IFN transcriptional profile between severe and moderate patients; thus they were grouped into a moderate/severe cohort (see Table S4 ). There were differences in the IFN response between mild/asymptomatic vs. moderate/severe patients for the following genes: GBP1, HERC6, IFIT3, OASL, SOCS1, and IFNA2, which were higher in moderate/severe than in mild/asymptomatic patients, and HERC6, OAS2, SIGLEC1, USP18, and CIITA, which had lower levels in moderate/severe patients ( Figure 4 , Table S5 ).

After describing the values for the z-score for each IFN gene, a variable selection method was applied, and CIITA and SOCS1 were found to be the most significant genes to build the decision tree (the correlation coefficient between the two was -0.33), which is represented in Figure   5A . With CITTA ≤-1.205 and SOCS1>0.45 the patient was classified as moderate/severe; otherwise, the patient was classified as mild/asymptomatic. The decision tree correctly classified 33/38 observations (AUC 80%).

Differences were found in the following genes: CXCL10, IFI27, IFI44L, SOCS1, CIITA, and IFNA2.

All of them presented higher values in severe patients that in mild ones except for CIITA, which was lower in severe patients than in mild/asymptomatic patients. (Figure 4 , Table S5 ).

A decision tree with CIITA and SOCS1 was also constructed for the adult population ( Figure 5B ).

With CITTA ≤-2.441 and SOCS1>0.491 the patient was classified as moderate/severe; otherwise, the patient was classified as mild/asymptomatic. This model, built for the adult population, correctly classified 27/34 observations (AUC 70%).

Regarding the significant differences between adults and children, moderate/severe children presented lower values of IFI27, LY6E, OAS1, OASL, SIGLEC1, and SPATS2L than did severe adults ( Figure 4 , Table S5 ). Figure 4 shows the distributions of both variables used to train the models, CIITA and SOCS1, for the cohorts in adults and children. No significant differences were found in CIITA and SOCS1 values when comparing moderate/severe adults vs. severe children (p=0. 236 and p=0.191, respectively) . However, children with mild disease had higher values for CIITA than did adults with mild disease (-0.3 vs. -1.4, p=0.009) . No differences were found in SOCS1 values for mild/asymptomatic patients when comparing adults and children (p=0.202).

Finally, the correlation between IFN genes and cytokines was analysed ( 

This is the first study to analyse the IFN and circulating cytokine profiles in children and adults with COVID-19. The importance of the role of interferon has been pointed up in other studies Hadjadj et al., 2020; Zhang et al., 2020) , and, related to its function, important data have come to light in this study. The main finding of the study is that the IFN response is different depending on the severity of the disease. Moreover, some circulating cytokine levels are tightly correlated with IFN gene expression, which could be useful as per future biomarkers and therapeutic targets.

First, from the general blood results, inflammatory biomarkers such as D-dimer, ferritin, LDH, and CRP were higher in severe cases than in milder disease, in line with findings in previous studies . Moreover, T lymphocytes seem to be generally decreased in COVID-J o u r n a l P r e -p r o o f 19, but they were significantly lower in severe cases, which suggests a diminished and delayed adaptive immune response to the virus (Bolouri et al., 2021; Chen et al., 2020) . Then, using the inflammatory biomarkers from the general blood test might be helpful to detect risky patients.

For instance, Gerotziafas et al. have validated a risk score, which includes biomarkers (Ddimer, CRP and lymphocytes count, among others) and clinical characteristics (obesity, male gender) to stratify patients into risk cohorts for severe disease (Gerotziafas et al., 2020) .

However, this study does not include paediatric population, nor immune factors. Including cytokines and IFN response as biomarkers would probably increase the sensitivity and specificity for detecting those vulnerable patients.

In examining the cytokine pattern in the different cohorts, we saw that patients with severe disease (adults and children) had an increased inflammatory cytokine profile with respect to patients with mild disease, which is consistent with other studies(Esteve-Sole et al., 2021) .

Interestingly, even if the described hyperinflammation in MIS-C differs from the severe acute COVID-19 inflammation in adults (Consiglio et al., 2020) , there was an increased inflammatory cytokine profile in both populations, although in children some inflammatory cytokines were higher than in adults. Higher values for IL-10 and TNF-alpha in children with MIS-C, in comparison to severe COVID-19 in adults, have previously been reported (Diorio et al., 2020) .

We have confirmed those results in our sample, since both IL-10 and TNF-alpha were higher in children with severe disease than in adults.

After correlating the general blood test results with the circulating cytokine profile, interesting correlations were found. IL-6, which has been described as elevated in severe disease (Galván-Román et al., 2021) and which is available in the regular blood test, was chosen to build a decision tree together with the lymphocyte count at admission. This correctly classified most patients; thus one of the main proposals emerging from the present study is to include them both in the regular screening of patients with COVID-19, since it allows physicians to discriminate severe from mild patients. Moreover, in cases with elevated IL-6, IL-6 blockers J o u r n a l P r e -p r o o f such as tocilizumab should be considered to diminish the inflammatory state and improve outcomes, as suggested in other studies (Gordon et al., 2021; Rubin et al., 2021) .

Interestingly, some of the cytokines, including CXCL10, IL-6, IL-8, TNF-alpha, INF-beta, and INFgamma, were correlated with the expression of IFN genes. The IFN pathway has been physiologically found to be a key point in the immune response (Calabrese et al., 2021; . In the present study, all the genes included in the IFN signature correlated with each other, and all the variables behaved similarly. Higher values of expression of these genes were related to more severe disease in both adults and children. This is consistent with the data reported before, in which SARS-CoV-2 infection was proven to cause an immune dysregulation with overproduction of cytokines (IL-1, IL-6, IL-2R, IL-10, and TNF-alpha) and a reduction in IFNγ expression in CD4 T cells, in severe cases Ong et al., 2020; Pedersen and Ho, 2020) . Therefore, anti-inflammatory treatments should be considered when it comes to severe diseases .

In contrast to the other IFN genes, which increase with the severity of the disease, low values of CIITA expression were linked to more severe disease and worse prognosis. CIITA is part of the IFN-stimulated immune response and it plays a role in the antiviral response as well as in antigen presentation. It can inhibit cathepsins, redirecting the virus into a degradation pathway in multivesicular bodies or lysosomes. Since it reduces viral transcription, assembly, and release, it plays a protective role (Bruchez et al., 2020; Butowt et al., 2020) . This could explain why patients from this cohort with lower CIITA values had more severe disease.

It is true that there are already some biomarkers and risk scores that might help to detect patients at high risk of developing a severe COVID-19(Bello-Chavolla et al., 2020; Gerotziafas et al., 2020) . However, this work adds valuable information, because it puts light on the physiopathology of the disease. Understanding the immune and inflammatory factors involved in SARS-CoV-2 might help to find other approaches and treatments. In concrete, CIITA might be J o u r n a l P r e -p r o o f a diagnostic tool, but it even might be a target in future research; since it seems a protective factor, increasing its activity might be useful for fighting against the disease. However, this is a preliminary data, thus these results need to be confirmed in further studies.

To sum up, there are different biomarkers that might be helpful to stratify patients with COVID-19 into risk cohorts. First, there is the lymphocyte count, since lymphopenia is related to severe disease. Secondly, a high cytokine response, including IL-6, may also be an indicator of severe disease, meaning that anti-inflammatory treatment such as tocilizumab should be considered. Finally, IFN response seems to be extremely important in managing SARS-CoV-2related infection. According to the findings of this study, the IFN pathway varies depending on the severity of the disease. Low CIITA values seem to be linked to more severe COVID-19; we propose as cut-off points for severe disease CITTA ≤-2.441 for adults and CITTA ≤-1.205 for children. In contrast, high values for SOCS1-and for all the other IFN genes were found in severe disease, and the proposed cut-off points are SOCS1>0.491 for adults and SOCS1>0.45 for children.

We reckon that these results need to be confirmed in additional studies with larger samples, but meanwhile these IFN genes may be taken into consideration when looking for therapeutic targets against SARS-CoV-2.

We acknowledge several limitations in this study. It is a single-centre study with a relatively small sample size. Therefore, the results may not be extrapolable to other populations.

Moreover, it would be extremely useful to examine the IFN expression and circulating cytokine profiles of the same patients at different time points of the disease, to better understand the overall immune response to it. Despite this, the results of the present study are promising, and further studies in line with this one should be performed so as to confirm them.

profile, were included. They were classified into different groups depending on the age and the severity of the disease (see Table S1 ). These different populations were analysed separately to compare their clinical characteristics and analytical patterns.

The study followed the Declaration of Helsinki recommendations and all subjects signed the informed consent form. The study was approved by the hospital ethics committee and institutional review board with the registry number PIC58-20.

The following demographic data were collected: age, gender, and comorbidities. Respiratory support used was also examined: oxygen, high flow nasal cannula, need for non-invasive ventilation (NIV), and mechanical ventilation (MV). Haemodynamic support requirement was analysed. The analytical biomarkers recorded were C-reactive protein, procalcitonin, midregional pro-adrenomedullin (MR-proADM), ferritin, D-dimer, lactate dehydrogenase, and blood counts. The length of stay in the hospital and in the PICU, as well as the mortality rate, were also recorded.

The sample analysis of IFN and cytokines was performed on admission, prior to any immunomodulatory treatment.

-IFN transcriptional profile: Whole blood samples were obtained in Tempus TM tubes (Applied Biosystems, California, USA) vigorously mixed and stored at -80ºC. Total RNA was then isolated using the Tempus Spin RNA Isolation Kit (Applied Biosystems, California, USA) following the manufacturer's recommendations. RNA was quantified by spectrophotometry (Nanodrop 2000, ThermoScientific, MA, USA) and assessed for J o u r n a l P r e -p r o o f quality using A260/A280 ratio. Gene expression of 32 genes was determined by Nanostring using the nCounter Elements Technology System according to the manufacturer's recommendations (NanoString Technologies, Seattle, WA). Twentyeight of the 32 genes have been previously described as "type-I IFN response genes" (IRG) (Kim et al., 2018) . In addition, we included 4 genes (CIITA, CXCL9, IFNA2, and STAT1) related to the IFN pathway but not strictly related with IFN-I response. Data were imported into nSolver analysis software (version 2.5) for quality checking and then were normalized, first to the geometric mean of six internal positive control sequences that are not native to any known organism, and second to the geometric mean of four housekeeping genes (ALAS1, HPRT1, TBP, and TUBB). A Z-score for each gene was calculated using the mean and standard deviation of a cohort of 26 healthy controls , which included 18 healthy adults and 8 healthy children, using the following equation for each gene:

The genes included were : CXCL10, DDX60, EPSTI1, GBP1, HERC5, HERC6, IFI27, IFI44,   IFI44L, IFI6, IFIT1, IFIT2, IFIT3, IFIT5, ISG15, LAMP3, LY6E, MX1, OAS1, OAS2, OAS3 Categorical variables were indicated as frequency (n) and percentage (%), whereas continuous variables were summarized as median and interquartile range (IQR) because they did not follow a normal distribution. Continuous variables were compared using the non-parametric

Kruskal-Wallis test. Probability values <0.05 were considered statistically significant.

Full correlation matrix between each pair of variables was computed. Since variables did not follow a normal distribution, correlations were analysed with Spearman's rank correlation coefficient. Once computed, the correlation matrix was plotted using a heat map where red indicates strong correlations (coefficient close to 1) and blue indicates negative correlations (coefficient close to -1).

With the objective of evaluating how well the patients could be classified based on the genetic variables collected, a machine learning supervised classifier was used. This type of method, when provided with a series of input variables and their known associated outcome class, automatically 'learns' the relationship between measured variables and outcome. Their output is then useful both to classify future cases and to quantify the relationships between inputs and clinical outcomes.

In this case, input variables were IFN gene variables values, and the class was the clinical outcome of the patients grouped in two classes: 1) asymptomatic and mild in one class and 2) moderate and severe in the other.

We used the classical CART algorithm as base classifier. This method maps the variables collected for each patient into their class in the shape of a decision tree. The tree is composed of nodes that branch to subsequent children nodes depending on the patient's variables. The J o u r n a l P r e -p r o o f tree is built in such a way that each branch separates the two classes as much as possible. This separation is measured as Shannon entropy, where a node with an entropy of zero means that the classification is perfect (either all or none of the patients developed a severe disease due to and an entropy of one is the worst possible mix (50%/50%). Once formed, each node contains a disease probability, calculated from the number of positive samples that reach the node, divided by the total number of samples. This approach has been widely used and validated in many applications and in other medical papers.

Due to the small sample size we wanted to avoid overfitting and sought to produce a simple, intuitive classifier that would explain the collected data in the simplest possible way.

Therefore, instead of directly applying the classifier using all input variables, we first applied a variable selection process. Variable selection is a process for selecting those variables which contribute most to the target outcome. These techniques allowed us to reduce overfitting, improve accuracy, and reduce training time, because having redundant or irrelevant data as input in classification models increases the probability of making a decision based on noise.

Variable importance is calculated as the decrease in node impurity weighted by the probability of reaching that node. The node probability can be calculated by the number of samples that reach the node, divided by the total number of samples. In this study we used the Extra Tree

Classifier. This algorithm for a given numerical attribute selects its cut-point fully at random, i.e., independently of the target variable. At each tree node, this is combined with a random choice of a certain number of attributes from which the best one is determined. Applying the variable selection algorithm, we obtained the list of variables sorted by importance in explaining the target variable. From the output of the variable selection algorithm, we chose only the two most relevant variables that did not correlate strongly; variables are correlated when the correlation coefficient is greater than 0.4 (in absolute value).

Finally, using the variables selected, the CART algorithm was used to build a tree classifier.

Since the goal of the study was to assess IFN response in children and evaluate the differences J o u r n a l P r e -p r o o f with respect to adults, different classifiers were built separately for children and adults, using the same variables. The classifiers were built using 10-fold cross-validation as recommended in TRIPOD (Transparent Reporting of a multivariable prediction model for Individual Prognosis or Diagnosis) guidance.

Once the tree classifiers were built, in order to make interpretability easier, in addition to the usual tree representation, we plotted the results in a bidimensional plot where each of the axes is one of the variables. The classification boundaries defined by the tree are shown with a red line and the patients are represented with a dot in different colours depending on their clinical outcome.

All statistical analysis was performed using Python version 3.7 (Python Foundation, USA).

The correlation coefficient was calculated using Spearman's rank correlation coefficient. From top to bottom and from left to right the variables plotted are : CXCL10, DDX60, EPSTI1, GBP1, HERC5,   HERC6, IFI27, IFI44, IFI44L, IFI6, IFIT1, IFIT2, IFIT3, IFIT5, ISG15, LAMP3, LY6E, OAS2, OAS3, OASL, RSAD2, RTP4, SIGLEC1, SOCS1, SPATS2L, USP18, CIITA, CXCL9, IFNA2,   STAT1, CXCL9, CXCL10, G-CSF, IFN-beta, IFN-gamma, IL-1RA, IL-2R alpha, IL-6, IL-8, IL-10, IL- 18, MCP-1, MCP-3, and TNF-alpha. A colour scale has been used to represent the correlations: red for positive correlations and blue for negative correlations. CXCL10, DDX60,   EPSTI1, GBP1, HERC5, HERC6, IFI27, IFI44, IFI44L, IFI6, IFIT1, IFIT2, IFIT3, IFIT5, ISG15,   LAMP3, LY6E, MX1, OAS1, OAS2, OAS3, OASL, RSAD2, RTP4, SIGLEC1, SOCS1, SPATS2L, USP18, CIITA, CXCL9, IFNA2, STAT1. 

Human genetic factors associated with susceptibility to SARS-CoV-2 infection and COVID-19 disease severity

Long COVID in children and adolescents

Auto-antibodies against type I IFNs in patients with life-threatening COVID-19. Science (80-. ). 370, eabd4585

Validation and repurposing of the MSL-COVID-19 score for prediction of severe COVID-19 using simple clinical predictors in a triage setting: The Nutri-CoV score

The COVID-19 immune landscape is dynamically and reversibly correlated with disease severity

MHC class II transactivator CIITA induces cell resistance to Ebola virus and SARS-like coronaviruses

Battle at the entrance gate: CIITA as a weapon to prevent the internalization of SARS-CoV-2 and Ebola viruses

Type I interferon, antiinterferon antibodies, and COVID-19

Clinical and immunological features of severe and moderate coronavirus disease 2019

The Immunology of Multisystem Inflammatory Syndrome in Children with COVID-19

Multisystem inflammatory syndrome in children and COVID-19 are distinct presentations of SARS-CoV-2

Similarities and differences between the immunopathogenesis of COVID-19-related pediatric multisystem inflammatory syndrome and Kawasaki disease

IL-6 serum levels predict severity and response to tocilizumab in COVID-19: An observational study

Risk factors for severe and critically ill COVID-19 patients: A review

Derivation and Validation of a Predictive Score for Disease Worsening

Impaired type I interferon activity and inflammatory responses in severe COVID-19 patients

Regulation of type I interferon responses

Development of a Validated Interferon Score Using NanoString Technology

Immunophenotyping of COVID-19 and influenza highlights the role of type I interferons in development of severe COVID-19

Systematic review of COVID-19 in children shows milder cases and a better prognosis than adults

Phenotyping long COVID

A Dynamic Immune Response Shapes COVID-19 Progression

Type I and Type III Interferons -Induction, Signaling, Evasion, and Application to Combat COVID-19

SARS-CoV-2: a storm is raging

Presenting Characteristics, Comorbidities, and Outcomes Among 5700 Patients Hospitalized With COVID-19 in the

Interleukin-6 Receptor Inhibition in Covid-19 -Cooling the Inflammatory Soup

COVID-19 in Children: Clinical Approach and Management

Clinical spectrum of COVID-19 and risk factors associated with severity in Spanish children

Presence of Genetic Variants Among Young Men With Severe COVID-19

COVID-19 Weekly Epidemiological Update, World Health Organization

Age-Related Differences in Immunological Responses to SARS-CoV-2

Multisystem inflammatory syndrome in children and adolescents with COVID-19

Inborn errors of type I IFN immunity in patients with lifethreatening COVID-19. Science (80-. ). 370, eabd4570

This study was supported by the projects PI18/00223, FI19/00208 and PI21/00211 to LA, integrated in the Plan Nacional de I+D+I and co-financed by the ISCIII -Subdirección General de Evaluación y Fomento de la Investigación Sanitaria -and the Fondo Europeo de Desarrollo Regional (FEDER), by Pla Estratègic de Recerca i Innovació en Salut (PERIS), Departament de