Genome-wide comparison of medieval and modern Mycobacterium leprae

Ancient Mycobacterium leprae genome reveals medieval English red squirrels as animal leprosy host

Leprosy, one of the oldest recorded diseases in human history, remains prevalent in Asia, Africa, and South America, with over 200,000 cases every year.1,2 Although ancient DNA (aDNA) approaches on the major causative agent, Mycobacterium leprae, have elucidated the disease's evolutionary history,3,4,5 the role of animal hosts and interspecies transmission in the past remains unexplored. Research has uncovered relationships between medieval strains isolated from archaeological human remains and modern animal hosts such as the red squirrel in England.6,7 However, the time frame, distribution, and direction of transmissions remains unknown. Here, we studied 25 human and 12 squirrel samples from two archaeological sites in Winchester, a medieval English city well known for its leprosarium and connections to the fur trade. We reconstructed four medieval M. leprae genomes, including one from a red squirrel, at a 2.2-fold average coverage. Our analysis revealed a phylogenetic placement of all strains on branch 3 as well as a close relationship between the squirrel strain and one newly reconstructed medieval human strain. In particular, the medieval squirrel strain is more closely related to some medieval human strains from Winchester than to modern red squirrel strains from England, indicating a yet-undetected circulation of M. leprae in non-human hosts in the Middle Ages. Our study represents the first One Health approach for M. leprae in archaeology, which is centered around a medieval animal host strain, and highlights the future capability of such approaches to understand the disease's zoonotic past and current potential.

Sequential trypsin and ProAlanase digestions unearth immunological protein biomarkers shrouded by skeletal collagen

This study investigates the efficacy of proteomic analysis of human remains to identify active infections in the past through the detection of pathogens and the host response to infection. We advance leprosy as a case study due to the sequestering of sufferers in leprosaria and the suggestive skeletal lesions that can result from the disease. Here we present a sequential enzyme extraction protocol, using trypsin followed by ProAlanase, to reduce the abundance of collagen peptides and in so doing increase the detection of non-collagenous proteins. Through our study of five individuals from an 11th to 18th century leprosarium, as well as four from a contemporaneous non-leprosy associated cemetery in Barcelona, we show that samples from 2 out of 5 leprosarium individuals extracted with the sequential digestion methodology contain numerous host immune proteins associated with modern leprosy. In contrast, individuals from the non-leprosy associated cemetery and all samples extracted with a trypsin-only protocol did not. Through this study, we advance a palaeoproteomic methodology to gain insights into the health of archaeological individuals and take a step toward a proteomics-based method to study immune responses in past populations.

Bioarchaeological investigation of individuals with suspected multibacillary leprosy from the mediaeval leprosarium of St Mary Magdalen, Winchester, Hampshire, UK

Introduction. We have examined four burials from the St Mary Magdalen mediaeval leprosarium cemetery in Winchester, Hampshire, UK. One (Sk.8) was a male child, two (Sk.45 and Sk.52) were adolescent females and the fourth (Sk.512) was an adult male. The cemetery was in use between the 10th and 12th centuries. All showed skeletal lesions of leprosy. Additionally, one of the two females (Sk.45) had lesions suggestive of multi-cystic tuberculosis and the second (Sk.52) of leprogenic odontodysplasia (LO), a rare malformation of the roots of the permanent maxillary incisors.Gap statement. Relatively little is known of the manifestations of lepromatous leprosy (LL) in younger individuals from the archaeological record.Aims and Methodology. To address this, we have used ancient DNA testing and osteological examination of the individuals, supplemented with X-ray and microcomputed tomography (micro-CT) scan as necessary to assess the disease status.Results and Conclusions. The presence of Mycobacterium leprae DNA was confirmed in both females, and genotyping showed SNP type 3I-1 strains but with a clear genotypic variation. We could not confirm Mycobacterium tuberculosis complex DNA in the female individual SK.45. High levels of M. leprae DNA were found within the pulp cavities of four maxillary teeth from the male child (Sk.8) with LO, consistent with the theory that the replication of M. leprae in alveolar bone may interfere with root formation at key stages of development. We report our biomolecular findings in these individuals and review the evidence this site has contributed to our knowledge of mediaeval leprosy.

Equine herpesvirus 4 infected domestic horses associated with Sintashta spoke-wheeled chariots around 4,000 years ago

Equine viral outbreaks have disrupted the socio-economic life of past human societies up until the late 19th century and continue to be of major concern to the horse industry today. With a seroprevalence of 60-80 per cent, equine herpesvirus 4 (EHV-4) is the most common horse pathogen on the planet. Yet, its evolutionary history remains understudied. Here, we screen the sequenced data of 264 archaeological horse remains to detect the presence of EHV-4. We recover the first ancient EHV-4 genome with 4.2× average depth-of-coverage from a specimen excavated in the Southeastern Urals and dated to the Early Bronze Age period, approximately 3,900 years ago. The recovery of an EHV-4 virus outside the upper respiratory tract not only points to an animal particularly infected but also highlights the importance of post-cranial bones in pathogen characterisation. Bayesian phylogenetic reconstruction provides a minimal time estimate for EHV-4 diversification to around 4,000 years ago, a time when modern domestic horses spread across the Central Asian steppes together with spoke-wheeled Sintashta chariots, or earlier. The analyses also considerably revise the diversification time of the two EHV-4 subclades from the 16th century based solely on modern data to nearly a thousand years ago. Our study paves the way for a robust reconstruction of the history of non-human pathogens and their impact on animal health.

Multiplex PCR-based RFLP assay for early identification of prevalent Mycobacterium leprae genotypes

Mycobacterium leprae is classified into four SNP genotypes and 16 subtypes (from 1A to 4P) that exhibit phylogeographical association reported from around the world. Among them, genotypes 1D and 3I represent more than 60% of M. leprae strains. Here, we report a new method for M. leprae genotyping which identifies the genotypes 1D and 3I by combining multiplex PCR amplification and restriction fragment length polymorphism (RFLP) of a M. leprae DNA amplicons using AgeI restriction enzyme. Agarose gel electrophoresis showed a deletion of 11 bp only among 3I genotypes by electrophoresis. When this multiplex PCR reaction is subjected to AgeI digestion, successful restriction digestion shows three bands for all the genotypes except 1D where only two bands were observed due to loss of restriction site. This method gives us the advantage of 1-step identification of the two most prevalent strains of M. leprae without using specialized equipments such as the Sanger sequencing system or quantitative PCR.

Professor David Minnikin Memorial Lecture: An era of the mycobacterial cell wall lipid biomarkers

This paper is dedicated to the memory of Professor David Ernest Minnikin (1939-2021). David was one of the key scientists who pioneered the field of Mycobacterium tuberculosis cell envelope research for over half a century. From the classification, identification, and extraction of the unusual lipids of the mycobacterial cell wall, to exploiting them as characteristic lipid biomarkers for sensitive detection, his ideas enlightened a whole world of possibilities within the tuberculosis (TB) field. In addition, his definition of the intricate models now forms a key milestone in our understanding of the M. tuberculosis cell envelope and has resolved many unanswered questions on the evolution of M. tuberculosis.

Benchmarking Metagenomic Classifiers on Simulated Ancient and Modern Metagenomic Data

