Retrospective pharmacogenetic study of psoriasis highlights the role of KLK7 in tumour necrosis factor signalling

BACKGROUND

Multiple treatment options are available for the management of psoriasis, but clinical response varies among individual patients and no biomarkers are available to facilitate treatment selection for improved patient outcomes.

OBJECTIVES

To utilize retrospective data to conduct a pharmacogenetic study to explore the potential genetic pathways associated with drug response in the treatment of psoriasis.

METHODS

We conducted a retrospective pharmacogenetic study using self-evaluated treatment response from 1942 genotyped patients with psoriasis. We examined 6 502 658 genetic markers to model their associations with response to six treatment options using linear regression, adjusting for cohort variables and demographic features. We further utilized an integrative approach incorporating epigenomics, transcriptomics and a longitudinal clinical cohort to provide biological implications for the topmost signals associated with drug response.

RESULTS

Two novel markers were revealed to be associated with treatment response: rs1991820 (P = 1.30 × 10-6) for anti-tumour necrosis factor (TNF) biologics; and rs62264137 (P = 2.94 × 10-6) for methotrexate, which was also associated with cutaneous mRNA expression levels of two known psoriasis-related genes KLK7 (P = 1.0 × 10-12) and CD200 (P = 5.4 × 10-6). We demonstrated that KLK7 expression was increased in the psoriatic epidermis, as shown by immunohistochemistry, as well as single-cell RNA sequencing, and its responsiveness to anti-TNF treatment was highlighted. By inhibiting the expression of KLK7, we further illustrated that keratinocytes have decreased proinflammatory responses to TNF.

CONCLUSIONS

Our study implicates the genetic regulation of cytokine responses in predicting clinical drug response and supports the association between pharmacogenetic loci and anti-TNF response, as shown here for KLK7.

GWAS meta-analysis of psoriasis identifies new susceptibility alleles impacting disease mechanisms and therapeutic targets

Psoriasis is a common, debilitating immune-mediated skin disease. Genetic studies have identified biological mechanisms of psoriasis risk, including those targeted by effective therapies. However, the genetic liability to psoriasis is not fully explained by variation at robustly identified risk loci. To move towards a saturation map of psoriasis susceptibility we meta-analysed 18 GWAS comprising 36,466 cases and 458,078 controls and identified 109 distinct psoriasis susceptibility loci, including 45 that have not been previously reported. These include susceptibility variants at loci in which the therapeutic targets IL17RA and AHR are encoded, and deleterious coding variants supporting potential new drug targets (including in STAP2, CPVL and POU2F3). We conducted a transcriptome-wide association study to identify regulatory effects of psoriasis susceptibility variants and cross-referenced these against single cell expression profiles in psoriasis-affected skin, highlighting roles for the transcriptional regulation of haematopoietic cell development and epigenetic modulation of interferon signalling in psoriasis pathobiology.

Shared Genetic Risk Factors for Multiple Sclerosis/Psoriasis Suggest Involvement of Interleukin-17 and Janus Kinase-Signal Transducers and Activators of Transcription Signaling

OBJECTIVE

Psoriasis and multiple sclerosis (MS) are complex immune diseases that are mediated by T cells and share multiple comorbidities. Previous studies have suggested psoriatic patients are at higher risk of MS; however, causal relationships between the two conditions remain unclear. Through epidemiology and genetics, we provide a comprehensive understanding of the relationship, and share molecular factors between psoriasis and MS.

METHODS

We used logistic regression, trans-disease meta-analysis and Mendelian randomization. Medical claims data were included from 30 million patients, including 141,544 with MS and 742,919 with psoriasis. We used genome-wide association study summary statistics from 11,024 psoriatic, 14,802 MS cases, and 43,039 controls for trans-disease meta-analysis, with additional summary statistics from 5 million individuals for Mendelian randomization.

RESULTS

Psoriatic patients have a significantly higher risk of MS (4,637 patients with both diseases; odds ratio [OR] 1.07, p = 1.2 × 10-5 ) after controlling for potential confounders. Using inverse variance and equally weighted trans-disease meta-analysis, we revealed >20 shared and opposing (direction of effect) genetic loci outside the major histocompatibility complex that showed significant genetic colocalization (in COLOC and COLOC-SuSiE v5.1.0). Co-expression analysis of genes from these loci further identified distinct clusters that were enriched among pathways for interleukin-17/tumor necrosis factor-α (OR >39, p < 1.6 × 10-3 ) and Janus kinase-signal transducers and activators of transcription (OR 35, p = 1.1 × 10-5 ), including genes, such as TNFAIP3, TYK2, and TNFRSF1A. Mendelian randomization found psoriasis as an exposure has a significant causal effect on MS (OR 1.04, p = 5.8 × 10-3 ), independent of type 1 diabetes (OR 1.05, p = 4.3 × 10-7 ), type 2 diabetes (OR 1.08, p = 2.3 × 10-3 ), inflammatory bowel disease (OR 1.11, p = 1.6 × 10-11 ), and vitamin D level (OR 0.75, p = 9.4 × 10-3 ).

INTERPRETATION

By investigating the shared genetics of psoriasis and MS, along with their modifiable risk factors, our findings will advance innovations in treatment for patients suffering from comorbidities. ANN NEUROL 2023;94:384-397.

Transethnic analysis of psoriasis susceptibility in South Asians and Europeans enhances fine-mapping in the MHC and genomewide

Because transethnic analysis may facilitate prioritization of causal genetic variants, we performed a genomewide association study (GWAS) of psoriasis in South Asians (SAS), consisting of 2,590 cases and 1,720 controls. Comparison with our existing European-origin (EUR) GWAS showed that effect sizes of known psoriasis signals were highly correlated in SAS and EUR (Spearman ρ = 0.78; p < 2 × 10-14). Transethnic meta-analysis identified two non-MHC psoriasis loci (1p36.22 and 1q24.2) not previously identified in EUR, which may have regulatory roles. For these two loci, the transethnic GWAS provided higher genetic resolution and reduced the number of potential causal variants compared to using the EUR sample alone. We then explored multiple strategies to develop reference panels for accurately imputing MHC genotypes in both SAS and EUR populations and conducted a fine-mapping of MHC psoriasis associations in SAS and the largest such effort for EUR. HLA-C*06 was the top-ranking MHC locus in both populations but was even more prominent in SAS based on odds ratio, disease liability, model fit and predictive power. Transethnic modeling also substantially boosted the probability that the HLA-C*06 protein variant is causal. Secondary MHC signals included coding variants of HLA-C and HLA-B, but also potential regulatory variants of these two genes as well as HLA-A and several HLA class II genes, with effects on both chromatin accessibility and gene expression. This study highlights the shared genetic basis of psoriasis in SAS and EUR populations and the value of transethnic meta-analysis for discovery and fine-mapping of susceptibility loci.

Epicardial adipose tissue volume is greater in men with severe psoriasis, implying an increased cardiovascular disease risk: A cross-sectional study

BACKGROUND

Patients with psoriasis have elevated risk of coronary artery disease.

