Genome-wide meta-analysis and fine-mapping prioritize potential causal variants and genes related to leprosy.

To date, genome-wide association studies (GWASs) have discovered 35 susceptible loci of leprosy; however, the cumulative effects of these loci can only partially explain the overall risk of leprosy, and the causal variants and genes within these loci remain unknown. Here, we conducted out new GWASs in two independent cohorts of 5007 cases and 4579 controls and then a meta-analysis in these newly generated and multiple previously published (2277 cases and 3159 controls) datasets were performed. Three novel and 15 previously reported risk loci were identified from these datasets, increasing the known leprosy risk loci of explained genetic heritability from 23.0 to 38.5%. A comprehensive fine-mapping analysis was conducted, and 19 causal variants and 14 causal genes were identified. Specifically, manual checking of epigenomic information from the Epimap database revealed that the causal variants were mainly located within the immune-relevant or immune-specific regulatory elements. Furthermore, by using gene-set, tissue, and cell-type enrichment analyses, we highlighted the key roles of immune-related tissues and cells and implicated the PD-1 signaling pathways in the pathogenetic mechanism of leprosy. Collectively, our study identified candidate causal variants and elucidated the potential regulatory and coding mechanisms for genes associated with leprosy.

Single-cell sequencing analysis reveals development and differentiation trajectory of Schwann cells manipulated by M. leprae.

BACKGROUND: M. leprae preferentially infects Schwann cells (SCs) in the peripheral nerves leading to nerve damage and irreversible disability. Knowledge of how M. leprae infects and interacts with host SCs is essential for understanding mechanisms of nerve damage and revealing potential new therapeutic strategies. METHODOLOGY/PRINCIPAL FINDINGS: We performed a time-course single-cell sequencing analysis of SCs infected with M. leprae at different time points, further analyzed the heterogeneity of SCs, subpopulations associated with M. leprae infection, developmental trajectory of SCs and validated by Western blot or flow cytometry. Different subpopulations of SCs exhibiting distinct genetic features and functional enrichments were present. We observed two subpopulations associated with M. leprae infection, a stem cell-like cell subpopulation increased significantly at 24 h but declined by 72 h after M. leprae infection, and an adipocyte-like cell subpopulation, emerged at 72 h post-infection. The results were validated and confirmed that a stem cell-like cell subpopulation was in the early stage of differentiation and could differentiate into an adipocyte-like cell subpopulation. CONCLUSIONS/SIGNIFICANCE: Our results present a systematic time-course analysis of SC heterogeneity after infection by M. leprae at single-cell resolution, provide valuable information to understand the critical biological processes underlying reprogramming and lipid metabolism during M. leprae infection of SCs, and increase understanding of the disease-causing mechanisms at play in leprosy patients as well as revealing potential new therapeutic strategies.

IL-23/IL-23R Promote Macrophage Pyroptosis and T Helper 1/T Helper 17 Cell Differentiation in Mycobacterial Infection.

Pathogen-induced epigenetic modifications can reshape anti-infection immune processes and control the magnitude of host responses. DNA methylation profiling has identified crucial aberrant methylation changes associated with diseases, thus providing biological insights into the roles of epigenetic factors in mycobacterial infection. In this study, we performed a genome-wide methylation analysis of skin biopsies from patients with leprosy and healthy controls. T helper 17 differentiation pathway was found to be significantly associated with leprosy through functional enrichment analysis. As a key gene in this pathway, IL-23R was found to be critical to mycobacterial immunity in leprosy, according to integrated analysis with DNA methylation, RNA sequencing, and GWASs. Functional analysis revealed that IL-23/IL-23R-enhanced bacterial clearance by activating caspase-1/GSDMD-mediated pyroptosis in a manner dependent on NLRP3 through signal transducer and activator of transcription 3 signaling in macrophages. Moreover, IL23/IL-23R promoted T helper 1 and T helper 17 cell differentiation and proinflammatory cytokine secretion, thereby increasing host bactericidal activity. IL-23R knockout attenuated the effects and increased susceptibility to mycobacterial infection mentioned earlier. These findings illustrate the biological functions of IL-23/IL-23R in modulating intracellular bacterial clearance in macrophages and further support their regulatory effects in T helper cell differentiation. Our study highlights that IL-23/IL-23R might serve as potential targets for the prevention and treatment of leprosy and other mycobacterial infections.

