Search

As the circadian clock regulates fundamental biological processes, disrupted clocks are often observed in patients and diseased tissues. Determining the circadian time of the patient or the tissue of focus is essential in circadian medicine and research. Here we present tauFisher, a computational pipeline that accurately predicts circadian time from a single transcriptomic sample by finding correlations between rhythmic genes within the sample. We demonstrate tauFisher's performance in adding timestamps to both bulk and single-cell transcriptomic samples collected from multiple tissue types and experimental settings. Application of tauFisher at a cell-type level in a single-cell RNAseq dataset collected from mouse dermal skin implies that greater circadian phase heterogeneity may explain the dampened rhythm of collective core clock gene expression in dermal immune cells compared to dermal fibroblasts. Given its robustness and generalizability across assay platforms, experimental setups, and tissue types, as well as its potential application in single-cell RNAseq data analysis, tauFisher is a promising tool that facilitates circadian medicine and research.

Single-cell technologies have become essential to driving discovery in both basic and translational investigative dermatology. Despite the multitude of available datasets, a central reference atlas of normal human skin, which can serve as a reference resource for skin cell types, cell states, and their molecular signatures, is still lacking. For any such atlas to receive broad acceptance, participation by many investigators during atlas construction is an essential prerequisite. As part of the Human Cell Atlas project, we have assembled a Skin Biological Network to build a consensus Human Skin Cell Atlas and outline a roadmap toward that goal. We define the drivers of skin diversity to be considered when selecting sequencing datasets for the atlas and list practical hurdles during skin sampling that can result in data gaps and impede comprehensive representation and technical considerations for tissue processing and computational analysis, the accounting for which should minimize biases in cell type enrichments and exclusions and decrease batch effects. By outlining our goals for Atlas 1.0, we discuss how it will uncover new aspects of skin biology.

The immunopathogenesis of psoriasis, a common chronic inflammatory disease of the skin, is incompletely understood. Here we demonstrate, using a combination of single cell and spatial RNA sequencing, IL-36 dependent amplification of IL-17A and TNF inflammatory responses in the absence of neutrophil proteases, which primarily occur within the supraspinous layer of the psoriatic epidermis. We further show that a subset of SFRP2+ fibroblasts in psoriasis contribute to amplification of the immune network through transition to a pro-inflammatory state. The SFRP2+ fibroblast communication network involves production of CCL13, CCL19 and CXCL12, connected by ligand-receptor interactions to other spatially proximate cell types: CCR2+ myeloid cells, CCR7+ LAMP3+ dendritic cells, and CXCR4 expressed on both CD8+ Tc17 cells and keratinocytes, respectively. The SFRP2+ fibroblasts also express cathepsin S, further amplifying inflammatory responses by activating IL-36G in keratinocytes. These data provide an in-depth view of psoriasis pathogenesis, which expands our understanding of the critical cellular participants to include inflammatory fibroblasts and their cellular interactions.

Palmoplantar skin is structurally and functionally unique, but the transcriptional programs driving this specialization are unclear. Here, we use bulk and single-cell RNA sequencing of human palm, sole, and hip skin to describe the distinguishing characteristics of palmoplantar and non-palmoplantar skin while also uncovering differences between palmar and plantar sites. Our approach reveals an altered immune environment in palmoplantar skin, with downregulation of diverse immunological processes and decreased immune cell populations. Further, we identify specific fibroblast populations that appear to orchestrate key differences in cell-cell communication in palm, sole, and hip. Dedicated keratinocyte analysis highlights major differences in basal cell fraction among the three sites and demonstrates the existence of two spinous keratinocyte populations constituting parallel, site-selective epidermal differentiation trajectories. In summary, this deep characterization of highly adapted palmoplantar skin contributes key insights into the fundamental biology of human skin and provides a valuable data resource for further investigation.

