Detection of Mycobacterium leprae DNA from archaeological skeletal remains in Japan using whole genome amplification and polymerase chain reaction

Evolutionary history of Mycobacterium leprae in the Pacific Islands

As one of the oldest known human diseases, leprosy or Hansen's disease remains a public health concern around the world with over 200 000 new cases in 2018. Most human leprosy cases are caused by Mycobacterium leprae, but a small number of cases are now known to be caused by Mycobacterium lepromatosis, a sister taxon of M. leprae. The global pattern of genomic variation in M. leprae is not well defined. Particularly, in the Pacific Islands, the origins of leprosy are disputed. Historically, it has been argued that leprosy arrived on the islands during nineteenth century colonialism, but some oral traditions and palaeopathological evidence suggest an older introduction. To address this, as well as investigate patterns of pathogen exchange across the Pacific Islands, we extracted DNA from 39 formalin-fixed paraffin-embedded biopsy blocks dating to 1992-2016. Using whole-genome enrichment and next-generation sequencing, we produced nine M. leprae genomes dating to 1998-2015 and ranging from 4-63× depth of coverage. Phylogenetic analyses indicate that these strains belong to basal lineages within the M. leprae phylogeny, specifically falling in branches 0 and 5. The phylogeographical patterning and evolutionary dating analysis of these strains support a pre-modern introduction of M. leprae into the Pacific Islands. This article is part of the theme issue 'Insights into health and disease from ancient biomolecules'.

Molecular epidemiology of leprosy: An update

Molecular epidemiology investigations are notoriously challenging in the leprosy field mainly because the inherent characteristics of the disease as well as its yet uncultivated causative agents, Mycobacterium leprae and M. lepromatosis. Despite significant developments in understanding the biology of leprosy bacilli through genomic approaches, the exact mechanisms of transmission is still unclear and the factors underlying pathological variation of the disease in different patients remain as major gaps in our knowledge about leprosy. Despite these difficulties, the last two decades have seen the development of genotyping procedures based on PCR-sequencing of target loci as well as by the genome-wide analysis of an increasing number of geographically diverse isolates of leprosy bacilli. This has provided a foundation for molecular epidemiology studies that are bringing a better understanding of strain evolution associated with ancient human migrations, and phylogeographical insights about the spread of disease globally. This review discusses the advantages and drawbacks of the main tools available for molecular epidemiological investigations of leprosy and summarizes various methods ranging from PCR-based genotyping to genome-typing techniques. We also describe their main applications in analyzing the short-range and long-range transmission of the disease. Finally, we summarise the current gaps and challenges that remain in the field of molecular epidemiology of leprosy.

Genomic Characterization of Mycobacterium leprae to Explore Transmission Patterns Identifies New Subtype in Bangladesh

Mycobacterium leprae, the causative agent of leprosy, is an unculturable bacterium with a considerably reduced genome (3.27 Mb) compared to homologues mycobacteria from the same ancestry. In 2001, the genome of M. leprae was first described and subsequently four genotypes (1-4) and 16 subtypes (A-P) were identified providing means to study global transmission patterns for leprosy. In order to understand the role of asymptomatic carriers we investigated M. leprae carriage as well as infection in leprosy patients (n = 60) and healthy household contacts (HHC; n = 250) from Bangladesh using molecular detection of the bacterial element RLEP in nasal swabs (NS) and slit skin smears (SSS). In parallel, to study M. leprae genotype distribution in Bangladesh we explored strain diversity by whole genome sequencing (WGS) and Sanger sequencing. In the studied cohort in Bangladesh, M. leprae DNA was detected in 33.3% of NS and 22.2% of SSS of patients with bacillary index of 0 whilst in HHC 18.0% of NS and 12.3% of SSS were positive. The majority of the M. leprae strains detected in this study belonged to genotype 1D (55%), followed by 1A (31%). Importantly, WGS allowed the identification of a new M. leprae genotype, designated 1B-Bangladesh (14%), which clustered separately between the 1A and 1B strains. Moreover, we established that the genotype previously designated 1C, is not an independent subtype but clusters within the 1D genotype. Intraindividual differences were present between the M. leprae strains obtained including mutations in hypermutated genes, suggesting mixed colonization/infection or in-host evolution. In summary, we observed that M. leprae is present in asymptomatic contacts of leprosy patients fueling the concept that these individuals contribute to the current intensity of transmission. Our data therefore emphasize the importance of sensitive and specific tools allowing post-exposure prophylaxis targeted at M. leprae-infected or -colonized individuals.

