The cellular architecture of the antimicrobial response network in human leprosy granulomas.

Granulomas are complex cellular structures composed predominantly of macrophages and lymphocytes that function to contain and kill invading pathogens. Here, we investigated the single-cell phenotypes associated with antimicrobial responses in human leprosy granulomas by applying single-cell and spatial sequencing to leprosy biopsy specimens. We focused on reversal reactions (RRs), a dynamic process whereby some patients with disseminated lepromatous leprosy (L-lep) transition toward self-limiting tuberculoid leprosy (T-lep), mounting effective antimicrobial responses. We identified a set of genes encoding proteins involved in antimicrobial responses that are differentially expressed in RR versus L-lep lesions and regulated by interferon-γ and interleukin-1ß. By integrating the spatial coordinates of the key cell types and antimicrobial gene expression in RR and T-lep lesions, we constructed a map revealing the organized architecture of granulomas depicting compositional and functional layers by which macrophages, T cells, keratinocytes and fibroblasts can each contribute to the antimicrobial response.

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.

Second-Strand Synthesis-Based Massively Parallel scRNA-Seq Reveals Cellular States and Molecular Features of Human Inflammatory Skin Pathologies.

High-throughput single-cell RNA-sequencing (scRNA-seq) methodologies enable characterization of complex biological samples by increasing the number of cells that can be profiled contemporaneously. Nevertheless, these approaches recover less information per cell than low-throughput strategies. To accurately report the expression of key phenotypic features of cells, scRNA-seq platforms are needed that are both high fidelity and high throughput. To address this need, we created Seq-Well S3 ("Second-Strand Synthesis"), a massively parallel scRNA-seq protocol that uses a randomly primed second-strand synthesis to recover complementary DNA (cDNA) molecules that were successfully reverse transcribed but to which a second oligonucleotide handle, necessary for subsequent whole transcriptome amplification, was not appended due to inefficient template switching. Seq-Well S3 increased the efficiency of transcript capture and gene detection compared with that of previous iterations by up to 10- and 5-fold, respectively. We used Seq-Well S3 to chart the transcriptional landscape of five human inflammatory skin diseases, thus providing a resource for the further study of human skin inflammation.

Single-cell profiling of prurigo nodularis demonstrates immune-stromal crosstalk driving profibrotic responses and reversal with nemolizumab.

BACKGROUND: Prurigo nodularis (PN) is a chronic neuroimmune skin disease characterized by bilaterally distributed pruritic hyperkeratotic nodules on extremities and trunk. Neuroimmune dysregulation and chronic scratching are believed to both induce and maintain the characteristic lesions. OBJECTIVES: This study sought to provide a comprehensive view of the molecular pathogenesis of PN at the single-cell level to identify and outline key pathologic processes and the cell types involved. Features that distinguish PN skin from the skin of patients with atopic dermatitis were of particular interest. We further aimed to determine the impact of the IL31RA antagonist, nemolizumab, and its specificity at the single-cell level. METHODS: Single-cell RNA-sequencing of skin from 15 healthy donors and nonlesional and lesional skin from 6 patients each with PN and atopic dermatitis, combined with spatial-sequencing using the 10x Visium platform. Integration with bulk RNA-sequencing data from patients treated with nemolizumab. RESULTS: This study demonstrates that PN is an inflammatory skin disease characterized by both keratinocyte proliferation and activation of profibrotic responses. This study also demonstrates that the COL11A1+ fibroblast subset is a major contributor to fibrosis and is predominantly found in the papillary dermis of PN skin. Activation of fibrotic responses is the main distinguishing feature between PN and atopic dermatitis skin. This study further shows the broad effect of nemolizumab on PN cell types, with a prominent effect driving COL11A1+ fibroblast and keratinocyte responses toward normal. CONCLUSIONS: This study provides a high-resolution characterization of the cell types and cellular processes activated in PN skin, establishing PN as a chronic fibrotic inflammatory skin disease. It further demonstrates the broad effect of nemolizumab on pathological processes in PN skin.