Psoriasis is an immune-mediated inflammatory skin disease typically characterized by erythematous and scaly plaques. It affects 3% of the Newfoundland population while only affecting 1.7% of the general Canadian population. Recent genome-wide association studies (GWAS) in psoriasis have identified more than 63 genetic susceptibility loci that individually have modest effects. Prior studies have shown that a genetic risk score (GRS) combining multiple loci can improve psoriasis disease prediction. However, these prior GRS studies have not fully explored the association of GRS with patient clinical characteristics. In this study, we calculated three types of GRS: one using all known GWAS SNPs (GRS-ALL), one using a subset of SNPs from the HLA region (GRS-HLA), and the last using non-HLA SNPs (GRS-noHLA). We examined the relationship between these GRS and a number of psoriasis features within a well characterized Newfoundland psoriasis cohort. We found that both GRS-ALL and GRS-HLA were significantly associated with early age of psoriasis onset, psoriasis severity, first presentation of psoriasis at the elbow or knee, and the total number of body locations affected, while only GRS-ALL was associated with a positive family history of psoriasis. GRS-noHLA was uniquely associated with genital psoriasis. These findings clarify the relationship of the HLA and non-HLA components of GRS with important clinical features of psoriasis.

The epidermis is the outermost layer of skin. Here, we used targeted lipid profiling to characterize the biogeographic alterations of human epidermal lipids across 12 anatomically distinct body sites, and we used single-cell RNA-Seq to compare keratinocyte gene expression at acral and nonacral sites. We demonstrate that acral skin has low expression of EOS acyl-ceramides and the genes involved in their synthesis, as well as low expression of genes involved in filaggrin and keratin citrullination (PADI1 and PADI3) and corneodesmosome degradation, changes that are consistent with increased corneocyte retention. Several overarching principles governing epidermal lipid expression were also noted. For example, there was a strong negative correlation between the expression of 18-carbon and 22-carbon sphingoid base ceramides. Disease-specific alterations in epidermal lipid gene expression and their corresponding alterations to the epidermal lipidome were characterized. Lipid biomarkers with diagnostic utility for inflammatory and precancerous conditions were identified, and a 2-analyte diagnostic model of psoriasis was constructed using a step-forward algorithm. Finally, gene coexpression analysis revealed a strong connection between lipid and immune gene expression. This work highlights (a) mechanisms by which the epidermis is uniquely adapted for the specific environmental insults encountered at different body surfaces and (b) how inflammation-associated alterations in gene expression affect the epidermal lipidome.

Proprotein convertase subtilisin/kexin type-9 (PCSK9) is a posttranslational regulator of the LDL receptor (LDLR). Recent studies have proposed a role for PCSK9 in regulating immune responses. Using RNA-Seq-based variant discovery, we identified a possible psoriasis-susceptibility locus at 1p32.3, located within PCSK9 (rs662145 C > T). This finding was verified in independently acquired genomic and RNA-Seq data sets. Single-cell RNA-Seq (scRNA-Seq) identified keratinocytes as the primary source of PCSK9 in human skin. PCSK9 expression, however, was not uniform across keratinocyte subpopulations. scRNA-Seq and IHC demonstrated an epidermal gradient of PCSK9, with expression being highest in basal and early spinous layer keratinocytes and lowest in granular layer keratinocytes. IL36G expression followed the opposite pattern, with expression highest in granular layer keratinocytes. PCSK9 siRNA knockdown experiments confirmed this inverse relationship between PCSK9 and IL36G expression. Other immune genes were also linked to PCSK9 expression, including IL27RA, IL1RL1, ISG20, and STX3. In both cultured keratinocytes and nonlesional human skin, homozygosity for PCSK9 SNP rs662145 C > T was associated with lower PCSK9 expression and higher IL36G expression, when compared with heterozygous skin or cell lines. Together, these results support PCSK9 as a psoriasis-susceptibility locus and establish a putative link between PCSK9 and inflammatory cytokine expression.

Acne affects 1 in 10 people globally, often resulting in disfigurement. The disease involves excess production of lipids, particularly squalene, increased growth of Cutibacterium acnes, and a host inflammatory response with foamy macrophages. By combining single-cell and spatial RNA sequencing as well as ultrahigh-resolution Seq-Scope analyses of early acne lesions on back skin, we identified TREM2 macrophages expressing lipid metabolism and proinflammatory gene programs in proximity to hair follicle epithelium expressing squalene epoxidase. We established that the addition of squalene induced differentiation of TREM2 macrophages in vitro, which were unable to kill C. acnes. The addition of squalene to macrophages inhibited induction of oxidative enzymes and scavenged oxygen free radicals, providing an explanation for the efficacy of topical benzoyl peroxide in the clinical treatment of acne. The present work has elucidated the mechanisms by which TREM2 macrophages and unsaturated lipids, similar to their involvement in atherosclerosis, may contribute to the pathogenesis of acne.