Taxonomic profiling of ancient metagenomic samples is challenging due to the accumulation of specific damage patterns on DNA over time. Although a number of methods for metagenome profiling have been developed, most of them have been assessed on modern metagenomes or simulated metagenomes mimicking modern metagenomes. Further, a comparative assessment of metagenome profilers on simulated metagenomes representing a spectrum of degradation depth, from the extremity of ancient (most degraded) to current or modern (not degraded) metagenomes, has not yet been performed. To understand the strengths and weaknesses of different metagenome profilers, we performed their comprehensive evaluation on simulated metagenomes representing human dental calculus microbiome, with the level of DNA damage successively raised to mimic modern to ancient metagenomes. All classes of profilers, namely, DNA-to-DNA, DNA-to-protein, and DNA-to-marker comparison-based profilers were evaluated on metagenomes with varying levels of damage simulating deamination, fragmentation, and contamination. Our results revealed that, compared to deamination and fragmentation, human and environmental contamination of ancient DNA (with modern DNA) has the most pronounced effect on the performance of each profiler. Further, the DNA-to-DNA (e.g., Kraken2, Bracken) and DNA-to-marker (e.g., MetaPhlAn4) based profiling approaches showed complementary strengths, which can be leveraged to elevate the state-of-the-art of ancient metagenome profiling.

Specialized active leprosy search strategies in an endemic area of the Brazilian Amazon identifies a hypermutated Mycobacterium leprae strain causing primary drug resistance

Introduction

Leprosy, an infectious disease caused by Mycobacterium leprae, remains a public health concern in endemic countries, particularly in Brazil. In this study, we conducted an active surveillance campaign in the hyperendemic city of Castanhal in the northeastern part of the state of Pará using clinical signs and symptoms combined with serological and molecular tools to diagnose new cases and to identify drug resistance of circulating M. leprae strains and their distribution in the community.

Methods

During an active surveillance of one week, we enrolled 318 individuals using three different strategies to enroll subjects for this study: (i) an active survey of previously treated cases from 2006 to 2016 found in the Brazil National Notifiable Disease Information System database (n = 23) and their healthy household contacts (HHC) (n = 57); (ii) an active survey of school children (SC) from two primary public schools in low-income neighborhoods (n = 178), followed by visits to the houses of these newly diagnosed SC (n = 7) to examine their HHC (n = 34) where we diagnosed additional new cases (n = 6); (iii) and those people who spontaneously presented themselves to our team or the local health center with clinical signs and/or symptoms of leprosy (n = 6) with subsequent follow-up of their HHC when the case was confirmed (n = 20) where we diagnosed two additional cases (n = 2). Individuals received a dermato-neurological examination, 5 ml of peripheral blood was collected to assess the anti-PGL-I titer by ELISA and intradermal earlobe skin scrapings were taken from HHC and cases for amplification of the M. leprae RLEP region by qPCR.

Results

Anti-PGL-I positivity was highest in the new leprosy case group (52%) followed by the treated group (40.9%), HHC (40%) and lowest in SC (24.6%). RLEP qPCR from SSS was performed on 124 individuals, 22 in treated cases, 24 in newly diagnosed leprosy cases, and 78 in HHC. We detected 29.0% (36/124) positivity overall in this sample set. The positivity in treated cases was 31.8% (7/22), while in newly diagnosed leprosy cases the number of positives were higher, 45.8% (11/23) and lower in HHC at 23.7% (18/76). Whole genome sequencing of M. leprae from biopsies of three infected individuals from one extended family revealed a hypermutated M. leprae strain in an unusual case of primary drug resistance while the other two strains were drug sensitive.

Discussion

This study represents the extent of leprosy in an active surveillance campaign during a single week in the city of Castanhal, a city that we have previously surveyed several times during the past ten years. Our results indicate the continuing high transmission of leprosy that includes fairly high rates of new cases detected in children indicating recent spread by multiple foci of infection in the community. An unusual case of a hypermutated M. leprae strain in a case of primary drug resistance was discovered. It also revealed a high hidden prevalence of overt disease and subclinical infection that remains a challenge for correct clinical diagnosis by signs and symptoms that may be aided using adjunct laboratory tests, such as RLEP qPCR and anti-PGL-I serology.

Ancient Clostridium DNA and variants of tetanus neurotoxins associated with human archaeological remains

The analysis of microbial genomes from human archaeological samples offers a historic snapshot of ancient pathogens and provides insights into the origins of modern infectious diseases. Here, we analyze metagenomic datasets from 38 human archaeological samples and identify bacterial genomic sequences related to modern-day Clostridium tetani, which produces the tetanus neurotoxin (TeNT) and causes the disease tetanus. These genomic assemblies had varying levels of completeness, and a subset of them displayed hallmarks of ancient DNA damage. Phylogenetic analyses revealed known C. tetani clades as well as potentially new Clostridium lineages closely related to C. tetani. The genomic assemblies encode 13 TeNT variants with unique substitution profiles, including a subgroup of TeNT variants found exclusively in ancient samples from South America. We experimentally tested a TeNT variant selected from an ancient Chilean mummy sample and found that it induced tetanus muscle paralysis in mice, with potency comparable to modern TeNT. Thus, our ancient DNA analysis identifies DNA from neurotoxigenic C. tetani in archaeological human samples, and a novel variant of TeNT that can cause disease in mammals.

Hi-plex deep amplicon sequencing for identification, high-resolution genotyping and multidrug resistance prediction of Mycobacterium leprae directly from patient biopsies by using Deeplex Myc-Lep

BACKGROUND

Expansion of antimicrobial resistance monitoring and epidemiological surveillance are key components of the WHO strategy towards zero leprosy. The inability to grow Mycobacterium leprae in vitro precludes routine phenotypic drug susceptibility testing, and only limited molecular tests are available. We evaluated a culture-free targeted deep sequencing assay, for mycobacterial identification, genotyping based on 18 canonical SNPs and 11 core variable-number tandem-repeat (VNTR) markers, and detection of rifampicin, dapsone and fluoroquinolone resistance-associated mutations in rpoB/ctpC/ctpI, folP1, gyrA/gyrB, respectively, and hypermutation-associated mutations in nth.

METHODS

The limit of detection (LOD) was determined using DNA of M. leprae reference strains and from 246 skin biopsies and 74 slit skin smears of leprosy patients, with genome copies quantified by RLEP qPCR. Sequencing results were evaluated versus whole genome sequencing (WGS) data of 14 strains, and versus VNTR-fragment length analysis (FLA) results of 89 clinical specimens.

FINDINGS

The LOD for sequencing success ranged between 80 and 3000 genome copies, depending on the sample type. The LOD for minority variants was 10%. All SNPs detected in targets by WGS were identified except in a clinical sample where WGS revealed two dapsone resistance-conferring mutations instead of one by Deeplex Myc-Lep, due to partial duplication of the sulfamide-binding domain in folP1. SNPs detected uniquely by Deeplex Myc-Lep were missed by WGS due to insufficient coverage. Concordance with VNTR-FLA results was 99.4% (926/932 alleles).

INTERPRETATION

Deeplex Myc-Lep may help improve the diagnosis and surveillance of leprosy. Gene domain duplication is an original putative drug resistance-related genetic adaptation in M. leprae.

FUNDING