Identification of shared loci associated with both Crohn's disease and leprosy in East Asians.

Genome-wide association studies (GWAS) of Crohn's disease (CD) in European and leprosy in Chinese population have shown that CD and leprosy share genetic risk loci. As these shared loci were identified through cross-comparisons across different ethnic populations, we hypothesized that meta-analysis of GWAS on CD and leprosy in East Asian populations would increase power to identify additional shared loci. We performed a cross-disease meta-analysis of GWAS data from CD (1621 cases and 4419 controls) and leprosy (2901 cases 3801 controls) followed by replication in additional datasets comprising 738 CD cases and 488 controls and 842 leprosy cases and 925 controls. We identified one novel locus at 7p22.3, rs77992257 in intron 2 of ADAP1, shared between CD and leprosy with genome-wide significance (P = 3.80 × 10-11) and confirmed 10 previously established loci in both diseases: IL23R, IL18RAP, IL12B, RIPK2, TNFSF15, ZNF365-EGR2, CCDC88B, LACC1, IL27, NOD2. Phenotype variance explained by the polygenic risk scores derived from Chinese leprosy data explained up to 5.28% of variance of Korean CD, supporting similar genetic structures between the two diseases. Although CD and leprosy shared a substantial number of genetic susceptibility loci in East Asians, the majority of shared susceptibility loci showed allelic effects in the opposite direction. Investigation of the genetic correlation using cross-trait linkage disequilibrium score regression also showed a negative genetic correlation between CD and leprosy (rg [SE] = -0.40[0.13], P = 2.6 × 10-3). These observations implicate the possibility that CD might be caused by hyper-sensitive reactions toward pathogenic stimuli.

The immune-suppressive landscape in lepromatous leprosy revealed by single-cell RNA sequencing.

Lepromatous leprosy (L-LEP), caused by the massive proliferation of Mycobacterium leprae primarily in macrophages, is an ideal disease model for investigating the molecular mechanism of intracellular bacteria evading or modulating host immune response. Here, we performed single-cell RNA sequencing of both skin biopsies and peripheral blood mononuclear cells (PBMCs) of L-LEP patients and healthy controls. In L-LEP lesions, we revealed remarkable upregulation of APOE expression that showed a negative correlation with the major histocompatibility complex II gene HLA-DQB2 and MIF, which encodes a pro-inflammatory and anti-microbial cytokine, in the subset of macrophages exhibiting a high expression level of LIPA. The exhaustion of CD8+ T cells featured by the high expression of TIGIT and LAG3 in L-LEP lesions was demonstrated. Moreover, remarkable enhancement of inhibitory immune receptors mediated crosstalk between skin immune cells was observed in L-LEP lesions. For PBMCs, a high expression level of APOE in the HLA-DRhighFBP1high monocyte subset and the expansion of regulatory T cells were found to be associated with L-LEP. These findings revealed the primary suppressive landscape in the L-LEP patients, providing potential targets for the intervention of intracellular bacteria caused persistent infections.

Validation study of HLA-B*13:01 as a biomarker of dapsone hypersensitivity syndrome in leprosy patients in Indonesia.

Leprosy is a stigmatizing, chronic infection which degenerates the nervous system and often leads to incapacitation. Multi-drug therapy which consists of dapsone, rifampicin and clofazimine has been effective to combat this disease. In Indonesia, especially in Papua Island, leprosy is still a problem. Furthermore, there had been higher reports of Dapsone Hypersensitivity Syndrome (DHS) which also challenges leprosy elimination in certain aspects. Globally, DHS has a prevalence rate of 1.4% and a fatality rate up to 13%. The aim of this study is to validate HLA-B*13:01, a previously discovered biomarker for DHS in the Chinese population, as a biomarker for DHS in the Papua population.This is a case-control study of 34 leprosy patients who presented themselves with DHS (case subjects) and 52 leprosy patients without DHS (control subjects). Patients were recruited from 2 provinces: Papua and West Papua. DNA was extracted from 3 ml blood specimens. HLA-B alleles were typed using the gold-standard sequence based typing method. Results were then analysed using logistic regression and risk assessment was carried out. The results of HLA-typing showed that HLA-B*13:01 was the most significant allele associated with DHS, with odds ratio = 233.64 and P-value = 7.11×10-9, confirming the strong association of HLA-B*13:01 to DHS in the Papua population. The sensitivity of this biomarker is 91.2% and specificity is 96.2%, with an area under the curve of 0.95. HLA-B*13:01 is validated as a biomarker for DHS in leprosy patients in Papua, Indonesia, and can potentially be a good predictor of DHS to help prevent this condition in the future.