Tape stripping is a minimally invasive, nonscarring method that can be utilized to assess gene expression in the skin but is infrequently used given technical constraints. By comparing different tape stripping technologies and full-thickness skin biopsy results of lesional and nonlesional psoriatic skin from the same patients, we demonstrate that tape stripping with optimized high-resolution transcriptomic profiling can be used to effectively assess and characterize inflammatory responses in the skin. Upon comparison with single-cell RNA-sequencing data from psoriatic full-thickness skin biopsies, we illustrate that tape-stripping efficiently captures the transcriptome of the upper layers of the epidermis with sufficient resolution to assess the molecular components of the feed-forward immune amplification pathway in psoriasis. Notably, nonlesional psoriatic skin sampled by tape stripping demonstrates activated, proinflammatory changes when compared to healthy control skin, suggesting a prepsoriatic state, which is not captured on full-thickness skin biopsy transcriptome profiling. This work illustrates an approach to assess inflammatory response in the epidermis by combining noninvasive sampling with high throughput RNA-sequencing, providing a foundation for biomarker discoveries and mechanism of action studies for inflammatory skin conditions.

Cutaneous lupus erythematosus (CLE) is a disfiguring and poorly understood condition frequently associated with systemic lupus. Previous studies suggest that nonlesional keratinocytes play a role in disease predisposition, but this has not been investigated in a comprehensive manner or in the context of other cell populations. To investigate CLE immunopathogenesis, normal-appearing skin, lesional skin, and circulating immune cells from lupus patients were analyzed via integrated single-cell RNA sequencing and spatial RNA sequencing. We demonstrate that normal-appearing skin of patients with lupus represents a type I interferon-rich, prelesional environment that skews gene transcription in all major skin cell types and markedly distorts predicted cell-cell communication networks. We also show that lupus-enriched CD16+ dendritic cells undergo robust interferon education in the skin, thereby gaining proinflammatory phenotypes. Together, our data provide a comprehensive characterization of lesional and nonlesional skin in lupus and suggest a role for skin education of CD16+ dendritic cells in CLE pathogenesis.

BackgroundA major issue with the current management of psoriasis is our inability to predict treatment response.ObjectiveOur aim was to evaluate the ability to use baseline molecular expression profiling to assess treatment outcome for patients with psoriasis.MethodsWe conducted a longitudinal study of 46 patients with chronic plaque psoriasis treated with anti-TNF agent etanercept, and molecular profiles were assessed in more than 200 RNA-seq samples.ResultsWe demonstrated correlation between clinical response and molecular changes during the course of the treatment, particularly for genes responding to IL-17A/TNF in keratinocytes. Intriguingly, baseline gene expressions in nonlesional, but not lesional, skin were the best marker of treatment response at week 12. We identified USP18, a known regulator of IFN responses, as positively correlated with Psoriasis Area and Severity Index (PASI) improvement (P = 9.8 × 10-4) and demonstrate its role in regulating IFN/TNF responses in keratinocytes. Consistently, cytokine gene signatures enriched in baseline nonlesional skin expression profiles had strong correlations with PASI improvement. Using this information, we developed a statistical model for predicting PASI75 (ie, 75% of PASI improvement) at week 12, achieving area under the receiver-operating characteristic curve value of 0.75 and up to 80% accurate PASI75 prediction among the top predicted responders.ConclusionsOur results illustrate feasibility of assessing drug response in psoriasis using nonlesional skin and implicate involvement of IFN regulators in anti-TNF responses.

Altered epidermal differentiation along with increased keratinocyte proliferation is a characteristic feature of psoriasis and pityriasis rubra pilaris (PRP). However, despite this large degree of overlapping clinical and histologic features, the molecular signatures these skin disorders share are unknown. Using global transcriptomic profiling, we demonstrate that plaque psoriasis and PRP skin lesions have high overlap, with all differentially expressed genes in PRP relative to normal skin having complete overlap with those in psoriasis. The major common pathway shared between psoriasis and PRP involves the phospholipases PLA2G2F, PLA2G4D, and PLA2G4E, which were found to be primarily expressed in the epidermis. Gene silencing each of the 3 PLA2s led to reduction in immune responses and epidermal thickness both in vitro and in vivo in a mouse model of psoriasis, establishing their proinflammatory roles. Lipidomic analyses demonstrated that PLA2s affect mobilization of a phospholipid-eicosanoid pool, which is altered in psoriatic lesions and functions to promote immune responses in keratinocytes. Taken together, our results highlight the important role of PLA2s as regulators of epidermal barrier homeostasis and inflammation, identify PLA2s as a shared pathogenic mechanism between PRP and psoriasis, and as potential therapeutic targets for both diseases.

