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.

Disease consequences of higher adiposity uncoupled from its adverse metabolic effects using Mendelian randomisation

Background

Some individuals living with obesity may be relatively metabolically healthy, whilst others suffer from multiple conditions that may be linked to adverse metabolic effects or other factors. The extent to which the adverse metabolic component of obesity contributes to disease compared to the non-metabolic components is often uncertain. We aimed to use Mendelian randomisation (MR) and specific genetic variants to separately test the causal roles of higher adiposity with and without its adverse metabolic effects on diseases.

Methods

We selected 37 chronic diseases associated with obesity and genetic variants associated with different aspects of excess weight. These genetic variants included those associated with metabolically 'favourable adiposity' (FA) and 'unfavourable adiposity' (UFA) that are both associated with higher adiposity but with opposite effects on metabolic risk. We used these variants and two sample MR to test the effects on the chronic diseases.

Results

MR identified two sets of diseases. First, 11 conditions where the metabolic effect of higher adiposity is the likely primary cause of the disease. Here, MR with the FA and UFA genetics showed opposing effects on risk of disease: coronary artery disease, peripheral artery disease, hypertension, stroke, type 2 diabetes, polycystic ovary syndrome, heart failure, atrial fibrillation, chronic kidney disease, renal cancer, and gout. Second, 9 conditions where the non-metabolic effects of excess weight (e.g. mechanical effect) are likely a cause. Here, MR with the FA genetics, despite leading to lower metabolic risk, and MR with the UFA genetics, both indicated higher disease risk: osteoarthritis, rheumatoid arthritis, osteoporosis, gastro-oesophageal reflux disease, gallstones, adult-onset asthma, psoriasis, deep vein thrombosis, and venous thromboembolism.

Conclusions

Our results assist in understanding the consequences of higher adiposity uncoupled from its adverse metabolic effects, including the risks to individuals with high body mass index who may be relatively metabolically healthy.

Funding

Diabetes UK, UK Medical Research Council, World Cancer Research Fund, National Cancer Institute.

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.

A high-resolution HLA reference panel capturing global population diversity enables multi-ancestry fine-mapping in HIV host response

Fine-mapping to plausible causal variation may be more effective in multi-ancestry cohorts, particularly in the MHC, which has population-specific structure. To enable such studies, we constructed a large (n = 21,546) HLA reference panel spanning five global populations based on whole-genome sequences. Despite population-specific long-range haplotypes, we demonstrated accurate imputation at G-group resolution (94.2%, 93.7%, 97.8% and 93.7% in admixed African (AA), East Asian (EAS), European (EUR) and Latino (LAT) populations). Applying HLA imputation to genome-wide association study data for HIV-1 viral load in three populations (EUR, AA and LAT), we obviated effects of previously reported associations from population-specific HIV studies and discovered a novel association at position 156 in HLA-B. We pinpointed the MHC association to three amino acid positions (97, 67 and 156) marking three consecutive pockets (C, B and D) within the HLA-B peptide-binding groove, explaining 12.9% of trait variance.

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.

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.

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.

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.

A rare coding allele in IFIH1 is protective for psoriatic arthritis

Objectives

Psoriatic arthritis (PsA) is an inflammatory arthritis associated with psoriasis. While many common risk alleles have been reported for association with PsA as well as psoriasis, few rare coding alleles have yet been identified.

Methods

To identify rare coding variation associated with PsA risk or protection, we genotyped 41 267 variants with the exome chip and investigated association within an initial cohort of 1980 PsA cases and 5913 controls. Genotype data for an independent cohort of 2234 PsA cases and 5708 controls was also made available, allowing for a meta-analysis to be performed with the discovery dataset.

Results

We identified an association with the rare variant rs35667974 (p=2.39x10⁻⁶, OR=0.47), encoding an Ile923Val amino acid change in the IFIH1 gene protein product. The association was reproduced in our independent cohort, which reached a high level of significance on meta-analysis with the discovery and replication datasets (p=4.67x10⁻¹⁰). We identified a strong association with IFIH1 when performing multiple-variant analysis (p=6.77x10⁻⁶), and found evidence of independent effects between the rare allele and the common PsA variant at the same locus.

Conclusion

For the first time, we report a rare coding allele in IFIH1 to be protective for PsA. This rare allele has also been identified to have the same direction of effect on type I diabetes and psoriasis. While this association further supports existing evidence for IFIH1 as a causal gene for PsA, mechanistic studies will need to be pursued to confirm that IFIH1 is indeed causal.

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.

