Search

S2. ASXL structure and evolution

S3. Domain architecture of human BAP1 and ASXL2 and the proteins/domains used in this study

S4. Refinement of the mechanism for ASXL stimulation of BAP deubiquitinase activity.

S1. BAP1 expression in mouse single cell datasets.

PDF file - 2.3MB

PDf file - 8.7MB

Purpose: Metastasis is responsible for the death of most cancer patients, yet few therapeutic agents are available which specifically target the molecular events that lead to metastasis. We recently showed that inactivating mutations in the tumor suppressor gene BAP1 are closely associated with loss of melanocytic differentiation in uveal melanoma (UM) and metastasis. The purpose of this study was to identify therapeutic agents that reverse the phenotypic effects of BAP1 loss in UM.Experimental Design:In silico screens were done to identify therapeutic compounds predicted to differentiate UM cells using Gene Set Enrichment Analysis and Connectivity Map databases. Valproic acid (VPA), trichostatin A, LBH-589, and suberoylanilide hydroxamic acid were evaluated for their effects on UM cells using morphologic evaluation, MTS viability assays, bromodeoxyuridine incorporation, flow cytometry, clonogenic assays, gene expression profiling, histone acetylation and ubiquitination assays, and a murine xenograft tumorigenicity model.Results: Histone deacetylase (HDAC) inhibitors induced morphologic differentiation, cell-cycle exit, and a shift to a differentiated, melanocytic gene expression profile in cultured UM cells. VPA inhibited the growth of UM tumors in vivo.Conclusions: These findings suggest that HDAC inhibitors may have therapeutic potential for inducing differentiation and prolonged dormancy of micrometastatic disease in UM. Clin Cancer Res; 18(2); 408–16. ©2011 AACR.

PDF file - 1.9MB

PDF file - 2.4MB

Purpose:Uveal melanoma is a primary malignancy of the eye with oncogenic mutations in GNAQ, GNA11, or CYSLTR2, and additional mutations in BAP1 (usually associated with LOH of Chr 3), SF3B1, or EIF1AX. There are other characteristic chromosomal alterations, but their significance is not clear.Experimental Design:To investigate genes driving chromosomal alterations, we integrated copy number, transcriptome, and mutation data from three cohorts and followed up key findings.Results:We observed significant enrichment of transcripts on chromosomes 1p, 3, 6, 8, and 16q and identified seven shared focal copy number alterations (FCNAs) on Chr 1p36, 2q37, 3, 6q25, 6q27, and 8q24. Integrated analyses revealed clusters of genes in focal copy number regions whose expression was associated with metastasis and worse overall survival. This included genes from Chr 1p36, 3p21, and 8q24.3. At Chr 6q27, we identified two tumors with homozygous deletion of PHF10/BAF45a and one with a frameshift mutation with concomitant loss of the wild-type allele. Downregulation of PHF10 in uveal melanoma cell lines and tumors altered a number of biological pathways including development and adhesion. These findings provide support for a role for PHF10 as a novel tumor suppressor at Chr 6q27.Conclusions:Integration of copy number, transcriptome, and mutation data revealed novel candidate genes playing a role in uveal melanoma pathogenesis and a potential tumor suppressor role for PHF10.

Purpose: Loss of chromosome 3 is strongly associated with metastasis in uveal melanoma and has been proposed as the basis for clinical prognostic testing. It is not known whether techniques that identify loss of heterozygosity for chromosome 3 predict metastasis more accurately than those that detect only numerical loss of chromosome 3 (monosomy 3).Experimental Design: Fifty-three uveal melanomas were analyzed by 28 single nucleotide polymorphisms (SNP) across chromosome 3. SNP was compared with fluorescence in situ hybridization (FISH) and array-based comparative genomic hybridization (aCGH) for metastasis prediction by sensitivity, specificity, and Kaplan-Meier survival analysis, using our validated gene expression-based classifier as a reference standard.Results: By Kaplan-Meier analysis, only the gene expression-based classifier (P = 0.001) and SNP-based detection of loss of heterozygosity for chromosome 3 (P = 0.04) were significantly associated with metastasis. Sensitivity and specificity were 95.2% and 80.8%, respectively, for SNP, 77.8% and 64.7%, respectively, for FISH, and 85.0% and 72.0%, respectively, for aCGH. Isodisomy 3 was identified by SNP but undetected by aCGH and FISH in three tumors.Conclusions: Prognostic tests based on SNP platforms, which detect both chromosomal homologues and their subregions, may be superior to techniques that only detect changes in chromosome number. These observations could have important implications for efforts to detect genetic alterations in cancer genomes with CGH-based approaches.