Many inflammatory skin diseases are characterized by altered epidermal differentiation. Whether this altered differentiation promotes inflammatory responses has been unknown. Here, we show that IRAK2, a member of the signaling complex downstream of IL-1 and IL-36, correlates positively with disease severity in both atopic dermatitis and psoriasis. Inhibition of epidermal IRAK2 normalizes differentiation and inflammation in two mouse models of psoriasis- and atopic dermatitis-like inflammation. Specifically, we demonstrate that IRAK2 ties together proinflammatory and differentiation-dependent responses and show that this function of IRAK2 is specific to keratinocytes and acts through the differentiation-associated transcription factor ZNF750. Taken together, our findings suggest that IRAK2 has a critical role in promoting feed-forward amplification of inflammatory responses in skin through modulation of differentiation pathways and inflammatory responses.

Granulomas are complex cellular structures composed predominantly of macrophages and lymphocytes that function to contain and kill invading pathogens. Here, we investigated the single-cell phenotypes associated with antimicrobial responses in human leprosy granulomas by applying single-cell and spatial sequencing to leprosy biopsy specimens. We focused on reversal reactions (RRs), a dynamic process whereby some patients with disseminated lepromatous leprosy (L-lep) transition toward self-limiting tuberculoid leprosy (T-lep), mounting effective antimicrobial responses. We identified a set of genes encoding proteins involved in antimicrobial responses that are differentially expressed in RR versus L-lep lesions and regulated by interferon-γ and interleukin-1β. By integrating the spatial coordinates of the key cell types and antimicrobial gene expression in RR and T-lep lesions, we constructed a map revealing the organized architecture of granulomas depicting compositional and functional layers by which macrophages, T cells, keratinocytes and fibroblasts can each contribute to the antimicrobial response.

Hidradenitis suppurativa (HS) is a chronic inflammatory disease characterized by abscesses, nodules, dissecting/ draining tunnels, and extensive fibrosis. Here, we integrate single-cell RNA sequencing, spatial transcriptomics, and immunostaining to provide an unprecedented view of the pathogenesis of chronic HS, characterizing the main cellular players and defining their interactions. We found a striking layering of the chronic HS infiltrate and identified the contribution of 2 fibroblast subtypes ([SFRP4.sup.+] and [CXCL13.sup.+]) in orchestrating this compartmentalized immune response. We further demonstrated the central role of the Hippo pathway in promoting extensive fibrosis in HS and provided preclinical evidence that the profibrotic fibroblast response in HS can be modulated through inhibition of this pathway. These data provide insights into key aspects of HS pathogenesis with broad therapeutic implications., Introduction Hidradenitis suppurativa (HS) is a chronic inflammatory disease of the skin that affects 1% of the general population (1). The disease is characterized by acute, recurrent inflammatory nodules and [...]

Psoriasis is a chronic inflammatory skin condition that has a fairly wide range of clinical presentations. Plaque psoriasis, which is the most common manifestation of psoriasis, is located on one end of the spectrum, dominated by adaptive immune responses, whereas the rarer pustular psoriasis lies on the opposite end, dominated by innate and autoinflammatory immune responses. In recent years, genetic studies have identified six genetic variants that predispose to pustular psoriasis, and these have highlighted the role of IL-36 cytokines as central to pustular psoriasis pathogenesis. In this review, we discuss the presentation and clinical subtypes of pustular psoriasis, contribution of genetic predisposing variants, critical role of the IL-36 family of cytokines in disease pathophysiology, and treatment perspectives for pustular psoriasis. We further outline the application of appropriate mouse models for the study of pustular psoriasis and address the outstanding questions and issues related to our understanding of the mechanisms involved in pustular psoriasis.

