Gene microarray analysis of expression profiles in Suberoyllanilide hyroxamic acid-treated Dendritic cells

Objective

The purpose of this study is to provide a further theoretical basis for the role of Suberoyllanilide hyroxamic acid (SAHA) affect on Dendritic cells (DCs).

Methods

We first downloaded the GSE74306 microarray data, which was about the effect of SAHA act on DCs, from the Gene Expression Omnibus database. Then we analyzed the differential expression genes (DEGs) between SAHA-treated DCs and SAHA-untreated DCs by limma package of R software; The Database for Annotation, Visualization and Integrated Discovery was used to analyze the Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways for these DEGs. The protein protein interaction (PPI) network was constructed by using STRING database, Cytoscape 3.6.1 software was used to dispose the PPI network for visualization. Finally, we determine the Hub genes in the PPI network according by the degree centrality and betweenness centrality, which were calculated by the CentScaPe 2.2 plug-in of Cytoscape 3.6.1 software.

Result

There were 551 DEGs between SAHA-treated DC cells and SAHA-untreated DC cells, including 357 upregulated genes and 194 downregulated genes. These DEGs genes were enriched in 115 Go terms (Biological Process, 51; Cellular Component, 35 and Molecular Function, 29) and a total of 16 pathways. Glutathione metabolic process, Glutathione metabolism pathway, Rheumatoid arthritis pathway and Systemic lupus erythematosus pathway were most significant function clusters. In the PPI network, Rad51, Src, and Eno2 were Hub genes.

Conclusion

The biological function and KEGG pathway enriched by DEGs may reveal the molecular mechanism of SAHA acting on DC cells. Its Hub genes, Src, Rad51 and Eno2, were expected to be new targets for SAHA therapeutic effects. However, it still need to be confirmed by the next more rigorous molecular biological experiments research.