MicroRNAs correlate with bacillary index and genes associated to cell death processes in leprosy

Mycobacterium leprae infects skin and peripheral nerves causing a broad of clinical forms. MicroRNAs (miRNAs) control immune mechanisms such as apoptosis, autophagy as well as to target genes leading to abnormal proliferation, metastasis, and invasion of cells. Herein we evaluated miRNAs expression for leprosy phenotypes in biopsies obtained from patients with and without reactions. We also correlated those miRNAs with both, bacillary index (BI) and genes involved in the micobacteria elimination process. Our results show a significant increase in the miR-125a-3p expression in paucibacillary (PB) patients vs multibacillary (MB) subjects (p = 0.007) and vs reversal reactions (RR) (p = 0.005), respectively. Likewise, there was a higher expression of miR-125a-3p in patients with erythema nodosum leprosum (ENL) vs MB without reactions (p = 0.002). Furthermore, there was a positive correlation between miR-125a-3p, miR-146b-5p and miR-132-5p expression and BI in patients with RR and ENL. These miRNAS were also correlated with genes such as ATG12 (miR-125a-3p), TNFRSF10A (miR-146b-5p), PARK2, CFLAR and STX7 (miR-132-5p). All together we underpin a role for these miRNAs in leprosy pathogenesis, implicating mechanisms such as apoptosis and autophagy in skin. The miR-125a-3p might have a distinct role associated with PB phenotype and ENL in MB patients.

Reanalysis and validation of the transcriptional pleural fluid signature in pleural tuberculosis

Introduction

Pleural tuberculosis (PlTB), the most common site of extrapulmonary TB, is characterized by a paucibacillary nature and a compartmentalized inflammatory response in the pleural cavity, both of which make diagnosis and management extremely challenging. Although transcriptional signatures for pulmonary TB have already been described, data obtained by using this approach for extrapulmonary tuberculosis and, specifically, for pleural tuberculosis are scarce and heterogeneous. In the present study, a set of candidate genes previously described in pulmonary TB was evaluated to identify and validate a transcriptional signature in clinical samples from a Brazilian cohort of PlTB patients and those with other exudative causes of pleural effusion.

Methods

As a first step, target genes were selected by a random forest algorithm with recursive feature elimination (RFE) from public microarray datasets. Then, peripheral blood (PB) and pleural fluid (PF) samples from recruited patients presenting exudative pleural effusion were collected during the thoracentesis procedure. Transcriptional analysis of the selected top 10 genes was performed by quantitative RT-PCR (RT-qPCR).

Results

Reanalysis of the public datasets identified a set of candidate genes (CARD17, BHLHE40, FCGR1A, BATF2, STAT1, BTN3A1, ANKRD22, C1QB, GBP2, and SEPTIN4) that demonstrated a global accuracy of 89.5% in discriminating pulmonary TB cases from other respiratory diseases. Our validation cohort consisted of PlTB (n = 35) patients and non-TB (n = 34) ones. The gene expressions of CARD17, GBP2, and C1QB in PF at diagnosis were significantly different between the two (PlTB and non-TB) groups (p < 0.0001). It was observed that the gene expressions of CARD17 and GBP2 were higher in PlTB PF than in non-TB patients. C1QB showed the opposite behavior, being higher in the non-TB PF. After anti-TB therapy, however, GBP2 gene expression was significantly reduced in PlTB patients (p < 0.001). Finally, the accuracy of the three above-cited highlighted genes in the PF was analyzed, showing AUCs of 91%, 90%, and 85%, respectively. GBP2 was above 80% (sensitivity = 0.89/specificity = 0.81), and CARD17 showed significant specificity (Se = 0.69/Sp = 0.95) in its capacity to discriminate the groups.

Conclusion

CARD17, GBP2, and C1QB showed promise in discriminating PlTB from other causes of exudative pleural effusion by providing accurate diagnoses, thus accelerating the initiation of anti-TB therapy.

MicroRNAs regulating macrophages infected with Leishmania L. (V.) Braziliensis isolated from different clinical forms of American tegumentary leishmaniasis

Background

Leishmaniasis is an infectious disease caused by protozoa of the genus Leishmania. There are still no vaccines, and therapeutic options are limited, indicating the constant need to understand the fine mechanisms of its pathophysiology. An approach that has been explored in leishmaniasis is the participation of microRNAs (miRNAs), a class of small non-coding RNAs that act, in most cases, to repress gene expression. miRNAs play a role in the complex and plastic interaction between the host and pathogens, either as part of the host's immune response to neutralize infection or as a molecular strategy employed by the pathogen to modulate host pathways to its own benefit.

