Skip to main content

Meta-analysis derived atopic dermatitis (MADAD) transcriptome defines a robust AD signature highlighting the involvement of atherosclerosis and lipid metabolism pathways



Atopic dermatitis (AD) is a common inflammatory skin disease with limited treatment options. Several microarray experiments have been conducted on lesional/LS and non-lesional/NL AD skin to develop a genomic disease phenotype. Although these experiments have shed light on disease pathology, inter-study comparisons reveal large differences in resulting sets of differentially expressed genes (DEGs), limiting the utility of direct comparisons across studies.


We carried out a meta-analysis combining 4 published AD datasets to define a robust disease profile, termed meta-analysis derived AD (MADAD) transcriptome.


This transcriptome enriches key AD pathways more than the individual studies, and associates AD with novel pathways, such as atherosclerosis signaling (IL-37, selectin E/SELE). We identified wide lipid abnormalities and, for the first time in vivo, correlated Th2 immune activation with downregulation of key epidermal lipids (FA2H, FAR2, ELOVL3), emphasizing the role of cytokines on the barrier disruption in AD. Key AD “classifier genes” discriminate lesional from nonlesional skin, and may evaluate therapeutic responses.


Our meta-analysis provides novel and powerful insights into AD disease pathology, and reinforces the concept of AD as a systemic disease.

Peer Review reports


Atopic dermatitis (AD) is the most common inflammatory skin disease (4–7 % prevalence in adults, and ~15 % in children), with a large unmet need for safer and more effective treatments [17]. Immune and barrier abnormalities characterize AD, with Th2/Th22 cytokine activation, increased hyperplasia and significant decreases in differentiation markers. These observations have led to two competing pathogenic hypotheses [1, 8, 9], although recent studies characterizing AD primarily as an immune-driven disease have shown reversal of barrier defects following specific and non-specific therapeutic interventions [1013].

Genomic expression profiling using gene-arrays and real time (RT)-PCR has been widely used to identify gene alterations in lesional (LS) and non-lesional (NL) AD compared to normal skin to better understand interactions between activation of cytokine pathways and epidermal abnormalities [6, 7, 10, 12, 14, 15]. Similar to other diseases, the AD phenotype/transcriptome can be defined genomically by differentially expressed genes (DEGs) between LS and NL skin [15]. A robust transcriptome was established as a powerful tool in identifying core psoriasis pathogenesis and evaluating the efficacy of targeted therapeutics at transcriptomic level [16, 17]. The high rates of placebo effect in AD patients contrasts with a worsening of disease phenotype at the transcriptomic level [10], reinforcing the importance of a robust disease transcriptome against which therapeutic effects can be evaluated [18]. Genomic profiling may also be used to predict therapeutic responses, as in juvenile idiopathic arthritis, in which profiling correctly predicted therapeutic responses at 6 months [19].

Nevertheless, high-throughput genomic analyses are vulnerable to multiple biases, including random noise, biological heterogeneity, and differences in experimental procedures (biopsies, hybridization, etc.), leading to remarkably little overlap between DEGs in similar scale studies [2022]. A meta-analysis approach that combines microarray data from independent yet similarly designed studies allows one to overcome these variations, ultimately increasing the power and reproducibility of the transcriptome [2325]. While several meta-analysis methods exist for combining microarray data from independent studies [26], the random-effects model (REM) has been established as one of the most suitable for heterogeneous studies [25, 27, 28].

We applied a REM meta-analysis model including 4 published AD microarray studies (including 97 samples; 54 LS and 43 NL; 41 paired) to determine core pathogenic elements and new disease associated genes [9, 10, 12, 15], resulting in the Meta-Analysis Derived AD (MADAD) disease transcriptome, a robust active disease signature of 595 DEGs, including 86 that were not previously detected by any individual study.


Sample collection

All samples were collected according to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) statement [29]. A total of 28 datasets were detected in Gene Expression Omnibus (GEO), but only datasets including LS and NL skin samples of AD patients were retained. Datasets run on platforms other than the HGU133Plus 2 chip, subject to treatments, or with non-randomly selected NL or LS samples (e.g. FLG homozygous/heterozygous loss of function mutation), and datasets without NL samples were excluded. When overlapping samples were found between datasets, only one copy was kept. Four microarray datasets satisfied the established criteria (GSE32924, GSE36842, GSE58558, GSE59294) [9, 10, 12, 15], including 97 samples (54 LS and 43 NL), which coincidently have been carried out by our group.

Pre-processing and expression analysis

Pre-processing and statistical analyses were performed using R and Bioconductor packages [30, 31]. Raw expression data were combined, summarized, and normalized using GCRMA [32]. Batch effects between datasets were adjusted for by the empirical Bayes method ComBat/SVA [3235]. Agreement of the individual studies raw microarray data was estimated by the Integrated Correlation Coefficient Analysis, which produces the general Integrated Correlation Coefficient (ICC), representing agreement between studies, and can be interpreted in the same way as Pearson correlation coefficient. The ICC was used to eliminate background noise prior to the analysis, by excluding genes with incoherent behavior across studies [36].

For each study, estimation of differences in expression levels of LS vs NL skin was performed using the mixed-effect framework of the limma package.


The formal random effects model (REM), described by Choi [27]. was applied to estimate the true effect size for each probe (see Additional file 1). These estimation and calculation steps were performed using the package GeneMeta. P-values were adjusted for multiple testing using the Benjamini-Hochberg procedure [37], with criteria for DEGs of absolute fold change (|FCH|) ≥ 2.0 and a false discovery rate (FDR) ≤ 0.05.


The MADAD transcriptome was subject to multiple downstream analysis methods. Integration-driven discovery/IDD-DEGs were defined as those not identified in any of the included studies. To explore functional annotations, overrepresentation analysis was carried out for BP GO-terms and KEGG pathways (both in DAVID) [38], Ingenuity Pathways (IPA –, as described) [25], and on previously reported gene-sets [39].

The normalized LS and NL expression data were subject to Weighted Gene Co-Expression Networks Analysis (WGCNA) to detect clusters (modules) of correlated genes and respective hub genes [40]. These modules were subject to trait correlation and corresponding gene-set overrepresentation analysis.

Meta Threshold Gradient Directed Regularization (MTGDR) method [41, 42] was used to select disease-associated genes while allowing for varied estimates of those genes across different experiments, as previously published [41, 42].

Primers and probes were designed for RT-PCR as previously described (see Additional file 1 and Additional file 2) [43]. Ribosomal protein, large, P0/RPLP0 normalized RT-PCR expression data was analyzed using a mixed-effect framework after log2-transformation.

