- Research
- Open access
- Published:
Bioinformatics analysis identifies a key gene HLA_DPA1 in severe influenza-associated immune infiltration
BMC Genomics volume 25, Article number: 257 (2024)
Abstract
Background
Severe influenza is a serious global health issue that leads to prolonged hospitalization and mortality on a significant scale. The pathogenesis of this infectious disease is poorly understood. Therefore, this study aimed to identify the key genes associated with severe influenza patients necessitating invasive mechanical ventilation.
Methods
The current study utilized two publicly accessible gene expression profiles (GSE111368 and GSE21802) from the Gene Expression Omnibus database. The research focused on identifying the genes exhibiting differential expression between severe and non-severe influenza patients. We employed three machine learning algorithms, namely the Least Absolute Shrinkage and Selection Operator regression model, Random Forest, and Support Vector Machine-Recursive Feature Elimination, to detect potential key genes. The key gene was further selected based on the diagnostic performance of the target genes substantiated in the dataset GSE101702. A single-sample gene set enrichment analysis algorithm was applied to evaluate the participation of immune cell infiltration and their associations with key genes.
Results
A total of 44 differentially expressed genes were recognized; among them, we focused on 10 common genes, namely PCOLCE2, HLA_DPA1, LOC653061, TDRD9, MPO, HLA_DQA1, MAOA, S100P, RAP1GAP, and CA1. To ensure the robustness of our findings, we employed overlapping LASSO regression, Random Forest, and SVM-RFE algorithms. By utilizing these algorithms, we were able to pinpoint the aforementioned 10 genes as potential biomarkers for distinguishing between both cases of influenza (severe and non-severe). However, the gene HLA_DPA1 has been recognized as a crucial factor in the pathological condition of severe influenza. Notably, the validation dataset revealed that this gene exhibited the highest area under the receiver operating characteristic curve, with a value of 0.891. The use of single-sample gene set enrichment analysis has provided valuable insights into the immune responses of patients afflicted with severe influenza that have further revealed a categorical correlation between the expression of HLA_DPA1 and lymphocytes.
Conclusion
The findings indicated that the HLA_DPA1 gene may play a crucial role in the immune-pathological condition of severe influenza and could serve as a promising therapeutic target for patients infected with severe influenza.
Background
Despite advances in biomedicine, the incidence of hospitalization and mortality rates elicited by influenza, a profoundly contagious respiratory disease, persistently exhibit an upward trend [1, 2]. The global prevalence of symptomatic flu is estimated to range from 10 to 20% annually, affecting a substantial portion of the population. However, the chronic manifestation of this disease afflicts approximately 3–5 million individuals worldwide. Tragically, influenza-related mortality rates vary between 290,000 and 650,000 deaths [3]. Furthermore, the clinical manifestations of influenza encompass a diverse array of symptoms, from acute upper respiratory tract infections to the development of severe pneumonia [4]. Conversely, some patients afflicted with severe influenza frequently exhibit respiratory dysfunction, as evidenced by reduced arterial pressure of oxygen to a fraction of the inspired oxygen ratio (≦ 200 mmHg). Consequently, these patients rely on IMV for respiratory support, and critical patients’ death rate reaches about 50–80% [5, 6]. The detailed mechanisms governing the pathological condition of severe influenza remain elusive.
Previous studies reported that immune cells and pathways are pivotal to the occurrence and progression of severe influenza [7, 8]. Reliable immunological biomarkers are urgently required to prevent and treat patients with severe influenza infections. Microarray technologies and bioinformatic analyses have been widely used to identify disease-specific biomarkers [7, 9]. However, due to the presence of sample heterogeneity and variations in sampling methods, as well as the utilization of diverse technology platforms and analysis strategies across individual studies, the execution of statistical analyses and the extraction of esteemed information pose significant challenges.
Hence, the integration of bioinformatics approaches together with expression profiling techniques presents an opportunity to obtain a comprehensive understanding of the molecular mechanisms underlying influenza infection. This approach can yield valuable insights and facilitate the development of novel molecular signatures. Here, we elucidated the key genes implicated in the requirement of IMV among influenza patients through bioinformatics analysis. Additionally, we sought to investigate the association among these genes and the levels of infiltrating distinct immune cells. The study design can be seen in Fig. 1.
Results
Cross‐platform normalization
The microarray platforms collectively identified 12,031 genes in the two patient samples. Before applying batch-effect removal techniques, the samples displayed clustering patterns influenced by batch effects along the two principal component (PC) axes with the highest variance. These axes were determined using gene expression values that had not been normalized (Fig. 2a). The principal component analysis (PCA) analysis conducted after normalization has validated the effective removal of batch effects (Fig. 2b), demonstrating the successful implementation of cross-platform normalization.
DEGs identification and functional analyses
A total of 44 DEGs between both types of influenza (severe and non-severe) samples were identified from the training dataset, including six downregulated DEGs and 38 upregulated DEGs (Fig. 3a, b). GO and KEGG enrichment analyses were employed to elucidate the specific biological roles played by the DEGs in severe influenza. GO analysis suggested that the DEGs were associated with the process of myeloid leukocyte activation, defence response to bacteria, and regulation of cytokine production (Fig. 3c and Supplementary File 1). KEGG enrichment analysis exhibited that the DEGs were predominantly involved in the pathways of transcriptional misregulation in cancer, neutrophil extracellular trap formation, and the IL − 17 signaling pathway (Fig. 3d and Supplementary File 2). In summary, the DEGs were mainly involved in immune and inflammatory responses.
Identification of the key gene for severe influenza
The ten common DEGs (PCOLCE2, HLA_DPA1, LOC653061, TDRD9, MPO, HLA_DQA1, MAOA, S100P, RAP1GAP, and CA1) that were obtained by overlapping genes from computing the three algorithms [LASSO regression (Fig. 4a, b), SVM-RFE algorithms (Fig. 4c, d), and RF (Fig. 4e, f) are candidate key genes for severe influenza (Fig. 4g). HLA_DPA1 and HLA_DQA1 expression was significantly lower in patients with severe influenza in the training dataset compared to patients with non-severe influenza. In contrast, PCOLCE2, TDRD9, MPO, MAOA, RAP1GAP, and S100P expression was higher in the severe influenza group compared to the non-severe influenza group in the training cohort (Fig. 5a-h), similar to the findings in the validation cohort (Fig. 6a-h). The expression of LOC653061 and CA1 was greater in severe influenza patients in contrast to non-severe patients in the training dataset (Fig. 5i, j), whereas it was comparable in the validation dataset. In the training dataset, HLA_DPA1 and LOC653061 genes exhibited the highest AUC of 0.788 as depicted in Fig. 7a, b, while others were below 0.7 (Fig. 7c-j). Conversely, in the validation dataset, the AUCs of HLA_DPA1 and PCOLCE2 were 0.891 and 0.838, respectively (Fig. 8a, b). The AUCs of all eight candidate genes were found to be less than 0.7. Thus, HLA_DPA1 was selected as a key gene in patients diagnosed with severe influenza needing IMV.
The severe influenza samples were categorized into two distinct groups by employing a division based on the median value of HLA_DPA1 expression: HLA_DPA1low (n = 44) and HLA_DPA1high (n = 45). Some genes (e.g., SPOCK2, ITGB7, GIMAP5, et al.) were upregulated, while others (e.g., PFKFB2, IRAK3, SIPA1L2, et al.) were downregulated in the group of HLA_DPA1high (Fig. 9a, b). The correlation between HLA_DPA1 and the other genes in the training dataset is shown in Fig. 9c. The median expression level of HLA_DPA1 from the training dataset in severe and non-severe influenza patients was 9.540 and 10.572, respectively.
Identification of the key gene via GSEA and GSVA analyses
In order to understand the possible functional importance of HLA_DPA1 in the pathogenesis of severe influenza, single-gene GSEA-KEGG pathway analysis was executed (Supplementary File 3), and the top six pathways enriched for HLA_DPA1 are presented in Fig. 10a. Overall, HLA_DPA1 was found to be involved in the pathological condition of severe influenza by regulating the immune or inflammatory responses such as KEGG_leishmania_infection, KEGG_Toll_like_receptor_signaling_pathway), carbohydrate and cofactor metabolism, and vitamin metabolism. The GSVA produced comparable outcomes (Fig. 10b).
Analysis of infiltration of immune cells
Significance variances in the numbers of specific immune cell populations in whole blood samples from individuals with HLA_DPA1low and HLA_DPA1high were compared using ssGSEA. This approach revealed noteworthy suppressed adaptive immune responses in patients with HLA_DPA1low. This suppression was characterized by reduced levels of CD8 + T-cells, B-cells, two T-cell subsets (Th1-cells and Th2-cells), tumor-infiltrating lymphocytes (TIL), T-cell co-stimulation, antigen-presenting cell (APC) co-stimulation, as well as elevated levels of regulatory T-cells (Treg) and APC co-inhibition (Fig. 10c and Supplementary file 4). Similarly, suppressive adaptive immune responses were observed in patients with severe influenza, which manifested as decreased levels of key lymphocyte populations, including activated CD8 + T cells, B cells, CD4 + T cells, and memory CD8 + T cells, B cells, and CD4 + T cells (Fig. 10d). In addition, subsequent correlation data exhibited a remarkable positive association between the expression of HLA_DPA1 and the abundance of these lymphocytes (Fig. 10e and Supplementary File 5).
Establishment of a key gene-based ceRNA network
A comprehensive analysis was executed by intersecting genes from the TargetScan, miRDB, and miRanda databases (Supplementary File 6) and via this approach, six miRNAs (hsa-miR-573, hsa-miR-1253, hsa-miR-877-3p, hsa-miR-429, hsa-miR-3182, and hsa-miR-22-5p) targeting HLA_DPA1 were screened. Based on starBase, three lncRNAs (LINC00689, LINC00940, and RP1-253P7.1) interacted with hsa-miR-877-3p. A ceRNA network comprising 5 nodes and 4 edges was established (Fig. 11).
qRT-PCR
The mRNA levels of HLA_DPA1 in blood samples from patients afflicted with severe and non-severe conditions were verified using qRT-PCR. This showed a significant reduction in the expression of HLA_DPA1 in patients afflicted with severe influenza compared to those who remain non-severe by infection (Fig. 12).
Discussion
Previous investigations have elucidated the host factors linked to the development of severe influenza. However, they have predominantly concentrated on a genetic event, genetic susceptibility [10,11,12]. Recently, transcriptomic investigations have documented comprehensive gene expression profiles pertaining to the host's response. The findings from these investigations suggest that the composition and functionality of gene sets deviate significantly among patients exhibiting different degrees of severity [13, 14]. Nonetheless, these findings were derived solely from a singular cohort study, thereby necessitating additional clinical validation and comprehensive functional analysis that needed to be explored. Thus, we have successfully recognized the key genes associated with severe influenza in the current study by integrating multiple datasets. Consequently, the outcomes obtained are anticipated to offer a more comprehensive understanding of the subject matter. Three distinct machine-learning methods were employed for the screening of potential key genes. The LASSO is a widely recognized regression analysis algorithm renowned for its distinctive variable selection and regularization features. These attributes are instrumental in mitigating the risk of overfitting and enhancing the accuracy of predictions [15]. The Support Vector Machine (SVM) is a well-established supervised machine learning approach that is commonly employed for classification and regression tasks. On the other hand, the Recursive Feature Elimination (RFE) algorithm is utilized to identify the most optimal combination of variables that maximizes the performance of the model [16]. Hence, the current investigation utilized the Support Vector Machine Recursive Feature Elimination (SVM-RFE) algorithm to ascertain feature biomarkers possessing exceptional discriminative capacity. The Random Forest technique is a widely used regression tree-based method that employs bootstrap aggregation and predictor randomization to attain notable predictive accuracy [17].
The candidate key genes obtained by overlapping the genes from the three algorithms exhibited higher reliability. Our study’s functional enrichment analysis displayed that DEGs between both influenza (severe and non-severe) cases were primarily associated with pathways with immune response and inflammation-related pathways. Moreover, the ICI analysis revealed a notable impairment in adaptive immune responses among patients afflicted with severe influenza, consistent with prior scientific findings [13, 18, 19]. Nguyen et al. [13] conducted a longitudinal study on patients hospitalized with acute influenza and found that a higher SOFA score was associated with lower adaptive-producing CD8 + T cell responses. Dunning et al. [18] reported that patients with the most severe illness exhibited a notable reduction in interferon (IFN)-related transcripts. The precise mechanisms responsible for inhibiting adaptive cellular immune responses during severe influenza infection remain poorly elucidated. The occurrence and progression of adaptive cellular immunosuppression may involve various mechanistic events, including directive killing, disruption of antigen presentation, apoptosis, abortive infection of primary human T cells, and T cell exhaustion or paralysis induced by viruses and cytokines [20,21,22].
From the candidate genes, HLA_DPA1 was selected as the key gene for patients with severe influenza requiring IMV, which showed the best differential performance in both the training and validation cohorts. Functional enrichment analysis suggested that HLA_DPA1 mainly participates in regulating immune and inflammatory pathways. HLA_DPA1 was significantly and positively associated with lymphocytes; thus, the patients with HLA_DPA1low often showed deficient adaptive immunity and were more likely to be classified as critically ill. HLA_DPA1 is a major histocompatibility complex (MHC) class II-related gene [23]. HLA-DP-restricted T-cells and antimicrobial immune responses have also been identified [24, 25]. HLA-DPA1 polymorphism is a major determinant of hepatitis B virus clearance [26, 27]. A previous study reported that downregulation of HLA_DPA1 is associated with immunosuppression and increased mortality in sepsis [28,29,30]. In the context of severe infection, some inflammatory mediators are possibly involved in the down-regulation of the gene expression of MHC II [31,32,33,34]. For example, interleukin-10 (IL-10) can reduce the membrane expression of MHC II in monocytes. This reduction is attributed to the internalization and sequestration of mature MHC II molecules within the intracellular compartments [31, 32]. In an in vitro study, transforming growth factor-1 (TFG-1) downregulates MHC II mRNA expression by suppressing transcription factor class II transactivator (CIITA) mRNA transcription, while prostaglandin E2 was found to suppress MHC II mRNA expression in macrophages [33, 34]. The downregulation of MHC II leads to defective antigen processing, presentation, and as well as the proliferation of lymphocytes [35, 36]. The immunosuppressive state of the immune system significantly impedes the patient's ability to eliminate the primary influenza virus infection and enhances vulnerability to subsequent opportunistic infections, thereby resulting in many detrimental clinical outcomes in patients afflicted with influenza infection.
The present study has several noteworthy constraints. First, we must recognize the complex pathology of severe influenza, which is not driven by a single gene. Nevertheless, it can be asserted with a certain degree of certainty that the HLA_DPA1 gene exerts a pivotal influence on the progression of severe influenza and therefore merits prioritization in subsequent investigations. Second, the sample size was comparatively small despite our efforts to retrieve all the online data. Hub gene-encoding protein tests revealed a correlation between hub genes and disease severity. Furthermore, it is noted that the association between hub genes and immune cells is based on statistical correlation rather than establishing a causal relationship. Lastly, identifying DEGs in patients with both types of influenza has shed light on potential host factors associated with the chronicity of infection. However, the specificity of these factors to severe influenza infection has yet to be determined. Additional cell culture and animal studies are necessary to investigate these hub genes' roles and underlying mechanisms in severe influenza.
Conclusions
In conclusion, the findings of our investigation declare that the HLA_DPA1 gene act as a crucial role in the immunopathological condition of severe influenza. Furthermore, because of the high discrimination potency and cost-efficient property of HLA_DPA1, its clinical assessment may provide an accurate and early diagnosis of severe influenza. Therefore, it is a promising candidate for targeted interventions for the management and prevention of severe influenza cases necessitating IMV.
Materials and methods
Data source
The National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database, accessible at (http://www.ncbi.nlm.nih.gov/geo) serves as a comprehensive repository for mRNA expression data pertaining to patients affected with influenza. The selection criteria employed in this study were as follows: I) Influenza infection was confirmed through the application of reverse transcription polymerase chain reaction (RT-PCR) methodology, which involved the analysis of respiratory tract samples; and ii) the disease severity classification was generally similar. In this investigation, the classification of severe influenza was established depending on the criterion of patients necessitating IMV; iii) Influenza patients were ≧ 16 years old, and intubated patients were included. Three datasets were obtained: GSE21802, GSE111368, and GSE101702. The GSE21802 microarray data consisted of blood samples obtained from 20 patients with severe influenza and 16 patients diagnosed with non-severe influenza, and the GSE111368 dataset comprised 69 samples of severe and 160 samples of non-severe influenza cases. The dataset GSE101702 included blood samples obtained from 107 individuals, consisting of 44 patients diagnosed with severe influenza and 63 with non-severe influenza. After the elimination of mRNA probes from the GSE21802 and GSE111368 datasets, the gene expression analysis was consolidated into a unified file, serving as the training dataset.
Data processing and screening of differentially expressed genes
The integration of genomic data batches to increase statistical power is often hindered by batch effects or unwanted variation in data caused by differences in technical factors across batches. To remove the batch effect from different platforms and batches, the R sva package (https://bioconductor.org/packages/sva/) was employed to mitigate batch effects. Before conducting cross-platform normalization, the expression values of individual datasets underwent log2 transformation. Expression values obtained from various platforms or sample batches were subjected to normalization via the ComBat method. Principal component analysis was executed to validate the successful removal of batch effects. We used specific criteria to identify DEGs among both types of influenza (severe and non-severe) cases. The threshold points for selection were set at a significance P < 0.05 level and a minimum log fold change (logFC) > 1. The experimental findings were graphically represented using a volcano plot.
Functional enrichment analyses
The enrichment analyses for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were executed via the R package 'clusterprofiler'. The significance threshold for these analyses was set with an adjusted FDR (false discovery rate) (FDR < 0.05) and P-value < 0.05. GO terms were categorized into three main classes: biological process (BP), molecular function (MF), and cellular component (CC). In this study, we presented the top 10 enriched terms.
Candidate key genes identification
Three machine learning algorithms, least absolute shrinkage and selection operator (LASSO), Random Forest (RF), and Support Vector Machine (SVM), were utilized in this study to detect significant diagnostic genes for severe influenza. The LASSO is a regression analysis algorithm, which is characterized by variable selection and regularization. It helps avoid overfitting and improves the prediction accuracy. RF uses different independent decision trees to predict the classification or regression. The SVM is a supervised machine learning technique widely used in classification and regression. The recursive feature elimination (RFE) algorithm is employed to acquire the optimal combination of variables that maximizes the performance of the model. Therefore, this study utilized the SVM-RFE algorithm to identify potent biomarkers with superior discriminative ability. Thus the candidate genes will have higher reliability as they are identified by overlapping genes via three algorithms. To validate their expression levels in severe influenza samples, the dataset GSE101702 was utilized.
Diagnostic performance examination
To assess the predictive efficiency of the candidate key genes for severe influenza, an ROC curve was plotted using the mRNA expression data obtained from patients diagnosed with influenza (severe and non-severe), sourced from both the training and validation datasets. The gene exhibiting the highest area under the ROC curve within the validation cohort was identified as a key gene.
Patients with key gene expression values above the median for all severe influenza patients were categorized as the genehigh group. In contrast, those with values below the median were assigned to the genelow group. The differential expression of the key gene was determined using analysis of an unpaired t-test, with a significance level of P < 0.05. A fold change (FC, log2) threshold of > 0.5 or < -0.5 was also applied.
Pathway evaluation by single-gene gene set enrichment analysis
The R GSEA package was utilized to conduct GSEA to identify the pathways linked to the key genes. This was achieved by assessing the correlations between the key genes and all other genes in the training dataset.
These genes were then ranked based on the strength of their correlative relationships. The “c2.cp.kegg.Hs.symbols” gene set was downloaded from the MSigDB database for GSEA analysis and an |NES|> 1, normalized p-value < 0.05, and FDR q-value < 0.25 denoted statistical significance. The genes were subsequently ranked according to the magnitude of their correlative associations. The gene set "c2.cp.kegg.Hs.symbols" was obtained from the Molecular Signatures Database (MSigDB) to conduct GSEA. Statistical significance was determined based on the criteria of an absolute Normalized Enrichment Score |NES|> 1, a normalized p-value < 0.05, and an FDR q-value < 0.25.
Single-gene gene set variation analysis of key genes
The GSVA analyses of key genes were executed using the R GSVA package, with the KEGG pathway gene set as the background. Using the Limma package, a comparison of the GSVA scores for marker genes between the low- and high-expression groups was conducted. Significance variations between groups were evaluated via a threshold of |t|> 2 and a level of significance (P < 0.05). A positive value of t > 0 indicated pathway activation in the high-expression group, while a negative value of t < 0 indicated pathway activation in the low-expression group.
Correlation between the key gene and infiltrating immune cells
The calculation of relative ICI levels in the training dataset was executed utilizing a ssGSEA algorithm. Immune cell enrichment levels were quantified using ssGSEA scores for each sample. Differential expression patterns of immune-infiltrating cells between the key genehigh and key genelow groups, and patients with both cases of influenza (severe and non-severe), were monitored via violin plots. The Spearman correlations between ICI and the key gene were assessed via the 'ggplot2' package in the R programming language.
Development of ceRNA network
The identification of miRNAs that interact with key genes was performed using the StarBase computational tool. The mRNA sequences of these genes were obtained from NCBI. Human miRNA sequences were acquired from miRbase. Subsequently, the TargetScan, miRDB, and miRanda databases were employed to forecast the target genes of miRNA. StarBase was used to conduct screening for interactions between mRNA-lncRNA. This facilitated the establishment of a comprehensive network involving mRNA, microRNA (miRNA), and lncRNA.
qRT-PCR
Total RNA content was extracted from a set of 10 paired severe and non-severe influenza samples by the reagent of TRIzol (Life Technologies, Carlsbad, CA, USA) as per the manufacturer's protocol guidelines. The reverse transcription process was executed via PrimeScript RT Master Mix (Takara in Tokyo, Japan). The resulting cDNA was amplified using the ABI 7700 system (Applied Biosystems in CA, USA). β-lactin was employed as housekeeping control to evaluate the relative expression levels. It was assessed by utilizing the 2-ΔΔCt method. The following primer sequences were used for the qRT-PCR:
-
Forward: 5’-CTGCCCAGAACAGATTACAGC-3’,
-
Reverse: 5’-ACAGTCTCCGTTGTCTCAGG-3’
Data analysis
The statistical analyses were executed by applying R software (version 4.2.0). Statistical analysis was performed using an unpaired t-test for variables that revealed a normal distribution. At the same time, the Mann–Whitney U test was utilized for variables that displayed a non-normal distribution. Spearman's correlation coefficient was employed to conduct the correlation analysis. Statistical significance was determined by assessing differences with a p < 0.05.
Availability of data and materials
Publicly available datasets were analyzed in this study. These data can be found in GSE111368 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE111368), GSE101702 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE101702), GSE70866 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE70866) and GSE10667 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi).
Abbreviations
- GEO:
-
Gene Expression Omnibus
- DEGs:
-
Differentially Expressed Genes
- LASSO:
-
Least Absolute Shrinkage and Selection Operator
- SVM-RFE:
-
Random Forest, and Support Vector Machine-Recursive Feature Elimination
- ssGSEA:
-
Single-sample Gene Set Enrichment Analysis
- GSVA:
-
Gene Set Variation Analysis
- ROC:
-
Receiver Operating Characteristic Curve
- AUC:
-
Area Under the Curve
References
Petrova VN, Russell CA. The evolution of seasonal influenza viruses. Nat Rev Microbiol. 2018;16(1):47–60.
Koszalka P, Subbarao K, Baz M. Preclinical and clinical developments for combination treatment of influenza. PLoS Pathog. 2022;18(5):e1010481.
Uyeki TM. Influenza. Ann Intern Med. 2017;167(5):ITC33–48.
Chan KKP, Hui DSC. Antiviral therapies for influenza. Curr Opin Infect Dis. 2023;36(2):124–31.
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.
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.
Tomic A, Pollard AJ, Davis MM. Systems immunology: revealing influenza immunological imprint. Viruses. 2021;13(5):948.
Malik G, Zhou Y. Innate immune sensing of influenza A virus. Viruses. 2020;12(7):755.
Ibrahim B, McMahon DP, Hufsky F, et al. A new era of virus bioinformatics. Virus Res. 2018;251:86–90.
Bryson KJ, Sives S, Lee HM, Borowska D, Smith J, Digard P, et al. Comparative analysis of different inbred chicken lines highlights how a hereditary inflammatory state affects susceptibility to avian influenza virus. Viruses. 2023;15(3):591.
Li J, Wagatsuma K, Sun Y, Sato I, Kawashima T, Saito T, et al. Factors associated with viral RNA shedding and evaluation of potential viral infectivity at returning to school in influenza outpatients after treatment with baloxavir marboxil and neuraminidase inhibitors during 2013/2014-2019/2020 seasons in Japan: an observational study. BMC Infect Dis. 2023;23(1):188.
Mizuguchi M, Shibata A, Kasai M, Hoshino A. Genetic and environmental risk factors of acute infection-triggered encephalopathy. Front Neurosci. 2023;17:1119708.
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.
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.
Rafique R, Islam SMR, Kazi JU. Machine learning in the prediction of cancer therapy. Comput Struct Biotechnol J. 2021;19:4003–17.
Rjoob K, Bond R, Finlay D, McGilligan V, Leslie SJ, Rababah A, Iftikhar A, Guldenring D, Knoery C, McShane A, et al. Machine learning and the electrocardiogram over two decades: time series and meta-analysis of the algorithms, evaluation metrics and applications. Artif Intell Med. 2022;132:102381.
Sarica A, Cerasa A, Quattrone A. Random forest algorithm for the classification of neuroimaging data in alzheimer’s disease: a systematic review. Front Aging Neurosci. 2017;9:329.
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.
Ludwig S, Pleschka S, Planz O. MEK inhibitors as novel host-targeted antivirals with a dual-benefit mode of action against hyperinflammatory respiratory viral diseases. Curr Opin Virol. 2023;59:101304.
Yunis J, Short KR, Yu D. Severe respiratory viral infections: T-cell functions diverging from immunity to inflammation. Trends Microbiol. 2023;31:644.
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.
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.
Snyder JA, Weston A, Tinkle SS, Demchuk E. Electrostatic potential on human leukocyte antigen: implications for putative mechanism of chronic beryllium disease. Environ Health Perspect. 2003;111(15):1827–34.
de Waal L, Yuksel S, Brandenburg AH, Langedijk JP, Sintnicolaas K, Verjans GM, et al. Identification of a common HLA-DP4-restricted T-cell epitope in the conserved region of the respiratory syncytial virus G protein. J Virol. 2004;78(4):1775–81.
Gaston JS, Life PF, van der Zee R, Jenner PJ, Colston MJ, Tonks S, et al. Epitope specificity and MHC restriction of rheumatoid arthritis synovial T cell clones which recognize a mycobacterial 65 kDa heat shock protein. Int Immunol. 1991;3(10):965–72.
An P, Winkler C, Guan L, O’Brien SJ, Zeng Z, Consortium HBVS. A common HLA-DPA1 variant is a major determinant of hepatitis B virus clearance in Han Chinese. J Infect Dis. 2011;203(7):943–7.
Hosaka T, Suzuki F, Kobayashi M, Fukushima T, Kawamura Y, Sezaki H, et al. HLA-DP genes polymorphisms associate with hepatitis B surface antigen kinetics and seroclearance during nucleot(s)ide analogue therapy. Liver Int. 2015;35(4):1290–302.
Chen Z, Chen R, Ou Y, Lu J, Jiang Q, Liu G, et al. Construction of an HLA classifier for early diagnosis, prognosis, and recognition of immunosuppression in sepsis by multiple transcriptome datasets. Front Physiol. 2022;13:870657.
Lu J, Chen R, Ou Y, Jiang Q, Wang L, Liu G, et al. Characterization of immune-related genes andimmune infiltration features for early diagnosis, prognosis and recognition of immunosuppression in sepsis. Int Immunopharmacol. 2022;107:108650.
Cazalis MA, Friggeri A, Cave L, Demaret J, Barbalat V, Cerrato E, et al. Decreased HLA-DR antigen-associated invariant chain (CD74) mRNA expression predicts mortality after septic shock. Crit Care. 2013;17(6):R287.
Kasraie S, Niebuhr M, Kopfnagel V, Dittrich-Breiholz O, Kracht M, Werfel T. Macrophages from patients with atopic dermatitis show a reduced CXCL10 expression in response to staphylococcal alpha-toxin. Allergy. 2012;67(1):41–9.
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.
Lee KS, Baek DW, Kim KH, Shin BS, Lee DH, Kim JW, et al. IL-10-dependent down-regulation of MHC class II expression level on monocytes by peritoneal fluid from endometriosis patients. Int Immunopharmacol. 2005;5(12):1699–712.
De Lerma BA, Procopio FA, Mortara L, Tosi G, Accolla RS. The MHC class II transactivator (CIITA) mRNA stability is critical for the HLA class II gene expression in myelomonocytic cells. Eur J Immunol. 2005;35(2):603–11.
Xu Y, McDonald J, Perloff E, Buttice G, Schreiber BM, Smith BD. Collagen and major histocompatibility class II expression in mesenchymal cells from CIITA hypomorphic mice. Mol Immunol. 2007;44(7):1709–21.
Mulder DJ, Pooni A, Mak N, Hurlbut DJ, Basta S, Justinich CJ. Antigen presentation and MHC class II expression by human esophageal epithelial cells: role in eosinophilic esophagitis. Am J Pathol. 2011;178(2):744–53.
Acknowledgements
The authors would like to express their gratitude to BMCSCI (http://www.bmcscience.com/) for the expert linguistic services provided.
Funding
This study is supported by the Nanjing medical science and technology development fund (NO. YKK22239) and Foundation for Advanced Talents (NO.08) (received by Liang Chen).
Author information
Authors and Affiliations
Contributions
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 revision.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Supplementary Material 1.
Detailed results of GO analysis.
Supplementary Material 2.
Detailed results of KEGG analysis.
Supplementary Material 3.
Detailed results of single-gene GSEA-KEGG pathway analysis for HLA_DPA1.
Supplementary Material 4.
Detailed results of ssGSEA for individuals with HLA_DPA1low and HLA_DPA1high.
Supplementary Material 5.
Correlation analysis of HLA_DPA1 expression level with the abundance of lymphocytes
Supplementary Material 6.
Predicted genes from the TargetScan, miRDB, and miRanda databases.
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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Chen, L., Hua, J. & He, X. Bioinformatics analysis identifies a key gene HLA_DPA1 in severe influenza-associated immune infiltration. BMC Genomics 25, 257 (2024). https://doi.org/10.1186/s12864-024-10184-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12864-024-10184-7