Identification of key genes and novel immune infiltration-associated biomarkers of sepsis

Sepsis is the major cause of mortality in the intensive care unit. The aim of this study was to identify the key prognostic biomarkers of abnormal expression and immune infiltration in sepsis. In this study, a total of 36 differentially expressed genes were identified to be mainly involved in a number of immune-related Gene Ontology terms and Kyoto Encyclopedia of Genes and Genomes pathways. The hub genes (MMP9 and C3AR1) were significantly related to the prognosis of sepsis patients. The immune infiltration analysis indicated a significant difference in the relative cell content of naive B cells, follicular Th cells, activated NK cells, eosinophils, neutrophils and monocytes between sepsis and normal controls. Weighted gene co-expression network analysis and a de-convolution algorithm that quantifies the cellular composition of immune cells were used to analyse the sepsis expression data from the Gene Expression Omnibus database and to identify modules related to differential immune cells. CEBPB is the key immune-related gene that may be involved in sepsis. Gene set enrichment analysis revealed that CEBPB is involved in the processes of T cell selection, B cell–mediated immunity, NK cell activation and pathways of T cells, B cells and NK cells. Therefore, CEBPB may play a key role in the biological and immunological processes of sepsis.     

Effect of dapagliflozin on diabetic patients with cardiovascular disease via MAPK signalling pathway

Clinical studies have shown that dapagliflozin can reduce cardiovascular outcome in patients with type 2 diabetes mellitus (T2DM), but the exact mechanism is unclear. In this study, we used the molecular docking and network pharmacology methods to explore the potential mechanism of dapagliflozin on T2DM complicated with cardiovascular diseases (CVD). Dapagliflozin's potential targets were predicted via the Swiss Target Prediction platform. The pathogenic targets of T2DM and CVD were screened by the Online Mendelian Inheritance in Man (OMIM) and Gene Cards databases. The common targets of dapagliflozin, T2DM and CVD were used to establish a protein-protein interaction (PPI) network; the potential protein functional modules in the PPI network were found out by MCODE. Metascape tool was used for Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway enrichment analysis. A potential protein functional module with the best score was obtained from the PPI network and 9 targets in the protein functional module all showed good binding properties when docking with dapagliflozin. The results of KEGG pathway enrichment analysis showed that the underlying mechanism mainly involved AGE-RAGE signalling pathway in diabetic complications, TNF signalling pathway and MAPK signalling pathway. Significantly, the MAPK signalling pathway was considered as the key pathway. In conclusion, we speculated that dapagliflozin played a therapeutic role in T2DM complicated with CVD mainly through MAPK signalling pathway. This study preliminarily reveals the possible mechanism of dapagliflozin in the treatment of T2DM complicated with CVD and provides a theoretical basis for future clinical research.     

生信分析筛选心衰氧化应激相关关键[]基因并预测候选药物

为探究氧化应激相关基因在心力衰竭发生发展中的作用,并发现核心基因进行靶基因药物预测。从GEO数据库下载GSE120895基因表达图谱,通过GEO2R筛选差异表达基因,将差异表达基因与GeneCard数据库中筛选的氧化应激相关基因取交集,得到心力衰竭氧化应激相关差异表达基因,利用R软件对差异表达基因进行GO及KEGG分析,利用Cytoscape进行PPI网络的模块以及关键基因的筛选。之后在GSE17800基因表达图谱中验证关键基因的表达,并针对关键基因进行相互作用药物预测。差异表达基因与氧化应激相关基因取交集后,共筛选出52个上调的氧化应激相关差异表达基因,在此基础上,筛选出ACTB,STAT3,FN1,EDN1,CAT共5个关键基因,在GSE17800基因表达图谱中验证后,针对4个关键基因预测了19个靶基因潜在药物。总之,本研究通过生物信息学方法鉴定关键基因,并预测潜在治疗药物,从而为了解心力衰竭的分子机制及其诊治方法提供新的见解。     

