Screening and identification of the core immune-related genes and immune cell infiltration in severe burns and sepsis

Severe burns often have a high mortality rate due to sepsis, but the genetic and immune crosstalk between them remains unclear. In the present study, the GSE77791 and GSE95233 datasets were analysed to identify immune-related differentially expressed genes (DEGs) involved in disease progression in both burns and sepsis. Subsequently, weighted gene coexpression network analysis (WGCNA), gene enrichment analysis, protein–protein interaction (PPI) network construction, immune cell infiltration analysis, core gene identification, coexpression network analysis and clinical correlation analysis were performed. A total of 282 common DEGs associated with burns and sepsis were identified. Kyoto Encyclopedia of Genes and Genomes pathway analysis identified the following enriched pathways in burns and sepsis: metabolic pathways; complement and coagulation cascades; legionellosis; starch and sucrose metabolism; and ferroptosis. Finally, six core DEGs were identified, namely, IL10, RETN, THBS1, FGF13, LCN2 and MMP9. Correlation analysis showed that some core DEGs were significantly associated with simultaneous dysregulation of immune cells. Of these, RETN upregulation was associated with a worse prognosis. The immune-related genes and dysregulated immune cells in severe burns and sepsis provide potential research directions for diagnosis and treatment.     

肾透明细胞癌患者铁死亡相关基因的预后作用

探讨铁死亡相关基因在肾透明细胞癌患者中的表达及其预后价值。通过TCGA数据库下载KIRC的相关测序数据与检索到的铁死亡相关基因取交集,进行铁死亡相关基因的差异分析。之后利用单变量和多变量Cox回归分析,筛选具有预后价值的基因,构建预测患者生存情况的风险评分模型,并对模型进行验证。对高低风险组进行GO与KEGG通路富集,探讨风险差异的可能原因;通过ssGSEA分析,评估高低风险组间的免疫浸润情况。在KIRC患者的肿瘤组织和正常组织中,共得到21个差异的铁死亡相关基因;通过单因素Cox回归分析,获得 28 个与KIRC预后相关的基因;之后进行Lasso回归与多因素Cox回归分析,结果显示有10个基因被纳入模型,计算公式为:风险值(Risk score)=(0.024 5)×ALOX5表达值+(0.126 0)×CBS表达值+(0.199 5)×CD44表达值+(0.218 3)×CHAC1表达值+(-0.295 9)×HMGCR表达值+(0.036 7)×MT1G表达值+(0.061 4)×SLC7A11表达值+(-0.080 7)×FDFT1表达值+(0.160 3)×PEBP1表达值+(-0.220 5)×GOT1表达值。生存状态图表明,高风险组死亡病例数多于低风险组;ROC曲线表明风险评分模型具备一定预测能力;K-M生存分析显示,高风险组总体生存率低于低风险组(P=5.73×10^-13)。GO与KEGG富集分析提示,高低风险组间免疫情况及IL-17信号通路存在显著差异;进一步的ssGSEA富集显示,高低风险组间大部分免疫细胞的评分存在显著差异。基于铁死亡相关基因的预后风险评分模型可用于KIRC的预后预测,针对铁死亡相关基因设计靶点可能是治疗KIRC的一种新选择。     

Prognostic value of tumour microenvironment-related genes by TCGA database in rectal cancer

Rectal cancer is a common malignant tumour and the progression is highly affected by the tumour microenvironment (TME). This study intended to assess the relationship between TME and prognosis, and explore prognostic genes of rectal cancer. The gene expression profile of rectal cancer was obtained from TCGA and immune/stromal scores were calculated by Estimation of Stromal and Immune cells in Malignant Tumors using Expression data (ESTIMATE) algorithm. The correlation between immune/stromal scores and survival time as well as clinical characteristics were evaluated. Differentially expressed genes (DEGs) were identified according to the stromal/immune scores, and the functional enrichment analyses were conducted to explore functions and pathways of DEGs. The survival analyses were conducted to clarify the DEGs with prognostic value, and the protein-protein interaction (PPI) network was performed to explore the interrelation of prognostic DEGs. Finally, we validated prognostic DEGs using data from the Gene Expression Omnibus (GEO) database by PrognoScan, and we verified these genes at the protein levels using the Human Protein Atlas (HPA) databases. We downloaded gene expression profiles of 83 rectal cancer patients from The Cancer Genome Atlas (TCGA) database. The Kaplan-Meier plot demonstrated that low-immune score was associated with worse clinical outcome (P = .034), metastasis (M1 vs. M0, P = .031) and lymphatic invasion (+ vs. -, P < .001). A total of 540 genes were screened as DEGs with 539 up-regulated genes and 1 down-regulated gene. In addition, 60 DEGs were identified associated with overall survival. Functional enrichment analyses and PPI networks showed that the DEGs are mainly participated in immune process, and cytokine-cytokine receptor interaction. Finally, 19 prognostic genes were verified by GSE17536 and GSE17537 from GEO, and five genes (ADAM23, ARHGAP20, ICOS, IRF4, MMRN1) were significantly different in tumour tissues compared with normal tissues at the protein level. In summary, our study demonstrated the associations between TME and prognosis as well as clinical characteristics of rectal cancer. Moreover, we explored and verified microenvironment-related genes, which may be the potential key prognostic genes of rectal cancer. Further clinical samples and functional studies are needed to validate this finding.     