OBJECTIVE

Do patients with severe psoriasis have larger epicardial adipose tissue volumes (EAT-V) that are associated with cardiovascular risk?

METHODS

For this cross-sectional study, we recruited dermatology patients with severe psoriasis and control patients without psoriasis or rheumatologic disease themselves or in a first-degree relative. Participants aged 34 to 55 years without known coronary artery disease or diabetes mellitus underwent computed tomography (CT); EAT-V was obtained from noncontrast CT heart images.

RESULTS

Twenty-five patients with psoriasis (14 men, 11 women) and 16 controls (5 men, 11 women) participated. Groups had no statistical difference in age, body mass index, various cardiovascular risk factors (except high-sensitivity C-reactive protein in men), CT-determined coronary artery calcium scores or plaque, or family history of premature cardiovascular disease. Mean EAT-V was greater in the psoriasis group compared to controls (P = .04). There was no statistically significant difference among women; however, male patients with psoriasis had significantly higher EAT-V than controls (P = .03), even when corrected for elevated high-sensitivity C-reactive protein (P = .05).

LIMITATIONS

A single-center convenience sample may not be representative.

CONCLUSION

Males with psoriasis without known coronary disease or diabetes had greater EAT-V than controls. EAT-V may be an early identifier of those at increased risk for cardiovascular events.

Large-Scale Imputation of KIR Copy Number and HLA Alleles in North American and European Psoriasis Case-Control Cohorts Reveals Association of Inhibitory KIR2DL2 With 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.

Causal Relationship and Shared Genetic Loci between Psoriasis and Type 2 Diabetes through Trans-Disease Meta-Analysis

Psoriasis and type 2 diabetes (T2D) are complex conditions with significant impacts on health. Patients with psoriasis have a higher risk of T2D (∼1.5 OR) and vice versa, controlling for body mass index; yet, there has been a limited study comparing their genetic architecture. We hypothesized that there are shared genetic components between psoriasis and T2D. Trans-disease meta-analysis was applied to 8,016,731 well-imputed genetic markers from large-scale meta-analyses of psoriasis (11,024 cases and 16,336 controls) and T2D (74,124 cases and 824,006 controls), adjusted for body mass index. We confirmed our findings in a hospital-based study (42,112 patients) and tested for causal relationships with multivariable Mendelian randomization. Mendelian randomization identified a causal relationship between psoriasis and T2D (P = 1.6 × 10^‒4, OR = 1.01) and highlighted the impact of body mass index. Trans-disease meta-analysis further revealed four genome-wide significant loci (P < 5 × 10^‒8) with evidence of colocalization and shared directions of effect between psoriasis and T2D not present in body mass index. The proteins coded by genes in these loci (ACTR2, ERLIN1, TRMT112, and BECN1) are connected through NF-κB signaling. Our results provide insight into the immunological components that connect immune-mediated skin conditions and metabolic diseases, independent of confounding factors.

Systemic evaluation of the relationship between psoriasis, psoriatic arthritis and osteoporosis: observational and Mendelian randomisation study

OBJECTIVES AND METHODS

With 432 513 samples from UK Biobank dataset, multivariable linear/logistic regression were used to estimate the relationship between psoriasis/psoriatic arthritis (PsA) and estimated bone mineral density (eBMD)/osteoporosis, controlling for potential confounders. Here, confounders were set in three ways: model0 (including age, height, weight, smoking and drinking), model1 (model0 +regular physical activity) and model2 (model1 +medication treatments). The eBMD was derived from heel ultrasound measurement. And 4904 patients with psoriasis and 847 patients with PsA were included in final analysis. Mendelian randomisation (MR) approach was used to evaluate the causal effect between them.

RESULTS

Lower eBMD were observed in patients with PsA than in controls in both model0 (β-coefficient=-0.014, p=0.0006) and model1 (β-coefficient=-0.013, p=0.002); however, the association disappeared when conditioning on treatment with methotrexate or ciclosporin (model2) (β-coefficient=-0.005, p=0.28), mediation analysis showed that 63% of the intermediary effect on eBMD was mediated by medication treatment (p<2E-16). Patients with psoriasis without arthritis showed no difference of eBMD compared with controls. Similarly, the significance of higher risk of osteopenia in patients with PsA (OR=1.27, p=0.002 in model0) could be eliminated by conditioning on medication treatment (p=0.244 in model2). Psoriasis without arthritis was not related to osteopenia and osteoporosis. The weighted Genetic Risk Score analysis found that genetically determined psoriasis/PsA were not associated with eBMD (p=0.24 and p=0.88). Finally, MR analysis showed that psoriasis/PsA had no causal effect on eBMD, osteoporosis and fracture.

CONCLUSIONS

The effect of PsA on osteoporosis was secondary (eg, medication) but not causal. Under this hypothesis, psoriasis without arthritis was not a risk factor for osteoporosis.

Quantifying Differences in Heritability among Psoriatic Arthritis (PsA), Cutaneous Psoriasis (PsC) and Psoriasis vulgaris (PsV)

Chronic plaque psoriasis and psoriatic arthritis are multifactorial inter-related diseases with strong genetic contributions. Better elucidation of the heritability of psoriatic disease subsets is important for identifying novel genes, risk stratification and potential clinical applications. In this study, we used two mixed-effect modelling methodologies to assess the additive contribution of common single nucleotide polymorphisms from genome-wide association studies to estimate the heritability of cutaneous psoriasis, psoriasis vulgaris and psoriatic arthritis. We found that cutaneous psoriasis and psoriatic arthritis both exhibit considerable heritability, with a greater contribution coming from cutaneous psoriasis.

Integrative Approach to Reveal Cell Type Specificity and Gene Candidates for Psoriatic Arthritis Outside the MHC

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.

Neutrophil Extracellular Traps Induce Human Th17 Cells: Effect of Psoriasis-Associated TRAF3IP2 Genotype

Psoriasis lesions are rich in IL-17-producing T cells as well as neutrophils, which release webs of DNA-protein complexes known as neutrophil extracellular traps (NETs). Because we and others have observed increased NETosis in psoriatic lesions, we hypothesized that NETs contribute to increased T helper type 17 (Th17) cells in psoriasis. After stimulating peripheral blood mononuclear cells with anti-CD3/CD28 beads for 7 days, we found significantly higher percentages of CD3⁺CD4⁺IL-17⁺ (Th17) cells in the presence versus absence of NETs, as assessed by flow cytometry, IL-17 ELISA, and IL17A/F and RORC mRNAs. Memory, but not naïve, T cells were competent and monocytes were required for CD3/CD28-mediated Th17 induction, with or without NETs. Th17 induction was enhanced by the T allele of rs33980500 (T/C), a psoriasis risk-associated variant in the TRAF3IP2 gene encoding the D10N variant of Act1, a key mediator of IL-17 signal transduction. Global transcriptome analysis of CD3/CD28-stimulated peripheral blood mononuclear cells by RNA sequencing confirmed the stimulatory effects of NETs, demonstrated NET-induced enhancement of cytokine gene expression, and verified that the effect of Act1 D10N was greater in the presence of NETs. Collectively, these results implicate NETs and the Act1 D10N variant in human Th17 induction from peripheral blood mononuclear cells, with ramifications for immunogenetic studies of psoriasis and other autoimmune diseases.