Supplementary Tables

PDF file - 65K

Supplementary Table S1 from Loss of Heterozygosity of Chromosome 3 Detected with Single Nucleotide Polymorphisms Is Superior to Monosomy 3 for Predicting Metastasis in Uveal Melanoma

Introduced mutations in BAP1 that are predicted to disrupt the helix structure of the BAP1-ULD domain dramatically disassociate its binding to the ASXL2-AB box. The GST-UCH association assay was performed with highly purified recombinant ULD-WT, L-P, or F-P mutant proteins and/or ASXL2-AB box proteins.

Deleterious mutations of the ubiquitin carboxy-terminal hydrolase BAP1 found in cancers are predicted to encode inactive truncated proteins, suggesting that loss of enzyme function is a primary tumorigenic mechanism. However, many tumors exhibit missense mutations or in-frame deletions or insertions, often outside the functionally critical UCH domain in this tumor suppressor protein. Thus, precisely how these mutations inactivate BAP1 is unknown. Here, we show how these mutations affect BAP1 interactions with the Polycomb group-like protein, ASXL2, using combinations of computational modeling technology, molecular biology, and in vitro reconstitution biochemistry. We found that the BAP1–ASXL2 interaction is direct and high affinity, occurring through the ASXH domain of ASXL2, an obligate partner for BAP1 enzymatic activity. The ASXH domain was the minimal domain for binding the BAP1 ULD domain, and mutations on the surfaces of predicted helices of ASXH abolished BAP1 association and stimulation of BAP1 enzymatic activity. The BAP1-UCH, BAP1-ULD, and ASXH domains formed a cooperative stable ternary complex required for deubiquitination. We defined four classes of alterations in BAP1 outside the UCH domain, each failing to productively recruit ASXH to the wild-type BAP1 catalytic site via the ULD, resulting in loss of BAP1 ubiquitin hydrolase activity. Our results indicate that many BAP1 mutations act allosterically to inhibit ASXH binding, thereby leading to loss of enzyme activity. Small-molecule approaches to reactivate latent wild-type UCH activity of these mutants might be therapeutically viable.Significance: Combined computational and biochemical approaches demonstrate that the BAP1–ASXL2 interaction is direct and high affinity and that many BAP1 mutations act allosterically to inhibit BAP1–ASXL2 binding. Cancer Res; 78(5); 1200–13. ©2017 AACR.

CCR Translation for This Article from Histone Deacetylase Inhibitors Induce Growth Arrest and Differentiation in Uveal Melanoma

PDF file - 572K

The deletion R666-H669 in the ULD domain of BAP1 corresponding to that observed in a uveal melanoma completely abolished ASXL2 binding in vivo. HEK 293 cells were co-transiently transfected with plasmids expressing WT or mutant BAP1 and ASXL2 (261-649aa) proteins. At 48 h post-transfection, the cells were collected and whole cell lysates were used for co-IP with ANTI-Flag M2 Affinity Gel. Western blot analyses were performed with BAP1 (rAb) and Flag (rAb).

Supplemental figure legends

PDF file - 17K

Supplementary Figures 1-8