EDCTP2 programme supported by the European Union (grant number RIA2017NIM-1847 -PEOPLE). EDCTP, R2Stop: Effect:Hope, The Mission To End Leprosy, the Flemish Fonds Wetenschappelijk Onderzoek.

Improving the extraction of ancient Yersinia pestis genomes from the dental pulp

Ancient DNA preserved in the dental pulp offers the opportunity to characterize the genome of some of the deadliest pathogens in human history. However, while DNA capture technologies help, focus sequencing efforts, and therefore, reduce experimental costs, the recovery of ancient pathogen DNA remains challenging. Here, we tracked the kinetics of ancient Yersinia pestis DNA release in solution during a pre-digestion of the dental pulp. We found that most of the ancient Y. pestis DNA is released within 60 min at 37°C in our experimental conditions. We recommend a simple pre-digestion as an economical procedure to obtain extracts enriched in ancient pathogen DNA, as longer digestion times release other types of templates, including host DNA. Combining this procedure with DNA capture, we characterized the genome sequences of 12 ancient Y. pestis bacteria from France dating to the second pandemic outbreaks of the 17^th and 18^th centuries Common Era.

Protocol, rationale and design of BE-PEOPLE (Bedaquiline enhanced exposure prophylaxis for LEprosy in the Comoros): a cluster randomized trial on effectiveness of rifampicin and bedaquiline as post-exposure prophylaxis of leprosy contacts

BACKGROUND

Leprosy is an ancient infectious disease with an annual global incidence of around 200,000 over the past decade. Since 2018, the World Health Organization (WHO) recommends single-dose rifampicin as post-exposure prophylaxis (SDR-PEP) for contacts of leprosy patients. The Post ExpOsure Prophylaxis for Leprosy (PEOPLE) trial evaluated PEP with a double dose of rifampicin in Comoros and Madagascar. Preliminary results of this trial show some reduction in leprosy incidence in intervention villages but a stronger regimen may be beneficial. The objective of the current Bedaquiline Enhanced ExpOsure Prophylaxis for LEprosy trial (BE-PEOPLE) is to explore effectiveness of a combination of bedaquiline and rifampicin as PEP.

METHODS

BE-PEOPLE is a cluster-randomized trial in which 44 clusters in Comoros will be randomized to two study arms. Door-to-door screening will be conducted annually during four years, leprosy patients identified will be offered standard of care treatment. Based on study arm, contacts aged five years and above and living within a 100-meter radius of an index case will either receive bedaquiline (400-800 mg) and rifampicin (150-600 mg) or only rifampicin (150-600 mg). Contacts aged two to four years will receive rifampicin only. Household contacts randomized to the bedaquiline plus rifampicin arm will receive a second dose four weeks later. Incidence rate ratios of leprosy comparing contacts who received either of the PEP regimens will be the primary outcome. We will monitor resistance to rifampicin and/or bedaquiline through molecular surveillance in all incident tuberculosis and leprosy patients nationwide. At the end of the study, we will assess anti-M. leprae PGL-I IgM seropositivity as a proxy for the population burden of M. leprae infection in 8 villages (17,000 individuals) that were surveyed earlier as part of the PEOPLE trial.

DISCUSSION

The COLEP trial on PEP in Bangladesh documented a reduction of 57% in incidence of leprosy among contacts treated with SDR-PEP after two years, which led to the WHO recommendation of SDR-PEP. Preliminary results of the PEOPLE trial show a lesser reduction in incidence. The BE-PEOPLE trial will explore whether reinforcing SDR-PEP with bedaquiline increases effectiveness and more rapidly reduces the incidence of leprosy, compared to SDR-PEP alone.

TRIAL REGISTRATION

NCT05597280. Protocol version 5.0 on 28 October 2022.

Rickettsia felis DNA recovered from a child who lived in southern Africa 2000 years ago

The Stone Age record of South Africa provides some of the earliest evidence for the biological and cultural origins of Homo sapiens. While there is extensive genomic evidence for the selection of polymorphisms in response to pathogen-pressure in sub-Saharan Africa, e.g., the sickle cell trait which provides protection against malaria, there is inadequate direct human genomic evidence for ancient human-pathogen infection in the region. Here, we analysed shotgun metagenome libraries derived from the sequencing of a Later Stone Age hunter-gatherer child who lived near Ballito Bay, South Africa, c. 2000 years ago. This resulted in the identification of ancient DNA sequence reads homologous to Rickettsia felis, the causative agent of typhus-like flea-borne rickettsioses, and the reconstruction of an ancient R. felis genome.

Mycobacterium leprae in Armadillo Tissues from Museum Collections, United States

We examined armadillos from museum collections in the United States using molecular assays to detect leprosy-causing bacilli. We found Mycobacterium leprae bacilli in samples from the United States, Bolivia, and Paraguay; prevalence was 14.8% in nine-banded armadillos. US isolates belonged to subtype 3I-2, suggesting long-term circulation of this genotype.

Analysis of Ancient Microbial DNA

The development of next-generation sequencing has led to a breakthrough in the analysis of ancient genomes, and the subsequent genomic analyses of ancient human skeletal remains have revolutionized our understanding of human evolution. This research led to the discovery of a new hominin lineage, and demonstrated multiple admixture events with more distantly related archaic human populations such as Neandertals and Denisovans over the last 100,000 years. Moreover, it has also yielded novel insights into the evolution of ancient pathogens. The analysis of ancient microbial genomes enables the study of their recent evolution, presently covering the last several millennia. These spectacular results have been obtained despite the degradation of DNA that takes place after the death of the host and increases with time. This cumulative degradation results in very short ancient DNA molecules, low in quantity, and highly prone to contamination by modern DNA molecules, especially from human and animal DNA present in reagents used in downstream biomolecular analyses. Finally, the minute amounts of ancient molecules are further diluted in environmental DNA from the soil microorganisms that colonize bones and teeth. Thus, ancient skeletal remains can share DNA profiles with environmental samples, and the identification of ancient microbial genomes among the more recent, presently poorly characterized, environmental microbiome is particularly challenging. Here, we describe the methods developed and/or in use in our laboratory to produce reliable and reproducible paleogenomic results from ancient skeletal remains that can be used to identify the presence of ancient microbiota.

Pathogenicity and virulence of Mycobacterium leprae

Leprosy is caused by Mycobacterium leprae (M. leprae) and M. lepromatosis, an obligate intracellular organism, and over 200,000 new cases occur every year. M. leprae parasitizes histiocytes (skin macrophages) and Schwann cells in the peripheral nerves. Although leprosy can be treated by multidrug therapy, some patients relapse or have a prolonged clinical course and/or experience leprosy reaction. These varying outcomes depend on host factors such as immune responses against bacterial components that determine a range of symptoms. To understand these host responses, knowledge of the mechanisms by which M. leprae parasitizes host cells is important. This article describes the characteristics of leprosy through bacteriology, genetics, epidemiology, immunology, animal models, routes of infection, and clinical findings. It also discusses recent diagnostic methods, treatment, and measures according to the World Health Organization (WHO), including prevention. Recently, the antibacterial activities of anti-hyperlipidaemia agents against other pathogens, such as M. tuberculosis and Staphylococcus aureus have been investigated. Our laboratory has been focused on the metabolism of lipids which constitute the cell wall of M. leprae. Our findings may be useful for the development of future treatments.