Lipid metabolism genes were defined as genes related to one of the four groups: Ceramides, Free Fatty Acids, Sphingolipids, and Cholesteryl Esters in the Gene Cards database ( We included all genes with a relevance score ≥ 10 [44, 45]. Pairwise Pearson correlations were calculated between the gene and patient specific deregulations. Multivariate correlations between sets of genes were calculated making use of the gene set specific μ-scores calculated by the muStat package [46].


Data collection

Two major repositories (GEO Omnibus and ArrayExpress) were queried to identify studies with expression profiles of LS and NL punch biopsies from AD patients. Four studies (GSE59294, GSE58558, GSE36842, GSE32924; including 97 samples, with 54 LS and 43 NL, 41 paired) met all inclusion criteria following the PRISMA guidelines (see Methods, Additional file 3: Figure E1 and Additional file 4: Table E1) [9, 10, 12, 15]. We only used data from chronic lesions in the analysis. No significant differences in disease severity (as measured by Scoring of AD/SCORAD and Eczema Area and Severity/EASI indices) or IgE levels were found across patients who all had moderate-to-severe AD (SCORAD > 25; EASI > 12).

Meta-analysis framework

The Venn diagram in Fig. 1a represents the overlap of DEGs identified by the individual studies, including 25 consensus DEGs (see Additional file 5: Table E2). Besides inter-study variation, sources of inconsistency in DEGs include choice of model, annotation, cut-off, and non-uniform pre-processing steps [47]. To combine results of individual studies and to aggregate robust DEGs with reliable effect sizes, we chose a meta-analysis approach [21, 27].

Fig. 1

a Approximate area proportional Venn diagram with overlaps of the differentially expressed genes of the four individually analyzed datasets under the same threshold (|FCH| ≥ 2, FDR ≤ 0.05). b Ingenuity Canonical Pathway overrepresentations compared between this MADAD transcriptome and the included dataset from Khattri et al. 2014 (GSE58558). The bars represent –log10 transformed Benjamini-Hochberg adjusted p-values termed –log(FDR), and the red threshold line indicates the FDR ≤ 0.05 cut-off. c Cutaneous localization of the MADAD transcriptome as compared to the epidermis and dermis related genes defined by Esaki et al. [63]

To address these sources of inconsistency, we applied a uniform pre-processing pipeline to the combined datasets, incorporating noise reduction techniques as proposed [23, 25, 36, 47] (see Methods, Additional file 3: Figure E1C, and Additional file 1). This pre-processing increased the average of the pairwise maximum canonical correlations between each pair of studies from 0.76 to 0.85.

To combine heterogeneous effect sizes across studies, a random effects model/REM was chosen over fixed effects model. The use of REM is supported by the sample distribution of the Cochran’s Q statistic [27], which departs substantially from the theoretical χ 23 -distribution under the assumption of homogeneity of study effects. Furthermore, a Kolmogorov-Smirnov-test also rejected the equality of these two distributions (p = 2.22−15; D = 0.235) (Additional file 6: Figure E2).

The Meta-analysis-derived AD (MADAD) transcriptome

Applying the meta-analysis approach to 4 studies and 97 samples, we identified a set of 595 DEGs (387 up- and 208 downregulated) using the classical |FCH| ≥ 2 and FDR ≤ 0.05 criteria, representing a robust profile defined as the meta-analysis derived AD (MADAD) transcriptome (see Additional file 7: Table E3 for the entire MADAD DEGs list). Among the highest dysregulated genes are key AD genes, including markers of general inflammation (MMP12), specific T helper activation (e.g. Th2/CCL18, Th1/IFN/CXCL10, Th17/PI3/elafin, Th17/Th22 S100A7/A8/A9), and markers of epidermal proliferation (KRT16, Mki67).

To show the robustness of the MADAD transcriptome versus a single study, we conducted an over-representation analysis (see Methods) of ingenuity canonical pathways (IPA) and of previously published immune and barrier gene-sets [12, 48]. Figure 1b illustrates the comparison of the largest available’ dataset (Khattri 2014 – GSE58558) and the MADAD transcriptome; similar results were obtained for the other three included datasets. Overall, the MADAD yields a more significant over-representation of key immune pathways, such as Granulocyte Adhesion and Diapedesis, T Cell Receptor Signaling and differentiation, iCOS-iCOSL Signaling, and IL-12, IL-17, and IL-9 signaling. Innate pathways (e.g. Role of Pattern Recognition Receptors in Recognition of Bacteria), which are associated with AD [49, 50], were also represented (Fig. 1b and Additional file 8: Tables E4 and Additional file 9: Table E5). Interestingly, Atherosclerosis Signaling was the second highest IPA pathway, and includes genes previously associated with vascular inflammation, such as IL-37, SERPINA1, S100A8, selectin E/SELE, lipoprotein lipase/LPL, and MMP1/3/9 [5162]. An equivalent analysis of previously reported immune and barrier gene-sets [12] (Additional file 10: Figure E3A) similarly showed increased sensitivity and representation in the MADAD transcriptome compared to the largest data set, including IFNα, IL-4, immune genes, cytokine-treated keratinocytes and epidermal differentiation gene subsets (Additional file 10: Figure E3A). Thus, the MADAD transcriptome provides a more robust AD-specific signal, consistent with known disease pathology.

Comparison of the MADAD transcriptome to the recently described epidermal and dermal layer-specific transcriptomes linked 223 and 231 DEGs to the epidermis and dermis, respectively (Fig. 1c, Additional file 7: Table E3) [63]. Of these, 55.3 % of the upregulated genes were epidermal, whereas 67.5 % of the downregulated genes were dermal (P = 2.6 × 10−5, Fisher’s exact test), as has been previously noted [63].

Of the top 25 up-regulated MADAD DEGs, 17 encode epidermal components, including key antimicrobial genes (DEFB4A, PI3/elafin, S100A9) (Additional file 7: Table E3). Top 25 up-regulated dermal genes include those related to collagen production (COL4A4, COL6A5) and inflammation (GZMB, OASL) [64, 65]. Among top down-regulated epidermal genes are structural/lipid-related genes (LEP, FABP7, ELOVL3), and genes linked to epidermal differentiation (CLDN8) [6668]. Several genes associated with the pathogenesis of AD (e.g. IL-22, OX40L and TSLP) and reported in the layer-specific AD transcriptomes [63] were not detected in the MADAD transcriptome, most likely due to low expression levels of cytokines on whole tissue microarrays, a known limitation of microarrays [24, 69], that cannot be overcome by the meta-analysis approach.

Integration-Driven Discovery Genes in the MADAD transcriptome

The MADAD comprises a subset of IDD-DEGs, which were not detected by any of the individual datasets using identical cutoffs [9, 10, 12, 15]. These 86 IDD-DEGs consist of 45 up- and 41 downregulated genes (Table 1). Using IPA and GOterm overrepresentation analyses on the IDD-DEGs, we found in addition to the expected inflammatory processes (e.g. Chemokine Signaling), several pathways related to lipid and fatty acid metabolic processes (Additional file 10: Figure E3B, Additional file 11: Table E6).

Table 1 Integration-Driven Discovery (IDD) Genes in the MADAD transcriptome, with indication of compartmental allocation as defined by Esaki et al. [63]

We further validated, using RT-PCR, several IDD-DEGs and top MADAD DEGs with plausible biologic relevance to AD, including immune (CCL8, CD1E, IL-37, IL-36G), structural (AQP5, a water channel) and lipid (FAR2, ELOVL3, FA2H) genes, which encode enzymes involved in fatty acid and ceramide metabolism (Fig. 2a). Immune-related genes (e.g. CCL8, LCP2, CD1E) and IL-36G, which was recently associated with psoriasis pathogenesis [70], were increased in AD LS skin. The only gene that showed increased expression in NL vs LS skin was IL-37, a negative immune regulator [71], consistent with past reports [72]. All structural and lipid metabolism genes showed downregulation in LS AD skin (P ≤ 0.05 for all), reinforcing the role of lipids and water channels in preserving barrier function in AD [73, 74].

Fig. 2

a Barplot representing means ± SEM of RT-PCR expression values for selected MADAD DEGs, normalized by hARP and log2 transformed. Underlined genes represent those not detected in any of the included datasets (IDD-genes). b Multivariate correlation between lipid and Th2-related genes dysregulations. Stars indicate the corresponding significance levels (+ ≤0.1, * ≤ 0.05, ** ≤0.01), while the correlations are shown from negative (−0.05, dark blue) to positive (0.05, dark red). The correlation data structure was assessed with Euclidean distance and complete agglomeration hierarchical clustering. Note that the right section gives a schematic overview of the gene set analysis and the corresponding correlation coefficients for the overall Th2 vs Lipid comparison and the sub-comparisons of C1 vs Lipid and C2 vs Lipid (see Additional file 15: Table E10 for corresponding BH adjusted p-values)

Suppression of lipid-related genes is coupled to increased Th2 activation

While we have previously shown inverse correlations between immune activation and terminal differentiation genes in AD lesions [11, 15], the relationship between immune activation and lipid metabolism genes has not been assessed in skin.

We thus investigated the relationship between expression of Th2-specific and epidermal lipid metabolism genes in the MADAD transcriptome. First we performed an unbiased database search for lipid related terms (see Methods) and chose genes with highest relevance scores. The correlation structure between lipid metabolism and immune genes in LS and NL skin is shown in Fig. 2b.

Unsupervised clustering identified two main sub-clusters in the Th2 gene-set and three main sub-clusters in the lipid metabolism gene-set. Cluster 3 consists only of SPTLC2, the only lipid metabolism gene with up-regulated expression in the MADAD transcriptome. This gene encodes serine palmitoyltransferase (SPT), a rate-limiting enzyme in sphingomyelin synthesis, whose elevated expression has been associated with increased barrier defects, including in AD [7577].

Multivariate u-statistics correlate the dysregulation of immune and lipid metabolism genes [46]. Overall, an inverse correlation of −0.46 (p = 0.003) was obtained between Th2 and lipid subsets, supporting a proposed model of Th2 cytokine effects on lipid suppression [78]. Among Th2 genes, cluster 1 showed the highest negative correlation with lipid metabolism genes (−0.49; p = 0.001), and includes key AD genes (CCL22, IL-7R, and IL-4R) [63, 7981]. Targeting IL-4R shows promise in early clinical trials as a possible therapeutic target for moderate-to-severe AD, and is now in phase 3 clinical trials for this disease (NCT02277743) [10, 18].

Weighted gene co-expression network analysis

Weighted gene co-expression network analysis (WGCNA) offers insights into disease pathogenesis by studying weighted co-expression of genes within tissue samples [40]. This technique requires a large sample size [40, 82], which we were able to apply here for the first time in AD (see Methods).

Using WGCNA, 21 distinct sub-networks were identified, and each was correlated with age, IgE level, and disease severity (measured by Scoring of Atopic Dermatitis/SCORAD) index (Fig. 3a-b, Additional file 12: Tables E7 and Additional file 13: Table E8). In LS skin, several networks showed significant positive correlations with SCORAD, including viral (M13) and innate immune response processes (M17), emphasizing cutaneous immune reactions to viral and/bacterial pathogens in LS AD skin. Proliferation and epidermal processes were also correlated with SCORAD (M10, M11), with a trend for negative correlations between SCORAD and structural epidermal modules (M4, M12). In NL skin, immune-related networks such as cytokine receptor signaling pathway and innate immune response modules (M4, M2) [63, 83] were also positively correlated with SCORAD, while extracellular matrix (ECM) organization (M10) and ECM-receptor interaction (M9) modules were negatively correlated (Fig. 3b). Interestingly, a trend for a positive correlation was observed between SCORAD and Staphylococcus aureus (S.aureus) infection in NL skin (M8), with possible clinical relevance, since AD patients are often colonized with S.aureus, even in NL skin [4, 84]. Similar but weaker correlations were seen with IgE, while minimal correlations where found with age (Fig. 3a-b).

Fig. 3

Module eigen-genes from weighted gene co-expression network analysis (WGCNA), correlated with Age, SCORAD and IgE, including annotation of top overrepresented BP GO-terms, KEGG Pathways or curated inflammatory skin diseases related gene-sets in (a) Lesional and (b) Non-lesional co-expression networks. Numbers in each cell represent the correlations with the respective trait and the associated p-values (see Additional file 16: Tables E11 and Additional file 17: Table E12 for corresponding BH adjusted p-values). c Radviz plots showing the separation of the lesional (LS) and non-lesional (NL) samples, as defined by the 19 genes signature identified by the MTGDR classifier (see Additional file 18: Table E13 for raw MTGDR results)

A robust genomic AD classifier

Because clinical trials with specific and broad therapeutics are being tested in AD [5, 85], establishing a reliable gene set to discriminate between diseased and normal-appearing skin in AD is crucial. Here, we applied the classification algorithm Meta Threshold Gradient Directed Regularization (MTGDR) [41, 42], to determine the smallest set of genes that distinguishes LS from NL AD skin (see Methods). The final model identifies 19 genes (Fig. 3c and Table 2), including both epidermal and dermal genes, emphasizing the importance of both compartments to the disease model.

Table 2 MTGDR robust classifier genes, compared to treatment response, and indication of compartmental allocation as defined by Esaki et al. [63]

To assess the translational validity of these 19 discriminating genes, we investigated the effect of various therapies on gene expression using previously published studies with Dupilumab (300 mg, 4 weeks of treatment) [10], Cyclosporin A (CsA; at 2 and 12 weeks) [12], and narrow band UVB/NB-UVB (12 weeks of treatment) [11], shown in Table 2. Reversal of disease phenotype to the NL state was observed in 18 of the 19 discriminating genes. Overall, a higher reversal was seen with CsA, a broad immune suppressant [86, 87], particularly after 12 weeks of treatment, with an average recovery of 94.64 % (see Table 2).

Among the genes down-regulated with treatment are inflammatory genes (S100A9, SELE) previously associated with AD [9, 59, 88]. Genes up-regulated with treatment (particularly with long-term CsA) include perilipin [89] and hydroxysteroid dehydrogenase [90], which are involved in lipid and steroid metabolism, respectively.

The LS and NL dysregulation in the MADAD correlates with the treatment change induced by CsA (r = 0.99; p = 2.5*10−14), dupilumab (r = 0.45; p = 0.189), and UBV (r = 0.90,p = 0.001) treatments, leading to restoration of gene expression similar to that seen in NL skin (Additional file 14: Figure E4).


A meta-analytic approach utilizes statistical processing and analysis to merge microarray studies from various populations and investigators, resulting in a single value that represents the estimated differential expression level of a gene between LS and NL skin. Combining multiple studies in a meta-analysis produces findings that more precisely reflect the differential expression of genes in a population, representing an accurate molecular characterization of a disease with increased power compared to individual analyses.

To address potential issues regarding meta-analysis application in gene expression studies, including laboratory effects, variations between probes and differential platforms, we planned this analysis beginning with data selection, through pre-processing and filtering, and finally to the meta-analysis model. This pipeline included the adjustment for study-related batch effects to minimize superficial inter-study discordance caused mainly by random noise and technical differences. Although in general this adjustment risk may confound true biological differences, we found no differences in disease severity, IgE and age across the four cohorts and thus feel confident that true biological differences are unlikely to be confounded. The resulting list of DEGs is presented here as a Meta-Analysis Derived Atopic Dermatitis (MADAD) transcriptome, utilizing 4 separate microarray studies for a total of 97 samples. It is our hope that the MADAD will prove helpful to investigators who may benefit from our robust characterization of the AD phenotype.

Similar to the MAD psoriasis transcriptome, our meta-analysis approach resulted in a more concise number of DEGs than that found in the individual studies [25]. However, this set of DEGs provides a more biologically relevant and powerful AD phenotype compared with previous studies [12], since important inflammatory and barrier pathways are more significantly enriched compared to individual studies, as is often observed with meta-analysis derived transcriptomes [91].

Our transcriptome is the first association of AD genomic fingerprinting with the atherosclerosis signaling pathway, which includes genes associated with vascular inflammation (SELE, IL-37, S100A8) [63, 77]. SELE has been independently associated with coronary heart disease and carotid artery atherosclerosis, and its expression in the vascular endothelium of the dermis of AD patients has also been observed [51, 88]. This supports the emerging view of AD as a systemic disease, which, like psoriasis, extends far beyond the skin [92, 93].

We have recently shown increased systemic immune activation in blood of moderate-to-severe AD patients among both skin homing and systemic T-cell subsets [94]. Furthermore, when comparing blood moderate-to-severe psoriasis and AD patients, we have shown that AD is associated with systemic activation and increased polar differentiation of effector and memory T-cell subsets, with higher and persistent activation within skin homing subsets. (Czarnowicki et al.-In press) AD patients also demonstrated higher levels of ICOS activation in circulating skin-homing subsets, (Czarnowicki- In press) consistent with the significant overrepresentation of ICOS signaling in our IPA analysis.

In large cohort studies, AD has also been recently associated with a variety of systemic diseases including inflammatory bowel disease [95, 96], Type 1 Diabetes Mellitus [97], and ADHD [98], providing further evidence for its systemic nature. Additionally, AD, like psoriasis, was recently shown to be associated with increased vascular inflammation using CT imaging [99101]. The association of AD with systemic involvement emphasizes the great unmet need among severe AD patients for systemic therapeutic approaches, which are now in clinical trials (NCT00638989) [10, 18].

An interesting association between barrier defects and vascular inflammation in AD may be represented by the gene SPTLC2, which encodes a SPT, the rate-limiting enzyme in de novo synthesis of sphingomyelin and ceramides, and which has also been shown to influence atherosclerosis. Plasma sphingomyelin level was found to be an independent risk factor for coronary artery disease and is associated with subclinical atherosclerosis in humans [102, 103]. In ApoE knockout mice, inhibition of SPT resulted in improved lipid profiles and prevented the development of atherosclerotic lesions [104]. Sphingomyelin is proposed to affect atherosclerosis by influencing lipid metabolism and regulating cell proliferation and apoptosis to modulate plaque growth and stability [105]. SPTLC2 increases with epidermal barrier abnormalities [75], as was observed in our MADAD transcriptome, and may be related to the previously unacknowledged systemic vascular inflammation in AD.

Lipid and differentiation abnormalities represent the hallmark of defective barrier function in AD [106, 107]. The Th2 cytokine effects on inhibition of terminal differentiation genes (e.g. filaggrin, loricrin) have been well documented by in vitro [1, 78] and in vivo studies [15, 108]. While the effects of Th2 cytokines on suppression of epidermal lipids have become recently available in model systems [78], our paper is the first to show negative correlations between expression of Th2 cytokines and epidermal lipids. It is established that AD LS skin shows alterations in lipid composition of the stratum corneum [109], with decreases in long-chain ceramides and free fatty acids in addition to disorganized lipid structure [107]. While Th1 cytokines, including TNFα and IFNγ, have been shown to induce ceramide synthesis, Th2 cytokines (including IL-4) inhibit the production of ceramides necessary for proper barrier function [76, 110]. Our findings support these in vitro models, and show significant negative correlations between increased Th2 cytokine production and decreased lipid expression.

We also identified for the first time as differentially expressed in AD several genes involved in lipid metabolism. These include FA2H, encoding protein essential to the de novo synthesis of specific ceramides that are critical in maintaining the permeability barrier of the epidermis [111, 112], and ELOVL3, encoding a protein involved in the elongation of long chain fatty acids and essential in prevention of transepidermal water loss [113, 114]. While not involved in fatty acid synthesis, aquaporin 5/AQP5 is a water channel that was found to be decreased in LS AD skin [115117].

The network analysis of MADAD DEGs shows a significant correlation between AD disease activity/SCORAD and the S. aureus module in NL skin, emphasizing the role of colonization and infection in the onset of immune activation in background skin. S. aureus colonization/infection has been shown to occur at significantly higher rates in both LS and NL AD skin compared to healthy controls [118120], and was shown to induce Th2 and Th22 immune activation [84, 121123]. In LS skin, viral and innate immune process modules were significantly positively correlated with SCORAD, emphasizing the association between reactions to external pathogens and active inflammation in AD [121].

The treatment effect on the 19 discriminatory genes highlights the diverse mechanistic effects of these therapies with several observations worth mentioning. First, MMP3, a marker of general inflammation, displays an impressive recovery with the targeted therapy dupilumab compared to the more nonspecific immune suppressant treatments. This inflammatory gene improved over 300 % in only 4 weeks of treatment with dupilumab, while long-term treatment with CsA and UVB did not reach 100 % recovery. This may suggest the ability of dupilumab to suppress immune dysregulation in AD in a shorter timeframe than with less specific agents. Another noteworthy gene is Selectin E, encoding a protein involved in leukocyte extravasation, which also showed higher levels of recovery in dupilumab compared to CsA and UVB therapy. The difference in these markers may be related to the mechanism of action of each drug; UVB has direct effects on keratinocytes and thus mainly mediates signals originating there, while CsA is a nonspecific inhibitor of T-cells, B-cells and related pathways. Dupilumab more specifically modulates IL-4/IL-13 signaling, which has been implicated in the pathogenesis of AD, and this may account for the impressive recovery seen in genes central to the inflammatory response generated by AD.


This meta-analysis provides a stable and robust AD transcriptome worthy of future biomarker selection and evaluation of treatment response. The value of a standardized AD transcriptome will only increase with the bench-to-bedside translational approach, as a standardized measure of treatment response that cannot be confounded with placebo effects. This meta-analysis provides an integrative model of AD, emphasizing both immune and barrier abnormalities and also highlighting the systemic nature of the disease.



Atopic dermatitis


Cyclosporin A


Differentially expressed genes (logFCH ≥ 1 and fdr ≤ 0.05)


False Discovery Rate


Fixed-effects model (special case of REM)


Integrated Correlation Coefficient


Integration-driven discovery


Immunoglobulin E


Ingenuity pathway analysis




Meta-analysis derived AD


Meta Threshold Gradient Directed Regularization




Preferred Reporting Items for Systematic Reviews and Meta-Analyses


Random-effects model


Quantitative Reverse Transcription PCR


Ribosomal protein, large, P0 (formerly termed human acidic ribosomal protein P0/hARP)

S. aureus:

Staphylococcus aureus


Scoring of Atopic Dermatitis


Serine palmitoyltransferase


Weighted gene co-expression network analysis


  1. 1.

    Boguniewicz M, Leung DYM. Atopic dermatitis: a disease of altered skin barrier and immune dysregulation. Immunol Rev. 2011;242(1):233–46.

    CAS  PubMed  PubMed Central  Google Scholar 

  2. 2.

    Leung DYM. New insights into atopic dermatitis: role of skin barrier and immune dysregulation. Allergol Int. 2013;62(2):151–61.

    CAS  PubMed  Google Scholar 

  3. 3.

    Dhingra N, Suarez-Farinas M, Fuentes-Duculan J, Gittler JK, Shemer A, Raz A, et al. Attenuated neutrophil axis in atopic dermatitis compared to psoriasis reflects TH17 pathway differences between these diseases. J Allergy Clin Immunol. 2013;132(2):498–501.e3.

    CAS  PubMed  Google Scholar 

  4. 4.

    Malajian D, & Guttman-Yassky E. New pathogenic and therapeutic paradigms in atopic dermatitis. Cytokine. 2015;73(2):311-318.

  5. 5.

    N & Guttman-Yassky E. The translational revolution and use of biologics in patients with inflammatory skin diseases. Journal of Allergy and Clinical Immunology. 2015;135(2):324-336.

  6. 6.

    Olsson M, Broberg A, Jernås M, Carlsson L, Rudemo M, Suurküla M, et al. Increased expression of aquaporin 3 in atopic eczema. Allergy. 2006;61(9):1132–7.

    CAS  PubMed  Google Scholar 

  7. 7.

    Guttman-Yassky E, Suárez-Fariñas M, Chiricozzi A, Nograles KE, Shemer A, Fuentes-Duculan J, et al. Broad defects in epidermal cornification in atopic dermatitis identified through genomic analysis. J Allergy Clin Immunol. 2009;124(6):1235–1244.e58.

    CAS  PubMed  Google Scholar 

  8. 8.

    Elias PM, Hatano Y, Williams ML. Basis for the barrier abnormality in atopic dermatitis: Outside-inside-outside pathogenic mechanisms. J Allergy Clin Immunol. 2008;121(6):1337–43.

    CAS  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Gittler JK, Shemer A, Suárez-Fariñas M, Fuentes-Duculan J, Gulewicz KJ, Wang CQF, et al. Progressive activation of T(H)2/T(H)22 cytokines and selective epidermal proteins characterizes acute and chronic atopic dermatitis. J Allergy Clin Immunol. 2012;130(6):1344–54.

    CAS  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Beck L, Thaçi D, Hamilton JD, Graham NM, Bieber T, Rocklin R, et al. Dupilumab Treatment in Adults with Moderate-to-Severe Atopic Dermatitis. N Engl J Med. 2014;371(2):130–9.

    PubMed  Google Scholar 

  11. 11.

    Tintle S, Shemer A, Suárez-Fariñas M, Fujita H, Gilleaudeau P, Sullivan-Whalen M, et al. Reversal of atopic dermatitis with narrow-band UVB phototherapy and biomarkers for therapeutic response. J Allergy Clin Immunol. 2011;128(3):583–93.e1–4.

    CAS  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Khattri S, Shemer A, Rozenblit M, Dhingra N, Czarnowicki T, Finney R, & Guttman-Yassky E. Cyclosporine in patients with atopic dermatitis modulates activated inflammatory pathways and reverses epidermal pathology. Journal of Allergy and Clinical Immunology. 2014;133(6):1626-634.

  13. 13.

    Jensen JM, Pfeiffer S, Witt M, Bräutigam M, Neumann C, Weichenthal M, et al. Different effects of pimecrolimus and betamethasone on the skin barrier in patients with atopic dermatitis. J Allergy Clin Immunol. 2009;123:1124–33.

    CAS  PubMed  Google Scholar 

  14. 14.

    Plager A, Leontovich AA, Henke SA, Davis MDP, McEvoy MT, Sciallis GF, et al. Early cutaneous gene transcription changes in adult atopic dermatitis and potential clinical implications. Exp Dermatol. 2007;16(1):28–36.

    CAS  PubMed  Google Scholar 

  15. 15.

    Suárez-Fariñas M, Tintle SJ, Shemer A, Chiricozzi A, Nograles K, Cardinale I, et al. Nonlesional atopic dermatitis skin is characterized by broad terminal differentiation defects and variable immune abnormalities. J Allergy Clin Immunol. 2011;127(4):954–64.e1–4.

    PubMed  PubMed Central  Google Scholar 

  16. 16.

    Zaba LC, Suárez-Fariñas M, Fuentes-Duculan J, Nograles KE, Guttman-Yassky E, Cardinale I, et al. Effective treatment of psoriasis with etanercept is linked to suppression of IL-17 signaling, not immediate response TNF genes. J Allergy Clin Immunol. 2009;124(5):1022–10.e1–395.

    CAS  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Krueger JG, Fretzin S, Suárez-Fariñas M, Haslett PA, Phipps KM, Cameron GS, et al. IL-17A is essential for cell activation and inflammatory gene circuits in subjects with psoriasis. J Allergy Clin Immunol. 2012;130(1):1–18.

    Google Scholar 

  18. 18.

    Hamilton JD, Suárez-Fariñas M, Dhingra N, Cardinale I, Li X, Kostic A, et al. Dupilumab improves the molecular signature in skin of patients with moderate-to-severe atopic dermatitis. J Allergy Clin Immunol. 2014;134(6):1293–300.

    CAS  PubMed  Google Scholar 

  19. 19.

    Jiang K, Frank MB, Chen Y, Osban J, Jarvis JN. Genomic characterization of remission in juvenile idiopathic arthritis. Arthritis Res Ther. 2013;15(4):R100.

    PubMed  PubMed Central  Google Scholar 

  20. 20.

    Cahan P, Rovegno F, Mooney D, Newman JC, St Laurent G, McCaffrey TA. Meta-analysis of microarray results: challenges, opportunities, and recommendations for standardization. Gene. 2007;401(1–2):12–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  21. 21.

    Ramasamy A, Mondry A, Holmes CC, Altman DG. Key issues in conducting a meta-analysis of gene expression microarray datasets. PLoS Med. 2008;5(9):e184.

    PubMed  PubMed Central  Google Scholar 

  22. 22.

    Zakharkin SO, Kim K, Mehta T, Chen L, Barnes S, Scheirer KE, et al. Sources of variation in Affymetrix microarray experiments. BMC Bioinformatics. 2005;6:214.

    PubMed  PubMed Central  Google Scholar 

  23. 23.

    Suárez-Fariñas M, Noggle S, Heke M, Hemmati-Brivanlou A, Magnasco MO. Comparing independent microarray studies: the case of human embryonic stem cells. BMC Genomics. 2005;6:99.

    PubMed  PubMed Central  Google Scholar 

  24. 24.

    Suárez-Fariñas M, Lowes MA, Zaba LC, Krueger JG. Evaluation of the psoriasis transcriptome across different studies by gene set enrichment analysis (GSEA). PLoS One. 2010;5(4):e10247.

    PubMed  PubMed Central  Google Scholar 

  25. 25.

    Tian S, Krueger JG, Li K, Jabbari A, Brodmerkel C, Lowes MA, et al. Meta-analysis derived (MAD) transcriptome of psoriasis defines the ‘core’ pathogenesis of disease. PLoS One. 2012;7(9):e44274.

    CAS  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Chang L-C, Lin H-M, Sibille E, Tseng GC. Meta-analysis methods for combining multiple expression profiles: comparisons, statistical characterization and an application guideline. BMC Bioinformatics. 2013;14(vi):368.

    PubMed  PubMed Central  Google Scholar 

  27. 27.

    Choi JK, Yu U, Kim S, Yoo OJ. Combining multiple microarray studies and modeling interstudy variation. Bioinformatics. 2003;19 Suppl 1:i84–90.

    PubMed  Google Scholar 

  28. 28.

    Tseng GC, Ghosh D, Feingold E. Comprehensive literature review and statistical considerations for microarray meta-analysis. Nucleic Acids Res. 2012;40(9):3785–99.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Liberati A, Altman DG, Tetzlaff J, Mulrow C, Gøtzsche PC, Ioannidis JPA, et al. The PRISMA statement for reporting systematic reviews and meta-analyses of studies that evaluate health care interventions: explanation and elaboration. PLoS Med. 2009;6(7):e1000100.

    PubMed  PubMed Central  Google Scholar 

  30. 30.

    Team R, Core. "R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria, 2012." (2014).

  31. 31.

    Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5(10):R80.

    PubMed  PubMed Central  Google Scholar 

  32. 32.

    Wu Z, Irizarry RA, Gentleman R, Martinez-Murillo F, Spencer F. A Model-Based Background Adjustment for Oligonucleotide Expression Arrays. J Am Stat Assoc. 2004;99(468):909–17.

    Google Scholar 

  33. 33.

    Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007;8(1):118–27.

    PubMed  Google Scholar 

  34. 34.

    Leek JT, Johnson WE, Parker HS, Jaffe AE, Storey JD. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. 2012;28(6):882–3.

    CAS  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Leek JT, Scharpf RB, Bravo HC, Simcha D, Langmead B, Johnson WE, et al. Tackling the widespread and critical impact of batch effects in high-throughput data. Nat Rev Genet. 2010;11:733–9.

    CAS  PubMed  Google Scholar 

  36. 36.

    Parmigiani G. A Cross-Study Comparison of Gene Expression Studies for the Molecular Classification of Lung Cancer. Clin Cancer Res. 2004;10(9):2922–7.

    CAS  PubMed  Google Scholar 

  37. 37.

    Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc. 1995;57(1):289–300.

    Google Scholar 

  38. 38.

    Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.

    CAS  Google Scholar 

  39. 39.

    Suárez-Fariñas M, Arbeit R, Jiang W, Ortenzio FS, Sullivan T, Krueger JG. Suppression of molecular inflammatory pathways by Toll-like receptor 7, 8, and 9 antagonists in a model of IL-23-induced skin inflammation. PLoS One. 2013;8(12):e84634.

    PubMed  PubMed Central  Google Scholar 

  40. 40.

    Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.

    PubMed  PubMed Central  Google Scholar 

  41. 41.

    Ma S, Huang J. Regularized gene selection in cancer microarray meta-analysis. BMC Bioinformatics. 2009;10:1.

    CAS  PubMed  PubMed Central  Google Scholar 

  42. 42.

    Tian S, Suarez-Farinas M. Multi-TGDR: a regularization method for multi-class classification in microarray experiments. PLoS One. 2013;8(11):e78302.

    PubMed  PubMed Central  Google Scholar 

  43. 43.

    Mitsui H, Suárez-Fariñas M, Belkin DA, Levenkova N, Fuentes-Duculan J, Coats I, et al. Combined use of laser capture microdissection and cDNA microarray analysis identifies locally expressed disease-related genes in focal regions of psoriasis vulgaris skin lesions. J Invest Dermatol. 2012;132(6):1615–26.

    CAS  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Rappaport N, Nativ N, Stelzer G, Twik M, Guan-Golan Y, Stein TI, et al. MalaCards: an integrated compendium for diseases and their annotation. Database (Oxford). 2013;2013:bat018.

    Google Scholar 

  45. 45.

    Appace. Solr. 2014.

    Google Scholar 

  46. 46.

    Wittkowski KM, Song T. Nonparametric methods for molecular biology. Methods Mol Biol. 2010;620:105–53.

    CAS  PubMed  PubMed Central  Google Scholar 

  47. 47.

    Irizarry RA, Warren D, Spencer F, Kim IF, Biswal S, Frank BC, et al. Multiple-laboratory comparison of microarray platforms. Nat Methods. 2005;2(5):345–50.

    CAS  PubMed  Google Scholar 

  48. 48.

    Lowes MA, Suárez-Fariñas M, Krueger JG. Immunology of psoriasis. Annu Rev Immunol. 2014;32:227–55.

    CAS  PubMed  PubMed Central  Google Scholar 

  49. 49.

    Panzer R, Blobel C, Fölster-Holst R, Proksch E. TLR2 and TLR4 expression in atopic dermatitis, contact dermatitis and psoriasis. Exp Dermatol. 2014;23(5):364–6.

    CAS  PubMed  Google Scholar 

  50. 50.

    Kuo I-H, Carpenter-Mendini A, Yoshida T, McGirt LY, Ivanov AI, Barnes KC, et al. Activation of epidermal toll-like receptor 2 enhances tight junction function: implications for atopic dermatitis and skin barrier repair. J Invest Dermatol. 2013;133(4):988–98.

    CAS  PubMed  Google Scholar 

  51. 51.

    Hwang SJ, Ballantyne CM, Sharrett AR, Smith LC, Davis CE, Gotto AM, et al. Circulating adhesion molecules VCAM-1, ICAM-1, and E-selectin in carotid atherosclerosis and incident coronary heart disease cases: the Atherosclerosis Risk In Communities (ARIC) study. Circulation. 1997;96(12):4219–25.

    CAS  PubMed  Google Scholar 

  52. 52.

    Kiechl S, Paré G, Barbalic M, Qi L, Dupuis J, Dehghan A, et al. Association of variation at the ABO locus with circulating levels of soluble intercellular adhesion molecule-1, soluble P-selectin, and soluble E-selectin: a meta-analysis. Circ Cardiovasc Genet. 2011;4(6):681–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  53. 53.

    Mead JR, Ramji DP. The pivotal role of lipoprotein lipase in atherosclerosis. Cardiovasc Res. 2002;55(2):261–9.

    CAS  PubMed  Google Scholar 

  54. 54.

    Chen Q, Razzaghi H, Demirci FY, Kamboh MI. Functional significance of lipoprotein lipase HindIII polymorphism associated with the risk of coronary artery disease. Atherosclerosis. 2008;200(1):102–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  55. 55.

    Goldberg IJ. Lipoprotein lipase and lipolysis: central roles in lipoprotein metabolism and atherogenesis. J Lipid Res. 1996;37(4):693–707.

    CAS  PubMed  Google Scholar 

  56. 56.

    Xu X, Wang L, Xu C, Zhang P, Yong F, Liu H, et al. Variations in matrix metalloproteinase-1, −3, and −9 genes and the risk of acute coronary syndrome and coronary artery disease in the Chinese Han population. Coron Artery Dis. 2013;24(4):259–65.

    PubMed  Google Scholar 

  57. 57.

    Abilleira S, Bevan S, Markus HS. The role of genetic variants of matrix metalloproteinases in coronary and carotid atherosclerosis. J Med Genet. 2006;43(12):897–901.

    CAS  PubMed  PubMed Central  Google Scholar 

  58. 58.

    Averill MM, Kerkhoff C, Bornfeldt KE. S100A8 and S100A9 in cardiovascular biology and disease. Arterioscler Thromb Vasc Biol. 2012;32(2):223–9.

    CAS  PubMed  Google Scholar 

  59. 59.

    Jin S, Park CO, Shin JU, Noh JY, Lee YS, Lee NR, et al. DAMP molecules S100A9 and S100A8 activated by IL-17A and house-dust mites are increased in atopic dermatitis. Exp Dermatol. 2014;23(12):938–41.

    CAS  PubMed  Google Scholar 

  60. 60.

    Ji Q, Zeng Q, Huang Y, Shi Y, Lin Y, Lu Z, et al. Elevated plasma IL-37, IL-18, and IL-18BP concentrations in patients with acute coronary syndrome. Mediators Inflamm. 2014;2014:165742.

    PubMed  PubMed Central  Google Scholar 

  61. 61.

    Inouye M, Ripatti S, Kettunen J, Lyytikäinen L-P, Oksala N, Laurila P-P, et al. Novel Loci for metabolic networks and multi-tissue expression studies reveal genes for atherosclerosis. PLoS Genet. 2012;8(8):e1002907.

    CAS  PubMed  PubMed Central  Google Scholar 

  62. 62.

    Tregouet DA, Barbaux S, Poirier O, Blankenberg S, Bickel C, Escolano S, et al. SELPLG gene polymorphisms in relation to plasma SELPLG levels and coronary artery disease. Ann Hum Genet. 2003;67(Pt 6):504–11.

    CAS  PubMed  Google Scholar 

  63. 63.

    Esaki H, Ewald DA, Ungar B, Rozenblit M. Identification of novel immune and barrier genes in atopic dermatitis by means of laser capture microdissection. J Allergy Clin Immunol. 2015;135(1):153–63.

    CAS  PubMed  PubMed Central  Google Scholar 

  64. 64.

    Choi UY, Kang JS, Hwang YS, Kim YJ. Oligoadenylate synthase-like (OASL) proteins: dual functions and associations with diseases. Experimental & molecular medicine. 2015;47(3):e144.

  65. 65.

    Wensink AC, Hack CE, Bovenschen N. Granzymes Regulate Proinflammatory Cytokine Responses. J Immunol. 2015;194(2):491–7.

    CAS  PubMed  Google Scholar 

  66. 66.

    Westerberg R, Mansson JE, Golozoubova V, Shabalina IG, Backlund EC, Tvrdik P, et al. ELOVL3 is an important component for early onset of lipid recruitment in brown adipose tissue. J Biol Chem. 2006;281(8):4958–68.

    CAS  PubMed  Google Scholar 

  67. 67.

    Turlo K, Leung CS, Seo JJ, Goulbourne CN, Adeyo O, Gin P, et al. Equivalent binding of wild-type lipoprotein lipase (LPL) and S447X-LPL to GPIHBP1, the endothelial cell LPL transporter. Biochim Biophys Acta. 2014;1841(7):963–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  68. 68.

    Goulbourne CN, Gin P, Tatar A, Nobumori C, Hoenger A, Jiang H, et al. The GPIHBP1-LPL complex is responsible for the margination of triglyceride-rich lipoproteins in capillaries. Cell Metab. 2014;19(5):849–60.

    CAS  PubMed  PubMed Central  Google Scholar 

  69. 69.

    Pan Y-S, Lee Y-S, Lee Y-L, Lee W-C, Hsieh S-Y. Differentially profiling the low-expression transcriptomes of human hepatoma using a novel SSH/microarray approach. BMC Genomics. 2006;7:131.

    PubMed  PubMed Central  Google Scholar 

  70. 70.

    Tortola L, Rosenwald E, Abel B, Blumberg H, Schäfer M, Coyle AJ, et al. Psoriasiform dermatitis is driven by IL-36-mediated DC-keratinocyte crosstalk. J Clin Invest. 2012;122(11):3965–76.

    CAS  PubMed  PubMed Central  Google Scholar 

  71. 71.

    Nold MF, Nold-Petry CA, Zepp JA, Palmer BE, Bufler P, Dinarello CA. IL-37 is a fundamental inhibitor of innate immunity. Nat Immunol. 2010;11(11):1014–22.

    CAS  PubMed  PubMed Central  Google Scholar 

  72. 72.

    da Rosa JC, Malajian D, Shemer A, Rozenblit M, Dhingra N, Czarnowicki T, et al. Patients with atopic dermatitis have attenuated and distinct contact hypersensitivity responses to common allergens in skin. Journal of Allergy and Clinical Immunology, 2015:135(3):712-20.

  73. 73.

    van Smeden J, Janssens M, Kaye ECJ, Caspers PJ, Lavrijsen AP, Vreeken RJ, et al. The importance of free fatty acid chain length for the skin barrier function in atopic eczema patients. Exp Dermatol. 2014;23:45–52.

    PubMed  Google Scholar 

  74. 74.

    Lee Y, Je Y-J, Lee S-S, Li ZJ, Choi D-K, Kwon Y-B, et al. Changes in transepidermal water loss and skin hydration according to expression of aquaporin-3 in psoriasis. Ann Dermatol. 2012;24(2):168–74.

    CAS  PubMed  PubMed Central  Google Scholar 

  75. 75.

    Stachowitz S, Alessandrini F, Abeck D, Ring J, Behrendt H. Permeability barrier disruption increases the level of serine palmitoyltransferase in human epidermis. J Invest Dermatol. 2002;119(5):1048–52.

    CAS  PubMed  Google Scholar 

  76. 76.

    Sawada E, Yoshida N, Sugiura A, Imokawa G. Th1 cytokines accentuate but Th2 cytokines attenuate ceramide production in the stratum corneum of human epidermal equivalents: an implication for the disrupted barrier mechanism in atopic dermatitis. J Dermatol Sci. 2012;68(1):25–35.

    CAS  PubMed  Google Scholar 

  77. 77.

    Tanno O, Ota Y, Kitamura N, Katsube T, Inoue S. Nicotinamide increases biosynthesis of ceramides as well as other stratum corneum lipids to improve the epidermal permeability barrier. Br J Dermatol. 2000;143(3):524–31.

    CAS  PubMed  Google Scholar 

  78. 78.

    Danso MO, van Drongelen V, Mulder A, van Esch J, Scott H, van Smeden J, et al. TNF-α and Th2 Cytokines Induce Atopic Dermatitis-Like Features on Epidermal Differentiation Proteins and Stratum Corneum Lipids in Human Skin Equivalents. J Invest Dermatol. 2014;134(7):1941–50.

    CAS  PubMed  Google Scholar 

  79. 79.

    Hirota T, Saeki H, Tomita K, Tanaka S, Ebe K, Sakashita M, et al. Variants of C-C motif chemokine 22 (CCL22) are associated with susceptibility to atopic dermatitis: case–control studies. PLoS One. 2011;6(11):e26987.

    CAS  PubMed  PubMed Central  Google Scholar 

  80. 80.

    Hussein YM, Shalaby SM, Nassar A, Alzahrani SS, Alharbi AS, Nouh M. Association between genes encoding components of the IL-4/IL-4 receptor pathway and dermatitis in children. Gene. 2014;545(2):276–81.

    CAS  PubMed  Google Scholar 

  81. 81.

    Gao P-S, Rafaels NM, Mu D, Hand T, Murray T, Boguniewicz M, et al. Genetic variants in thymic stromal lymphopoietin are associated with atopic dermatitis and eczema herpeticum. J Allergy Clin Immunol. 2010;125(6):1403–1407.e4.

    CAS  PubMed  PubMed Central  Google Scholar 

  82. 82.

    Langfelder P, Horvath S. WGCNA package FAQ. 2014 [Online]. Available: [Accessed: 03-Mar-2015].

    Google Scholar 

  83. 83.

    Hon KL, Ching GK, Ng PC, Leung TF. Exploring CCL18, eczema severity and atopy. Pediatr Allergy Immunol. 2011;22(7):704–7.

    PubMed  Google Scholar 

  84. 84.

    Niebuhr M, Mamerow D, Heratizadeh A, Satzger I, Werfel T. Staphylococcal alpha-toxin induces a higher T cell proliferation and interleukin-31 in atopic dermatitis. Int Arch Allergy Immunol. 2011;156(4):412–5.

    CAS  PubMed  Google Scholar 

  85. 85.

    Guttman-Yassky E, Dhingra N, Leung DYM. New era of biologic therapeutics in atopic dermatitis. Expert Opin Biol Ther. 2013;13(4):549–61.

    CAS  PubMed  Google Scholar 

  86. 86.

    Rozenblit M, Suarez-Farinas M, Shemer A, Khattri S, Gilleaudeau P, Sullivan-Whalen M, et al. Residual genomic profile after cyclosporine treatment may offer insights into atopic dermatitis reoccurrence. J Allergy Clin Immunol. 2014;134(4):955–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  87. 87.

    Brandt C, Pavlovic V, Radbruch A, Worm M, Baumgrass R. Low-dose cyclosporine A therapy increases the regulatory T cell population in patients with atopic dermatitis. Allergy. 2009;64(11):1588–96.

    CAS  PubMed  Google Scholar 

  88. 88.

    Marinović Kulišić S, Lipozenčić J, Gregorović G, Gregurek G, Lacković G, Jurić-Lekić G. Expression of e-selectin in the skin of patients with atopic dermatitis: morphometric study. Acta Dermatovenerol Croat. 2010;18(3):141–5.

    PubMed  Google Scholar 

  89. 89.

    Tsao C-H, Shiau M-Y, Chuang P-H, Chang Y-H, Hwang J. Interleukin-4 regulates lipid metabolism by inhibiting adipogenesis and promoting lipolysis. J Lipid Res. 2014;55(3):385–97.

    CAS  PubMed  PubMed Central  Google Scholar 

  90. 90.

    Morgan SA, McCabe EL, Gathercole LL, Hassan-Smith ZK, Larner DP, Bujalska IJ, et al. 11β-HSD1 is the major regulator of the tissue-specific effects of circulating glucocorticoid excess. Proc Natl Acad Sci U S A. 2014;111(24):E2482–91.

    CAS  PubMed  PubMed Central  Google Scholar 

  91. 91.

    Pérez-Montarelo D, Hudson NJ, Fernández AI, Ramayo-Caldas Y, Dalrymple BP, Reverter A. Porcine tissue-specific regulatory networks derived from meta-analysis of the transcriptome. PLoS One. 2012;7(9):e46159.

    PubMed  PubMed Central  Google Scholar 

  92. 92.

    Voiculescu VM, Lupu M, Papagheorghe L, Giurcaneanu C, Micu E. Psoriasis and Metabolic Syndrome–scientific evidence and therapeutic implications. Journal of medicine and life. 2014;7(4):468.

  93. 93.

    Oliveira Mde F, Rocha Bde O, Duarte GV. Psoriasis: classical and emerging comorbidities. An Bras Dermatol. 2015;90(1):9–20.

    PubMed  Google Scholar 

  94. 94.

    Czarnowicki T, Gonzalez J, Shemer A, Malajian D, Xu H, Zheng X, et al. Severe atopic dermatitis is characterized by selective expansion of circulating TH2/TC2 and TH22/TC22, but not TH17/TC17, cells within the skin-homing T-cell population. J Allergy Clin Immunol. 2015.

  95. 95.

    Niwa Y, Sumi H, Akamatsu H. An association between ulcerative colitis and atopic dermatitis, diseases of impaired superficial barriers. J Invest Dermatol. 2004;123(5):999–1000.

    CAS  PubMed  Google Scholar 

  96. 96.

    Myrelid P, Dufmats M, Lilja I, Grinn C, Lannerstad O, Sjödahl R. Atopic manifestations are more common in patients with Crohn disease than in the general population. Scand J Gastroenterol. 2004;39(8):731–6.

    CAS  PubMed  Google Scholar 

  97. 97.

    Cardwell CR, Shields MD, Carson DJ, Patterson CC. A meta-analysis of the association between childhood type 1 diabetes and atopic disease. Diabetes Care. 2003;26(9):2568–74.

    PubMed  Google Scholar 

  98. 98.

    Schmitt J, Buske-Kirschbaum A, Roessner V. Is atopic disease a risk factor for attention-deficit/hyperactivity disorder? A systematic review. Allergy. 2010;65(12):1506–24.

    CAS  PubMed  Google Scholar 

  99. 99.

    7 th International Congress of Psoriasis: from Gene to Clinic The Queen Elizabeth II Conference Centre, London, U.K. 11th-13th December 2014. Br. J. Dermatol. 2014;171:6, e105–e173.

  100. 100.

    Dave J, Ahlman MA, Lockshin BN, Bluemke DA, Mehta NN. Vascular inflammation in psoriasis localizes to the arterial wall using a novel imaging technique. J Am Acad Dermatol. 2014;70(6):1137–8.

    PubMed  PubMed Central  Google Scholar 

  101. 101.

    Rose S, Sheth NH, Baker JF, Ogdie A, Raper A, Saboury B, et al. A comparison of vascular inflammation in psoriasis, rheumatoid arthritis, and healthy subjects by FDG-PET/CT: a pilot study. Am J Cardiovasc Dis. 2013;3(4):273–8.

    PubMed  PubMed Central  Google Scholar 

  102. 102.

    Jiang XC, Paultre F, Pearson TA, Reed RG, Francis CK, Lin M, et al. Plasma sphingomyelin level as a risk factor for coronary artery disease. Arterioscler Thromb Vasc Biol. 2000;20(12):2614–8.

    CAS  PubMed  Google Scholar 

  103. 103.

    Nelson JC, Jiang X-C, Tabas I, Tall A, Shea S. Plasma sphingomyelin and subclinical atherosclerosis: findings from the multi-ethnic study of atherosclerosis. Am J Epidemiol. 2006;163(10):903–12.

    PubMed  Google Scholar 

  104. 104.

    Park T-S, Panek RL, Mueller SB, Hanselman JC, Rosebury WS, Robertson AW, et al. Inhibition of sphingomyelin synthesis reduces atherogenesis in apolipoprotein E-knockout mice. Circulation. 2004;110(22):3465–71.

    CAS  PubMed  Google Scholar 

  105. 105.

    Tabas I. Sphingolipids and atherosclerosis: a mechanistic connection? A therapeutic opportunity? Circulation. 2004;110(22):3400–1.

    PubMed  Google Scholar 

  106. 106.

    van Smeden J, Janssens M, Gooris GS, Bouwstra JA. The important role of stratum corneum lipids for the cutaneous barrier function. Biochim Biophys Acta. 2014;1841(3):295–313.

    PubMed  Google Scholar 

  107. 107.

    Janssens M, van Smeden J, Gooris GS, Bras W, Portale G, Caspers PJ, et al. Increase in short-chain ceramides correlates with an altered lipid organization and decreased barrier function in atopic eczema patients. J Lipid Res. 2012;53(12):2755–66.

    CAS  PubMed  PubMed Central  Google Scholar 

  108. 108.

    Fujita H, Shemer A, Suárez-Fariñas M, Johnson-Huang LM, Tintle S, Cardinale I, et al. Lesional dendritic cells in patients with chronic atopic dermatitis and psoriasis exhibit parallel ability to activate T-cell subsets. J Allergy Clin Immunol. 2011;128(3):574–82.e1–12.

    CAS  PubMed  Google Scholar 

  109. 109.

    Sääf AM, Tengvall-Linder M, Chang HY, Adler AS, Wahlgren C-FF, Scheynius A, et al. Global expression profiling in atopic eczema reveals reciprocal expression of inflammatory and lipid genes. PLoS One. 2008;3(12):e4017.

    PubMed  PubMed Central  Google Scholar 

  110. 110.

    Hatano Y, Terashi H, Arakawa S, Katagiri K. Interleukin-4 suppresses the enhancement of ceramide synthesis and cutaneous permeability barrier functions induced by tumor necrosis factor-alpha and interferon-gamma in human epidermis. J Invest Dermatol. 2005;124(4):786–92.

    CAS  PubMed  Google Scholar 

  111. 111.

    Hama H. Fatty acid 2-Hydroxylation in mammalian sphingolipid biology. Biochim Biophys Acta. 2010;1801(4):405–14.

    CAS  PubMed  Google Scholar 

  112. 112.

    Pappas A. Epidermal surface lipids. Dermatoendocrinol. 2009;1(2):72–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  113. 113.

    Westerberg R, Tvrdik P, Unden AB, Mansson JE, Norlen L, Jakobsson A, et al. Role for ELOVL3 and fatty acid chain length in development of hair and skin function. J Biol Chem. 2004;279(7):5621–9.

    CAS  PubMed  Google Scholar 

  114. 114.

    Fehrenschild D, Galli U, Breiden B, Bloch W, Schettina P, Brodesser S, et al. TCF/Lef1-mediated control of lipid metabolism regulates skin barrier function. J Invest Dermatol. 2012;132(2):337–45.

    CAS  PubMed  Google Scholar 

  115. 115.

    Blaydon DC, Lind LK, Plagnol V, Linton KJ, Smith FJ, Wilson NJ, et al. Mutations in AQP5, encoding a water-channel protein, cause autosomal-dominant diffuse nonepidermolytic palmoplantar keratoderma. Am J Hum Genet. 2013;93(2):330–5.

    CAS  PubMed  PubMed Central  Google Scholar 

  116. 116.

    Boury-Jamot M, Daraspe J, Bonte F, Perrier E, Schnebert S, Dumas M, et al. Skin aquaporins: function in hydration, wound healing, and skin epidermis homeostasis. Handb Exp Pharmacol. 2009;190:205–17.

    CAS  PubMed  Google Scholar 

  117. 117.

    Nejsum LN, Kwon TH, Jensen UB, Fumagalli O, Frokiaer J, Krane CM, et al. Functional requirement of aquaporin-5 in plasma membranes of sweat glands. Proc Natl Acad Sci U S A. 2002;99(1):511–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  118. 118.

    Breuer K, HAussler S, Kapp A, Werfel T. Staphylococcus aureus: colonizing features and influence of an antibacterial treatment in adults with atopic dermatitis. Br J Dermatol. 2002;147(1):55–61.

    CAS  PubMed  Google Scholar 

  119. 119.

    Higaki S, Morohashi M, Yamagishi T, Hasegawa Y. Comparative study of staphylococci from the skin of atopic dermatitis patients and from healthy subjects. Int J Dermatol. 1999;38(4):265–9.

    CAS  PubMed  Google Scholar 

  120. 120.

    Park H-Y, Kim C-R, Huh I-S, Jung M-Y, Seo E-Y, Park J-H, et al. Staphylococcus aureus Colonization in Acute and Chronic Skin Lesions of Patients with Atopic Dermatitis. Ann Dermatol. 2013;25(4):410–6.

    PubMed  PubMed Central  Google Scholar 

  121. 121.

    Matsui K, Nishikawa A. Peptidoglycan from Staphylococcus aureus induces T(H)2 immune response in mice. J Investig Allergol Clin Immunol. 2012;22(2):80–6.

    CAS  PubMed  Google Scholar 

  122. 122.

    Warner JA, McGirt LY, Beck LA. Biomarkers of Th2 polarity are predictive of staphylococcal colonization in subjects with atopic dermatitis. Br J Dermatol. 2009;160(1):183–5.

    CAS  PubMed  Google Scholar 

  123. 123.

    Niebuhr M, Scharonow H, Gathmann M, Mamerow D, Werfel T. Staphylococcal exotoxins are strong inducers of IL-22: A potential role in atopic dermatitis. J Allergy Clin Immunol. 2010;126(6):1176–83.e4.

    CAS  PubMed  Google Scholar 

Download references


The authors thank Joel Correa da Rosa, Sandra Garcet, and Troels Marstrand for helpful discussions on the analysis pipeline. This work was supported by LEO Pharma AS and the Danish Ministry of Science and Higher Education.


This work was supported by the danish ministry of scienceand higher education in terms of partly funding of a PhD fellowship (

Author information



Corresponding author

Correspondence to Mayte Suárez-Fariñas.

Additional information

Competing interests

DAE is partly funded by the Danish Ministry of Science and Higher Education and LEO Pharma AS. TL is employed by LEO Pharma AS. All other authors have no competing interests to disclose.

Authors’ contributions

DAE carried out the main statistical and bioinformatics analysis, and prepared the manuscript. DM prepared the manuscript and carried out biological assessment and evaluations of the statistical results. JGK participated in the design of the study and manuscript preparation. CTW advised and participated in the bioinformatics analysis. TW and ST carried out the classification analysis. TL participated in the design of the study and manuscript preparation. EGY participated in the design of the study and directed manuscript preparation. MSF designed the study, directed manuscript preparation, and directed the statistical analysis. All authors have read and approved the final manuscript.

Authors’ information

Not applicable.

Availability of data and materials

Not applicable.

Additional files

Additional file 1:

Supplementary materials and methods. (DOCX 111 kb)

Additional file 2: Table E9.

Primer IDs for RT-PCR analysis. (XLSX 37 kb)

Additional file 3: Figure E1.

A) Prisma. B) MADAD Datasets. C) Meta-Analysis Workflow. (PDF 608 kb)

