Psoriasis is a common inflammatory skin disease that has been reported to be associated with obesity. We aimed to investigate a possible causal relationship between body mass index (BMI) and psoriasis.

Psoriasis vulgaris (PsV) risk is strongly associated with variation within the major histocompatibility complex (MHC) region, but its genetic architecture has yet to be fully elucidated. Here, we conducted a large-scale fine-mapping study of PsV risk in the MHC region in 9,247 PsV-affected individuals and 13,589 controls of European descent by imputing class I and II human leukocyte antigen (HLA) genes from SNP genotype data. In addition, we imputed sequence variants for MICA, an MHC HLA-like gene that has been associated with PsV, to evaluate association at that locus as well. We observed that HLA-C(∗)06:02 demonstrated the lowest p value for overall PsV risk (p = 1.7 × 10(-364)). Stepwise analysis revealed multiple HLA-C(∗)06:02-independent risk variants in both class I and class II HLA genes for PsV susceptibility (HLA-C(∗)12:03, HLA-B amino acid positions 67 and 9, HLA-A amino acid position 95, and HLA-DQα1 amino acid position 53; p < 5.0 × 10(-8)), but no apparent risk conferred by MICA. We further evaluated risk of two major clinical subtypes of PsV, psoriatic arthritis (PsA; n = 3,038) and cutaneous psoriasis (PsC; n = 3,098). We found that risk heterogeneity between PsA and PsC might be driven by HLA-B amino acid position 45 (Pomnibus = 2.2 × 10(-11)), indicating that different genetic factors underlie the overall risk of PsV and the risk of specific PsV subphenotypes. Our study illustrates the value of high-resolution HLA and MICA imputation for fine mapping causal variants in the MHC.

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.

We have shown that autocrine proliferation of human keratinocytes (KCs) is strongly dependent upon amphiregulin (AREG), whereas blockade of heparin-binding EGF-like growth factor (HB-EGF) inhibits KC migration in scratch wound assays. Here we demonstrate that expression of soluble HB-EGF (sHB-EGF) or full-length transmembrane HB-EGF (proHB-EGF), but not proAREG, results in profound increases in KC migration and invasiveness in monolayer culture. Coincident with these changes, HB-EGF significantly decreases mRNA expression of several epithelial markers including keratins 1, 5, 10, and 14 while increasing expression of markers of cellular motility including SNAI1, ZEB1, COX-2, and MMP1. Immunostaining revealed HB-EGF-induced expression of the mesenchymal protein vimentin and decreased expression of E-cadherin, as well as nuclear translocation of β-catenin. Suggestive of a trade-off between KC motility and proliferation, overexpression of HB-EGF also reduced KC growth by >90%. We also show that HB-EGF is strongly induced in regenerating epidermis after partial-thickness wounding of human skin. Taken together, our data suggest that expression of HB-EGF in human KCs triggers a migratory and invasive phenotype with many features of epithelial-mesenchymal transition (EMT), which may be beneficial in the context of cutaneous wound healing.

Merkel cell carcinoma (MCC) is an aggressive cutaneous neuroendocrine tumor with high mortality rates. Merkel cell polyomavirus (MCPyV), identified in the majority of MCCs, may drive tumorigenesis via viral T antigens. However, the mechanisms underlying pathogenesis in MCPyV-negative MCCs remain poorly understood. To nominate genes contributing to the pathogenesis of MCPyV-negative MCCs, we performed DNA microarray analysis on 30 MCCs. The MCPyV status of MCCs was determined by PCR for viral DNA and RNA. A total of 1,593 probe sets were differentially expressed between MCPyV-negative and MCPyV-positive MCCs, with significant differential expression defined as at least a 2-fold change in either direction and a P-value 0.05. MCPyV-negative tumors showed decreased RB1 expression, whereas MCPyV-positive tumors were enriched for immune response genes. Validation studies included immunohistochemistry demonstration of decreased RB protein expression in MCPyV-negative tumors and increased peritumoral CD8+ T lymphocytes surrounding MCPyV-positive tumors. In conclusion, our data suggest that loss of RB1 expression may have an important role in the tumorigenesis of MCPyV-negative MCCs. Functional and clinical validation studies are needed to determine whether this tumor-suppressor pathway represents an avenue for targeted therapy.