Methods

Monocyte-derived macrophages from healthy subjects were infected with isolates of three clinical forms of L. braziliensis: cutaneous (CL), mucosal (ML), and disseminated (DL) leishmaniasis. We compared the expression of miRNAs that take part in the TLR/NFkB pathways. Correlations with parasite load as well as immune parameters were analyzed.

Results

miRNAs -103a-3p, -21-3p, 125a-3p -155-5p, -146a-5p, -132- 5p, and -147a were differentially expressed in the metastatic ML and DL forms, and there was a direct correlation between miRNAs -103a-3p, -21-3p, -155-5p, -146a-5p, -132-5p, and -9-3p and parasite load with ML and DL isolates. We also found a correlation between the expression of miR-21-3p and miR-146a-5p with the antiapoptotic gene BCL2 and the increase of viable cells, whereas miR-147a was indirectly correlated with CXCL-9 levels.

Conclusion

The expression of miRNAs is strongly correlated with the parasite load and the inflammatory response, suggesting the participation of these molecules in the pathogenesis of the different clinical forms of L. braziliensis.

Gene expression patterns associated with multidrug therapy in multibacillary leprosy

Multidrug therapy (MDT) has been successfully used in the treatment of leprosy. However, although patients are cured after the completion of MDT, leprosy reactions, permanent disability, and occasional relapse/reinfection are frequently observed in patients. The immune system of multibacillary patients (MB) is not able to mount an effective cellular immune response against M. leprae. Consequently, clearance of bacilli from the body is a slow process and after 12 doses of MDT not all MB patients reduce bacillary index (BI). In this context, we recruited MB patients at the uptake and after 12-month of MDT. Patients were stratified according to the level of reduction of the BI after 12 doses MDT. A reduction of at least one log in BI was necessary to be considered a responder patient. We evaluated the pattern of host gene expression in skin samples with RNA sequencing before and after MDT and between samples from patients with or without one log reduction in BI. Our results demonstrated that after 12 doses of MDT there was a reduction in genes associated with lipid metabolism, inflammatory response, and cellular immune response among responders (APOBEC3A, LGALS17A, CXCL13, CXCL9, CALHM6, and IFNG). Also, by comparing MB patients with lower BI reduction versus responder patients, we identified high expression of CDH19, TMPRSS4, PAX3, FA2H, HLA-V, FABP7, and SERPINA11 before MDT. From the most differentially expressed genes, we observed that MDT modulates pathways related to immune response and lipid metabolism in skin cells from MB patients after MDT, with higher expression of genes like CYP11A1, that are associated with cholesterol metabolism in the group with the worst response to treatment. Altogether, the data presented contribute to elucidate gene signatures and identify differentially expressed genes associated with MDT outcomes in MB patients.

The Type I Interferon Pathway Is Upregulated in the Cutaneous Lesions and Blood of Multibacillary Leprosy Patients With Erythema Nodosum Leprosum

In leprosy patients, acute inflammatory episodes, known as erythema nodosum leprosum (ENL), are responsible for high morbidity and tissue damage that occur during the course of Mycobacterium leprae infection. In a previous study, we showed evidence implicating DNA-sensing via TLR9 as an important inflammatory pathway in ENL. A likely important consequence of TLR9 pathway activation is the production of type I interferons (IFN-I) by plasmacytoid dendritic cells (pDCs), also implicated in the pathogenesis of several chronic inflammatory diseases. In this study, we investigated whether the IFN-I pathway is activated during ENL. Blood samples and skin lesions from multibacillary patients diagnosed with ENL were collected and the expression of genes of the IFN-I pathway and interferon-stimulated genes were compared with samples collected from non-reactional multibacillary (NR) patients. Whole blood RNAseq analysis suggested higher activation of the IFN-I pathway in ENL patients, confirmed by RT-qPCR. Likewise, significantly higher mRNA levels of IFN-I-related genes were detected in ENL skin biopsies when compared to NR patient lesions. During thalidomide administration, the drug of choice for ENL treatment, a decrease in the mRNA and protein levels of some of these genes both in the skin and blood was observed. Indeed, in vitro assays showed that thalidomide was able to block the secretion of IFN-I by peripheral blood mononuclear cells in response to M. leprae sonicate or CpG-A, a TLR9 ligand. Finally, the decreased frequencies of peripheral pDCs in ENL patients, along with the higher TLR9 expression in ENL pDCs and the enrichment of CD123⁺ cells in ENL skin lesions, suggest the involvement of these cells as IFN-I producers in this type of reaction. Taken together, our data point to the involvement of the pDC/type I IFN pathway in the pathogenesis of ENL, opening new avenues in identifying biomarkers for early diagnosis and new therapeutic targets for the better management of this reactional episode.

