Skip to main content

Co-expression network analysis identifies potential candidate hub genes in severe influenza patients needing invasive mechanical ventilation



Influenza is a contagious disease that affects people of all ages and is linked to considerable mortality during epidemics and occasional outbreaks. Moreover, effective immunological biomarkers are needed for elucidating aetiology and preventing and treating severe influenza. Herein, we aimed to evaluate the key genes linked with the disease severity in influenza patients needing invasive mechanical ventilation (IMV). Three gene microarray data sets (GSE101702, GSE21802, and GSE111368) from blood samples of influenza patients were made available by the Gene Expression Omnibus (GEO) database. The GSE101702 and GSE21802 data sets were combined to create the training set. Hub indicators for IMV patients with severe influenza were determined using differential expression analysis and Weighted correlation network analysis (WGCNA) from the training set. The receiver operating characteristic curve (ROC) was also used to evaluate the hub genes from the test set's diagnostic accuracy. Different immune cells' infiltration levels in the expression profile and their correlation with hub gene markers were examined using single-sample gene set enrichment analysis (ssGSEA).


In the present study, we evaluated a total of 447 differential genes. WGCNA identified eight co-expression modules, with the red module having the strongest correlation with IMV patients. Differential genes were combined to obtain 3 hub genes (HLA-DPA1, HLA-DRB3, and CECR1). The identified genes were investigated as potential indicators for patients with severe influenza who required IMV using the least absolute shrinkage and selection operator (LASSO) approach. The ROC showed the diagnostic value of the three hub genes in determining the severity of influenza. Using ssGSEA, it has been revealed that the expression of key genes was negatively correlated with neutrophil activation and positively associated with adaptive cellular immune response.


We evaluated three novel hub genes that could be linked to the immunopathological mechanism of severe influenza patients who require IMV treatment and could be used as potential biomarkers for severe influenza prevention and treatment.

Peer Review reports


Influenza is a contagious respiratory disease caused by the influenza virus. Despite advances in medical technology, influenza continues to cause many hospitalizations and deaths [1, 2]. Symptomatic influenza affects 10–20% of the global population yearly, with 3–5 million cases of severe disease and 290,000–650,000 deaths reported [3]. Furthermore, influenza is a disease with a wide range of clinical manifestations, from a self-limited upper respiratory tract infection to severe pneumonia [4, 5]. With a mortality rate as high as 50–80%, critically ill patients with influenza frequently have severe respiratory failure (e.g., arterial pressure of oxygen/fraction of inspiration oxygen 200 mmHg) and require invasive mechanical ventilation (IMV) [6,7,8]. Moreover, the detailed pathogenesis of severe influenza is still unknown.

Previous studies have revealed that immunological cells and immune pathways have been implicated in the onset and progression of severe influenza [9, 10]. Hence, it is needed to evaluate promising immunological biomarkers for diagnosing and treating patients with severe influenza. Bioinformatics analysis has been frequently used to screen disease-specific biomarkers since the emergence of microarray technologies [11, 12]. However, it is challenging to conduct statistical evaluations and extract relevant information due to sample heterogeneity and sampling discrepancies, different technology platforms, and technique utilization in individual studies [13,14,15]. Hence, integrated bioinformatics tools that provide comprehensive and valuable information have been employed to investigate the molecular pathophysiology of influenza infection and evaluate novel biomarkers. Weighted correlation network analysis, also known as weighted gene co-expression network analysis (WGCNA), is a tool utilized for indicating gene interaction patterns across models, which can characterize the relationship between genes and interaction pathways based on the endogeneity of the gene set and the link between the gene set and the phenotype [16, 17]. Therefore, this promising method is increasingly being utilized to identify highly synergistic gene sets, potential biomarker genes, and therapeutic targets.

Herein, we aimed to evaluate the key genes linked with the severe influenza patients needing invasive mechanical ventilation (IMV) by WGCNA and their association with the infiltration levels of distinctive immune cells via single-sample gene set enrichment analysis (ssGSEA).


The stable DEGs detection and functional enrichment between IMV and NIMV group

Differentiation analysis was employed to examine the gene expression profiling from the training set. According to the criteria, between IMV and NIMV patients with influenza, a total of 447 DEGs (containing 261 upregulated and 186 downregulated genes) were found (Fig. 1).

Fig. 1
figure 1

