Skip to main content
  • Research article
  • Open access
  • Published:

Pathway mapping of leukocyte transcriptome in influenza patients reveals distinct pathogenic mechanisms associated with progression to severe infection



Influenza infections produce a spectrum of disease severity, ranging from a mild respiratory illness to respiratory failure and death. The host-response pathways associated with the progression to severe influenza disease are not well understood.


To gain insight into the disease mechanisms associated with progression to severe infection, we analyzed the leukocyte transcriptome in severe and moderate influenza patients and healthy control subjects. Pathway analysis on differentially expressed genes was performed using a topology-based pathway analysis tool that takes into account the interaction between multiple cellular pathways. The pathway profiles between moderate and severe influenza were then compared to delineate the biological mechanisms underpinning the progression from moderate to severe influenza.


107 patients (44 severe and 63 moderate influenza patients) and 52 healthy control subjects were included in the study. Severe influenza was associated with upregulation in several neutrophil-related pathways, including pathways involved in neutrophil differentiation, migration, degranulation and neutrophil extracellular trap (NET) formation. The degree of upregulation in neutrophil-related pathways were significantly higher in severely infected patients compared to moderately infected patients. Severe influenza was also associated with downregulation in immune response pathways, including pathways involved in antigen presentation such as CD4+ T-cell co-stimulation, CD8+ T cell and Natural Killer (NK) cells effector functions. Apoptosis pathways were also downregulated in severe influenza patients compare to moderate and healthy controls.


These findings showed that there are changes in gene expression profile that may highlight distinct pathogenic mechanisms associated with progression from moderate to severe influenza infection.

Peer Review reports


Infection with influenza virus is a significant cause of morbidity and mortality worldwide and is seen as one of the major global health threats of our time [1]. Influenza can cause a mild respiratory illness but can also lead to respiratory failure and death [2, 3]. Virulence factors (of different virus strains) contribute significantly to this variability in infection severity [4, 5]. However, it is increasingly being recognized that host factors also play an important role in contributing to infection severity [6,7,8,9,10]. Recently, circulating leukocytes have been identified as the main host factors linked to infection severity, as revealed by transcriptomics studies [11,12,13]. Transcriptomics studies capture global gene-expression changes expressed by circulating leukocytes and findings from these studies showed that influenza host response displayed distinctively changes across the full range of mild, moderate and severe infection [12,13,14,15,16]. However, further advances are impeded by a major limitation in these studies. Transcriptomics studies provide gene-level analysis; they offer limited insight into the underlying biological pathways affected by influenza infection. A pathway-level analysis approach will be more informative since it provides a considerably greater amount of biologically relevant information and thus will allow a better understanding of the pathogenic mechanisms linked to disease progression.

In this study, we aimed to gain insight into influenza pathogenesis on a pathway level. We used a topology-based pathway analysis technique to study the leukocytes-mediated host response of patients with influenza infection. Topology-based techniques have the advantage over other analytical techniques in that they incorporate the interactions of different genes and proteins, potentially revealing far more biologically relevant information than standard “gene-set” pathway analysis techniques. A recent study showed that topology-based pathway analysis can better model the biological phenomena and more realistically reflect the host response [17]. Here, we applied topology-based pathway analysis on blood transcriptomics data obtained from patients with moderate and severe influenza infection. Our objective was to delineate, from a pathway perspective, the host response mechanisms associated with disease progression from moderate to severe influenza infection.


Setting and patients

We performed a multicentre prospective study to recruit adult patients (> 18 years) fulfilling the World Health Organization criteria of influenza-like illness (fever of 38 °C or higher, cough and illness onset within the last 10 days). Only patients with a positive test for influenza virus were included. Severe influenza illness was defined as a severe influenza pneumonitis that requires mechanical ventilation. Moderate influenza illness was defined as a significant symptomatic disease (resulting consultation to an emergency department) but did not require mechanical ventilation support. Two physicians independently assigned the patients to groups (moderate vs. severity), based on the following criteria (1) whether mechanical ventilation was used, (2) that influenza virus was confirmed on PCR of airway samples, (3), that the clinical features are consistent with influenza illness. Patients infected with other viruses than influenza were excluded from the study. Healthy control subjects were also included. Study protocol was approved by the institutional review board of each participating institution. Informed written consents were obtained from all study participants.

Data collection

Demographic data included age, gender, comorbidities including chronic respiratory disease and standard laboratory test were obtained. Moreover, airway sample (nasopharyngeal swab, throat sample) and 2.5 ml peripheral blood sample in PAXgene tubes were collected. For those admitted in intensive care unit (ICU) and were under mechanical ventilation, additional respiratory samples were obtained. Virus testing consisted of nucleic acid PCR in order to detect influenza A, influenza B, respiratory syncytial virus, rhinovirus, parainfluenza virus and metapneumovirus.

RNA extraction, normalisation and microarray analysis

