Psoriasis is a common immune mediated disorder that affects the skin, nails, and joints. To identify psoriasis susceptibility loci, we genotyped 438,670 SNPs in 1,409 European ancestry psoriasis cases and 1,436 controls. Twenty-one promising SNPs were followed-up in 5,048 psoriasis cases and 5,041 controls. Our results provide strong support for the association of at least seven genetic loci and psoriasis (each with p < 5×10−8 overall). Loci with confirmed association encode HLA-C, three genes involved in IL-23 signaling (IL23A, IL23R, IL12B), two genes that act downstream of TNF-α and regulate NF-κB signaling (TNIP1, TNFAIP3), and two genes involved in the modulation of Th2 immune responses (IL4, IL13). Although the proteins encoded in these loci are known to interact biologically, we found no evidence for epistasis between associated SNPs. Our results expand the catalog of genetic loci implicated in psoriasis susceptibility and suggest priority targets for study in other auto-immune disorders.

Summary To gain further insight into the genetic architecture of psoriasis, we conducted a meta-analysis of three genome-wide association studies (GWAS) and two independent datasets genotyped on the Immunochip, involving 10,588 cases and 22,806 controls in total. We identified 15 new disease susceptibility regions, increasing the number of psoriasis-associated loci to 36 for Caucasians. Conditional analyses identified five independent signals within previously known loci. The newly identified shared disease regions encompassed a number of genes whose products regulate T-cell function (e.g. RUNX3, TAGAP and STAT3). The new psoriasis-specific regions were notable for candidate genes whose products are involved in innate host defense, encoding proteins with roles in interferon-mediated antiviral responses (DDX58), macrophage activation (ZC3H12C), and NF-κB signaling (CARD14 and CARM1). These results portend a better understanding of shared and distinctive genetic determinants of immune-mediated inflammatory disorders and emphasize the importance of the skin in innate and acquired host defense.

To increase our understanding of psoriasis, we utilized RNA-seq to assay the transcriptomes of lesional psoriatic and normal skin. We sequenced polyadenylated RNA-derived cDNAs from 92 psoriatic and 82 normal punch biopsies, generating an average of ~38 million single-end 80-bp reads per sample. Comparison of 42 samples examined by both RNA-seq and microarray revealed marked differences in sensitivity, with transcripts identified only by RNA-seq having much lower expression than those also identified by microarray. RNA-seq identified many more differentially expressed transcripts enriched in immune system processes. Weighted gene co-expression network analysis (WGCNA) revealed multiple modules of coordinately expressed epidermal differentiation genes, overlapping significantly with genes regulated by the long non-coding RNA TINCR, its target gene, staufen-1 (STAU1), the p63 target gene ZNF750, and its target KLF4. Other coordinately expressed modules were enriched for lymphoid and/or myeloid signature transcripts and genes induced by IL-17 in keratinocytes. Dermally-expressed genes were significantly down-regulated in psoriatic biopsies, most likely due to expansion of the epidermal compartment. These results demonstrate the power of WGCNA to elucidate gene regulatory circuits in psoriasis, and emphasize the influence of tissue architecture in both differential expression and co-expression analysis.

Psoriasis is a complex disease of skin with a prevalence of about 2%. We conducted the largest meta-analysis of genome-wide association studies (GWAS) for psoriasis to date, including data from eight different Caucasian cohorts, with a combined effective sample size >39,000 individuals. We identified 16 additional psoriasis susceptibility loci achieving genome-wide significance, increasing the number of identified loci to 63 for European-origin individuals. Functional analysis highlighted the roles of interferon signalling and the NFκB cascade, and we showed that the psoriasis signals are enriched in regulatory elements from different T cells (CD8+ T-cells and CD4+ T-cells including TH0, TH1 and TH17). The identified loci explain ∼28% of the genetic heritability and generate a discriminatory genetic risk score (AUC=0.76 in our sample) that is significantly correlated with age at onset (p=2 × 10−89). This study provides a comprehensive layout for the genetic architecture of common variants for psoriasis.

Psoriasis is a multifactorial skin disease characterized by epidermal hyperproliferation and chronic inflammation, the most common form of which is psoriasis vulgaris (PsV). We present a genome-wide association analysis of 2,339,118 SNPs in 472 psoriasis patients and 1,146 controls from Germany, with follow-up of the 147 most significant SNPs in 2,746 PsV cases and 4,140 controls from three independent replication panels. We identified an association at TRAF3IP2 on 6q21 and genotyped two SNPs at this locus in two additional replication panels (the combined discovery and replication panels consisted of 6,487 cases and 8,037 controls; combinded P = 2.36×10−10 for rs13210247 and combined P = 1.24×10−16 for rs33980500). About 15% of psoriasis cases develope psoriatic arthritis (PsA). A stratified analysis of our datasets including only PsA cases (1,922 cases compared to 8,037 controls, P=4.57×10−12 for rs33980500) suggested that TRAF3IP2 represents a shared susceptibility for PsV and PsA. TRAF3IP2 encodes a protein involved in IL-17 signaling and which interacts with memebers of the Rel/NF-κB transcription factor family.