CARD14-associated papulosquamous eruption (CAPE) is a proposed term that encompasses features ranging from psoriasis to pityriasis rubra pilaris (PRP) in association with CARD14 mutations. The early onset of the disease, prominent facial involvement, family history of an autosomal dominant trait, and poor response to conventional treatment are characteristics of CAPE that distinguish it from classical psoriasis and PRP. We describe the clinical features, family history, and response to therapy in three unrelated children with CAPE and compare these characteristics with those of previously described pediatric patients. Testing for CARD14 mutations in children with early onset of features of psoriasis or pityriasis rubra pilaris and resistance to conventional therapy should be considered.

Dermatologists stand at the gateway of individualization of classification, treatment, and outcomes of acral melanoma patients. The acral melanoma genetic landscape differs in vital ways from that of other cutaneous melanomas. These differences have important implications in understanding pathogenesis, treatment, and prognosis. The selection of molecularly targeted therapy must be adapted for acral melanoma. It is also critical to recognize that tumor development is far more complex than an isolated event, reliably treated by a medication acting on a single target. Tumors exhibit intratumor genetic heterogeneity, metastasis may have different genetic or epigenetic features than primary tumors, and tumor resistance may develop because of the activation of alternative genetic pathways. Microenvironmental, immune, and epigenetic events contribute and sustain tumors in complex ways. Treatment strategies with multiple targets are required to effectively disrupt the tumor ecosystem. This review attempts to translate the current molecular understanding of acral melanoma into digestible concepts relevant to the practice of dermatology. The focus is tumor genetics defining potentially treatable cancer pathways, contextualized within the relevant pathologic and molecular features., Graphical abstract

Funder: Libor Fund grant from the UK Department of Health, by the British Lung Foundation and by the Asmarley Foundation, Funder: UK Medical Research Council, Pleural mesothelioma is an aggressive malignancy with limited effective therapies. In order to identify therapeutic targets, we integrated SNP genotyping, sequencing and transcriptomics from tumours and low-passage patient-derived cells. Previously unrecognised deletions of SUFU locus (10q24.32), observed in 21% of 118 tumours, resulted in disordered expression of transcripts from Hedgehog pathways and the T-cell synapse including VISTA. Co-deletion of Interferon Type I genes and CDKN2A was present in half of tumours and was a predictor of poor survival. We also found previously unrecognised deletions in RB1 in 26% of cases and show sub-micromolar responses to downstream PLK1, CHEK1 and Aurora Kinase inhibitors in primary mesothelioma cells. Defects in Hippo pathways that included RASSF7 amplification and NF2 or LATS1/2 mutations were present in 50% of tumours and were accompanied by micromolar responses to the YAP1 inhibitor Verteporfin. Our results suggest new therapeutic avenues in mesothelioma and indicate targets and biomarkers for immunotherapy.

Purpose: Uveal melanoma is a primary malignancy of the eye with oncogenic mutations in GNAQ, GNA11, or CYSLTR2, and additional mutations in BAP1 (usually associated with LOH of Chr 3), SF3B1, or EIF1AX. There are other characteristic chromosomal alterations, but their significance is not clear. Experimental Design: To investigate genes driving chromosomal alterations, we integrated copy number, transcriptome, and mutation data from three cohorts and followed up key findings. Results: We observed significant enrichment of transcripts on chromosomes 1p, 3, 6, 8, and 16q and identified seven shared focal copy number alterations (FCNAs) on Chr 1p36, 2q37, 3, 6q25, 6q27, and 8q24. Integrated analyses revealed clusters of genes in focal copy number regions whose expression was associated with metastasis and worse overall survival. This included genes from Chr 1p36, 3p21, and 8q24.3. At Chr 6q27, we identified two tumors with homozygous deletion of PHF10/BAF45a and one with a frameshift mutation with concomitant loss of the wild-type allele. Downregulation of PHF10 in uveal melanoma cell lines and tumors altered a number of biological pathways including development and adhesion. These findings provide support for a role for PHF10 as a novel tumor suppressor at Chr 6q27. Conclusions: Integration of copy number, transcriptome, and mutation data revealed novel candidate genes playing a role in uveal melanoma pathogenesis and a potential tumor suppressor role for PHF10.