The psoriasis-protective TYK2 I684S variant impairs IL-12 stimulated pSTAT4 response in skin-homing CD4+ and CD8+ memory T-cells

Tyrosine kinase 2 (TYK2) belongs to the Janus kinase (JAK) family of tyrosine kinases, which transmit signals from activated cytokine receptors. GWAS have consistently implicated TYK2 in psoriasis susceptibility. We performed an in-depth association analysis of TYK2 using GWAS and resequencing data. Strong genetic association of three nonsynonymous variants in the exonic regions of the TYK2 gene (rs34536443, rs12720356, and rs2304256) were found. rs12720356 encoding I684S is predicted to be deleterious based on its location in the pseudokinase domain. We analyzed PBMCs from 29 individuals representing the haplotypes containing each of the significantly associated signals. STAT4 phosphorylation was evaluated by phospho-flow cytometry after CD3/CD28 activation of cells followed by IL-12 stimulation. Individuals carrying the protective I684S variant manifested significantly reduced p-STAT4 levels in CD4 + CD25 + CD45RO+ (mean Stimulation Index (S.I.) 48.08, n = 10) and CD8 + CD25 + CD45RO + cells (S.I. 55.71, n = 10), compared to controls homozygous for the ancestral haplotype (S.I. 68.19, n = 10 (p = 0.002) and 76.76 n = 10 (p = 0.0008) respectively). Reduced p-STAT4 levels were also observed in skin-homing, cutaneous lymphocyte associated antigen (CLA)-positive CD4 and CD8 cells from I684S carriers. No significant changes in p-STAT4 for the psoriasis-associated variant rs34536443 was found. These data establish the functional significance of the TYK2 I684S variant in psoriasis susceptibility.

The Act1 D10N missense variant impairs CD40 signaling in human B-cells

The TRAF3IP2 gene resides within one of at least 63 psoriasis susceptibility loci and encodes Act1, an adapter protein involved in IL-17 receptor and CD40 signaling pathways. TRAF3IP2 is distinctive (among <10% of candidate susceptibility genes) in that a strongly disease-associated variant encodes a missense SNP predicted to be functionally relevant (SNP rs33980500 C/T encoding Act1 pD10N). As assessed by flow cytometry, Act1 protein was expressed at the highest levels in monocytes, with lower levels in T-cells and B-cells. However, monocytes, T-cells and B-cells failed to respond to IL-17A stimulation of PBMC, as measured by flow cytometric determination of NF-κB phospho-p65. As an alternative stimulus, we treated PBMCs with trimerized recombinant human CD40L and assessed p65, p38 and Erk phosphorylation in CD19⁺ B-cells as a function of D10N genotype. The increase of phosphorylated p65, p38 and Erk was well-correlated across individuals, and CD40L-induced phosphorylation of p65, p38, and Erk was significantly attenuated in B-cells from Act1 D10N homozygotes, compared to heterozygotes and nullizygotes. Our results indicate that the Act1 D10N variant is a relevant genetic determinant of CD40L responsiveness in human B-cells, with the risk allele being associated with lower B-cell responses in an acute signaling context.

Exome-wide association study reveals novel psoriasis susceptibility locus at TNFSF15 and rare protective alleles in genes contributing to type I IFN signalling

Psoriasis is a common inflammatory skin disorder for which multiple genetic susceptibility loci have been identified, but few resolved to specific functional variants. In this study, we sought to identify common and rare psoriasis-associated gene-centric variation. Using exome arrays we genotyped four independent cohorts, totalling 11 861 psoriasis cases and 28 610 controls, aggregating the dataset through statistical meta-analysis. Single variant analysis detected a previously unreported risk locus at TNFSF15 (rs6478108; P = 1.50 × 10-8, OR = 1.10), and association of common protein-altering variants at 11 loci previously implicated in psoriasis susceptibility. We validate previous reports of protective low-frequency protein-altering variants within IFIH1 (encoding an innate antiviral receptor) and TYK2 (encoding a Janus kinase), in each case establishing a further series of protective rare variants (minor allele frequency < 0.01) via gene-wide aggregation testing (IFIH1: pburden = 2.53 × 10-7, OR = 0.707; TYK2: pburden = 6.17 × 10-4, OR = 0.744). Both genes play significant roles in type I interferon (IFN) production and signalling. Several of the protective rare and low-frequency variants in IFIH1 and TYK2 disrupt conserved protein domains, highlighting potential mechanisms through which their effect may be exerted.

IFN-γ and TNF-α synergism may provide a link between psoriasis and inflammatory atherogenesis

Chronic inflammation is a critical component of atherogenesis, however, reliable human translational models aimed at characterizing these mechanisms are lacking. Psoriasis, a chronic inflammatory skin disease associated with increased susceptibility to atherosclerosis, provides a clinical human model that can be utilized to investigate the links between chronic inflammation and atherosclerosis development. We sought to investigate key biological processes in psoriasis skin and human vascular tissue to identify biological components that may promote atherosclerosis in chronic inflammatory conditions. Using a bioinformatics approach of human skin and vascular tissue, we determined IFN-γ and TNF-α are the dominant pro-inflammatory signals linking atherosclerosis and psoriasis. We then stimulated primary aortic endothelial cells and ex-vivo atherosclerotic tissue with IFN-γ and TNF-α and found they synergistically increased monocyte and T-cell chemoattractants, expression of adhesion molecules on the endothelial cell surface, and decreased endothelial barrier integrity in vitro, therefore increasing permeability. Our data provide strong evidence of synergism between IFN-γ and TNF- α in inflammatory atherogenesis and provide rationale for dual cytokine antagonism in future studies.

Serum bone-turnover biomarkers are associated with the occurrence of peripheral and axial arthritis in psoriatic disease: a prospective cross-sectional comparative study

Background

A recent systematic review identified four candidate serum-soluble bone-turnover biomarkers (dickkopf-1, Dkk-1; macrophage-colony stimulating factor, M-CSF; matrix metalloproteinase-3, MMP-3; osteoprotegerin, OPG) showing possible association with psoriatic arthritis (PsA). We aimed to: (i) confirm and determine if these four biomarkers are associated with PsA; (ii) differentiate psoriasis cases with and without arthritis; and (iii) differentiate PsA cases with and without axial arthritis.

Methods

A prospective cross-sectional comparative two-centre study recruited 200 patients with psoriasis without arthritis (PsC), 127 with PsA without axial arthritis (pPsA), 117 with PsA with axial arthritis (psoriatic spondyloarthritis, PsSpA), 157 with ankylosing spondylitis (AS) without psoriasis, and 50 matched healthy controls (HC). Serum biomarker concentrations were measured using ELISA. Multivariable regression and receiver operating characteristic analyses were performed.

Results