Risk factor for <em>Mycobacterium leprae</em> detection in household contacts with leprosy patients: a study in Papua, East Indonesia

BACKGROUND In the era of leprosy eradication, Jayapura is still one of the biggest leprosy pockets in Papua, Indonesia. The trend for leprosy case detection rate has remained relatively stable over recent years. This study was aimed to detect Mycobacterium leprae in household contacts and to evaluate the associated factors with the detection. METHODS This cross-sectional study recruited household contacts of leprosy patients who were diagnosed consecutively from March to August 2015 in Hamadi Point of Care, Jayapura. The leprosy patients were diagnosed using polymerase chain reaction (PCR). For each leprosy patient, up to four household contacts that had no symptom were included. Every household contact received screening through DNA detection of M. leprae extracted from nasal swab specimens and examined using PCR. Factors for bacteria detection included intensity, time duration and number of contacts living together in the same house, and random blood glucose levels were evaluated. Bivariate analysis was used to associate them with M. leprae detection in household contacts. RESULTS From 107 household contacts of 35 patients who had leprosy, M. leprae was detected in 19.6%. Household contacts with leprosy patients for >1 year was a risk factor for detection (OR = 12.45; 95% CI = 1.595–97.20; p = 0.002). Blood glucose (p = 0.444), ethnic (p = 0.456), sleeping proximity to leprosy case (p = 0.468) and relatives (p = 0.518) give no effect to M. leprae detection in household contacts. CONCLUSIONS Among the various risk factors studied, duration of living together with the patient significantly increased the risk of M. leprae transmission.

A rapid and sensitive system for recovery of nucleic acids from Mycobacteria sp. on archived glass slides

BACKGROUND

The field of diagnostics continues to advance rapidly with a variety of novel approaches, mainly dependent upon high technology platforms. Nonetheless much diagnosis, particularly in developing countries, still relies upon traditional methods such as microscopy. Biological material, particularly nucleic acids, on archived glass slides is a potential source of useful information both for diagnostic and epidemiological purposes. There are significant challenges faced when examining archived samples in order that an adequate amount of amplifiable DNA can be obtained. Herein, we describe a model system to detect low numbers of bacterial cells isolated from glass slides using (laser capture microscopy) LCM coupled with PCR amplification of a suitable target.

RESULTS

Mycobacterium smegmatis was used as a model organism to provide a proof of principle for a method to recover bacteria from a stained sample on a glass slide using a laser capture system. Ziehl-Neelsen (ZN) stained cells were excised and catapulted into tubes. Recovered cells were subjected to DNA extraction and pre-amplified with multiple displacement amplification (MDA). This system allowed a minimum of 30 catapulted cells to be detected following a nested real-time PCR assay, using rpoB specific primers. The combination of MDA and nested real-time PCR resulted in a 30-fold increase in sensitivity for the detection of low numbers of cells isolated using LCM.

CONCLUSIONS

This study highlights the potential of LCM coupled with MDA as a tool to improve the recovery of amplifiable nucleic acids from archived glass slides. The inclusion of the MDA step was essential to enable downstream amplification. This platform should be broadly applicable to a variety of diagnostic applications and we have used it as a proof of principle with a Mycobacterium sp. model system.

Quantitative evaluation of maxillary bone deformation by computed tomography in patients with leprosy

Background

Facial deformation as a sequela of leprosy is caused not only by a saddle nose but also by regression of the maxilla, as well documented in paleopathological observations of excavated skeletal remains of patients with leprosy. However, maxillary changes in living patients have been evaluated only by the subjective visual grading. Here, we attempted to evaluate maxillary bone deformation in patients with leprosy using three-dimensional computed tomography (3D-CT).

Methods

Three-dimensional images centered on the maxilla were reconstructed using multiplanar reconstruction methods in former patients with leprosy (n = 10) and control subjects (n = 5); the anterior-posterior length of the maxilla (M_A-P) was then measured. The difference between the M_A-P of the patients and those of controls was evaluated after compensating for individual skull size. These findings were also compared with those from previous paleopathological studies.

Findings

Three former patients with lepromatous leprosy showed marked atrophy of the maxilla at the prosthion (-8.6, -11.1 and -17.9 mm) which corresponded with the visual appearance of the maxillary deformity, and these results were consistent with paleopathological findings of excavated skeletal remains. Additionally, the precise bone defects of the maxilla could be individually calculated for accurate reconstructive surgery.

Interpretation

We have successfully illustrated maxillary bone deformities in living patients with leprosy. This study also confirmed the maxillary regression described in paleopathological studies.

The sequelae of leprosy are diverse and include facial deformity. A characteristic facial deformity is saddle nose, and the common clinical approach to this facial deformity has been rhinoplasty. Paleopathological studies from excavated skeletal remains have revealed that leprosy also causes deformation of the maxillary bone, which supports the external nose; however, these changes have not been evaluated in living patients with leprosy. In this study, we reconstructed three-dimensional and stereoscopic images of the cranium from CT scan data and evaluated maxillary deformities both qualitatively and quantitatively. Out of the ten former patients with leprosy, three had severe deformity of the maxilla including: atrophy of the anterior nasal spine, resorption of the medial part of the alveolar process, loss of sharpness of the piriform aperture margins, and atrophy of the nasal turbinate and the septum. These changes were consistent with previous findings from paleopathological studies. We also developed a novel method for measuring the anterior-posterior length of the maxilla using 3D-CT images, with compensation for individual differences in cranial size. The calculated values are indicative of the severity of the maxillary defects as well as the thickness of bone that needs to be reconstructed to correct the deformed maxilla.

Inability of 'Whole Genome Amplification' to Improve Success Rates for the Biomolecular Detection of Tuberculosis in Archaeological Samples

We assessed the ability of whole genome amplification (WGA) to improve the efficiency of downstream polymerase chain reactions (PCRs) directed at ancient DNA (aDNA) of members of the Mycobacterium tuberculosis complex (MTBC). Using extracts from a variety of bones and a tooth from human skeletons with or without lesions indicative of tuberculosis, from multiple time periods, we obtained inconsistent results. We conclude that WGA does not provide any advantage in studies of MTBC aDNA. The sporadic nature of our results are probably due to the fact that WGA is itself a PCR-based procedure which, although designed to deal with fragmented DNA, might be inefficient with the low concentration of templates in an aDNA extract. As such, WGA is subject to similar, if not the same, restrictions as PCR when applied to aDNA.

Detection of Mycobacterium ulcerans subsp. shinshuense DNA from a water channel in familial Buruli ulcer cases in Japan

AIM

To determine if Mycobacterium ulcerans subsp. shinshuense is present in the common living environment in a case of familial concurrent occurrence of Buruli ulcer in Japan.

MATERIALS & METHODS

Environmental samples were collected primarily from a water channel in the patient's residence. A combination of whole-genome amplification, touchdown PCR and DNA sequencing was used to detect M. ulcerans subsp. shinshuense DNA.

RESULTS

M. ulcerans subsp. shinshuense DNA was detected in a crayfish sampled in July, but not in other samples including a crayfish sampled in October.

CONCLUSION

These findings support an association between contaminated aquatic environments and Buruli ulcer in Japan. The data also suggest dynamic seasonal appearance of the pathogen in the environment may contribute to the seasonal variation of Buruli ulcer occurrence in Japan.

Molecular analysis of ancient caries

An 84 base pair sequence of the Streptococcus mutans virulence factor, known as dextranase, has been obtained from 10 individuals from the Bronze Age to the Modern Era in Europe and from before and after the colonization in America. Modern samples show four polymorphic sites that have not been found in the ancient samples studied so far. The nucleotide and haplotype diversity of this region have increased over time, which could be reflecting the footprint of a population expansion. While this segment has apparently evolved according to neutral evolution, we have been able to detect one site that is under positive selection pressure both in present and past populations. This study is a first step to study the evolution of this microorganism, analysed using direct evidence obtained from ancient remains.

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.

Paleopathological evidence and detection of Mycobacterium leprae DNA from archaeological skeletal remains of Nabe-kaburi (head-covered with iron pots) burials in Japan

The Nabe-kaburi is a unique burial method, the purpose of which is shrouded in mystery. The burials were performed during the 15^th to 18^th centuries in eastern Japan, and involved covering the heads of the deceased with iron pots or mortars. The identification of leprosy-specific osteological lesions among some of the excavated remains has led to the suggestion that Nabe-kaburi burials were a reflection of the social stigma against certain infectious diseases, such as leprosy, tuberculosis or syphilis. However, molecular evidence for the presence of disease has been lacking. The goal of this study was to detect Mycobacterium leprae (M. leprae) DNA in archaeological human skeletal remains from Nabe-kaburi burials. The paleopathological data from three Nabe-kaburi burials were re-evaluated before small samples were taken from affected and control areas. DNA was extracted and used as a template to target the M. leprae-specific DNA using a combination of whole genome amplification, PCR analysis and DNA sequencing. M. leprae DNA fragments were detected in the two sets of skeletal remains that had also shown paleopathological evidence of leprosy. These findings provide definitive evidence that some of the Nabe-kaburi burials were performed for people affected by leprosy. Demonstration of the presence of M. leprae DNA, combined with archeological and anthropological examinations, will aid in solving the mystery of why Nabe-kaburi burials were performed in medieval Japan.

Comprehensive analysis of prokaryotes in environmental water using DNA microarray analysis and whole genome amplification

The microflora in environmental water consists of a high density and diversity of bacterial species that form the foundation of the water ecosystem. Because the majority of these species cannot be cultured in vitro, a different approach is needed to identify prokaryotes in environmental water. A novel DNA microarray was developed as a simplified detection protocol. Multiple DNA probes were designed against each of the 97,927 sequences in the DNA Data Bank of Japan and mounted on a glass chip in duplicate. Evaluation of the microarray was performed using the DNA extracted from one liter of environmental water samples collected from seven sites in Japan. The extracted DNA was uniformly amplified using whole genome amplification (WGA), labeled with Cy3-conjugated 16S rRNA specific primers and hybridized to the microarray. The microarray successfully identified soil bacteria and environment-specific bacteria clusters. The DNA microarray described herein can be a useful tool in evaluating the diversity of prokaryotes and assessing environmental changes such as global warming.

Palaeopathology and genes: investigating the genetics of infectious diseases in excavated human skeletal remains and mummies from past populations

The aim of this paper is to review the use of genetics in palaeomicrobiology, and to highlight the importance of understanding past diseases. Palaeomicrobiology is the study of disease pathogens in skeletal and mummified remains from archaeological contexts. It has revolutionarised our understanding of health in the past by enabling a deeper knowledge of the origins and evolution of many diseases that have shaped us as a species. Bacterial diseases explored include tuberculosis, leprosy, bubonic plague, typhoid, syphilis, endemic and epidemic typhus, trench fever, and Helicobacter pylori. Viral diseases discussed include influenza, hepatitis B, human papilloma virus (HPV), human T-cell lymphotrophic virus (HTLV-1) and human immunodeficiency virus (HIV). Parasitic diseases investigated include malaria, leishmaniasis, Chagas' disease, roundworm, whipworm, pinworm, Chinese liver fluke, fleas and lice. Through a better understanding of disease origins and their evolution, we can place into context how many infectious diseases are changing over time, and so help us estimate how they may change in the future.

Detection and strain typing of ancient Mycobacterium leprae from a medieval leprosy hospital