Circulating tumor DNA (ctDNA) sequencing studies could provide novel insights into the molecular pathology of cancer in sub-Saharan Africa. ctDNA was readily detected in 15 blood samples collected in Ghana at the time of suspected breast cancer. Genomic alterations previously associated with unfavorable prognostic outcomes were observed in the majority of patients. This supports the use of liquid biopsies for diagnosis, surveillance and clinical management of breast cancer in Ghana.

BAP1 is a ubiquitin hydrolase whose deubiquitinase activity is mediated by polycomb group-like protein ASXL2. Cancer-related mutations/deletions of BAP1 lead to loss-of-function either by directly targeting the catalytic (UCH) or ULD domains of BAP1, the latter disrupts binding to ASXL2, an obligate partner for BAP1 enzymatic activity. However, the biochemical and biophysical properties of the domains involved in forming the enzymatically active complex are unknown. Here we investigate the molecular dynamics, kinetics and stoichiometry of these interactions. We demonstrate that the BAP1 and ASXL2 domain/proteins or protein complexes produced in either bacteria or baculovirus are structurally and functionally active. The interaction between BAP1 and ASXL2 is direct, specific, and stable to in vitro biochemical and biophysical manipulations as detected by isothermal titration calorimetry, GST association, and optical biosensor assays. Association of the ASXL2-AB box greatly stimulates BAP1 deubiquitinase activity. A stable ternary complex can be formed comprised of the BAP1-UCH, BAP1-ULD, and ASXL2-AB domains. Binding of the BAP1-ULD domain to the ASXL2-AB box is rapid, with fast association and slow dissociation rates. Stoichiometric analysis revealed that one molecule of the ULD domain directly interacts with one molecule of the AB Box. Real-time kinetics analysis of ULD/AB protein complex to the UCH domain of BAP1, based on SPR, indicated that formation of the ULD/AB complex with the UCH domain is a single-step event with fast association and slow dissociation rates. These structural and dynamic parameters implicate the possibility for future small-molecule approaches to reactivate latent wild-type UCH activity in BAP-mutant malignancies.

BackgroundMalignant pleural mesothelioma (MPM) is an aggressive malignancy with limited effective therapies.MethodsIn order to identify therapeutic targets, we integrated SNP genotyping, sequencing and transcriptomics from tumours and low-passage patient-derived cells.ResultsPreviously unrecognised losses of SUFU locus (10q24.32), observed in 21% of 118 tumours, resulted in disordered expression of transcripts from Hedgehog pathways and the T-cell synapse including VISTA. Co-deletion of Interferon Type I genes and CDKN2A was present in half of tumours and was a predictor of poor survival. We also found previously unrecognised deletions in RB1 in 26% of cases and show sub-micromolar responses to downstream PLK1, CHEK1 and Aurora Kinase inhibitors in primary MPM cells. Defects in Hippo pathways that included RASSF7 amplification and NF2 or LATS1/2 mutations were present in 50% of tumours and were accompanied by micromolar responses to the YAP1 inhibitor Verteporfin.ConclusionsOur results suggest new therapeutic avenues in MPM and provide targets and biomarkers for immunotherapy.

BAP1 is an ubiquitin hydrolase whose deubiquitinase activity is mediated by polycomb group-like protein ASXL2. Cancer-related BAP1 mutations/deletions lead to loss-of-function by targeting the catalytic ubiquitin C-terminal hydrolase (UCH) or UCH37-like domain (ULD) domains of BAP1, and the latter disrupts binding to ASXL2, an obligate partner for BAP1 enzymatic activity. However, the biochemical and biophysical properties of domains involved in forming the enzymatically active complex are unknown. Here, we report the molecular dynamics, kinetics, and stoichiometry of these interactions. We demonstrate that interactions between BAP1 and ASXL2 are direct, specific, and stable to biochemical and biophysical manipulations as detected by isothermal titration calorimetry (ITC), GST association, and optical biosensor assays. Association of the ASXL2-AB box greatly stimulates BAP1 activity. A stable ternary complex is formed, comprised of the BAP1-UCH, BAP1-ULD, and ASXL2-AB domains. Stoichiometric analysis revealed that one molecule of the ULD domain directly interacts with one molecule of the AB box. Real-time kinetic analysis of the ULD/AB protein complex to the BAP1-UCH domain, based on surface plasmon resonance, indicated that formation of the ULD/AB complex with the UCH domain is a single-step event with fast association and slow dissociation rates. In vitro experiments validated in cells that the ASXL-AB box directly regulates BAP1 activity. Implications: Collectively, these data elucidate molecular interactions between specific protein domains regulating BAP1 deubiquitinase activity, thus establishing a foundation for small-molecule approaches to reactivate latent wild-type BAP1 catalytic activity in BAP1-mutant cancers.