Examining pathogen DNA recovery across the remains of a 14th century Italian friar (Blessed Sante) infected with Brucella melitensis

OBJECTIVE

To investigate variation in ancient DNA recovery of Brucella melitensis, the causative agent of brucellosis, from multiple tissues belonging to one individual MATERIALS: 14 samples were analyzed from the mummified remains of the Blessed Sante, a 14 ^th century Franciscan friar from central Italy, with macroscopic diagnosis of probable brucellosis.

METHODS

Shotgun sequencing data from was examined to determine the presence of Brucella DNA.

RESULTS

Three of the 14 samples contained authentic ancient DNA, identified as belonging to B. melitensis. A genome (23.81X depth coverage, 0.98 breadth coverage) was recovered from a kidney stone. Nine of the samples contained reads classified as B. melitensis (7-169), but for many the data quality was insufficient to withstand our identification and authentication criteria.

CONCLUSIONS

We identified significant variation in the preservation and abundance of B. melitensis DNA present across multiple tissues, with calcified nodules yielding the highest number of authenticated reads. This shows how greatly sample selection can impact pathogen identification.

SIGNIFICANCE

Our results demonstrate variation in the preservation and recovery of pathogen DNA across tissues. This study highlights the importance of sample selection in the reconstruction of infectious disease burden and highlights the importance of a holistic approach to identifying disease.

LIMITATIONS

Study focuses on pathogen recovery in a single individual.

SUGGESTIONS FOR FURTHER RESEARCH

Further analysis of how sampling impacts aDNA recovery will improve pathogen aDNA recovery and advance our understanding of disease in past peoples.

Historic and Prehistoric Epidemics: An Overview of Sources Available for the Study of Ancient Pathogens

Since life on earth developed, parasitic microbes have thrived. Increases in host numbers, or the conquest of a new species, provide an opportunity for such a pathogen to enjoy, before host defense systems kick in, a similar upsurge in reproduction. Outbreaks, caused by "endemic" pathogens, and epidemics, caused by "novel" pathogens, have thus been creating chaos and destruction since prehistorical times. To study such (pre)historic epidemics, recent advances in the ancient DNA field, applied to both archeological and historical remains, have helped tremendously to elucidate the evolutionary trajectory of pathogens. These studies have offered new and unexpected insights into the evolution of, for instance, smallpox virus, hepatitis B virus, and the plague-causing bacterium Yersinia pestis. Furthermore, burial patterns and historical publications can help in tracking down ancient pathogens. Another source of information is our genome, where selective sweeps in immune-related genes relate to past pathogen attacks, while multiple viruses have left their genomes behind for us to study. This review will discuss the sources available to investigate (pre)historic diseases, as molecular knowledge of historic and prehistoric pathogens may help us understand the past and the present, and prepare us for future epidemics.

Metagenomic Research of Infectious Diseases in Archaeological Contexts: Evidence from the Hospital Real de Todos-os-Santos (Portugal)

Syphilis is one of the most exciting diseases explored in paleopathology and, therefore, tracing back its origin and development has provided a prolific debate. The combination of paleopathological data with historical sources, iconography, and archaeological contexts were the primary sources used to reconstruct its historical path. However, there are some limitations to paleopathological diagnosis due to the nature of bone reaction to stimuli. In addition, historical sources are subjected to a bias of social and cultural nature and the knowledge of those who wrote them. Hence, ancient DNA analysis offers the possibility of acquiring proof of cause by identifying pathogens in an organism. We undertook a metagenomic study of a skeleton exhumed from the Royal Hospital of All Saints (Portugal), renowned for treating syphilis from the 16th century onwards. The skeleton had previously been diagnosed with syphilis according to paleopathological analysis. However, the metagenomics analysis showed no presence of the pathogen associated with syphilis (i.e., Treponema pallidum) but revealed pathogenic microorganisms related to respiratory diseases (pneumonia), nonspecific bone infections (osteomyelitis), and oral bacterial pathologies as well as Hansen’s disease (also known as leprosy). The results are exciting and demand a reappraisal of the observed bone changes, recontextualizing their characterization as syphilis related. They prove that past reconstruction of health and disease diagnoses based on assessing human osteological remains of known context (such as a syphilitic hospital) may bias interpretations and, therefore, caution is recommended, not forgetting that the absence of evidence is not evidence of absence (in this case of syphilis) in life.

Creating communities of care: Sex estimation and mobility histories of adolescents buried in the cemetery of St. Mary Magdalen leprosarium (Winchester, England)

Objectives

This study examines the biological sex and geographical origins of adolescents buried at the St Mary Magdalen leprosarium (Winchester, UK). The data are combined with archaeological and palaeopathological evidence to broaden the understanding of mobility and its relationship to leprosy and leprosaria in Medieval England.

Materials and Methods

Nineteen individuals (~10–25 at death) with skeletal lesions diagnostic of leprosy were analyzed using standard osteological methods. Amelogenin peptides were extracted from five individuals whose biological sex could not be assessed from macroscopic methods. Enamel samples were analyzed to produce ⁸⁷Sr/⁸⁶Sr and δ¹⁸O values to explore mobility histories.

Results

Amelogenin peptides revealed three males and two females. Tooth enamel samples provided an ⁸⁷Sr/⁸⁶Sr ratio range from 0.7084 to 0.7103 (mean 0.7090, ±0.0012, 2σ). δ¹⁸O_P values show a wide range of 15.6‰–19.3‰ (mean 17.8 ± 1.6‰ 2σ), with corresponding δ¹⁸O_DW values ranging from −9.7‰ to −4.1‰ (mean −6.3 ± 2.4‰ 2σ).

Discussion

Amelogenin peptide data reveal the presence of adolescent females with bone changes of leprosy, making them the youngest confirmed females with leprosy in the archaeological record. Results also show at least 12 adolescents were local, and seven were from further afield, including outside Britain. Since St. Mary Magdalen was a leprosarium, it is possible that these people traveled there specifically for care. Archaeological and palaeopathological data support the notion that care was provided at this facility and that leprosy stigma, as we understand it today, may not have existed in this time and place.

Stone Age Yersinia pestis genomes shed light on the early evolution, diversity, and ecology of plague

The bacterium Yersinia pestis has caused numerous historically documented outbreaks of plague and research using ancient DNA could demonstrate that it already affected human populations during the Neolithic. However, the pathogen’s genetic diversity, geographic spread, and transmission dynamics during this early period of Y. pestis evolution are largely unexplored. Here, we describe a set of ancient plague genomes up to 5,000 y old from across Eurasia. Our data demonstrate that two genetically distinct forms of Y. pestis evolved in parallel and were both distributed across vast geographic distances, potentially occupying different ecological niches. Interpreted within the archeological context, our results suggest that the spread of plague during this period was linked to increased human mobility and intensification of animal husbandry.