High-abundance mRNAs in Apis mellifera: comparison between nurses and foragers

A large volume of honey bee (Apis mellifera) tag-seq was obtained to identify differential gene expression via Solexa/lllumina Digital Gene Expression tag profiling (DGE) based on next generation sequencing. In total, 4,286,250 (foragers) and 3,422,327 (nurses) clean tags were sequenced, 24,568 (foragers) and 13,134 (nurses) distinct clean tags could not be match to the reference database, and 7508 and 6875 mapped genes were detected in foragers and nurses respectively. 7045 genes were found differentially expressed between foragers and nurses. Of those genes, 1621genes had significantly different expression, that is, they showed an expression ratio (foragers/nurses) of more than 2 and FDR (False Discovery Rate) of less than 0.001. We identified 101 genes that were uniquely expressed in foragers, and 9 genes that were only expressed in nurses. We performed the Gene Ontology (GO) category and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, and found 415 genes with annotation terms linked to the GO cellular component category. 200 components of KEGG pathways were obtained, including 21 signaling pathways. The PPAR signaling pathway was the most highly enriched, with the lowest Q-value.     

Bioinformatics analysis reveals disturbance mechanism of MAPK signaling pathway and cell cycle in Glioblastoma multiforme

Background & objectives

To analyze the reversal gene pairs and identify featured reversal genes related to mitogen-activated protein kinases (MAPK) signaling pathway and cell cycle in Glioblastoma multiforme (GBM) to reveal its pathogenetic mechanism.

Methods

We downloaded the gene expression profile GSE4290 from the Gene Expression Omnibus database, including 81 gene chips of GBM and 23 gene chips of controls. The t test was used to analyze the DEGs (differentially expressed genes) between 23 normal and 81 GBM samples. Then some perturbing metabolic pathways, including MAPK (mitogen-activated protein kinases) and cell cycle signaling pathway, were extracted from KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway database. Cancer genes were obtained from the database of Cancer Gene Census. The reversal gene pairs between DEGs and cancer genes were further analyzed in MAPK and cell cycle signaling pathway.

Results

A total 8523 DEGs were obtained including 4090 up-regulated and 4433 down-regulated genes. Among them, ras-related protein rab-13(RAB13), neuroblastoma breakpoint family member 10 (NBPF10) and disks large homologue 4 (DLG4) were found to be involved in GBM for the first time. We obtained MAPK and cell cycle signaling pathways from KEGG database. By analyzing perturbing mechanism in these two pathways, we identified several reversal gene pairs, including NRAS (neuroblastoma RAS) and CDK2 (cyclin-dependent kinase 2), CCND1 (cyclin D1) and FGFR (fibroblast growth factor receptor). Further analysis showed that NRAS and CDK2 were positively related with GBM. However, FGFR2 and CCND1 were negatively related with GBM.

Interpretation & conclusions

These findings suggest that newly identified DEGs and featured reversal gene pairs participated in MAPK and cell cycle signaling pathway may provide a new therapeutic line of approach to GBM.     

Transcriptomic Analysis of Insulin-Sensitive Tissues from Anti-Diabetic Drug Treated ZDF Rats,a T2DM Animal Model

Gene expression changes have been associated with type 2 diabetes mellitus (T2DM); however, the alterations are not fully understood. We investigated the effects of anti-diabetic drugs on gene expression in Zucker diabetic fatty (ZDF) rats using oligonucleotide microarray technology to identify gene expression changes occurring in T2DM. Global gene expression in the pancreas, adipose tissue, skeletal muscle, and liver was profiled from Zucker lean control (ZLC) and anti-diabetic drug treated ZDF rats compared with those in ZDF rats. We showed that anti-diabetic drugs regulate the expression of a large number of genes. We provided a more integrated view of the diabetic changes by examining the gene expression networks. The resulting sub-networks allowed us to identify several biological processes that were significantly enriched by the anti-diabetic drug treatment, including oxidative phosphorylation (OXPHOS), systemic lupus erythematous, and the chemokine signaling pathway. Among them, we found that white adipose tissue from ZDF rats showed decreased expression of a set of OXPHOS genes that were normalized by rosiglitazone treatment accompanied by rescued blood glucose levels. In conclusion, we suggest that alterations in OXPHOS gene expression in white adipose tissue may play a role in the pathogenesis and drug mediated recovery of T2DM through a comprehensive gene expression network study after multi-drug treatment of ZDF rats.     