In each sample, whole blood RNA was extracted from PAXgene tubes as per manufacturer’s protocol (QIAGEN PreAnalytiX – Blood RNA version 2; 2015). After checking RNA integrity on Bioanalyzer (Agilent Technologies; Waldbronn, Germany), 100 ng of total RNA was applied for Cy3-labelling reaction using the one color Quick Amp labelling protocol (Agilent Technologies; Waldbronn, Germany). Labelled cRNA was hybridized to Agilent 8x60k Human V3 (Design ID: 072363) microarrays for 16 h at 68 °C and scanned using the Agilent DNA Microarray Scanner. Results were then analysed using the R software package (version 3.1.3). Pre-processing steps included background correction, adding an offset of 50, quantile normalization and annotation using the limma package and Agi4x44PreProcess packages. Multi-group comparisons and identification of differentially expressed probe sets were performed using limma with Benjamini and Hochberg correction for multiple testing. Differentially expressed probesets (some genes are represented by several probesets) were identified based on an adjusted p-value of < 0.05 and exhibiting more than a two-fold difference in expression levels ([log2] > 1). Full dataset of the expression data is available at the National Centre for Biotechnology Information Gene Expression Omnibus (GEO accession number GSE101702).

Statistical analysis

Continuous variables were expressed as mean and standard deviation (SD) and were compared by using Mann-Whitney U test. Categorical variables were expressed as numbers and percentages and were compared by the chi-square test or the Fisher exact test as appropriate.

Pathway analysis

Three lists of differentially expressed genes were generated by comparing three different phenotypes (Fig. 1):

  1. I.

    Moderate influenza vs. healthy controls

  2. II.

    Severe influenza vs. healthy controls

  3. III.

    Severe influenza vs. moderate influenza

Fig. 1
figure 1

Flow chart and study scheme. Flow chart shows study design and analysis workflow

Differentially expressed (DE) probesets between these comparison groups were identified based on an adjusted p-value of < 0.05 and exhibiting more than a two-fold difference in expression levels ([log2] > 1). Changes in expression levels were presented as fold changes for probesets of a given gene.

For pathway analysis, probeset lists were collapsed to gene lists (some genes are represented by multiple probesets that may be differentially expressed) and these three DE gene lists were used as input for the pathway analysis software, MetaCore. MetaCore is a topology-based analysis software (Clarivate Analytics, Philadelphia, U.S.A.) designed for functional analysis of transcriptomic data. Pathway analysis consists in examining the intersection between our lists of differentially expressed genes and prebuilt canonical pathways. A first statistical analysis is performed that take into account the number of genes differentially expressed in each comparison, the number of genes that intersect the prebuilt maps and the number of genes in the database. A False Discovery Rate (FDR) adjustment is apply for multiple testing. An FDR of 5% was used as the cut-off to determine whether a pathway was statistically overrepresented in the gene list. Adjusted p-value are expressed in -log(p-value) and ranked by statistical significance. Finally, these statistically overrepresented pathways were re-organized under broad biological themes (e.g. “Interferon response”, “Neutrophils”, “Immune response” and “Cell cycle”) [17,18,19].


Clinical data

One hundred seven patients were recruited with laboratory-confirmed influenza infection (either nasopharyngeal swab or bronchoalveolar lavage (BAL)). There was no viral coinfection observed in the cohort. Microarray analysis was performed for all 107 samples. The moderate group and the severe group, as defined by clinical criteria (see Methods), consisted of 63 and 44 patients respectively. Patient assignment matches perfectly between the 2 experts. A group of 52 healthy control subjects without any medical illnesses had also been enrolled (Fig. 1). Baseline characteristics are described in Table 1. Age, gender, cell counts did not differ significantly between the groups. Eighty five percent (n = 91) of the patients had at least one comorbidity and the proportion was not significantly different between moderate and severe cases (n = 51 (81%), n = 40 (91%), p value: 0.18). In the moderate group, 45 patients (71%) required hospitalization and 7 (11%) required ICU admission. Severity disease was associated with a longer length of stay, (1.4 days vs 26 days; p value< 0,0001). The hospital mortality rate in the severe cases of influenza pneumonitis was 20% (9/44 patients).

Table 1 Demographics and clinical characteristics of patients

Gene expression profile in severe influenza illness differs from moderate influenza illness

Influenza infection was associated with significant changes in gene expression. Compared to healthy control subjects, 994 transcripts from unique genes were found to be differentially expressed in severe influenza illness of which 535 were up-regulated and 459 down-regulated. Similarly, 252 transcripts from unique genes were differentially expressed in moderate influenza illness compared to healthy controls subjects of which 185 were up-regulated and 67 were down-regulated. Finally, the comparison between severe and moderate influenza illness revealed 211 transcripts (from unique genes) that were differentially expressed, of which 103 were up-regulated and 108 were down-regulated (Fig. 2a). Severe and moderate influenza illness share commonly expressed genes (147 up-regulated and 62 down-regulated). However, 388 genes were found to be up-regulated only in the severe influenza group and 38 genes only in the moderate influenza group. Conversely, 397 genes were found to be down-regulated only in the severe influenza group and 5 genes only in the moderate influenza group (Fig. 2b, c).