Additional file 4: Table E1.

All datasets intially considered before applying inclusion/exclusion criteria. (DOCX 65 kb)

Additional file 5: Table E2.

Consensus of re-analyized included datasets (by Probes). (XLSX 60 kb)

Additional file 6: Figure E2.

QQ plot. Comparison of parameters. (PDF 1005 kb)

Additional file 7: Table E3.

ALL MADAD DEGs (FCH ≥ 2, FDR ≤ 0.05). (XLSX 171 kb)

Additional file 8: Table E4.

MADAD – Ingenuity Pathway Analysis. (XLSX 33 kb)

Additional file 9: Table E5.

Khattri 2014 (GSE58558) – Ingenuity Pathway Analysis. (XLSX 26 kb)

Additional file 10: Figure E3.

A) Inflammatory Skin Disease Gene Sets. B) MADAD IDD-DEGs IPA Canonical Pathways. (PDF 85 kb)

Additional file 11: Table E6.

IPA Pathways and GO terms IDD genes. (XLSX 9 kb)

Additional file 12: Table E7.

Overrepresented Hub30 Genes and module names. (LS) (XLSX 38 kb)

Additional file 13: Table E8.

Overrepresented Hub30 Genes and module names. (NL) (XLSX 47 kb)

Additional file 14: Figure E4.

logFCHs correlations of MADAD and treatments A) Dupilumab (300mg), B) Cyclosporin (W12), and C) UVB. (PDF 42 kb)

Additional file 15: Table E10.

Adjusted P-values for the significance indicators in Figure 2B. (XLSX 8 kb)

Additional file 16: Table E11.

BH adjusted p-values WGCNA LS Modules. (XLSX 9 kb)

Additional file 17: Table E12.

BH adjusted p-values WGCNA NL Modules. (XLSX 9 kb)

Additional file 18: Table E13.

MTGDR results. (XLSX 9 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Ewald, D.A., Malajian, D., Krueger, J.G. et al. Meta-analysis derived atopic dermatitis (MADAD) transcriptome defines a robust AD signature highlighting the involvement of atherosclerosis and lipid metabolism pathways. BMC Med Genomics 8, 60 (2015).

Download citation


  • Atopic dermatitis
  • Meta-analysis
  • Transcriptome
  • Atherosclerosis
  • Expression analysis