Eighteenth-century genomes show that mixed infections were common at time of peak tuberculosis in Europe

Integrative and comparative genomic analyses of mammalian macrophage responses to intracellular mycobacterial pathogens

Mycobacterium tuberculosis, the causative agent of human tuberculosis (hTB), is a close evolutionary relative of Mycobacterium bovis, which causes bovine tuberculosis (bTB), one of the most damaging infectious diseases to livestock agriculture. Previous studies have shown that the pathogenesis of bTB disease is comparable to hTB disease, and that the bovine and human alveolar macrophage (bAM and hAM, respectively) transcriptomes are extensively reprogrammed in response to infection with these intracellular mycobacterial pathogens. In this study, a multi-omics integrative approach was applied with functional genomics and GWAS data sets across the two primary hosts (Bos taurus and Homo sapiens) and both pathogens (M. bovis and M. tuberculosis). Four different experimental infection groups were used: 1) bAM infected with M. bovis, 2) bAM infected with M. tuberculosis, 3) hAM infected with M. tuberculosis, and 4) human monocyte-derived macrophages (hMDM) infected with M. tuberculosis. RNA-seq data from these experiments 24 h post-infection (24 hpi) was analysed using three computational pipelines: 1) differentially expressed genes, 2) differential gene expression interaction networks, and 3) combined pathway analysis. The results were integrated with high-resolution bovine and human GWAS data sets to detect novel quantitative trait loci (QTLs) for resistance to mycobacterial infection and resilience to disease. This revealed common and unique response macrophage pathways for both pathogens and identified 32 genes (12 bovine and 20 human) significantly enriched for SNPs associated with disease resistance, the majority of which encode key components of the NF-κB signalling pathway and that also drive formation of the granuloma.

Metagenomic and paleopathological analyses of a historic documented collection explore ancient dental calculus as a diagnostic tool

Dental calculus is a microbial biofilm that contains biomolecules from oral commensals and pathogens, including those potentially related to cause of death (CoD). To assess the utility of calculus as a diagnostically informative substrate, in conjunction with paleopathological analysis, calculus samples from 39 individuals in the Smithsonian Institution's Robert J. Terry Collection with CoDs of either syphilis or tuberculosis were assessed via shotgun metagenomic sequencing for the presence of Treponema pallidum subsp. pallidum and Mycobacterium tuberculosis complex (MTBC) DNA. Paleopathological analysis revealed that frequencies of skeletal lesions associated with these diseases were partially inconsistent with diagnostic criteria. Although recovery of T. p. pallidum DNA from individuals with a syphilis CoD was elusive, MTBC DNA was identified in at least one individual with a tuberculosis CoD. The authenticity of MTBC DNA was confirmed using targeted quantitative PCR assays, MTBC genome enrichment, and in silico bioinformatic analyses; however, the lineage of the MTBC strain present could not be determined. Overall, our study highlights the utility of dental calculus for molecular detection of tuberculosis in the archaeological record and underscores the effect of museum preparation techniques and extensive handling on pathogen DNA preservation in skeletal collections.

Exploring the Complexity of the Human Respiratory Virome through an In Silico Analysis of Shotgun Metagenomic Data Retrieved from Public Repositories

BACKGROUND

Respiratory viruses significantly impact global morbidity and mortality, causing more disease in humans than any other infectious agent. Beyond pathogens, various viruses and bacteria colonize the respiratory tract without causing disease, potentially influencing respiratory diseases' pathogenesis. Nevertheless, our understanding of respiratory microbiota is limited by technical constraints, predominantly focusing on bacteria and neglecting crucial populations like viruses. Despite recent efforts to improve our understanding of viral diversity in the human body, our knowledge of viral diversity associated with the human respiratory tract remains limited.

METHODS

Following a comprehensive search in bibliographic and sequencing data repositories using keyword terms, we retrieved shotgun metagenomic data from public repositories (n = 85). After manual curation, sequencing data files from 43 studies were analyzed using EVEREST (pipEline for Viral assEmbly and chaRactEriSaTion). Complete and high-quality contigs were further assessed for genomic and taxonomic characterization.

RESULTS

Viral contigs were obtained from 194 out of the 868 FASTQ files processed through EVEREST. Of the 1842 contigs that were quality assessed, 8% (n = 146) were classified as complete/high-quality genomes. Most of the identified viral contigs were taxonomically classified as bacteriophages, with taxonomic resolution ranging from the superkingdom level down to the species level. Captured contigs were spread across 25 putative families and varied between RNA and DNA viruses, including previously uncharacterized viral genomes. Of note, airway samples also contained virus(es) characteristic of the human gastrointestinal tract, which have not been previously described as part of the lung virome. Additionally, by performing a meta-analysis of the integrated datasets, ecological trends within viral populations linked to human disease states and their biogeographical distribution along the respiratory tract were observed.

CONCLUSION

By leveraging publicly available repositories of shotgun metagenomic data, the present study provides new insights into viral genomes associated with specimens from the human respiratory tract across different disease spectra. Further studies are required to validate our findings and evaluate the potential impact of these viral communities on respiratory tract physiology.

Evolution and emergence of Mycobacterium tuberculosis

Tuberculosis (TB) remains one of the deadliest infectious diseases in human history, prevailing even in the 21st century. The causative agents of TB are represented by a group of closely related bacteria belonging to the Mycobacterium tuberculosis complex (MTBC), which can be subdivided into several lineages of human- and animal-adapted strains, thought to have shared a last common ancestor emerged by clonal expansion from a pool of recombinogenic Mycobacterium canettii-like tubercle bacilli. A better understanding of how MTBC populations evolved from less virulent mycobacteria may allow for discovering improved TB control strategies and future epidemiologic trends. In this review, we highlight new insights into the evolution of mycobacteria at the genus level, describing different milestones in the evolution of mycobacteria, with a focus on the genomic events that have likely enabled the emergence and the dominance of the MTBC. We also review the recent literature describing the various MTBC lineages and highlight their particularities and differences with a focus on host preferences and geographic distribution. Finally, we discuss on putative mechanisms driving the evolution of tubercle bacilli and mycobacteria in general, by taking the mycobacteria-specific distributive conjugal transfer as an example.

Detection of a historic reservoir of bedaquiline/clofazimine resistance-associated variants in Mycobacterium tuberculosis

BACKGROUND

Drug resistance in tuberculosis (TB) poses a major ongoing challenge to public health. The recent inclusion of bedaquiline into TB drug regimens has improved treatment outcomes, but this advance is threatened by the emergence of strains of Mycobacterium tuberculosis (Mtb) resistant to bedaquiline. Clinical bedaquiline resistance is most frequently conferred by off-target resistance-associated variants (RAVs) in the mmpR5 gene (Rv0678), the regulator of an efflux pump, which can also confer cross-resistance to clofazimine, another TB drug.

METHODS

We compiled a dataset of 3682 Mtb genomes, including 180 carrying variants in mmpR5, and its immediate background (i.e. mmpR5 promoter and adjacent mmpL5 gene), that have been associated to borderline (henceforth intermediate) or confirmed resistance to bedaquiline. We characterised the occurrence of all nonsynonymous mutations in mmpR5 in this dataset and estimated, using time-resolved phylogenetic methods, the age of their emergence.

