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.     

牛痘病毒感染人巨噬细胞基因表达谱的RNA-Seq分析

[目的] 利用RNA-Seq,分析人巨噬细胞在牛痘病毒（VACV）感染前后基因表达的变化,探索牛痘病毒与宿主细胞相互作用的机制。[方法] 用牛痘病毒感染人巨噬细胞,采用RNA-Seq比较感染组和对照组的差异表达基因,并进行KEGG、GO以及STRING网络分析。[结果] 感染组与对照组相比,筛选出显著性差异表达基因4796个,其中上调表达2416个,下调表达2380个。KEGG富集分析表明差异基因主要参与新陈代谢、信号转导、免疫系统和感染疾病等通路。GO功能注释显示这些基因富集在细胞功能调控、物质代谢和免疫调控等生命过程。运用STRING在线蛋白互作数据库,对NOD样信号通路进行分析,鉴定出以JUN、CHUK、IL1B和PYCARD 为核心的调控基因。[结论] 牛痘病毒感染可以诱导人巨噬细胞基因差异性表达,涉及的生物学过程有很多,对免疫相关的信号通路进行深入的分析,发现C型凝集素受体信号通路（C-type lectin receptor signaling pathway）、Toll样受体信号通路以及NOD样通路等多条通路可能参与调控牛痘病毒感染引起的炎症反应,该研究为解析牛痘病毒的感染机制和免疫逃逸机制以及其在临床上治疗感染性疾病和癌症提供了新思路。     

三棱内酯B调控人冠状动脉内皮细胞基因表达谱的生物信息学分析

目的:分析三棱内酯B在人冠状动脉内皮细胞中的表达谱数据集,寻找三棱内酯B调控血管内皮功能的关键作用靶点。方法:基于GEO公共数据库,下载原始表达谱数据集(GSE44598),经过差异基因筛选,功能注释,通路富集,信号通路网络以及基因互作网络分析,找出三棱内酯B对人冠状动脉内皮细胞基因表达谱产生影响的关键基因和信号通路。结果:同对照组相比,三棱内酯B给药组共有5224个基因有显著性差异,包括2628个上调基因和2596个下调基因。基因功能注释和信号通路富集分析表明,差异基因主要参与了细胞周期过程。网络分析显示,MAPK信号通路、细胞周期通路以及PLCG2,PRKACA和ADCY4等为关键信号通路和基因。结论:三棱内酯B通过影响PLCG2,PRKACA,ADCY4等基因的表达,参与MAPK和细胞周期等信号通路,从而调节人冠状内皮细胞的功能。这些关键基因和信号通路是三棱内酯B在心血管疾病治疗应用中潜在的作用靶点。     

All-trans retinoic acid (ATRA) induces miR-23a expression,decreases CTSC expression and granzyme B activity leading to impaired NK cell cytotoxicity

NK cell is an innate immune system lymphocyte lineage with natural cytotoxicity. Its optimal use in the clinic requires in vitro expansion and activation. Cytokines and encounter with target cells activate NK cells and induce proliferation, and this could depend on the presence of other immune cells. Here we activated PBMCs during 5 days with IL-2, with IL-2 plus the tumor cell line K562 and with the lymphoblastoid cell line R69 and perform integrated analyses of microRNA and mRNA expression profiles of purified NK cells. The samples cluster depending on the stimuli and not on the donor, indicating that the pattern of NK cell stimulation is acutely well conserved between individuals. Regulation of mRNA expression is tighter than that of miRNA expression. All stimuli induce a common preserved genetic remodeling. In addition, encounter with target cells mainly activates pathways related to metabolism. Different target cells induce different NK cell remodeling which affects cytokine response and cytotoxicity, supporting the notion that encounter with different target cells significantly changing the activation pattern. We validate our analysis by showing that activation down regulates miR-23a, which is a negative regulator of cathepsin C (CTSC) mRNA, a gene up regulated by all stimuli. The peptidase CTSC activates the granzymes, the main effector proteases involved in NK cell cytotoxicity. All-trans retinoic acid (ATRA), which induces miR-23a expression, decreases CTSC expression and granzyme B activity leading to impaired NK cell cytotoxicity in an in vivo mouse model.     

BACL Is a Novel Brain-Associated,Non-NKC-Encoded Mammalian C-Type Lectin-Like Receptor of the CLEC2 Family