Evaluation of Prospective HLA-B*13:01 Screening to Prevent Dapsone Hypersensitivity Syndrome in Patients With Leprosy.

Importance: Dapsone hypersensitivity syndrome (DHS) is the most serious adverse reaction associated with dapsone administration and one of the major causes of death in patients with leprosy, whose standard treatment includes multidrug therapy (MDT) with dapsone, rifampicin, and clofazimine. Although the HLA-B*13:01 polymorphism has been identified as the genetic determinant of DHS in the Chinese population, no studies to date have been done to evaluate whether prospective HLA-B*13:01 screening could prevent DHS by identifying patients who should not receive dapsone. Objective: To evaluate the clinical use of prospective HLA-B*13:01 screening for reduction of the incidence of DHS by excluding dapsone from the treatment for patients with HLA-B*13:01-positive leprosy. Design, Setting, and Participants: A prospective cohort study was conducted from February 15, 2015, to April 30, 2018, in 21 provinces throughout China. A total of 1539 patients with newly diagnosed leprosy were enrolled who had not received dapsone previously. After excluding patients who had a history of allergy to sulfones or glucose-6-phosphate dehydrogenase deficiency, 1512 individuals underwent HLA-B*13:01 genotyping. All of the patients were followed up weekly for the first 8 weeks after treatment to monitor for adverse events. Exposures: Patients who were HLA-B*13:01 carriers were instructed to eliminate dapsone from their treatment regimens, and noncarrier patients received standard MDT. Main Outcomes and Measures: The primary outcome was the incidence of DHS. The historical incidence rate of DHS (1.0%) was used as a control. Results: Among 1512 patients (1026 [67.9%] men, 486 [32.1%] women; mean [SD] age, 43.1 [16.2] years), 261 (17.3%) were identified as carriers of the HLA-B*13:01 allele. A total of 714 adverse events in 384 patients were observed during the follow-up period. Dapsone hypersensitivity syndrome did not develop in any of the 1251 patients who were HLA-B*13:01-negative who received dapsone, while approximately 13 patients would be expected to experience DHS, based on the historical incidence rate of 1.0% per year (P = 2.05 × 10-5). No significant correlation was found between other adverse events, including dermatologic or other events, and HLA-B*13:01 status. Conclusions and Relevance: Prospective HLA-B*13:01 screening and subsequent elimination of dapsone from MDT for patients with HLA-B*13:01-positive leprosy may significantly reduce the incidence of DHS in the Chinese population.

Development and evaluation of a droplet digital PCR assay for the diagnosis of paucibacillary leprosy in skin biopsy specimens.

BACKGROUND: The reduced amounts of Mycobacterium leprae (M. leprae) among paucibacillary (PB) patients reflect the need to further optimize methods for leprosy diagnosis. An increasing number of reports have shown that droplet digital polymerase chain reaction (ddPCR) is a promising tool for diagnosis of infectious disease among samples with low copy number. To date, no publications have investigated the utility of ddPCR in the detection of M. leprae. The aim of this study was to develop and evaluate a ddPCR assay for the diagnosis of PB leprosy. METHODOLOGY: The two most sensitive DNA targets for detection of M. leprae were selected from electronic databases for assessment of sensitivity and specificity by quantitative polymerase chain reaction (qPCR) and ddPCR. Control patients (n = 59) suffering from other dermatological diseases were used to define the cut-off of the duplex ddPCR assay. For comparative evaluation, qPCR and ddPCR assays were performed in 44 PB patients and 68 multibacillary (MB) patients. PRINCIPAL FINDINGS: M. leprae-specific repetitive element (RLEP) and groEL (encoding the 65 kDa molecular chaperone GroEL) were used to develop the ddPCR assay by systematically analyzing specificity and sensitivity. Based on the defined cut-off value, the ddPCR assay showed greater sensitivity in detecting M. leprae DNA in PB patients compared with qPCR (79.5% vs 36.4%), while both assays have a 100% sensitivity in MB patients. CONCLUSIONS/SIGNIFICANCE: We developed and evaluated a duplex ddPCR assay for leprosy diagnosis in skin biopsy samples from leprosy patients. While still costly, ddPCR might be a promising diagnostic tool for detection of PB leprosy.