Killer cell immunoglobulin-like receptors (KIR) regulate immune responses in NK and CD8+ T cells via interaction with HLA ligands. KIR genes, including KIR2DS1, KIR3DL1, and KIR3DS1 have previously been implicated in psoriasis susceptibility. However, these previous studies were constrained to small sample sizes, in part due to the time and expense required for direct genotyping of KIR genes. Here, we implemented KIR*IMP to impute KIR copy number from single-nucleotide polymorphisms (SNPs) on chromosome 19 in the discovery cohort (n=11,912) from the PAGE consortium, University of California San Francisco, and the University of Dundee, and in a replication cohort (n=66,357) from Kaiser Permanente Northern California. Stratified multivariate logistic regression that accounted for patient ancestry and high-risk HLA alleles revealed that KIR2DL2 copy number was significantly associated with psoriasis in the discovery cohort (p ≤ 0.05). The KIR2DL2 copy number association was replicated in the Kaiser Permanente replication cohort. This is the first reported association of KIR2DL2 copy number with psoriasis and highlights the importance of KIR genetics in the pathogenesis of psoriasis.

The skin serves as the primary interface between our body and the external environment and acts as a barrier against entry of physical agents, chemicals, and microbes. Keratinocytes make up the main cellular constitute of the outermost layer of the skin, contributing to the formation of the epidermis, and they are crucial for maintaining the integrity of this barrier. Beyond serving as a physical barrier component, keratinocytes actively participate in maintaining tissue homeostasis, shaping, amplifying, and regulating immune responses in skin. Keratinocytes act as sentinels, continuously monitoring changes in the environment, and, through microbial sensing, stretch, or other physical stimuli, can initiate a broad range of inflammatory responses via secretion of various cytokines, chemokines, and growth factors. This diverse function of keratinocytes contributes to the highly variable clinical manifestation of skin immune responses. In this Review, we highlight the highly diverse functions of epidermal keratinocytes and their contribution to various immune-mediated skin diseases.

Hidradenitis suppurativa (HS) is a debilitating chronic inflammatory skin disease characterized by chronic abscess formation and development of multiple draining sinus tracts in the groin, axillae, and perineum. Using proteomic and transcriptomic approaches, we characterized the inflammatory responses in HS in depth, revealing immune responses centered on IFN-γ, IL-36, and TNF, with lesser contribution from IL-17A. We further identified B cells and plasma cells, with associated increases in immunoglobulin production and complement activation, as pivotal players in HS pathogenesis, with Bruton's tyrosine kinase (BTK) and spleen tyrosine kinase (SYK) pathway activation as a central signal transduction network in HS. These data provide preclinical evidence to accelerate the path toward clinical trials targeting BTK and SYK signaling in moderate-to-severe HS.

BackgroundAlthough multiple studies have assessed molecular changes in chronic atopic dermatitis (AD) lesions, little is known about the transition from acute to chronic disease stages, and the factors and mechanisms that shape chronic inflammatory activity.ObjectivesWe sought to assess the global transcriptome changes that characterize the progression from acute to chronic stages of AD.MethodsWe analyzed transcriptome changes in paired nonlesional skin, acute and chronic AD lesions from 11 patients and 38 healthy controls by RNA-sequencing, and conducted in vivo and histological assays to evaluate findings.ResultsOur data demonstrate that approximately 74% of the genes dysregulated in acute lesions remain or are further dysregulated in chronic lesions, whereas only 34% of the genes dysregulated in chronic lesions are altered already in the acute stage. Nonlesional AD skin exhibited enrichment of TNF, TH1, TH2, and TH17 response genes. Acute lesions showed marked dendritic-cell signatures and a prominent enrichment of TH1, TH2, and TH17 responses, along with increased IL-36 and thymic stromal lymphopoietin expression, which were further heightened in chronic lesions. In addition, genes involved in skin barrier repair, keratinocyte proliferation, wound healing, and negative regulation of T-cell activation showed a significant dysregulation in the chronic versus acute comparison. Furthermore, our data show progressive changes in vasculature and maturation of dendritic-cell subsets with chronicity, with FOXK1 acting as immune regulator.ConclusionsOur results show that the changes accompanying the transition from nonlesional to acute to chronic inflammation in AD are quantitative rather than qualitative, with chronic AD having heightened TH2, TH1, TH17, and IL36 responses and skin barrier repair mechanisms. These findings provide novel insights and highlight underappreciated pathways in AD pathogenesis that may be amenable to therapeutic targeting.