Deleterious mutations of the ubiquitin carboxy-terminal hydrolase BAP1 found in cancers are predicted to encode inactive truncated proteins, suggesting that loss of enzyme function is a primary tumorigenic mechanism. However, many tumors exhibit missense mutations or in-frame deletions or insertions, often outside the functionally critical UCH domain in this tumor suppressor protein. Thus, precisely how these mutations inactivate BAP1 is unknown. Here, we show how these mutations affect BAP1 interactions with the Polycomb group-like protein, ASXL2, using combinations of computational modeling technology, molecular biology, and in vitro reconstitution biochemistry. We found that the BAP1–ASXL2 interaction is direct and high affinity, occurring through the ASXH domain of ASXL2, an obligate partner for BAP1 enzymatic activity. The ASXH domain was the minimal domain for binding the BAP1 ULD domain, and mutations on the surfaces of predicted helices of ASXH abolished BAP1 association and stimulation of BAP1 enzymatic activity. The BAP1-UCH, BAP1-ULD, and ASXH domains formed a cooperative stable ternary complex required for deubiquitination. We defined four classes of alterations in BAP1 outside the UCH domain, each failing to productively recruit ASXH to the wild-type BAP1 catalytic site via the ULD, resulting in loss of BAP1 ubiquitin hydrolase activity. Our results indicate that many BAP1 mutations act allosterically to inhibit ASXH binding, thereby leading to loss of enzyme activity. Small-molecule approaches to reactivate latent wild-type UCH activity of these mutants might be therapeutically viable. Significance: Combined computational and biochemical approaches demonstrate that the BAP1–ASXL2 interaction is direct and high affinity and that many BAP1 mutations act allosterically to inhibit BAP1–ASXL2 binding. Cancer Res; 78(5); 1200–13. ©2017 AACR.

Purpose: Analysis of cell free DNA could provide a rapid and non-invasive approach to detect cancer and provide new molecular insights in many African countries where expert pathology is lacking. Hence, we tested whether whole-genome sequencing of cfDNA (WGS-cfDNA) could identify somatic alterations that drive breast cancer. Methods: We conducted a pilot on 15 Ghanaian women (median age 49.5 years) recruited as part of the Ghana Breast Health Study. cfDNA was extracted and subjected to WGS at 30x and 0.1x. ichorCNA software was used to predict copy number alterations and ctDNA fractions. Results: We found extensive amplification and deletion of multiple chromosomal regions including those with oncogenes and tumor suppressor genes associated with breast cancer. Similar copy number alterations for selected breast cancer genes were observed with 0.1x and 30x cfDNA-WGS with increasing concordance between the two instruments as the ctDNA fraction increases. We observed a high frequency (>50%) of copy number gain in 3/5 regions and potential target genes for the amplification (chr8p11-12 [ZNF703] n=8, 53.3%; chr8q24.2 [MYC] n=9, 60%; chr19q12 [CCNE1] n=9, 60%), which were in agreement to previous observations among African-American (AA) ancestry compared to European-American (EA) ancestry in TCGA datasets. Conclusion: Our data provided evidence that ctDNA-based genomic studies are possible and ctDNA analysis could be a tool for future molecular oncology studies in Africa for cancer etiology, surveillance and clinical trials. Citation Format: Samuel Ahuno, Anna-Lisa Doebley, Thomas Ahearn, Joel Yarney, Nicholas Titiloye, Nancy Hamel, Ernest Adjei, Joe-Nat Clegg-Lamptey, Lawrence Edusei, Baffour Awuah, Xiaoyu Song, Verne Vanderpuye, Mustapha Abubakar, Maire Duggan, Daniel Stover, Kofi Nyarko, John Bartlet, Francis Aitpillah, Daniel Ansong, Kevin Gardner, Anne Bowcock, Carlos Caldas, William Foulkes, Seth Wiafe, Wiafe-Addai, Montserrat Garcia-Closas, Alexander Kwarteng, Gavin Ha, Jonine Figueroa, Paz Polak, On Behalf Of Ghana Breast Health Study Team. Studying Ghanian Cancer Genomes Using Cell-free DNA [abstract]. In: Proceedings of the 9th Annual Symposium on Global Cancer Research; Global Cancer Research and Control: Looking Back and Charting a Path Forward; 2021 Mar 10-11. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2021;30(7 Suppl):Abstract nr 81.