Amino Acid Variants of HLA-DRB1 Confer Susceptibility to Dapsone Hypersensitivity Syndrome in Addition to HLA-B*13:01.

Dapsone hypersensitivity syndrome is a rare yet severe adverse drug reaction caused by dapsone, a principal drug in multidrug therapy for leprosy. HLA-B*13:01 has been identified as a strong risk factor of dapsone hypersensitivity syndrome; however, its low positive predictive value indicated that additional genetic variants may be involved in the disease development. To discover contributing genetic variants within HLA loci in addition to HLA-B*13:01, we performed a high-coverage next-generation sequencing (NGS)-based HLA typing analysis in 103 dapsone-hypersensitive and 857 dapsone-tolerant HLA-B*13:01-positive leprosy patients in a Chinese population. Five amino acid variants in high linkage disequilibrium of HLA-DRB1 were significantly associated with dapsone hypersensitivity syndrome (positions 133, 142, -17, 11, and 13). DRB1*16:02 and DRB1*15:01 tagged by these risk-conferring amino acid residues were associated at a nominal significance level. This study identifies five amino acid variants within HLA-DRB1 that are in high linkage disequilibrium and significantly associated with dapsone hypersensitivity syndrome in a Chinese population.

A pathway-based association analysis identified FMNL1-MAP3K14 as susceptibility genes for leprosy.

The nuclear transcription factor-κB (NF-κB) plays a pivotal role in controlling both innate and adaptive immunity and regulates the expressions of many immunological mediators. Abundant evidences have showed the importance of NF-κB pathway in the host immune responses against Mycobacterium leprae in the development of leprosy. However, no particular association study between leprosy and NF-κB pathway-related gene polymorphisms was reported. Here, we performed a large-scale and two-stage candidate association study to investigate the association between 94 NF-κB pathway-related genes and leprosy. Our results showed that rs58744688 was significantly associated with leprosy (P = 7.57 × 10-7 , OR = 1.12) by combining the previous genomewide association data sets and four independent validation sample series, consisting of a total of 4631 leprosy cases and 6413 healthy controls. This founding implicated that MAP3K14 and FMNL1 were susceptibility genes for leprosy, which suggested the involvement of macrophage targeting and NF-κB pathway in the development of leprosy.

Genome-Wide Analysis of Protein-Coding Variants in Leprosy.

Although genome-wide association studies have greatly advanced our understanding of the contribution of common noncoding variants to leprosy susceptibility, protein-coding variants have not been systematically investigated. We carried out a three-stage genome-wide association study of protein-coding variants in Han Chinese, of whom were 7,048 leprosy patients and 14,398 were healthy control subjects. Seven coding variants of exome-wide significance were discovered, including two rare variants: rs145562243 in NCKIPSD (P = 1.71 × 10-9, odds ratio [OR] = 4.35) and rs149308743 in CARD9 (P = 2.09 × 10-8, OR = 4.75); three low-frequency variants: rs76418789 in IL23R (P = 1.03 × 10-10, OR = 1.36), rs146466242 in FLG (P = 3.39 × 10-12, OR = 1.45), and rs55882956 in TYK2 (P = 1.04 × 10-6, OR = 1.30); and two common variants: rs780668 in SLC29A3 (P = 2.17 × 10-9, OR = 1.14) and rs181206 in IL27 (P = 1.08 × 10-7, OR = 0.83). Discovered protein-coding variants, particularly low-frequency and rare ones, showed involvement of skin barrier and endocytosis/phagocytosis/autophagy, in addition to known innate and adaptive immunity, in the pathogenesis of leprosy, highlighting the merits of protein-coding variant studies for complex diseases.