Monogenic etiologies of persistent human papillomavirus infections: A comprehensive systematic review.

PURPOSE: Persistent human papillomavirus infection (PHPVI) causes cutaneous, anogenital, and mucosal warts. Cutaneous warts include common warts, Treeman syndrome, and epidermodysplasia verruciformis, among others. Although more reports of monogenic predisposition to PHPVI have been published with the development of genomic technologies, genetic testing is rarely incorporated into clinical assessments. To encourage broader molecular testing, we compiled a list of the various monogenic etiologies of PHPVI. METHODS: We conducted a systematic literature review to determine the genetic, immunological, and clinical characteristics of patients with PHPVI. RESULTS: The inclusion criteria were met by 261 of 40,687 articles. In 842 patients, 83 PHPVI-associated genes were identified, including 42, 6, and 35 genes with strong, moderate, and weak evidence for causality, respectively. Autosomal recessive inheritance predominated (69%). PHPVI onset age was 10.8 ± 8.6 years, with an interquartile range of 5 to 14 years. GATA2,IL2RG,DOCK8, CXCR4, TMC6, TMC8, and CIB1 are the most frequently reported PHPVI-associated genes with strong causality. Most genes (74 out of 83) belong to a catalog of 485 inborn errors of immunity-related genes, and 40 genes (54%) are represented in the nonsyndromic and syndromic combined immunodeficiency categories. CONCLUSION: PHPVI has at least 83 monogenic etiologies and a genetic diagnosis is essential for effective management.

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.

Genome-Wide Association Study in Acute Tubulointerstitial Nephritis.

SIGNIFICANCE STATEMENT: Polymorphisms of HLA genes may confer susceptibility to acute tubulointerstitial nephritis (ATIN), but small sample sizes and candidate gene design have hindered their investigation. The first genome-wide association study of ATIN identified two significant loci, risk haplotype DRB1*14-DQA1*0101-DQB1*0503 (DR14 serotype) and protective haplotype DRB1*1501-DQA1*0102-DQB1*0602 (DR15 serotype), with amino acid position 60 in the peptide-binding groove P10 of HLA-DR ß 1 key. Risk alleles were shared among different causes of ATIN and HLA genotypes associated with kidney injury and immune therapy response. HLA alleles showed the strongest association. The findings suggest that a genetically conferred risk of immune dysregulation is part of the pathogenesis of ATIN. BACKGROUND: Acute tubulointerstitial nephritis (ATIN) is a rare immune-related disease, accounting for approximately 10% of patients with unexplained AKI. Previous elucidation of the relationship between genetic factors that contribute to its pathogenesis was hampered because of small sample sizes and candidate gene design. METHODS: We undertook the first two-stage genome-wide association study and meta-analysis involving 544 kidney biopsy-defined patients with ATIN and 2346 controls of Chinese ancestry. We conducted statistical fine-mapping analysis, provided functional annotations of significant variants, estimated single nucleotide polymorphism (SNP)-based heritability, and checked genotype and subphenotype correlations. RESULTS: Two genome-wide significant loci, rs35087390 of HLA-DQA1 ( P =3.01×10 -39 ) on 6p21.32 and rs2417771 of PLEKHA5 on 12p12.3 ( P =2.14×10 -8 ), emerged from the analysis. HLA imputation using two reference panels suggested that HLA-DRB1*14 mainly drives the HLA risk association . HLA-DRB1 residue 60 belonging to pocket P10 was the key amino acid position. The SNP-based heritability estimates with and without the HLA locus were 20.43% and 10.35%, respectively. Different clinical subphenotypes (drug-related or tubulointerstitial nephritis and uveitis syndrome) seemed to share the same risk alleles. However, the HLA risk genotype was associated with disease severity and response rate to immunosuppressive therapy. CONCLUSIONS: We identified two candidate genome regions associated with susceptibility to ATIN. The findings suggest that a genetically conferred risk of immune dysregulation is involved in the pathogenesis of ATIN.