The circadian clock coordinates a variety of immune responses with signals from the external environment to promote survival. We investigated the potential reciprocal relationship between the circadian clock and skin inflammation. We treated mice topically with the Toll-like receptor 7 (TLR7) agonist imiquimod (IMQ) to activate IFN-sensitive gene (ISG) pathways and induce psoriasiform inflammation. IMQ transiently altered core clock gene expression, an effect mirrored in human patient psoriatic lesions. In mouse skin 1 d after IMQ treatment, ISGs, including the key ISG transcription factor IFN regulatory factor 7 (Irf7), were more highly induced after treatment during the day than the night. Nuclear localization of phosphorylated-IRF7 was most prominently time-of-day dependent in epidermal leukocytes, suggesting that these cell types play an important role in the diurnal ISG response to IMQ. Mice lacking Bmal1 systemically had exacerbated and arrhythmic ISG/Irf7 expression after IMQ. Furthermore, daytime-restricted feeding, which affects the phase of the skin circadian clock, reverses the diurnal rhythm of IMQ-induced ISG expression in the skin. These results suggest a role for the circadian clock, driven by BMAL1, as a negative regulator of the ISG response, and highlight the finding that feeding time can modulate the skin immune response. Since the IFN response is essential for the antiviral and antitumor effects of TLR activation, these findings are consistent with the time-of-day-dependent variability in the ability to fight microbial pathogens and tumor initiation and offer support for the use of chronotherapy for their treatment.

Oil-secreting sebaceous glands (SGs) are critical for proper skin function; however, it remains unclear how different factors act together to modulate SG stem cells. Here, we provide functional evidence that each SG lobe is serviced by its own dedicated stem cell population. Upon ablating Notch signaling in different skin subcompartments, we find that this pathway exerts dual counteracting effects on SGs. Suppressing Notch in SG progenitors traps them in a hybrid state where stem and differentiation features become intermingled. In contrast, ablating Notch outside of the SG stem cell compartment indirectly drives SG expansion. Finally, we report that a K14:K5→K14:K79 keratin shift occurs during SG differentiation. Deleting K79 destabilizes K14 in sebocytes, and attenuates SGs and eyelid meibomian glands, leading to corneal ulceration. Altogether, our findings demonstrate that SGs integrate diverse signals from different niches and suggest that mutations incurred within one stem cell compartment can indirectly influence another.

Atopic dermatitis (AD) affects up to 20% of children and adults worldwide. To gain a deeper understanding of the pathophysiology of AD, we conducted a large-scale transcriptomic study of AD with deeply sequenced RNA-sequencing samples using long (126-bp) paired-end reads. In addition to the comparisons against previous transcriptomic studies, we conducted in-depth analysis to obtain a high-resolution view of the global architecture of the AD transcriptome and contrasted it with that of psoriasis from the same cohort. By using 147 RNA samples in total, we found striking correlation between dysregulated genes in lesional psoriasis and lesional AD skin with 81% of AD dysregulated genes being shared with psoriasis. However, we described disease-specific molecular and cellular features, with AD skin showing dominance of IL-13 pathways, but with near undetectable IL-4 expression. We also demonstrated greater disease heterogeneity and larger proportion of dysregulated long noncoding RNAs in AD, and illustrated the translational impact, including skin-type classification and drug-target prediction. This study is by far the largest study comparing the AD and psoriasis transcriptomes using RNA sequencing and demonstrating the shared inflammatory components, as well as specific discordant cytokine signatures of these two skin diseases.

The present paradigm of psoriasis pathogenesis revolves around the IL-23/IL-17A axis. Dual-secreting Th17 T cells presumably are the predominant sources of the psoriasis phenotype-driving cytokines, IL-17A and IL-22. We thus conducted a meta-analysis of independently acquired RNA-seq psoriasis datasets to explore the relationship between the expression of IL17A and IL22. This analysis failed to support the existence of dual secreting IL-17A/IL-22 Th17 cells as a major source of these cytokines. However, variable relationships amongst the expression of psoriasis susceptibility genes and of IL17A, IL22, and IL23A were identified. Additionally, to shed light on gene expression relationships in psoriasis, we applied a machine learning nonlinear dimensionality reduction strategy (t-SNE) to display the entire psoriasis transcriptome as a 2-dimensonal image. This analysis revealed a variety of gene clusters, relevant to psoriasis pathophysiology but failed to support a relationship between IL17A and IL22. These results support existing theories on alternative sources of IL-17A and IL-22 in psoriasis such as a Th22 cells and non-T cell populations.