A large-scale genome-wide association and meta-analysis identified four novel susceptibility loci for leprosy.

Leprosy, a chronic infectious disease, results from the uncultivable pathogen Mycobacterium leprae (M. leprae), and usually progresses to peripheral neuropathy and permanent progressive deformity if not treated. Previously published genetic studies have identified 18 gene/loci significantly associated with leprosy at the genome-wide significant level. However as a complex disease, only a small proportion of leprosy risk could be explained by those gene/loci. To further identify more susceptibility gene/loci, we hereby performed a three-stage GWAS comprising 8,156 leprosy patients and 15,610 controls of Chinese ancestry. Four novel loci were identified including rs6807915 on 3p25.2 (P=1.94 × 10-8, OR=0.89), rs4720118 on 7p14.3 (P=3.85 × 10-10, OR=1.16), rs55894533 on 8p23.1 (P=5.07 × 10-11, OR=1.15) and rs10100465 on 8q24.11 (P=2.85 × 10-11, OR=0.85). Altogether, these findings have provided new insight and significantly expanded our understanding of the genetic basis of leprosy.

Fine-mapping analysis revealed complex pleiotropic effect and tissue-specific regulatory mechanism of TNFSF15 in primary biliary cholangitis, Crohn's disease and leprosy.

Genetic polymorphism within the 9q32 locus is linked with increased risk of several diseases, including Crohn's disease (CD), primary biliary cholangitis (PBC) and leprosy. The most likely disease-causing gene within 9q32 is TNFSF15, which encodes the pro-inflammatory cytokine TNF super-family member 15, but it was unknown whether these disparate diseases were associated with the same genetic variance in 9q32, and how variance within this locus might contribute to pathology. Using genetic data from published studies on CD, PBC and leprosy we revealed that bearing a T allele at rs6478108/rs6478109 (r(2) = 1) or rs4979462 was significantly associated with increased risk of CD and decreased risk of leprosy, while the T allele at rs4979462 was associated with significantly increased risk of PBC. In vitro analyses showed that the rs6478109 genotype significantly affected TNFSF15 expression in cells from whole blood of controls, while functional annotation using publicly-available data revealed the broad cell type/tissue-specific regulatory potential of variance at rs6478109 or rs4979462. In summary, we provide evidence that variance within TNFSF15 has the potential to affect cytokine expression across a range of tissues and thereby contribute to protection from infectious diseases such as leprosy, while increasing the risk of immune-mediated diseases including CD and PBC.

Discovery of six new susceptibility loci and analysis of pleiotropic effects in leprosy.

Genome-wide association studies (GWAS) have led to the discovery of several susceptibility loci for leprosy with robust evidence, providing biological insight into the role of host genetic factors in mycobacterial infection. However, the identified loci only partially explain disease heritability, and additional genetic risk factors remain to be discovered. We performed a 3-stage GWAS of leprosy in the Chinese population using 8,313 cases and 16,017 controls. Besides confirming all previously published loci, we discovered six new susceptibility loci, and further gene prioritization analysis of these loci implicated BATF3, CCDC88B and CIITA-SOCS1 as new susceptibility genes for leprosy. A systematic evaluation of pleiotropic effects demonstrated a high tendency for leprosy susceptibility loci to show association with autoimmunity and inflammatory diseases. Further analysis suggests that molecular sensing of infection might have a similar pathogenic role across these diseases, whereas immune responses have discordant roles in infectious and inflammatory diseases.

An association study of TOLL and CARD with leprosy susceptibility in Chinese population.