Tissue Specific Dual RNA-Seq Defines Host-Parasite Interplay in Murine Visceral Leishmaniasis Caused by Leishmania donovani and Leishmania infantum

Visceral leishmaniasis is associated with hepato-splenomegaly and altered immune and hematological parameters in both preclinical animal models and humans. We studied mouse experimental visceral leishmaniasis caused by Leishmania infantum and Leishmania donovani in BALB/c mice using dual RNA-seq to investigate the transcriptional response of host and parasite in liver and spleen. We identified only 4 species-specific parasite expressed genes (SSPEGs; log2FC >1, FDR <0.05) in the infected spleen, and none in the infected liver. For the host transcriptome, we found 789 differentially expressed genes (DEGs; log2FC >1, FDR <0.05) in the spleen that were common to both infections, with IFNγ signaling and complement and coagulation cascade pathways highly enriched, and an additional 286 and 186 DEGs that were selective to L. donovani and L. infantum infection, respectively. Among those, there were network interactions between genes of amino acid metabolism and PPAR signaling in L. donovani infection and increased IL1β and positive regulation of fatty acid transport in L. infantum infection, although no pathway enrichment was observed. In the liver, there were 1,939 DEGs in mice infected with either L. infantum or L. donovani in comparison to uninfected mice, and the most enriched pathways were IFNγ signaling, neutrophil mediated immunity, complement and coagulation, cytokine-chemokine responses, and hemostasis. Additionally, 221 DEGs were selective in L. donovani and 429 DEGs in L. infantum infections. These data show that the host response for these two visceral leishmaniasis infection models is broadly similar, and ∼10% of host DEGs vary in infections with either parasite species. IMPORTANCE Visceral leishmaniasis (VL) is caused by two species of Leishmania parasites, L. donovani in the Old World and L. infantum in the New World and countries bordering the Mediterranean. Although cardinal features such as hepato-splenomegaly and alterations in blood and immune function are evident, clinical presentation may vary by geography, with for example severe bleeding often associated with VL in Brazil. Although animal models of both L. donovani and L. infantum have been widely used to study disease pathogenesis, a direct side-by-side comparison of how these parasites species impact the infected host and/or how they might respond to the stresses of mammalian infection has not been previously reported. Identifying common and distinct pathways to pathogenesis will be important to ensure that new therapeutic or prophylactic approaches will be applicable across all forms of VL.

Development and validation of a multiplex real-time qPCR assay using GMP-grade reagents for leprosy diagnosis

Leprosy is a chronic dermato-neurological disease caused by Mycobacterium leprae, an obligate intracellular bacterium. Timely detection is a challenge in leprosy diagnosis, relying on clinical examination and trained health professionals. Furthermore, adequate care and transmission control depend on early and reliable pathogen detection. Here, we describe a qPCR test for routine diagnosis of leprosy-suspected patients. The reaction simultaneously amplifies two specific Mycobacterium leprae targets (16S rRNA and RLEP), and the human 18S rRNA gene as internal control. The limit of detection was estimated to be 2.29 copies of the M. leprae genome. Analytical specificity was evaluated using a panel of 20 other skin pathogenic microorganisms and Mycobacteria, showing no cross-reactivity. Intra- and inter-operator C_p variation was evaluated using dilution curves of M. leprae DNA or a synthetic gene, and no significant difference was observed between three operators in two different laboratories. The multiplex assay was evaluated using 97 patient samples with clinical and histopathological leprosy confirmation, displaying high diagnostic sensitivity (91%) and specificity (100%). Validation tests in an independent panel of 50 samples confirmed sensitivity and specificity of 97% and 98%, respectively. Importantly, assay performance remained stable for at least five months. Our results show that the newly developed multiplex qPCR effectively and specifically detects M. leprae DNA in skin samples, contributing to an efficient diagnosis that expedites the appropriate treatment.