MMP-3 concentrations were significantly higher and M-CSF significantly lower in each arthritis disease group compared with HC (p ≤ 0.02). MMP-3 concentrations were significantly higher (adjusted odds ratio, OR_adj 1.02 per ng/ml increase in concentration; p = 0.0004) and M-CSF significantly lower (OR_adj 0.44 per ng/ml increase; p = 0.01) in PsA (pPsA and PsSpA combined) compared with PsC. Dkk-1 concentrations were significantly higher (OR_adj 1.22 per ng/mL increase; p = 0.01), and OPG concentrations significantly lower (OR_adj 0.20 per ng/mL increase; p = 0.02) in patients with axial arthritis (PsSpA and AS combined) than in those without (pPsA). Furthermore, Dkk-1 concentrations were significantly higher along a spectrum of increasing axial arthritis; Dkk-1 concentrations were higher in AS compared with PsSpA (OR_adj 1.18 per ng/mL increase; p = 0.02). Receiver operating characteristic analysis showed MMP-3 to be the best single biomarker for differentiating PsA from PsC (AUC 0.70 for a cut-off of 14.51 ng/mL; sensitivity 0.76, specificity 0.60).

Conclusions

MMP-3 and M-CSF are biomarkers for the presence of arthritis in psoriatic disease, and could therefore be used to screen for PsA in psoriasis cohorts. Dkk-1 and OPG are biomarkers of axial arthritis; they could therefore be used to screen for the presence of axial disease in PsA cases, and help differentiate PsSpA from AS. High concentrations of Dkk-1 in AS and PsSpA compared with HC, support previous reports that Dkk-1 is dysfunctional in the spondyloarthritides.

Psoriasis-Associated Late Cornified Envelope (LCE) Proteins Have Antibacterial Activity

Terminally differentiating epidermal keratinocytes express a large number of structural and antimicrobial proteins that are involved in the physical barrier function of the stratum corneum and provide innate cutaneous host defense. Late cornified envelope (LCE) genes, located in the epidermal differentiation complex on chromosome 1, encode a family of 18 proteins of unknown function, whose expression is largely restricted to epidermis. Deletion of two members, LCE3B and LCE3C (LCE3B/C-del), is a widely-replicated psoriasis risk factor that interacts with the major psoriasis-psoriasis risk gene HLA-C*06. Here we performed quantitative trait locus analysis, utilizing RNA-seq data from human skin and found that LCE3B/C-del was associated with a markedly increased expression of LCE3A, a gene directly adjacent to LCE3B/C-del. We confirmed these findings in a 3-dimensional skin model using primary keratinocytes from LCE3B/C-del genotyped donors. Functional analysis revealed that LCE3 proteins, and LCE3A in particular, have defensin-like antimicrobial activity against a variety of bacterial taxa at low micromolar concentrations. No genotype-dependent effect was observed for the inside-out or outside-in physical skin barrier function. Our findings identify an unknown biological function for LCE3 proteins and suggest a role in epidermal host defense and LCE3B/C-del-mediated psoriasis risk.

Large scale meta-analysis characterizes genetic architecture for common psoriasis associated variants

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 T_H0, T_H1 and T_H17). 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.

Endogenous Glucocorticoid Deficiency in Psoriasis Promotes Inflammation and Abnormal Differentiation

The factors involved in maintaining a localized inflammatory state in psoriatic skin remain poorly understood. Here, we demonstrate through metabolomic and transcriptomic profiling marked suppression of glucocorticoid biosynthesis in the epidermis of psoriatic skin leading to localized deficiency of cortisol. Utilizing a 3D human epidermis model, we demonstrate that glucocorticoid biosynthesis is suppressed by proinflammatory cytokines and that glucocorticoid deficiency promotes inflammatory responses in keratinocytes. Finally, we show in vitro and in vivo that treatment with topical glucocorticoids leads to rapid restoration of glucocorticoid biosynthesis gene expression coincident with normalization of epidermal differentiation and suppression of inflammatory responses. Taken together, our data suggest that localized glucocorticoid deficiency in psoriatic skin interferes with epidermal differentiation and promotes a sustained and localized inflammatory response. This may shed new light on the mechanism of action of topical steroids, and demonstrates the critical role of endogenous steroid in maintaining both inflammatory and differentiation homeostasis in the epidermis.

A gene network regulated by the transcription factor VGLL3 as a promoter of sex-biased autoimmune diseases

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.

Sebaceous Gland Atrophy in Psoriasis: An Explanation for Psoriatic Alopecia?

In a transcriptome study of lesional psoriatic skin (PP) versus normal skin, we found a coexpressed gene module (N5) enriched 11.5-fold for lipid biosynthetic genes. We also observed fewer visible hairs in PP skin, compared with uninvolved nonlesional psoriatic skin or normal skin (P < 0.0001). To ask whether these findings might be due to abnormalities of the pilosebaceous unit, we carried out three-dimensional morphometric analysis of paired PP and nonlesional psoriatic skin biopsies. Sebaceous glands were markedly atrophic in PP versus nonlesional psoriatic skin (91% average reduction in volume, P = 0.031). Module N5 genes were strongly downregulated in PP versus normal skin (fold change < 0.25, 44.4-fold) and strongly upregulated in sebaceous hyperplasia (fold change > 4, 54.1-fold). The intersection of PP-downregulated and sebaceous hyperplasia-upregulated gene lists generated a gene expression signature consisting solely of module N5 genes, whose expression in PP versus normal skin was inversely correlated with the signature of IL17-stimulated keratinocytes. Despite loss of visible hairs, morphometry identified elongated follicles in PP versus nonlesional psoriatic skin (average 1.7 vs. 1.2 μm, P = 0.020). These results document sebaceous gland atrophy in nonscalp psoriasis, identify a cytokine-regulated set of sebaceous gland signature genes, and suggest that loss of visible hair in PP skin may result from abnormal sebaceous gland function.

Investigating the Causal Relationship of C-Reactive Protein with 32 Complex Somatic and Psychiatric Outcomes: A Large-Scale Cross-Consortium Mendelian Randomization Study

BACKGROUND

C-reactive protein (CRP) is associated with immune, cardiometabolic, and psychiatric traits and diseases. Yet it is inconclusive whether these associations are causal.

METHODS AND FINDINGS