Systematic analysis of alternative splicing signature unveils prognostic predictor for kidney renal clear cell carcinoma

There is growing evidence that alternative splicing (AS) plays an important role in cancer development. However, a comprehensive analysis of AS signatures in kidney renal clear cell carcinoma (KIRC) is lacking and urgently needed. It remains unclear whether AS acts as diagnostic biomarkers in predicting the prognosis of KIRC patients. In the work, gene expression and clinical data of KIRC were obtained from The Cancer Genome Atlas (TCGA), and profiles of AS events were downloaded from the SpliceSeq database. The RNA sequence/AS data and clinical information were integrated, and we conducted the Cox regression analysis to screen survival-related AS events and messenger RNAs (mRNAs). Correlation between prognostic AS events and gene expression were analyzed using the Pearson correlation coefficient. Protein-protein interaction analysis was conducted for the prognostic AS-related genes, and a potential regulatory network was built using Cytoscape (version 3.6.1). Meanwhile, functional enrichment analysis was conducted. A prognostic risk score model is then established based on seven hub genes (KRT222, LENG8, APOB, SLC3A1, SCD5, AQP1, and ADRA1A) that have high performance in the risk classification of KIRC patients. A total 46,415 AS events including 10,601 genes in 537 patients with KIRC were identified. In univariate Cox regression analysis, 13,362 survival associated AS events and 8,694 survival-specific mRNAs were detected. Common 3,105 genes were screen by overlapping 13,362 survival associated AS events and 8,694 survival-specific mRNAs. The Pearson correlation analysis suggested that 13 genes were significantly correlated with AS events (Pearson correlation coefficient >0.8 or <−0.8). Then, We conducted multivariate Cox regression analyses to select the potential prognostic AS genes. Seven genes were identified to be significantly related to OS. A prognostic model based on seven genes was constructed. The area under the ROC curve was 0.767. In the current study, a robust prognostic prediction model was constructed for KIRC patients, and the findings revealed that the AS events could act as potential prognostic biomarkers for KIRC.     

Discovery of core gene families associated with liver metastasis in colorectal cancer and regulatory roles in tumor cell immune infiltration

In this study, we aimed to uncover genes that drive the pathogenesis of liver metastasis in colorectal cancer (CRC), and identify effective genes that could serve as potential therapeutic targets for treating with colorectal liver metastasis patients based on two GEO datasets. Several bioinformatics approaches were implemented. First, differential expression analysis screened out key differentially expressed genes (DEGs) across the two GEO datasets. Based on gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, we identified the enrichment functions and pathways of the DEGs that were associated with liver metastasis in CRC. Second, immune infiltration analysis identified key immune signature gene sets associated with CRC liver metastasis, among which two key immune gene families (CD and CCL) identified as key DEGs were filtered by protein-protein interaction (PPI) network. Some of the members in these gene families were associated with disease free survival (DFS) or overall survival (OS) in two subtypes of CRC, namely COAD and READ. Finally, functional enrichment analysis of the two gene families and their neighboring genes revealed that they were closely associated with cytokine, leukocyte proliferation and chemotaxis. These results are valuable in comprehending the pathogenesis of liver metastasis in CRC, and are of seminal importance in understanding the role of immune tumor infiltration in CRC. Our study also identified potentially effective therapeutic targets for liver metastasis in CRC including CCL20, CCL24 and CD70.     

Five candidate biomarkers associated with the diagnosis and prognosis of cervical cancer