Membrane-Tethered Intracellular Domain of Amphiregulin Promotes Keratinocyte Proliferation

The epidermal growth factor receptor (EGFR) and its ligands are essential regulators of epithelial biology, which are often amplified in cancer cells. We have previously shown that shRNA-mediated silencing of one of these ligands, amphiregulin (AREG), results in keratinocyte growth arrest that cannot be rescued by soluble extracellular EGFR ligands. To further explore the functional importance of specific AREG domains, we stably transduced keratinocytes expressing tetracycline-inducible AREG-targeted shRNA with lentiviruses expressing silencing-proof, membrane-tethered AREG cytoplasmic and extracellular domains (AREG-CTD and AREG-ECD), as well as full-length AREG precursor (proAREG). Here we show that growth arrest of AREG-silenced keratinocytes occurs in G2/M and is significantly restored by proAREG and AREG-CTD but not by AREG-ECD. Moreover, the AREG-CTD was sufficient to normalize cell cycle distribution profiles and expression of mitosis-related genes. Our findings uncover an important role of the AREG-CTD in regulating cell division, which may be relevant to tumor resistance to EGFR-directed therapies.

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

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.

Psoriasis drug development and GWAS interpretation through in silico analysis of transcription factor binding sites

Background

Psoriasis is a cytokine-mediated skin disease that can be treated effectively with immunosuppressive biologic agents. These medications, however, are not equally effective in all patients and are poorly suited for treating mild psoriasis. To develop more targeted therapies, interfering with transcription factor (TF) activity is a promising strategy.

Methods

Meta-analysis was used to identify differentially expressed genes (DEGs) in the lesional skin from psoriasis patients (n = 237). We compiled a dictionary of 2935 binding sites representing empirically-determined binding affinities of TFs and unconventional DNA-binding proteins (uDBPs). This dictionary was screened to identify “psoriasis response elements” (PREs) overrepresented in sequences upstream of psoriasis DEGs.

Results

PREs are recognized by IRF1, ISGF3, NF-kappaB and multiple TFs with helix-turn-helix (homeo) or other all-alpha-helical (high-mobility group) DNA-binding domains. We identified a limited set of DEGs that encode proteins interacting with PRE motifs, including TFs (GATA3, EHF, FOXM1, SOX5) and uDBPs (AVEN, RBM8A, GPAM, WISP2). PREs were prominent within enhancer regions near cytokine-encoding DEGs (IL17A, IL19 and IL1B), suggesting that PREs might be incorporated into complex decoy oligonucleotides (cdODNs). To illustrate this idea, we designed a cdODN to concomitantly target psoriasis-activated TFs (i.e., FOXM1, ISGF3, IRF1 and NF-kappaB). Finally, we screened psoriasis-associated SNPs to identify risk alleles that disrupt or engender PRE motifs. This identified possible sites of allele-specific TF/uDBP binding and showed that PREs are disproportionately disrupted by psoriasis risk alleles.

Conclusions

We identified new TF/uDBP candidates and developed an approach that (i) connects transcriptome informatics to cdODN drug development and (ii) enhances our ability to interpret GWAS findings. Disruption of PRE motifs by psoriasis risk alleles may contribute to disease susceptibility.

Electronic supplementary material

The online version of this article (doi:10.1186/s40169-015-0054-5) contains supplementary material, which is available to authorized users.

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.

Cellular dissection of psoriasis for transcriptome analyses and the post-GWAS era

Background

Genome-scale studies of psoriasis have been used to identify genes of potential relevance to disease mechanisms. For many identified genes, however, the cell type mediating disease activity is uncertain, which has limited our ability to design gene functional studies based on genomic findings.

Methods

We identified differentially expressed genes (DEGs) with altered expression in psoriasis lesions (n = 216 patients), as well as candidate genes near susceptibility loci from psoriasis GWAS studies. These gene sets were characterized based upon their expression across 10 cell types present in psoriasis lesions. Susceptibility-associated variation at intergenic (non-coding) loci was evaluated to identify sites of allele-specific transcription factor binding.

Results

Half of DEGs showed highest expression in skin cells, although the dominant cell type differed between psoriasis-increased DEGs (keratinocytes, 35%) and psoriasis-decreased DEGs (fibroblasts, 33%). In contrast, psoriasis GWAS candidates tended to have highest expression in immune cells (71%), with a significant fraction showing maximal expression in neutrophils (24%, P < 0.001). By identifying candidate cell types for genes near susceptibility loci, we could identify and prioritize SNPs at which susceptibility variants are predicted to influence transcription factor binding. This led to the identification of potentially causal (non-coding) SNPs for which susceptibility variants influence binding of AP-1, NF-κB, IRF1, STAT3 and STAT4.