Previous genome-wide association studies (GWASs) identified multiple susceptibility loci that have highlighted the important role of TLR (Toll-like receptor) and CARD (caspase recruitment domain) genes in leprosy. A large three-stage candidate gene-based association study of 30 TLR and 47 CARD genes was performed in the leprosy samples of Chinese Han. Of 4363 SNPs investigated, eight SNPs showed suggestive association (P < 0.01) in our previously published GWAS datasets (Stage 1). Of the eight SNPs, rs2735591 and rs4889841 showed significant association (P < 0.001) in an independent series of 1504 cases and 1502 controls (Stage 2), but only rs2735591 (next to BCL10) showed significant association in the second independent series of 938 cases and 5827 controls (Stage 3). Rs2735591 showed consistent association across the three stages (P > 0.05 for heterogeneity test), significant association in the combined validation samples (Pcorrected = 5.54 × 10(-4) after correction for 4363 SNPs tested) and genome-wide significance in the whole GWAS and validation samples (P = 1.03 × 10(-9), OR = 1.24). In addition, we demonstrated the lower expression of BCL10 in leprosy lesions than normal skins and a significant gene connection between BCL10 and the eight previously identified leprosy loci that are associated with NFκB, a major regulator of downstream inflammatory responses, which provides further biological evidence for the association. We have discovered a novel susceptibility locus on 1p22, which implicates BCL10 as a new susceptibility gene for leprosy. Our finding highlights the important role of both innate and adaptive immune responses in leprosy.

Identification of IL18RAP/IL18R1 and IL12B as leprosy risk genes demonstrates shared pathogenesis between inflammation and infectious diseases.

Of eight leprosy susceptibility loci identified by genome-wide association studies, five have been implicated in Crohn disease, suggesting a common genetic fingerprint between leprosy and inflammatory bowel disease (IBD). Here, we conducted a multiple-stage genetic association study of 133 IBD susceptibility loci in multiple leprosy samples (totaling 4,971 leprosy cases and 5,503 controls) from a Chinese population and discovered two associations at rs2058660 on 2q12.1 (p = 4.57 × 10(-19); odds ratio [OR] = 1.30) and rs6871626 on 5q33.3 (p = 3.95 × 10(-18); OR = 0.75), implicating IL18RAP/IL18R1 and IL12B as susceptibility genes for leprosy. Our study reveals the important role of IL12/IL18-mediated transcriptional regulation of IFN-γ production in leprosy, and together with previous findings, it demonstrates the shared genetic susceptibility between infectious and inflammatory diseases.

Identification of two new loci at IL23R and RAB32 that influence susceptibility to leprosy.

We performed a genome-wide association study with 706 individuals with leprosy and 5,581 control individuals and replicated the top 24 SNPs in three independent replication samples, including a total of 3,301 individuals with leprosy and 5,299 control individuals from China. Two loci not previously associated with the disease were identified with genome-wide significance: rs2275606 (combined P = 3.94 × 10(-14), OR = 1.30) on 6q24.3 and rs3762318 (combined P = 3.27 × 10(-11), OR = 0.69) on 1p31.3. These associations implicate IL23R and RAB32 as new susceptibility genes for leprosy. Furthermore, we identified evidence of interaction between the NOD2 and RIPK2 loci, which is consistent with the biological association of the proteins encoded by these genes (NOD2-RIPK2 complex) in activating the NF-κB pathway as a part of the host defense response to infection. Our findings have expanded the biological functions of IL23R by uncovering its involvement in infectious disease susceptibility and suggest a potential involvement of autophagocytosis in leprosy pathogenesis. The IL23R association supports previous observations of the marked overlap of susceptibility genes for leprosy and Crohn's disease, implying common pathogenesis mechanisms.

Gene expression analysis of leprosy by using a multiplex branched DNA assay.

Leprosy is caused by Mycobacterium leprae, and the global registered prevalence of leprosy at the beginning of 2009 stood at 213,036 cases. It has long been thought that leprosy has a strong genetic risk. Recently, we have identified significant associations (P < 1.00 × 10(-10)) between SNPs in the genes CCDC122, C13orf31, NOD2, TNFSF15, HLA-DR and RIPK2 and a trend towards an association (P = 5.10 × 10(-5)) with a SNP in LRRK2. Here, we investigated the expression of these seven genes in formalin-fixed, paraffin-embedded skin tissues of leprosy and matched normal tissues using branched DNA technology. This technology allows for direct measurement of targeted mRNA within cellular lysate using a 96-well plate format in a time frame compared to a reporter gene assay. The clear upregulation of all seven genes was found in leprosy tissues compared to normal tissues, which further supports our genome-wide association study results.