Novel insights into atopic dermatitis.

Recent research into the pathophysiology and treatment of atopic dermatitis (AD) has shown notable progress. An increasing number of aspects of the immune system are being implicated in AD, including the epithelial barrier, TH2 cytokines, and mast cells. Major advances in therapeutics were made in biologic cytokine and receptor antagonists and among Janus kinase inhibitors. We focus on these areas and address new insights into AD epidemiology, biomarkers, endotypes, prevention, and comorbidities. Going forward, we expect future mechanistic insights and therapeutic advances to broaden physicians' ability to diagnose and manage AD patients, and perhaps to find a cure for this chronic condition.

Serial KinderMiner (SKiM) discovers and annotates biomedical knowledge using co-occurrence and transformer models.

BACKGROUND: The PubMed archive contains more than 34 million articles; consequently, it is becoming increasingly difficult for a biomedical researcher to keep up-to-date with different knowledge domains. Computationally efficient and interpretable tools are needed to help researchers find and understand associations between biomedical concepts. The goal of literature-based discovery (LBD) is to connect concepts in isolated literature domains that would normally go undiscovered. This usually takes the form of an A-B-C relationship, where A and C terms are linked through a B term intermediate. Here we describe Serial KinderMiner (SKiM), an LBD algorithm for finding statistically significant links between an A term and one or more C terms through some B term intermediate(s). The development of SKiM is motivated by the observation that there are only a few LBD tools that provide a functional web interface, and that the available tools are limited in one or more of the following ways: (1) they identify a relationship but not the type of relationship, (2) they do not allow the user to provide their own lists of B or C terms, hindering flexibility, (3) they do not allow for querying thousands of C terms (which is crucial if, for instance, the user wants to query connections between a disease and the thousands of available drugs), or (4) they are specific for a particular biomedical domain (such as cancer). We provide an open-source tool and web interface that improves on all of these issues. RESULTS: We demonstrate SKiM's ability to discover useful A-B-C linkages in three control experiments: classic LBD discoveries, drug repurposing, and finding associations related to cancer. Furthermore, we supplement SKiM with a knowledge graph built with transformer machine-learning models to aid in interpreting the relationships between terms found by SKiM. Finally, we provide a simple and intuitive open-source web interface ( https://skim.morgridge.org ) with comprehensive lists of drugs, diseases, phenotypes, and symptoms so that anyone can easily perform SKiM searches. CONCLUSIONS: SKiM is a simple algorithm that can perform LBD searches to discover relationships between arbitrary user-defined concepts. SKiM is generalized for any domain, can perform searches with many thousands of C term concepts, and moves beyond the simple identification of an existence of a relationship; many relationships are given relationship type labels from our knowledge graph.

MTMR3 risk alleles enhance Toll Like Receptor 9-induced IgA immunity in IgA nephropathy.

Multiple genome-wide association studies (GWASs) have reproducibly identified the MTMR3/HORMAD2/LIF/OSM locus to be associated with IgA nephropathy (IgAN). However, the causal variant(s), implicated gene(s), and altered mechanisms remain poorly understood. Here, we performed fine-mapping analyses based on GWAS datasets encompassing 2762 IgAN cases and 5803 control individuals, and identified rs4823074 as the candidate causal variant that intersects the MTMR3 promoter in B-lymphoblastoid cells. Mendelian randomization studies suggested the risk allele may modulate disease susceptibility by affecting serum IgA levels through increased MTMR3 expression. Consistently, elevated MTMR3 expression in peripheral blood mononuclear cells was observed in patients with IgAN. Further mechanistic studies in vitro demonstrated that MTMR3 increased IgA production dependent upon its phosphatidylinositol 3-phosphate binding domain. Moreover, our study provided the in vivo functional evidence that Mtmr3-/- mice exhibited defective Toll Like Receptor 9-induced IgA production, glomerular IgA deposition, as well as mesangial cell proliferation. RNA-seq and pathway analyses showed that MTMR3 deficiency resulted in an impaired intestinal immune network for IgA production. Thus, our results support the role of MTMR3 in IgAN pathogenesis by enhancing Toll Like Receptor 9-induced IgA immunity.