RESULTS

We identified eight cases where RAVs were present in the genomes of strains collected prior to the use of bedaquiline in TB treatment regimes. Phylogenetic reconstruction points to multiple emergence events and circulation of RAVs in mmpR5, some estimated to predate the introduction of bedaquiline. However, epistatic interactions can complicate bedaquiline drug-susceptibility prediction from genetic sequence data. Indeed, in one clade, Ile67fs (a RAV when considered in isolation) was estimated to have emerged prior to the antibiotic era, together with a resistance reverting mmpL5 mutation.

CONCLUSIONS

The presence of a pre-existing reservoir of Mtb strains carrying bedaquiline RAVs prior to its clinical use augments the need for rapid drug susceptibility testing and individualised regimen selection to safeguard the use of bedaquiline in TB care and control.

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.

Bioarchaeological and molecular evidence of tuberculosis in human skeletal remains from 18th-19th century orthodox cemeteries in Irkutsk, Eastern Siberia

In this study, we tested the skeletal human remains from the 18th - early 19th century Orthodox cemeteries in Irkutsk, Eastern Siberia, for tuberculosis-associated morphological alterations and Mycobacterium tuberculosis DNA. The morphologically studied bone collection included 591 individuals of mainly Caucasian origin. The molecular methods (IS6110-PCR and spoligotyping) suggested that at least four individuals (out of 15 TB-suspected, DNA-tested) were positive for the presence of M. tuberculosis DNA. All of them were males (3 maturus, 1 maturus senilis). Two of them date back to the second and third quarters of the 18th century, another to the last quarter of the 18th century, and the last one to the second half of the 19th century. The combined molecular analysis cautiously suggested presence of different strains and at least some of them represented not the currently predominant in Siberia Beijing genotype (M. tuberculosis East-Asian lineage) but strains of European origin. In conclusion, this study presented bioarchaeological and molecular evidence of tuberculosis in human skeletal remains from 18th-19th century Orthodox cemeteries in Irkutsk, Eastern Siberia. The samples are not M. bovis and represent human M. tuberculosis sensu stricto. Their precise phylogenetic identity is elusive but evokes the European/Russian origin of at least some isolates.

More than one millennium (2nd-16th century CE) of the White Plague in the Carpathian Basin - New cases, expanding knowledge

The causative agent of tuberculosis is still a widespread pathogen, which caused the death of ca. 1.6 million people globally in 2021. The paleopathological study of human remains revealed the antiquity of the disease and its continuous presence throughout the history of humankind. The Carpathian Basin has always been a biocultural melting pot, since it has seen several migrations over the centuries, and served as a location of admixture and interaction for numerous populations of different cultures. Thus, this geographical territory is ideal for the examination of the coevolutionary processes of hosts and their pathogens. We aimed to reveal the spatial and temporal distribution of tuberculosis cases excavated inside the borders of Hungary between the 2nd and 16th centuries CE. We established a comprehensive database by collecting 114 already published cases and introducing 39 new cases. The involved cases include those that have been confirmed by different molecular methods, as well as possible infections that were identified based on the presence of macromorphological and radiological alterations. The progress of future molecular and paleopathological studies can be facilitated by our dataset, as it presents spatial and temporal information concerning the spread of the disease in the Carpathian Basin, as well as the biological profile and detailed paleopathological description of lesions illustrated by photo- and radiographs.

Hit or miss - A metagenomic evaluation of intra-bone variability of host pathogen load in tuberculosis-infected human remains

Many sampling protocols have been established to successfully retrieve human DNA from archaeological remains, however the systematic detection of ancient pathogens remains challenging. Here, we present a first assessment of the intra-bone variability of metagenomic composition in human skeletal remains and its effect on the sampling success for Mycobacterium tuberculosis (MTB) and human endogenous DNA. For this purpose, four bone samples from published peer-reviewed studies with PCR-based evidence for ancient MTB DNA were selected. Two bone samples of a Neolithic individual from Halberstadt, Germany and two ribs of two 18th-century Hungarian church mummies were sampled at multiple locations for equal amounts, followed by DNA extraction and library construction. Shotgun sequencing data was generated for taxonomic profiling as well as quantitative and qualitative evaluation of MTB and human endogenous DNA. Despite low variance in microbial diversity within and across samples, intra-bone variability of up to 36.45- and 62.88-fold for authentic ancient MTB and human reads, respectively, was detected. This study demonstrates the variable sampling success for MTB and human endogenous DNA within single skeletal samples despite relatively consistent microbial composition and highlights how a multisampling approach can facilitate the detection of hotspots with highly concentrated pathogen and human endogenous DNA.

Tuberculosis in prehistory in Eastern Central Europe (Hungary) - Chronological and geographical distribution

Tuberculosis (TB) has long been a major scourge of humankind. Paleopathological and paleomicrobiological studies have revealed the past presence of the disease on a large spatial and temporal scale. The antiquity of the disease has extensively been studied in the Carpathian Basin, given its dynamic population and cultural changes since prehistory. These studies, however, have mainly focused on the populations living during the Common Era. The aim of this paper is to present the published and the recently discovered cases of prehistoric TB, from the Neolithic (6000-4500/4400 BCE) to the Bronze Age (2600/2500-800 BCE) Central Carpathian Basin (Hungary). We summarize 18 published cases and present new cases dating to the Neolithic period and introduce 3 newly discovered Bronze Age cases of TB. Despite extensive research, TB has not yet been identified from the Copper and Iron Ages in the Carpathian Basin. Considering the state of TB research, and supplemented by our prehistoric dataset, the spatio-temporal pattern of the disease can be further elucidated, thus advancing future molecular and paleopathological studies. Our dataset offers comprehensive spatial and temporal information on the spread of the disease in the Carpathian Basin, along with a detailed biological profile of the demonstrated cases and extensive paleopathological descriptions of the observed lesions, complemented by photographic evidence. This invaluable resource paves the way for enhanced understanding and progress in the field.

Tuberculosis in mummies - New findings, perspectives and limitations

The molecular analysis of ancient pathogen DNA represents a unique opportunity for the study of infectious diseases in ancient human remains. Among other diseases, paleogenetic studies have been successful in detecting tuberculous DNA in ancient human remains. In the beginning of ancient DNA (aDNA) studies, the presence of tuberculosis (TB) DNA was assessed using a PCR-based assay targeting specific regions of the Mycobacterium tuberculosis (MTB) complex, such as the repetitive element IS6110. The advent of high-throughput sequencing has enabled the reconstruction of full ancient TB genomes in the field of paleomicrobiology. However, despite the numerous paleopathological and PCR-based studies on the presence of tuberculosis in historic human remains, full genome wide reconstructions are still limited to well-preserved specimens with low environmental contamination and connected with extensive screening efforts. This has led to some controversies regarding the evolutionary history of its causative agent Mycobacterium tuberculosis. In this context, mummies have been shown to be a good source for the detection of MTB complex DNA due to a low exposure to environmental influences and the overall good state of preservation of hard and soft tissues in the human remains. Here, we present the major findings on the presence of TB infections in the 18th century naturally mummified human remains from Vác, Hungary and the current status of the detection of MTB complex DNA in mummified human remains. The future perspectives of detecting tuberculosis in mummies will be discussed in the light of methodological aspects, as well as ethical and curational challenges.