Fig. 2
figure 2

Differentially expressed genes in moderate and severe influenza. a Break down of statistically significant differentially expressed genes. The table showed the total number of differentially expressed genes in the three comparisons and the number of upregulated or downregulated genes. b Venn diagrams to indicate overlap of up-regulated genes. c Venn diagrams to indicate overlap of down-regulated genes. The Venn diagrams showed evidence that severe and moderate infection share common characteristics, but also have a distinctive gene expression profile

Unsupervised principal component analysis (PCA) was performed using normalized log2 gene-expression levels (Additional file 1: Fig. S1). Based on two principal components, the analysis showed a separation in the gene expression between severe influenza, moderate influenza and healthy control subjects. We note that gender didn’t seem to be associated with a separation in gene expression profile (Additional file 1: Fig. S1B).

Taken together, these findings suggest that severe and moderate infection share common characteristics, but they also have a distinctive gene expression profiles, indicating severe and moderate infection might be associated with distinct host response. This was confirmed by subsequent pathway analysis, which revealed four biological themes (Interferon response, Neutrophils, Immune response and Cell cycle) of which three (Neutrophils, Immune response and Cell cycle) are differentially expressed in the severe group compare to the moderate group (Fig. 3, Additional file 1: Tables S1, S2, S3 and Additional files 2, 3 and 4).

Fig. 3
figure 3

Top 10 pathways ranked by statistical significance. Top 10 pathways ranked by p-values (vertical bars) in three conditions. Vertical axis denotes statistical significance as measured by minus logarithm of p-values. Blue bars are downregulated pathways; red bars are upregulated pathways. a Moderate influenza patients compared to healthy controls. Upregulated pathways correspond to activation of interferon and neutrophil themes. Downregulated pathway corresponds to the immune response theme. b Severe influenza patients compared to healthy controls. Upregulated pathways correspond to activation of interferon and neutrophil themes. Downregulated pathway corresponds to the immune response theme. c Severe influenza patients compared to moderate influenza patients. Upregulated pathways correspond to activation of neutrophil and cell cycle (delayed apoptosis) themes. Downregulated pathway corresponds to the immune response theme

Interferon-related pathway genes

As expected, the most statistically significant pathways (in both moderate and severe groups) were those involved in antiviral response and interferon signalling. Notably, the interferon stimulated gene (ISG) IFI27 was the most up-regulated transcript in almost all infected patients. Other transcripts from interferon-stimulated genes, with anti-viral functions, were also up-regulated in both moderate and severe groups compared to controls including IFI44 (inhibits viral replication, p <  0.0001), IFIT1 and IFIT2 (inhibits viral mRNA translation initiation, p <  0.0001), MX1 (anti-influenza, p <  0.0001), OAS3 (degrades viral RNA, p <  0.0001). Conversely, we found only a small number of interferon stimulated genes (ISGs from the INTERFEROM database [20]) for the comparison between severe and moderate illness (Additional file 1: Figure S2). However, when we compared individual interferon-related genes, we observed a trend in downregulation of these genes in patients with severe illness (although these differences were statistically not significant) (Additional file 1: Figure S3). These findings showed that influenza infection up-regulates interferon-stimulated genes and other antiviral genes but the degree of up-regulation was not significantly different between severe and moderate illness.

Neutrophil-related pathway genes