Purpose: Cervical cancer (CC) is one of the most general gynecological malignancies and is associated with high morbidity and mortality. We aimed to select candidate genes related to the diagnosis and prognosis of CC.Methods: The mRNA expression profile datasets were downloaded. We also downloaded RNA-sequencing gene expression data and related clinical materials from TCGA, which included 307 CC samples and 3 normal samples. Differentially expressed genes (DEGs) were obtained by R software. GO function analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of DEGs were performed in the DAVID dataset. Using machine learning, the optimal diagnostic mRNA biomarkers for CC were identified. We used qRT-PCR and Human Protein Atlas (HPA) database to exhibit the differences in gene and protein levels of candidate genes.Results: A total of 313 DEGs were screened from the microarray expression profile datasets. DNA methyltransferase 1 (DNMT1), Chromatin Assembly Factor 1, subunit B (CHAF1B), Chromatin Assembly Factor 1, subunit A (CHAF1A), MCM2, CDKN2A were identified as optimal diagnostic mRNA biomarkers for CC. Additionally, the GEPIA database showed that the DNMT1, CHAF1B, CHAF1A, MCM2 and CDKN2A were associated with the poor survival of CC patients. HPA database and qRT-PCR confirmed that these genes were highly expressed in CC tissues.Conclusion: The present study identified five DEmRNAs, including DNMT1, CHAF1B, CHAF1A, MCM2 and Kinetochore-related protein 1 (KNTC1), as potential diagnostic and prognostic biomarkers of CC.     

Identification of prognostic biomarkers for breast cancer brain metastases based on the bioinformatics analysis

PurposeThe prognosis of breast cancer (BC) patients who develop into brain metastases (BMs) is very poor. Thus, it is of great significance to explore the etiology of BMs in BC and identify the key genes involved in this process to improve the survival of BC patients with BMs.Patients and methodsThe gene expression data and the clinical information of BC patients were downloaded from TCGA and GEO database. Differentially expressed genes (DEGs) in TCGA-BRCA and GSE12276 were overlapped to find differentially expressed metastatic genes (DEMGs). The protein-protein interaction (PPI) network of DEMGs was constructed via STRING database. ClusterProfiler R package was applied to perform the gene ontology (GO) enrichment analysis of DEMGs. The univariate Cox regression analysis and the Kaplan-Meier (K-M) curves were plotted to screen DEMGs associated with the overall survival and the metastatic recurrence survival, which were identified as the key genes associated with the BMs in BC. The immune infiltration and the expressions of immune checkpoints for BC patients with brain relapses and BC patients with other relapses were analyzed respectively. The correlations among the expressions of key genes and the differently infiltrated immune cells or the differentially expressed immune checkpoints were calculated. The gene set enrichment analysis (GSEA) of each key gene was conducted to investigate the potential mechanisms of key genes involved in BC patients with BMs. Moreover, CTD database was used to predict the drug-gene interaction network of key genes.ResultsA total of 154 DEGs were identified in BC patients at M0 and M1 in TCGA database. A total of 667 DEGs were identified in BC patients with brain relapses and with other relapses. By overlapping these DEGs, 17 DEMGs were identified, which were enriched in the cell proliferation related biological processes and the immune related molecular functions. The univariate Cox regression analysis and the Kaplan-Meier curves revealed that CXCL9 and GPR171 were closely associated with the overall survival and the metastatic recurrence survival and were identified as key genes associated with BMs in BC. The analyses of immune infiltration and immune checkpoint expressions showed that there was a significant difference of the immune microenvironment between brain relapses and other relapses in BC. GSEA indicated that CXCL9 and GPR171 may regulate BMs in BC via the immune-related pathways.ConclusionOur study identified the key genes associated with BMs in BC patients and explore the underlying mechanisms involved in the etiology of BMs in BC. These findings may provide a promising approach for the treatments of BC patients with BMs.     

Identification of potential core genes and pathways predicting pathogenesis in head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) is the most common subtype of head and neck cancer; however, its pathogenesis and potential therapeutic targets remain largely unknown. In the present study, we analyzed three gene expression profiles and screened differentially expressed genes (DEGs) between HNSCC and normal tissues. The DEGs were subjected to gene ontology (GO), Kyoto encyclopedia of genes and genomes (KEGG), protein–protein interaction (PPI), and survival analyses, while the connectivity map (CMap) database was used to predict candidate small molecules that may reverse the biological state of HNSCC. Finally, we measured the expression of the most relevant core gene in vitro and examined the effect of the top predicted potential drug against the proliferation of HNSCC cell lines. Among the 208 DEGs and ten hub genes identified, CDK1 and CDC45 were associated with unfavorable HNSCC prognosis, and three potential small molecule drugs for treating HNSCC were identified. Increased CDK1 expression was confirmed in HNSCC cells, and menadione, the top predicted potential drug, exerted significant inhibitory effects against HNSCC cell proliferation and markedly reversed CDK1 expression. Together, the findings of the present study suggest that the ten hub genes and pathways identified may be closely related to HNSCC pathogenesis. In particular, CDK1 and CDC45 overexpression could be reliable biomarkers for predicting unfavorable prognosis in patients with HNSCC, while the new candidate small molecules identified by CMap analysis provide new avenues for the development of potential drugs to treat HNSCC.     