Parallel signatures of Mycobacterium tuberculosis and human Y-chromosome phylogeography support the Two Layer model of East Asian population history

The Two Layer hypothesis is fast becoming the favoured narrative describing East Asian population history. Under this model, hunter-gatherer groups who initially peopled East Asia via a route south of the Himalayas were assimilated by agriculturalist migrants who arrived via a northern route across Eurasia. A lack of ancient samples from tropical East Asia limits the resolution of this model. We consider insight afforded by patterns of variation within the human pathogen Mycobacterium tuberculosis (Mtb) by analysing its phylogeographic signatures jointly with the human Y-chromosome. We demonstrate the Y-chromosome lineages enriched in the traditionally hunter-gatherer groups associated with East Asia's first layer of peopling to display deep roots, low long-term effective population size, and diversity patterns consistent with a southern entry route. These characteristics mirror those of the evolutionarily ancient Mtb lineage 1. The remaining East Asian Y-chromosome lineage is almost entirely absent from traditionally hunter-gatherer groups and displays spatial and temporal characteristics which are incompatible with a southern entry route, and which link it to the development of agriculture in modern-day China. These characteristics mirror those of the evolutionarily modern Mtb lineage 2. This model paves the way for novel host-pathogen coevolutionary research hypotheses in East Asia.

Metagenomic analysis reveals mixed Mycobacterium tuberculosis infection in a 18th century Hungarian midwife

The Vác Mummy Collection comprises 265 well documented mummified individuals from the late 16th to the early 18th century that were discovered in 1994 inside a crypt in Vác, Hungary. This collection offers a unique opportunity to study the relationship between humans and pathogens in the pre-antibiotic era, as previous studies have shown a high proportion of tuberculosis (TB) infections among the individuals. In this study, we recovered ancient DNA with shotgun sequencing from a rib bone sample of a 18th century midwife. This individual is part of the collection and shows clear skeletal changes that are associated with tuberculosis and syphilis. To provide molecular proof, we applied a metagenomic approach to screen for ancient pathogen DNA. While we were unsuccessful to recover any ancient Treponema pallidum DNA, we retrieved high coverage ancient TB DNA and identified a mixed infection with two distinct TB strains by detailed single-nucleotide polymorphism and phylogenetic analysis. Thereby, we have obtained comprehensive results demonstrating the long-time prevalence of mixed infections with the sublineages L4.1.2.1/Haarlem and L4.10/PGG3 within the local community in preindustrial Hungary and put them in context of sociohistorical factors.

Biocultural perspectives of infectious diseases and demographic evolution: Tuberculosis and its comorbidities through history

Anthropologists recognize the importance of conceptualizing health in the context of the mutually evolving nature of biology and culture through the biocultural approach, but biocultural anthropological perspectives of infectious diseases and their impacts on humans (and vice versa) through time are relatively underrepresented. Tuberculosis (TB) has been a constant companion of humans for thousands of years and has heavily influenced population health in almost every phase of cultural and demographic evolution. TB in human populations has been dramatically influenced by behavior, demographic and epidemiological shifts, and other comorbidities through history. This paper critically discusses TB and some of its major comorbidities through history within a biocultural framework to show how transitions in human demography and culture affected the disease-scape of TB. In doing so, I address the potential synthesis of biocultural and epidemiological transition theory to better comprehend the mutual evolution of infectious diseases and humans.

Recent progress on methods for estimating and updating large phylogenies

With the increased availability of sequence data and even of fully sequenced and assembled genomes, phylogeny estimation of very large trees (even of hundreds of thousands of sequences) is now a goal for some biologists. Yet, the construction of these phylogenies is a complex pipeline presenting analytical and computational challenges, especially when the number of sequences is very large. In the past few years, new methods have been developed that aim to enable highly accurate phylogeny estimations on these large datasets, including divide-and-conquer techniques for multiple sequence alignment and/or tree estimation, methods that can estimate species trees from multi-locus datasets while addressing heterogeneity due to biological processes (e.g. incomplete lineage sorting and gene duplication and loss), and methods to add sequences into large gene trees or species trees. Here we present some of these recent advances and discuss opportunities for future improvements. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.

Evolution of Mycobacterium tuberculosis drug resistance in the genomic era

Mycobacterium tuberculosis has acquired drug resistance to all drugs that have been used against it, including those only recently introduced into clinical practice. Compared to other bacteria, it has a well conserved genome due to its role as an obligate human pathogen that has adapted to a niche over five to ten thousand years. These features facilitate reconstruction and dating of M. tuberculosis phylogenies, giving key insights into how resistance has been acquired and spread globally. Resistance to each new drug has occurred within five to ten years of clinical use and has occurred even more rapidly with recently introduced drugs. In most cases, resistance-conferring mutations come with a fitness cost, but this can be overcome by compensatory mutations which restore fitness to that of wild-type bacteria. It is likely that M. tuberculosis acquires drug resistance while maintaining limited genomic variability due the generation of low frequency within-host variation, combined with ongoing purifying selection causing loss of variants without a clear fitness advantage. However, variants that do confer an advantage, such as drug resistance, can increase in prevalence amongst all bacteria within a host and become the dominant clone. These resistant strains can then be transmitted leading to primary drug resistant infection in a new host. As many countries move towards genomic methods for diagnosis of M. tuberculosis infection and drug resistance, it is important to be aware of the implications for the evolution of resistance. Currently, understanding of resistance-conferring mutations is incomplete, and some targeted genetic diagnostics create their own selective pressures. We discuss an example where a rifampicin resistance-conferring mutation which was not routinely covered by standard testing became dominant. Finally, resistance to new drugs such as bedaquiline and delamanid is caused by individually rare mutations occurring across a large mutational genomic target that have been detected over a short time, and do not provide statistical power for genotype-phenotype correlation – in contrast to longer-established drugs that form the backbone of drug-sensitive antituberculosis therapy. Therefore, we need a different approach to identify resistance-conferring mutations of new drugs before their resistance becomes widespread, abrogating their usefulness.

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.

M. tuberculosis microvariation is common and is associated with transmission: analysis of three years prospective universal sequencing in England

BACKGROUND

The prevalence, association with disease status, and public health impact of infection with mixtures of M. tuberculosis strains is unclear, in part due to limitations of existing methods for detecting mixed infections.

METHODS

We developed an algorithm to identify mixtures of M. tuberculosis strains using next generation sequencing data, assessing performance using simulated sequences. We identified mixed M. tuberculosis strains when there was at least one mixed nucleotide position, and where both the mixture's components were present in similar isolates from other individuals, compatible with transmission of the component strains. We determined risk factors for mixed infection among isolations of M. tuberculosis in England using logistic regression. We used survival analyses to assess the association between mixed infection and putative transmission.

FINDINGS

6,560 isolations of TB were successfully sequenced in England 2016-2018. Of 3,691 (56%) specimens for which similar sequences had been isolated from at least two other individuals, 341 (9.2%) were mixed. Mixed infection was more common in lineages other than Lineage 4. Among the 1,823 individuals with pulmonary infection with Lineage 4 M. tuberculosis, mixed infection was associated with significantly increased risk of subsequent isolation of closely related organisms from a different individual (HR 1.43, 95% CI 1.05,1.94), indicative of transmission.

INTERPRETATION

Mixtures of transmissible strains occur in at least 5% of tuberculosis infections in England; when present in pulmonary disease, such mixtures are associated with an increased risk of tuberculosis transmission.

FUNDING

Public Health England; NIHR Health Protection Research Units; European Union.

The Microbiome as Part of the Contemporary View of Tuberculosis Disease

The study of the microbiome has changed our overall perspective on health and disease. Although studies of the lung microbiome have lagged behind those on the gastrointestinal microbiome, there is now evidence that the lung microbiome is a rich, dynamic ecosystem. Tuberculosis is one of the oldest human diseases, it is primarily a respiratory infectious disease caused by strains from the Mycobacterium tuberculosis Complex. Even today, during the COVID-19 pandemic, it remains one of the principal causes of morbidity and mortality worldwide. Tuberculosis disease manifests itself as a dynamic spectrum that ranges from asymptomatic latent infection to life-threatening active disease. The review aims to provide an overview of the microbiome in the tuberculosis setting, both in patients’ and animal models. We discuss the relevance of the microbiome and its dysbiosis, and how, probably through its interaction with the immune system, it is a significant factor in tuberculosis’s susceptibility, establishment, and severity.

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.

OUP accepted manuscript

Escherichia coli has a rich history as biology's ‘rock star’, driving advances across many fields. In the wild, E. coli resides innocuously in the gut of humans and animals but is also a versatile pathogen commonly associated with intestinal and extraintestinal infections and antimicrobial resistance—including large foodborne outbreaks such as the one that swept across Europe in 2011, killing 54 individuals and causing approximately 4000 infections and 900 cases of haemolytic uraemic syndrome. Given that most E. coli are harmless gut colonizers, an important ecological question plaguing microbiologists is what makes E. coli an occasionally devastating pathogen? To address this question requires an enhanced understanding of the ecology of the organism as a commensal. Here, we review how our knowledge of the ecology and within-host diversity of this organism in the vertebrate gut has progressed in the 137 years since E. coli was first described. We also review current approaches to the study of within-host bacterial diversity. In closing, we discuss some of the outstanding questions yet to be addressed and prospects for future research.

Genomic analysis of an outbreak of bovine tuberculosis in a man-made multi-host species system: A call for action on wildlife in Brazil

We report on a 15-year-long outbreak of bovine tuberculosis (bTB) in wildlife from a Brazilian safari park. A timeline of diagnostic events and whole-genome sequencing (WGS) of 21 Mycobacterium bovis isolates from deer and llamas were analyzed. Accordingly, from 2003 to 2018, at least 16 animals, from eight species, died due to TB, which is likely an underestimated number. In three occasions since 2013, the deer presented positive tuberculin tests, leading to the park closure and culling of all deer. WGS indicated that multiple M. bovis strains were circulating, with at least three founding introductions since the park inauguration in 1977. Using a previously sequenced dataset of 71 M. bovis genomes from cattle, we found no recent transmission events between nearby farms and the park based on WGS. Lastly, by discussing socio-economic and environmental factors escaping current regulatory gaps that were determinant of this outbreak, we pledge for the development of a plan to report and control bTB in wildlife in Brazil.

A 3,000-year-old, basal S. enterica lineage from Bronze Age Xinjiang suggests spread along the Proto-Silk Road

Salmonella enterica (S. enterica) has infected humans for a long time, but its evolutionary history and geographic spread across Eurasia is still poorly understood. Here, we screened for pathogen DNA in 14 ancient individuals from the Bronze Age Quanergou cemetery (XBQ), Xinjiang, China. In 6 individuals we detected S. enterica. We reconstructed S. enterica genomes from those individuals, which form a previously undetected phylogenetic branch basal to Paratyphi C, Typhisuis and Choleraesuis-the so-called Para C lineage. Based on pseudogene frequency, our analysis suggests that the ancient S. enterica strains were not host adapted. One genome, however, harbors the Salmonella pathogenicity island 7 (SPI-7), which is thought to be involved in (para)typhoid disease in humans. This offers first evidence that SPI-7 was acquired prior to the emergence of human-adapted Paratyphi C around 1,000 years ago. Altogether, our results show that Salmonella enterica infected humans in Eastern Eurasia at least 3,000 years ago, and provide the first ancient DNA evidence for the spread of a pathogen along the Proto-Silk Road.

Multiple Merger Genealogies in Outbreaks of Mycobacterium tuberculosis

The Kingman coalescent and its developments are often considered among the most important advances in population genetics of the last decades. Demographic inference based on coalescent theory has been used to reconstruct the population dynamics and evolutionary history of several species, including Mycobacterium tuberculosis (MTB), an important human pathogen causing tuberculosis. One key assumption of the Kingman coalescent is that the number of descendants of different individuals does not vary strongly, and violating this assumption could lead to severe biases caused by model misspecification. Individual lineages of MTB are expected to vary strongly in reproductive success because 1) MTB is potentially under constant selection due to the pressure of the host immune system and of antibiotic treatment, 2) MTB undergoes repeated population bottlenecks when it transmits from one host to the next, and 3) some hosts show much higher transmission rates compared with the average (superspreaders).

Here, we used an approximate Bayesian computation approach to test whether multiple-merger coalescents (MMC), a class of models that allow for large variation in reproductive success among lineages, are more appropriate models to study MTB populations. We considered 11 publicly available whole-genome sequence data sets sampled from local MTB populations and outbreaks and found that MMC had a better fit compared with the Kingman coalescent for 10 of the 11 data sets. These results indicate that the null model for analyzing MTB outbreaks should be reassessed and that past findings based on the Kingman coalescent need to be revisited.

New insights into human immunity from ancient genomics

Population genetic studies have clearly indicated that immunity and host defense are among the functions most frequently subject to natural selection, and increased our understanding of the biological relevance of the corresponding genes and their contribution to variable immune traits and diseases. Herein, we will focus on some recently studied forms of human adaptation to infectious agents, including hybridization with now-extinct hominins, such as Neanderthals and Denisovans, and admixture between modern human populations. These studies, which are partly enabled by the technological advances in the sequencing of DNA from ancient remains, provide new insight into the sources of immune response variation in contemporary humans, such as the recently reported link between Neanderthal heritage and susceptibility to severe COVID-19 disease. Furthermore, ancient DNA analyses, in both humans and pathogens, allow to measure the action of natural selection on immune genes across time and to reconstruct the impact of past epidemics on the evolution of human immunity.

Human ancient DNA analyses reveal the high burden of tuberculosis in Europeans over the last 2,000 years

Tuberculosis (TB), usually caused by Mycobacterium tuberculosis bacteria, is the first cause of death from an infectious disease at the worldwide scale, yet the mode and tempo of TB pressure on humans remain unknown. The recent discovery that homozygotes for the P1104A polymorphism of TYK2 are at higher risk to develop clinical forms of TB provided the first evidence of a common, monogenic predisposition to TB, offering a unique opportunity to inform on human co-evolution with a deadly pathogen. Here, we investigate the history of human exposure to TB by determining the evolutionary trajectory of the TYK2 P1104A variant in Europe, where TB is considered to be the deadliest documented infectious disease. Leveraging a large dataset of 1,013 ancient human genomes and using an approximate Bayesian computation approach, we find that the P1104A variant originated in the common ancestors of West Eurasians ∼30,000 years ago. Furthermore, we show that, following large-scale population movements of Anatolian Neolithic farmers and Eurasian steppe herders into Europe, P1104A has markedly fluctuated in frequency over the last 10,000 years of European history, with a dramatic decrease in frequency after the Bronze Age. Our analyses indicate that such a frequency drop is attributable to strong negative selection starting ∼2,000 years ago, with a relative fitness reduction on homozygotes of 20%, among the highest in the human genome. Together, our results provide genetic evidence that TB has imposed a heavy burden on European health over the last two millennia.

Methods for Detecting Mycobacterial Mixed Strain Infections-A Systematic Review

Mixed strain infection (MSI) refers to the concurrent infection of a susceptible host with multiple strains of a single pathogenic species. Known to occur in humans and animals, MSIs deserve special consideration when studying transmission dynamics, evolution, and treatment of mycobacterial diseases, notably tuberculosis in humans and paratuberculosis (or Johne's disease) in ruminants. Therefore, a systematic review was conducted to examine how MSIs are defined in the literature, how widespread the phenomenon is across the host species spectrum, and to document common methods used to detect such infections. Our search strategy identified 121 articles reporting MSIs in both humans and animals, the majority (78.5%) of which involved members of the Mycobacterium tuberculosis complex, while only a few (21.5%) examined non-tuberculous mycobacteria (NTM). In addition, MSIs exist across various host species, but most reports focused on humans due to the extensive amount of work done on tuberculosis. We reviewed the strain typing methods that allowed for MSI detection and found a few that were commonly employed but were associated with specific challenges. Our review notes the need for standardization, as some highly discriminatory methods are not adapted to distinguish between microevolution of one strain and concurrent infection with multiple strains. Further research is also warranted to examine the prevalence of NTM MSIs in both humans and animals. In addition, it is envisioned that the accurate identification and a better understanding of the distribution of MSIs in the future will lead to important information on the epidemiology and pathophysiology of mycobacterial diseases.

Verification of tuberculosis infection among Vác mummies (18th century CE, Hungary) based on lipid biomarker profiling with a new HPLC-HESI-MS approach

Tuberculosis (TB) was a large burden of infections that peaked during the 19th century in Europe. Mummies from the 18th century CE, discovered in the crypt of a church at Vác, Hungary, had high TB prevalence, as revealed by amplification of key fragments of TB DNA and genome-wide TB analysis. Complementary methods are needed to confirm these diagnoses and one approach uses the identification of specific lipid biomarkers, such as TB mycocerosic acids (MCs). Previously, MC derivatives were profiled by specialised gas chromatography-mass spectrometry (GC-MS), so an alternative more direct approach has been developed. Underivatized MCs are extracted and analysed by high-performance liquid chromatography linked to a mass spectrometer, in heated electrospray ionisation mode (HPLC-HESI-MS). The method was validated using representatives of the Mycobacterium tuberculosis complex and other mycobacteria and tested on six Vác mummy cases, previously considered positive for TB infection. Analysing both rib and soft tissue samples, four out of six cases gave profiles of main C32 and major C29 and C39 mycocerosates correlating well with those of M. tuberculosis. Multidisciplinary methods are needed in the diagnosis of ancient tuberculosis; this new protocol accesses important confirmatory evidence, as demonstrated by the confirmation of TB in the Vác mummies.

Successful expansion of Mycobacterium tuberculosis Latin American and Mediterranean sublineage (L4.3/LAM) in Tunisia mainly driven by a single, long-established clonal complex

OBJECTIVES

To explore the evolutionary history of Mycobacterium tuberculosis Latin American and Mediterranean (L4.3/LAM) sublineage in Tunisia, where it predominates.

METHODS

High-resolution genotyping of 252 L4.3/LAM clinical strains was undertaken, and whole-genome sequencing was performed on 31 representative isolates.

RESULTS

Genotyping data coupled with Bayesian analyses split the Tunisian L4.3/LAM strain collection into two divergent entities (65.07% vs 34.92%): a major subpopulation, dominated by a single clonal complex (CC), TUN4.3_CC1 (94.51%); and a minor subpopulation, dominated by TUN4.3_CC2 (42.04%). TUN4.3_CC1 is clearly thriving in Tunisia, accounting for 61.5% of the L4.3/LAM sublineage. TUN4.3_CC1 displayed higher mean allelic richness compared with TUN4.3_CC2 and predominated throughout the entire region, indicating a long-established history. The very low proportion of drug resistance among TUN4.3_CC1 isolates is indicative of their intrinsic ability to spread successfully in the host population. Genomic analyses further confirmed the clear genetic separation between the two main CCs (pairwise fixation index 0.56), and suggested the relatively ancient origin of TUN4.3_CC1. Consistent with its successful expansion, TUN4.3_CC1 showed reduced mean pairwise genetic distance between genomes.

CONCLUSIONS

These findings link the successful expansion of L4.3/LAM in Tunisia to a single long-established clone.

Metagenomics of the modern and historical human oral microbiome with phylogenetic studies on Streptococcus mutans and Streptococcus sobrinus

We have recently developed bioinformatic tools to accurately assign metagenomic sequence reads to microbial taxa: SPARSE for probabilistic, taxonomic classification of sequence reads; EToKi for assembling and polishing genomes from short-read sequences; and GrapeTree, a graphic visualizer of genetic distances between large numbers of genomes. Together, these methods support comparative analyses of genomes from ancient skeletons and modern humans. Here, we illustrate these capabilities with 784 samples from historical dental calculus, modern saliva and modern dental plaque. The analyses revealed 1591 microbial species within the oral microbiome. We anticipated that the oral complexes of Socransky et al., which were defined in 1998, would predominate among taxa whose frequencies differed by source. However, although some species discriminated between sources, we could not confirm the existence of the complexes. The results also illustrate further functionality of our pipelines with two species that are associated with dental caries, Streptococcus mutans and Streptococcus sobrinus. They were rare in historical dental calculus but common in modern plaque, and even more common in saliva. Reconstructed draft genomes of these two species from metagenomic samples in which they were abundant were combined with modern public genomes to provide a detailed overview of their core genomic diversity. This article is part of the theme issue 'Insights into health and disease from ancient biomolecules'.

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

The population genomics of within-host Mycobacterium tuberculosis

Recent progress in genomic sequencing from patient samples has allowed for the first detailed insight into the within-host genetic diversity of Mycobacterium tuberculosis (M.TB), revealing remarkably low levels of variation. While this has often been attributed to low mutation rates, other factors have been described, including resistance evolution (i.e., selective sweeps), widespread purifying and background selection, and, more recently, progeny skew. Here we review recent findings pertaining to the processes governing the evolutionary dynamics of M.TB, discuss their implications for improving our understanding of this important human pathogen, and make recommendations for future work. Significantly, this emerging evolutionary framework involving the joint estimation of demographic, selective, and reproductive processes is forming a new paradigm for the study of within-host pathogen evolution that will be widely applicable across organisms.

'TB or not TB': the conundrum of pre-European contact tuberculosis in the Pacific

Tuberculosis (TB) is a major global health threat, infecting one-third of the world's population. Despite this prominence, the age, origin and spread of the disease have been topics of contentious debate. Molecular studies suggest that Mycobacterium tuberculosis 'sensu stricto', the most common strain of TB infecting humans today, originated in Africa and from there spread into Europe and Asia. The M. tuberculosis strains most commonly found across the Pacific and the Americas today are most closely related to European strains, supporting a hypothesis that the disease only reached these regions relatively recently via European sailors or settlers. However, this hypothesis is inconsistent with palaeopathological evidence of TB-like lesions in human remains from across the Pacific that predate European contact. Similarly, genetic evidence from pre-European South American mummies challenges the notion of a European introduction of the disease into the Pacific. Here, we review the complex evidence for the age and origin of TB in the Pacific, and discuss key gaps in our knowledge and how these may be addressed. This article is part of the theme issue 'Insights into health and disease from ancient biomolecules'.

Oldest evidence of tuberculosis in Argentina: A multidisciplinary investigation in an adult male skeleton from Saujil, Tinogasta, Catamarca (905-1030 CE)

The Mycobacterium tuberculosis complex (MTC) has affected South American populations since ca. 200 years BCE. In Argentina, possible cases date from ca. 1000-1400 Common Era (CE). This paper describes the oldest (905-1030 CE) confirmed case of tuberculosis (TB) in a young adult male from Lomitas de Saujil (Tinogasta, Catamarca, Argentina). Osteolytic lesions on the bodies of the lower spine were macroscopically and radiographically identified. Bilateral new bone formation was seen on the visceral vertebral third of several ribs and in long bones, compatible with hypertrophic osteoarthropathy. Representative rib and hand bones gave profiles for MTC-specific C₂₇-C₃₂ mycocerosic acid lipid biomarkers; these were strongest in one heavily-lesioned lower rib, which also had MTC-diagnostic C₇₆-C₈₉ mycolic acids and positive amplification of MTC-typical IS6110 aDNA fragments. During the first millennium CE, the intense social interaction, the spatial circumscription of villages among the pre-Hispanic societies in the mesothermal valleys of Catamarca and the fluid contacts with the Eastern lowlands, valleys and puna, were factors likely to favor disease transmission. It is proposed that TB arrived from northern Chile and dispersed towards the northeast into the Yocavil valley, where several cases of TB infection were macroscopically identified for a later chronology.

A seventeenth-century Mycobacterium tuberculosis genome supports a Neolithic emergence of the Mycobacterium tuberculosis complex

BACKGROUND

Although tuberculosis accounts for the highest mortality from a bacterial infection on a global scale, questions persist regarding its origin. One hypothesis based on modern Mycobacterium tuberculosis complex (MTBC) genomes suggests their most recent common ancestor followed human migrations out of Africa approximately 70,000 years before present. However, studies using ancient genomes as calibration points have yielded much younger dates of less than 6000 years. Here, we aim to address this discrepancy through the analysis of the highest-coverage and highest-quality ancient MTBC genome available to date, reconstructed from a calcified lung nodule of Bishop Peder Winstrup of Lund (b. 1605-d. 1679).

RESULTS

A metagenomic approach for taxonomic classification of whole DNA content permitted the identification of abundant DNA belonging to the human host and the MTBC, with few non-TB bacterial taxa comprising the background. Genomic enrichment enabled the reconstruction of a 141-fold coverage M. tuberculosis genome. In utilizing this high-quality, high-coverage seventeenth-century genome as a calibration point for dating the MTBC, we employed multiple Bayesian tree models, including birth-death models, which allowed us to model pathogen population dynamics and data sampling strategies more realistically than those based on the coalescent.

CONCLUSIONS

The results of our metagenomic analysis demonstrate the unique preservation environment calcified nodules provide for DNA. Importantly, we estimate a most recent common ancestor date for the MTBC of between 2190 and 4501 before present and for Lineage 4 of between 929 and 2084 before present using multiple models, confirming a Neolithic emergence for the MTBC.

The past, present and future of ancient bacterial DNA

Groundbreaking studies conducted in the mid-1980s demonstrated the possibility of sequencing ancient DNA (aDNA), which has allowed us to answer fundamental questions about the human past. Microbiologists were thus given a powerful tool to glimpse directly into inscrutable bacterial history, hitherto inaccessible due to a poor fossil record. Initially plagued by concerns regarding contamination, the field has grown alongside technical progress, with the advent of high-throughput sequencing being a breakthrough in sequence output and authentication. Albeit burdened with challenges unique to the analysis of bacteria, a growing number of viable sources for aDNA has opened multiple avenues of microbial research. Ancient pathogens have been extracted from bones, dental pulp, mummies and historical medical specimens and have answered focal historical questions such as identifying the aetiological agent of the black death as Yersinia pestis. Furthermore, ancient human microbiomes from fossilized faeces, mummies and dental plaque have shown shifts in human commensals through the Neolithic demographic transition and industrial revolution, whereas environmental isolates stemming from permafrost samples have revealed signs of ancient antimicrobial resistance. Culminating in an ever-growing repertoire of ancient genomes, the quickly expanding body of bacterial aDNA studies has also enabled comparisons of ancient genomes to their extant counterparts, illuminating the evolutionary history of bacteria. In this review we summarize the present avenues of research and contextualize them in the past of the field whilst also pointing towards questions still to be answered.

Citywide Transmission of Multidrug-resistant Tuberculosis Under China's Rapid Urbanization: A Retrospective Population-based Genomic Spatial Epidemiological Study

BACKGROUND

Population movement could extend multidrug-resistant tuberculosis (MDR-TB) transmission and complicate its global prevalence. We sought to identify the high-risk populations and geographic sites of MDR-TB transmission in Shenzhen, the most common destination for internal migrants in China.

METHODS

We performed a population-based, retrospective study in patients diagnosed with MDR-TB in Shenzhen during 2013-2017. By defining genomic clusters with a threshold of 12-single-nucleotide polymorphism distance based on whole-genome sequencing of their clinical strains, the clustering rate was calculated to evaluate the level of recent transmission. Risk factors were identified by multivariable logistic regression. To further delineate the epidemiological links, we invited the genomic-clustered patients to an in-depth social network investigation.

RESULTS

In total, 105 (25.2%) of the 417 enrolled patients with MDR-TB were grouped into 40 genome clusters, suggesting recent transmission of MDR strains. The adjusted risk for student to have a clustered strain was 4.05 (95% confidence interval, 1.06-17.0) times greater than other patients. The majority (70%, 28/40) of the genomic clusters involved patients who lived in different districts, with residences separated by an average of 8.76 kilometers. Other than household members, confirmed epidemiological links were also identified among classmates and workplace colleagues.

CONCLUSIONS

These findings demonstrate that local transmission of MDR-TB is a serious problem in Shenzhen. While most transmission occurred between people who lived distant from each other, there was clear evidence that transmission occurred in schools and workplaces, which should be included as targeted sites for active case finding.The average residential distance between genomic-clustered cases was more than 8 kilometers, while schools and workplaces, identified as sites of transmission in this study, deserve increased vigilance for targeted case finding of multidrug-resistant tuberculosis.

Exploration of space to achieve scientific breakthroughs

Living organisms adapt to changing environments using their amazing flexibility to remodel themselves by a process called evolution. Environmental stress causes selective pressure and is associated with genetic and phenotypic shifts for better modifications, maintenance, and functioning of organismal systems. The natural evolution process can be used in complement to rational strain engineering for the development of desired traits or phenotypes as well as for the production of novel biomaterials through the imposition of one or more selective pressures. Space provides a unique environment of stressors (e.g., weightlessness and high radiation) that organisms have never experienced on Earth. Cells in the outer space reorganize and develop or activate a range of molecular responses that lead to changes in cellular properties. Exposure of cells to the outer space will lead to the development of novel variants more efficiently than on Earth. For instance, natural crop varieties can be generated with higher nutrition value, yield, and improved features, such as resistance against high and low temperatures, salt stress, and microbial and pest attacks. The review summarizes the literature on the parameters of outer space that affect the growth and behavior of cells and organisms as well as complex colloidal systems. We illustrate an understanding of gravity-related basic biological mechanisms and enlighten the possibility to explore the outer space environment for application-oriented aspects. This will stimulate biological research in the pursuit of innovative approaches for the future of agriculture and health on Earth.

Advances in the molecular detection of tuberculosis in pre-contact Andean South America

Andean paleopathological research has significantly enhanced knowledge about the geographical distribution and evolution of tuberculosis (TB) in pre-Columbian South America. In this paper, we review the history and progress of research on ancient tuberculosis (TB) in the Andean region, focusing on the strengths and limitations of current approaches for the molecular detection of ancient pathogens, with special attention to TB. As a case study, we describe a molecular screening approach for the detection of ancient Mycobacterium tuberculosis in individuals from Late Intermediate Period (1000-1400 CE) contexts at the site of Huari, Peru. We evaluate 34 commingled human vertebrae and combine morphological assessments of pathology with high throughput sequencing and a non-selective approach to ancient pathogen DNA screening. Our method enabled the simultaneous detection of ancient M. tuberculosis DNA and an evaluation of the environmental microbial composition of each sample. Our results show that despite the dominance of environmental DNA, molecular signatures of M. tuberculosis were identified in eight vertebrae, six of which had no observable skeletal pathology classically associated tuberculosis infection. This screening approach will assist in the identification of candidate samples for downstream genomic analyses. The method permits higher resolution disease identification in cases where pathology may be absent, or where the archaeological context may necessitate a broad differential diagnosis based on morphology alone.

Global Distribution and Evolution of Mycobacterium bovis Lineages

Mycobacterium bovis is the main causative agent of zoonotic tuberculosis in humans and frequently devastates livestock and wildlife worldwide. Previous studies suggested the existence of genetic groups of M. bovis strains based on limited DNA markers (a.k.a. clonal complexes), and the evolution and ecology of this pathogen has been only marginally explored at the global level. We have screened over 2,600 publicly available M. bovis genomes and newly sequenced four wildlife M. bovis strains, gathering 1,969 genomes from 23 countries and at least 24 host species, including humans, to complete a phylogenomic analyses. We propose the existence of four distinct global lineages of M. bovis (Lb1, Lb2, Lb3, and Lb4) underlying the current disease distribution. These lineages are not fully represented by clonal complexes and are dispersed based on geographic location rather than host species. Our data divergence analysis agreed with previous studies reporting independent archeological data of ancient M. bovis (South Siberian infected skeletons at ∼2,000 years before present) and indicates that extant M. bovis originated between 715 and 3,556 years BP, with later emergence in the New World and Oceania, likely influenced by trades among countries.

TbD1 deletion as a driver of the evolutionary success of modern epidemic Mycobacterium tuberculosis lineages

Mycobacterium tuberculosis (Mtb) strains are classified into different phylogenetic lineages (L), three of which (L2/L3/L4) emerged from a common progenitor after the loss of the MmpS6/MmpL6-encoding Mtb-specific deletion 1 region (TbD1). These TbD1-deleted “modern” lineages are responsible for globally-spread tuberculosis epidemics, whereas TbD1-intact “ancestral” lineages tend to be restricted to specific geographical areas, such as South India and South East Asia (L1) or East Africa (L7). By constructing and characterizing a panel of recombinant TbD1-knock-in and knock-out strains and comparison with clinical isolates, here we show that deletion of TbD1 confers to Mtb a significant increase in resistance to oxidative stress and hypoxia, which correlates with enhanced virulence in selected cellular, guinea pig and C3HeB/FeJ mouse infection models, the latter two mirroring in part the development of hypoxic granulomas in human disease progression. Our results suggest that loss of TbD1 at the origin of the L2/L3/L4 Mtb lineages was a key driver for their global epidemic spread and outstanding evolutionary success.

Mycobacterium tuberculosis (Mtb) modern strains emerged from a common progenitor after the loss of Mtb-specific deletion 1 region (TbD1). Here, the authors show that deletion of TbD1 correlates with enhanced Mtb virulence in animal models, mirroring the development of hypoxic granulomas in human disease progression.

Comparative genomics confirms a rare melioidosis human-to-human transmission event and reveals incorrect phylogenomic reconstruction due to polyclonality

Human-to-human transmission of the melioidosis bacterium, Burkholderia pseudomallei, is exceedingly rare, with only a handful of suspected cases documented to date. Here, we used whole-genome sequencing (WGS) to characterize one such unusual B. pseudomallei transmission event, which occurred between a breastfeeding mother with mastitis and her child. Two strains corresponding to multilocus sequence types (STs)-259 and -261 were identified in the mother's sputum from both the primary culture sweep and in purified colonies, confirming an unusual polyclonal infection in this patient. In contrast, primary culture sweeps of the mother's breast milk and the child's cerebrospinal fluid and blood samples contained only ST-259, indicating monoclonal transmission to the child. Analysis of purified ST-259 isolates showed no genetic variation between mother and baby isolates, providing the strongest possible evidence of B. pseudomallei human-to-human transmission, probably via breastfeeding. Next, phylogenomic analysis of all isolates, including the mother's mixed ST-259/ST-261 sputum sample, was performed to investigate the effects of mixtures on phylogenetic inference. Inclusion of this mixture caused a dramatic reduction in the number of informative SNPs, resulting in branch collapse of ST-259 and ST-261 isolates, and several instances of incorrect topology in a global B. pseudomallei phylogeny, resulting in phylogenetic incongruence. Although phylogenomics can provide clues about the presence of mixtures within WGS datasets, our results demonstrate that this methodology can lead to phylogenetic misinterpretation if mixed genomes are not correctly identified and omitted. Using current bioinformatic tools, we demonstrate a robust method for bacterial mixture identification and strain parsing that avoids these pitfalls.

Resistance Sniffer: An online tool for prediction of drug resistance patterns of Mycobacterium tuberculosis isolates using next generation sequencing data

The effective control of multidrug resistant tuberculosis (MDR-TB) relies upon the timely diagnosis and correct treatment of all tuberculosis cases. Whole genome sequencing (WGS) has great potential as a method for the rapid diagnosis of drug resistant Mycobacterium tuberculosis (Mtb) isolates. This method overcomes most of the problems that are associated with current phenotypic drug susceptibility testing. However, the application of WGS in the clinical setting has been deterred by data complexities and skill requirements for implementing the technologies as well as clinical interpretation of the next generation sequencing (NGS) data. The proposed diagnostic application was drawn upon recent discoveries of patterns of Mtb clade-specific genetic polymorphisms associated with antibiotic resistance. A catalogue of genetic determinants of resistance to thirteen anti-TB drugs for each phylogenetic clade was created. A computational algorithm for the identification of states of diagnostic polymorphisms was implemented as an online software tool, Resistance Sniffer (http://resistance-sniffer.bi.up.ac.za/), and as a stand-alone software tool to predict drug resistance in Mtb isolates using complete or partial genome datasets in different file formats including raw Illumina fastq read files. The program was validated on sequenced Mtb isolates with data on antibiotic resistance trials available from GMTV database and from the TB Platform of South African Medical Research Council (SAMRC), Pretoria. The program proved to be suitable for probabilistic prediction of drug resistance profiles of individual strains and large sequence data sets.

Evolution of Mycobacterium tuberculosis complex lineages and their role in an emerging threat of multidrug resistant tuberculosis in Bamako, Mali

In recent years Bamako has been faced with an emerging threat from multidrug resistant TB (MDR-TB). Whole genome sequence analysis was performed on a subset of 76 isolates from a total of 208 isolates recovered from tuberculosis patients in Bamako, Mali between 2006 and 2012. Among the 76 patients, 61(80.3%) new cases and 15(19.7%) retreatment cases, 12 (16%) were infected by MDR-TB. The dominant lineage was the Euro-American lineage, Lineage 4. Within Lineage 4, the Cameroon genotype was the most prevalent genotype (n = 20, 26%), followed by the Ghana genotype (n = 16, 21%). A sub-clade of the Cameroon genotype, which emerged ~22 years ago was likely to be involved in community transmission. A sub-clade of the Ghana genotype that arose approximately 30 years ago was an important cause of MDR-TB in Bamako. The Ghana genotype isolates appeared more likely to be MDR than other genotypes after controlling for treatment history. We identified a clade of four related Beijing isolates that included one MDR-TB isolate. It is a major concern to find the Cameroon and Ghana genotypes involved in community transmission and MDR-TB respectively. The presence of the Beijing genotype in Bamako remains worrying, given its high transmissibility and virulence.

HOPS: automated detection and authentication of pathogen DNA in archaeological remains

High-throughput DNA sequencing enables large-scale metagenomic analyses of complex biological systems. Such analyses are not restricted to present-day samples and can also be applied to molecular data from archaeological remains. Investigations of ancient microbes can provide valuable information on past bacterial commensals and pathogens, but their molecular detection remains a challenge. Here, we present HOPS (Heuristic Operations for Pathogen Screening), an automated bacterial screening pipeline for ancient DNA sequences that provides detailed information on species identification and authenticity. HOPS is a versatile tool for high-throughput screening of DNA from archaeological material to identify candidates for genome-level analyses.

Revisiting the tuberculosis and leprosy cross-immunity hypothesis: Expanding the dialogue between immunology and paleopathology

OBJECTIVE

Our primary objective is to re-visit the tuberculosis and leprosy cross-immunity. hypothesis through the careful integration of immunology and paleopathology.

METHODS

Using an integrated theoretical analysis that evaluates clinical literature on human innate immunological responses, paleomicrobiology, bioarchaeology, and paleopathology, we develop a multifactorial model.

RESULTS

Past populations do not represent homogeneous immunological landscapes, and therefore it is likely that leprosy in Medieval Europe did not uniformly decline due to cross-immunity.

CONCLUSIONS

We recommend that bioarchaeological reconstructions of past disease experience take into consideration models that include variation in immune function based on past environments and social contexts. This provides a unique opportunity to conduct comprehensive analyses on complex immunological processes.

SIGNIFICANCE

Extrapolating results from experimental immunology to larger populations elucidates complexities of disease cross-immunity and highlights the importance of synthesizing archaeological, social, paleopathological and biological data as a means of understanding disease in the past.

LIMITATIONS

All extrapolations from data produced from in vitro studies to past populations, using living donors, pose significant limitations where, among other factors, the full reconstruction of past environmental and social contexts can frequently be sparse or incomplete.

SUGGESTIONS FOR FUTURE RESEARCH

To reduce the limitations of integrating experimental immunology with bioarchaeological reconstructions (i.e. how to use skeletal samples to reconstruct inflammatory phenotypes), we propose that osteoimmunology, or the study of the interplay between immune cells and bone cells, should be considered a vital discipline and perhaps the foundation for the expansion of paleoimmunology.

The molecular clock of Mycobacterium tuberculosis

The molecular clock and its phylogenetic applications to genomic data have changed how we study and understand one of the major human pathogens, Mycobacterium tuberculosis (MTB), the etiologic agent of tuberculosis. Genome sequences of MTB strains sampled at different times are increasingly used to infer when a particular outbreak begun, when a drug-resistant clone appeared and expanded, or when a strain was introduced into a specific region. Despite the growing importance of the molecular clock in tuberculosis research, there is a lack of consensus as to whether MTB displays a clocklike behavior and about its rate of evolution. Here we performed a systematic study of the molecular clock of MTB on a large genomic data set (6,285 strains), covering different epidemiological settings and most of the known global diversity. We found that sampling times below 15-20 years were often insufficient to calibrate the clock of MTB. For data sets where such calibration was possible, we obtained a clock rate between 1x10-8 and 5x10-7 nucleotide changes per-site-per-year (0.04-2.2 SNPs per-genome-per-year), with substantial differences between clades. These estimates were not strongly dependent on the time of the calibration points as they changed only marginally when we used epidemiological isolates (sampled in the last 40 years) or three ancient DNA samples (about 1,000 years old) to calibrate the tree. Additionally, the uncertainty and the discrepancies in the results of different methods were sometimes large, highlighting the importance of using different methods, and of considering carefully their assumptions and limitations.