The investigation of biological systems involving all organs of the body including the skin is in an era of big data. This requires heavy-duty computational tools, and novel statistical methods. Microarrays have allowed the interrogation of thousands of common genetic markers in thousands of individuals from the same population (termed genome wide association studies or GWAS) to reveal common variation associated with disease or phenotype. These markers are usually single nucleotide polymorphisms (SNPs) that are relatively common in the population. In the case of dermatological diseases such as alopecia areata, vitiligo, psoriasis and atopic dermatitis, common variants have been identified that are associated with disease, and these provide insights into biological pathways and reveal possible novel drug targets. Other skin phenotypes such as acne, color and skin cancers are also being investigated with GWAS. Analyses of such large GWAS datasets require a consideration of a number of statistical issues including the testing of multiple markers, population substructure, and ultimately a requirement for replication. There are also issues regarding the missing heritability of disease that cannot be entirely explained with current GWAS approaches. Next generation sequencing technologies such as exome and genome sequencing of similar patient cohorts will reveal additional variants contributing to disease susceptibility. However, the data generated with these approaches will be orders of magnitude greater than that those generated with arrays, with concomitant challenges in the identification of disease causing variants.

// Sing Y. Moorcraft 1 , Thomas Jones 2 , Brian A. Walker 1 , George Ladas 2 , Eleftheria Kalaitzaki 1 , Lina Yuan 1 , Ruwaida Begum 1 , Zakaria Eltahir 1 , Andrew Wotherspoon 1 , Angeles Montero-Fernandez 2 , Larissa S. Teixeira Mendes 1 , David Gonzalez de Castro 1 , Sanna Hulkki Wilson 1 , Paula Proszek 1 , Ye M. To 1 , Eliza Hawkes 1 , Amitesh Roy 1 , David Cunningham 1 , Sheela Rao 1 , David Watkins 1 , Naureen Starling 1 , Anne M. Bowcock 3 and Ian Chau 1 1 The Royal Marsden NHS Foundation Trust, London and Sutton, United Kingdom 2 The Royal Brompton and Harefield NHS Foundation Trust, London, United Kingdom 3 National Heart and Lung Institute, Imperial College, London, United Kingdom Correspondence to: Ian Chau, email: ian.chau@rmh.nhs.uk Keywords: colorectal cancer, heterogeneity, metastasectomy, pulmonary metastases, RAS Received: October 06, 2016 Accepted: March 29, 2017 Published: April 11, 2017 ABSTRACT This study aimed to molecularly characterise colorectal pulmonary metastases (PM) and investigate whether their molecular profiles were concordant with those of the primary tumour. Clinical data and archival formalin fixed paraffin embedded tissue samples were retrospectively collected from patients who underwent ≥ 1 pulmonary metastasectomies for colorectal cancer between 1997–2012. Primary tumour and metastatic samples were analysed using a targeted capture sequencing panel of 46 cancer-associated genes. The 5-year progression-free and overall survival rates for the 81 patients in this study were 32% (95% CI 22–42%) and 77% (95% CI 66–85%) respectively. Fifty-four patients had samples available from ≥ 1 PM, and sequencing data were successfully obtained from 33 PM from 24 patients. The most frequently mutated genes were APC (71%), KRAS (58%) and TP53 (46%). Seventy-three percent of the 15 patients with matched primary and PM samples and 6 of the 7 patients (86%) with data from ≥ 2 PM had concordant molecular profiles. The concordance for KRAS and NRAS was 100%. At our institutions, patients with resectable colorectal PM had a favourable prognosis. RAS mutations were commonly detected in PM and the molecular profiles of colorectal PM were highly concordant with the primary tumour.