Conclusions

These findings underscore the role of innate immunity in psoriasis and highlight neutrophils as a cell type linked with pathogenetic mechanisms. Assignment of candidate cell types to genes emerging from GWAS studies provides a first step towards functional analysis, and we have proposed an approach for generating hypotheses to explain GWAS hits at intergenic loci.

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.

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.

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.

Polymorphisms of the IL12B and IL23R genes are associated with psoriasis

Psoriasis is a common inflammatory and hyperproliferative skin disease with a multifactorial genetic basis. A recent study reported that psoriasis was associated with the IL12B haplotype rs3212227 (3'-untranslated region)-rs6887695 (60 kb, 5') and the IL23R haplotype rs7530511 (L310P)-rs11209026 (Q381R). We examined these four single-nucleotide polymorphisms (SNPs) for association with psoriasis in two groups of North American and German Caucasians: (1) 1,810 cases and 2,522 controls; and (2) 509 pedigrees. Both IL12B markers showed highly significant association with psoriasis in the case-control (rs3212227, odds ratio (OR)=1.62, P=1.7 x 10(-15); rs6887695, OR=1.49, P=2.7 x 10(-15)) and in the family-based analysis (rs3212227, P=2.2 x 10(-3); rs6887695, P=1.7 x 10(-3)). The IL23R SNPs also showed significant association in the cases and controls (rs7530511, OR=1.22, P=3.9 x 10(-3); rs11209026, OR=1.40, P=3.8 x 10(-4)). For both genes, common risk haplotypes were identified whose statistical significance approached (IL23R) or exceeded (IL12B) genome-wide criteria. We found no statistical evidence for interactions of these haplotypes with HLA-Cw6. Our results confirm associations between IL12B and IL23R and psoriasis in Caucasians, and provide a genetic basis for the clinical association between psoriasis and Crohn's disease.

Sequence and haplotype analysis supports HLA-C as the psoriasis susceptibility 1 gene

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.

Visualizing human leukocyte antigen class II risk haplotypes in human systemic lupus erythematosus

Human leukocyte antigen (HLA) class I and class II alleles are implicated as genetic risk factors for many autoimmune diseases. However, the role of the HLA loci in human systemic lupus erythematosus (SLE) remains unclear. Using a dense map of polymorphic microsatellites across the HLA region in a large collection of families with SLE, we identified three distinct haplotypes that encompassed the class II region and exhibited transmission distortion. DRB1 and DQB1 typing of founders showed that the three haplotypes contained DRB1*1501/ DQB1*0602, DRB1*0801/ DQB1*0402, and DRB1*0301/DQB1*0201 alleles, respectively. By visualizing ancestral recombinants, we narrowed the disease-associated haplotypes containing DRB1*1501 and DRB1*0801 to an approximately 500-kb region. We conclude that HLA class II haplotypes containing DRB1 and DQB1 alleles are strong risk factors for human SLE.

Genetic linkage and transmission disequilibrium of marker haplotypes at chromosome 1q41 in human systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease characterized by the production of autoantibodies to a wide range of self-antigens. Recent genome screens have implicated numerous chromosomal regions as potential SLE susceptibility loci. Among these, the 1q41 locus is of particular interest, because evidence for linkage has been found in several independent SLE family collections. Additionally, the 1q41 locus appears to be syntenic with a susceptibility interval identified in the NZM2410 mouse model for SLE. Here, we report the results of genotyping of 11 microsatellite markers within the 1q41 region in 210 SLE sibpair and 122 SLE trio families. These data confirm the modest evidence for linkage at 1q41 in our family collection (LOD = 1.21 at marker D1S2616). Evidence for significant linkage disequilibrium in this interval was also found. Multiple markers in the region exhibit transmission disequilibrium, with the peak single marker multiallelic linkage disequilibrium noted at D1S490 (pedigree disequilibrium test [PDT] global P value = 0.0091). Two- and three-marker haplotypes from the 1q41 region similarly showed strong transmission distortion in the collection of 332 SLE families. The finding of linkage together with significant transmission disequilibrium provides strong evidence for a susceptibility locus at 1q41 in human SLE.