We performed Mendelian randomization (MR) analyses using two genetic risk scores (GRSs) as instrumental variables (IVs). The first GRS consisted of four single nucleotide polymorphisms (SNPs) in the CRP gene (GRSCRP), and the second consisted of 18 SNPs that were significantly associated with CRP levels in the largest genome-wide association study (GWAS) to date (GRSGWAS). To optimize power, we used summary statistics from GWAS consortia and tested the association of these two GRSs with 32 complex somatic and psychiatric outcomes, with up to 123,865 participants per outcome from populations of European ancestry. We performed heterogeneity tests to disentangle the pleiotropic effect of IVs. A Bonferroni-corrected significance level of less than 0.0016 was considered statistically significant. An observed p-value equal to or less than 0.05 was considered nominally significant evidence for a potential causal association, yet to be confirmed. The strengths (F-statistics) of the IVs were 31.92-3,761.29 and 82.32-9,403.21 for GRSCRP and GRSGWAS, respectively. CRP GRSGWAS showed a statistically significant protective relationship of a 10% genetically elevated CRP level with the risk of schizophrenia (odds ratio [OR] 0.86 [95% CI 0.79-0.94]; p < 0.001). We validated this finding with individual-level genotype data from the schizophrenia GWAS (OR 0.96 [95% CI 0.94-0.98]; p < 1.72 × 10-6). Further, we found that a standardized CRP polygenic risk score (CRPPRS) at p-value thresholds of 1 × 10-4, 0.001, 0.01, 0.05, and 0.1 using individual-level data also showed a protective effect (OR < 1.00) against schizophrenia; the first CRPPRS (built of SNPs with p < 1 × 10-4) showed a statistically significant (p < 2.45 × 10-4) protective effect with an OR of 0.97 (95% CI 0.95-0.99). The CRP GRSGWAS showed that a 10% increase in genetically determined CRP level was significantly associated with coronary artery disease (OR 0.88 [95% CI 0.84-0.94]; p < 2.4 × 10-5) and was nominally associated with the risk of inflammatory bowel disease (OR 0.85 [95% CI 0.74-0.98]; p < 0.03), Crohn disease (OR 0.81 [95% CI 0.70-0.94]; p < 0.005), psoriatic arthritis (OR 1.36 [95% CI 1.00-1.84]; p < 0.049), knee osteoarthritis (OR 1.17 [95% CI 1.01-1.36]; p < 0.04), and bipolar disorder (OR 1.21 [95% CI 1.05-1.40]; p < 0.007) and with an increase of 0.72 (95% CI 0.11-1.34; p < 0.02) mm Hg in systolic blood pressure, 0.45 (95% CI 0.06-0.84; p < 0.02) mm Hg in diastolic blood pressure, 0.01 ml/min/1.73 m2 (95% CI 0.003-0.02; p < 0.005) in estimated glomerular filtration rate from serum creatinine, 0.01 g/dl (95% CI 0.0004-0.02; p < 0.04) in serum albumin level, and 0.03 g/dl (95% CI 0.008-0.05; p < 0.009) in serum protein level. However, after adjustment for heterogeneity, neither GRS showed a significant effect of CRP level (at p < 0.0016) on any of these outcomes, including coronary artery disease, nor on the other 20 complex outcomes studied. Our study has two potential limitations: the limited variance explained by our genetic instruments modeling CRP levels in blood and the unobserved bias introduced by the use of summary statistics in our MR analyses.

CONCLUSIONS

Genetically elevated CRP levels showed a significant potentially protective causal relationship with risk of schizophrenia. We observed nominal evidence at an observed p < 0.05 using either GRSCRP or GRSGWAS-with persistence after correction for heterogeneity-for a causal relationship of elevated CRP levels with psoriatic osteoarthritis, rheumatoid arthritis, knee osteoarthritis, systolic blood pressure, diastolic blood pressure, serum albumin, and bipolar disorder. These associations remain yet to be confirmed. We cannot verify any causal effect of CRP level on any of the other common somatic and neuropsychiatric outcomes investigated in the present study. This implies that interventions that lower CRP level are unlikely to result in decreased risk for the majority of common complex outcomes.

The EGF receptor ligand amphiregulin controls cell division via FoxM1

Epidermal growth factor receptor (EGFR) is central to epithelial cell physiology, and deregulated EGFR signaling has an important role in a variety of human carcinomas. Here we show that silencing of the EGF-related factor amphiregulin (AREG) markedly inhibits the expansion of human keratinocytes through mitotic failure and accumulation of cells with ⩾ 4n DNA content. RNA-sequencing-based transcriptome analysis revealed that tetracycline-mediated AREG silencing significantly altered the expression of 2331 genes, 623 of which were not normalized by treatment with EGF. Interestingly, genes irreversibly upregulated by suppression of AREG overlapped with genes involved in keratinocyte differentiation. Moreover, a significant proportion of the irreversibly downregulated genes featured upstream binding sites recognized by forkhead box protein M1 (FoxM1), a key transcription factor in the control of mitosis that is widely dysregulated in cancer. The downregulation of FoxM1 and its target genes preceded mitotic arrest. Constitutive expression of FoxM1 in AREG knockdown cells normalized cell proliferation, reduced the number of cells with ⩾ 4n DNA content and rescued expression of FoxM1 target genes. These results demonstrate that AREG controls G2/M progression and cytokinesis in keratinocytes via activation of a FoxM1-dependent transcriptional program, suggesting new avenues for treatment of epithelial cancer.

Widespread non-additive and interaction effects within HLA loci modulate the risk of autoimmune diseases

Human leukocyte antigen (HLA) genes confer strong risk for autoimmune diseases on a log-additive scale. Here we speculated that differences in autoantigen binding repertoires between a heterozygote’s two expressed HLA variants may result in additional non-additive risk effects. We tested non-additive disease contributions of classical HLA alleles in patients and matched controls for five common autoimmune diseases: rheumatoid arthritis (RA, N_cases=5,337), type 1 diabetes (T1D, N_cases=5,567), psoriasis vulgaris (N_cases=3,089), idiopathic achalasia (N_cases=727), and celiac disease (N_cases=11,115). In four out of five diseases, we observed highly significant non-additive dominance effects (RA: P=2.5×10¹²; T1D: P=2.4×10⁻¹⁰; psoriasis: P=5.9×10⁻⁶; celiac disease: P=1.2×10⁻⁸⁷). In three of these diseases, the dominance effects were explained by interactions between specific classical HLA alleles (RA: P=1.8×10⁻³; T1D: P=8.6×10²⁷; celiac disease: P=6.0×10⁻¹⁰⁰). These interactions generally increased disease risk and explained moderate but significant fractions of phenotypic variance (RA: 1.4%, T1D: 4.0%, and celiac disease: 4.1%, beyond a simple additive model).

Proteogenomic analysis of psoriasis reveals discordant and concordant changes in mRNA and protein abundance

BACKGROUND

Psoriasis is a chronic disease characterized by the development of scaly red skin lesions and possible co-morbid conditions. The psoriasis lesional skin transcriptome has been extensively investigated, but mRNA levels do not necessarily reflect protein abundance. The purpose of this study was therefore to compare differential expression patterns of mRNA and protein in psoriasis lesions.

METHODS

Lesional (PP) and uninvolved (PN) skin samples from 14 patients were analyzed using high-throughput complementary DNA sequencing (RNA-seq) and liquid chromatography-tandem mass spectrometry (LC-MS/MS).

RESULTS