In this issue of Cancer Research, Wang and colleagues identify a large number of regulatory sites within the genomes of non–small cell lung cancers with a global scan for open chromatin (assaying for transposase-accessible chromatin with sequencing). They show that this type of profiling might substitute RNA sequencing in classifying lung cancer samples and in making predictions about prognosis. They also show experimentally that genome editing of some regulatory sites upregulates the expression of GSTM1 and GSTT1, which are required for detoxification of carcinogens and whose low expression levels are associated with lung cancer risk. See related article by Wang et al., p. 4840

Psoriasis (PS) is a common inflammatory and incurable skin disease affecting 2–3% of the human population. Although genome-wide association studies implicate more than 60 loci, the full complement of genetic factors leading to disease is not known. Rare, highly penetrant, gain-of-function, dominantly acting mutations within the human caspase recruitment domain family, member 14 (CARD14) gene lead to the development of PS and psoriatic arthritis (PSA) (a familial p.G117S and de-novo p.E138A alteration). These residues are conserved in mouse and orthologous Knock-In (KI) mutations within Card14 were created. The Card14tm.1.1Sun allele (G117S) resulted in no clinically or histologically evident phenotype of the skin or joints in young adult or old mice. However, mice carrying the Card14tm2.1Sun mutant allele (E138A) were runted and developed thick, white, scaly skin soon after birth, dying within two weeks or less. The skin hyperplasia and inflammation was remarkable similarity to human PS at the clinical, histological, and transcriptomic levels. For example, the skin was markedly acanthotic and exhibited orthokeratotic hyperkeratosis with minimal inflammation and no pustules and transcripts affecting critical pathways of epidermal differentiation and components of the IL17 axis (IL23, IL17A, IL17C, TNF and IL22) were altered. Similar changes were seen in a set of orthologous microRNAs previously associated with PS suggesting conservation across species. Crossing the Card14tm2.1Sun/WT mice to C57BL/6NJ, FVB/NJ, CBA/J, C3H/HeJ, and 129S1/SvImJ generated progeny with epidermal acanthosis and marked orthokeratotic hyperkeratosis regardless of the hybrid strain. Of these hybrid lines, only the FVB;B6N(129S4) mice survived to 250 days of age or older and has led to recombinant inbred lines homozygous for Card14E138A that are fecund and have scaly skin disease. This implicates that modifiers of PS severity exist in mice, as in the familial forms of the disease in humans.

Inherited and de novo mutations in the CARD14 gene promote the development of psoriasis, an inflammatory disease of the skin. Caspase recruitment domain-containing protein 14 (CARD14) is a member of the CARMA protein family that includes the structurally related CARD11 adaptor that mediates NF-κB activation by antigen receptors. We investigated the mechanism by which CARD14 mutation in psoriasis activates NF-κB. In contrast with wild-type CARD14, CARD14E138A and CARD14G117S psoriasis mutants interacted constitutively with BCL10 and MALT1, and triggered BCL10- and MALT1-dependent activation of NF-κB in keratinocytes. These alterations disrupted the inhibitory effect of the CARD14 linker region (LR) on NF-κB activation by facilitating BCL10 binding. Therefore, psoriasis mutations activated CARD14 by a mechanism analogous to oncogenic CARD11 mutations in non-Hodgkin B cell lymphomas. CARD14E138A also stimulated MALT1 paracaspase activity and activated both ERK1/2 and p38α MAP kinases. Inhibition of MALT1 with mepazine reduced CARD14E138A-induced expression of specific psoriasis-associated transcripts in keratinocytes. Our results establish the mechanism whereby gain-of-function CARD14 variants, which induce psoriatic disease in affected individuals, activate pro-inflammatory signalling.

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.