Natural Killer Gene Complex (NKC)–encoded C-type lectin-like receptors (CTLRs) are expressed on various immune cells including T cells, NK cells and myeloid cells and thereby contribute to the orchestration of cellular immune responses. Some NKC-encoded CTLRs are grouped into the C-type lectin family 2 (CLEC2 family) and interact with genetically linked CTLRs of the NKRP1 family. While many CLEC2 family members are expressed by hematopoietic cells (e.g. CD69 (CLEC2C)), others such as the keratinocyte-associated KACL (CLEC2A) are specifically expressed by other tissues. Here we provide the first characterization of the orphan gene CLEC2L. In contrast to other CLEC2 family members, CLEC2L is conserved among mammals and located outside of the NKC. We show that CLEC2L-encoded CTLRs are expressed as non-glycosylated, disulfide-linked homodimers at the cell surface. CLEC2L expression is fairly tissue-restricted with a predominant expression in the brain. Thus CLEC2L-encoded CTLRs were designated BACL (brain-associated C-type lectin). Combining in situ hybridization and immunohistochemistry, we show that BACL is expressed by neurons in the CNS, with a pronounced expression by Purkinje cells. Notably, the CLEC2L locus is adjacent to another orphan CTLR gene (KLRG2), but reporter cell assays did neither indicate interaction of BACL with the KLRG2 ectodomain nor with human NK cell lines or lymphocytes. Along these lines, growth of BACL-expressing tumor cell lines in immunocompetent mice did not provide evidence for an immune-related function of BACL. Altogether, the CLEC2L gene encodes a homodimeric cell surface CTLR that stands out among CLEC2 family members by its conservation in mammals, its biochemical properties and the predominant expression in the brain. Future studies will have to reveal insights into the functional relevance of BACL in the context of its neuronal expression.     

Mechanisms of NK Cell-Macrophage Bacillus anthracis Crosstalk: A Balance between Stimulation by Spores and Differential Disruption by Toxins

NK cells are important immune effectors for preventing microbial invasion and dissemination, through natural cytotoxicity and cytokine secretion. Bacillus anthracis spores can efficiently drive IFN-γ production by NK cells. The present study provides insights into the mechanisms of cytokine and cellular signaling that underlie the process of NK-cell activation by B. anthracis and the bacterial strategies to subvert and evade this response. Infection with non-toxigenic encapsulated B. anthracis induced recruitment of NK cells and macrophages into the mouse draining lymph node. Production of edema (ET) or lethal (LT) toxin during infection impaired this cellular recruitment. NK cell depletion led to accelerated systemic bacterial dissemination. IFN-γ production by NK cells in response to B. anthracis spores was: i) contact-dependent through RAE-1-NKG2D interaction with macrophages; ii) IL-12, IL-18, and IL-15-dependent, where IL-12 played a key role and regulated both NK cell and macrophage activation; and iii) required IL-18 for only an initial short time window. B. anthracis toxins subverted both NK cell essential functions. ET and LT disrupted IFN-γ production through different mechanisms. LT acted both on macrophages and NK cells, whereas ET mainly affected macrophages and did not alter NK cell capacity of IFN-γ secretion. In contrast, ET and LT inhibited the natural cytotoxicity function of NK cells, both in vitro and in vivo. The subverting action of ET thus led to dissociation in NK cell function and blocked natural cytotoxicity without affecting IFN-γ secretion. The high efficiency of this process stresses the impact that this toxin may exert in anthrax pathogenesis, and highlights a potential usefulness for controlling excessive cytotoxic responses in immunopathological diseases. Our findings therefore exemplify the delicate balance between bacterial stimulation and evasion strategies. This highlights the potential implication of the crosstalk between host innate defences and B. anthracis in initial anthrax control mechanisms.     

Correlation of zinc finger protein 2, a prognostic biomarker,with immune infiltrates in liver cancer

Purpose: The expression and clinical value of zinc finger protein 2 gene (ZIC2) in hepatocellular carcinoma (HCC) were analyzed by mining gene information from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases.Methods: Gene chip data sets were retrieved from GEO and TCGA and screened for differentially expressed genes in HCC. Gene expression profile interaction analysis (GEPIA) and Kaplan–Meier curves were used to analyze the relationship between differentially expressed genes (DEGs) and survival and prognosis in patients with HCC. Moreover, the Genecards database was used to extract ZIC2-related proteins and to analyze the physiological process of protein enrichment. Furthermore, the relationships between ZIC2 gene and tumor cell immune invasion and that between immune cell infiltration and the 5-year survival rate were studied using the tumor immune evaluation resource (TIMER) database.Results: Datasets from GEO and TCGA revealed that ZIC2 was differentially expressed in HCC tissues and normal tissues (P<0.05). High ZIC2 expression was associated with overall survival (OS) and progress-free survival in HCC patients. Overall, 25 ZIC2 related proteins, including Gli3, PRKDC, and rnf180 were identified and protein enrichment analysis indicated these were associated with four types of cell components, six types of cell functions, and eight types of biological processes. ZIC2 was positively correlated with immune infiltration cells in patients with HCC, and higher expression of ZIC2 mRNA CD4+T cells is associated with a better 5-year survival.Conclusion: ZIC2 gene may be used as an immune response marker in liver cancer to predict the prognosis of HCC.     