The Histone Methyltransferase SETDB2 Modulates Tissue Inhibitors of Metalloproteinase-Matrix Metalloproteinase Activity During Abdominal Aortic Aneurysm Development.

OBJECTIVE: To determine macrophage-specific alterations in epigenetic enzyme function contributing to the development of abdominal aortic aneurysms (AAAs). BACKGROUND: AAA is a life-threatening disease, characterized by pathologic vascular remodeling driven by an imbalance of matrix metalloproteinases and tissue inhibitors of metalloproteinases (TIMPs). Identifying mechanisms regulating macrophage-mediated extracellular matrix degradation is of critical importance to developing novel therapies. METHODS: The role of SET Domain Bifurcated Histone Lysine Methyltransferase 2 (SETDB2) in AAA formation was examined in human aortic tissue samples by single-cell RNA sequencing and in a myeloid-specific SETDB2 deficient murine model induced by challenging mice with a combination of a high-fat diet and angiotensin II. RESULTS: Single-cell RNA sequencing of human AAA tissues identified SETDB2 was upregulated in aortic monocyte/macrophages and murine AAA models compared with controls. Mechanistically, interferon-ß regulates SETDB2 expression through Janus kinase/signal transducer and activator of transcription signaling, which trimethylates histone 3 lysine 9 on the TIMP1-3 gene promoters thereby suppressing TIMP1-3 transcription and leading to unregulated matrix metalloproteinase activity. Macrophage-specific knockout of SETDB2 ( Setdb2f/fLyz2Cre+ ) protected mice from AAA formation with suppression of vascular inflammation, macrophage infiltration, and elastin fragmentation. Genetic depletion of SETDB2 prevented AAA development due to the removal of the repressive histone 3 lysine 9 trimethylation mark on the TIMP1-3 gene promoter resulting in increased TIMP expression, decreased protease activity, and preserved aortic architecture. Lastly, inhibition of the Janus kinase/signal transducer and activator of the transcription pathway with an FDA-approved inhibitor, Tofacitinib, limited SETDB2 expression in aortic macrophages. CONCLUSIONS: These findings identify SETDB2 as a critical regulator of macrophage-mediated protease activity in AAAs and identify SETDB2 as a mechanistic target for the management of AAAs.

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.

DDX58(RIG-I)-related disease is associated with tissue-specific interferon pathway activation.

BACKGROUND: Singleton-Merten syndrome (SGMRT) is a rare immunogenetic disorder that variably features juvenile open-angle glaucoma (JOAG), psoriasiform skin rash, aortic calcifications and skeletal and dental dysplasia. Few families have been described and the genotypic and phenotypic spectrum is poorly defined, with variants in DDX58 (DExD/H-box helicase 58) being one of two identified causes, classified as SGMRT2. METHODS: Families underwent deep systemic phenotyping and exome sequencing. Functional characterisation with in vitro luciferase assays and in vivo interferon signature using bulk and single cell RNA sequencing was performed. RESULTS: We have identified a novel DDX58 variant c.1529A>T p.(Glu510Val) that segregates with disease in two families with SGMRT2. Patients in these families have widely variable phenotypic features and different ethnic background, with some being severely affected by systemic features and others solely with glaucoma. JOAG was present in all individuals affected with the syndrome. Furthermore, detailed evaluation of skin rash in one patient revealed sparse inflammatory infiltrates in a unique distribution. Functional analysis showed that the DDX58 variant is a dominant gain-of-function activator of interferon pathways in the absence of exogenous RNA ligands. Single cell RNA sequencing of patient lesional skin revealed a cellular activation of interferon-stimulated gene expression in keratinocytes and fibroblasts but not in neighbouring healthy skin. CONCLUSIONS: These results expand the genotypic spectrum of DDX58-associated disease, provide the first detailed description of ocular and dermatological phenotypes, expand our understanding of the molecular pathogenesis of this condition and provide a platform for testing response to therapy.