Systematic profiling of invasion-related gene signature predicts prognostic features of lung adenocarcinoma

Due to the high heterogeneity of lung adenocarcinoma (LUAD), molecular subtype based on gene expression profiles is of great significance for diagnosis and prognosis prediction in patients with LUAD. Invasion-related genes were obtained from the CancerSEA database, and LUAD expression profiles were downloaded from The Cancer Genome Atlas. The ConsensusClusterPlus was used to obtain molecular subtypes based on invasion-related genes. The limma software package was used to identify differentially expressed genes (DEGs). A multi-gene risk model was constructed by Lasso-Cox analysis. A nomogram was also constructed based on risk scores and meaningful clinical features. 3 subtypes (C1, C2 and C3) based on the expression of 97 invasion-related genes were obtained. C3 had the worst prognosis. A total of 669 DEGs were identified among the subtypes. Pathway enrichment analysis results showed that the DEGs were mainly enriched in the cell cycle, DNA replication, the p53 signalling pathway and other tumour-related pathways. A 5-gene signature (KRT6A, MELTF, IRX5, MS4A1 and CRTAC1) was identified by using Lasso-Cox analysis. The training, validation and external independent cohorts proved that the model was robust and had better prediction ability than other lung cancer models. The gene expression results showed that the expression levels of MS4A1 and KRT6A in tumour tissues were higher than in normal tissues, while CRTAC1 expression in tumour tissues was lower than in normal tissues. The 5-gene signature prognostic stratification system based on invasion-related genes could be used to assess prognostic risk in patients with LUAD.     

Identification of key candidate genes and pathways in multiple myeloma by integrated bioinformatics analysis

Multiple myeloma (MM) is a common hematologic malignancy for which the underlying molecular mechanisms remain largely unclear. This study aimed to elucidate key candidate genes and pathways in MM by integrated bioinformatics analysis. Expression profiles GSE6477 and GSE47552 were obtained from the Gene Expression Omnibus database, and differentially expressed genes (DEGs) with p < .05 and [logFC] > 1 were identified. Functional enrichment, protein–protein interaction network construction and survival analyses were then performed. First, 51 upregulated and 78 downregulated DEGs shared between the two GSE datasets were identified. Second, functional enrichment analysis showed that these DEGs are mainly involved in the B cell receptor signaling pathway, hematopoietic cell lineage, and NF-kappa B pathway. Moreover, interrelation analysis of immune system processes showed enrichment of the downregulated DEGs mainly in B cell differentiation, positive regulation of monocyte chemotaxis and positive regulation of T cell proliferation. Finally, the correlation between DEG expression and survival in MM was evaluated using the PrognoScan database. In conclusion, we identified key candidate genes that affect the outcomes of patients with MM, and these genes might serve as potential therapeutic targets.     

人肝细胞癌预后不良相关基因的生物信息学分析及其临床意义

目的:筛选肝细胞癌(HCC)预后不良相关基因,并探讨其临床意义。方法:在基因表达综合数据库(GEO)中获取符合分析条件的肝细胞癌全基因组表达谱数据并分析得到差异表达基因(DEGs),再运用生物学信息注释及可视化数据库 (DAVID) 和蛋白相互作用数据库 (String) 分别进行功能富集分析和蛋白质互作用网络的构建。利用癌症基因组图谱数据库(TCGA)和Cox比例风险回归模型对相关差异基因进行预后分析。结果:找到一个符合条件的人类HCC数据库 (GSE84402),共筛选出1141个差异表达基因(DEGs),其中上调基因720个,下调基因421个。基因功能富集分析和蛋白质互作用分析结果显示CDK1、CDC6、CCNA2、CHEK1、CENPE 、PIK3R1、RACGAP1、BIRC5、KIF11和CYP2B6为HCC预后的关键基因。TCGA数据库和Cox回归模型分析显示CDC6、PIK3R1、RACGAP1和KIF11的表达升高,CENPE的表达降低与HCC预后不良密切相关。结论:CDC6、CENPE、PIK3R1、RACGAP1和KIF11可能和HCC的预后不良相关,可作为未来HCC预后研究的参考标志物。     

