Multiple Transcriptome Data Analysis Reveals Biologically Relevant Atopic Dermatitis Signature Genes and Pathways

Several studies have identified genes that are differentially expressed in atopic dermatitis (AD) compared to normal skin. However, there is also considerable variation in the list of differentially expressed genes (DEGs) reported by different groups and the exact cause of AD is still not fully understood. Using a rank-based approach, we analyzed gene expression data from five different microarray studies, comprising a total of 127 samples and more than 250,000 transcripts. A total of 89 AD gene expression signatures ‘89ADGES’, including FLG gene, were identified to show dysregulation consistently across these studies. Using a Support Vector Machine, we showed that the ‘89ADGES’ discriminates AD from normal skin with 98% predictive accuracy. Functional annotation of these genes implicated their roles in immune responses (e.g., betadefensin, microseminoprotein), keratinocyte differentiation/epidermal development (e.g., FLG, CORIN, AQP, LOR, KRT16), inflammation (e.g., IL37, IL27RA, CCL18) and lipid metabolism (e.g., AKR1B10, FAD7, FAR2). Subsequently, we validated a subset of signature genes using quantitative PCR in a mouse model. Using a bioinformatic approach, we identified keratinocyte pathway over-represented (P = <0.0006) among the 89 signature genes. Keratinocytes are known to play a major role in barrier function due to their location in the epidermis. Our result suggests that besides immune- mediated pathway, skin barrier pathways such as the keratinocyte differentiation pathway play a key role in AD pathogenesis. A better understanding of the role of keratinocytes in AD will be important for developing novel “barrier therapy” for this disease.