Psoriasis, an immune-mediated, inflammatory disease of the skin and joints, provides an ideal system for expression quantitative trait locus (eQTL) analysis, because it has a strong genetic basis and disease-relevant tissue (skin) is readily accessible. To better understand the role of genetic variants regulating cutaneous gene expression, we identified 841 cis-acting eQTLs using RNA extracted from skin biopsies of 53 psoriatic individuals and 57 healthy controls. We found substantial overlap between cis-eQTLs of normal control, uninvolved psoriatic, and lesional psoriatic skin. Consistent with recent studies and with the idea that control of gene expression can mediate relationships between genetic variants and disease risk, we found that eQTL SNPs are more likely to be associated with psoriasis than are randomly selected SNPs. To explore the tissue specificity of these eQTLs and hence to quantify the benefits of studying eQTLs in different tissues, we developed a refined statistical method for estimating eQTL overlap and used it to compare skin eQTLs to a published panel of lymphoblastoid cell line (LCL) eQTLs. Our method accounts for the fact that most eQTL studies are likely to miss some true eQTLs as a result of power limitations and shows that ∼70% of cis-eQTLs in LCLs are shared with skin, as compared with the naive estimate of < 50% sharing. Our results provide a useful method for estimating the overlap between various eQTL studies and provide a catalog of cis-eQTLs in skin that can facilitate efforts to understand the functional impact of identified susceptibility variants on psoriasis and other skin traits.

We carried out a meta-analysis of two recent psoriasis genome-wide association studies with a combined discovery sample of 1,831 affected individuals (cases) and 2,546 controls. One hundred and two loci selected based on P value rankings were followed up in a three-stage replication study including 4,064 cases and 4,685 controls from Michigan, Toronto, Newfoundland and Germany. In the combined meta-analysis, we identified three new susceptibility loci, including one at NOS2 (rs4795067, combined P = 4 × 10⁻¹¹), one at FBXL19 (rs10782001, combined P = 9 × 10⁻¹⁰) and one near PSMA6-NFKBIA (rs12586317, combined P = 2 × 10⁻⁸). All three loci were also associated with psoriatic arthritis (rs4795067, combined P = 1 × 10⁻⁵; rs10782001, combined P = 4 × 10⁻⁸; and rs12586317, combined P = 6 × 1⁻⁵) and purely cutaneous psoriasis (rs4795067, combined P = 1 × 10⁻⁸; rs10782001, combined P = 2 × 10⁻⁶; and rs12586317, combined P = 1 × 10⁻⁶). We also replicated a recently identified association signal near RNF114 (rs495337, combined P = 2 × 10⁻⁷).

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.

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 psoriasis cases and 1,436 controls of European ancestry. We followed up 21 promising SNPs 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 combined P < 5 x 10(-8)). Loci with confirmed association include HLA-C, three genes involved in IL-23 signaling (IL23A, IL23R, IL12B), two genes that act downstream of TNF-alpha and regulate NF-kappaB 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.

We recently conducted a large association analysis to compare the genetic profiles between patients with psoriatic arthritis (PsA) and cutaneous-only psoriasis (PsC). Despite including over 7,000 genotyped patients, only the MHC achieved genome-wide significance. In this study, we hypothesized that appropriate epigenomic elements (H3K27ac marks for active enhancers) can guide us to reveal valuable information about the loci with suggestive evidence of association. Our aim is to investigate these loci and explore how they may lead to the development of PsA. We evaluated this potential by investigating the genes connected with these loci from the perspective of pharmacogenomics and gene expression. We illustrated that markers with suggestive evidence of association outside the MHC region are enriched in H3K27ac marks for osteoblast and chondrogenic differentiated cells; using pharmacogenomics resources, we showed that genes near these markers are targeted by existing drugs used to treat psoriatic arthritis. Significantly, six of the ten suggestive significant loci overlapping the regulatory elements encompass genes differentially expressed (FDR < 5%) in differentiated osteoblasts, including genes participating in the Wnt signaling such as RUNX1, FUT8, and CTNNAL1. Our approach shows that epigenomic information can be used as cost-effective approach to enhance the inferences for GWAS results, especially in situations when few genome-wide significant loci are available. Our results also point the way to more directed investigations comparing the genetics of PsA and PsC.

To increase our understanding of psoriasis, we used high-throughput complementary DNA sequencing (RNA-seq) to assay the transcriptomes of lesional psoriatic and normal skin. We sequenced polyadenylated RNA-derived complementary DNAs 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 coexpression network analysis (WGCNA) revealed multiple modules of coordinately expressed epidermal differentiation genes, overlapping significantly with genes regulated by the long noncoding 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 downregulated in psoriatic biopsies, most likely because of expansion of the epidermal compartment. These results show the power of WGCNA to elucidate gene regulatory circuits in psoriasis, and emphasize the influence of tissue architecture in both differential expression and coexpression analysis.