The bacterial pathogen Yersinia pestis gave rise to devastating outbreaks throughout human history, and ancient DNA evidence has shown it afflicted human populations as far back as the Neolithic. Y. pestis genomes recovered from the Eurasian Late Neolithic/Early Bronze Age (LNBA) period have uncovered key evolutionary steps that led to its emergence from a Yersinia pseudotuberculosis-like progenitor; however, the number of reconstructed LNBA genomes are too few to explore its diversity during this critical period of development. Here, we present 17 Y. pestis genomes dating to 5,000 to 2,500 y BP from a wide geographic expanse across Eurasia. This increased dataset enabled us to explore correlations between temporal, geographical, and genetic distance. Our results suggest a nonflea-adapted and potentially extinct single lineage that persisted over millennia without significant parallel diversification, accompanied by rapid dispersal across continents throughout this period, a trend not observed in other pathogens for which ancient genomes are available. A stepwise pattern of gene loss provides further clues on its early evolution and potential adaptation. We also discover the presence of the flea-adapted form of Y. pestis in Bronze Age Iberia, previously only identified in in the Caucasus and the Volga regions, suggesting a much wider geographic spread of this form of Y. pestis. Together, these data reveal the dynamic nature of plague’s formative years in terms of its early evolution and ecology.

Geographically dispersed zoonotic tuberculosis in pre-contact South American human populations

Previous ancient DNA research has shown that Mycobacterium pinnipedii, which today causes tuberculosis (TB) primarily in pinnipeds, infected human populations living in the coastal areas of Peru prior to European colonization. Skeletal evidence indicates the presence of TB in several pre-colonial South and North American populations with minimal access to marine resources— a scenario incompatible with TB transmission directly from infected pinnipeds or their tissues. In this study, we investigate the causative agent of TB in ten pre-colonial, non-coastal individuals from South America. We reconstruct M. pinnipedii genomes (10- to 15-fold mean coverage) from three contemporaneous individuals from inland Peru and Colombia, demonstrating the widespread dissemination of M. pinnipedii beyond the coast, either through human-to-human and/or animal-mediated routes. Overall, our study suggests that TB transmission in the pre-colonial era Americas involved a more complex transmission pathway than simple pinniped-to-human transfer.

A probable case of leprosy from colonial period St. Vincent and the Grenadines, Southeastern Caribbean

OBJECTIVE

To document and differentially diagnose facial pathology found in an isolated skull from St. Vincent and the Grenadines, southeastern Caribbean. To directly date this individual using radiocarbon dating.

MATERIALS

Isolated skull recovered from Petite Mustique Island.

METHODS

Describe facial pathology occurring in this individual and compare with known diseases or disease processes that impact the craniofacial complex.

RESULTS

Features of the rhinomaxillary syndrome are present, indicating a diagnosis of leprosy. Dating places the time of death to the late 18th or early 19th centuries.

CONCLUSIONS

Analysis of the rhinomaxillary syndrome produces a diagnosis of early-stage leprosy in an individual that correlates with the apparent attempt to locate a leprosarium on Petite Mustique Island in the first decade of the 19th century.

SIGNIFICANCE

Location and time corroborate historical records of at least one attempt to locate a leprosarium on Petite Mustique Island. Only directly dated individual with leprosy in the western hemisphere and possibly the earliest yet recorded.

LIMITATIONS

This is an isolated find that is archaeologically unprovenienced.

SUGGESTIONS FOR FURTHER RESEARCH

Professional archaeological survey of Petite Mustique.

Genomics of Ancient Pathogens: First Advances and Prospects

Paleogenomics is one of the urgent and promising areas of interdisciplinary research in the today’s world science. New genomic methods of ancient DNA (aDNA) analysis, such as next generation sequencing (NGS) technologies, make it possible not only to obtain detailed genetic information about historical and prehistoric human populations, but also to study individual microbial and viral pathogens and microbiomes from different ancient and historical objects. Studies of aDNA of pathogens by reconstructing their genomes have so far yielded complete sequences of the ancient pathogens that played significant role in the history of the world: Yersiniapestis (plague), Variola virus (smallpox), Vibriocholerae (cholera), HBV (hepatitis B virus), as well as the equally important endemic human infectious agents: Mycobacteriumtuberculosis (tuberculosis), Mycobacteriumleprae (leprosy), and Treponemapallidum (syphilis). Genomic data from these pathogens complemented the information previously obtained by paleopathologists and allowed not only to identify pathogens from the past pandemics, but also to recognize the pathogen lineages that are now extinct, to refine chronology of the pathogen appearance in human populations, and to reconstruct evolutionary history of the pathogens that are still relevant to public health today. In this review, we describe state-of-the-art genomic research of the origins and evolution of many ancient pathogens and viruses and examine mechanisms of the emergence and spread of the ancient infections in the mankind history.

A new perspective on ancient Mitis group streptococcal genetics

Mitis group Streptococcus are human obligate bacteria residing in the nasopharynx and oral cavity. They comprise both commensal and pathogenic species with the most well-known being Streptococcus pneumoniae – a leading cause of meningitis and pneumonia. A primary difference between the commensal and pathogenic species is the presence of the polysaccharide capsule – a major virulence factor in S. pneumoniae, also present in other commensal species. Our current understanding of the evolutionary divergence of the pathogenic and commensal species has been inferred from extant strains. Ancient genomes can further elucidate streptococcal evolutionary history. We extracted streptococcal genome reads from a 5700-year-old ancient metagenome and worked towards characterizing them. Due to excessive within- and between-species recombination common among streptococci we were unable to parse individual species. Further, the composite reads of the ancient metagenome do not fit within the diversity of any specific extant species. Using a capsular gene database and AT-content analysis we determined that this ancient metagenome is missing polysaccharide synthesis genes integral to streptococcal capsule formation. The presence of multiple zinc metalloproteases suggests that adaptation to host IgA1 had begun and the presence of other virulence factors further implies development of close host–microbe interactions, though the absence of a capsule suggests an inability to cause invasive disease. The presence of specific virulence factors such as pneumolysin implies stable maintenance of such genes through streptococcal evolution that may strengthen their value as anti-pneumococcal vaccine antigens, while maintaining awareness of their potential presence in commensal species. Following from Jensen et al.’s initial analysis we provide historical context for this long time human nasopharyngeal resident, the Mitis group Streptococcus.

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-DR^highFBP1^high 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.

Ancient Metagenomic Studies: Considerations for the Wider Scientific Community

Like modern metagenomics, ancient metagenomics is a highly data-rich discipline, with the added challenge that the DNA of interest is degraded and, depending on the sample type, in low abundance. This requires the application of specialized measures during molecular experiments and computational analyses. Furthermore, researchers often work with finite sample sizes, which impedes optimal experimental design and control of confounding factors, and with ethically sensitive samples necessitating the consideration of additional guidelines. In September 2020, early career researchers in the field of ancient metagenomics met (Standards, Precautions & Advances in Ancient Metagenomics 2 [SPAAM2] community meeting) to discuss the state of the field and how to address current challenges. Here, in an effort to bridge the gap between ancient and modern metagenomics, we highlight and reflect upon some common misconceptions, provide a brief overview of the challenges in our field, and point toward useful resources for potential reviewers and newcomers to the field.

Human macrophage polarization in the response to Mycobacterium leprae genomic DNA