We identified 4122 differentially expressed genes (DEGs) along with 748 differentially expressed proteins (DEPs). Global shifts in mRNA were modestly correlated with changes in protein abundance (r = 0.40). We identified similar numbers of increased and decreased DEGs, but 4-fold more increased than decreased DEPs. Ribosomal subunit and translation proteins were elevated within lesions, without a corresponding shift in mRNA expression (RPL3, RPS8, RPL11). We identified 209 differentially expressed genes/proteins (DEGPs) with corresponding trends at the transcriptome and proteome levels. Most DEGPs were similarly altered in at least one other skin disease. Psoriasis-specific and non-specific DEGPs had distinct cytokine-response patterns, with only the former showing disproportionate induction by IL-17A in cultured keratinocytes.

CONCLUSIONS

Our findings reveal global imbalance between the number of increased and decreased proteins in psoriasis lesions, consistent with heightened translation. This effect could not have been discerned from mRNA profiling data alone. High-confidence DEGPs were identified through transcriptome-proteome integration. By distinguishing between psoriasis-specific and non-specific DEGPs, our analysis uncovered new functional insights that would otherwise have been overlooked.

Enhanced meta-analysis and replication studies identify five new psoriasis susceptibility loci

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.

Association of CLEC16A with human common variable immunodeficiency disorder and role in murine B cells

Common variable immunodeficiency disorder (CVID) is the most common symptomatic primary immunodeficiency in adults, characterized by B-cell abnormalities and inadequate antibody response. CVID patients have considerable autoimmune comorbidity and we therefore hypothesized that genetic susceptibility to CVID may overlap with autoimmune disorders. Here, in the largest genetic study performed in CVID to date, we compare 778 CVID cases with 10,999 controls across 123,127 single-nucleotide polymorphisms (SNPs) on the Immunochip. We identify the first non-HLA genome-wide significant risk locus at CLEC16A (rs17806056, P=2.0 × 10(-9)) and confirm the previously reported human leukocyte antigen (HLA) associations on chromosome 6p21 (rs1049225, P=4.8 × 10(-16)). Clec16a knockdown (KD) mice showed reduced number of B cells and elevated IgM levels compared with controls, suggesting that CLEC16A may be involved in immune regulatory pathways of relevance to CVID. In conclusion, the CLEC16A associations in CVID represent the first robust evidence of non-HLA associations in this immunodeficiency condition.

Analysis of long non-coding RNAs highlights tissue-specific expression patterns and epigenetic profiles in normal and psoriatic skin

BACKGROUND

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.

RESULTS

We used RNA-seq data from 99 lesional psoriatic, 27 uninvolved psoriatic, and 90 normal skin biopsies, and applied computational approaches to identify and characterize expressed lncRNAs. We detect 2,942 previously annotated and 1,080 novel lncRNAs which are expected to be skin specific. Notably, over 40% of the novel lncRNAs are differentially expressed and the proportions of differentially expressed transcripts among protein-coding mRNAs and previously-annotated lncRNAs are lower in psoriasis lesions versus uninvolved or normal skin. We find that many lncRNAs, in particular those that are differentially expressed, are co-expressed with genes involved in immune related functions, and that novel lncRNAs are enriched for localization in the epidermal differentiation complex. We also identify distinct tissue-specific expression patterns and epigenetic profiles for novel lncRNAs, some of which are shown to be regulated by cytokine treatment in cultured human keratinocytes.

CONCLUSIONS

Together, our results implicate many lncRNAs in the immunopathogenesis of psoriasis, and our results provide a resource for lncRNA studies in other autoimmune diseases.

Transcriptome analysis of psoriasis in a large case-control sample: RNA-seq provides insights into disease mechanisms

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.

Genetic associations of psoriasis in a Pakistani population

BACKGROUND

Genetic predisposition to psoriasis, an inflammatory skin disease affecting 0·2-4% of the world population, is well established. Thus far, 41 psoriasis susceptibility loci reach genome-wide significance (P ≤ 5 × 10(-8) ). Identification of genetic susceptibility loci in diverse populations will help understand the underlying biology of psoriasis susceptibility.

OBJECTIVES

The primary objective of this study was to examine psoriasis susceptibility associations previously reported in Chinese and caucasian populations in a Pakistani cohort.

METHODS

Blood samples and phenotype data were collected from psoriasis cases and controls in Islamabad, Pakistan. DNA was isolated and genotypes of selected susceptibility markers were determined. The data were analysed using χ(2) tests or logistic regression for psoriasis association.

RESULTS

HLA-Cw6 showed the strongest association [odds ratio (OR) 2·43, P = 2·3 × 10(-12) ]. HLA-Cw1 showed marginally significant association (OR 1·66, P = 0·049), suggesting that the HLA-Cw1-B46 risk haplotype may be present in the Pakistani population. Three other loci (IL4/IL13, NOS2, TRAF3IP2) showed nominally significant association (P < 0·05).

CONCLUSIONS

HLA-Cw6 is strongly associated with psoriasis susceptibility in the Pakistani population, as has been found in every other population studied. In addition, HLA-Cw1 showed marginal association, reflecting the relative geographical proximity and thus likely genetic relatedness to other populations in which the HLA-Cw1-B46 haplotype is known to be associated. A larger cohort and a denser marker set will be required for further analysis of psoriasis associations in the South Asian population.

Susceptibility-associated genetic variation at IL12B enhances Th1 polarization in psoriasis

The IL12B gene encodes the common p40 subunit of IL-12 and IL-23, cytokines with key roles in Th1 and Th17 biology, respectively, and genetic variation in this region significantly influences risk of psoriasis. Here, we demonstrate that a psoriasis-associated risk haplotype at the IL12B locus leads to increased expression of IL12B by monocytes and correlated with increased serum levels of IL-12, IFN-γ and the IFN-γ induced chemokine, CXCL10. In contrast, serum IL-23 levels were decreased in risk carriers when compared with non-carriers. We further demonstrate that IL-12 is increased in psoriatic skin and that risk carriers manifest a skewing of the inflammatory network toward stronger IFN-γ responses. Taken together, our data demonstrate that the risk variant in IL12B associates with its increased expression and predisposes to stronger Th1 polarization through deviation of the local inflammatory environment toward increased IL-12/IFN-γ at the expense of IL-23/IL-17 responses.

Identification of 15 new psoriasis susceptibility loci highlights the role of innate immunity

To gain further insight into the genetic architecture of psoriasis, we conducted a meta-analysis of 3 genome-wide association studies (GWAS) and 2 independent data sets genotyped on the Immunochip, including 10,588 cases and 22,806 controls. We identified 15 new susceptibility loci, increasing to 36 the number associated with psoriasis in European individuals. We also identified, using conditional analyses, five independent signals within previously known loci. The newly identified loci shared with other autoimmune diseases include candidate genes with roles in regulating T-cell function (such as RUNX3, TAGAP and STAT3). Notably, they included candidate genes whose products are involved in innate host defense, including interferon-mediated antiviral responses (DDX58), macrophage activation (ZC3H12C) and nuclear factor (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.

Rare and common variants in CARD14, encoding an epidermal regulator of NF-kappaB, in psoriasis

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.

Combined analysis of genome-wide association studies for Crohn disease and psoriasis identifies seven shared susceptibility loci

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 × 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 × 10(-8), 11q13 near PRDX5, p(rs694739) = 3.71 × 10(-09), 22q11 at YDJC, p(rs181359) = 8.02 × 10(-10)). In addition, we identified one susceptibility locus for CD (16p13 near SOCS1, p(rs4780355) = 4.99 × 10(-8)). Refinement of association signals identified shared genome-wide 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.

