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.

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.

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.

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.

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.

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.

Genome-wide Association Analysis of Psoriatic Arthritis and Cutaneous Psoriasis Reveals Differences in Their Genetic Architecture.

Psoriasis vulgaris (PsV) is a common inflammatory and hyperproliferative skin disease. Up to 30% of people with PsV eventually develop psoriatic arthritis (PsA), an inflammatory musculoskeletal condition. To discern differences in genetic risk factors for PsA and cutaneous-only psoriasis (PsC), we carried out a genome-wide association study (GWAS) of 1,430 PsA case subjects and 1,417 unaffected control subjects. Meta-analysis of this study with three other GWASs and two targeted genotyping studies, encompassing a total of 9,293 PsV case subjects, 3,061 PsA case subjects, 3,110 PsC case subjects, and 13,670 unaffected control subjects of European descent, detected 10 regions associated with PsA and 11 with PsC at genome-wide (GW) significance. Several of these association signals (IFNLR1, IFIH1, NFKBIA for PsA; TNFRSF9, LCE3C/B, TRAF3IP2, IL23A, NFKBIA for PsC) have not previously achieved GW significance. After replication, we also identified a PsV-associated SNP near CDKAL1 (rs4712528, odds ratio [OR] = 1.16, p = 8.4 × 10(-11)). Among identified psoriasis risk variants, three were more strongly associated with PsC than PsA (rs12189871 near HLA-C, p = 5.0 × 10(-19); rs4908742 near TNFRSF9, p = 0.00020; rs10888503 near LCE3A, p = 0.0014), and two were more strongly associated with PsA than PsC (rs12044149 near IL23R, p = 0.00018; rs9321623 near TNFAIP3, p = 0.00022). The PsA-specific variants were independent of previously identified psoriasis variants near IL23R and TNFAIP3. We also found multiple independent susceptibility variants in the IL12B, NOS2, and IFIH1 regions. These results provide insights into the pathogenetic similarities and differences between PsC and PsA.

Fine mapping major histocompatibility complex associations in psoriasis and its clinical subtypes.

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

SLC44A2 single nucleotide polymorphisms, isoforms, and expression: Association with severity of Meniere's disease?

SLC44A2 was discovered as the target of an antibody that causes hearing loss. Knockout mice develop age related hearing loss, loss of sensory cells and spiral ganglion neurons. SLC44A2 has polymorphic sites implicated in human disease. Transfusion related acute lung injury (TRALI) is linked to rs2288904 and genome wide association studies link rs2288904 and rs9797861 to venous thromboembolism (VTE), coronary artery disease and stroke. Here we report linkage disequilibrium of rs2288904 with rs3087969 and the association of these SLC44A2 SNPs with Meniere's disease severity. Tissue-specific isoform expression differences suggest that the N-terminal domain is linked to different functions in different cell types. Heterozygosity at rs2288904 CGA/CAA and rs3087969 GAT/GAC showed a trend for association with intractable Meniere's disease compared to less severe disease and to controls. The association of SLC44A2 SNPs with VTE suggests that thrombi affecting cochlear vessels could be a factor in Meniere's disease.

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.

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.

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.

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.

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.

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 x 10(-6), OR = 0.66), located within c6orf10, a potential downstream effecter of TNF-alpha, and one indicated by rs13437088 (p = 9 x 10(-6), 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 x 10(-47), 6 x 10(-8), and 3 x 10(-7), 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.

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.

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.

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.