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.

Psoriasis is a common inflammatory skin disease with a prevalence of 2% to 3% in Caucasians1. In a genome-wide search for copy number variants (CNV) using a sample pooling approach we have identified a deletion comprising LCE3B and LCE3C, members of the late cornified envelope (LCE) gene cluster2. The absence of LCE3B and LCE3C (LCE3C-LCE3B-del) is significantly associated (p=1.38E-08) with risk of psoriasis in 2,831 samples from Spain, The Netherlands, Italy and the USA, and in a family-based study (p=5.4E-04). LCE3C-LCE3B-del is tagged by rs4112788 (r2=0.93), which is also strongly associated with psoriasis (p<6.6E-09). LCE3C-LCE3B-del shows epistatic effects with the HLA-Cw6 allele on the development of psoriasis in Dutch samples, and multiplicative effects in the other samples. LCE expression can be induced in normal epidermis by skin barrier disruption and is strongly expressed in psoriatic lesions, suggesting that compromised skin barrier function plays a role in psoriasis susceptibility.

Previous studies have narrowed the interval containing PSORS1, the psoriasis-susceptibility locus in the major histocompatibility complex (MHC), to an approximately 300-kb region containing HLA-C and at least 10 other genes. In an effort to identify the PSORS1 gene, we cloned and completely sequenced this region from both chromosomes of five individuals. Two of the sequenced haplotypes were associated with psoriasis (risk), and the other eight were clearly unassociated (nonrisk). Comparison of sequence of the two risk haplotypes identified a 298-kb region of homology, extending from just telomeric of HLA-B to the HCG22 gene, which was flanked by clearly nonhomologous regions. Similar haplotypes cloned from unrelated individuals had nearly identical sequence. Combinatorial analysis of exonic variations in the known genes of the candidate interval revealed that HCG27, PSORS1C3, OTF3, TCF19, HCR, STG, and HCG22 bore no alleles unique to risk haplotypes among the 10 sequenced haplotypes. SPR1 and SEEK1 both had messenger RNA alleles specific to risk haplotypes, but only HLA-C and CDSN yielded protein alleles unique to risk. The risk alleles of HLA-C and CDSN (HLA-Cw6 and CDSN*TTC) were genotyped in 678 families with early-onset psoriasis; 620 of these families were also typed for 34 microsatellite markers spanning the PSORS1 interval. Recombinant haplotypes retaining HLA-Cw6 but lacking CDSN*TTC were significantly associated with psoriasis, whereas recombinants retaining CDSN*TTC but lacking HLA-Cw6 were not associated, despite good statistical power. By grouping recombinants with similar breakpoints, the most telomeric quarter of the 298-kb candidate interval could be excluded with high confidence. These results strongly suggest that HLA-Cw6 is the PSORS1 risk allele that confers susceptibility to early-onset psoriasis.

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 common inflammatory disorder of the skin and other organs. We have determined that mutations in CARD14, encoding a nuclear factor of kappa light chain enhancer in B cells (NF-kB) activator within skin epidermis, account for PSORS2. Here, we describe fifteen additional rare missense variants in CARD14, their distribution in seven psoriasis cohorts (>6,000 cases and >4,000 controls), and their effects on NF-kB activation and the transcriptome of keratinocytes. There were more CARD14 rare variants in cases than in controls (burden test p value = 0.0015). Some variants were only seen in a single case, and these included putative pathogenic mutations (c.424G>A [p.Glu142Lys] and c.425A>G [p.Glu142Gly]) and the generalized-pustular-psoriasis mutation, c.413A>C (p.Glu138Ala); these three mutations lie within the coiled-coil domain of CARD14. The c.349G>A (p.Gly117Ser) familial-psoriasis mutation was present at a frequency of 0.0005 in cases of European ancestry. CARD14 variants led to a range of NF-kB activities; in particular, putative pathogenic variants led to levels >2.5× higher than did wild-type CARD14. Two variants (c.511C>A [p.His171Asn] and c.536G>A [p.Arg179His]) required stimulation with tumor necrosis factor alpha (TNF-α) to achieve significant increases in NF-kB levels. Transcriptome profiling of wild-type and variant CARD14 transfectants in keratinocytes differentiated probably pathogenic mutations from neutral variants such as polymorphisms. Over 20 CARD14 polymorphisms were also genotyped, and meta-analysis revealed an association between psoriasis and rs11652075 (c.2458C>T [p.Arg820Trp]; p value = 2.1 × 10(-6)). In the two largest psoriasis cohorts, evidence for association increased when rs11652075 was conditioned on HLA-Cw*0602 (PSORS1). These studies contribute to our understanding of the genetic basis of psoriasis and illustrate the challenges faced in identifying pathogenic variants in common disease.

In a 12.5 cM genome-wide scan for psoriasis susceptibility loci, recombination-based tests revealed linkage to the HLA region (Zmax = 3.52), as well as suggestive linkage to two novel regions: chromosome 16q (60-83.1 cM from pter, Zmax = 2.50), and chromosome 20p (7.5-25 cM from pter, Zmax = 2.62). All three regions yielded P values < or = 0.01 by non-parametric analysis. Recombination-based and allele sharing methods also confirmed a previous report of a dominant susceptibility locus on distal chromosome 17q (108.2 cM from pter, Zmax = 2.09, GENEHUNTER P = 0.0056). We could not confirm a previously reported locus on distal chromosome 4q; however, a broad region of unclear significance was identified proximal to this proposed locus (153.6-178.4 cM from pter, Zmax = 1.01). Taken together with our recent results demonstrating linkage to HLA-B and -C, this genome-wide scan identifies a psoriasis susceptibility locus at HLA, confirms linkage to 17q, and recommends two novel genomic regions for further scrutiny. One of these regions (16q) overlaps with a recently-identified susceptibility locus for Crohn's disease. Psoriasis is much more common in patients with Crohn's disease than in controls, suggesting that an immunomodulatory locus capable of influencing both diseases may reside in this region.

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.