Genome-wide analysis reveals the active roles of keratinocytes in oral mucosal adaptive immune response

Abstract

To elucidate the roles of oral keratinocytes in the adaptive immune response of oral mucosa, global gene expression analysis was performed by microarray technique and integrating computational methods, including hierarchical clustering, biological process Gene Ontology analysis, Kyoto Encyclopedia of Genes and Genomes pathway analysis, self-organizing maps (SOMs) and biological association network analysis (BAN). Raw data from microarray experiments were uploaded to the Gene Expression Omnibus Database, http://www.ncbi.nlm.nih.gov/geo/ (GEO accession GSE28035). We identified 666 differentially expressed genes in the early stage (48 h) and 993 in the late stage (96 h) of the oral mucosal adaptive immune response. The analysis revealed that oral keratinocytes exerted diverse biological functions in different stages of immune response. Specifically, in 48 h the differentially expressed genes encompassed an array of biological ontology associated with immune response, such as antigen processing and presentation, and positive regulation of T-cell-mediated cytotoxicity. Several pathways which have been reported to be critical in inflammation, including mitogen-activated protein kinase pathway, were activated. Furthermore, after BAN construction, some putative hub genes and networks such as interleukin-1α and its subnetwork were recognized. Taken together, these results give substantial evidence to support the active roles of keratinocytes in the oral mucosal adaptive immune response.