Systematic analysis on multiple Gene Expression Omnibus data sets reveals fierce immune response in hepatitis B virus‐related acute liver failure

Acute liver failure (ALF) caused by hepatitis B virus (HBV) is common type of liver failure in the world, with high morbidity and mortality rates. However, the prevalence, genetic background and factors determining the development of HBV‐related ALF are rarely studied. In this study, we examined three Gene Expression Omnibus (GEO) data sets by bioinformatics analysis to identify differentially expressed genes (DEGs), key biological processes and pathways. Immune infiltration analysis showed high immune cells infiltration in HBV‐related ALF tissue. We then confirmed natural killer cells and macrophages infiltration in clinical samples by immunohistochemistry assay, implying these cells play a significant role in HBV‐ALF. We found 1277 genes were co‐up‐regulated and that 1082 genes were co‐down‐regulated in the 3 data sets. Inflammation‐related pathways were enriched in the co‐up‐regulated genes and synthetic metabolic pathways were enriched in the co‐down‐regulated genes. WGCNA also revealed a key module enriching in immune inflammation response and identified 10 hub genes, differentially expressed in an independent data set. In conclusion, we identified fierce immune inflammatory response to elucidate the immune‐driven mechanism of HBV‐ALF and 10 hub genes based on gene expression profiles.     

Gene Expression Profiling of the Response to Interferon Beta in Epstein-Barr-Transformed and Primary B Cells of Patients with Multiple Sclerosis

The effects of interferon-beta (IFN-β), one of the key immunotherapies used in multiple sclerosis (MS), on peripheral blood leukocytes and T cells have been extensively studied. B cells are a less abundant leukocyte type, and accordingly less is known about the B cell-specific response to IFN-β. To identify gene expression changes and pathways induced by IFN-β in B cells, we studied the in vitro response of human Epstein Barr-transformed B cells (lymphoblast cell lines-LCLs), and validated our results in primary B cells. LCLs were derived from an MS patient repository. Whole genome expression analysis identified 115 genes that were more than two-fold differentially up-regulated following IFN-β exposure, with over 50 previously unrecognized as IFN-β response genes. Pathways analysis demonstrated that IFN-β affected LCLs in a similar manner to other cell types by activating known IFN-β canonical pathways. Additionally, IFN-β increased the expression of innate immune response genes, while down-regulating many B cell receptor pathway genes and genes involved in adaptive immune responses. Novel response genes identified herein, NEXN, DDX60L, IGFBP4, and HAPLN3, B cell receptor pathway genes, CD79B and SYK, and lymphocyte activation genes, LAG3 and IL27RA, were validated as IFN-β response genes in primary B cells. In this study new IFN-β response genes were identified in B cells, with possible implications to B cell-specific functions. The study''s results emphasize the applicability of LCLs for studies of human B cell drug response. The usage of LCLs from patient-based repositories may facilitate future studies of drug response in MS and other immune-mediated disorders with a B cell component.     

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.     

Single-cell RNA sequencing reveals the sustained immune cell dysfunction in the pathogenesis of sepsis secondary to bacterial pneumonia

Sepsis is a leading cause of mortality in intensive care unit worldwide, it's accompanied by immune cell dysfunction induced by multiple factors. However, little is known about the specific alterations in immune cells in the dynamic pathogenesis of sepsis secondary to bacterial pneumonia. Here, we used single cell RNA sequencing (scRNA-seq) to profile peripheral blood mononuclear cells (PBMCs) in a healthy control and two patients with sepsis secondary to bacterial pneumonia, including acute, stable and recovery stage. We analyzed the quantity and function of immune cells. During disease course, interferon gamma response was upregulated; T/NK cell subtypes presented activation and exhaustion properties, which might be driven by monocytes through IL-1β signaling pathways; The proportion of plasma cells was increased, which might be driven by NK cells through IFN signaling pathways; Additionally, interferon gamma response was upregulated to a greater degree in sepsis secondary to pneumonia induced by SARS-COV-2 compared with that induced by influenza virus and bacteria.     

Immunogenomics Analysis Reveals that TP53 Mutations Inhibit Tumor Immunity in Gastric Cancer