Infection with Mycobacterium leprae, the causative organism of leprosy, is still endemic in numerous parts of the world including the southwestern United States. The broad variation of symptoms in the leprosy disease spectrum range from the milder tuberculoid leprosy (paucibacillary) to the more severe and disfiguring lepromatous leprosy (multibacillary). The established thinking in the health community is that host response, rather than M. leprae strain variation, is the reason for the range of disease severity. More recent discoveries suggest that macrophage polarization also plays a significant role in the spectrum of leprosy disease but to what degree it contributes is not fully established. In this study, we aimed to analyze if different strains of M. leprae elicit different transcription responses in human macrophages, and to examine the role of macrophage polarization in these responses. Genomic DNA from three different strains of M. leprae DNA (Strains NHDP, Br4923, and Thai-53) were used to stimulate human macrophages under three polarization conditions (M1, M1-activated, and M2). Transcriptome analysis revealed a large number of differentially expressed (DE) genes upon stimulation with DNA from M. leprae strain Thai-53 compared to strains NHDP and Br4923, independent of the macrophage polarization condition. We also found that macrophage polarization affects the responses to M. leprae DNA, with up-regulation of numerous interferon stimulated genes. These findings provide a deeper understanding of the role of macrophage polarization in the recognition of M. leprae DNA, with the potential to improve leprosy treatment strategies.

One Health Approaches to Trace Mycobacterium leprae's Zoonotic Potential Through Time

Hansen's disease (leprosy), mainly caused by infection with Mycobacterium leprae, has accompanied humanity for thousands of years. Although currently rare in Europe, there are over 200,000 new infections annually in South East Asia, Africa, and South America. Over the years many disciplines - palaeopathology, ancient DNA and other ancient biomolecules, and history - have contributed to a better understanding of leprosy's past, in particular its history in medieval Europe. We discuss their contributions and potential, especially in relation to the role of inter-species transmission, an unexplored phenomenon in the disease's history. Here, we explore the potential of interdisciplinary approaches that understand disease as a biosocial phenomenon, which is a product of both infection with M. leprae and social behaviours that facilitate transmission and spread. Genetic evidence of M. leprae isolated from archaeological remains combined with systematic zooarchaeological and historical analysis would not only identify when and in what direction transmission occurred, but also key social behaviours and motivations that brought species together. In our opinion, this combination is crucial to understand the disease's zoonotic past and current potential.

Molecular archeology of human viruses

The evolution of human-virus associations is usually reconstructed from contemporary patterns of genomic diversity. An intriguing, though still rarely implemented, alternative is to search for the genetic material of viruses in archeological and medical archive specimens to document evolution as it happened. In this chapter, we present lessons from ancient DNA research and incorporate insights from virology to explore the potential range of applications and likely limitations of archeovirological approaches. We also highlight the numerous questions archeovirology will hopefully allow us to tackle in the near future, and the main expected roadblocks to these avenues of research.

Mycobacterium leprae diversity and population dynamics in medieval Europe from novel ancient genomes

BACKGROUND

Hansen's disease (leprosy), widespread in medieval Europe, is today mainly prevalent in tropical and subtropical regions with around 200,000 new cases reported annually. Despite its long history and appearance in historical records, its origins and past dissemination patterns are still widely unknown. Applying ancient DNA approaches to its major causative agent, Mycobacterium leprae, can significantly improve our understanding of the disease's complex history. Previous studies have identified a high genetic continuity of the pathogen over the last 1500 years and the existence of at least four M. leprae lineages in some parts of Europe since the Early Medieval period.

RESULTS

Here, we reconstructed 19 ancient M. leprae genomes to further investigate M. leprae's genetic variation in Europe, with a dedicated focus on bacterial genomes from previously unstudied regions (Belarus, Iberia, Russia, Scotland), from multiple sites in a single region (Cambridgeshire, England), and from two Iberian leprosaria. Overall, our data confirm the existence of similar phylogeographic patterns across Europe, including high diversity in leprosaria. Further, we identified a new genotype in Belarus. By doubling the number of complete ancient M. leprae genomes, our results improve our knowledge of the past phylogeography of M. leprae and reveal a particularly high M. leprae diversity in European medieval leprosaria.

CONCLUSIONS

Our findings allow us to detect similar patterns of strain diversity across Europe with branch 3 as the most common branch and the leprosaria as centers for high diversity. The higher resolution of our phylogeny tree also refined our understanding of the interspecies transfer between red squirrels and humans pointing to a late antique/early medieval transmission. Furthermore, with our new estimates on the past population diversity of M. leprae, we gained first insights into the disease's global history in relation to major historic events such as the Roman expansion or the beginning of the regular transatlantic long distance trade. In summary, our findings highlight how studying ancient M. leprae genomes worldwide improves our understanding of leprosy's global history and can contribute to current models of M. leprae's worldwide dissemination, including interspecies transmissions.

Salmonella enterica from a soldier from the 1652 siege of Barcelona (Spain) supports historical transatlantic epidemic contacts

Ancient pathogen genomics is an emerging field allowing reconstruction of past epidemics. The demise of post-contact American populations may, at least in part, have been caused by paratyphoid fever brought by Europeans. We retrieved genome-wide data from two Spanish soldiers who were besieging the city of Barcelona in 1652, during the Reapers' War. Their ancestry derived from the Basque region and Sardinia, respectively, (at that time, this island belonged to the Spanish kingdom). Despite the proposed plague epidemic, we could not find solid evidence for the presence of the causative plague agent in these individuals. However, we retrieved from one individual a substantial fraction of the Salmonella enterica serovar Paratyphi C lineage linked to paratyphoid fever in colonial period Mexico. Our results support a growing body of evidence that Paratyphi C enteric fever was more prevalent in Europe and the Americas in the past than it is today.

Simultaneous Determination of Mycobacterium leprae Drug Resistance and Single-Nucleotide Polymorphism Genotype by Use of Nested Multiplex PCR with Amplicon Sequencing

Mycobacterium leprae is the predominant cause of leprosy worldwide, and its genotypes can be classified into four single-nucleotide polymorphism (SNP) types and 16 subtypes. Determining M. leprae drug resistance and genotype is typically done by PCR and Sanger DNA sequencing, which require substantial effort. Here, we describe a rapid method involving multiplex PCR in combination with nested amplification and next-generation sequence analysis that allows simultaneous determination of M. leprae drug resistance and SNP genotype directly from clinical specimens. We used this method to analyze clinical samples from two paucibacillary, nine multibacillary, and six type-undetermined leprosy patients. Regions in folP1, rpoB, gyrA, and gyrB that determine drug resistance and those for 84 SNP-InDels in the M. leprae genome were amplified from clinical samples and their sequences determined. The results showed that seven samples were subtype 1A, three were 1D, and seven were 3K. Three samples of the subtype 3K had folp1 mutation. The method may allow more rapid genetic analyses of M. leprae in clinical samples.

The Recovery, Interpretation and Use of Ancient Pathogen Genomes

The ability to sequence genomes from ancient biological material has provided a rich source of information for evolutionary biology and engaged considerable public interest. Although most studies of ancient genomes have focused on vertebrates, particularly archaic humans, newer technologies allow the capture of microbial pathogens and microbiomes from ancient and historical human and non-human remains. This coming of age has been made possible by techniques that allow the preferential capture and amplification of discrete genomes from a background of predominantly host and environmental DNA. There are now near-complete ancient genome sequences for three pathogens of considerable historical interest — pre-modern bubonic plague (Yersinia pestis), smallpox (Variola virus) and cholera (Vibrio cholerae) — and for three equally important endemic human disease agents — Mycobacterium tuberculosis (tuberculosis), Mycobacterium leprae (leprosy) and Treponema pallidum pallidum (syphilis). Genomic data from these pathogens have extended earlier work by paleopathologists. There have been efforts to sequence the genomes of additional ancient pathogens, with the potential to broaden our understanding of the infectious disease burden common to past populations from the Bronze Age to the early 20^th century. In this review we describe the state-of-the-art of this rapidly developing field, highlight the contributions of ancient pathogen genomics to multidisciplinary endeavors and describe some of the limitations in resolving questions about the emergence and long-term evolution of pathogens.