Although analysis pipelines have been developed to use RNA-seq to identify long non-coding RNAs (lncRNAs), inference of their biological and pathological relevance remains a challenge. As a result, most transcriptome studies of autoimmune disease have only assessed protein-coding transcripts.

The ability to obtain gene expression profiles from human disease specimens provides an opportunity to identify relevant gene pathways, but is limited by the absence of data sets spanning a broad range of conditions. Here, we analyzed publicly available microarray data from 16 diverse skin conditions in order to gain insight into disease pathogenesis. Unsupervised hierarchical clustering separated samples by disease as well as common cellular and molecular pathways. Disease-specific signatures were leveraged to build a multi-disease classifier, which predicted the diagnosis of publicly and prospectively collected expression profiles with 93% accuracy. In one sample, the molecular classifier differed from the initial clinical diagnosis and correctly predicted the eventual diagnosis as the clinical presentation evolved. Finally, integration of IFN-regulated gene programs with the skin database revealed a significant inverse correlation between IFN-β and IFN-γ programs across all conditions. Our study provides an integrative approach to the study of gene signatures from multiple skin conditions, elucidating mechanisms of disease pathogenesis. In addition, these studies provide a framework for developing tools for personalized medicine toward the precise prediction, prevention, and treatment of disease on an individual level.

Psoriasis has been linked to cardiometabolic diseases, but epidemiological findings are inconsistent. We investigated the association between psoriasis and cardiometabolic outcomes in a German cross-sectional study (n=4,185) and a prospective cohort of German Health Insurance beneficiaries (n=1,811,098). A potential genetic overlap was explored using genome-wide data from >22,000 coronary artery disease and >4,000 psoriasis cases, and with a dense genotyping study of cardiometabolic risk loci on 927 psoriasis cases and 3,717 controls. After controlling for major confounders, in the cross-sectional analysis psoriasis was significantly associated with type 2 diabetes (T2D, adjusted odds ratio (OR)=2.36; 95% confidence interval CI=1.26-4.41) and myocardial infarction (MI, OR=2.26; 95% CI=1.03-4.96). In the longitudinal study, psoriasis slightly increased the risk for incident T2D (adjusted relative risk (RR)=1.11; 95% CI=1.08-1.14) and MI (RR=1.14; 95% CI=1.06-1.22), with highest risk increments in systemically treated psoriasis, which accounted for 11 and 17 excess cases of T2D and MI per 10,000 person-years. Except for weak signals from within the major histocompatibility complex, there was no evidence of genetic risk loci shared between psoriasis and cardiometabolic traits. Our findings suggest that psoriasis, in particular severe psoriasis, increases the risk for T2D and MI, and that the genetic architecture of psoriasis and cardiometabolic traits is largely distinct.

Psoriasis is a chronic autoimmune disease with complex genetic architecture. Previous genome-wide association studies (GWAS) and a recent meta-analysis using Immunochip data have uncovered 36 susceptibility loci. Here, we extend our previous meta-analysis of European ancestry by refined genotype calling and imputation and by the addition of 5,033 cases and 5,707 controls. The combined analysis, consisting of over 15,000 cases and 27,000 controls, identifies five new psoriasis susceptibility loci at genome-wide significance (P<5 × 10(-8)). The newly identified signals include two that reside in intergenic regions (1q31.1 and 5p13.1) and three residing near PLCL2 (3p24.3), NFKBIZ (3q12.3) and CAMK2G (10q22.2). We further demonstrate that NFKBIZ is a TRAF3IP2-dependent target of IL-17 signalling in human skin keratinocytes, thereby functionally linking two strong candidate genes. These results further integrate the genetics and immunology of psoriasis, suggesting new avenues for functional analysis and improved therapies.

Keratinocyte (KC) hyper-proliferation and epidermal thickening are characteristic features of psoriasis lesions, but the specific contributions of KCs to plaque formation are not fully understood. This study used RNA-seq to investigate the transcriptome of primary monolayer KC cultures grown from lesional (PP) and non-lesional (PN) biopsies of psoriasis patients and control subjects (NN). Whole skin biopsies from the same subjects were evaluated concurrently. RNA-seq analysis of whole skin identified a larger number of psoriasis-increased differentially expressed genes (DEGs), but analysis of KC cultures identified more PP- and PN-decreased DEGs. These latter DEG sets overlapped more strongly with genes near loci identified by psoriasis genome-wide association studies and were enriched for genes associated with epidermal differentiation. Consistent with this, the frequency of AP-1 motifs was elevated in regions upstream of PN-KC-decreased DEGs. A subset of these genes belonged to the same co-expression module, mapped to the epidermal differentiation complex, and exhibited differentiation-dependent expression. These findings demonstrate a decreased differentiation gene signature in PP/PN-KCs that had not been identified by pre-genomic studies of patient-derived monolayers. This may reflect intrinsic defects limiting psoriatic KC differentiation capacity, which may contribute to compromised barrier function in normal-appearing uninvolved psoriatic skin.