Relationship between genetically proxied vitamin D and psoriasis risk: a Mendelian randomization study.

BACKGROUND: Observational research suggests that vitamin D levels affect psoriasis. However, observational studies are prone to potential confounding or reverse causation, which complicates interpreting the data and drawing causal conclusions. AIM: To apply Mendelian randomization (MR) methods to comprehensively assess a potential association between vitamin D and psoriasis. METHODS: Genetic variants strongly associated with 25-hydroxyvitamin D (25OHD) in genome-wide association study (GWAS) data from 417 580 and 79 366 individuals from two independent studies served as instrumental variables (used as the discovery and replication datasets, respectively). As the outcome variable, we used GWAS data of psoriasis (13 229 people in the case group, 21 543 in the control group). We used (i) biologically validated genetic instruments, and (ii) polygenic genetic instruments to assess the relationship between genetically proxied vitamin D and psoriasis. We carried out inverse-variance weighted (IVW) MR analyses for the primary analysis. In sensitivity analyses, we used robust MR approaches. RESULTS: MR analyses of both the discovery and replication datasets did not show an effect of 25OHD on psoriasis. Neither the IVW MR analysis of the biologically validated instruments [discovery dataset: odds ratio (OR) 0.99; 95% confidence interval (CI) 0.88-1.12, P = 0.873; replication dataset: OR 0.98, 95% CI 0.66-1.46, P = 0.930] nor that of the polygenic genetic instruments (discovery dataset: OR 1.00, 95% CI 0.81-1.22, P = 0.973; replication dataset: OR 0.94, 95% CI 0.64-1.38, P = 0.737) revealed an impact of 25OHD on psoriasis. CONCLUSION: The present MR study did not support the hypothesis that vitamin D levels, measured by 25OHD, affect psoriasis. This study was conducted on Europeans, so the conclusions may not be applicable to all ethnicities.

Circadian control of interferon-sensitive gene expression in murine skin.

The circadian clock coordinates a variety of immune responses with signals from the external environment to promote survival. We investigated the potential reciprocal relationship between the circadian clock and skin inflammation. We treated mice topically with the Toll-like receptor 7 (TLR7) agonist imiquimod (IMQ) to activate IFN-sensitive gene (ISG) pathways and induce psoriasiform inflammation. IMQ transiently altered core clock gene expression, an effect mirrored in human patient psoriatic lesions. In mouse skin 1 d after IMQ treatment, ISGs, including the key ISG transcription factor IFN regulatory factor 7 (Irf7), were more highly induced after treatment during the day than the night. Nuclear localization of phosphorylated-IRF7 was most prominently time-of-day dependent in epidermal leukocytes, suggesting that these cell types play an important role in the diurnal ISG response to IMQ. Mice lacking Bmal1 systemically had exacerbated and arrhythmic ISG/Irf7 expression after IMQ. Furthermore, daytime-restricted feeding, which affects the phase of the skin circadian clock, reverses the diurnal rhythm of IMQ-induced ISG expression in the skin. These results suggest a role for the circadian clock, driven by BMAL1, as a negative regulator of the ISG response, and highlight the finding that feeding time can modulate the skin immune response. Since the IFN response is essential for the antiviral and antitumor effects of TLR activation, these findings are consistent with the time-of-day-dependent variability in the ability to fight microbial pathogens and tumor initiation and offer support for the use of chronotherapy for their treatment.

Netherton syndrome subtypes share IL-17/IL-36 signature with distinct IFN-α and allergic responses.