A Small Protein but with Diverse Roles: A Review of EsxA in Mycobacterium-Host Interaction

As a major effector of the ESX-1 secretion system, EsxA is essential for the virulence of pathogenic mycobacteria, such as Mycobacterium tuberculosis (Mtb) and Mycobacterium marinum (Mm). EsxA possesses an acidic pH-dependent membrane permeabilizing activity and plays an essential role by mediating mycobacterial escape from the phagosome and translocation to the cytosol for intracellular replication. Moreover, EsxA regulates host immune responses as a potent T-cell antigen and a strong immunoregulator. EsxA interacts with multiple cellular proteins and stimulates several signal pathways, such as necrosis, apoptosis, autophagy, and antigen presentation. Interestingly, there is a co-dependency in the expression and secretion of EsxA and other mycobacterial factors, which greatly increases the complexity of dissecting the precise roles of EsxA and other factors in mycobacterium-host interaction. In this review, we summarize the current understandings of the roles and functions of EsxA in mycobacterial infection and discuss the challenges and future directions.

Patterns of transmission and horizontal gene transfer in the Dioscorea sansibarensis leaf symbiosis revealed by whole-genome sequencing

Leaves of the wild yam species Dioscorea sansibarensis display prominent forerunner or "drip" tips filled with extracellular bacteria of the species Orrella dioscoreae.¹ This species of yam is native to Madagascar and tropical Africa and reproduces mainly asexually through aerial bulbils and underground tubers, which also contain a small population of O. dioscoreae.^2,3 Despite apparent vertical transmission, the genome of O. dioscoreae does not show any of the hallmarks of genome erosion often found in hereditary symbionts (e.g., small genome size and accumulation of pseudogenes).^4-6 We investigated here the range and distribution of leaf symbiosis between D. sansibarensis and O. dioscoreae using preserved leaf samples from herbarium collections that were originally collected from various locations in Africa. We recovered DNA from the extracellular symbiont in all samples, showing that the symbiosis is widespread throughout continental Africa and Madagascar. Despite the degraded nature of this DNA, we constructed 17 symbiont genomes using de novo methods without relying on a reference. Phylogenetic and genomic analyses revealed that horizontal transmission of symbionts and horizontal gene transfer have shaped the evolution of the symbiont. These mechanisms could help explain lack of signs of reductive genome evolution despite an obligate host-associated lifestyle. Furthermore, phylogenetic analysis of D. sansibarensis based on plastid genomes revealed a strong geographical clustering of samples and provided evidence that the symbiosis originated at least 13 mya, earlier than previously estimated.³.

Additive Uncorrelated Relaxed Clock Models for the Dating of Genomic Epidemiology Phylogenies

Phylogenetic dating is one of the most powerful and commonly used methods of drawing epidemiological interpretations from pathogen genomic data. Building such trees requires considering a molecular clock model which represents the rate at which substitutions accumulate on genomes. When the molecular clock rate is constant throughout the tree then the clock is said to be strict, but this is often not an acceptable assumption. Alternatively, relaxed clock models consider variations in the clock rate, often based on a distribution of rates for each branch. However, we show here that the distributions of rates across branches in commonly used relaxed clock models are incompatible with the biological expectation that the sum of the numbers of substitutions on two neighboring branches should be distributed as the substitution number on a single branch of equivalent length. We call this expectation the additivity property. We further show how assumptions of commonly used relaxed clock models can lead to estimates of evolutionary rates and dates with low precision and biased confidence intervals. We therefore propose a new additive relaxed clock model where the additivity property is satisfied. We illustrate the use of our new additive relaxed clock model on a range of simulated and real data sets, and we show that using this new model leads to more accurate estimates of mean evolutionary rates and ancestral dates.

Reconstruction of ancient microbial genomes from the human gut

Loss of gut microbial diversity^1–6 in industrial populations is associated with chronic diseases⁷, underscoring the importance of studying our ancestral gut microbiome. However, relatively little is known about the composition of pre-industrial gut microbiomes. Here we performed a large-scale de novo assembly of microbial genomes from palaeofaeces. From eight authenticated human palaeofaeces samples (1,000–2,000 years old) with well-preserved DNA from southwestern USA and Mexico, we reconstructed 498 medium- and high-quality microbial genomes. Among the 181 genomes with the strongest evidence of being ancient and of human gut origin, 39% represent previously undescribed species-level genome bins. Tip dating suggests an approximate diversification timeline for the key human symbiont Methanobrevibacter smithii. In comparison to 789 present-day human gut microbiome samples from eight countries, the palaeofaeces samples are more similar to non-industrialized than industrialized human gut microbiomes. Functional profiling of the palaeofaeces samples reveals a markedly lower abundance of antibiotic-resistance and mucin-degrading genes, as well as enrichment of mobile genetic elements relative to industrial gut microbiomes. This study facilitates the discovery and characterization of previously undescribed gut microorganisms from ancient microbiomes and the investigation of the evolutionary history of the human gut microbiota through genome reconstruction from palaeofaeces.

Ancient microbiomes from palaeofaeces are more similar to non-industrialized than industrialized human gut microbiomes regardless of geography, but 39% of their de novo reconstructed genomes represent previously undescribed microbial species.

Mass burial genomics reveals outbreak of enteric paratyphoid fever in the Late Medieval trade city Lübeck

Medieval Europe was repeatedly affected by outbreaks of infectious diseases, some of which reached epidemic proportions. A Late Medieval mass burial next to the Heiligen-Geist-Hospital in Lübeck (present-day Germany) contained the skeletal remains of more than 800 individuals who had presumably died from infectious disease. From 92 individuals, we screened the ancient DNA extracts for the presence of pathogens to determine the cause of death. Metagenomic analysis revealed evidence of Salmonella enterica subsp. enterica serovar Paratyphi C, suggesting an outbreak of enteric paratyphoid fever. Three reconstructed S. Paratyphi C genomes showed close similarity to a strain from Norway (1200 CE). Radiocarbon dates placed the disease outbreak in Lübeck between 1270 and 1400 cal CE, with historical records indicating 1367 CE as the most probable year. The deceased were of northern and eastern European descent, confirming Lübeck as an important trading center of the Hanseatic League in the Baltic region.

•

Salmonella enterica Paratyphi C detected in remains from a mass burial in Lübeck

•

Outbreak of enteric paratyphoid fever likely occurred in 1367 CE

•

Pathogen genomes showed close similarity to a strain from Norway (1200 CE)

Development of a novel loop-mediated isothermal amplification assay for rapid detection of Mycobacterium leprae in clinical samples

BACKGROUND

Sensitive and definitive diagnostic tests are required for timely treatment of leprosy and to control its transmission.

AIM

In the present study, we report the development of loop-mediated isothermal amplification assay using six primers targeting the RLEP gene sequence uniquely present in Mycobacterium leprae.

METHODS

