Assessment of the genes and molecular mechanisms of B cells activation through systems biology approaches

Two classes of T helper lymphocytes, Th1 and Th2, have different roles in B cell activation based on specific cytokines. To understand the difference of molecular mechanisms of B cell activation, the microarray dataset of B cells co-cultured with type 1 and 2 T helper, Be1 and Be2, were investigated. After quality assessment, using the GEO2R tool, the GSE84948 dataset was re-analyzed. Genes with adjusted p-value ⩽ 0.05 were assumed as differentially expressed (DE). The protein-protein interaction (PPI) networks were constructed using CluePedia plugin of Cytoscape, and analyzed by NetworkAnalyzer tool and MCODE plugin. Using ClueGO plugin of Cytoscape software, gene ontology (GO) analysis was performed. The comparison of Be1 and Be2 cells with naive B cells revealed 8742 and 8748 DE genes, respectively. The topology analysis of PPI networks predicted central genes. Among these, Jak3, Actrt3, and Pik3cb genes were determined as crucial genes in Be1 network. Prkx, Smarca4, and Jak2 genes were defined as critical genes in Be2 PPI network. GO analysis with PPI networks genes resulted in the promotion of immune system activation. In conclusion, we explored holistic methods for molecular assay of the difference between B cell activation mechanisms with Th1 and Th2.