Illustrates the training set gene expression profiling. A A heatmap of the top 50 DEGs. Upregulated genes are seen in red, while downregulated genes are highlighted in blue. B The DEGs volcanic plot. Upregulated genes are highlighted in red, while downregulated genes are indicated in green

The biological activities and signal cascades of DEGs co-related with severe influenza patients with IMV were studied using GO and KEGG analyses. Based on the GO enrichment analysis, DEGs were mostly involved in antigen processing and presentation (e.g., antigen processing and presentation of exogenous peptide antigen via MHC class II, MHC class II protein complex assembly, and peptide antigen assembly with MHC class II protein complex) and cell–cell adhesion (e.g., regulation of T cell activation), as shown in Fig. 2A. According to the KEGG signaling cascades analyses, DEGs were abundant in pathways linked to infection and cellular differentiation (e.g., Th1, Th2, and Th17 cell differentiation, Staphylococcus aureus infection, hematopoietic cell lineage, and human T-cell leukaemia virus-1 infection) (Fig. 2B). These findings demonstrated the immunological and inflammation-related pathways are associated with influenza severity in patients requiring IMV for severe influenza.

Fig. 2
figure 2

The DEGs from the training set were analyzed for functional enrichment Analysis of; A GO enrichment, and B Analysis of KEGG enrichment

WGCNA and hub gene analysis

The training set was utilized to build the co-expression network by WGCNA analysis to evaluate the important modules linked with influenza severity. The significance of the diagnosis and evaluation of hub gene expression level was examined using a soft-thresholding power of 6 (scale-free R2 = 0.85) and a cut height of 0.25. We then narrowed it down to eight modules, as shown in Fig. 3A-C. The link between module eigengene (ME) values and sample traits was used to quantify the link between the modules and clinical sample traits, which was observed using heat map profiling. The red module was the most closely related to illness severity (cor = 0.53, P = 7e-12) (Fig. 3D). Subsequently, three intersecting genes were determined based on DEG intersections (Fig. 3E). LASSO analysis confirmed 3 hub genes as follows: CECR1, HLA-DPA1, and HLA-DRB3 (Fig. 3F-G).

Fig. 3
figure 3

Establishing a weighted gene co-expression network analysis (WGCNA) and screening for Hub genes. A Scale-free fit index and mean connectivity analysis for various soft-thresholding powers. The red line denotes the point when the correlation coefficient is 0.9, and the soft-thresholding power (β) is 6. B Connectivity distribution histogram and scale-free topology check when β = 6. C Gene dendrograms and clustering of module eigengenes using a dissimilarity metric (1-TOM) (the red line represents a cut height of 0.25). D Module-trait correlations between module eigengenes and sample traits were assessed. The correlation coefficient and P value are displayed in each cell. E Venn diagram showing where the DEGs and the red module overlap. F LASSO regression's partial likelihood deviance with changing log(l) is shown in tenfold cross-validations. Using the minimum criterion (lambda.min) and 1 standard error of the minimum criterion (1-SE criteria), dotted vertical lines were created at the ideal values. G The tenfold cross-LASSO validation coefficient profiles for three hub genes are shown

Evaluation of hub gene expression and diagnostic value

Boxplots were utilized to identify the expression of the three hub genes. Figure 4A shows that IMV patients have considerably lower CECR1, HLA-DPA1, and HLA-DRB3 expression levels than NIMV patients from the training set. Moreover, similar results were observed in the patients from the test set (Fig. 4C). ROC analysis from the training set showed that area under the ROC curve (AUC) was 0.862 [95% confidence interval (CI) 0.800–0.921] for CECR1, 0.866 (95% CI 0.805–0.919) for HLA-DPA1, and 0.821 (95% CI 0.748–0.886) for HLA-DRB3, accordingly (Fig. 4B). Meanwhile, the AUC was 0.700–0.800 for the three hub genes from the test set, suggesting a moderate diagnostic accuracy (Fig. 4D).

Fig. 4
figure 4

Confirmation of hub genes. A Boxplots in the training set were applied to validate the hub genes' expression levels. B ROC analysis in the training set was used to validate the diagnostic utility of the hub genes. C Boxplots in the test set were used to verify the hub genes' expression levels, and (D) ROC analysis in the test set was utilized to establish the hub genes' diagnostic value

The link between hub genes and ICI