Tissue punch samples (n = 50) and slit aspirates (n = 50) from confirmed cases of leprosy (M. leprae positive by quantitative polymerase chain reaction), reporting at the Department of Dermatology, Safdarjung Hospital, New Delhi, were analyzed using newly developed closed tube loop-mediated isothermal amplification assay. The sensitivity and specificity; positive predictive value, negative predictive value and accuracy were calculated using MedCalc statistical software.

RESULTS

The loop-mediated isothermal amplification assay specifically amplified M. leprae genomic DNA with an analytical sensitivity of 100 fg. About 47 Out of the 50 quantitative polymerase chain reactions confirmed M. leprae positive tissue samples, 47 were positive by loop-mediated isothermal amplification assay (sensitivity 94%; 95% confidence interval 83.5%-98.8%) while only 31/50 were positive by histopathology (sensitivity 62%; 95% confidence interval 47.2%-75.4%) . Using slit aspirate samples of these 50 patients, 42 were positive by both quantitative polymerase chain reaction and loop-mediated isothermal amplification assay (sensitivity 84%; 95% confidence interval 70.9%-92.8%) while only 23/50 (sensitivity 46%; 95% confidence interval 31.8%-60.7%) were positive by microscopy.

LIMITATIONS

In the present study, the leprosy patient cohort was not uniform, as it comprised a lower number of paucibacillary cases (22%) compared to multibacillary (78%) cases.

CONCLUSION

Loop-mediated isothermal amplification assay established here provides a rapid and accurate diagnostic test for leprosy in terms of sensitivity and specificity. The assay is simple to perform in comparison with other molecular techniques (polymerase chain reaction/quantitative polymerase chain reaction) and has potential for field applicability.

Culture-free genome-wide locus sequence typing (GLST) provides new perspectives on Trypanosoma cruzi dispersal and infection complexity

Analysis of genetic polymorphism is a powerful tool for epidemiological surveillance and research. Powerful inference from pathogen genetic variation, however, is often restrained by limited access to representative target DNA, especially in the study of obligate parasitic species for which ex vivo culture is resource-intensive or bias-prone. Modern sequence capture methods enable pathogen genetic variation to be analyzed directly from host/vector material but are often too complex and expensive for resource-poor settings where infectious diseases prevail. This study proposes a simple, cost-effective 'genome-wide locus sequence typing' (GLST) tool based on massive parallel amplification of information hotspots throughout the target pathogen genome. The multiplexed polymerase chain reaction amplifies hundreds of different, user-defined genetic targets in a single reaction tube, and subsequent agarose gel-based clean-up and barcoding completes library preparation at under 4 USD per sample. Our study generates a flexible GLST primer panel design workflow for Trypanosoma cruzi, the parasitic agent of Chagas disease. We successfully apply our 203-target GLST panel to direct, culture-free metagenomic extracts from triatomine vectors containing a minimum of 3.69 pg/μl T. cruzi DNA and further elaborate on method performance by sequencing GLST libraries from T. cruzi reference clones representing discrete typing units (DTUs) TcI, TcIII, TcIV, TcV and TcVI. The 780 SNP sites we identify in the sample set repeatably distinguish parasites infecting sympatric vectors and detect correlations between genetic and geographic distances at regional (< 150 km) as well as continental scales. The markers also clearly separate TcI, TcIII, TcIV and TcV + TcVI and appear to distinguish multiclonal infections within TcI. We discuss the advantages, limitations and prospects of our method across a spectrum of epidemiological research.

An investigation of irreproducibility in maximum likelihood phylogenetic inference

Phylogenetic trees are essential for studying biology, but their reproducibility under identical parameter settings remains unexplored. Here, we find that 3515 (18.11%) IQ-TREE-inferred and 1813 (9.34%) RAxML-NG-inferred maximum likelihood (ML) gene trees are topologically irreproducible when executing two replicates (Run1 and Run2) for each of 19,414 gene alignments in 15 animal, plant, and fungal phylogenomic datasets. Notably, coalescent-based ASTRAL species phylogenies inferred from Run1 and Run2 sets of individual gene trees are topologically irreproducible for 9/15 phylogenomic datasets, whereas concatenation-based phylogenies inferred twice from the same supermatrix are reproducible. Our simulations further show that irreproducible phylogenies are more likely to be incorrect than reproducible phylogenies. These results suggest that a considerable fraction of single-gene ML trees may be irreproducible. Increasing reproducibility in ML inference will benefit from providing analyses’ log files, which contain typically reported parameters (e.g., program, substitution model, number of tree searches) but also typically unreported ones (e.g., random starting seed number, number of threads, processor type).

Replicate runs of maximum likelihood phylogenetic analyses can generate different tree topologies due to differences in parameters, such as random seeds. Here, Shen et al. demonstrate that replicate runs can generate substantially different tree topologies even with identical data and parameters.

Variola virus genome sequenced from an eighteenth-century museum specimen supports the recent origin of smallpox

Smallpox, caused by the variola virus (VARV), was a highly virulent disease with high mortality rates causing a major threat for global human health until its successful eradication in 1980. Despite previously published historic and modern VARV genomes, its past dissemination and diversity remain debated. To understand the evolutionary history of VARV with respect to historic and modern VARV genetic variation in Europe, we sequenced a VARV genome from a well-described eighteenth-century case from England (specimen P328). In our phylogenetic analysis, the new genome falls between the modern strains and another historic strain from Lithuania, supporting previous claims of larger diversity in early modern Europe compared to the twentieth century. Our analyses also resolve a previous controversy regarding the common ancestor between modern and historic strains by confirming a later date around the seventeenth century. Overall, our results point to the benefit of historic genomes for better resolution of past VARV diversity and highlight the value of such historic genomes from around the world to further understand the evolutionary history of smallpox as well as related diseases. This article is part of the theme issue 'Insights into health and disease from ancient biomolecules'.

Comparison of target enrichment strategies for ancient pathogen DNA

In ancient DNA research, the degraded nature of the samples generally results in poor yields of highly fragmented DNA; targeted DNA enrichment is thus required to maximize research outcomes. The three commonly used methods - array-based hybridization capture and in-solution capture using either RNA or DNA baits - have different characteristics that may influence the capture efficiency, specificity and reproducibility. Here we compare their performance in enriching pathogen DNA of Mycobacterium leprae and Treponema pallidum from 11 ancient and 19 modern samples. We find that in-solution approaches are the most effective method in ancient and modern samples of both pathogens and that RNA baits usually perform better than DNA baits.

A systematic investigation of human DNA preservation in medieval skeletons

Ancient DNA (aDNA) analyses necessitate the destructive sampling of archaeological material. Currently, the cochlea, part of the osseous inner ear located inside the petrous pyramid, is the most sought after skeletal element for molecular analyses of ancient humans as it has been shown to yield high amounts of endogenous DNA. However, destructive sampling of the petrous pyramid may not always be possible, particularly in cases where preservation of skeletal morphology is of top priority. To investigate alternatives, we present a survey of human aDNA preservation for each of ten skeletal elements in a skeletal collection from Medieval Germany. Through comparison of human DNA content and quality we confirm best performance of the petrous pyramid and identify seven additional sampling locations across four skeletal elements that yield adequate aDNA for most applications in human palaeogenetics. Our study provides a better perspective on DNA preservation across the human skeleton and takes a further step toward the more responsible use of ancient materials in human aDNA studies.

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'.