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.

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.

Disruption of LRRK2 in Zebrafish leads to hyperactivity and weakened antibacterial response.

As a protein with complex domain structure and roles in kinase, GTPase and scaffolding, LRRK2 is believed to be an important orchestration node leading to several cascades of signal transduction rather than one specific pathway. LRRK2 variants were found to be associated with Parkinson's disease, Crohn's disease and leprosy. Here we disrupt LRRK2 in zebrafish and found hyperactivity rather than hypoactivity in adult zebrafish mutants. By RNA-seq we found genes involved in infectious disease and immunological disease were notably affected. Functional studies also revealed a weakened antibacterial response in LRRK2 mutant. This mutant can be further explored for revealing molecular mechanisms and modeling of LRRK2 related diseases.

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.

Genome-wide analysis of the genetic regulation of gene expression in human neutrophils.

Neutrophils are an abundant immune cell type involved in both antimicrobial defence and autoimmunity. The regulation of their gene expression, however, is still largely unknown. Here we report an eQTL study on isolated neutrophils from 114 healthy individuals of Chinese ethnicity, identifying 21,210 eQTLs on 832 unique genes. Unsupervised clustering analysis of these eQTLs confirms their role in inflammatory responses and immunological diseases but also indicates strong involvement in dermatological pathologies. One of the strongest eQTL identified (rs2058660) is also the tagSNP of a linkage block reported to affect leprosy and Crohn's disease in opposite directions. In a functional study, we can link the C allele with low expression of the ß-chain of IL18-receptor (IL18RAP). In neutrophils, this results in a reduced responsiveness to IL-18, detected both on the RNA and protein level. Thus, the polymorphic regulation of human neutrophils can impact beneficial as well as pathological inflammatory responses.

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.

Chromosome 2p14 is linked to susceptibility to leprosy.

BACKGROUND: A genetic component to the etiology of leprosy is well recognized but the mechanism of inheritance and the genes involved are yet to be fully established. METHODOLOGY: A genome-wide single nucleotide polymorphism (SNP) based linkage analysis was carried out using 23 pedigrees, each with 3 to 7 family members affected by leprosy. Multipoint parametric and non-parametric linkage analyses were performed using MERLIN 1.1.1. PRINCIPAL FINDINGS: Genome-wide significant evidence for linkage was identified on chromosome 2p14 with a heterogeneity logarithm of odds (HLOD) score of 3.51 (rs1106577) under a recessive model of inheritance, while suggestive evidence was identified on chr.4q22 (HLOD 2.92, rs1349350, dominant model), chr. 8q24 (HLOD 2.74, rs1618523, recessive model) and chr.16q24 (HLOD 1.93, rs276990 dominant model). Our study also provided moderate evidence for a linkage locus on chromosome 6q24-26 by non-parametric linkage analysis (rs6570858, LOD 1.54, p = 0.004), overlapping a previously reported linkage region on chromosome 6q25-26. CONCLUSION: A genome-wide linkage analysis has identified a new linkage locus on chromosome 2p14 for leprosy in Pedigrees from China.

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.

Evidence for an association of HLA-DRB1*15 and DRB1*09 with leprosy and the impact of DRB1*09 on disease onset in a Chinese Han population.

BACKGROUND: Human leukocyte antigens (HLAs) have been proposed to modulate the immune response to Mycobacterium leprae. The association of HLA-DRB1 with leprosy has been reported in several populations, but not in a Chinese population. METHODS: The polymerase chain reaction-sequence-specific oligonucleotide probe with Luminex100 (PCR-SSOP-Luminex) method was used to genotype HLA-DRB1 alleles in 305 leprosy patients and 527 healthy control individuals. RESULTS: The HLA-DRB1*15 allele was significantly more prevalent among leprosy patients than healthy controls, whereas the frequency of the HLA-DRB1*09 allele was lower among leprosy patients, especially those with early-onset disease. CONCLUSION: HLA-DRB1 alleles are associated with leprosy susceptibility in a Chinese population. The HLA-DRB1*09 allele was found to be protective exclusively in a subset of early-onset leprosy patients.