Skin aging is associated with intrinsic processes that compromise the structure of the extracellular matrix while promoting loss of functional and regenerative capacity. These processes are accompanied by a large-scale shift in gene expression, but underlying mechanisms are not understood and conservation of these mechanisms between humans and mice is uncertain.

EphA2 is a receptor tyrosine kinase (RTK) that triggers keratinocyte differentiation upon activation and subsequent downregulation by ephrin-A1 ligand. The objective of this study was to determine whether the EphA2/ephrin-A1 signaling axis was altered in psoriasis, an inflammatory skin condition in which keratinocyte differentiation is abnormal. Microarray analysis of skin biopsies from psoriasis patients revealed increased mRNA transcripts for several members of this RTK family in plaques, including the EphA1, EphA2, and EphA4 subtypes prominently expressed by keratinocytes. Of these, EphA2 showed the greatest upregulation, a finding that was confirmed by quantitative reverse-transcriptase-PCR, immunohistochemistry (IHC), and ELISA. In contrast, psoriatic lesions exhibited reduced ephrin-A ligand immunoreactivity. Exposure of primary keratinocytes induced to differentiate in high calcium or a three-dimensional (3D) raft culture of human epidermis to a combination of growth factors and cytokines elevated in psoriasis increased EphA2 mRNA and protein expression while inducing S100A7 and disrupting differentiation. Pharmacological delivery of a soluble ephrin-A1 peptidomimetic ligand led to a reduction in EphA2 expression and ameliorated proliferation and differentiation in raft cultures exposed to EGF and IL-1α. These findings suggest that ephrin-A1-mediated downregulation of EphA2 supports keratinocyte differentiation in the context of cytokine perturbation.

The tumor necrosis factor-alpha-induced protein 3 (TNFAIP3) gene has been associated with psoriasis, rheumatoid arthritis, type 1 diabetes mellitus, systemic lupus erythematosus, and celiac disease. TNFAIP3 encodes A20, a tumor necrosis factor (TNF)-α-inducible zinc finger protein thought to limit NF-κB-mediated immune responses. In this study, we report association of response of psoriasis to TNF blockers with two TNFAIP3 single-nucleotide polymorphisms (rs2230926 in exon 7 and rs610604 in intron 3) and their haplotypes. Treatment response was self-evaluated using a 0-5 visual analog scale in 433 psoriasis patients who received TNF blockers. Confirmation was sought in 199 psoriasis and psoriatic arthritis patients from Toronto who were followed up prospectively. Response variables were dichotomized separately in the two cohorts, yielding similar proportions of good responses. Whereas significant associations were observed only for the Michigan cohort, fixed-effects meta-analysis retained significant association between dosage of the G allele of rs610604 and good combined response to all TNF blockers (odds ratio (OR) = 1.50, P(corr) = 0.050) and etanercept (OR = 1.64, P(corr) = 0.016). The rs2230926 T-rs610604 G haplotype was similarly associated. By demonstrating an association with therapeutic response, these results provide a clinically relevant functional correlate to the recently described genetic association between psoriasis and TNFAIP3.

Autoimmune diseases affect 7.5% of the US population, and they are among the leading causes of death and disability. A notable feature of many autoimmune diseases is their greater prevalence in females than in males, but the underlying mechanisms of this have remained unclear. Through the use of high-resolution global transcriptome analyses, we demonstrated a female-biased molecular signature associated with susceptibility to autoimmune disease and linked this to extensive sex-dependent co-expression networks. This signature was independent of biological age and sex-hormone regulation and was regulated by the transcription factor VGLL3, which also had a strong female-biased expression. On a genome-wide level, VGLL3-regulated genes had a strong association with multiple autoimmune diseases, including lupus, scleroderma and Sjögren's syndrome, and had a prominent transcriptomic overlap with inflammatory processes in cutaneous lupus. These results identified a VGLL3-regulated network as a previously unknown inflammatory pathway that promotes female-biased autoimmunity. They demonstrate the importance of studying immunological processes in females and males separately and suggest new avenues for therapeutic development.