Identification of potential biomarkers associated with immune infiltration in the esophageal carcinoma tumor microenvironment

Tumor immune cell infiltration was significantly correlated with the progression and the effect of immunotherapy in cancers including esophageal carcinoma (ESCA). However, no biomarkers were identified which were associated with immune infiltration in ESCA. In the present study, a total of 128 common differentially expressed genes (DEGs) were identified between esophageal squamous cell carcinomas (ESCC) and esophageal adenocarcinomas (EAC). The results of gene ontology (GO) enrichment and Reactome pathway analysis displayed that the up-regulated DEGs were mainly involved in the regulation of extracellular matrix (ECM), while the down-regulated DEGs were mainly involved in the regulation of cornification and keratinocyte differentiation. The most significant module of up-regulated DEGs was selected by Molecular Complex Detection (MCODE). Top ten similar genes of COL1A2 were explored, then validation and the prognostic analysis of these genes displayed that COL1A2, COL1A1, COL3A1, ZNF469 and Periostin (POSTN) had the prognostic value which were up-regulated in ESCA. The expressions of COL1A2 and its four similar genes were mainly correlated with infiltrating levels of macrophages and dendritic cells (DCs) and showed strong correlations with diverse immune marker sets in ESCA. To summarize, COL1A2 and its four similar genes were identified as the potential biomarkers associated with immune infiltration in ESCA. These genes might be applied to immunotherapy for ESCA.     

Single-cell RNA sequencing reveals transcriptional profiles of monocytes in HBV-infected pregnant women during mid-pregnancy

There is a growing body of evidence that innate immunity also plays an important role in the progression of hepatitis B virus (HBV) infection. However, there is less study on systematically elucidating the characteristics of innate immunity in HBV-infected pregnant women. We compared the features of peripheral blood mononuclear cells in three healthy pregnant women and three HBV-infected pregnant women by single-cell RNA sequencing. 10 DEGs were detected between groups and monocytes were the main expression source of most of the DEGs, which involved in the inflammatory response, apoptosis and immune regulation. Meanwhile, qPCR and ELISA were performed to verify above genes. Monocytes displayed immune response defect, reflecting poor ability of response to IFN. In addition, eight clusters were identified in monocytes. We identified molecular drivers in monocytes subpopulations.TNFSF10+ monocytes, MT1G+ monocytes and TUBB1+ monocytes were featured with different gene expression pattern and biological function.TNFSF10+ monocytes and MT1G+ monocytes were characterized by high levels of inflammation response.TNFSF10+ monocytes, MT1G+ monocytes and TUBB1+ monocytes showed decreased response to IFN. Our results dissects alterations in monocytes related to the immune response of HBV-infected pregnant women and provides a rich resource for fully understanding immunopathogenesis and developing effective preventing HBV intrauterine infection strategies.     

Ferroptosis-related gene CHAC1 is a valid indicator for the poor prognosis of kidney renal clear cell carcinoma

To evaluate the validity of CHAC1 for predicting the prognosis of kidney renal clear cell carcinoma (KIRC) and to explore its therapeutic potential for KIRC, we conducted several bioinformatic analyses using the sequencing data and clinical information derived from online databases. We found CHAC1 is down-regulated in KIRC samples when compared with normal samples but up-regulated in KIRC samples with relatively higher malignancy and later stages. Univariate cox analysis and multivariate cox regression analysis were conducted and the results revealed up-regulated CHAC1 is an independent risk factor for poor prognosis of KIRC. Further, the nomogram model based on the result of multivariate cox regression analysis was constructed and effectively predicted patients' 1-year, 3-year and 5-year survival respectively. The correlation analyses showed CHAC1 is associated with the immune pathway markers of memory B cell, natural killer cell and type1 T helper cell as well as the checkpoint genes like ADORA2A, CD200, CD44, CD70, HHLA2, NRP1, PDCD1LG2 and TNFRSF18. Furthermore, experiments in vitro indicated CHAC1 could induce cell death in KIRC cell lines but had limited influence on cell migration and cell invasion. In conclusion, CHAC1 is found a valid indicator for poor prognosis of kidney renal clear cell carcinoma.