Although immunotherapy continues to demonstrate efficacy in a variety of refractory cancers, currently, no any immunotherapeutic strategy is clinically used for gastric cancer (GC) except its microsatellite instable subtype. Thus, it is important to identify molecular biomarkers for predicting the responders to GC immunotherapy. TP53 mutations frequently occur in GC and are associated with unfavorable clinical outcomes in GC. We performed a comprehensive characterization of the associations between TP53 mutations and immune activities in GC based on two large-scale GC cancer genomics data. We compared expression and enrichment levels of 787 immune-related genes and 23 immune gene-sets among TP53-mutated GCs, TP53‐wildtype GCs, and normal tissue, and explored the correlations between p53-mediated pathways and immune activities in GC. Strikingly, almost all analyzed immune gene-sets were significantly downregulated in enrichment levels in TP53-mutated GCs compared to TP53‐wildtype GCs. These less active immune pathways and cell types in TP53-mutated GCs included 15 immune cell types and function, tumor-infiltrating lymphocytes, regulatory T cells, immune checkpoint, cytokine and cytokine receptor, human leukocyte antigen, pro‐inflammatory, and parainflammation. Moreover, we identified a number of p53-mediated pathways and proteins that were significantly associated with immune activities in GC. Furthermore, we demonstrated that the TP53 mutation itself could result in the depressed immune activities in GC and other cancer types. We revealed that chromosomal instability was an important mechanism for the depressed tumor immunity in TP53-mutated cancers. Finally, we showed that immune cell infiltration and immune activities were likely positively associated with survival prognosis in GC. Our findings suggest that p53 may play an important role in activating tumor immunity in GC and other cancer types and that the TP53 mutation status could be useful in stratifying cancer patients responsive to a certain immunotherapy.     

Identifying significant genes and functionally enriched pathways in familial hypercholesterolemia using integrated gene co-expression network analysis

Familial hypercholesterolemia (FH) is a monogenic lipid disorder which promotes atherosclerosis and cardiovascular diseases. Owing to the lack of sufficient published information, this study aims to identify the potential genetic biomarkers for FH by studying the global gene expression profile of blood cells. The microarray expression data of FH patients and controls was analyzed by different computational biology methods like differential expression analysis, protein network mapping, hub gene identification, functional enrichment of biological pathways, and immune cell restriction analysis. Our results showed the dysregulated expression of 115 genes connected to lipid homeostasis, immune responses, cell adhesion molecules, canonical Wnt signaling, mucin type O-glycan biosynthesis pathways in FH patients. The findings from expanded protein interaction network construction with known FH genes and subsequent Gene Ontology (GO) annotations have also supported the above findings, in addition to identifying the involvement of dysregulated thyroid hormone and ErbB signaling pathways in FH patients. The genes like CSNK1A1, JAK3, PLCG2, RALA, and ZEB2 were found to be enriched under all GO annotation categories. The subsequent phenotype ontology results have revealed JAK3I, PLCG2, and ZEB2 as key hub genes contributing to the inflammation underlying cardiovascular and immune response related phenotypes. Immune cell restriction findings show that above three genes are highly expressed by T-follicular helper CD4⁺ T cells, naïve B cells, and monocytes, respectively. These findings not only provide a theoretical basis to understand the role of immune dysregulations underlying the atherosclerosis among FH patients but may also pave the way to develop genomic medicine for cardiovascular diseases.     

NK cell activation and tumor infiltration are involved in the antitumor mechanism of Virulizin

Previous studies have demonstrated antitumor efficacy of Virulizin in several human tumor xenograft models and a critical role for macrophages in the antitumor mechanism of Virulizin. Although there is growing support for an immune stimulatory mechanism of action for Virulizin, the details remain to be elucidated. The aim of this study was to determine whether infiltration of natural killer (NK) cells into xenografted tumors is altered by Virulizin treatment, and whether such alterations contribute to the antitumor activity of Virulizin. Immunohistochemical analysis demonstrated that xenografted tumors from Virulizin-treated mice had an increase in infiltration of F4/80⁺ (macrophages) and NK1.1⁺ (NK) cells. The increase in NK1.1⁺ cell infiltration occurred at an early stage of Virulizin treatment, which correlated with an early sign of apoptosis. In addition, Virulizin resulted in an increase in the number of NK cells in the spleens, and NK cells isolated from the spleen exhibited increased cytotoxicity to tumor cells in vitro. In NK cell–deficient SCID-beige mice, the antitumor activity of Virulizin was compromised, providing additional support to the hypothesis that NK cells are necessary for inhibition of tumor growth by Virulizin. Finally, depletion of macrophages resulted in the loss of Virulizin-induced increase in NK1.1⁺ cell infiltration into xenografted tumors, suggesting the involvement of macrophages in NK cell infiltration into tumors. Taken together, these results strongly support a mechanism in which Virulizin stimulates a sustained expansion and infiltration of NK cells and macrophages into tumors with subsequent activation of NK cells that is responsible for the observed antitumor activity.