2型糖尿病模型大鼠肝脏脂代谢相关基因的筛选及生物信息学分析

肝的脂肪代谢异常和胰岛素抵抗（insulin resistance,IR）对促进2型糖尿病（type 2 diabetes mellitus,T2DM）的发生与发展具有显著影响。但此过程复杂,参与调控基因目前尚未完全清楚。有研究表明,脂肪酸分解、氨基酸代谢、肝糖原合成等生物过程对糖尿病的形成具有促进作用。为了阐明这一调控机制,本文通过基因芯片技术研究GK（Goto-Kakizaki）大鼠和WKY（Wistar-Kyoto）大鼠肝差异基因对肝的脂肪代谢和胰岛素抵抗的影响,探讨可引起2型糖尿病发病的分子机制。从基因表达数据库（GEO）获取GSE13271基因表达谱,并对原始数据进行标准化处理。通过GO（Gene Ontology）、KEGG（Kyoto Encyclopedia of Genes and Genomes enrichment）、String和Cytoscape软件对差异表达基因进行功能分析。结果从GK和WKY大鼠中分别获得179和278个差异基因,同时从排名前10的路径中筛选出21个差异基因(Aldh1a1, Cyp2c22, bp2,Fabp7,Cyp4a3, Acot1, Acot2,Hsd17b2, Ech1, Hmgcl,Bdh1, Crot, Pex11a, Cpt1a, Hadhb, Gda, Elovl2, Prodh, Agpat3, Sardh, Pigu),将这些基因与前10个的GO term取交集。最终得到10个显著差异基因(Aldh1a1, Fabp2, Acot1, Acot2, Ech1, Hmgcl, Bdh1, Crot, Cpt1a, Hadhb),功能分析结果显示,肝组织相关基因通过一系列生物过程对肝的脂肪代谢和胰岛素抵抗产生调节作用,从而也为临床糖尿病的治疗以及新作用靶点的发现提供更多参考依据。     

Differential gene expression in pancreatic tissues of streptozocin-induced diabetic rats and genetically-diabetic mice in response to hypoglycemic dipeptide cyclo (His-Pro) treatment

Diabetic studies are mostly interested in gene expression in the pancreas, the site of insulin secretion that regulates blood glucose levels. However, a single gene approach has been ruled out for many years in discovering new genes or the molecular networks involved in the induction process of diabetes. To understand the molecular mechanisms by which cyclo (His-Pro) (CHP) affects amelioration of diabetes mellitus, we performed gene expression profiling in the pancreatic tissues of two diabetic animal models, streptozocin (STZ)-induced diabetic rats (T1DM) and genetically-diabetic (C57BL/6J ob/ob) mice (T2DM). To understand the healing process of these diabetic rodents, we examined the effects of CHP on various gene expression in pancreatic tissues of both animal models. Our microarray analysis revealed that a total of 1,175 genes were down-regulated and 629 genes were up-regulated in response to STZ treatment, and the altered expression levels of numerous genes were restored to normal state upon CHP treatment. In particular, 476 genes showed significantly altered gene expression upon CHP treatment. In a functional classification, 7,198 genes were counted as differentially expressed in pancreatic tissues of STZ- and CHP-treated rats compared with control, whereas 1,534 genes were restored to normal states by CHP treatment. Microarray data demonstrated for the first time that overexpression of the genes encoding IL-1 receptor, lipid metabolic enzymes (e.g. Mte1, Ptdss1, and Sult2a1), myo-inositol oxygenase, glucagon, and somatostatin as well as down-regulation of olfactory receptor 984 and mitochondrial ribosomal protein, which are highly linked to T1DM etiology. In genetically-diabetic mice, 4,384 genes were altered in gene expression by more than 2-fold compared to the control mice, when counted differentially expressed. In genetically-diabetic mice, 4,384 genes altered in expression by higher than 2-fold were counted as differentially expressed genes in pancreatic tissues of CHP-treated mice. On the other hand, 2,140 genes were up-regulated and 2,244 genes were down-regulated by CHP treatment. The results of the microarray analysis revealed that up-regulation of IL-2, IL12a, and leptin receptor and down-regulation of PIK3 played important physiological roles in the onset of T2DM. In conclusion, we hypothesize that CHP accelerates alterations of gene expression in ameliorating diabetes and antagonizes those that induces the disease.     

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.     