We also identified neutrophils-related pathways as the most strongly associated with severe disease. These pathways were involved in key neutrophil process including neutrophils differentiation, degranulation and neutrophil extracellular traps (NETs) formation (Fig. 4). Furthermore, transcripts from genes involved in driving neutrophil differentiation from myeloid cell lineages were significantly up-regulated in severe influenza. Among them, RETN (resistin), also known as C/EBP-ε, is a critical transcription factor and was highly up-regulated in severe influenza patients (p <  0.0001) [21, 22]. RETN induce myeloid/granulocyte-specific genes expression such as PRTN3 (proteinase 3) MPO (myeloperoxidase), LCN2 (lipocalin 2) and LTF (lactotransferrin), which promote cell differentiation into mature neutrophils. Other transcripts from genes encoding for specific neutrophils granules were also more highly expressed in severe patients (MMP-9 (matrix metallopeptidase-9, p <  0.0001), HP (haptoglobin, p <  0.0001), OLM4 (olfactomedin, p <  0.0001). We noted that CD177 was the most abundant transcript in severe influenza illness compared to moderate illness and healthy control subjects (p <  0.0001). CD177 is a specific neutrophils protein that plays a role in neutrophils adhesion and transendothelial migration [23].

Fig. 4
figure 4

Histogram of neutrophil-related significant changes in gene expression between severe influenza, moderate influenza illness and healthy controls. Y-axis shows normalised log2 expression levels. * indicate p < 0.001, adjusted for multiple testing by Bonferroni method. ns denotes non-significant. HC denotes healthy control. a Genes encoding proteins involved in neutrophil extracellular trap formation. Expression differences are shown for (the strongest regulated) probesets of the individual gene. b Genes encoding proteins involved in neutrophil migration. Expression differences are shown for (the strongest regulated) probesets of the individual gene. c Genes encoding components of neutrophil granules. Expression differences are shown for (the strongest regulated) probesets of the individual gene. Neutrophils-related genes were upregulated in patients with severe influenza illness compared to moderate influenza illness and healthy control subjects

NETosis was also one of the most up-regulated pathway in severe influenza illness. Among transcripts from genes involved in NETs, all were significantly more expressed in severe influenza compared to moderate influenza (Fig. 3c). Furthermore, among the 30 most up-regulated transcripts in the severe group compared to the moderate group, 25 transcripts (83%) are from neutrophils related genes (Additional file 1: Table S4).

These findings suggest that, at the pathway level, neutrophil differentiation, activation, migration and finally NET formation are associated with a more severe disease in influenza patient.

Immune response-related pathway genes

Compared to healthy control subjects, data from influenza patients showed a downregulation of the immune response pathways. Among the top 10 ranking pathways of biological significance in the moderate group, one pathway was associated with downregulation of CD8+ T cell recruitment and cytotoxicity. In the severe influenza group, there was a stronger effect of down-regulation in immune pathways; six pathways from the top 10 ranking pathways of biological significance were down-regulated (p <  0.0001). These pathways are involved in both innate and adaptive immune response: CD4+ T cell co-stimulation, MHC class II expression and antigen presentation, CD8+ cytotoxicity and NK cell cytotoxicity (Fig. 3).

Comparing transcriptomes between severe and moderate influenza illness, we observed that the most down-regulated transcripts were from genes coding for Major histocompatibility complex (MHC) class II (HLA-DRB5, HLA-DRA1, MHC class II alfa and beta chain, HLA-DRB1, HLA-DRB4, HLA-DRB3, HLA-DRB) (p <  0.0001). Transcripts involved in the CD8+ T cell recruitment pathway also was found to be differentially expressed between moderate and severe illness but involved only a small number of transcripts (CX3CR1 downregulation, IL-18R1, ELANE, MMP-9 and SLPI upregulation) (Additional file 1: Table S3). These findings are summarised in Figs. 5 and 6.

Fig. 5
figure 5

Histogram of MHC class II significant changes in gene expression between severe influenza compared to moderate influenza illness and healthy controls. Expression differences are shown for (the strongest regulated) probesets of the individual gene. Y-axis shows normalised log2 expression levels. * indicate p < 0.001, adjusted for multiple testing by Bonferroni method. ns denotes non-significant. HC denotes healthy control MHC class II were downregulated in patients with severe influenza illness compared to moderate influenza illness and healthy control subjects

Fig. 6
figure 6

Radar chart: gene expression profile in moderate and severe influenza illness. Most representative differentially expressed probesets of the individual gene in the top three principal biological themes were computed in a radar chart. Level of expression was normalized to healthy controls. Expression differences are shown for (the strongest regulated) probesets of the individual gene. In moderate infection, the radar chart showed an upregulation of the interferon and neutrophil genes. In severe infection, besides an upregulation of the interferon and neutrophil genes, we observed a down regulation of the MHC class II genes

Cell cycle pathway genes

Cell cycle was also a distinct biological theme that we observed in severe patients. BIRC5 (baculoviral IAP repeat containing 5) also known as survivin was up regulated in severe patients and is known to play a role in apoptosis inhibition (p <  0.0001). CCNB2 (cyclin B2) was also up-regulated and encodes an essential component of the cell cycle regulatory machinery (p < 0.0001). Moreover, CDKN1C (Cyclin Dependant Kinase Inhibitor 1C) is a negative regulator of cell proliferation. We observed a downregulation of CDKN1C transcript in the severe influenza cases compared to moderate influenza cases (p < 0.0001). These findings suggest that cell proliferation and decrease in apoptosis could be associated with severity in influenza patients.

External validation

We also performed a gene-based comparison with an external dataset (GEO 111368) [12]. This external dataset has a similar study design to ours, which included 109 adult patients with laboratory-confirmed influenza infection and 130 healthy participants. We found strong similarity in the upregulated genes between the two studies: 17 upregulated genes were shared between two studies (70%) (Although only 3 (15%) downregulated genes were shared between the two studies). Notably, the upregulated genes encode proteins involved in neutrophil functions and the downregulated genes encode proteins involved in the immune response. Among the most differentially expressed genes, neutrophil activation and reduced immune response were the dominant biological themes found in both dataset (Additional file 1: Table S5).


A large number of studies have focused on the immunological mechanism underpinning host responses against influenza infection. These studies have identified risk factors for progression to severe influenza infection including viral factors (subtype of influenza virus, viral load, mutation in viral genome) and host factors (genetic susceptibility, pre-existing medical conditions) [7,8,9,10]. Host factors play a protective role against influenza infection, but also may contribute to immunopathology, leading to tissue damage, organ failure and disease severity.

We herein report a study of systemic host factors in a large cohort of patients presenting with influenza illness and well characterized different level of severity. Pathway mapping of leukocyte transcriptome in influenza patients has revealed distinct pathogenic mechanisms associated with severe infection. These mechanisms were represented by four biological functional themes associated with influenza illness (Neutrophils, Immune Response, Cell cycle and Interferon response). Three of the biological themes (Neutrophils, Immune Response, Cell cycle) were more highly associated with severe infection. Of these three themes, the neutrophils theme was the most represented in severe illness. These results were consistent with an unsupervised analysis previously published by our team [13]. In that analysis, a weighted gene co-expression analysis (WGCNA) was performed to identify disease modules associated with infection severity. The results of that analysis are in keeping with the main findings of the current paper, which include (1) the neutrophil module displayed the highest increase in modular expression as influenza severity progress from moderate to severe form, (2), cell cycle module was upregulated, (3) immune response module revealed broad downregulation in gene expression of key genes involved in innate and adaptive immunity. Overall, the pathway analysis presented in this paper extends the previous WGCNA analysis by providing additional insights on the biological pathways associated with severe influenza infection.

The neutrophils theme consisted of several key pathogenic pathways, including neutrophils differentiation, neutrophils migration and degranulation and neutrophil extracellular traps formation. These findings are consistent with previous studies highlighting the complex dual role of neutrophils in influenza infection. Neutrophils are key cells in innate immune responses that can be protective during influenza infection [24, 25]. However, mice and human studies showed that lung damage can also be associated with neutrophils infiltration, matrix metalloproteinase 9 and myeloperoxidase activities [26,27,28]. Interestingly, alveolar damage was less severe in neutrophils depleted mice infected with influenza A H1N1 illustrating the dual roles of neutrophils (protective vs. harmful) in influenza pneumonitis [29]. NETs were also strongly induced in neutrophils from infected mice lungs and were responsible for an increase in endothelial damage. Moreover, it has been shown that NET’s formation is dependent on redox enzymes activation such as myeloperoxidase and superoxide dismutase and blocking these enzymes can reduce alveolar damages [30]. In line with these findings, extracellular histones, major components of NET’s, play an important role in causing lung injury in mice infected by influenza, most probably mediated by cytotoxicity and thrombus formation after platelet binding [31]. Taken together, these data support the possible role of neutrophils in severe influenza immunopathology and suggest that more studies are needed to explore whether neutrophil can be a potential target of host-directed therapy [12].

Interferon has been identified as a critical key host response against influenza infection. First, inborn errors on immunity in IRF7 and IRF9 have been identified as responsible for life threatening pulmonary influenza. These patients cell fail to amplify type I and III IFNs and to control viral replication [7, 8]. Second, previous study in severely infected patients revealed an increase of neutrophils transcripts while interferon-related transcripts were down-regulated [12]. Finally, impaired INF production is responsible for reduced immune responses in lungs leading to acute lung injury [8, 10]. We found that interferon genes were not differentially expressed between moderate and severe infection. The lack of a further up-regulated interferon response observed in the severe group is consistent with the hypothesis that a reduced interferon response may contribute to compromised host response. Therefore, it is possible that an insufficient interferon response may be linked with a deregulated neutrophil activation. Whether this is indeed true will require further mechanistic study (e.g. in animal model).

Neutrophilia is well known marker of bacterial infection. In this study, super-imposed bacterial infections may explain the upregulation of neutrophil pathway observed in severe infection. Therefore, additional studies are required to investigate if the neutrophil signature was associated with bacterial co-infection.

Finally, a previous in vitro study suggests that H1N1 influenza virus could replicate in and be released from neutrophils. It has also been shown that influenza virus itself can activate neutrophils and induce oxidative burst [32]. Whether these mechanisms could explain our observations is uncertain and will require additional studies in the future.

Our findings also suggested significant changes in the immune response at a transcriptomic level, including, impairment of NK cells cytotoxicity, CD4 T cells co-stimulation, T cell CD8 recruitment and MHC II antigen presentation. Effective antigen presentation is critical to develop a robust adaptive immune response to the virus, and without this response virus clearance cannot occur effectively. Effective antigen presentation requires the cooperative interaction between MHC II molecules, CD74 protein from the endoplasmic reticulum (ER) and T-cell-receptor (TCR)-CD3 complex. This interaction is followed by binding of the CD4 co-receptor to the MHC II complex and this binding is a critical step in initiating the signalling pathways leading to T-cell activation and differentiation. In the current study, we observed that many transcripts involved in these processes were down-regulated in patients with severe influenza illness compare to moderate patients and healthy control. These findings are consistent with previous reports, showing that NK cell cytotoxicity and specific CD8+ T-cell response are essential for protection against severe influenza disease [33, 34]. In line with these findings, it has been shown that memory cross reactive CD8+ T Cell may provide protection in case of infection. Thus, patients a lower count of pre-existing cross reactive CD8+ T cell could develop more severe illness [33, 34]. Interestingly, ARG1 (arginase 1) was one of the most abundant transcripts found in severe influenza patients (Additional file 1: Fig. S5). ARG1 is known to be stored in granules of neutrophils. Once released and activated, ARG1 can degrade extra cellular arginine resulting in inhibition of T-cell proliferation [35].

CD8+ T-cell seems to play dual roles in influenza pathogenesis. If CD8+ cytotoxicity and the production of pro-inflammatory cytokines (INF-γ TNF-α) are critical for the efficiency of infection resolution in mice, they may also contribute to immunopathology and lung injury [36]. Our current study does not provide sufficient data to determine whether decreased CD8+ cytotoxicity could have contributed to immunopathology.

Apoptosis is also a well-known defence mechanism against viruses leading to the inhibition of the replication and the spread of virus into host [37]. Our data are consistent with previous work showing severe influenza infection results in the dysregulation of apoptosis in infected cells [38].

In summary, our transcriptomic data showed that, increase of neutrophils activity, decrease of antigen presentation and cytotoxicity (CD8+ T-cell and NK cells) are associated with severity while interferon pathway activation was common in all infected patients (Fig. 6).

Some limitations have to be acknowledged. An important limitation lies in the nature of the data we used. Our transcriptomic data requires further validation on both protein and functional levels. Also, it is important to keep in mind that our observations suggest biological and clinical association. However, establishing a causal link is not possible to do in humans but will require controlled animal experiments. Moreover, the immune response is a dynamic process and different pathways may be involved at different time points during an influenza infection. Also, patients arrive in hospital with different delay after onset. Although we found that the time elapsed since symptom did not statistically significantly impact on gene-expression levels, we cannot confidently exclude its potential confounding effect. Therefore, this issue should be clarified in future studies.


This study provides evidence that influenza infection severity is associated with changes in gene expression profile. It allows hypothesis-generation suggesting that excessive neutrophil activation, impaired adaptive immune function and apoptosis may be important host factors in mediating the progression to severe influenza disease. Further studies are needed to confirm these findings and to explore host mechanisms in particular in neutrophil-related host response.

Availability of data and materials

Full dataset of the gene-expression data is available at the National Centre for Biotechnology Information Gene Expression Omnibus (GEO accession number GSE101702). Additional data can be obtained by contacting our data manager by email:



Bronchoalveolar Lavage


Differentially Expressed


Gene Expression Omnibus


Intensive Care Unit


Interferon Stimulated Gene


Major Histocompatibility Complex


Neutrophil Extracellular Trap


Natural Killer


  1. WHO. WHO Influenza Fact Sheets 2016, Available from:

  2. Bermejo-Martin JF, Ortiz de Lejarazu R, Pumarola T, Rello J, Almansa R, Ramírez P, et al. Th1 and Th17 hypercytokinemia as early host response signature in severe pandemic influenza. Crit Care. 2009;13(6):R201.

    Article  Google Scholar 

  3. Mauad T, Hajjar LA, Callegari GD, da Silva LFF, Schout D, Galas FRBG, et al. Lung pathology in fatal novel human influenza a (H1N1) infection. Am J Respir Crit Care Med. 2010 Jan;181(1):72–9.

    Article  Google Scholar 

  4. Nieto A, Pozo F, Vidal-García M, Omeñaca M, Casas I, Falcón A. Identification of rare PB2-D701N mutation from a patient with severe influenza: contribution of the PB2-D701N mutation to the pathogenicity of human influenza. Front Microbiol. 2017;3(8):575.

    Article  Google Scholar 

  5. Park J-E, Ryu Y. Transmissibility and severity of influenza virus by subtype. Infect Genet Evol. 2018;65:288–92.

    Article  Google Scholar 

  6. Keynan Y, Malik S, Fowke KR. The role of polymorphisms in host immune genes in determining the severity of respiratory illness caused by pandemic H1N1 influenza. Public Health Genomics. 2013;16(1–2):9–16.

    Article  CAS  Google Scholar 

  7. Ciancanelli MJ, Huang SXL, Luthra P, Garner H, Itan Y, Volpi S, et al. Life-threatening influenza and impaired interferon amplification in human IRF7 deficiency. Science. 2015;348(6233):448–53.

    Article  CAS  Google Scholar 

  8. Hernandez N, Melki I, Jing H, Habib T, Huang SSY, Danielson J, et al. Life-threatening influenza pneumonitis in a child with inherited IRF9 deficiency. J Exp Med. 2018;215(10):2567–85.

    Article  CAS  Google Scholar 

  9. Allen EK, Randolph AG, Bhangale T, Dogra P, Ohlson M, Oshansky CM, et al. SNP-mediated disruption of CTCF binding at the IFITM3 promoter is associated with risk of severe influenza in humans. Nat Med. 2017;23(8):975–83.

    Article  CAS  Google Scholar 

  10. Ciancanelli MJ, Abel L, Zhang S-Y, Casanova J-L. Host genetics of severe influenza: from mouse Mx1 to human IRF7. Curr Opin Immunol. 2016;38:109–20.

    Article  CAS  Google Scholar 

  11. Zaas AK, Chen M, Varkey J, Veldman T, Hero AO, Lucas J, et al. Gene expression signatures diagnose influenza and other symptomatic respiratory viral infections in humans. Cell Host Microbe. 2009;6(3):207–17.

    Article  CAS  Google Scholar 

  12. Investigators MOSAIC, Dunning J, Blankley S, Hoang LT, Cox M, Graham CM, et al. Progression of whole-blood transcriptional signatures from interferon-induced to neutrophil-associated patterns in severe influenza. Nat Immunol. 2018;19(6):625–35.

    Article  Google Scholar 

  13. Tang BM, Shojaei M, Teoh S, Meyers A, Ho J, Blake T, et al. Neutrophils-related host factors are associated with disease progression and fatality in patients with severe influenza infection. Nat Commun. 2019;10(1):3422.

    Article  Google Scholar 

  14. Parnell GP, McLean AS, Booth DR, Armstrong NJ, Nalos M, Huang SJ, et al. A distinct influenza infection signature in the blood transcriptome of patients with severe community-acquired pneumonia. Crit Care. 2012;16(4):R157.

    Article  Google Scholar 

  15. Ramilo O, Allman W, Chung W, Mejias A, Ardura M, Glaser C, et al. Gene expression patterns in blood leukocytes discriminate patients with acute infections. Blood. 2007;109(5):2066–77.

    Article  CAS  Google Scholar 

  16. Tang BM, Shojaei M, Parnell GP, Huang S, Nalos M, Teoh S, et al. A novel immune biomarker IFI27 discriminates between influenza and bacteria in patients with suspected respiratory infection. Eur Respir J. 2017;49(6):1602098.

  17. Mitrea C, Taghavi Z, Bokanizad B, Hanoudi S, Tagett R, Donato M, et al. Methods and approaches in the topology-based analysis of biological pathways. Front Physiol. 2013;4:278.

    Article  Google Scholar 

  18. Jin L, Zuo X-Y, Su W-Y, Zhao X-L, Yuan M-Q, Han L-Z, et al. Pathway-based analysis tools for complex diseases: a review. Genomics Proteomics Bioinformatics. 2014;12(5):210–20.

    Article  Google Scholar 

  19. Tarca AL, Bhatti G, Romero R. A Comparison of Gene Set Analysis Methods in Terms of Sensitivity, Prioritization and Specificity. Chen L, editor. PLoS ONE. 2013 8(11):e79217.

  20. Rusinova I, Forster S, Yu S, Kannan A, Masse M, Cumming H, et al. INTERFEROME v2.0: an updated database of annotated interferon-regulated genes. Nucleic Acids Res. 2012;41(D1):D1040–6.

    Article  Google Scholar 

  21. Friedman AD. Transcriptional regulation of granulocyte and monocyte development. Oncogene. 2002;21(21):3377.

    Article  CAS  Google Scholar 

  22. Gombart AF, Kwok SH, Anderson KL, Yamaguchi Y, Torbett BE, Koeffler HP. Regulation of neutrophil and eosinophil secondary granule gene expression by transcription factors C/EBP and PU.1. Blood. 2003;101(8):3265–73.

    Article  CAS  Google Scholar 

  23. Bai M, Grieshaber-Bouyer R, Wang J, Schmider AB, Wilson ZS, Zeng L, et al. CD177 modulates human neutrophil migration through activation-mediated integrin and chemoreceptor regulation. Blood. 2017;130(19):2092–100.

    Article  CAS  Google Scholar 

  24. Tate MD, Deng Y-M, Jones JE, Anderson GP, Brooks AG, Reading PC. Neutrophils ameliorate lung injury and the development of severe disease during influenza infection. J Immunol. 2009;183(11):7441–50.

    Article  CAS  Google Scholar 

  25. White MR, Tecle T, Crouch EC, Hartshorn KL. Impact of neutrophils on antiviral activity of human bronchoalveolar lavage fluid. Am J Physiol-Lung Cell Mol Physiol. 2007;293(5):L1293–9.

    Article  CAS  Google Scholar 

  26. Damjanovic D, Small C-L, Jeyananthan M, McCormick S, Xing Z. Immunopathology in influenza virus infection: uncoupling the friend from foe. Clin Immunol. 2012;144(1):57–69.

    Article  CAS  Google Scholar 

  27. Pechous RD. With friends like these: the complex role of neutrophils in the progression of severe pneumonia. Front Cell Infect Microbiol. 2017;7:160.

    Article  Google Scholar 

  28. Drescher B, Bai F. Neutrophil in viral infections, friend or foe? Virus Res. 2013;171(1):1–7.

    Article  CAS  Google Scholar 

  29. Tumpey TM, Garcia-Sastre A, Taubenberger JK, Palese P, Swayne DE, Pantin-Jackwood MJ, et al. Pathogenicity of influenza viruses with genes from the 1918 pandemic virus: functional roles of alveolar macrophages and neutrophils in limiting virus replication and mortality in mice. J Virol. 2005;79(23):14933–44.

    Article  CAS  Google Scholar 

  30. Narasaraju T, Yang E, Samy RP, Ng HH, Poh WP, Liew A-A, et al. Excessive neutrophils and neutrophil extracellular traps contribute to acute lung injury of influenza pneumonitis. Am J Pathol. 2011;179(1):199–210.

    Article  CAS  Google Scholar 

  31. Ashar HK, Mueller NC, Rudd JM, Snider TA, Achanta M, Prasanthi M, et al. The role of extracellular histones in influenza virus pathogenesis. Am J Pathol. 2018;188(1):135–48.

    Article  CAS  Google Scholar 

  32. Zhang Z, Huang T, Yu F, Liu X, Zhao C, Chen X, et al. Infectious progeny of 2009 a (H1N1) influenza virus replicated in and released from human neutrophils. Sci Rep. 2015;5:17809.

    Article  CAS  Google Scholar 

  33. Sridhar S, Begom S, Bermingham A, Hoschler K, Adamson W, Carman W, et al. Cellular immune correlates of protection against symptomatic pandemic influenza. Nat Med. 2013;19(10):1305–12.

    Article  CAS  Google Scholar 

  34. Quinones-Parra S, Grant E, Loh L, Nguyen THO, Campbell K-A, Tong SYC, et al. Preexisting CD8+ T-cell immunity to the H7N9 influenza a virus varies across ethnicities. Proc Natl Acad Sci. 2014;111(3):1049–54.

    Article  CAS  Google Scholar 

  35. Rodriguez PC, Quiceno DG, Ochoa AC. L-arginine availability regulates T-lymphocyte cell-cycle progression. Blood. 2007;109(4):1568–73.

    Article  CAS  Google Scholar 

  36. Rygiel TP, Rijkers ESK, de Ruiter T, Stolte EH, van der Valk M, Rimmelzwaan GF, et al. Lack of CD200 enhances pathological T cell responses during influenza infection. J Immunol. 2009;183(3):1990–6.

    Article  CAS  Google Scholar 

  37. Sano R, Reed JC. ER stress-induced cell death mechanisms. Biochim Biophys Acta BBA - Mol Cell Res. 2013;1833(12):3460–70.

    Article  CAS  Google Scholar 

  38. Parnell G, McLean A, Booth D, Huang S, Nalos M, Tang B. Aberrant Cell Cycle and Apoptotic Changes Characterise Severe Influenza A Infection – A Meta-Analysis of Genomic Signatures in Circulating Leukocytes. Chen Z, editor. PLoS ONE. 2011;6(3):e17186.

Download references


Not Applicable.


This work was supported by the Nepean Institute of Critical Care Education and Research and the Marie Bashir Institute for Infectious Diseases and Biosecurity. Klaus Schughart was supported by intra-mural grants from the Helmholtz-Association (Program Infection and Immunity) and a start-up grant from the University of Memphis Tennessee Health Science Center. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Authors and Affiliations




Conceived the study; BT and KS. Performed the study; MS, AM, JH, TB, YK, AP, AK, AK, MN, AM. Analysed the data; KS, BT, EJ, YZ. Wrote the manuscript; YZ, EJ, BT. Review of manuscript; all co-authors. All authors have read and approved the manuscript.

Corresponding author

Correspondence to Yoann Zerbib.

Ethics declarations

Ethics approval and consent to participate

Approval for the study was obtained from Human Research Ethics Committee, Nepean Blue Mountains Local Health District. Informed written consent was obtained from all study participants.

Consent for publication

No applicable as no personal data was used in this manuscript.

Competing interests

There is no conflict of interest to declare.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Figure S1.

Principal component analysis. Figure S2. Number of differentially expressed genes and proportion of interferon-stimulated gene for the three comparisons. Figure S3. Boxplots of Interferon-stimulated genes-expression. Table S1. Top 10 pathways ranked by statistical significance - moderate influenza patients compared to healthy controls. Table S2. Top 10 pathways ranked by statistical significance - severe influenza patients compare to healthy controls. Table S3. Top 10 pathways ranked by statistical significance - severe influenza patients compare to moderate influenza patients. Table S4. Top 30 differentially expressed genes (ranked by expression levels) between moderate and severe infections. Table S5. a Top 20 most upregulated genes in both datasets. b Top 20 most downregulated genes in both datasets.

Additional file 2.

List of differentially expressed genes moderate vs healthy controls.

Additional file 3.

List of differentially expressed genes severe vs healthy controls.

Additional file 4.

List of differentially expressed genes severe vs moderate.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zerbib, Y., Jenkins, E.K., Shojaei, M. et al. Pathway mapping of leukocyte transcriptome in influenza patients reveals distinct pathogenic mechanisms associated with progression to severe infection. BMC Med Genomics 13, 28 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: