Ancient genomes reveal a high diversity of Mycobacterium leprae in medieval Europe.

Studying ancient DNA allows us to retrace the evolutionary history of human pathogens, such as Mycobacterium leprae, the main causative agent of leprosy. Leprosy is one of the oldest recorded and most stigmatizing diseases in human history. The disease was prevalent in Europe until the 16th century and is still endemic in many countries with over 200,000 new cases reported annually. Previous worldwide studies on modern and European medieval M. leprae genomes revealed that they cluster into several distinct branches of which two were present in medieval Northwestern Europe. In this study, we analyzed 10 new medieval M. leprae genomes including the so far oldest M. leprae genome from one of the earliest known cases of leprosy in the United Kingdom-a skeleton from the Great Chesterford cemetery with a calibrated age of 415-545 C.E. This dataset provides a genetic time transect of M. leprae diversity in Europe over the past 1500 years. We find M. leprae strains from four distinct branches to be present in the Early Medieval Period, and strains from three different branches were detected within a single cemetery from the High Medieval Period. Altogether these findings suggest a higher genetic diversity of M. leprae strains in medieval Europe at various time points than previously assumed. The resulting more complex picture of the past phylogeography of leprosy in Europe impacts current phylogeographical models of M. leprae dissemination. It suggests alternative models for the past spread of leprosy such as a wide spread prevalence of strains from different branches in Eurasia already in Antiquity or maybe even an origin in Western Eurasia. Furthermore, these results highlight how studying ancient M. leprae strains improves understanding the history of leprosy worldwide.

Dapsone Resistance in Leprosy Patients Originally from American Samoa, United States, 2010-2012.

Skin biopsies from US leprosy patients were tested for mutations associated with drug resistance. Dapsone resistance was found in 4 of 6 biopsies from American Samoa patients. No resistance was observed in patients from other origins. The high rate of dapsone resistance in patients from American Samoa warrants further investigation.

Intrapatient comparison of Mycobacterium leprae by VNTR analysis in nasal secretions and skin biopsy in a Brazilian leprosy endemic region.

Background: This study compares the strains of genotypes of M. leprae from nasal secretions (NS) and skin biopsy (SB) in the same patient, supplementing conventional epidemiology to gain insight into the infection of leprosy in Fortaleza, Brazil. Methods: The sample consisted of 38 newly diagnosed leprosy patients attending the National Reference Center of Dermatology Dona Libania (CDERM), in Fortaleza, who tested positive for M. leprae by PCR in DNA extracts of nasal secretions. DNA was also extracted from skin biopsy (SB) scrapings of each patient and used for multiplex PCR amplification of M. leprae VNTR loci. The number of repeats at 15 loci were determined by the fragment length analysis method. Results: Locus VNTR genotypes were achieved in 38 NS, and in 38 SB specimens. M. leprae strains differed in their genotypes in paired specimens in all but two of 38 patients. The genotype similarity in the remainder ranged from 53% to 87%. Conclusion: M. leprae 15 VNTR loci genotypes of paired nasal and biopsy skin samples from five patients were identical, while as many as seven loci differed in the 33 other patients. When the NS and biopsy genotypes were pooled and compared, it was found that there was a great variability among different VNTR markers. It is important to investigate other molecular markers suitable for typing genetic variations of the bacilli.

Mycobacterium leprae genomes from a British medieval leprosy hospital: towards understanding an ancient epidemic.

BACKGROUND: Leprosy has afflicted humankind throughout history leaving evidence in both early texts and the archaeological record. In Britain, leprosy was widespread throughout the Middle Ages until its gradual and unexplained decline between the 14th and 16th centuries. The nature of this ancient endemic leprosy and its relationship to modern strains is only partly understood. Modern leprosy strains are currently divided into 5 phylogenetic groups, types 0 to 4, each with strong geographical links. Until recently, European strains, both ancient and modern, were thought to be exclusively type 3 strains. However, evidence for type 2 strains, a group normally associated with Central Asia and the Middle East, has recently been found in archaeological samples in Scandinavia and from two skeletons from the medieval leprosy hospital (or leprosarium) of St Mary Magdalen, near Winchester, England. RESULTS: Here we report the genotypic analysis and whole genome sequencing of two further ancient M. leprae genomes extracted from the remains of two individuals, Sk14 and Sk27, that were excavated from 10th-12th century burials at the leprosarium of St Mary Magdalen. DNA was extracted from the surfaces of bones showing osteological signs of leprosy. Known M. leprae polymorphisms were PCR amplified and Sanger sequenced, while draft genomes were generated by enriching for M. leprae DNA, and Illumina sequencing. SNP-typing and phylogenetic analysis of the draft genomes placed both of these ancient strains in the conserved type 2 group, with very few novel SNPs compared to other ancient or modern strains. CONCLUSIONS: The genomes of the two newly sequenced M. leprae strains group firmly with other type 2F strains. Moreover, the M. leprae strain most closely related to one of the strains, Sk14, in the worldwide phylogeny is a contemporaneous ancient St Magdalen skeleton, vividly illustrating the epidemic and clonal nature of leprosy at this site. The prevalence of these type 2 strains indicates that type 2F strains, in contrast to later European and associated North American type 3 isolates, may have been the co-dominant or even the predominant genotype at this location during the 11th century.

Probable zoonotic leprosy in the southern United States.

BACKGROUND: In the southern region of the United States, such as in Louisiana and Texas, there are autochthonous cases of leprosy among native-born Americans with no history of foreign exposure. In the same region, as well as in Mexico, wild armadillos are infected with Mycobacterium leprae. METHODS: Whole-genome resequencing of M. leprae from one wild armadillo and three U.S. patients with leprosy revealed that the infective strains were essentially identical. Comparative genomic analysis of these strains and M. leprae strains from Asia and Brazil identified 51 single-nucleotide polymorphisms and an 11-bp insertion-deletion. We genotyped these polymorphic sites, in combination with 10 variable-number tandem repeats, in M. leprae strains obtained from 33 wild armadillos from five southern states, 50 U.S. outpatients seen at a clinic in Louisiana, and 64 Venezuelan patients, as well as in four foreign reference strains. RESULTS: The M. leprae genotype of patients with foreign exposure generally reflected their country of origin or travel history. However, a unique M. leprae genotype (3I-2-v1) was found in 28 of the 33 wild armadillos and 25 of the 39 U.S. patients who resided in areas where exposure to armadillo-borne M. leprae was possible. This genotype has not been reported elsewhere in the world. CONCLUSIONS: Wild armadillos and many patients with leprosy in the southern United States are infected with the same strain of M. leprae. Armadillos are a large natural reservoir for M. leprae, and leprosy may be a zoonosis in the region. (Funded by the National Institute of Allergy and Infectious Diseases and others.).

Distribution of Mycobacterium leprae strains among cases in a rural and urban population of Maharashtra, India.

The elimination of leprosy continues to be a challenge, with the disease remaining endemic in several countries. India accounts for the highest number of cases, and the identification of child cases indicates recent transmission. Genetic markers, like variable-number tandem repeats (VNTRs) and single-nucleotide polymorphisms (SNPs), have been identified to track transmission of the pathogen Mycobacterium leprae. They were used to describe M. leprae strains detected in 48 skin biopsy specimens from leprosy patients in the state of Maharashtra in western India in rural and urban areas near Mumbai. Ninety-three percent of strains across both settings belonged to the SNP type 1D, with three of SNP type 1B being identified in patients living within 3 km of each other. The VNTR profiles of the Maharashtra strains clustered with those from Southern India reported previously and a few other Asian strains, indicating that the Indian strains are genotypically conserved at the level of many VNTR loci. Taken together, SNP and VNTR markers are sufficiently reliable and suitable for both localized and broad geographical genotype associations. VNTR profiles of additional cases may aid in distinguishing the SNP type 1B and 1D strains.

Molecular drug susceptibility testing and genotyping of Mycobacterium leprae strains from South America.

Possible drug resistance in Mycobacterium leprae strains from Venezuela and three other South American countries was surveyed by molecular methods. None of the 230 strains from new leprosy cases exhibited drug resistance-associated mutations. However, two of the three strains from relapsed cases contained dapsone resistance mutations, and one strain also harbored a rifampin resistance mutation. Single nucleotide polymorphism analysis of these strains revealed five subtypes: 3I (73.8%), 4P (11.6%), 1D (6.9%), 4N (6%), and 4O (1.7%).

What is the evidence that the putative Mycobacterium lepromatosis species causes diffuse lepromatous leprosy?

Han et al. have made a retrospective isolation of DNA from two patients with fatal Lucio's phenomenon. This DNA does have some molecular differences to M. leprae and may constitute a variant of M. leprae. However the experiments and data needed to confirm that this is a new leprosy-causing species have not yet been done. We have outlined the work that does need to be done. For the moment the assertion that 'M. lepromatosis' is a new leprosy-causing species is not proven.

Molecular typing of Mycobacterium leprae strains from northern India using short tandem repeats.

BACKGROUND & OBJECTIVES: Due to the inability to cultivate Mycobacterium leprae in vitro and most cases being paucibacillary, it has been difficult to apply classical genotyping methods to this organism. The objective of this study was therefore, to analyze the diversity among M. leprae strains from Uttar Pradesh, north India, by targeting ten short tandem repeats (STRs) as molecular markers. METHODS: Ninety specimens including 20 biopsies and 70 slit scrappings were collected in TE buffer from leprosy patients, who attended the OPD of National JALMA Institute for Leprosy and Other Mycobacterial Diseases, Tajganj, Agra, and from villages of Model Rural Health Research Unit (MRHRU) at Ghatampur, Kanpur, Uttar Pradesh. DNA was extracted from these specimens and ten STRs loci were amplified by using published and in-house designed primers. The copy numbers were determined by electrophoretic mobility as well as sequence analysis. Phylogenetic analysis was done on variable number of tandem repeats (VNTRs) data sets using start software. RESULTS: Diversity was observed in the cross-sectional survey of isolates obtained from 90 patients. Allelic index for different loci was found to vary from 0.7 to 0.8 except for rpoT for which allelic index was 0.186. Similarity in fingerprinting profiles observed in specimens from the cases from same house or nearby locations indicated a possible common source of infection. Such analysis was also found to be useful in discriminating the relapse from possible reinfection. INTERPRETATION & CONCLUSIONS: This study led to identification of STRs eliciting polymorphism in north Indian strains of M. leprae. The data suggest that these STRs can be used to study the sources and transmission chain in leprosy, which could be very important in monitoring of the disease dynamics in high endemic foci.

Clinicohistopathological correlation in leprosy.

BACKGROUND: Leprosy is a chronic, infectious disease caused by Mycobacterium leprae. It is classified into five groups based on clinical, histological, microbiological and immunological criteria (Ridley and Jopling Classification) . However, a great variation has been observed in the interpretation of histopathological examination ok skin biopsies and clinical presentation of the disease. OBJECTIVE: To correlate clinical diagnosis with histopathological diagnosis of leprosy patients in Nepal. METHODS: A retrospective hospital-based study was conducted among patients with all clinical types of leprosy, classified as per the Ridley-Jopling classification. Skin biopsies were taken from active lesions in all patients and were stained with Hematoxylin and Eosin stain and modified Fite-Ferraco stain for identification of Mycobacterium leprae. The histopathological findings were compared with clinical diagnoses. RESULTS: A total 156 patients were studied, out of which 84 (53.8%) males and 72 (46.1%) females between 8 and 86 years of age. The majority of patients 33 (23.57%) were in the age group of 21-30 years and least affected was children below 10 years 1(0.007%).Overall coincidence of clinical and histopathological diagnoses of classification was seen in 115 cases (80.4%). The maximum correlation (95.2%) was noted in LL patients (p value 0.000049) followed by BT(89.74%), TT (73.2%),BL(72.4%), BB(64.7%). CONCLUSION: Leprosy still continues to be one of the common infectious disease in Nepal and skin biopsy is a useful tool in confirming the clinical diagnosis of leprosy as well as for the therapeutic guide.

Minimally invasive sampling to identify leprosy patients with a high bacterial burden in the Union of the Comoros.