Leprosy is a chronic dermato-neurological disease caused by Mycobacterium leprae, an obligate intracellular bacterium. Diagnosis of leprosy often relies on skin examinations for clinical signs, bacilli staining from skin smears and invasive skin biopsies. However, the spectrum of clinical manifestations and, often, low bacilli numbers can hinder accurate diagnosis. Timely detection is a challenge in leprosy diagnosis, relying on clinical examination and requiring trained health professionals. Proper intervention for adequate care and transmission control depends on early and reliable pathogen detection. Quantitative PCR methods for detecting bacterial DNA are more sensitive and could aid in differentially diagnosing leprosy from other dermatological conditions. In this work, we present a new multiplex PCR that was assessed for quality control standards, and the data indicate that the assay is stable and reproducible. The results presented here are the basis of a novel and robust tool with potential to increase the accuracy of leprosy diagnosis in routine or reference laboratories.

A new paradigm for leprosy diagnosis based on host gene expression

Transcriptional profiling is a powerful tool to investigate and detect human diseases. In this study, we used bulk RNA-sequencing (RNA-Seq) to compare the transcriptomes in skin lesions of leprosy patients or controls affected by other dermal conditions such as granuloma annulare, a confounder for paucibacillary leprosy. We identified five genes capable of accurately distinguishing multibacillary and paucibacillary leprosy from other skin conditions. Indoleamine 2,3-dioxygenase 1 (IDO1) expression alone was highly discriminatory, followed by TLR10, BLK, CD38, and SLAMF7, whereas the HS3ST2 and CD40LG mRNA separated multi- and paucibacillary leprosy. Finally, from the main differentially expressed genes (DEG) and enriched pathways, we conclude that paucibacillary disease is characterized by epithelioid transformation and granuloma formation, with an exacerbated cellular immune response, while multibacillary leprosy features epithelial-mesenchymal transition with phagocytic and lipid biogenesis patterns in the skin. These findings will help catalyze the development of better diagnostic tools and potential host-based therapeutic interventions. Finally, our data may help elucidate host-pathogen interplay driving disease clinical manifestations.

Putative pathogen-selected polymorphisms in the PKLR gene are associated with mycobacterial susceptibility in Brazilian and African populations