基于单细胞转录组的多级别胶质瘤异质性及免疫微环境分析揭示了潜在的预后生物标志物

脑胶质瘤(glioma)是中枢神经系统最常见的内在肿瘤,具有发病率高、预后较差等特点。本研究旨在鉴定多形性胶质母细胞瘤(glioblastoma multiforme,GBM)和低级别胶质瘤(lower-grade gliomas, LGG)之间的差异表达基因(differentially expressed genes, DEGs),以探讨不同级别胶质瘤的预后影响因素。从NCBI基因表达综合数据库中收集了胶质瘤的单细胞转录组测序数据,其中包括来自3个数据集的共29 097个细胞样本。对于不同分级的人脑胶质瘤进行分析,经过滤得到21 071个细胞,通过基因本体分析、京都基因与基因组百科全书途径分析,从差异表达基因中筛选出70个基因,我们通过查阅文献,聚焦到delta样典型Notch配体3 (delta like canonical Notch ligand 3,DLL3)这个基因。基于TCGA的基因表达谱交互分析(gene expression profiling interactive analysis, GEPIA)数据库用于探索LGG和GBM中DLL3基因的表达差异,采用基因表达...     

MiR-128-3p accelerates cardiovascular calcification and insulin resistance through ISL1-dependent Wnt pathway in type 2 diabetes mellitus rats

Vascular calcification is highly prevalent in patients with type 2 diabetes mellitus (T2DM), one of the most common chronic diseases with high morbidity and mortality. In recent years, microRNAs have been widely reported as potential biomarkers for the diagnosis and treatment of T2DM. We hypothesized that miR-128-3p is associated with cardiovascular calcification and insulin resistance (IR) in rats with T2DM by targeting ISL1 via the Wnt pathway. Microarray analysis was adopted to identify differentially expressed genes related to T2DM. T2DM models were induced in rats. Blood samples from normal and T2DM rats were used to detect islet β-cell function, islet sensitivity, and calcium content. Next, islet tissues were obtained to identify the expression of miR-128-3p, ISL1, and the Wnt signaling pathway- and apoptosis-related genes. Finally, apoptosis of islet β-cells was determined by flow cytometry. Through microarray analysis of GSE27382 and GSE23343, ISL1 was found to be downregulated in T2DM. In blood samples from T2DM rats, basic biochemical indicators, IR, and calcium content were increased, and islet sensitivity and islet β-cell function were decreased. Furthermore, upregulation of miR-128-3p and ISL1 gene silencing promoted the expression of Wnt-1, β-catenin, GSK-3β, and Bax and the phosphorylation of β-catenin and GSK-3β, inhibited c-fos, PDX-1, and Bcl-2 expression, and enhanced cell apoptosis. The key findings of our study demonstrate that miR-128-3p aggravates cardiovascular calcification and IR in T2DM rats by downregulating ISL1 through the activation of the Wnt pathway. Thus, miR-128-3p may serve as a potential target for the treatment of T2DM.