We extended our previous GWAS for psoriasis with a a multistage replication study including 8,312 cases and 12,919 controls from China as well as 3,293 cases, 4,188 controls from Germany and the USA, and 254 nuclear families from the USA. We identified 6 new susceptibility loci associated to psoriasis in Chinese, containing candidate genes ERAP1, PTTG1, CSMD1, GJB2, SERPINB8, ZNF816A (PCombined<5×10−8) and replicated one locus 5q33.1 (TNIP1/ANXA6) previously reported (PCombined=3.8×10−21) in European studies. Two of these loci showed evidence for association evidence in the German study, at ZNF816A and GJB2 with P=3.6×10−3 and P=7.9×10−3, respectively. ERAP1 and ZNF816A were preferentially associated with Type I (early onset) psoriasis in Chinese Han population (test for heterogeneity P=6.5×10−3 and P=1.5×10−3, respectively). Comparisons with previous GWAS of psoriasis highlight the heterogeneity of disease susceptibility between Chinese and European populations. Our study identifies new genetic susceptibility factors and suggests new biological pathways in psoriasis.

To further elucidate molecular alterations in psoriasis, we performed a gene expression study on 58 paired lesional and uninvolved psoriatic and 64 control skin samples. Comparison of involved psoriatic (PP) and normal (NN) skin identified 1,326 differentially regulated transcripts encoding 918 unique genes (548 up- and 369 down-regulated), of which over 600 were to our knowledge unreported, including S100A7A, THRSP and ELOVL3. Strongly up-regulated genes included SERPINB4, PI3, DEFB4 and several S100 family members. Strongly down-regulated genes included Wnt-inhibitory factor-1 (WIF1), betacellulin (BTC), and CCL27. Enriched gene ontology categories included immune response, defense response, and keratinocyte differentiation. Biological processes regulating fatty acid and lipid metabolism were enriched in the down-regulated gene set. Comparison of the psoriatic transcriptome to the transcriptomes of cytokine-stimulated cultured keratinocytes (IL-17, IL-22, IL-1α, IFN-γ, TNF-α and oncostatin-M) revealed surprisingly little overlap, with the cytokine stimulated keratinocyte expression representing only 2.5%, 0.7%, 1.5%, 5.6%, 5.0% and 1.9% of the lesional psoriatic dysregulated transcriptome, respectively. This comprehensive analysis of differentially regulated transcripts in psoriasis provides additional insight into the pathogenic mechanisms involved and emphasizes the need for more complex yet tractable experimental models of psoriasis.

Psoriasis (PS) and Crohn disease (CD) have been shown to be epidemiologically, pathologically, and therapeutically connected, but little is known about their shared genetic causes. We performed meta-analyses of five published genome-wide association studies on PS (2,529 cases and 4,955 controls) and CD (2,142 cases and 5,505 controls), followed up 20 loci that showed strongest evidence for shared disease association and, furthermore, tested cross-disease associations for previously reported PS and CD risk alleles in additional 6,115 PS cases, 4,073 CD cases, and 10,100 controls. We identified seven susceptibility loci outside the human leukocyte antigen region (9p24 near JAK2, 10q22 at ZMIZ1, 11q13 near PRDX5, 16p13 near SOCS1, 17q21 at STAT3, 19p13 near FUT2, and 22q11 at YDJC) shared between PS and CD with genome-wide significance (p < 5 3 10 À8 ) and confirmed four already established PS and CD risk loci (IL23R, IL12B, REL, and TYK2). Three of the shared loci are also genome-wide significantly associated with PS alone (10q22 at ZMIZ1, p rs1250544 ¼ 3.53 3 10 À8 , 11q13 near PRDX5, p rs694739 ¼ 3.71 3 10 À09 , 22q11 at YDJC, p rs181359 ¼ 8.02 3 10 À10 ). In addition, we identified one susceptibility locus for CD (16p13 near SOCS1, p rs4780355 ¼ 4.99 3 10 À8 ). Refinement of association signals identified shared genomewide significant associations for exonic SNPs at 10q22 (ZMIZ1) and in silico expression quantitative trait locus analyses revealed that the associations at ZMIZ1 and near SOCS1 have a potential functional effect on gene expression. Our results show the usefulness of joint analyses of clinically distinct immune-mediated diseases and enlarge the map of shared genetic risk loci.