BACKGROUND: Netherton syndrome (NS) is a rare recessive skin disorder caused by loss-of-function mutations in SPINK5 encoding the protease inhibitor LEKTI (lymphoepithelial Kazal-type-related inhibitor). NS patients experience severe skin barrier defects, display inflammatory skin lesions, and have superficial scaling with atopic manifestations. They present with typical ichthyosis linearis circumflexa (NS-ILC) or scaly erythroderma (NS-SE). OBJECTIVE: We used a combination of several molecular profiling methods to comprehensively characterize the skin, immune cells, and allergic phenotypes of NS-ILC and NS-SE patients. METHODS: We studied a cohort of 13 patients comprising 9 NS-ILC and 4 NS-SE. RESULTS: Integrated multiomics revealed abnormal epidermal proliferation and differentiation and IL-17/IL-36 signatures in lesion skin and in blood in both NS endotypes. Although the molecular profiles of NS-ILC and NS-SE lesion skin were very similar, nonlesion skin of each disease subtype displayed distinctive molecular features. Nonlesion and lesion NS-SE epidermis showed activation of the type I IFN signaling pathway, while lesion NS-ILC skin differed from nonlesion NS-ILC skin by increased complement activation and neutrophil infiltration. Serum cytokine profiling and immunophenotyping of circulating lymphocytes showed a TH2-driven allergic response in NS-ILC, whereas NS-SE patients displayed mainly a TH9 axis with increased CCL22/MDC and CCL17/TARC serum levels. CONCLUSIONS: This study confirms IL-17/IL-36 as the predominant signaling axes in both NS endotypes and unveils molecular features distinguishing NS-ILC and NS-SE. These results identify new therapeutic targets and could pave the way for precision medicine of NS.

Cytokine responses in nonlesional psoriatic skin as clinical predictor to anti-TNF agents.

BACKGROUND: A major issue with the current management of psoriasis is our inability to predict treatment response. OBJECTIVE: Our aim was to evaluate the ability to use baseline molecular expression profiling to assess treatment outcome for patients with psoriasis. METHODS: We conducted a longitudinal study of 46 patients with chronic plaque psoriasis treated with anti-TNF agent etanercept, and molecular profiles were assessed in more than 200 RNA-seq samples. RESULTS: We demonstrated correlation between clinical response and molecular changes during the course of the treatment, particularly for genes responding to IL-17A/TNF in keratinocytes. Intriguingly, baseline gene expressions in nonlesional, but not lesional, skin were the best marker of treatment response at week 12. We identified USP18, a known regulator of IFN responses, as positively correlated with Psoriasis Area and Severity Index (PASI) improvement (P = 9.8 × 10-4) and demonstrate its role in regulating IFN/TNF responses in keratinocytes. Consistently, cytokine gene signatures enriched in baseline nonlesional skin expression profiles had strong correlations with PASI improvement. Using this information, we developed a statistical model for predicting PASI75 (ie, 75% of PASI improvement) at week 12, achieving area under the receiver-operating characteristic curve value of 0.75 and up to 80% accurate PASI75 prediction among the top predicted responders. CONCLUSIONS: Our results illustrate feasibility of assessing drug response in psoriasis using nonlesional skin and implicate involvement of IFN regulators in anti-TNF responses.

Transcriptomic characterization of prurigo nodularis and the therapeutic response to nemolizumab.

BACKGROUND: Prurigo nodularis (PN) is a debilitating, difficult-to-treat, intensely pruritic, chronic inflammatory skin disease characterized by hyperkeratotic skin nodules. The pathogenesis of PN is not well understood but is believed to involve cross talk between sensory nerve fibers, immune cells, and the epidermis. It is centered around the neuroimmune cytokine IL-31, driving an intractable itch-scratch cycle. OBJECTIVE: We sought to provide a comprehensive view of the transcriptomic changes in PN skin and characterize the mechanism of action of the anti-IL-31 receptor inhibitor nemolizumab. METHOD: RNA sequencing of biopsy samples obtained from a cohort of patients treated with the anti-IL-31 receptor inhibitor nemolizumab and taken at baseline and week 12. Generation and integration of patient data with RNA-Seq data generated from reconstructed human epidermis stimulated with IL-31 and other proinflammatory cytokines. RESULTS: Our results demonstrate that nemolizumab effectively decreases IL-31 responses in PN skin, leading to effective suppression of downstream inflammatory responses including TH2/IL-13 and TH17/IL-17 responses. This is accompanied by decreased keratinocyte proliferation and normalization of epidermal differentiation and function. Furthermore, our results demonstrate how transcriptomic changes associated with nemolizumab treatment correlate with improvement in lesions, pruritus, stabilization of extracellular matrix remodeling, and processes associated with cutaneous nerve function. CONCLUSION: These data demonstrate a broad response to IL-31 receptor inhibition with nemolizumab and confirm the critical upstream role of IL-31 in PN pathogenesis.