Heterogeneity of inflammatory and cytokine networks in chronic plaque psoriasis

The clinical features of psoriasis, characterized by sharply demarcated scaly erythematous plaques, are typically so distinctive that a diagnosis can easily be made on these grounds alone. However, there is great variability in treatment response between individual patients, and this may reflect heterogeneity of inflammatory networks driving the disease. In this study, whole-genome transcriptional profiling was used to characterize inflammatory and cytokine networks in 62 lesional skin samples obtained from patients with stable chronic plaque psoriasis. We were able to stratify lesions according to their inflammatory gene expression signatures, identifying those associated with strong (37% of patients), moderate (39%) and weak inflammatory infiltrates (24%). Additionally, we identified differences in cytokine signatures with heightened cytokine-response patterns in one sub-group of lesions (IL-13-strong; 50%) and attenuation of these patterns in a second sub-group (IL-13-weak; 50%). These sub-groups correlated with the composition of the inflammatory infiltrate, but were only weakly associated with increased risk allele frequency at some psoriasis susceptibility loci (e.g., REL, TRAF3IP2 and NOS2). Our findings highlight variable points in the inflammatory and cytokine networks known to drive chronic plaque psoriasis. Such heterogeneous aspects may shape clinical course and treatment responses, and can provide avenues for development of personalized treatments.

Genome-wide meta-analysis of psoriatic arthritis identifies susceptibility locus at REL

Psoriatic arthritis (PsA) is a chronic inflammatory musculoskeletal disease affecting up to 30% of psoriasis vulgaris (PsV) cases and approximately 0.25% to 1% of the general population. To identify common susceptibility loci, we performed a meta-analysis of three imputed genome-wide association studies (GWAS) on psoriasis, stratified for PsA. A total of 1,160,703 SNPs were analyzed in the discovery set consisting of 535 PsA cases and 3,432 controls from Germany, the United States and Canada. We followed up two SNPs in 1,931 PsA cases and 6,785 controls comprising six independent replication panels from Germany, Estonia, the United States and Canada. In the combined analysis, a genome-wide significant association was detected at 2p16 near the REL locus encoding c-Rel (rs13017599, P=1.18×10⁻⁸, OR=1.27, 95% CI=1.18–1.35). The rs13017599 polymorphism is known to associate with rheumatoid arthritis (RA), and another SNP near REL (rs702873) was recently implicated in PsV susceptibility. However, conditional analysis indicated that rs13017599, rather than rs702873, accounts for the PsA association at REL. We hypothesize that c-Rel, as a member of the Rel/NF-κB family, is associated with PsA in the context of disease pathways that involve other identified PsA and PsV susceptibility genes including TNIP1, TNFAIP3 and NFκBIA.

Genome-wide expression profiling of five mouse models identifies similarities and differences with human psoriasis

Development of a suitable mouse model would facilitate the investigation of pathomechanisms underlying human psoriasis and would also assist in development of therapeutic treatments.

However, while many psoriasis mouse models have been proposed, no single model recapitulates all features of the human disease, and standardized validation criteria for psoriasis mouse models have not been widely applied. In this study, whole-genome transcriptional profiling is used to compare gene expression patterns manifested by human psoriatic skin lesions with those that occur in five psoriasis mouse models (K5-Tie2, imiquimod, K14-AREG, K5-Stat3C and K5-TGFbeta1). While the cutaneous gene expression profiles associated with each mouse phenotype exhibited statistically significant similarity to the expression profile of psoriasis in humans, each model displayed distinctive sets of similarities and differences in comparison to human psoriasis. For all five models, correspondence to the human disease was strong with respect to genes involved in epidermal development and keratinization. Immune and inflammation-associated gene expression, in contrast, was more variable between models as compared to the human disease. These findings support the value of all five models as research tools, each with identifiable areas of convergence to and divergence from the human disease. Additionally, the approach used in this paper provides an objective and quantitative method for evaluation of proposed mouse models of psoriasis, which can be strategically applied in future studies to score strengths of mouse phenotypes relative to specific aspects of human psoriasis.

Gene expression in skin and lymphoblastoid cells: Refined statistical method reveals extensive overlap in cis-eQTL signals

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.

Genome-wide association analysis identifies three psoriasis susceptibility loci

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⁻⁷).

Genome-wide association study identifies a psoriasis susceptibility locus at TRAF3IP2

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 PsV cases 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; combined P = 2.36 × 10⁻¹⁰ for rs13210247 and combined P = 1.24 × 10⁻¹⁶ for rs33980500). About 15% of psoriasis cases develop 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⁻¹² 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 members of the Rel/NF-κB transcription factor family.

Association analyses identify six new psoriasis susceptibility loci in the Chinese population

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

Comparison of MHC class I risk haplotypes in Thai and Caucasian psoriatics shows locus heterogeneity at PSORS1

Earlier studies have shown that psoriasis in Japan and Thailand is associated with two different major histocompatibility complex (MHC) haplotypes - those bearing HLA-Cw6 and those bearing HLA-Cw1 and HLA-B46. In an independent case-control sample from Thailand, we confirmed the association of psoriasis with both haplotypes. No association was seen in Thai HLA-Cw1 haplotypes lacking HLA-B46, nor was HLA-Cw1 associated with psoriasis in a large Caucasian sample. To assess whether these risk haplotypes share a common origin, we sequenced genomic DNA from a Thai HLA-Cw1-B46 homozygote across the ∼300 kb MHC risk interval, and compared it with sequence of a HLA-Cw6-B57 risk haplotype. Three small regions of homology were found, but these regions share equivalent sequence similarity with one or more clearly non-risk haplotypes, and they contain no polymorphism alleles unique to all risk haplotypes. Differences in psoriasis phenotype were also observed, including lower risk of disease, greater nail involvement, and later age at onset in HLA-Cw1-B46 carriers compared with HLA-Cw6 carriers. These findings suggest locus heterogeneity at PSORS1 (psoriasis susceptibility 1), the major psoriasis susceptibility locus in the MHC, with HLA-Cw6 imparting risk in both Caucasians and Asians, and an allele other than HLA-Cw1 on the HLA-Cw1-B46 haplotype acting as an additional risk variant in East Asians.

Assessment of the psoriatic transcriptome in a large sample: additional regulated genes and comparisons with in vitro models

To further elucidate molecular alterations in psoriasis, we performed a gene expression study of 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 (549 up- and 369 downregulated), of which over 600 are to our knowledge previously unreported, including S100A7A, THRSP, and ELOVL3. Strongly upregulated genes included SERPINB4, PI3, DEFB4, and several S100-family members. Strongly downregulated genes included Wnt-inhibitory factor-1 (WIF1), beta-cellulin (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-1alpha, IFN-gamma, TNF-alpha, and OSM) showed 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.