Pyruvate kinase (PK), encoded by the PKLR gene, is a key player in glycolysis controlling the integrity of erythrocytes. Due to Plasmodium selection, mutations for PK deficiency, which leads to hemolytic anemia, are associated with resistance to malaria in sub-Saharan Africa and with susceptibility to intracellular pathogens in experimental models. In this case-control study, we enrolled 4,555 individuals and investigated whether PKLR single nucleotide polymorphisms (SNPs) putatively selected for malaria resistance are associated with susceptibility to leprosy across Brazil (Manaus-North; Salvador-Northeast; Rondonópolis-Midwest and Rio de Janeiro-Southeast) and with tuberculosis in Mozambique. Haplotype T/G/G (rs1052176/rs4971072/rs11264359) was associated with leprosy susceptibility in Rio de Janeiro (OR = 2.46, p = 0.00001) and Salvador (OR = 1.57, p = 0.04), and with tuberculosis in Mozambique (OR = 1.52, p = 0.07). This haplotype downregulates PKLR expression in nerve and skin, accordingly to GTEx, and might subtly modulate ferritin and haptoglobin levels in serum. Furthermore, we observed genetic signatures of positive selection in the HCN3 gene (xpEHH>2 -recent selection) in Europe but not in Africa, involving 6 SNPs which are PKLR/HCN3 eQTLs. However, this evidence was not corroborated by the other tests (FST, Tajima's D and iHS). Altogether, we provide evidence that a common PKLR locus in Africans contribute to mycobacterial susceptibility in African descent populations and also highlight, for first, PKLR as a susceptibility gene for leprosy and TB.

Interferon-lambda 3 and 4 Polymorphisms Increase Sustained Virological Responses and Regulate Innate Immunity in Antiviral Therapy With Pegylated Interferon-Alpha

Sustained virologic response (SVR) in chronic hepatitis C (CHC) treatment denotes that the host genetics controls the immune response and unequivocally contribute to viral clearance or disease severity. In this context, single nucleotide polymorphisms (SNPs) in the locus of interferon lambda 3 and 4 genes (IFNL3/4) have been important genetic markers of responsiveness to CHC as prognostic markers for the pegylated-Interferon-alpha/ribavirin (Peg-IFN-α/RBV). Here, we analyzed 12 SNPs at the IFNL3/4 region in 740 treatment-naïve patients with CHC infected with hepatitis C virus (HCV) genotypes 1, 2, or 3 treated with Peg-IFN-α/RBV. Individually, rs12979860-CC, rs8109886-CC, or rs8099917-TT were predictive markers of SVR, while rs12979860-CC demonstrated the stronger effect. Besides, the genotypic combination of these three predictors' genotypes, CC/CC/TT, increased the rate of SVR. Serum levels of cytokines and gene expression analysis on the genes IFNL3, IFNL4, IFNA1, and some of the IFN-stimulated genes (ISGs) were measured in a subgroup of 24 treated patients and 24 healthy volunteers. An antagonist effect was highlighted between the expression of IFNL3/4 and IFNA1 mRNA among patients. Besides, a prominent production of the pro-inflammatory chemokines CCL4 and CXCL10 was observed at a 12-week treatment follow-up. Lower serum levels of these chemokines were detected in patients with an rs12979860-CC genotype associated with the better treatment outcome. Also, lower expression levels of the IFI6, IFI16, IRF9 genes were observed among rs12979860-CC individuals. In conclusion, a combination of the genotypes at the IFNL3/4 locus can act as a better marker for the prognosis for virological responses in an admixed Brazilian population presenting the modulating effect over innate immunity and inflammation that are controlling the outcome of the viral infection, but also other infectious diseases. This study is registered on the ClinicalTrials.gov platform (accession number NCT01889849 and NCT01623336).

Large-Scale Gene Expression Signatures Reveal a Microbicidal Pattern of Activation in Mycobacterium leprae-Infected Monocyte-Derived Macrophages With Low Multiplicity of Infection

Leprosy is a disease with a clinical spectrum of presentations that is also manifested in diverse histological features. At one pole, lepromatous lesions (L-pole) have phagocytic foamy macrophages heavily parasitized with freely multiplying intracellular Mycobacterium leprae. At the other pole, the presence of epithelioid giant cells and granulomatous formation in tuberculoid lesions (T-pole) lead to the control of M. leprae replication and the containment of its spread. The mechanism that triggers this polarization is unknown, but macrophages are central in this process. Over the past few years, leprosy has been studied using large scale techniques to shed light on the basic pathways that, upon infection, rewire the host cellular metabolism and gene expression. M. leprae is particularly peculiar as it invades Schwann cells in the nerves, reprogramming their gene expression leading to a stem-like cell phenotype. This modulatory behavior exerted by M. leprae is also observed in skin macrophages. Here, we used live M. leprae to infect (10:1 multiplicity of infection) monocyte-derived macrophages (MDMs) for 48 h and analyzed the whole gene expression profile using microarrays. In this model, we observe an intense upregulation of genes consistent with a cellular immune response, with enriched pathways including peptide and protein secretion, leukocyte activation, inflammation, and cellular divalent inorganic cation homeostasis. Among the most differentially expressed genes (DEGs) are CCL5/RANTES and CYP27B1, and several members of the metallothionein and metalloproteinase families. This is consistent with a proinflammatory state that would resemble macrophage rewiring toward granulomatous formation observed at the T-pole. Furthermore, a comparison with a dataset retrieved from the Gene Expression Omnibus of M. leprae-infected Schwann cells (MOI 100:1) showed that the patterns among the DEGs are highly distinct, as the Schwann cells under these conditions had a scavenging and phagocytic gene profile similar to M2-like macrophages, with enriched pathways rearrangements in the cytoskeleton, lipid and cholesterol metabolism and upregulated genes including MVK, MSMO1, and LACC1/FAMIN. In summary, macrophages may have a central role in defining the paradigmatic cellular (T-pole) vs. humoral (L-pole) responses and it is likely that the multiplicity of infection and genetic polymorphisms in key genes are gearing this polarization.

Presence of Senescent and Memory CD8+ Leukocytes as Immunocenescence Markers in Skin Lesions of Elderly Leprosy Patients

Leprosy is an infectious disease that remains endemic in approximately 100 developing countries, where about 200,000 new cases are diagnosed each year. Moreover, multibacillary leprosy, the most contagious form of the disease, has been detected at continuously higher rates among Brazilian elderly people. Due to the so-called immunosenescence, characterized by several alterations in the quality of the immune response during aging, this group is more susceptible to infectious diseases. In view of such data, the purpose of our work was to investigate if age-related alterations in the immune response could influence the pathogenesis of leprosy. As such, we studied 87 individuals, 62 newly diagnosed and untreated leprosy patients distributed according to the age range and to the clinical forms of the disease and 25 healthy volunteers, who were studied as controls. The frequency of senescent and memory CD8⁺ leukocytes was assessed by immunofluorescence of biopsies from cutaneous lesions, while the serum levels of IgG anti-CMV antibodies were analyzed by chemiluminescence and the gene expression of T cell receptors' inhibitors by RT-qPCR. We noted an accumulation of memory CD8⁺ T lymphocytes, as well as reduced CD8⁺CD28⁺ cell expression in skin lesions from elderly patients, when compared to younger people. Alterations in LAG3 and PDCD1 gene expression in cutaneous lesions of young MB patients were also observed, when compared to elderly patients. Such data suggest that the age-related alterations of T lymphocyte subsets can facilitate the onset of leprosy in elderly patients, not to mention other chronic inflammatory diseases.

Ultra-sensitive detection of Mycobacterium leprae: DNA extraction and PCR assays

Leprosy urgently needs a precise and early diagnostic tool. The sensitivity of the direct (bacilli staining, Mycobacterium leprae DNA) and indirect (antibody levels, T cell assays) diagnostics methods vary based on the clinical form. Recently, PCR-based M. leprae DNA detection has been shown to differentially diagnose leprosy from other dermatological conditions. However, accuracy can still be improved, especially for use with less invasive clinical samples. We tested different commercial DNA extraction kits: DNeasy Blood & Tissue, QIAamp DNA Microbiome, Maxwell 16 DNA Purification, PowerSoil DNA Isolation; as well as in-house phenol-chloroform and Trizol/FastPrep methods. Extraction was performed on M. leprae-infected mouse footpads and different clinical samples of leprosy patients (skin biopsies and scrapings, lesion, oral and nasal swabs, body hair, blood on FTA cards, peripheral whole blood). We observed that the Microbiome kit was able to enrich for mycobacterial DNA, most likely due the enzymatic digestion cocktail along with mechanical disruption involved in this method. Consequently, we had a significant increase in sensitivity in skin biopsies from paucibacillary leprosy patients using a duplex qPCR targeting 16S rRNA (M. leprae) and 18S rRNA (mammal) in the StepOnePlus system. Our data showed that the presence of M. leprae DNA was best detected in skin biopsies and skin scrapings, independent of the extraction method or the clinical form. For multibacillary patients, detection of M. leprae DNA in nasal swabs indicates the possibility of having a much less invasive sample that can be used for the purposes of DNA sequencing for relapse analysis and drug resistance monitoring. Overall, DNA extracted with the Microbiome kit presented the best bacilli detection rate for paucibacillary cases, indicating that investments in extraction methods with mechanical and DNA digestion should be made.

Population Genomics of Mycobacterium leprae Reveals a New Genotype in Madagascar and the Comoros

Human settlement of Madagascar traces back to the beginning of the first millennium with the arrival of Austronesians from Southeast Asia, followed by migrations from Africa and the Middle East. Remains of these different cultural, genetic, and linguistic legacies are still present in Madagascar and other islands of the Indian Ocean. The close relationship between human migration and the introduction and spread of infectious diseases, a well-documented phenomenon, is particularly evident for the causative agent of leprosy, Mycobacterium leprae. In this study, we used whole-genome sequencing (WGS) and molecular dating to characterize the genetic background and retrace the origin of the M. leprae strains circulating in Madagascar (n = 30) and the Comoros (n = 3), two islands where leprosy is still considered a public health problem and monitored as part of a drug resistance surveillance program. Most M. leprae strains (97%) from Madagascar and Comoros belonged to a new genotype as part of branch 1, closely related to single nucleotide polymorphism (SNP) type 1D, named 1D-Malagasy. Other strains belonged to the genotype 1A (3%). We sequenced 39 strains from nine other countries, which, together with previously published genomes, amounted to 242 genomes that were used for molecular dating. Specific SNP markers for the new 1D-Malagasy genotype were used to screen samples from 11 countries and revealed this genotype to be restricted to Madagascar, with the sole exception being a strain from Malawi. The overall analysis thus ruled out a possible introduction of leprosy by the Austronesian settlers and suggests a later origin from East Africa, the Middle East, or South Asia.