TIGAR: An Improved Bayesian Tool for Transcriptomic Data Imputation Enhances Gene Mapping of Complex Traits.

The transcriptome-wide association studies (TWASs) that test for association between the study trait and the imputed gene expression levels from cis-acting expression quantitative trait loci (cis-eQTL) genotypes have successfully enhanced the discovery of genetic risk loci for complex traits. By using the gene expression imputation models fitted from reference datasets that have both genetic and transcriptomic data, TWASs facilitate gene-based tests with GWAS data while accounting for the reference transcriptomic data. The existing TWAS tools like PrediXcan and FUSION use parametric imputation models that have limitations for modeling the complex genetic architecture of transcriptomic data. Therefore, to improve on this, we employ a nonparametric Bayesian method that was originally proposed for genetic prediction of complex traits, which assumes a data-driven nonparametric prior for cis-eQTL effect sizes. The nonparametric Bayesian method is flexible and general because it includes both of the parametric imputation models used by PrediXcan and FUSION as special cases. Our simulation studies showed that the nonparametric Bayesian model improved both imputation R2 for transcriptomic data and the TWAS power over PrediXcan when ≥1% cis-SNPs co-regulate gene expression and gene expression heritability ≤0.2. In real applications, the nonparametric Bayesian method fitted transcriptomic imputation models for 57.8% more genes over PrediXcan, thus improving the power of follow-up TWASs. We implement both parametric PrediXcan and nonparametric Bayesian methods in a convenient software tool "TIGAR" (Transcriptome-Integrated Genetic Association Resource), which imputes transcriptomic data and performs subsequent TWASs using individual-level or summary-level GWAS data.

Single-cell Transcriptomics Reveals Dynamic Role of Smooth Muscle Cells and Enrichment of Immune Cell Subsets in Human Abdominal Aortic Aneurysms.

OBJECTIVE: To determine cell-specific gene expression profiles that contribute to development of abdominal aortic aneurysms (AAAs). BACKGROUND: AAAs represent the most common pathological aortic dilation leading to the fatal consequence of aortic rupture. Both immune and structural cells contribute to aortic degeneration, however, gene specific alterations in these cellular subsets are poorly understood. METHODS: We performed single-cell RNA sequencing (scRNA-seq) analysis of AAAs and control tissues. AAA-related changes were examined by comparing gene expression profiles as well as detailed receptor-ligand interactions. An integrative analysis of scRNA-seq data with large genome-wide association study data was conducted to identify genes critical for AAA development. RESULTS: Using scRNA-seq we provide the first comprehensive characterization of the cellular landscape in human AAA tissues. Unbiased clustering analysis of transcriptional profiles identified seventeen clusters representing 8 cell lineages. For immune cells, clustering analysis identified 4 T-cell and 5 monocyte/macrophage subpopulations, with distinct transcriptional profiles in AAAs compared to controls. Gene enrichment analysis on immune subsets identified multiple pathways only expressed in AAA tissue, including those involved in mitochondrial dysfunction, proliferation, and cytokine secretion. Moreover, receptor-ligand analysis defined robust interactions between vascular smooth muscle cells and myeloid populations in AAA tissues. Lastly, integrated analysis of scRNA-seq data with genome-wide association study studies determined that vascular smooth muscle cell expression of SORT1 is critical for maintaining normal aortic wall function. CONCLUSIONS: Here we provide the first comprehensive evaluation of single-cell composition of the abdominal aortic wall and reveal how the gene expression landscape is altered in human AAAs.