A genome-wide meta-analysis of palmoplantar pustulosis implicates Th2 responses and cigarette smoking in disease pathogenesis ============================================================================================================================= * Ariana Hernandez-Cordero * Laurent Thomas * Alice Smail * Zhao Qin Lim * Jake R Saklatvala * Raymond Chung * Charles J Curtis * Patrick Baum * Sudha Visvanathan * A David Burden * Hywel L Cooper * Giles Dunnill * Christopher EM Griffiths * Nick J Levell * Richard Parslew * Nick J Reynolds * Shyamal Wahie * Richard B Warren * Andrew Wright * The APRICOT and PLUM Study Team * Michael Simpson * Kristian Hveem * Jonathan N Barker * Nick Dand * Mari Loset * Catherine H Smith * Francesca Capon ## ABSTRACT **Background** Palmoplantar pustulosis (PPP) is an inflammatory skin disorder that mostly affects smokers and manifests with painful pustular eruptions on the palms and soles. While the disease can present with concurrent plaque psoriasis, TNF and IL-17/IL-23 inhibitors show limited efficacy. There is therefore a pressing need to uncover PPP disease drivers and therapeutic targets. **Objectives** To identify genetic determinants of PPP and investigate whether cigarette smoking contributes to disease pathogenesis. **Methods** We performed a genome-wide association meta-analysis of three North-European cohorts (n=1,456 PPP cases and 402,050 controls). We then used the scGWAS program to investigate the cell-type specificity of the resulting association signals. We undertook genetic correlation analyses to examine the similarities between PPP and other immune-mediated diseases. Finally, we applied Mendelian randomization to analyze the causal relationship between cigarette smoking and PPP. **Results** We found that PPP is not associated with the main genetic determinants of plaque psoriasis. Conversely, we identified genome-wide significant associations with the *FCGR3A/FCGR3B* and *CCHCR1* loci. We also observed 13 suggestive (*P*<5X10-6) susceptibility regions, including the *IL4/IL13* interval. Accordingly, we demonstrated a significant genetic correlation between PPP and Th2-mediated diseases like atopic dermatitis and ulcerative colitis. We also found that genes mapping to PPP-associated intervals were preferentially expressed in dendritic cells and enriched for T-cell activation pathways. Finally, we undertook a Mendelian randomization analysis, which supported a causal role of cigarette smoking in PPP. **Conclusions** The first genome-wide association study of PPP points to a pathogenic role for deregulated Th2 responses and cigarette smoking. **Clinical implications** The results of the first PPP GWAS support the therapeutic potential of agents that inhibit Th2 responses and target inflammatory pathways activated by cigarette smoking. **Capsule** The genetic analysis of ∼1,400 PPP cases and 400,000 healthy controls points to a causal role of abnormal Th2 responses and cigarette smoking. This supports the therapeutic utility of Th2 inhibition. Key words * palmoplantar pustulosis * genome-wide association study * Th2 * cigarette smoking ## Introduction Palmar plantar pustulosis (PPP) is a severe inflammatory skin disorder presenting with neutrophil-filled pustules on the palms and soles1. The lesions are painful, disabling and stigmatizing, so that PPP has a profound impact on quality of life. This is often compounded by comorbidities such as type 1 diabetes, psoriatic arthritis and Graves’ disease2. At the same time, PPP remains very difficult to manage, with treatment options limited by poor efficacy and/or toxicity3. As PPP is traditionally classified as a pustular variant of psoriasis, clinical trials have tested the efficacy of agents used to good effect in plaque psoriasis (TNF, interleukin [IL]-17 and IL-23 blockers) or generalized pustular psoriasis (IL-1 and IL-36 blockers). These therapeutics, however, showed limited efficacy in PPP is limited4-7. These disappointing results reflect a poor understanding of the pathways underlying disease pathogenesis. The analysis of candidate genes and bulk RNA-sequencing profiles has documented a significant up-regulation of Th17 responses in PPP8-10. This was confirmed in a recent single-cell study carried out by our group, which also demonstrated an abnormal activation of Th2 pathways10. It is, however, unclear whether these alterations are a cause or a secondary manifestation of the disease. PPP disproportionately affects women11, 12 and smokers 11, 12, with evidence of increased disease severity in current/heavy smokers vs former/never smokers13-15. While it has been hypothesized that nicotine alters the expression of IL-36 cytokines14 and anti-inflammatory nicotinic acetylcholine receptors16, a causal role of cigarette smoking has not been established in PPP. Genome-wide association studies (GWAS) and Mendelian randomization approaches have unique potential to reveal pathogenic pathways and disease-promoting exposures, especially when applied to large population biobanks17. This promise, however, has not been realized in PPP, as the disease is poorly annotated in the UK Biobank (<50 cases reported) and its rarity has hindered the ascertainment of adequately powered datasets. Here, we have addressed these issues by combining three independent Northern European cohorts, including a total of 1,456 cases and 402,050 controls. By undertaking a genome-wide association meta-analysis in this extended dataset, we identified two genome-wide significant and 13 suggestive susceptibility loci. We then investigated our association results by genetic correlation analysis and Mendelian randomization, implicating Th2 responses and cigarette smoking in the pathogenesis of PPP. ## Methods ### Study cohorts This work was undertaken according to the principles of the Declaration of Helsinki. Ethical approval was obtained from the Norwegian Data Protection Authority and the Regional Committee for Medical and Health Research Ethics in Central Norway (Reference number 2015/586); London Bridge Research Ethics Committee, UK (ref: 16/LO/2190) and the Coordinating Ethics Committee of the Hospital District of Helsinki and Uusimaa (ref: HUS/990/2017). Written informed consent was obtained from all study participants. Three patient cohorts, originating from the UK, Norway and Finland, were examined. The UK dataset included 288 affected individuals of European descent, who were recruited in specialist dermatology centers as part of the APRICOT clinical trial4 and its sister study PLUM. All were diagnosed by dermatologists, based on clinical examination and consensus criteria18. A total of 7,321 unrelated individuals from the English Longitudinal Study of Aging were analyzed as controls. The Norwegian cases (n=225) were ascertained from the Trøndelag Health Study (HUNT), a population-based study of adult residents from Trøndelag County19. Participant IDs were linked to regional- and national health registries, which enabled the identification of PPP cases through the L40.3 (pustulosis palmaris et plantaris) ICD-10 code. HUNT participants who were not affected by PPP or any other form of psoriasis (n=64,050) were analysed as controls. The Finnish sample was also a population-based cohort. It was derived from the 8th data release of the FinnGen study, a partnership aiming to analyse genome and health data from Finnish biobank participants20. Cases (n=969) were ascertained based on the L40.3 ICD-10 code and the remaining FinnGen participants (n=330,975) were analysed as controls. Key demographic and clinical information for the PPP cases in the three cohorts is summarised in Table E1. ### Effect of disease associated SNPs on gene expression Following genotyping, association testing and meta-analysis (all described in the Supplementary methods), the potential effects of PPP associated SNPs on gene expression were investigated. Summary data-based Mendelian randomization was therefore implemented with SMR v1.3.121. Summary statistics for expression quantitative trait loci (eQTL) identified in non-sun exposed skin were retrieved from the GTEx database (V7 release, [https://gtexportal.org/home/datasets](https://gtexportal.org/home/datasets)) and examined in conjunction with the PPP meta-analysis results. The genome-wide threshold for statistical significance was set at *P*<9.6x10-6 (0.05/5,189 probes). To determine whether gene expression and disease phenotype were influenced by a shared variant (pleiotropy) or multiple variants in linkage disequilibrium (LD) with each other (linkage), a HEIDI test (heterogeneity in dependent instruments) was undertaken. Any SNP generating a statistically significant SMR p-value and a HEIDI p-value >0.05 was deemed to have pleiotropic effects on PPP and gene expression. ### scRNA-seq assisted GWAS analysis (scGWAS) We used scGWAS_r122 to determine if the PPP-associated genes identified in the GWAS are concordantly activated in a particular cell type. We first calculated gene-based p-values by processing the meta-analysis summary statistics with MAGMA23 (Multi-marker Analysis of GenoMic Annotation). The programme aggregates the SNP data available for each gene and implements a multiple regression analysis to test their joint association with the phenotype of interest. Next, we obtained two reference scRNA-seq datasets: i) the peripheral blood mononuclear cells (PBMC) profiles provided by Jia et al (n=1 donor)22 and ii) a subset of the healthy PBMC profiles generated by McCluskey et al (n=3 donors)10, which we re-analysed as described in the Supplementary Methods. We applied the Box-Cox transformation to the original distribution of -log10 (p) from MAGMA and the distribution of log (CPM+1) from the scRNA-seq datasets. Finally, we processed the two normalised datasets with scGWAS. ### Genetic covariance and pathway enrichment analysis We used GNOVA (GeNetic cOVariance Analyzer)24 to investigate pairwise genetic correlations between PPP and traits of interest (plaque psoriasis and four diseases where the involvement of Th2 pathways is well recognised: allergic rhinitis, asthma, atopic dermatitis (AD), and ulcerative colitis (UC)). We retrieved summary statistics from studies shared by the UK biobank (allergic rhinitis; dataset: ukb-b-16499) or published in the literature (all other diseases25-28). We removed low-frequency (minor allele frequency <5%) and palindromic SNPs. We then calculated genetic correlations, using the European 1000 Genomes Phase 329 (v5) dataset as a reference panel for LD estimation. We set the threshold for statistical significance at *P*<0.01 (0.05/5 traits). For pathway enrichment analyses, we processed GWAS summary statistics with MAGMA to obtain gene-based p-values. We ranked the output gene list by p-value and analysed it with GSEA30, using Gene Ontology Biological Process terms as a reference dataset. Network analysis and visualization was then undertaken with Cytoscape31 3.10.1. ### Mendelian randomization To investigate the causality of cigarette smoking in PPP, Mendelian randomization was implemented with TwoSampleMR v0.5.732. The genetic instrument recapitulating the effects of the exposure was derived from a GWAS for smoking initiation undertaken in 3,383,199 individuals33. Clumping of associated variants was performed to only keep one SNP per LD block (LD window=500kb, r2=0.01). If a lead SNP from the exposure GWAS was not present in the outcome dataset (our PPP meta-analysis), a proxy with r2>0.5 was identified with LDproxy v5.6.334, using the European 1000 Genomes Phase 3 dataset as a reference. Finally, palindromic SNPs were removed and the power of the resulting genetic instrument (n=240 markers) was validated by calculating the F-statistic. As F was 38.1, we were able to exclude weak instrument bias (observed when F<10) and implement Mendelian randomization with the inverse variance weighted and weighted median methods. To validate the robustness of our findings, additional tests were performed to examine heterogeneity (inverse variance weighted method), horizontal pleiotropy (Egger intercept) and direction of effect (r2 exposure vs r2 outcome). The effects of individual SNPs were also examined to investigate sensitivity. ## Results ### Genome-wide meta-analysis To identify genetic determinants of PPP we undertook a genome-wide association scan of 288 UK cases and 7,321 unaffected controls. We then used METAL35 to meta-analyze our results with those obtained in two population-based cohorts analyzed by the HUNT study (225 cases and 64,050 controls) and the FinnGen consortium (969 cases and 330,975 controls). The total dataset included 403,506 individuals (1,456 cases and 402,050 controls), genotyped for ∼9.2M SNPs (Table E1). The distribution of P-values did not deviate from that expected in the absence of population stratification (λ=1.05) (Figure E1). The meta-analysis yielded two genome-wide significant loci, spanning *CCHCR1/POU5F1* (*P*=2.9x10-11) and *FCGR3A/3B* (1.6 x10-8) (Figure 1; Figure E2). Of note, *CCHCR1* maps in close proximity of *HLA-C*0602*, the major genetic determinant of plaque psoriasis36. To exclude the possibility that the association we had observed was driven by co-morbid plaque psoriasis (reported in 7-28% of our PPP cases) we examined LD conservation around the *CCHR1* locus. By interrogating population data generated in Finnish and British individuals, we found that the lead SNP in this region was not in LD with rs4406273 (r2<0.15) (Figure E3), a well-known *HLA-C*0602* proxy37. Accordingly, rs4406273 was not associated with PPP in the FinnGen or UK dataset (*P*>0.05 for both; no data was available for HUNT where the SNP had not been typed or imputed). Thus, the *CCHCR1* association is independent of *HLA-C*0602*, which is consistent with LD patterns previously described in other datasets38. ![Figure 1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/01/17/2024.01.17.24301406/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2024/01/17/2024.01.17.24301406/F1) Figure 1: Manhattan plot showing the association signals detected in the GWAS meta-analysis. We next undertook conditional analysis at the *CCHCR1/POU5F1* and *FCGR3A/FCGR3B* loci. We did not detect any secondary association signals, although this may be due to limited statistical power. Beyond *CCHCR1/POU5F1* and *FCGR3A/FCGR3B*, a further 13 regions showed suggestive evidence for association, with P-values <5x10-6 observed in proximity of immune genes such as *IL4, HLA-DRA* and *TNFSF15* (Table 1, Table E2). View this table: [Table 1:](http://medrxiv.org/content/early/2024/01/17/2024.01.17.24301406/T1) Table 1: Summary of genome-wide significant and suggestive association signals ### Analysis of genes and cell types underlying association signals While the lead SNPs for the *CCHCR1/POU5F1* and *FCGR3A/FCGR3B* loci were in LD with coding variants (rs130068 and rs76714703, respectively; r2>0.75 for both) neither of these missense changes was classified as damaging by pathogenicity predictors (CADD39 and MutationTaster40). We therefore examined the overlap between the two association signals and skin expression quantitative trait loci (eQTL) identified by the GTEX consortium41. Using the SMR program, we found evidence of pleiotropy for SNP rs74538320 in the chromosome 1 locus, indicating that the variant affects PPP risk and *FCGR3B* expression (pSMR=5.3x10-6; pHEIDI>0.05). Conversely, the analysis of the *CCHCR1/POU5F1* locus identified a disease-associated SNP (rs1265079) that is unlikely to directly influence gene expression, but may be in LD with a variant regulating *CCHR1* levels (pSMR=4.4x10-6; pHEIDI=4.5x10-5). We ran a parallel analysis of GTEX blood eQTL, but neither of the above SNPs generated a significant SMR p-value. Of note, we previously showed that circulating immune cells are likely to play a role in the immune pathogenesis of PPP10. To further investigate our GWAS results in the light of these observations and determine whether the effects of susceptibility alleles were mediated by specific immune populations, we undertook scRNA-seq assisted GWAS analysis22 (scGWAS). scGWAS derives biological networks that are simultaneously enriched with disease associated genes and genes that are transcriptionally active in a specific cell type. Here, we observed that networks centered around Transmembrane Immune Signaling Adaptor *TYROBP* and *FOXO1* were enriched for genes that are associated with PPP and expressed in dendritic cells (Figure 2). ![Figure 2:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/01/17/2024.01.17.24301406/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2024/01/17/2024.01.17.24301406/F2) Figure 2: PPP-associated genes are enriched in cell networks (gene modules) that are preferentially expressed in denritic cells. The plots show the results obtained using two scRNA-seq reference dasets for PBMC gene expression. These were retrieved from the scGWAS package (left) and the literature (right, see Methods). The plots show the GWAS (x-axis) and scRNA-seq (y-axis) enrichment scores for each module (dots). Real modules are shown in grey with their confidence intervals indicated by a red circle. Virtual modules, which show the null distribution of module scores, are plotted in blue with confidence intervals in black. The vertical and horizontal dashed lines indicate the enrichment significance threshold (z=1.96) and significant real modules are highlighted in red. The z1 value illustrates the signficance of the red module enrichment. Taken together, these post-GWAS analyses indicate that the PPP associated alleles affect immune pathways in skin and circulating dendritic cells. ### Genetic correlation between PPP and Th2-mediated diseases To further investigate the pathophysiological relevance of our findings, we used GNOVA to explore genetic correlations between PPP and other immune traits of interest. Given that PPP can present with concurrent plaque psoriasis1,13, we first investigated this condition, examining publicly available GWAS summary statistics28. However, we did not find any evidence of correlation between PPP and psoriasis. As we and others have reported Th2 activation in PPP skin8,10, we next examined four Th2 mediated disorders (allergy, asthma, AD and UC). We found that PPP shows a positive correlation with AD (r=0.49; *P*=6.1x10-9) and an inverse one with UC (r=-0.20; *P*=1.5x10-4). To validate these observations, we undertook pathway enrichment analyses on the GWAS summary statistics generated in PPP, AD, UC and psoriasis. Using GSEA, we identified the 50 most significantly enriched Gene Ontology Biological Processes for each disease (False Discovery Rate <0.05). As expected for the analysis of four immune-mediated disorders, we observed a substantial overlap between the enriched pathways. Importantly, however, the largest number of shared pathways was observed between PPP and UC (Jaccard similarity index=0.33), and PPP and AD (Jaccard similarity index=0.25). A closer inspection of the results revealed that most of the overlapping Gene Ontology processes relate to T-cell differentiation and activation. Conversely, the overlap between PPP and psoriasis was more limited (Jaccard similarity index=0.22) and mostly underpinned by pathways involved in antigen processing/presentation and T-cell cytotoxicity (Figure 3). Taken together, these observations support a shared genetic basis between PPP and Th2 mediated diseases. ![Figure 3:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/01/17/2024.01.17.24301406/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2024/01/17/2024.01.17.24301406/F3) Figure 3: Overlap between GO Biological Process pathways identified from PPP, psoriasis (PsO), AD and UC GWAS summary statistics. **A)** Venn diagrams illustrating the overlap between the 50 most significantly enriched GO Biological Process pathways identified by GSEA in each condition. **B)** Cytoscape network analysis of the 50 most significantly enriched pathways in PPP. Pathways were clustered using the MCL clusterMaker alogrithm and annotated using AutoAnnotate. Each GO pathway is coloured to indicate overlap with AD, PsO and UC; pathways identified in PPP but showing no overlap with AD, PsO or UC are shown in grey. ### Mendelian randomization indicates a causal role of cigarette smoking in PPP Given the high prevalence of smokers among PPP patients1,11,13, we next used Mendelian Randomization to determine whether cigarette smoking may play a causal role in the disease. We first derived a genetic instrument (n=240 independent SNPs associated with smoking initiation33) recapitulating the exposure to smoking. We then implemented Mendelian randomization using an inverse variance weighted method. This showed a causal influence of the exposure (smoking) on the outcome (PPP) (*P*=2.9x10-4), which was also confirmed with the weighted median method (*P*=1.79x10-3). We tested the robustness of our findings with sensitivity tests measuring heterogeneity between variant causal estimates, horizontal pleiotropy (the phenomenon whereby a SNP may affects the outcome through pathways unrelated to the exposure) and reverse causation. We found no evidence for heterogeneity between variant effects (*P*>0.05), with a leave-one-out test also demonstrating that no SNP had a disproportionate impact on the association between exposure and outcome (Figure E4). We also excluded the presence of horizontal pleiotropy (*P*>0.05) and reverse causation (correct causal direction P-value< 0.001). These results demonstrate that our genetic instrument did not violate key mendelian randomization assumptions, confirming a causal influence of cigarette smoking in PPP. ## Discussion The aim of this study was to identify novel genetic determinants of PPP through the analysis of an extended dataset. We therefore ascertained three independent case resources and obtained a patient sample of unprecedented size (n=1,456), which enabled us to undertake the first GWAS of PPP. Our analysis identified two genome-wide significant association signals. The first maps to the *FCGR* cluster on chromosome 1q23, which spans five genes encoding immunoglobulin gamma receptors. Co-localization with skin eQTLs suggested that the association is driven by SNPs downregulating the expression of *FCGR3B*. Interestingly, variation in *FCGR3B* copy number has been associated with susceptibility to systemic lupus erythematosus and rheumatoid arthritis42,43. Disease risk is specifically driven by reduced *FCGR3B* copy number, which is thought to impair immunoglobulin binding and clearance of immune complexes44. While it is tempting to speculate that a similar mechanism may be at play in PPP, an in-depth dissection of the *FCGR* locus will be required to disentangle the effects of copy number variants and gene eQTLs. The role of the various cell types expressing *FCGR3B* (including T cells, natural killer cells, macrophages and granulocytes) will also require further investigation. Nonetheless, the association between PPP and *FCGR* genes points to an involvement of antibody-mediated pathways, suggesting an autoimmune component for the pathogenesis of PPP. The second association signal maps to chromosome 6p21. While the lead SNP was located within a *POU5F1* intron, our analysis suggested that it may tag an eQTL for the neighbouring *CCHCR1* gene. The latter encodes a coiled coil protein which has been implicated in microtubule assembly, centrosome maturation, mRNA metabolism and cell proliferation45. While a *CCHCR1* coding haplotype has been repeatedly associated with psoriasis susceptibility46, the relevant SNPs are not in LD (r2<0.5) with the PPP risk allele identified here. This is in keeping with the lack of association observed for the *HLA-C*0602* proxy rs4406273. Of note, we also obtained negative findings (*P*>0.05) for the genomic region encoding *IL36RN*, the main genetic determinant of generalised pustular psoriasis47. Thus, PPP is genetically distinct from co-morbid and phenotypically related forms of skin inflammation. The inspection of suggestive loci revealed that SNPs spanning the *IL4/IL3* gene region, which encodes key Th2 cytokines, are associated with PPP (*P*=6.4x10-7 for the lead SNP) and AD27. In this context, our GNOVA analysis identified a significant correlation between PPP and AD, suggesting the existence of shared genetic risk. Gene set enrichment analysis also showed that the genes that are associated with both conditions preferentially map to T-cell differentiation and T cell activation pathways. This is in keeping with the notion of a shared T-cell mediated aetiology. The latter is also consistent with the enrichment of PPP associated genes within expression modules centred around *TYROBP* and *FOXO1*. Indeed, both genes are required for T cell activation by dendritic cells48,49. Thus, our findings point to a causal role of abnormal Th2 responses in PPP. This observation has significant implications, given that Th2 pathways can be targeted by emerging therapeutics, such as IL-4/IL-13 blockers and JAK1 inhibitors50. Interestingly, case reports published in recent months suggest that both class of drugs could have therapeutic efficacy in PPP51,52. To complement our genetic findings, we used Mendelian randomization to show that cigarette smoking has a causal influence on PPP. Interestingly, cigarette smoking is thought to be protective in UC53, a condition that showed a negative correlation with PPP in the GNOVA analysis. This suggests that the shared genetic determinants between PPP and UC regulate inflammatory pathways that can be influenced by cigarette smoking. Our Mendelian randomization results are in keeping with reports that disease severity correlates with cigarette pack years14. In this context, the further observation that PPP severity is worse in current versus former smokers13, suggests that targeting inflammatory pathways influenced by smoking may be an attractive therapeutic strategy. While ours is by some margin, the largest genetic study of PPP, it is still limited by the relatively small size of the patient sample. The different approaches to case ascertainment (recruitment through Dermatology specialist centres for the UK dataset and analysis of ICD-10 codes for the HUNT and Finngen studies) may have also generated heterogeneity within the case sample, affecting statistical power. Thus, the establishment of larger collaborative consortia will be required to identify further genetic determinants of PPP and validate the involvement of pathways that can be targeted by existing therapeutics. ## Supporting information Supplementary Material [[supplements/301406_file02.docx]](pending:yes) ## Data Availability The FinnGen data are publicly available. The UK data will be submitted to the GWAS catalogue. ## Funding This work was supported by NIHR BioResource Centre Maudsley, National Institute for Health Research Maudsley Biomedical Research Centre (BRC) at South London and Maudsley NHS Foundation Trust and Institute of Psychiatry, Psychology and Neuroscience (IoPPN), King’s College London. We gratefully acknowledge capital equipment funding from the Maudsley Charity (Grant Ref. 980) and Guy’s and St Thomas’s Charity (Grant Ref. STR130505). The work was also supported by the NIHR Manchester Biomedical Research Centre (NIHR203308). The Trøndelag Health Study (The HUNT Study) is a collaboration between HUNT Research Centre (Faculty of Medicine and Health Sciences, NTNU, Norwegian University of Science and Technology), Trøndelag County Council, Central Norway Regional Health Authority, and the Norwegian Institute of Public Health. The genotyping in HUNT was financed by the National Institutes of Health; University of Michigan; the Research Council of Norway; the Liaison Committee for Education, Research and Innovation in Central Norway; and the Joint Research Committee between St Olavs hospital and the Faculty of Medicine and Health Sciences, NTNU. The APRICOT trial was funded by the Efficacy and Mechanism Evaluation (EME) Programme, an MRC and NIHR partnership (grant EME 13/50/17 to C.H.S., F.C., J.N.B., C.E.M.G., N.J.R. and A.D.B). This study was supported by Boehringer-Ingelheim and by the Psoriasis Association (grants BSTOP50/5 to CHS and PhD studentship ST3/20 to AHC). ZQL was supported by the Talent Development Fund of KK Women’s and Children’s Hospital, Singapore. NJR is a NIHR Senior Investigator and is also supported by NIHR Newcastle Biomedical Research Centre, NIHR Newcastle In Vitro Diagnostics Co-operative and NIHR Newcastle Patient Safety Research Collaboration. The views expressed are those of the author(s) and not necessarily those of the NHS, Department of Health or King’s College London. ## Conflict of interests F.C. has received grants and consultancy fees from Boehringer-Ingelheim. A.D.B. has received consultancy fees from Boehringer-Ingelheim. C.E.M.G. has received research grants and/or honoraria from Abbvie, Almirall, Anaptysbio, Boehringer-Ingelheim, Bristol Meyers Squibb, Evelo, GSK, Inmagene, Janssen, Lilly, ONO Pharmaceuticals, Novartis, Pfizer and UCB. P.B. and S.V. are Boehringer-Ingelheim employees. S.W. has had non-financial support (sponsorship to attend dermatology conferences) from Janssen, Abbvie, Novartis, Almirall and UCB. J.N.B. declares paid activities as an advisor and speaker for AbbVie, Amgen, Boehringer-Ingelheim, Bristol Myers Squibb, Johnson & Johnson, Lilly, Novartis. ## Acknowledgements We want to acknowledge the participants and investigators of the FinnGen study. In addition to A. David Burden, Hywell L. Cooper, Giles Dunnil, Christopher E. M. Griffiths, Nick J. Levell, Richard Parslew, Nick J Reynolds, Shyamal Wahie, Richard B. Warren, Andrew Wright, Jonathan N. Barker, Catherine H. Smith, and Francesca Capon, who are authors, the following APRICOT and PLUM study team members facilitated patient recruitment and data processing for the APRICOT clinical trial and the PLUM study: Thamir Abraham (Peterborough City Hospital, Peterborough), Muhmad Ali (Worthing Hospital, Worthing, Suzannah August (Poole Hospital NHS Foundation Trust, Poole), David Baudry (Guy’s Hospital, London), Gabrielle Becher (West Glasgow Ambulatory Care Hospital, Glasgow), Anthony Bewley (Whipps Cross Hospital, Leytonstone), Victoria Brown (St Albans City Hospital, St Albans), Victoria Cornelius (Imperial College London, London), Sharizan Ghaffar (Ninewells Hospital and Medical School, Dundee), John Ingram (University Hospital of Wales, Cardiff), Svetlana Kavakleiva (Royal Lancaster Infirmary, Lancaster), Susan Kelly (The Royal Shrewsbury Hospital, Shrewsbury), Mohsen Khorshid (Basildon Hospital, Essex), Helen Lachmann (Royal Free Hospital, London), Effie Ladoyanni (Russells Hall Hospital, Dudley), Helen McAteer (Psoriasis Association, Northampton), John McKenna (Leicester Royal Infirmary, Leicester), Freya Meynell (Guy’s Hospital, London), Prakash Patel (Guy’s Hospital, London), Andrew Pink (Guy’s Hospital, London), Kingsley Powell (Guy’s Hospital, London), Angela Pushparajah (Guy’s Hospital, London), Catriona Sinclair (Broomfield Hospital, Essex) and Rachel Wachsmuth (Royal Devon and Exeter NHS Foundation Trust, Exeter). ## Abbreviations AD : atopic dermatitis; eQTL : expression quantitative trait locus; GWAS : genome-wide association study; LD : linkage disequilibrium; PBMC : peripheral blood mononuclear cells; PPP : palmoplantar pustulosis; UC : ulcerative colitis. * Received January 17, 2024. * Revision received January 17, 2024. * Accepted January 17, 2024. * © 2024, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. 1.1. Griffiths CEM, 2. Barker JN, 3. Bleiker T, 4. Chalmers RJ, 5. Creamer D Burden AD, Kirby B. Psoriasis and related disorders. In: Griffiths CEM, Barker JN, Bleiker T, Chalmers RJ, Creamer D, editors. Rook’s Textbook of Dermatology. Chichester: Wiley-Blackwell; 2016. 2. 2.Kim DH, Lee JY, Cho SI, Jo SJ. Risks of Comorbidities in Patients With Palmoplantar Pustulosis vs Patients With Psoriasis Vulgaris or Pompholyx in Korea. JAMA Dermatol 2022; 158:650–60. 3. 3.Obeid G, Do G, Kirby L, Hughes C, Sbidian E, Le Cleach L. Interventions for chronic palmoplantar pustulosis: abridged Cochrane systematic review and GRADE assessments. Br J Dermatol 2021; 184:1023–32. 4. 4.Cro S, Cornelius VR, Pink AE, Wilson R, Pushpa-Rajah A, Patel P, et al. Anakinra for palmoplantar pustulosis: results from a randomized, double-blind, multicentre, two staged, adaptive placebo controlled trial (APRICOT). Br J Dermatol 2021. 5. 5.Fukasawa T, Yamashita T, Enomoto A, Yoshizaki-Ogawa A, Miyagawa K, Sato S, Yoshizaki A. Optimal treatments and outcome measures of palmoplantar pustulosis: A systematic review and network meta-analysis-based comparison of treatment efficacy. J Eur Acad Dermatol Venereol 2023. 6. 6.Mrowietz U, Bachelez H, Burden AD, Rissler M, Sieder C, Orsenigo R, Chaouche-Teyara K. Secukinumab for moderate-to-severe palmoplantar pustular psoriasis: Results of the 2PRECISE study. J Am Acad Dermatol 2019; 80:1344–52. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jaad.2019.01.066&link_type=DOI) 7. 7.Mrowietz U, Burden AD, Pinter A, Reich K, Schakel K, Baum P, et al. Spesolimab, an Anti-Interleukin-36 Receptor Antibody, in Patients with Palmoplantar Pustulosis: Results of a Phase IIa, Multicenter, Double-Blind, Randomized, Placebo-Controlled Pilot Study. Dermatol Ther (Heidelb) 2021. 8. 8.Baum P, Visvanathan S, Garcet S, Roy J, Schmid R, Bossert S, et al. Pustular psoriasis: Molecular pathways and effects of spesolimab in generalized pustular psoriasis. J Allergy Clin Immunol 2022; 149:1402–12. 9. 9.Bissonnette R, Nigen S, Langley RG, Lynde CW, Tan J, Fuentes-Duculan J, Krueger JG. Increased expression of IL-17A and limited involvement of IL-23 in patients with palmo-plantar (PP) pustular psoriasis or PP pustulosis; results from a randomised controlled trial. J Eur Acad Dermatol Venereol 2014; 28:1298–305. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/jdv.12272&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24112799&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F17%2F2024.01.17.24301406.atom) 10. 10.McCluskey D, Benzian-Olsson N, Mahil SK, Hassi NK, Wohnhaas CT Apricot, et al. Single-cell analysis implicates TH17-to-TH2 cell plasticity in the pathogenesis of palmoplantar pustulosis. J Allergy Clin Immunol 2022. 11. 11.Assan F, Husson B, Hegazy S, Seneschal J, Aubin F, Mahe E, et al. Palmoplantar pustulosis and acrodermatitis continua of Hallopeau: demographic and clinical comparative study in a large multicentre cohort. J Eur Acad Dermatol Venereol 2022; 36:1578–83. 12. 12.Twelves S, Mostafa A, Dand N, Burri E, Farkas K, Wilson R, et al. Clinical and genetic differences between pustular psoriasis subtypes J Allergy Clin Immunol 2019; 143:1021–6. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F17%2F2024.01.17.24301406.atom) 13. 13.Benzian-Olsson N, Dand N, Chaloner C, Bata-Csorgo Z, Borroni R, Burden AD, et al. Association of Clinical and Demographic Factors With the Severity of Palmoplantar Pustulosis. JAMA Dermatol 2020. 14. 14.Kobayashi K, Kamekura R, Kato J, Kamiya S, Kamiya T, Takano K, et al. Cigarette Smoke Underlies the Pathogenesis of Palmoplantar Pustulosis via an IL-17A-Induced Production of IL-36gamma in Tonsillar Epithelial Cells. J Invest Dermatol 2021; 141:1533–41 e4. 15. 15.Sarikaya Solak S, Kara Polat A, Kilic S, Oguz Topal I, Saricaoglu H, Karadag AS, et al. Clinical characteristics, quality of life and risk factors for severity in palmoplantar pustulosis: a cross-sectional, multicentre study of 263 patients. Clin Exp Dermatol 2022; 47:63–71. 16. 16.Hagforsen E, Edvinsson M, Nordlind K, Michaelsson G. Expression of nicotinic receptors in the skin of patients with palmoplantar pustulosis. Br J Dermatol 2002; 146:383–91. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1046/j.1365-2133.2002.04640.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11952537&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F17%2F2024.01.17.24301406.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000174889600004&link_type=ISI) 17. 17.Sanderson E, Glymour MM, Holmes MV, Kang H, Morrison J, Munafo MR, et al. Mendelian randomization. Nat Rev Methods Primers 2022; 2. 18. 18.Navarini AA, Burden AD, Capon F, Mrowietz U, Puig L, Koks S, et al. European Consensus Statement on Phenotypes of Pustular Psoriasis. J Eur Acad Dermatol Venereol 2017:1792–9. 19. 19.Brumpton BM, Graham S, Surakka I, Skogholt AH, Loset M, Fritsche LG, et al. The HUNT study: A population-based cohort for genetic research. Cell Genom 2022; 2:100193. 20. 20.Kurki MI, Karjalainen J, Palta P, Sipila TP, Kristiansson K, Donner KM, et al. FinnGen provides genetic insights from a well-phenotyped isolated population. Nature 2023; 613:508–18. 21. 21.Zhu Z, Zhang F, Hu H, Bakshi A, Robinson MR, Powell JE, et al. Integration of summary data from GWAS and eQTL studies predicts complex trait gene targets. Nat Genet 2016; 48:481–7. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3538&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27019110&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F17%2F2024.01.17.24301406.atom) 22. 22.Jia P, Hu R, Yan F, Dai Y, Zhao Z. scGWAS: landscape of trait-cell type associations by integrating single-cell transcriptomics-wide and genome-wide association studies. Genome Biol 2022; 23:220. 23. 23.de Leeuw CA, Mooij JM, Heskes T, Posthuma D. MAGMA: generalized gene-set analysis of GWAS data. PLoS Comput Biol 2015; 11:e1004219. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pcbi.1004219&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25885710&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F17%2F2024.01.17.24301406.atom) 24. 24.Lu Q, Li B, Ou D, Erlendsdottir M, Powles RL, Jiang T, et al. A Powerful Approach to Estimating Annotation-Stratified Genetic Covariance via GWAS Summary Statistics. Am J Hum Genet 2017; 101:939–64. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2017.11.001&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F17%2F2024.01.17.24301406.atom) 25. 25.de Lange KM, Moutsianas L, Lee JC, Lamb CA, Luo Y, Kennedy NA, et al. Genomewide association study implicates immune activation of multiple integrin genes in inflammatory bowel disease. Nat Genet 2017; 49:256–61. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3760.&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28067908&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F17%2F2024.01.17.24301406.atom) 26. 26.Ferreira MAR, Mathur R, Vonk JM, Szwajda A, Brumpton B, Granell R, et al. Genetic Architectures of Childhood- and Adult-Onset Asthma Are Partly Distinct. Am J Hum Genet 2019; 104:665–84. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2019.02.022&link_type=DOI) 27. 27.Sliz E, Huilaja L, Pasanen A, Laisk T, Reimann E, Magi R, et al. Uniting biobank resources reveals novel genetic pathways modulating susceptibility for atopic dermatitis. J Allergy Clin Immunol 2022; 149:1105–12 e9. 28. 28.Tsoi LC, Spain SL, Knight J, Ellinghaus E, Stuart PE, Capon F, et al. Identification of 15 new psoriasis susceptibility loci highlights the role of innate immunity. Nat Genet 2012; 44:1341–8. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2467&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23143594&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F17%2F2024.01.17.24301406.atom) 29. 29.The 1000 Genomes Project Consortium, Auton A, Brooks LD, Durbin RM, Garrison EP, Kang HM, et al. A global reference for human genetic variation. Nature 2015; 526:68–74. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature15393&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26432245&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F17%2F2024.01.17.24301406.atom) 30. 30.Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A 2005; 102:15545–50. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTAyLzQzLzE1NTQ1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDEvMTcvMjAyNC4wMS4xNy4yNDMwMTQwNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 31. 31.Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 2003; 13:2498–504. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjEwOiIxMy8xMS8yNDk4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDEvMTcvMjAyNC4wMS4xNy4yNDMwMTQwNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 32. 32.Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife 2018; 7. 33. 33.Saunders GRB, Wang X, Chen F, Jang SK, Liu M, Wang C, et al. Genetic diversity fuels gene discovery for tobacco and alcohol use. Nature 2022; 612:720–4. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-022-05477-4&link_type=DOI) 34. 34.Machiela MJ, Chanock SJ. LDlink: a web-based application for exploring populationspecific haplotype structure and linking correlated alleles of possible functional variants. Bioinformatics 2015; 31:3555–7. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btv402&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26139635&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F17%2F2024.01.17.24301406.atom) 35. 35.Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics 2010; 26:2190–1. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btq340&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20616382&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F17%2F2024.01.17.24301406.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000281738900017&link_type=ISI) 36. 36.Tsoi LC, Stuart PE, Tian C, Gudjonsson JE, Das S, Zawistowski M, et al. Large scale meta-analysis characterizes genetic architecture for common psoriasis associated variants. Nat Commun 2017; 8:15382. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncomms15382&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F17%2F2024.01.17.24301406.atom) 37. 37.Stuart PE, Tejasvi T, Shaiq PA, Kullavanijaya P, Qamar R, Raja GK, et al. A single SNP surrogate for genotyping HLA-C*06:02 in diverse populations. J Invest Dermatol 2015; 135:1177–80. 38. 38.Veal CD, Capon F, Allen MH, Heath EK, Evans JC, Jones A, et al. Family-based analysis using a dense single-nucleotide polymorphism-based map defines genetic variation at PSORS1, the major psoriasis-susceptibility locus. Am J Hum Genet 2002; 71:554–64. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/342289&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12148091&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F17%2F2024.01.17.24301406.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000177489600009&link_type=ISI) 39. 39.Rentzsch P, Witten D, Cooper GM, Shendure J, Kircher M. CADD: predicting the deleteriousness of variants throughout the human genome. Nucleic Acids Res 2019; 47:D886–D94. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gky1016&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F17%2F2024.01.17.24301406.atom) 40. 40.Schwarz JM, Rodelsperger C, Schuelke M, Seelow D. MutationTaster evaluates disease-causing potential of sequence alterations. Nat Methods 2010; 7:575–6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nmeth0810-575&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20676075&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F17%2F2024.01.17.24301406.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000280500000014&link_type=ISI) 41. 41.Consortium G. The GTEx Consortium atlas of genetic regulatory effects across human tissues. Science 2020; 369:1318–30. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNjkvNjUwOS8xMzE4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDEvMTcvMjAyNC4wMS4xNy4yNDMwMTQwNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 42. 42.Mueller M, Barros P, Witherden AS, Roberts AL, Zhang Z, Schaschl H, et al. Genomic pathology of SLE-associated copy-number variation at the FCGR2C/FCGR3B/FCGR2B locus. Am J Hum Genet 2013; 92:28–40. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2012.11.013&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23261299&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F17%2F2024.01.17.24301406.atom) 43. 43.Rahbari R, Zuccherato LW, Tischler G, Chihota B, Ozturk H, Saleem S, et al. Understanding the Genomic Structure of Copy-Number Variation of the Low-Affinity Fcgamma Receptor Region Allows Confirmation of the Association of FCGR3B Deletion with Rheumatoid Arthritis. Hum Mutat 2017; 38:390–9. 44. 44.Nossent JC, Becker-Merok A, Rischmueller M, Lester S. Susceptibility for Lupus Nephritis by Low Copy Number of the FCGR3B Gene Is Linked to Increased Levels of Pathogenic Autoantibodies. Autoimmune Dis 2013; 2013:750814. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23864940&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F17%2F2024.01.17.24301406.atom) 45. 45.Tervaniemi MH, Katayama S, Skoog T, Siitonen HA, Vuola J, Nuutila K, et al. Intracellular signalling pathways and cytoskeletal functions converge on the psoriasis candidate gene CCHCR1 expressed at P-bodies and centrosomes. BMC Genomics 2018; 19:432. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12864-018-4810-y&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29866042&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F17%2F2024.01.17.24301406.atom) 46. 46.Asumalahti K, Veal C, Laitinen T, Suomela S, Allen M, Elomaa O, et al. Coding haplotype analysis supports HCR as the putative susceptibility gene for psoriasis at the MHC PSORS1 locus. Hum Mol Genet 2002; 11:589–97. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/11.5.589&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11875053&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F17%2F2024.01.17.24301406.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000174281300012&link_type=ISI) 47. 47.Hussain S, Berki DM, Choon SE, Burden AD, Allen MH, Arostegui JI, et al. IL36RN mutations define a severe auto-inflammatory phenotype of generalized pustular psoriasis. Journal of Allergy and Clinical Immunology 2015; 135:1067–70. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jaci.2014.09.043&link_type=DOI) 48. 48.Dong G, Wang Y, Xiao W, Pacios Pujado S, Xu F, Tian C, et al. FOXO1 regulates dendritic cell activity through ICAM-1 and CCR7. J Immunol 2015; 194:3745–55. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6MTA6IjE5NC84LzM3NDUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wMS8xNy8yMDI0LjAxLjE3LjI0MzAxNDA2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 49. 49.Graham DB, Akilesh HM, Gmyrek GB, Piccio L, Gilfillan S, Sim J, et al. ITAM signaling in dendritic cells controls T helper cell priming by regulating MHC class II recycling. Blood 2010; 116:3208–18. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTI6ImJsb29kam91cm5hbCI7czo1OiJyZXNpZCI7czoxMToiMTE2LzE3LzMyMDgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wMS8xNy8yMDI0LjAxLjE3LjI0MzAxNDA2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 50. 50.Bieber T. Atopic dermatitis: an expanding therapeutic pipeline for a complex disease. Nat Rev Drug Discov 2022; 21:21–40. 51. 51.Gleeson D, Barker J, Capon F, Pink AE, Woolf RT, Smith CH, Mahil SK. Are Janus kinase inhibitors an effective treatment for palmoplantar pustulosis? A critically appraised topic. Br J Dermatol 2023; 188:471–3. 52. 52.Zheng YX, Chen XB, Wang ZY, Cai SQ, Zheng M, Koh LF, et al. Efficacy of dupilumab in palmoplantar pustulosis treatment highlights the role of Th2 inflammation. Allergy 2024. 53. 53.Piovani D, Pansieri C, Kotha SRR, Piazza AC, Comberg CL, Peyrin-Biroulet L, et al. Ethnic Differences in the Smoking-related Risk of Inflammatory Bowel Disease: A Systematic Review and Meta-analysis. J Crohns Colitis 2021; 15:1658–78. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ecco-jcc/jjab047&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F01%2F17%2F2024.01.17.24301406.atom)