Construction of an HLA Classifier for Early Diagnosis, Prognosis, and Recognition of Immunosuppression in Sepsis by Multiple Transcriptome Datasets

Background: Sepsis is a clinical syndrome, due to a dysregulated inflammatory response to infection. Accumulating evidence shows that human leukocyte antigen (HLA) genes play a key role in the immune responses to sepsis. Nevertheless, the effects of HLA genes in sepsis have still not been comprehens...

Бүрэн тодорхойлолт

Номзүйн дэлгэрэнгүй
Үндсэн зохиолчид: Zhen Chen, Rui Chen, Yangpeng Ou, Jianhai Lu, Qianhua Jiang, Genglong Liu, Liping Wang, Yayun Liu, Zhujiang Zhou, Ben Yang, Liuer Zuo
Формат: Өгүүллэг
Хэл сонгох:English
Хэвлэсэн: Frontiers Media S.A. 2022-05-01
Цуврал:Frontiers in Physiology
Нөхцлүүд:
Онлайн хандалт:https://www.frontiersin.org/articles/10.3389/fphys.2022.870657/full
_version_ 1828730451283410944
author Zhen Chen
Rui Chen
Yangpeng Ou
Jianhai Lu
Qianhua Jiang
Genglong Liu
Liping Wang
Yayun Liu
Zhujiang Zhou
Ben Yang
Liuer Zuo
author_facet Zhen Chen
Rui Chen
Yangpeng Ou
Jianhai Lu
Qianhua Jiang
Genglong Liu
Liping Wang
Yayun Liu
Zhujiang Zhou
Ben Yang
Liuer Zuo
author_sort Zhen Chen
collection DOAJ
description Background: Sepsis is a clinical syndrome, due to a dysregulated inflammatory response to infection. Accumulating evidence shows that human leukocyte antigen (HLA) genes play a key role in the immune responses to sepsis. Nevertheless, the effects of HLA genes in sepsis have still not been comprehensively understood.Methods: A systematical search was performed in the Gene Expression Omnibus (GEO) and ArrayExpress databases from inception to 10 September 2021. Random forest (RF) and modified Lasso penalized regression were conducted to identify hub genes in multi-transcriptome data, thus we constructed a prediction model, namely the HLA classifier. ArrayExpress databases, as external validation, were utilized to evaluate its diagnostic, prognostic, and predictive performance. Immune cell infiltration score was calculated via CIBERSORTx tools and single-sample gene set enrichment analysis (ssGSEA). Gene set variation analysis (GSVA) and ssGSEA were conducted to determine the pathways that are significantly enriched in different subgroups. Next, we systematically correlated the HLA classifier with immunological characteristics from multiple perspectives, such as immune-related cell infiltration, pivotal molecular pathways, and cytokine expression. Finally, quantitative real-time polymerase chain reaction (qRT-PCR) was conducted to validate the expression level of HLA genes in clinical samples.Results: A total of nine datasets comprising 1,251 patients were included. Based on RF and modified Lasso penalized regression in multi-transcriptome datasets, five HLA genes (B2M, HLA-DQA1, HLA-DPA1, TAP1, and TAP2) were identified as hub genes, which were used to construct an HLA classifier. In the discovery cohort, the HLA classifier exhibited superior diagnostic value (AUC = 0.997) and performed better in predicting mortality (AUC = 0.716) than clinical characteristics or endotypes. Encouragingly, similar results were observed in the ArrayExpress databases. In the E-MTAB-7581 dataset, the use of hydrocortisone in the HLA high-risk subgroup (OR: 2.84, 95% CI 1.07–7.57, p = 0.037) was associated with increased risk of mortality, but not in the HLA low-risk subgroup. Additionally, immune infiltration analysis by CIBERSORTx and ssGSEA revealed that B cells, activated dendritic cells, NK cells, T helper cells, and infiltrating lymphocytes (ILs) were significantly richer in HLA low-risk phenotypes, while Tregs and myeloid-derived suppressor cells (MDSCs) were more abundant in HLA high-risk phenotypes. The HLA classifier was significantly negatively correlated with B cells, activated dendritic cells, NK cells, T helper cells, and ILs, yet was significantly positively correlated with Tregs and MDSCs. Subsequently, molecular pathways analysis uncovered that cytokine-cytokine receptor (CCR) interaction, human leukocyte antigen (HLA), and antigen-presenting cell (APC) co-stimulation were significantly enriched in HLA low-risk endotypes, which was significantly negatively correlated with the HLA classifier in multi-transcriptome data. Finally, the expression levels of several cytokines (IL-10, IFNG, TNF) were significantly different between the HLA subgroups, and the ratio of IL-10/TNF was significantly positively correlated with HLA score in multi-transcriptome data. Results of qRT-PCR validated the higher expression level of B2M as well as lower expression level of HLA-DQA1, HLA-DPA1, TAP1, and TAP2 in sepsis samples compared to control sample.Conclusion: Based on five HLA genes, a diagnostic and prognostic model, namely the HLA classifier, was established, which is closely correlated with responses to hydrocortisone and immunosuppression status and might facilitate personalized counseling for specific therapy.
first_indexed 2024-04-12T17:30:27Z
format Article
id doaj.art-9ee5491b51e54b9194d9cb0c36f23ba0
institution Directory Open Access Journal
issn 1664-042X
language English
last_indexed 2024-04-12T17:30:27Z
publishDate 2022-05-01
publisher Frontiers Media S.A.
record_format Article
series Frontiers in Physiology
spelling doaj.art-9ee5491b51e54b9194d9cb0c36f23ba02022-12-22T03:23:09ZengFrontiers Media S.A.Frontiers in Physiology1664-042X2022-05-011310.3389/fphys.2022.870657870657Construction of an HLA Classifier for Early Diagnosis, Prognosis, and Recognition of Immunosuppression in Sepsis by Multiple Transcriptome DatasetsZhen Chen0Rui Chen1Yangpeng Ou2Jianhai Lu3Qianhua Jiang4Genglong Liu5Liping Wang6Yayun Liu7Zhujiang Zhou8Ben Yang9Liuer Zuo10Department of Intensive Care Unit, Shunde Hospital, Southern Medical University (The First People’s Hospital of Shunde), Foshan, ChinaDepartment of Medical Intensive Care Unit, General Hospital of Southern Theater Command, Guangzhou, ChinaDepartment of Oncology, Huizhou Third People’s Hospital, Guangzhou Medical University, Huizhou, ChinaDepartment of Intensive Care Unit, Shunde Hospital, Southern Medical University (The First People’s Hospital of Shunde), Foshan, ChinaDepartment of Intensive Care Unit, Shunde Hospital, Southern Medical University (The First People’s Hospital of Shunde), Foshan, ChinaDepartment of Pathology, The Third Affiliated Hospital of Guangdong Medical University (Longjiang Hospital of Shunde District), Foshan, ChinaDepartment of Otorhinolaryngology Head and Neck Surgery, The First Affiliated Hospital of Hainan Medical University, Haikou, ChinaDepartment of Endocrinology, GuiYang Huaxi District People’s Hospital, Guiyang, ChinaDepartment of Intensive Care Unit, Shunde Hospital, Southern Medical University (The First People’s Hospital of Shunde), Foshan, ChinaDepartment of Burn Surgery, Huizhou Municipal Central Hospital, Huizhou, ChinaDepartment of Intensive Care Unit, Shunde Hospital, Southern Medical University (The First People’s Hospital of Shunde), Foshan, ChinaBackground: Sepsis is a clinical syndrome, due to a dysregulated inflammatory response to infection. Accumulating evidence shows that human leukocyte antigen (HLA) genes play a key role in the immune responses to sepsis. Nevertheless, the effects of HLA genes in sepsis have still not been comprehensively understood.Methods: A systematical search was performed in the Gene Expression Omnibus (GEO) and ArrayExpress databases from inception to 10 September 2021. Random forest (RF) and modified Lasso penalized regression were conducted to identify hub genes in multi-transcriptome data, thus we constructed a prediction model, namely the HLA classifier. ArrayExpress databases, as external validation, were utilized to evaluate its diagnostic, prognostic, and predictive performance. Immune cell infiltration score was calculated via CIBERSORTx tools and single-sample gene set enrichment analysis (ssGSEA). Gene set variation analysis (GSVA) and ssGSEA were conducted to determine the pathways that are significantly enriched in different subgroups. Next, we systematically correlated the HLA classifier with immunological characteristics from multiple perspectives, such as immune-related cell infiltration, pivotal molecular pathways, and cytokine expression. Finally, quantitative real-time polymerase chain reaction (qRT-PCR) was conducted to validate the expression level of HLA genes in clinical samples.Results: A total of nine datasets comprising 1,251 patients were included. Based on RF and modified Lasso penalized regression in multi-transcriptome datasets, five HLA genes (B2M, HLA-DQA1, HLA-DPA1, TAP1, and TAP2) were identified as hub genes, which were used to construct an HLA classifier. In the discovery cohort, the HLA classifier exhibited superior diagnostic value (AUC = 0.997) and performed better in predicting mortality (AUC = 0.716) than clinical characteristics or endotypes. Encouragingly, similar results were observed in the ArrayExpress databases. In the E-MTAB-7581 dataset, the use of hydrocortisone in the HLA high-risk subgroup (OR: 2.84, 95% CI 1.07–7.57, p = 0.037) was associated with increased risk of mortality, but not in the HLA low-risk subgroup. Additionally, immune infiltration analysis by CIBERSORTx and ssGSEA revealed that B cells, activated dendritic cells, NK cells, T helper cells, and infiltrating lymphocytes (ILs) were significantly richer in HLA low-risk phenotypes, while Tregs and myeloid-derived suppressor cells (MDSCs) were more abundant in HLA high-risk phenotypes. The HLA classifier was significantly negatively correlated with B cells, activated dendritic cells, NK cells, T helper cells, and ILs, yet was significantly positively correlated with Tregs and MDSCs. Subsequently, molecular pathways analysis uncovered that cytokine-cytokine receptor (CCR) interaction, human leukocyte antigen (HLA), and antigen-presenting cell (APC) co-stimulation were significantly enriched in HLA low-risk endotypes, which was significantly negatively correlated with the HLA classifier in multi-transcriptome data. Finally, the expression levels of several cytokines (IL-10, IFNG, TNF) were significantly different between the HLA subgroups, and the ratio of IL-10/TNF was significantly positively correlated with HLA score in multi-transcriptome data. Results of qRT-PCR validated the higher expression level of B2M as well as lower expression level of HLA-DQA1, HLA-DPA1, TAP1, and TAP2 in sepsis samples compared to control sample.Conclusion: Based on five HLA genes, a diagnostic and prognostic model, namely the HLA classifier, was established, which is closely correlated with responses to hydrocortisone and immunosuppression status and might facilitate personalized counseling for specific therapy.https://www.frontiersin.org/articles/10.3389/fphys.2022.870657/fullsepsisHLA genesimmune infiltrationimmunosuppressionmodel
spellingShingle Zhen Chen
Rui Chen
Yangpeng Ou
Jianhai Lu
Qianhua Jiang
Genglong Liu
Liping Wang
Yayun Liu
Zhujiang Zhou
Ben Yang
Liuer Zuo
Construction of an HLA Classifier for Early Diagnosis, Prognosis, and Recognition of Immunosuppression in Sepsis by Multiple Transcriptome Datasets
Frontiers in Physiology
sepsis
HLA genes
immune infiltration
immunosuppression
model
title Construction of an HLA Classifier for Early Diagnosis, Prognosis, and Recognition of Immunosuppression in Sepsis by Multiple Transcriptome Datasets
title_full Construction of an HLA Classifier for Early Diagnosis, Prognosis, and Recognition of Immunosuppression in Sepsis by Multiple Transcriptome Datasets
title_fullStr Construction of an HLA Classifier for Early Diagnosis, Prognosis, and Recognition of Immunosuppression in Sepsis by Multiple Transcriptome Datasets
title_full_unstemmed Construction of an HLA Classifier for Early Diagnosis, Prognosis, and Recognition of Immunosuppression in Sepsis by Multiple Transcriptome Datasets
title_short Construction of an HLA Classifier for Early Diagnosis, Prognosis, and Recognition of Immunosuppression in Sepsis by Multiple Transcriptome Datasets
title_sort construction of an hla classifier for early diagnosis prognosis and recognition of immunosuppression in sepsis by multiple transcriptome datasets
topic sepsis
HLA genes
immune infiltration
immunosuppression
model
url https://www.frontiersin.org/articles/10.3389/fphys.2022.870657/full
work_keys_str_mv AT zhenchen constructionofanhlaclassifierforearlydiagnosisprognosisandrecognitionofimmunosuppressioninsepsisbymultipletranscriptomedatasets
AT ruichen constructionofanhlaclassifierforearlydiagnosisprognosisandrecognitionofimmunosuppressioninsepsisbymultipletranscriptomedatasets
AT yangpengou constructionofanhlaclassifierforearlydiagnosisprognosisandrecognitionofimmunosuppressioninsepsisbymultipletranscriptomedatasets
AT jianhailu constructionofanhlaclassifierforearlydiagnosisprognosisandrecognitionofimmunosuppressioninsepsisbymultipletranscriptomedatasets
AT qianhuajiang constructionofanhlaclassifierforearlydiagnosisprognosisandrecognitionofimmunosuppressioninsepsisbymultipletranscriptomedatasets
AT genglongliu constructionofanhlaclassifierforearlydiagnosisprognosisandrecognitionofimmunosuppressioninsepsisbymultipletranscriptomedatasets
AT lipingwang constructionofanhlaclassifierforearlydiagnosisprognosisandrecognitionofimmunosuppressioninsepsisbymultipletranscriptomedatasets
AT yayunliu constructionofanhlaclassifierforearlydiagnosisprognosisandrecognitionofimmunosuppressioninsepsisbymultipletranscriptomedatasets
AT zhujiangzhou constructionofanhlaclassifierforearlydiagnosisprognosisandrecognitionofimmunosuppressioninsepsisbymultipletranscriptomedatasets
AT benyang constructionofanhlaclassifierforearlydiagnosisprognosisandrecognitionofimmunosuppressioninsepsisbymultipletranscriptomedatasets
AT liuerzuo constructionofanhlaclassifierforearlydiagnosisprognosisandrecognitionofimmunosuppressioninsepsisbymultipletranscriptomedatasets