CYBB-Mediated Ferroptosis Associated with Immunosuppression in Mycobacterium leprae-Infected Monocyte-Derived Macrophages.

Mycobacterium leprae-infected macrophages preferentially exhibit the regulatory M2 phenotype in vitro, which helps the immune escape unabated growth of M leprae in host cells. The mechanism that triggers macrophage polarization is still unknown. In this study, we performed single-cell RNA sequencing to determine the initial responses of human monocyte-derived macrophages against M leprae infection of 4 healthy individuals and found an increase in a major alternative-activated macrophage type that overexpressed NEAT1, CCL2, and CD163. Importantly, further functional analysis showed that ferroptosis was positively correlated with M2 polarization of macrophages, and in vitro experiments have shown that inhibition of ferroptosis promotes the survival of M leprae within macrophages. In addition, further joint analysis of our results with mutisequencing data from patients with leprosy and in vitro validation identified that CYBB was the pivotal molecule for ferroptosis that could promote the M2 polarization of M leprae-infected macrophages, resulting in the immune escape and unabated growth of pathogenic bacteria. Overall, our results suggest that M leprae facilitated its survival by inducing CYBB-mediated macrophage ferroptosis leading to its alternative activation and might reveal the potential for a new therapeutic strategy of leprosy.

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.

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.

Massively Parallel Sequencingof the Filaggrin Gene Reveals an Association Between FLG Loss-of-function Mutations and Leprosy.

Filaggrin, encoded by the FLG gene, plays a crucial role in the barrier function of epidermis, but the association between FLG loss-of-function mutations and infectious skin diseases has not been systematically studied. FLG coding sequences from 945 patients with leprosy and 916 healthy controls were captured and enriched using an array-based high-throughput system, and subjected to next-generation sequencing. The loss-of-function mutations found were further validated by Sanger sequencing. A total of 21 loss-of-function mutations were found in 945 patients with leprosy, with a carrier rate of 17.53%, while the prevalence of these mutations in 916 healthy controls was 14.77%, which was significantly lower than in patients. Two individual FLG loss-of-function mutations (K4022X and Q1790X) were found to be significantly associated with leprosy. These results suggest a possible role for filaggrin in defending against leprosy pathogens.

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.

Prediction of leprosy in the Chinese population based on a weighted genetic risk score.

Genome wide association studies (GWASs) have revealed multiple genetic variants associated with leprosy in the Chinese population. The aim of our study was to utilize the genetic variants to construct a risk prediction model through a weighted genetic risk score (GRS) in a Chinese set and to further assess the performance of the model in identifying higher-risk contact individuals in an independent set. The highest prediction accuracy, with an area under the curve (AUC) of 0.743 (95% confidence interval (CI): 0.729-0.757), was achieved with a GRS encompassing 25 GWAS variants in a discovery set that included 2,144 people affected by leprosy and 2,671 controls. Individuals in the high-risk group, based on genetic factors (GRS > 28.06), have a 24.65 higher odds ratio (OR) for developing leprosy relative to those in the low-risk group (GRS≤18.17). The model was then applied to a validation set consisting of 1,385 people affected by leprosy and 7,541 individuals in contact with leprosy, which yielded a discriminatory ability with an AUC of 0.707 (95% CI: 0.691-0.723). When a GRS cut-off value of 22.38 was selected with the optimal sensitivity and specificity, it was found that 39.31% of high risk contact individuals should be screened in order to detect leprosy in 64.9% of those people affected by leprosy. In summary, we developed and validated a risk model for the prediction of leprosy that showed good discrimination capabilities, which may help physicians in the identification of patients coming into contact with leprosy and are at a higher-risk of developing this condition.

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.

Tuberculosis risk-associated single nucleotide polymorphisms do not show association with leprosy in Chinese population.

OBJECTIVE: Leprosy and tuberculosis (TB) are chronic granulomatous infectious diseases. As well as pathogen and environmental factors, host genetic factors make a substantial contribution to susceptibility to both diseases. More importantly, leprosy and TB also have pathogenic mechanisms and clinical features in common. In this study, the genetic association between leprosy and TB was investigated in a Chinese Han population. METHODS: A genetic association study that included 46 TB susceptibility single nucleotide polymorphisms (SNPs) was performed, involving 1150 leprosy cases and 1150 controls from the Chinese Han population. The Sequenom MassARRAY system was used. RESULTS: No significant association was found between the 46 SNPs and leprosy. Therefore, according to the present study, there is no shared susceptibility locus between leprosy and TB in the Chinese Han population. CONCLUSIONS: Although leprosy and TB have a number of similar characteristics, no shared susceptibility loci were found in the Chinese Han population. Thus, this study demonstrated that the genetic basis of the pathogenesis of the two diseases may vary greatly.

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.