The ssGSEA algorithm was employed to evaluate the variances in ICI between IMV and NIMV groups. The results obtained from sGSEA revealed a considerably elevated infiltration of neutrophils and dendritic cells, a lower infiltration of activated CD8 T cells, memory CD8 T cells, memory CD4 T cells, and natural killer cells in the IMV group. Immunological cell correlation analysis with hub genes suggested that the three hub genes were positively linked with activated CD8 T-cells, T follicular helper cells, memory CD8 T-cells, and memory CD4 T-cells, while negatively associated with activated dendritic cells. In addition, HLA − DPA1 and HLA − DRB3 were negatively related to neutrophils (Fig. 5A, B).

Fig. 5
figure 5

Analysis of the immunological landscape concerning disease severity. The distribution of immune cells in the IMV and NIMV groups is depicted in an (A) heatmap and (B) violin plot. C The connection between immune cell infiltration and hub genes


Influenza is a contagious respiratory disease caused by the influenza virus. Despite advances in medical technology, influenza continues to cause many hospitalizations and deaths [1,2,3]. In recent decades, transcriptomics studies revealed that the structure of gene sets and their roles differed throughout the broad spectrum of patients with varying degrees of severity [13,14,15, 18, 19]. By combining numerous datasets and employing systematic bioinformatics tools, we identified three key genes related to severe influenza patients requiring IMV. WGCNA offers numerous advantages over other bioinformatics methods because the analysis focuses on the link between clinical features and co-expression modules, resulting in more complete data with high reliability and biological significance.

In this study, the functional enrichment analysis indicated that the DEGs between IMV and NIMV patients were mostly linked with immunological and inflammatory cascades, consistent with previous findings [18,19,20,21]. Moreover, the ICI showed two opposing results. First, the adaptive cellular immune response to influenza viruses was suppressed, as evidenced by a significant drop in active and memory CD8 + T cells in IMV patients, despite an increase in activated dendritic cells and Th1 cells. Second, the neutrophils were found to be activated. A similar study by Nguyen et al. [19] reported comparable results. They revealed a relationship between increased Sequential Organ Failure Assessment (SOFA) and lower adaptive interferon (IFN)-γ producing CD8 + T cell responses in people hospitalized with acute influenza. According to Jake Dunning et al., patients with the most severe illnesses also have a higher proportion of transcripts associated with neutrophils and fewer transcripts connected to IFN- γ [21]. Although multiple studies [22, 23] have shown a correlation between lymphopenia and severe influenza, the mechanism by which the adaptive cellular immune response is inhibited in severe influenza remains largely unknown [24,25,26]. When the lungs get infected with influenza, neutrophils are among the first cells attracted there and perform a defensive function [27]. In severe influenza infection, higher numbers of circulating and lung neutrophils are associated with, but not dependent on, bacterial co-infection [28, 29]. A significant neutrophil influx was associated with severe lung inflammation, whereas neutrophil reduction had little effect on viral clearance and caused only minor lung pathology [27,28,29]. Myeloperoxidase and neutrophil extracellular traps (NETs), which have been recognized as endothelial-damaging agents, are released by neutrophils [30, 31]. NETs produced by influenza infection may compromise lung function rather than contribute to bacterial killing or protect against secondary bacterial infection. According to recent studies, neutrophil-dependent tissue damage causes mortality regardless of viral load, and DNase neutralization of NETs enhances longevity [32].

Furthermore, poor outcomes in septic patients have previously been associated with decreased expression of MHC class II-related genes, including HLA-DPA1 and HLA-DRB3 [33,34,35]. Several inflammatory mediators may be responsible for the downregulation of MHC II mRNA expression [36, 37]. Rare accounts, nonetheless, did concentrate on influenza in the literature. We hypothesized that decreased MHC II expression results in impaired antigen processing and presentation, abnormal lymphocyte proliferation, and impaired viral clearance, all of which affect the development of severe influenza. The CECR1 (Cat Eye Syndrome Chromosome Region 1) gene, found on chromosome 22q11, also produces the ADA2 protein. According to previous research, ADA2 deficiency is frequently linked to lymphopenia (including CD8 memory cells, T cells, and follicular T cells) [38, 39]. It is also linked to increased myeloperoxidase expression and up-regulation of neutrophil-expressed gene transcripts, which can lead to the formation of NET [40]. In severe influenza, these could cause endothelium damage and severe illness. The current study indicated that the expression of three key genes was considerably lower in severe influenza patients requiring IMV compared to NIMV patients.