Nine burials excavated from the Magdalen Hill Archaeological Research Project (MHARP) in Winchester, UK, showing skeletal signs of lepromatous leprosy (LL) have been studied using a multidisciplinary approach including osteological, geochemical and biomolecular techniques. DNA from Mycobacterium leprae was amplified from all nine skeletons but not from control skeletons devoid of indicative pathology. In several specimens we corroborated the identification of M. leprae with detection of mycolic acids specific to the cell wall of M. leprae and persistent in the skeletal samples. In five cases, the preservation of the material allowed detailed genotyping using single-nucleotide polymorphism (SNP) and multiple locus variable number tandem repeat analysis (MLVA). Three of the five cases proved to be infected with SNP type 3I-1, ancestral to contemporary M. leprae isolates found in southern states of America and likely carried by European migrants. From the remaining two burials we identified, for the first time in the British Isles, the occurrence of SNP type 2F. Stable isotope analysis conducted on tooth enamel taken from two of the type 3I-1 and one of the type 2F remains revealed that all three individuals had probably spent their formative years in the Winchester area. Previously, type 2F has been implicated as the precursor strain that migrated from the Middle East to India and South-East Asia, subsequently evolving to type 1 strains. Thus we show that type 2F had also spread westwards to Britain by the early medieval period.

Current status of leprosy: epidemiology, basic science and clinical perspectives

Leprosy has affected humans for millennia and remains an important health problem worldwide, as evidenced by nearly 250 000 new cases detected every year. It is a chronic infectious disorder, caused by Mycobacterium leprae, that primarily affects the skin and peripheral nerves. Recent advances in basic science have improved our knowledge of the disease. Variation in the cellular immune response is the basis of a range of clinical manifestations. The introduction of multidrug therapy has significantly contributed to a decrease in the prevalence of the disease. However, leprosy control activities, including monitoring and prevention programs, must be maintained.

Mycobacterium leprae: genes, pseudogenes and genetic diversity

Leprosy, which has afflicted human populations for millenia, results from infection with Mycobacterium leprae, an unculturable pathogen with an exceptionally long generation time. Considerable insight into the biology and drug resistance of the leprosy bacillus has been obtained from genomics. M. leprae has undergone reductive evolution and pseudogenes now occupy half of its genome. Comparative genomics of four different strains revealed remarkable conservation of the genome (99.995% identity) yet uncovered 215 polymorphic sites, mainly single nucleotide polymorphisms, and a handful of new pseudogenes. Mapping these polymorphisms in a large panel of strains defined 16 single nucleotide polymorphism-subtypes that showed strong geographical associations and helped retrace the evolution of M. leprae.

Multiple loci variable number tandem repeat (VNTR) analysis (MLVA) of Mycobacterium leprae isolates amplified from European archaeological human remains with lepromatous leprosy

Molecular typing methods based on polymorphisms in single nucleotides and short tandem repeat motifs have been developed as epidemiological typing tools for Mycobacterium leprae. We have used a variable number tandem repeat method based on three variable loci to identify strain variation in archaeological cases of lepromatous leprosy. The panel of polymorphic loci used revealed unique profiles in five cases of leprosy, including those with identical SNP type and subtype. These were also different from profiles of three previously studied lepromatous skeletons. Whilst examination with SNP typing provides evidence for disease origins, dissemination and phylogeny, tandem repeat typing may be useful for studying cases from within a defined area or community where SNP types may be identical due to geographical constraints. We envisage the technique may be useful in studying contemporaneous burials such as those associated with leprosaria and will prove invaluable in authentication of ancient DNA analyses.

Museums and disease: using tissue archive and museum samples to study pathogens

Molecular studies of archival and fossil samples have traditionally focused on the nucleic acids derived from the host species. However, there has recently been an increase in ancient DNA research on the identification and characterization of infectious agents within the hosts. The study of pathogens from the past provides great opportunities for discovering the causes of historical infection events, characterizing host-microorganism co-evolution and directly investigating the evolution of specific pathogens. Several research teams have been able to isolate and characterize a variety of different bacterial, parasite and viral microorganisms. However, this emerging field is not without obstacles. The diagenetic processes that make ancient DNA research generally difficult are also impediments to ancient pathogen research and perhaps more so given that their DNA may represent an even rarer proportion of the remaining nucleic acids in a fossil sample than host DNA. However, studies performed under controlled conditions and following stringent ancient DNA protocols can and have yielded reliable and often surprising results. This article reviews the advantages, problems, and failures of ancient microbiological research.