The World Health Organization (WHO) endorsed diagnosis of leprosy (also known as Hansen's disease) entirely based on clinical cardinal signs, without microbiological confirmation, which may lead to late or misdiagnosis. The use of slit skin smears is variable, but lacks sensitivity. In 2017-2018 during the ComLep study, on the island of Anjouan (Union of the Comoros; High priority country according to WHO, 310 patients were diagnosed with leprosy (paucibacillary = 159; multibacillary = 151), of whom 263 were sampled for a skin biopsy and fingerstick blood, and 260 for a minimally-invasive nasal swab. In 74.5% of all skin biopsies and in 15.4% of all nasal swabs, M. leprae DNA was detected. In 63.1% of fingerstick blood samples, M. leprae specific antibodies were detected with the quantitative αPGL-I test. Results show a strong correlation of αPGL-I IgM levels in fingerstick blood and RLEP-qPCR positivity of nasal swabs, with the M. leprae bacterial load measured by RLEP-qPCR of skin biopsies. Patients with a high bacterial load (≥50,000 bacilli in a skin biopsy) can be identified with combination of counting lesions and the αPGL-I test. To our knowledge, this is the first study that compared αPGL-I IgM levels in fingerstick blood with the bacterial load determined by RLEP-qPCR in skin biopsies of leprosy patients. The demonstrated potential of minimally invasive sampling such as fingerstick blood samples to identify high bacterial load persons likely to be accountable for the ongoing transmission, merits further evaluation in follow-up studies.

Molecular epidemiology of Mycobacterium leprae as determined by structure-neighbor clustering.

It has proven challenging to investigate the molecular epidemiology of Mycobacterium leprae, the causative agent of leprosy, due to difficulties with culturing of the organism and a lack of genetic heterogeneity between strains. Recently, a cost-effective panel of variable-number tandem-repeat (VNTR) markers has been developed. Use of this panel allows some of those limitations to be overcome and has allowed the genotyping of 475 M. leprae strains from six different countries. In the present report, we provide a comprehensive analysis of the relationships among the strains in order to investigate the patterns of transmission and migration of M. leprae. We find phylogenetic analysis to be inadequate and have developed an alternative method, structure-neighbor clustering, which assigns isolates with the most similar genotypes to the same groups and, subsequently, subgroups, without inferring how the strains descended from a common ancestor. We validate the approach by using simulated data and detecting expected epidemiological relationships from experimental data. Our results suggest that most M. leprae strains from a given country cluster together and that the occasional isolates assigned to different clusters are a consequence of migration. We found three genetically distinguishable populations among isolates from the Philippines, as well as evidence for the significant influx of strains to that nation from India. We also report that reference strain TN originated from the Philippines and not from India, as was previously believed. Lastly, analysis of isolates from the same families and villages suggests that most community infections originate from a common source or person-to-person transmission but that infection from independent sources does occur with measurable frequency.

Molecular epidemiology and transmission dynamics of leprosy among multicase families and case-contact pairs.

BACKGROUND: Human-to-human transmission of Mycobacterium leprae among household contacts of active leprosy cases is significant, and surveillance of household contacts is vital to interrupting the transmission chain for this disease. This study was conducted to identify similarities in M. leprae strains, based on genomic single nucleotide polymorphisms (SNPs), among cases and their household contacts and in multicase families in order to decipher possible associations, transmission links, various clinical conditions of index cases that enhance person-to-person transmission, and timelines for transmission patterns. METHODS: PCR for M. leprae DNA detection (amplification of the Rlep gene) and SNP subtyping of M. leprae strains was performed for 61 index cases and one of their household contacts. Additionally, we studied six families with multiple cases of leprosy, to understand timelines of infectivity and its relation to severity of the disease in the index cases. RESULTS: Index cases with lepromatous (LL) and borderline lepromatous (BL) leprosy, together with a positive bacteriological index (BI) for M. leprae, result in a higher percentage of their contacts subclinically infected with M. leprae, with odds ratios (OR) of 6.6 (95% confidence interval (CI) 1.6-27.6) for BL and LL, and 7.07 (CI 1.41-35.41) for BI-positive index cases. 75% of the case-contact pairs had a similar SNP subtype of M. leprae. The timeline of infection in multicase families revealed that contacts were infected during the BI-positive period of the index case. CONCLUSION: Using molecular methods, we determined that positivity for M. leprae DNA in contacts of index leprosy cases was attributed to clinical characteristics of leprosy in the index cases. LL and BL forms of leprosy, together with positive BI, contributed to dissemination of infection to household contacts. In conclusion, we found a relationship between SNP subtypes within index case-contact pairs. This method can help decipher the transmission patterns and identify individuals at risk of contracting leprosy.

Genotyping of Mycobacterium leprae for understanding the distribution and transmission of leprosy in endemic provinces of China.

OBJECTIVES: Understanding the nature of Mycobacterium leprae transmission is vital to implement better control strategies for leprosy elimination. The present study expands the knowledge of county-level strain diversity, distribution, and transmission patterns of leprosy in endemic provinces of China. METHODS: We genetically characterized 290 clinical isolates of M. leprae from four endemic provinces using variable number tandem repeats (VNTR) and single nucleotide polymorphisms (SNPs). Attained genetic profiles and cluster consequences were contrasted with geographical and migration features of leprosy at county levels. RESULTS: Considering the allelic variability of 17 VNTR loci by the discriminatory index, (GTA)9, (AT)17, (AT)15, (TA)18, (TTC)21, and (TA)10 are reported to be more highly polymorphic than other loci. The VNTR profile generated the low-density clustering pattern in the counties of Sichuan and Yunnan, whereas clusters have been observed from the isolates from Huayuan (N = 6), Yongding (N = 3), Zixing (N = 3), Chenxi (N = 2) and Zhongfang (N = 2) counties of Hunan, and Zhijin (N = 3), Anlong (N = 2), Zhenning (N = 2), and Xixiu (N = 2) counties of Guizhou. In some clusters, people's social relations have been observed between villages. From the 290 clinical isolates, the most predominantly reported SNP was 3K (278, 95.8%), followed by SNP 1D (10, 3.4%), which are typically observed to be predominant in China. We also detected the novel SNP 3J (2, 0.8%), which has not yet been reported in China. CONCLUSION: The clustering pattern of M. leprae indicates the transmission of leprosy still persists at county levels, suggesting that there is a need to implement better approaches for tracing the close contacts of leprosy patients.

Comparative sequence analysis of Mycobacterium leprae and the new leprosy-causing Mycobacterium lepromatosis.

Mycobacterium lepromatosis is a newly discovered leprosy-causing organism. Preliminary phylogenetic analysis of its 16S rRNA gene and a few other gene segments revealed significant divergence from Mycobacterium leprae, a well-known cause of leprosy, that justifies the status of M. lepromatosis as a new species. In this study we analyzed the sequences of 20 genes and pseudogenes (22,814 nucleotides). Overall, the level of matching of these sequences with M. leprae sequences was 90.9%, which substantiated the species-level difference; the levels of matching for the 16S rRNA genes and 14 protein-encoding genes were 98.0% and 93.1%, respectively, but the level of matching for five pseudogenes was only 79.1%. Five conserved protein-encoding genes were selected to construct phylogenetic trees and to calculate the numbers of synonymous substitutions (dS values) and nonsynonymous substitutions (dN values) in the two species. Robust phylogenetic trees constructed using concatenated alignment of these genes placed M. lepromatosis and M. leprae in a tight cluster with long terminal branches, implying that the divergence occurred long ago. The dS and dN values were also much higher than those for other closest pairs of mycobacteria. The dS values were 14 to 28% of the dS values for M. leprae and Mycobacterium tuberculosis, a more divergent pair of species. These results thus indicate that M. lepromatosis and M. leprae diverged approximately 10 million years ago. The M. lepromatosis pseudogenes analyzed that were also pseudogenes in M. leprae showed nearly neutral evolution, and their relative ages were similar to those of M. leprae pseudogenes, suggesting that they were pseudogenes before divergence. Taken together, the results described above indicate that M. lepromatosis and M. leprae diverged from a common ancestor after the massive gene inactivation event described previously for M. leprae.

Rapid variable-number tandem-repeat genotyping for Mycobacterium leprae clinical specimens.

Mycobacterium leprae is the noncultivable pathogen of leprosy. Since the genome sequence of an isolate of M. leprae has become available, multiple-locus variable-number tandem-repeat (VNTR) analysis (MLVA) has been explored as a tool for strain typing and identification of chains of transmission of leprosy. In order to discover VNTRs and develop methods transferable to clinical samples, MLVA was applied to a global collection of M. leprae isolates derived from leprosy patients and propagated in armadillo hosts. PCR amplification, agarose gel electrophoresis, and sequencing methods were applied to DNA extracts from these infected armadillo tissues (n = 21). We identified polymorphisms in 15 out of 25 short-tandem-repeat (STR) loci previously selected by in silico analyses of the M. leprae genome. We then developed multiplex PCR for amplification of these 15 loci in four separate PCRs suitable for fluorescent fragment length analysis and demonstrated STR profiles highly concordant with those from the sequencing methods. Subsequently, we extended this method to DNA extracts from human clinical specimens, such as skin biopsy specimens (n = 30). With these techniques, mapping of multiple loci and differentiation of genotypes have been possible using total DNA extracts from limited amounts of clinical samples at a reduced cost and with less time. These practical methods are therefore available and applicable to answer focused epidemiological questions and to allow monitoring of the transmission of M. leprae in different countries where leprosy is endemic.

Population-based molecular epidemiology of leprosy in Cebu, Philippines.

To address the persisting problem of leprosy in Cebu, Philippines, we compiled a database of more than 200 patients who attend an established referral skin clinic. We described the patient characteristics in conventional demographic parameters and also applied multiple-locus variable-number tandem-repeat (VNTR) analysis (MLVA) and single nucleotide polymorphism (SNP) typing for Mycobacterium leprae in biopsied skin lesion samples. These combined approaches revealed that transmission is ongoing, with the affected including the young Cebuano population under 40 years of age in both crowded cities and rural areas of the island. The emergence of multicase families (MCF) is indicative of infection unconstrained by standard care measures. For the SNPs, we designed a low-cost PCR-restriction fragment length polymorphism typing method. MLVA in M. leprae was highly discriminatory in this population yet could retain broad groups, as defined by the more stable SNPs, implying temporal marker stability suitable for interpreting population structures and evolution. The majority of isolates belong to an Asian lineage (SNP type 1), and the rest belong to a putative postcolonial lineage (SNP type 3). Specific alleles at two VNTR loci, (GGT)5 and 21-3, were highly associated with SNP type 3 in this population. MLVA identified M. leprae genotype associations for patients with known epidemiological links such as in MCFs and in some villages. These methods provide a molecular database and a rational framework for targeted approaches to search and confirm leprosy transmission in various scenarios.

Quantitative real-time PCR analysis of Mycobacterium leprae DNA and mRNA in human biopsy material from leprosy and reactional cases.

Mycobacterium leprae, the causative agent of leprosy, is uncultivable in defined media. Development of new diagnostic tools which do not depend on growth of bacteria is needed for the early detection of M. leprae and for monitoring the effectiveness of chemotherapy. We used a real-time PCR-based assay to quantify the copy number of bacterial DNA and hsp18 mRNA from 47 leprosy patients using paraffin-embedded biopsy samples. The assay used was specific, sensitive and reproducible. The applicability of this approach in monitoring the chemotherapy of leprosy was examined. A reduction in DNA and mRNA during chemotherapy was observed and hsp18 mRNA could not be detected in patients who underwent 2 years of multidrug therapy (MDT). However, a considerable amount of M. leprae DNA could be detected even after 2 years of MDT. A significant amount of hsp18 mRNA was found in reactional cases as well. This raises important questions regarding the role of bacterial antigens in leprosy reactions and the rationale of omitting antibiotics in the treatment of reactional cases. Results in this study show that real-time PCR could be a better tool for the careful monitoring of bacillary DNA and mRNA in lesions, which will help to improve diagnosis, disease progression and the treatment regimen.