This study has some limitations. First, the sample size is still relatively small, which might limit the accuracy of the obtained data. Second, we should classify the association of hub genes and immune cells identified in this study as just a statistical correlation rather than causality. Finally, despite identifying several DEGs in IMV and NIMV patients, it is unknown if these host factors are unique to severe influenza infection. Hence, in vivo and in vitro studies are necessary to identify the function of identified hub genes and the underlying mechanism in severe influenza.


In conclusion, we evaluated three novel hub genes that could be linked to the immunopathological mechanism of severe influenza patients that require IMV treatment and could be used as candidate biomarkers for severe influenza prevention and treatment.


Data source

The NCBI-GEO ( was used to obtain the mRNA expression patterns of influenza patients' blood samples. The selection criteria included: i) influenza infection confirmed by RT-PCR using respiratory tract samples; ii) identical disease severity categorization; and iii) influenza patients ≥ 16 years of age and intubated. In this study, we identified three data sets such as GSE101702, GSE21802, and GSE111368. The GSE21802 microarray data contained 20 whole-blood samples from IMV patients and 16 whole-blood samples from non-IMV patients (NIMV), whereas the GSE101702 included 44 whole-blood samples from the IMV group and 63 whole-blood samples from the NIMV group. The microarray data of GSE111368 included 69 whole-blood samples from the IMV group and 160 whole-blood samples from the NIMV group. The training data for the GSE101702 and GSE21802 were combined to identify differentially expressed genes (DEGs) to construct WGCNA. The hub genes were validated using GSE111368 as the test set.

Evaluation of DEGs

The 'limma' and 'GEOquery' packages of R software (version 4.2.0) were used to normalize and annotate the data from the training set (GSE101702 + GSE21802), with DEG screening criteria of adjusted P-value < 0.05 and log fold change (logFC) > 0.5. The data was plotted on a volcano, and the top 50 DEGs were plotted on a heatmap.

Functional enrichment analysis

The R package 'cluster-profile' was utilized to conduct GO enrichment and KEGG pathway analyses. With an adjusted P < 0.05, GO terms or KEGG cascades were considered statistically significant. Biological process, cellular component, and molecular function (abbreviated as BP, CC, and MF, respectively) were GO terms' three aspects.

Generation of WGCNA

With the 'WGCNA' package of R software, a WGCNA was generated for the expression profile data of the training set. Next, genes were selected with the top 25% absolute deviation from the median. The 'goodSampleGenes' function was used to ensure that the data was accurate. With the help of the 'pickSoftThreshold' function, an ideal soft threshold (β) was chosen and confirmed. The data from the matrix were then changed into an adjacency matrix, clustered to find modules based on the topological overlap. A hierarchical clustering dendrogram was constructed after completing the module eigengene (ME) calculation and combining related modules in the ME-based clustering tree. Gene significance (GS) and module significance (MS) were evaluated using modules and phenotypic data to identify the GS and clinical information and examine the relationship between modules and models. In addition, for each gene, the module membership (MM) was determined to examine the GS in the module.

Hub gene evaluation and confirmation

Candidate hub genes were chosen based on their significant inter-module interaction. The absolute GS values of genes having biological importance are frequently greater. The criteria were used to screen candidate hub genes (absolute value of GS > 0.20; absolute value of MM > 0.80). The DEGs were then intersected with the potential hub genes using the R software's 'glmnet' package to conduct LASSO analysis to find the final hub genes.

Box plots were used to compare the levels of hub gene expression in the IMV and NIMV groups. The expression of hub genes was assessed through receiver operating characteristic curves (ROCs) to differentiate IMV from NIMV patients.

Evaluation of immune cell infiltration (ICI) and its association with hub genes

The ssGSEA algorithm was used to calculate the relative infiltration levels of immune cells in the training set samples. Differential expression levels of immune invading cells were visualized using violin plots. The 'ggplot2' package was used to visualize the Spearman correlations for immune infiltrating cells with hub genes.

Availability of data and materials

Publicly available datasets were analyzed in this study. These data can be found in GSE101702 (, GSE21802 ( and GSE111368 (



Invasive mechanical ventilation


Noninvasive mechanical ventilation


Gene Expression Omnibus


Receiver operating characteristic curve


Single-sample gene set enrichment analysis


Weighted gene co-expression network analysis


Least absolute shrinkage and selection operator


Differentially expressed genes


Sequential Organ Failure Assessment


Confidence interval


  1. Koszalka P, Subbarao K, Baz M. Preclinical and clinical developments for combination treatment of influenza. PLoS Pathog. 2022;18(5):e1010481.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Petrova VN, Russell CA. The evolution of seasonal influenza viruses. Nat Rev Microbiol. 2018;16(1):47–60.

    Article  CAS  PubMed  Google Scholar 

  3. Lozano R, Naghavi M, Foreman K, et al. Global and regional mortality from 235 causes of death for 20 age groups in 1990 and 2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2012;380(9859):2095–128.

    Article  PubMed  Google Scholar 

  4. Uyeki TM. Influenza. Annals of internal medicine. 2017;167(5):ITC33–48.

    Article  PubMed  Google Scholar 

  5. Influenza DD. Indian J Pediatr. 2020;87(10):828–32.

    Article  Google Scholar 

  6. Dominguez-Cherit G, De la Torre A, Rishu A, et al. Influenza A (H1N1pdm09)-Related Critical Illness and Mortality in Mexico and Canada, 2014. Crit Care Med. 2016;44(10):1861–70.

    Article  PubMed  Google Scholar 

  7. Oliveira EC, Marik PE, Colice G. Influenza pneumonia: a descriptive study. Chest. 2001;119(6):1717–23.

    Article  CAS  PubMed  Google Scholar 

  8. Chen L, Han X, Li Y, Zhang C, Xing X. Flu-IV score: a predictive tool for assessing the risk of invasive mechanical ventilation in patients with influenza-related pneumonia. BMC Pulm Med. 2022;22(1):47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Malik G, Zhou Y. Innate Immune Sensing of Influenza A Virus. Viruses. 2020;12(7):755.

    Article  CAS  PubMed Central  Google Scholar 

  10. Tomic A, Pollard AJ, Davis MM. Systems Immunology: Revealing Influenza Immunological Imprint. Viruses. 2021;13(5):948.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Ibrahim B, McMahon DP, Hufsky F, et al. A new era of virus bioinformatics. Virus Res. 2018;251:86–90.

    Article  CAS  PubMed  Google Scholar 

  12. Zou J, Wang E. Cancer Biomarker Discovery for Precision Medicine: New Progress. Curr Med Chem. 2019;26(42):7655–71.

    Article  CAS  PubMed  Google Scholar 

  13. Mathew NR, Jayanthan JK, Smirnov IV, et al. Single-cell BCR and transcriptome analysis after influenza infection reveals spatiotemporal dynamics of antigen-specific B cells. Cell Rep. 2021;35(12):109286.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Tang BM, Shojaei M, Teoh S, et al. Neutrophils-related host factors associated with severe disease and fatality in patients with influenza infection. Nat Commun. 2019;10(1):3422.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Parnell GP, McLean AS, Booth DR, 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  PubMed  PubMed Central  Google Scholar 

  16. Hu G, Grover CE, Arick MA, Liu M, Peterson DG, Wendel JF. Homoeologous gene expression and co-expression network analyses and evolutionary inference in allopolyploids. Brief Bioinform. 2021;22(2):1819–35.

    Article  CAS  PubMed  Google Scholar 

  17. Nomiri S, Karami H, Baradaran B, et al. Exploiting systems biology to investigate the gene modules and drugs in ovarian cancer: A hypothesis based on the weighted gene co-expression network analysis. Biomed Pharmacother. 2022;146:112537.

    Article  CAS  PubMed  Google Scholar 

  18. Tang BM, Shojaei M, Parnell GP, 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.

    Article  CAS  PubMed  Google Scholar 

  19. Nguyen THO, Koutsakos M, van de Sandt CE, et al. Immune cellular networks underlying recovery from influenza virus infection in acute hospitalized patients. Nat Commun. 2021;12(1):2691.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Zerbib Y, Jenkins EK, 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. 2020;13(1):28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Dunning J, Blankley S, Hoang LT, 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  CAS  PubMed  PubMed Central  Google Scholar 

  22. Lalueza A, Folgueira D, Diaz-Pedroche C, et al. Severe lymphopenia in hospitalized patients with influenza virus infection as a marker of a poor outcome. Infect Dis. 2019;51(7):543–6.

    Article  Google Scholar 

  23. Fox A, Le NM, Horby P, et al. Severe pandemic H1N1 2009 infection is associated with transient NK and T deficiency and aberrant CD8 responses. PLoS ONE. 2012;7(2):e31535.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Gonzalez Y, Juarez E, Carranza C, Sada E, Pedraza-Sanchez S, Torres M. Diminished effector and memory CD8+ circulating T lymphocytes in patients with severe influenza caused by the AH1N1 pdm09 virus. Virology. 2017;500:139–48.

    Article  CAS  PubMed  Google Scholar 

  25. Boonnak K, Vogel L, Feldmann F, Feldmann H, Legge KL, Subbarao K. Lymphopenia associated with highly virulent H5N1 virus infection due to plasmacytoid dendritic cell-mediated apoptosis of T cells. J Immunol. 2014;192(12):5906–12.

    Article  CAS  PubMed  Google Scholar 

  26. Yu J, Li H, Jia J, et al. Pandemic influenza A (H1N1) virus causes abortive infection of primary human T cells. Emerg Microbes Infect. 2022;11(1):1191–204.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Tate MD, Brooks AG, Reading PC. The role of neutrophils in the upper and lower respiratory tract during influenza virus infection of mice. Respir Res. 2008;9:57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Narasaraju T, Yang E, Samy RP, 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  PubMed  PubMed Central  Google Scholar 

  29. Tumpey TM, Garcia-Sastre A, Taubenberger JK, 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  PubMed  PubMed Central  Google Scholar 

  30. Yipp BG, Kubes P. NETosis: how vital is it? Blood. 2013;122(16):2784–94.

    Article  CAS  PubMed  Google Scholar 

  31. Narayana Moorthy A, Narasaraju T, Rai P, et al. In vivo and in vitro studies on the roles of neutrophil extracellular traps during secondary pneumococcal pneumonia after primary pulmonary influenza infection. Front Immunol. 2013;4:56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Siegler BH, Altvater M, Thon JN, et al. Postoperative abdominal sepsis induces selective and persistent changes in CTCF binding within the MHC-II region of human monocytes. PLoS ONE. 2021;16(5):e0250818.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Siegler BH, Uhle F, Lichtenstern C, et al. Impact of human sepsis on CCCTC-binding factor associated monocyte transcriptional response of Major Histocompatibility Complex II components. PLoS ONE. 2018;13(9):e0204168.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Lukaszewicz AC, Grienay M, Resche-Rigon M, et al. Monocytic HLA-DR expression in intensive care patients: interest for prognosis and secondary infection prediction. Crit Care Med. 2009;37(10):2746–52.

    Article  CAS  PubMed  Google Scholar 

  35. Fumeaux T, Pugin J. Role of interleukin-10 in the intracellular sequestration of human leukocyte antigen-DR in monocytes during septic shock. Am J Respir Crit Care Med. 2002;166(11):1475–82.

    Article  PubMed  Google Scholar 

  36. Cazalis MA, Friggeri A, Cave L, et al. Decreased HLA-DR antigen-associated invariant chain (CD74) mRNA expression predicts mortality after septic shock. Crit Care. 2013;17(6):R287.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Pillai PS, Molony RD, Martinod K, et al. Mx1 reveals innate pathways to antiviral resistance and lethal influenza disease. Science. 2016;352(6284):463–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Pinto B, Deo P, Sharma S, Syal A, Sharma A. Expanding spectrum of DADA2: a review of phenotypes, genetics, pathogenesis and treatment. Clin Rheumatol. 2021;40(10):3883–96.

    Article  PubMed  Google Scholar 

  39. Caorsi R, Penco F, Grossi A, et al. ADA2 deficiency (DADA2) as an unrecognised cause of early onset polyarteritis nodosa and stroke: a multicentre national study. Ann Rheum Dis. 2017;76(10):1648–56.

    Article  CAS  PubMed  Google Scholar 

  40. Belot A, Wassmer E, Twilt M, et al. Mutations in CECR1 associated with a neutrophil signature in peripheral blood. Pediatr Rheumatol Online J. 2014;12:44.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


The authors would like to express their gratitude to BMCSCI ( for the expert linguistic services provided.


This study is founded by the Liyang City 2020 Science and Technology Project (LC2020002).

Author information

Authors and Affiliations



LC: conception, design, data collection, data analysis, manuscript writing, manuscript revision, and supervision; JH: conception, design, data collection; XP H: conception, design, data collection, data analysis, manuscript writing, manuscript revision. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Xiaopu He.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Chen, L., Hua, J. & He, X. Co-expression network analysis identifies potential candidate hub genes in severe influenza patients needing invasive mechanical ventilation. BMC Genomics 23, 703 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Severe influenza
  • Invasive mechanical ventilation
  • Co-expression network analysis
  • Hub gene