Psoriasis vulgaris is a common, chronic inflammatory skin disease with a complex etiology involving genetic risk factors and environmental triggers. Here we describe the many known genetic predispositions of psoriasis with respect to immune genes and their encoded pathways in psoriasis susceptibility. These genes span an array of functions that involve antigen presentation (HLA-Cw6, ERAP1, ERAP2, MICA), the IL-23 axis (IL12Bp40, IL23Ap19, IL23R, JAK2, TYK2), T-cell development and T-cells polarization (RUNX1, RUNX3, STAT3, TAGAP, IL4, IL13), innate immunity (CARD14, c-REL, TRAF3IP2, DDX58, IFIH1), and negative regulators of immune responses (TNIP1, TNFAIP3, NFKBIA, ZC3H12C, IL36RN, SOCS1). The contribution of some of these gene products to psoriatic disease has also been revealed in recent years through targeting of key immune components, such as the Th17/IL-23 axis which has been highly successful in disease treatment. However, many of the genetic findings involve immune genes with less clear roles in psoriasis pathogenesis. This is particularly the case for those genes involved in innate immunity and negative regulation of immune specific pathways. It is possible that risk alleles of these genes decrease the threshold for the initial activation of the innate immune response. This could then lead to the onslaught of the pathogenic adaptive immune response known to be active in psoriatic skin. However, precisely how these various genes affect immunobiology need to be determined and some are speculated upon in this review. These novel genetic findings also open opportunities to explore novel therapeutic targets and potentially the development of personalized medicine, as well as discover new biology of human skin disease.

Previous studies have identified 41 independent genome-wide significant psoriasis susceptibility loci. After our first psoriasis genome-wide association study, we designed a custom genotyping array to fine-map eight genome-wide significant susceptibility loci known at that time (IL23R, IL13, IL12B, TNIP1, MHC, TNFAIP3, IL23A and RNF114) enabling genotyping of 2269 single-nucleotide polymorphisms (SNPs) in the eight loci for 2699 psoriasis cases and 2107 unaffected controls of European ancestry. We imputed these data using the latest 1000 Genome reference haplotypes, which included both indels and SNPs, to increase the marker density of the eight loci to 49 239 genetic variants. Using stepwise conditional association analysis, we identified nine independent signals distributed across six of the eight loci. In the major histocompatibility complex (MHC) region, we detected three independent signals at rs114255771 (P = 2.94 × 10(-74)), rs6924962 (P = 3.21 × 10(-19)) and rs892666 (P = 1.11 × 10(-10)). Near IL12B we detected two independent signals at rs62377586 (P = 7.42 × 10(-16)) and rs918518 (P = 3.22 × 10(-11)). Only one signal was observed in each of the TNIP1 (rs17728338; P = 4.15 × 10(-13)), IL13 (rs1295685; P = 1.65 × 10(-7)), IL23A (rs61937678; P = 1.82 × 10(-7)) and TNFAIP3 (rs642627; P = 5.90 × 10(-7)) regions. We also imputed variants for eight HLA genes and found that SNP rs114255771 yielded a more significant association than any HLA allele or amino-acid residue. Further analysis revealed that the HLA-C*06-B*57 haplotype tagged by this SNP had a significantly higher odds ratio than other HLA-C*06-bearing haplotypes. The results demonstrate allelic heterogeneity at IL12B and identify a high-risk MHC class I haplotype, consistent with the existence of multiple psoriasis effectors in the MHC.

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; p5.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.