Evidence for altered Wnt signaling in psoriatic skin

The Wnt gene family encodes a set of highly conserved secreted signaling proteins that have major roles in embryogenesis and tissue homeostasis. Yet the expression of this family of important mediators in psoriasis, a disease characterized by marked changes in keratinocyte growth and differentiation, is incompletely understood. We subjected 58 paired biopsies from lesional and uninvolved psoriatic skin and 64 biopsies from normal skin to global gene expression profiling. WNT5A transcripts were upregulated fivefold in lesional skin, accompanied by increased Wnt-5a protein levels. Notably, WNT5A mRNA was markedly induced by IL-1alpha, tumor necrosis factor-alpha, IFN-gamma, and transforming growth factor-alpha in cultured keratinocytes. Frizzled 2 (FZD2) and FZD5, which encode receptors for Wnt5A, were also increased in lesional psoriatic skin. In contrast, expression of WIF1 mRNA, encoding a secreted antagonist of the Wnt proteins, was downregulated >10-fold in lesional skin, along with decreased WNT inhibitory factor (WIF)-1 immunostaining. Interestingly, pathway analysis along with reduced AXIN2 expression and lack of nuclear translocation of beta-catenin indicated a suppression of canonical Wnt signaling in lesional skin. The results of our study suggest a shift away from canonical Wnt signaling toward noncanonical pathways driven by interactions between Wnt-5a and its cognate receptors in psoriasis, accompanied by impaired homeostatic inhibition of Wnt signaling by WIF-1 and dickkopf.

Molecular dissection of psoriasis: integrating genetics and biology

Psoriasis is a common and debilitating disease of the skin, nails, and joints, with an acknowledged but complex genetic basis. Early genome-wide linkage studies of psoriasis focused on segregation of microsatellite markers in families; however, the only locus consistently identified resided in the major histocompatibility complex. Subsequently, several groups mapped this locus to the vicinity of HLA-C, and two groups have reported HLA-Cw6 itself to be the major susceptibility allele. More recently, the development of millions of single-nucleotide polymorphisms, coupled with the development of high-throughput genotyping platforms and a comprehensive map of human haplotypes, has made possible a genome-wide association approach using cases and controls rather than families. Taking advantage of these developments, we participated in a collaborative genome-wide association study of psoriasis involving thousands of cases and controls. Initial analysis of these data revealed and/or confirmed association between psoriasis and seven genetic loci-HLA-C, IL12B, IL23R, IL23A, IL4/IL13, TNFAIP3, and TNIP1-and ongoing studies are revealing additional loci. Here, we review the epidemiology, immunopathology, and genetics of psoriasis, and present a disease model integrating its genetics and immunology.

Genetic evidence for involvement of the IL23 pathway in Thai psoriatics

A recent genome-wide association analysis of psoriasis identified IL12B and IL23R as significantly associated with psoriasis. Here we report association test results of a Thai cohort consisting of 206 psoriasis cases and 114 controls. The IL23R SNPs rs7530511 and rs11209026, and IL12B SNPs rs3212227 and rs6887695 were genotyped using Taqman assays. Data were analyzed using a logistic regression model for linear trend of association. One of the IL23R markers, rs7530511, was marginally significant (P = 0.017). The other IL23R marker, rs11209026, was not polymorphic. One of the IL12B markers, rs3212227, showed significant association with psoriasis (OR = 1.64, P = 0.0058) while the other, rs6887695, did not (OR = 1.29, P = 0.12). Haplotype analysis of the two IL12B SNPs yielded highly significant association (P = 0.00081, OR = 1.73). These results showed that IL12B is an important genetic factor in psoriasis pathogenesis in the Thai population, with an association strong enough to yield significant confirmatory evidence using a modest sample size. Together with previously reported evidence for IL12B association in Caucasian, Japanese, and Chinese psoriatics, our results support the hypothesis that genes encoding components of the IL23-mediated inflammatory pathway are important determinants of psoriasis pathogenesis across multiple racial groups.

Multiple Loci within the major histocompatibility complex confer risk of psoriasis

Psoriasis is a common inflammatory skin disease characterized by thickened scaly red plaques. Previously we have performed a genome-wide association study (GWAS) on psoriasis with 1,359 cases and 1,400 controls, which were genotyped for 447,249 SNPs. The most significant finding was for SNP rs12191877, which is in tight linkage disequilibrium with HLA-Cw*0602, the consensus risk allele for psoriasis. However, it is not known whether there are other psoriasis loci within the MHC in addition to HLA-C. In the present study, we searched for additional susceptibility loci within the human leukocyte antigen (HLA) region through in-depth analyses of the GWAS data; then, we followed up our findings in an independent Han Chinese 1,139 psoriasis cases and 1,132 controls. Using the phased CEPH dataset as a reference, we imputed the HLA-Cw*0602 in all samples with high accuracy. The association of the imputed HLA-Cw*0602 dosage with disease was much stronger than that of the most significantly associated SNP, rs12191877. Adjusting for HLA-Cw*0602, there were two remaining association signals: one demonstrated by rs2073048 (p = 2×10⁻⁶, OR = 0.66), located within c6orf10, a potential downstream effecter of TNF-alpha, and one indicated by rs13437088 (p = 9×10⁻⁶, OR = 1.3), located 30 kb centromeric of HLA-B and 16 kb telomeric of MICA. When HLA-Cw*0602, rs2073048, and rs13437088 were all included in a logistic regression model, each of them was significantly associated with disease (p = 3×10⁻⁴⁷, 6×10⁻⁸, and 3×10⁻⁷, respectively). Both putative loci were also significantly associated in the Han Chinese samples after controlling for the imputed HLA-Cw*0602. A detailed analysis of HLA-B in both populations demonstrated that HLA-B*57 was associated with an increased risk of psoriasis and HLA-B*40 a decreased risk, independently of HLA-Cw*0602 and the C6orf10 locus, suggesting the potential pathogenic involvement of HLA-B. These results demonstrate that there are at least two additional loci within the MHC conferring risk of psoriasis.

Psoriasis (Ps) is a chronic inflammatory disease of the skin, affecting approximately 2% of Europeans. The HLA-C gene, located within the major histocompatibility complex (MHC) region on chromosome 6, is the major genetic determinant of psoriasis. However, multiple susceptibility genes within MHC are also hypothesized. Recently, we carried out a genome-wide association scan (GWAS) on psoriasis with 1,359 patients and 1,400 healthy controls, which identified seven psoriasis loci in the human genome and confirmed the effect of HLA-C. This dataset contains densely distributed genetic variations, single nucleotide polymorphisms (SNPs), which were then further analyzed in search for additional susceptibility genes within the MHC region. Using the SNP data, we imputed in all samples the HLA-C risk allele with high accuracy. Adjusting for the HLA-C, two additional loci, one near C6orf10 and one near HLA-B/MICA, have significant associations with psoriasis, which were also observed in an independent Han Chinese dataset, suggesting that within the MHC there are at least three genes moderating susceptibility to psoriasis.