1
|
Derelle R, von Wachsmann J, Mäklin T, Hellewell J, Russell T, Lalvani A, Chindelevitch L, Croucher NJ, Harris SR, Lees JA. Seamless, rapid, and accurate analyses of outbreak genomic data using split k-mer analysis. Genome Res 2024; 34:1661-1673. [PMID: 39406504 PMCID: PMC11529842 DOI: 10.1101/gr.279449.124] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2024] [Accepted: 09/16/2024] [Indexed: 11/01/2024]
Abstract
Sequence variation observed in populations of pathogens can be used for important public health and evolutionary genomic analyses, especially outbreak analysis and transmission reconstruction. Identifying this variation is typically achieved by aligning sequence reads to a reference genome, but this approach is susceptible to reference biases and requires careful filtering of called genotypes. There is a need for tools that can process this growing volume of bacterial genome data, providing rapid results, but that remain simple so they can be used without highly trained bioinformaticians, expensive data analysis, and long-term storage and processing of large files. Here we describe split k-mer analysis (SKA2), a method that supports both reference-free and reference-based mapping to quickly and accurately genotype populations of bacteria using sequencing reads or genome assemblies. SKA2 is highly accurate for closely related samples, and in outbreak simulations, we show superior variant recall compared with reference-based methods, with no false positives. SKA2 can also accurately map variants to a reference and be used with recombination detection methods to rapidly reconstruct vertical evolutionary history. SKA2 is many times faster than comparable methods and can be used to add new genomes to an existing call set, allowing sequential use without the need to reanalyze entire collections. With an inherent absence of reference bias, high accuracy, and a robust implementation, SKA2 has the potential to become the tool of choice for genotyping bacteria. SKA2 is implemented in Rust and is freely available as open-source software.
Collapse
Affiliation(s)
- Romain Derelle
- NIHR Health Protection Research Unit in Respiratory Infections, National Heart and Lung Institute, Imperial College London, London W21PG, United Kingdom
| | - Johanna von Wachsmann
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton CB10 1SD, United Kingdom
| | - Tommi Mäklin
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton CB10 1SD, United Kingdom
- Department of Mathematics and Statistics, University of Helsinki, Helsinki 00014, Finland
| | - Joel Hellewell
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton CB10 1SD, United Kingdom
| | - Timothy Russell
- Centre for Mathematical Modelling of Infectious Diseases, London School of Hygiene & Tropical Medicine, London WC1E 7HT, United Kingdom
| | - Ajit Lalvani
- NIHR Health Protection Research Unit in Respiratory Infections, National Heart and Lung Institute, Imperial College London, London W21PG, United Kingdom
| | - Leonid Chindelevitch
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London W12 0BZ, United Kingdom
| | - Nicholas J Croucher
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London W12 0BZ, United Kingdom
| | - Simon R Harris
- Bill and Melinda Gates Foundation, Westminster, London SW1E 6AJ, United Kingdom
| | - John A Lees
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton CB10 1SD, United Kingdom;
| |
Collapse
|
2
|
Wood AJ, Benton CH, Delahay RJ, Marion G, Palkopoulou E, Pooley CM, Smith GC, Kao RR. The utility of whole-genome sequencing to identify likely transmission pairs for pathogens with slow and variable evolution. Epidemics 2024; 48:100787. [PMID: 39197305 DOI: 10.1016/j.epidem.2024.100787] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2024] [Revised: 07/03/2024] [Accepted: 08/14/2024] [Indexed: 09/01/2024] Open
Abstract
Pathogen whole-genome sequencing (WGS) has been used to track the transmission of infectious diseases in extraordinary detail, especially for pathogens that undergo fast and steady evolution, as is the case with many RNA viruses. However, for other pathogens evolution is less predictable, making interpretation of these data to inform our understanding of their epidemiology more challenging and the value of densely collected pathogen genome data uncertain. Here, we assess the utility of WGS for one such pathogen, in the "who-infected-whom" identification problem. We study samples from hosts (130 cattle, 111 badgers) with confirmed infection of M. bovis (causing bovine Tuberculosis), which has an estimated clock rate as slow as ∼0.1-1 variations per year. For each potential pathway between hosts, we calculate the relative likelihood that such a transmission event occurred. This is informed by an epidemiological model of transmission, and host life history data. By including WGS data, we shrink the number of plausible pathways significantly, relative to those deemed likely on the basis of life history data alone. Despite our uncertainty relating to the evolution of M. bovis, the WGS data are therefore a valuable adjunct to epidemiological investigations, especially for wildlife species whose life history data are sparse.
Collapse
Affiliation(s)
- A J Wood
- Roslin Institute, University of Edinburgh, United Kingdom
| | - C H Benton
- Animal & Plant Health Agency, United Kingdom
| | - R J Delahay
- Animal & Plant Health Agency, United Kingdom
| | - G Marion
- Biomathematics and Statistics Scotland, United Kingdom
| | | | - C M Pooley
- Biomathematics and Statistics Scotland, United Kingdom
| | - G C Smith
- Animal & Plant Health Agency, United Kingdom
| | - R R Kao
- Roslin Institute, University of Edinburgh, United Kingdom; Royal (Dick) School of Veterinary Studies, University of Edinburgh, United Kingdom.
| |
Collapse
|
3
|
Wang L, Lim DR, Phelan JE, Reyes LT, Palparan AG, Sanchez MGC, Asuncion LAA, Tujan MAA, Medado IAP, Daga CMA, Llames JHS, Mitarai S, Murase Y, Morishige Y, Ang CF, Malabad JCM, Montoya JC, Hafalla JC, Campino S, Hibberd ML, Ama CG, Basilio RP, La Paz EMCD, Clark TG. Whole genome sequencing analysis of Mycobacterium tuberculosis reveals circulating strain types and drug-resistance mutations in the Philippines. Sci Rep 2024; 14:19602. [PMID: 39179783 PMCID: PMC11344074 DOI: 10.1038/s41598-024-70471-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2024] [Accepted: 08/16/2024] [Indexed: 08/26/2024] Open
Abstract
The Philippines is a high-incidence country for tuberculosis, with the increasing prevalence of multi- (MDR-TB) and extensively-drug (XDR-TB) resistant Mycobacterium tuberculosis strains posing difficulties to disease control. Understanding the genetic diversity of circulating strains can provide insights into underlying drug resistance mutations and transmission dynamics, thereby assisting the design of diagnostic tools, including those using next generation sequencing (NGS) platforms. By analysing genome sequencing data of 732 isolates from Philippines drug-resistance survey collections spanning from 2011 to 2019, we found that the majority belonged to lineages L1 (531/732; 72.5%) and L4 (European-American; n = 174; 23.8%), with the Manila strain (L1.2.1.2.1) being the most prominent (475/531). Approximately two-thirds of isolates were found to be at least MDR-TB (483/732; 66.0%), and potential XDR-TB genotypic resistance was observed (3/732; 0.4%), highlighting an emerging problem in the country. Genotypic resistance was highly concordant with laboratory drug susceptibility testing. By finding isolates with (near-)identical genomic variation, five major clusters containing a total of 114 isolates were identified: all containing either L1 or L4 isolates with at least MDR-TB resistance and spanning multiple years of collection. Closer inspection of clusters revealed transmission in prisons, some involving isolates with XDR-TB, and mutations linked to third-line drug bedaquiline. We have also identified previously unreported mutations linked to resistance for isoniazid, rifampicin, ethambutol, and fluoroquinolones. Overall, this study provides important insights into the genetic diversity, transmission and circulating drug resistance mutations of M. tuberculosis in the Philippines, thereby informing clinical and surveillance decision-making, which is increasingly using NGS platforms.
Collapse
Affiliation(s)
- Linfeng Wang
- Department of Infection Biology, Faculty of Infectious and Tropical Diseases, London School of Hygiene & Tropical Medicine, London, WC1E 7HT, UK
| | - Dodge R Lim
- National Tuberculosis Reference Laboratory, Research Institute for Tropical Medicine, Department of Health, Filinvest, Alabang, Muntinlupa City, Philippines
| | - Jody E Phelan
- Department of Infection Biology, Faculty of Infectious and Tropical Diseases, London School of Hygiene & Tropical Medicine, London, WC1E 7HT, UK
| | - Lorenzo T Reyes
- National Tuberculosis Reference Laboratory, Research Institute for Tropical Medicine, Department of Health, Filinvest, Alabang, Muntinlupa City, Philippines
| | - Alma G Palparan
- National Tuberculosis Reference Laboratory, Research Institute for Tropical Medicine, Department of Health, Filinvest, Alabang, Muntinlupa City, Philippines
| | - Maria Guileane C Sanchez
- National Tuberculosis Reference Laboratory, Research Institute for Tropical Medicine, Department of Health, Filinvest, Alabang, Muntinlupa City, Philippines
| | - Louella Abigail A Asuncion
- National Tuberculosis Reference Laboratory, Research Institute for Tropical Medicine, Department of Health, Filinvest, Alabang, Muntinlupa City, Philippines
| | - Ma Angelica A Tujan
- National Tuberculosis Reference Laboratory, Research Institute for Tropical Medicine, Department of Health, Filinvest, Alabang, Muntinlupa City, Philippines
| | - Inez Andrea P Medado
- Molecular Biology Laboratory, Research Institute for Tropical Medicine, Department of Health, Filinvest, Alabang, Muntinlupa City, Philippines
| | - Chona Mae A Daga
- Department of Epidemiology and Biostatistics, Research Institute for Tropical Medicine, Department of Health, Filinvest, Alabang, Muntinlupa City, Philippines
| | - Jo-Hannah S Llames
- Philippine Genome Center, University of Philippines Diliman, Quezon City, Philippines
| | - Satoshi Mitarai
- Department of Mycobacterium Reference and Research, Research Institute of Tuberculosis, Japan Anti-Tuberculosis Association, Tokyo, Japan
| | - Yoshiro Murase
- Department of Mycobacterium Reference and Research, Research Institute of Tuberculosis, Japan Anti-Tuberculosis Association, Tokyo, Japan
| | - Yuta Morishige
- Department of Mycobacterium Reference and Research, Research Institute of Tuberculosis, Japan Anti-Tuberculosis Association, Tokyo, Japan
| | - Concepcion F Ang
- College of Medicine, University of the Philippines, Manila, Philippines
| | - John Carlo M Malabad
- Department of Science and Technology, DOST Compound, Bicutan, Taguig City, Philippine, Philippines
| | - Jaime C Montoya
- College of Medicine, University of the Philippines, Manila, Philippines
| | - Julius C Hafalla
- Department of Infection Biology, Faculty of Infectious and Tropical Diseases, London School of Hygiene & Tropical Medicine, London, WC1E 7HT, UK
| | - Susana Campino
- Department of Infection Biology, Faculty of Infectious and Tropical Diseases, London School of Hygiene & Tropical Medicine, London, WC1E 7HT, UK
| | - Martin L Hibberd
- Department of Infection Biology, Faculty of Infectious and Tropical Diseases, London School of Hygiene & Tropical Medicine, London, WC1E 7HT, UK
| | - Cecilia G Ama
- National Tuberculosis Reference Laboratory, Research Institute for Tropical Medicine, Department of Health, Filinvest, Alabang, Muntinlupa City, Philippines
| | - Ramon P Basilio
- National Tuberculosis Reference Laboratory, Research Institute for Tropical Medicine, Department of Health, Filinvest, Alabang, Muntinlupa City, Philippines
| | - Eva Maria Cutiongoco-De La Paz
- Philippine Genome Center, University of Philippines Diliman, Quezon City, Philippines
- College of Medicine, University of the Philippines, Manila, Philippines
| | - Taane G Clark
- Department of Infection Biology, Faculty of Infectious and Tropical Diseases, London School of Hygiene & Tropical Medicine, London, WC1E 7HT, UK.
- Department of Infectious Disease Epidemiology, Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London, WC1E 7HT, UK.
| |
Collapse
|
4
|
Weimann A, Dinan AM, Ruis C, Bernut A, Pont S, Brown K, Ryan J, Santos L, Ellison L, Ukor E, Pandurangan AP, Krokowski S, Blundell TL, Welch M, Blane B, Judge K, Bousfield R, Brown N, Bryant JM, Kukavica-Ibrulj I, Rampioni G, Leoni L, Harrison PT, Peacock SJ, Thomson NR, Gauthier J, Fothergill JL, Levesque RC, Parkhill J, Floto RA. Evolution and host-specific adaptation of Pseudomonas aeruginosa. Science 2024; 385:eadi0908. [PMID: 38963857 DOI: 10.1126/science.adi0908] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2023] [Accepted: 05/02/2024] [Indexed: 07/06/2024]
Abstract
The major human bacterial pathogen Pseudomonas aeruginosa causes multidrug-resistant infections in people with underlying immunodeficiencies or structural lung diseases such as cystic fibrosis (CF). We show that a few environmental isolates, driven by horizontal gene acquisition, have become dominant epidemic clones that have sequentially emerged and spread through global transmission networks over the past 200 years. These clones demonstrate varying intrinsic propensities for infecting CF or non-CF individuals (linked to specific transcriptional changes enabling survival within macrophages); have undergone multiple rounds of convergent, host-specific adaptation; and have eventually lost their ability to transmit between different patient groups. Our findings thus explain the pathogenic evolution of P. aeruginosa and highlight the importance of global surveillance and cross-infection prevention in averting the emergence of future epidemic clones.
Collapse
Affiliation(s)
- Aaron Weimann
- Victor Phillip Dahdaleh Heart & Lung Research Institute, University of Cambridge, Cambridge, UK
- University of Cambridge Molecular Immunity Unit, MRC Laboratory of Molecular Biology, Cambridge, UK
- Cambridge Centre for AI in Medicine, University of Cambridge, Cambridge, UK
- Department of Veterinary Medicine, University of Cambridge, Cambridge, UK
| | - Adam M Dinan
- Victor Phillip Dahdaleh Heart & Lung Research Institute, University of Cambridge, Cambridge, UK
- University of Cambridge Molecular Immunity Unit, MRC Laboratory of Molecular Biology, Cambridge, UK
- Cambridge Centre for AI in Medicine, University of Cambridge, Cambridge, UK
| | - Christopher Ruis
- Victor Phillip Dahdaleh Heart & Lung Research Institute, University of Cambridge, Cambridge, UK
- University of Cambridge Molecular Immunity Unit, MRC Laboratory of Molecular Biology, Cambridge, UK
- Cambridge Centre for AI in Medicine, University of Cambridge, Cambridge, UK
- Department of Veterinary Medicine, University of Cambridge, Cambridge, UK
| | - Audrey Bernut
- Laboratory of Pathogens and Host Immunity (LPHI), UMR5235, CNRS/Université de Montpellier, Montpellier, France
| | - Stéphane Pont
- Laboratory of Pathogens and Host Immunity (LPHI), UMR5235, CNRS/Université de Montpellier, Montpellier, France
| | - Karen Brown
- Victor Phillip Dahdaleh Heart & Lung Research Institute, University of Cambridge, Cambridge, UK
- University of Cambridge Molecular Immunity Unit, MRC Laboratory of Molecular Biology, Cambridge, UK
- Cambridge Centre for Lung Infection, Royal Papworth Hospital, Cambridge, UK
| | - Judy Ryan
- Victor Phillip Dahdaleh Heart & Lung Research Institute, University of Cambridge, Cambridge, UK
- University of Cambridge Molecular Immunity Unit, MRC Laboratory of Molecular Biology, Cambridge, UK
| | - Lúcia Santos
- Department of Physiology, Bioscience Institute, University College Cork, Cork, Ireland
| | - Louise Ellison
- University of Cambridge Molecular Immunity Unit, MRC Laboratory of Molecular Biology, Cambridge, UK
| | - Emem Ukor
- University of Cambridge Molecular Immunity Unit, MRC Laboratory of Molecular Biology, Cambridge, UK
- Cambridge Centre for Lung Infection, Royal Papworth Hospital, Cambridge, UK
| | - Arun P Pandurangan
- Victor Phillip Dahdaleh Heart & Lung Research Institute, University of Cambridge, Cambridge, UK
- Department of Biochemistry, University of Cambridge, Cambridge, UK
- Department of Medicine, University of Cambridge, Cambridge, UK
| | - Sina Krokowski
- Victor Phillip Dahdaleh Heart & Lung Research Institute, University of Cambridge, Cambridge, UK
- University of Cambridge Molecular Immunity Unit, MRC Laboratory of Molecular Biology, Cambridge, UK
| | - Tom L Blundell
- Victor Phillip Dahdaleh Heart & Lung Research Institute, University of Cambridge, Cambridge, UK
- Department of Biochemistry, University of Cambridge, Cambridge, UK
- Department of Medicine, University of Cambridge, Cambridge, UK
| | - Martin Welch
- Department of Biochemistry, University of Cambridge, Cambridge, UK
| | - Beth Blane
- Department of Medicine, University of Cambridge, Cambridge, UK
| | - Kim Judge
- Wellcome Sanger Institute, Hinxton, UK
| | - Rachel Bousfield
- Department of Medicine, University of Cambridge, Cambridge, UK
- Cambridge University Hospitals Trust, Cambridge, UK
| | | | | | - Irena Kukavica-Ibrulj
- Institut de Biologie Intégrative et des Systèmes (IBIS), Université Laval, Québec City, Québec, Canada
| | - Giordano Rampioni
- Department of Science, University Roma Tre, Rome, Italy
- IRCCS Fondazione Santa Lucia, Rome, Italy
| | - Livia Leoni
- Department of Science, University Roma Tre, Rome, Italy
| | - Patrick T Harrison
- Department of Physiology, Bioscience Institute, University College Cork, Cork, Ireland
| | - Sharon J Peacock
- Department of Medicine, University of Cambridge, Cambridge, UK
- Cambridge University Hospitals Trust, Cambridge, UK
| | - Nicholas R Thomson
- Wellcome Sanger Institute, Hinxton, UK
- London School of Hygiene and Tropical Medicine, London, UK
| | - Jeff Gauthier
- Institut de Biologie Intégrative et des Systèmes (IBIS), Université Laval, Québec City, Québec, Canada
| | - Jo L Fothergill
- Department of Clinical Infection, Microbiology and Immunology, University of Liverpool, Liverpool, UK
| | - Roger C Levesque
- Institut de Biologie Intégrative et des Systèmes (IBIS), Université Laval, Québec City, Québec, Canada
| | - Julian Parkhill
- Department of Veterinary Medicine, University of Cambridge, Cambridge, UK
| | - R Andres Floto
- Victor Phillip Dahdaleh Heart & Lung Research Institute, University of Cambridge, Cambridge, UK
- University of Cambridge Molecular Immunity Unit, MRC Laboratory of Molecular Biology, Cambridge, UK
- Cambridge Centre for AI in Medicine, University of Cambridge, Cambridge, UK
- Cambridge Centre for Lung Infection, Royal Papworth Hospital, Cambridge, UK
- Department of Medicine, University of Cambridge, Cambridge, UK
- Cambridge University Hospitals Trust, Cambridge, UK
| |
Collapse
|
5
|
Pereira AC, Lourenço J, Themudo G, Botelho A, Cunha MV. Population structure and history of Mycobacterium bovis European 3 clonal complex reveal transmission across ecological corridors of unrecognized importance in Portugal. Microbiol Spectr 2024; 12:e0382923. [PMID: 38771094 PMCID: PMC11218495 DOI: 10.1128/spectrum.03829-23] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2023] [Accepted: 04/22/2024] [Indexed: 05/22/2024] Open
Abstract
Mycobacterium bovis causes animal tuberculosis in livestock and wildlife, with an impact on animal health and production, wildlife management, and public health. In this work, we sampled a multi-host tuberculosis community from the official hotspot risk area of Portugal over 16 years, generating the largest available data set in the country. Using phylogenetic and ecological modeling, we aimed to reconstruct the history of circulating lineages across the livestock-wildlife interface to inform intervention and the implementation of genomic surveillance within the official eradication plan. We find evidence for the co-circulation of M. bovis European 1 (Eu1), Eu2, and Eu3 clonal complexes, with Eu3 providing sufficient temporal signal for further phylogenetic investigation. The Eu3 most recent common ancestor (bovine) was dated in the 1990s, subsequently transitioning to wildlife (red deer and wild boar). Isolate clustering based on sample metadata was used to inform phylogenetic inference, unravelng frequent transmission between two clusters that represent an ecological corridor of previously unrecognized importance in Portugal. The latter was associated with transmission at the livestock-wildlife interface toward locations with higher temperature and precipitation, lower agriculture and road density, and lower host densities. This is the first analysis of M. bovis Eu3 complex in Iberia, shedding light on background ecological factors underlying long-term transmission and informing where efforts could be focused within the larger hotspot risk area of Portugal. IMPORTANCE Efforts to strengthen surveillance and control of animal tuberculosis (TB) are ongoing worlwide. Here, we developed an eco-phylodynamic framework based on discrete phylogenetic approaches informed by M. bovis whole-genome sequence data representing a multi-host transmission system at the livestock-wildlife interface, within a rich ecological landscape in Portugal, to understand transmission processes and translate this knowledge into disease management benefits. We find evidence for the co-circulation of several M. bovis clades, with frequent transmission of the Eu3 lineage among cattle and wildlife populations. Most transition events between different ecological settings took place toward host, climate and land use gradients, underscoring animal TB expansion and a potential corridor of unrecognized importance for M. bovis maintenance. Results stress that animal TB is an established wildlife disease without ecological barriers, showing that control measures in place are insufficient to prevent long-distance transmission and spillover across multi-host communities, demanding new interventions targeting livestock-wildlife interactions.
Collapse
Affiliation(s)
- André C. Pereira
- Centre for Ecology, Evolution and Environmental Changes (cE3c) & CHANGE—Global Change and Sustainability Institute, Faculdade de Ciências, Universidade de Lisboa, Lisbon, Portugal
- Biosystems & Integrative Sciences Institute (BioISI), Faculdade de Ciências, Universidade de Lisboa, Lisbon, Portugal
| | - José Lourenço
- Biosystems & Integrative Sciences Institute (BioISI), Faculdade de Ciências, Universidade de Lisboa, Lisbon, Portugal
| | - Gonçalo Themudo
- Centre for Ecology, Evolution and Environmental Changes (cE3c) & CHANGE—Global Change and Sustainability Institute, Faculdade de Ciências, Universidade de Lisboa, Lisbon, Portugal
| | - Ana Botelho
- INIAV, I.P.—National Institute for Agrarian and Veterinary Research, Oeiras, Portugal
| | - Mónica V. Cunha
- Centre for Ecology, Evolution and Environmental Changes (cE3c) & CHANGE—Global Change and Sustainability Institute, Faculdade de Ciências, Universidade de Lisboa, Lisbon, Portugal
- Biosystems & Integrative Sciences Institute (BioISI), Faculdade de Ciências, Universidade de Lisboa, Lisbon, Portugal
| |
Collapse
|
6
|
Liu Z, Li X, Xiong H, Zhou Q, Yi H, Wu K, Zhou Z, Lu Y, Zhu Y, Zhou L, Zhang M, Gao J, Gao J, Chen S, Wang X, Wang W. Genomic and spatial analysis reveal the transmission dynamics of tuberculosis in areas with high incidence of Zhejiang, China: A prospective cohort study. INFECTION, GENETICS AND EVOLUTION : JOURNAL OF MOLECULAR EPIDEMIOLOGY AND EVOLUTIONARY GENETICS IN INFECTIOUS DISEASES 2024; 121:105603. [PMID: 38723983 DOI: 10.1016/j.meegid.2024.105603] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/01/2023] [Revised: 03/28/2024] [Accepted: 05/05/2024] [Indexed: 05/14/2024]
Abstract
In the mountainous, rural regions of eastern China, tuberculosis (TB) remains a formidable challenge; however, the long-term molecular epidemiological surveillance in these regions is limited. This study aimed to investigate molecular and spatial epidemiology of TB in two mountainous, rural counties of Zhejiang Province, China, from 2015 to 2021, to elucidate the recent transmission and drug-resistance profiles. The predominant Lineage 2 (L2) Beijing family accounted for 80.1% of total 532 sequenced Mycobacterium tuberculosis (Mtb) strains, showing consistent prevalence over seven years. Gene mutations associated with drug resistance were identified in 19.4% (103/532) of strains, including 47 rifampicin or isoniazid-resistant strains, eight multi-drug-resistant (MDR) strains, and five pre-extensively drug-resistant (pre-XDR) strains. Genomic clustering revealed 53 distinct clusters with an overall transmission clustering rate of 23.9% (127/532). Patients with a history of retreatment and those infected with L2 strains had a higher risk of recent transmission. Spatial and epidemiological analysis unveiled significant transmission hotspots, especially in densely populated urban areas, involving various public places such as medical institutions, farmlands, markets, and cardrooms. The study emphasizes the pivotal role of Beijing strains and urban-based TB transmission in the western mountainous regions in Zhejiang, highlighting the urgent requirement for specific interventions to mitigate the impact of TB in these unique communities.
Collapse
Affiliation(s)
- Zhengwei Liu
- The Institute of TB Control, Zhejiang Provincial Center for Disease Control and Prevention, Hangzhou, Zhejiang 310051, China
| | - Xiangchen Li
- Key Laboratory of Precision Medicine in Diagnosis and Monitoring Research of Zhejiang Province, Hangzhou, Zhejiang 310020, China
| | - Haiyan Xiong
- School of Public Health, Fudan University, Shanghai 200032, China; Key Laboratory of Public Health Safety of Ministry of Education, Fudan University, Shanghai 200032, China
| | - Qingrong Zhou
- Center for Disease Control and Prevention of Jiangshan City, Zhejiang 324100, China
| | - Huaiming Yi
- Center for Disease Control and Prevention of Changshan County, Zhejiang 324200, China
| | - Kunyang Wu
- The Institute of TB Control, Zhejiang Provincial Center for Disease Control and Prevention, Hangzhou, Zhejiang 310051, China
| | - Zonglei Zhou
- School of Public Health, Fudan University, Shanghai 200032, China
| | - Yewei Lu
- Key Laboratory of Precision Medicine in Diagnosis and Monitoring Research of Zhejiang Province, Hangzhou, Zhejiang 310020, China
| | - Yelei Zhu
- The Institute of TB Control, Zhejiang Provincial Center for Disease Control and Prevention, Hangzhou, Zhejiang 310051, China
| | - Lin Zhou
- The Institute of TB Control, Zhejiang Provincial Center for Disease Control and Prevention, Hangzhou, Zhejiang 310051, China
| | - Mingwu Zhang
- The Institute of TB Control, Zhejiang Provincial Center for Disease Control and Prevention, Hangzhou, Zhejiang 310051, China
| | - Junshun Gao
- Key Laboratory of Precision Medicine in Diagnosis and Monitoring Research of Zhejiang Province, Hangzhou, Zhejiang 310020, China
| | - Junli Gao
- Key Laboratory of Precision Medicine in Diagnosis and Monitoring Research of Zhejiang Province, Hangzhou, Zhejiang 310020, China
| | - Songhua Chen
- The Institute of TB Control, Zhejiang Provincial Center for Disease Control and Prevention, Hangzhou, Zhejiang 310051, China.
| | - Xiaomeng Wang
- The Institute of TB Control, Zhejiang Provincial Center for Disease Control and Prevention, Hangzhou, Zhejiang 310051, China.
| | - Weibing Wang
- School of Public Health, Fudan University, Shanghai 200032, China; Key Laboratory of Public Health Safety of Ministry of Education, Fudan University, Shanghai 200032, China.
| |
Collapse
|
7
|
Austin RM, Honap TP, Mann AE, Hübner A, DeGaglia CMS, Warinner C, Zuckerman MK, Hofman CA. Metagenomic and paleopathological analyses of a historic documented collection explore ancient dental calculus as a diagnostic tool. Sci Rep 2024; 14:14720. [PMID: 38926415 PMCID: PMC11208530 DOI: 10.1038/s41598-024-64818-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2024] [Accepted: 06/13/2024] [Indexed: 06/28/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Rita M Austin
- Frontiers in Evolutionary Zoology Research Group, Natural History Museum of Oslo, University of Oslo, Oslo, 0562, Norway.
- Department of Anthropology, National Museum of Natural History, Smithsonian Institution, Washington, DC, 20560, USA.
- Department of Anthropology, University of Oklahoma, Norman, OK, 73019, USA.
- Laboratories of Molecular Anthropology and Microbiome Research, University of Oklahoma, Norman, OK, 73019, USA.
| | - Tanvi P Honap
- Department of Anthropology, University of Oklahoma, Norman, OK, 73019, USA
- Laboratories of Molecular Anthropology and Microbiome Research, University of Oklahoma, Norman, OK, 73019, USA
| | - Allison E Mann
- Department of Biological Sciences, Clemson University, Clemson, SC, 29634, USA
| | - Alexander Hübner
- Department Archaeogenetics, Max-Planck-Institute for Evolutionary Anthropology, Leipzig, 04103, Germany
| | | | - Christina Warinner
- Department of Anthropology, Harvard University, Cambridge, MA, 02138, USA
| | - Molly K Zuckerman
- Department of Anthropology and Middle Eastern Cultures, Mississippi State University, Mississippi State, MS, 39762, USA.
| | - Courtney A Hofman
- Department of Anthropology, National Museum of Natural History, Smithsonian Institution, Washington, DC, 20560, USA.
- Department of Anthropology, University of Oklahoma, Norman, OK, 73019, USA.
- Laboratories of Molecular Anthropology and Microbiome Research, University of Oklahoma, Norman, OK, 73019, USA.
| |
Collapse
|
8
|
Bainomugisa A, Pandey S, O'Connor B, Syrmis M, Whiley D, Sintchenko V, Coin LJ, Marais BJ, Coulter C. Sustained transmission over two decades of a previously unrecognised MPT64 negative Mycobacterium tuberculosis strain in Queensland, Australia: a whole genome sequencing study. THE LANCET REGIONAL HEALTH. WESTERN PACIFIC 2024; 47:101105. [PMID: 39022748 PMCID: PMC11253042 DOI: 10.1016/j.lanwpc.2024.101105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/09/2024] [Revised: 04/09/2024] [Accepted: 05/16/2024] [Indexed: 07/20/2024]
Abstract
Background MPT64 is a key protein used for Mycobacterium tuberculosis (MTB) complex strain identification. We describe protracted transmission of an MPT64 negative MTB strain in Queensland, Australia, and explore genomic factors related to its successful spread. Methods All MPT64 negative strains identified between 2002 and 2022 by the Queensland Mycobacteria Reference Laboratory, and an additional 2 isolates from New South Wales (NSW), were whole genome sequenced. Bayesian modelling and phylogeographical analyses were used to assess their evolutionary history and transmission dynamics. Protein structural modelling to understand the putative functional effects of the mutated gene coding for MPT64 protein was performed. Findings Forty-three MPT64 negative isolates were sequenced, belonging to a single MTB cluster of Lineage 4.1.1.1 strains. Combined with a UK dataset of the same lineage, molecular dating estimated 1990 (95% HPD 1987-1993) as the likely time of strain introduction into Australia. Although the strain has spread over a wide geographic area and new cases linked to the cluster continue to arise, phylodynamic analysis suggest the outbreak peaked around 2003. All MPT64 negative strains had a frame shift mutation (delAT, p.Val216fs) within the MPT64 gene, which confers two major structural rearrangements at the C-terminus of the protein. Interpretation This study uncovered the origins of an MPT64 negative MTB outbreak in Australia, providing a richer understanding of its biology and transmission dynamics, as well as guidance for clinical diagnosis and public health action. The potential spread of MPT64 negative strains undermines the diagnostic utility of the MPT64 immunochromatographic test. Funding This study was funded from an operational budget provided to the Queensland Mycobacterium Reference Laboratory by Pathology Queensland, Queensland Department of Health.
Collapse
Affiliation(s)
- Arnold Bainomugisa
- Queensland Mycobacterium Reference Laboratory, Brisbane, Queensland, Australia
| | - Sushil Pandey
- Queensland Mycobacterium Reference Laboratory, Brisbane, Queensland, Australia
| | - Bridget O'Connor
- Public Health Intelligence Branch, Department of Health, Brisbane, Queensland, Australia
| | - Melanie Syrmis
- Queensland Mycobacterium Reference Laboratory, Brisbane, Queensland, Australia
| | - David Whiley
- University of Queensland Centre for Clinical Research, Brisbane, Queensland, Australia
| | - Vitali Sintchenko
- Centre for Infectious Diseases and Microbiology-Public Health, Westmead Hospital, Western Sydney Local Health District, Sydney, New South Wales, Australia
- NSW Mycobacterium Reference Laboratory, Centre for Infectious Diseases and Microbiology Laboratory Services, Institute of Clinical Pathology and Medical Research, NSW Health Pathology - Western, Sydney, New South Wales, Australia
- Sydney Infectious Diseases Institute (Sydney ID), The University of Sydney, Sydney, New South Wales, Australia
| | - Lachlan J.M. Coin
- Department of Microbiology and Immunology, University of Melbourne, Peter Doherty Institute for Infection and Immunity, Melbourne, Victoria 3000, Australia
| | - Ben J. Marais
- Centre for Infectious Diseases and Microbiology-Public Health, Westmead Hospital, Western Sydney Local Health District, Sydney, New South Wales, Australia
- Sydney Infectious Diseases Institute (Sydney ID), The University of Sydney, Sydney, New South Wales, Australia
| | - Christopher Coulter
- Queensland Mycobacterium Reference Laboratory, Brisbane, Queensland, Australia
- Communicable Diseases Branch, Department of Health, Brisbane, Queensland, Australia
| |
Collapse
|
9
|
Dixit A, Ektefaie Y, Kagal A, Freschi L, Karyakarte R, Lokhande R, Groschel M, Tornheim JA, Gupte N, Pradhan NN, Paradkar MS, Deshmukh S, Kadam D, Schito M, Engelthaler DM, Gupta A, Golub J, Mave V, Farhat M. Drug resistance and epidemiological success of modern Mycobacterium tuberculosis lineages in western India. J Infect Dis 2024:jiae240. [PMID: 38819323 DOI: 10.1093/infdis/jiae240] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2023] [Revised: 04/24/2024] [Accepted: 05/03/2024] [Indexed: 06/01/2024] Open
Abstract
BACKGROUND Transmission is contributing to the slow decline of tuberculosis (TB) incidence globally. Drivers of TB transmission in India, the country estimated to carry a quarter of the World's burden, are not well studied. We conducted a genomic epidemiology study to compare epidemiological success, host factors and drug resistance (DR) among the four major Mycobacterium tuberculosis (Mtb) lineages (L1-4) circulating in Pune, India. METHODS We performed whole-genome sequencing (WGS) of Mtb sputum culture-positive isolates from participants in two prospective cohort studies and predicted genotypic susceptibility using a validated random forest model. We used maximum likelihood estimation to build phylogenies. We compared lineage specific phylogenetic and time-scaled metrics to assess epidemiological success. RESULTS Of the 642 isolates that underwent WGS, 612 met sequence quality criteria. Most isolates belonged to L3 (44.6%). The majority (61.1%) of multidrug-resistant isolates belonged to L2 (P < 0.001). In molecular dating, L2 demonstrated a higher rate and more recent resistance acquisition. We measured higher clustering, and time-scaled haplotypic density (THD) for L4 and L2 compared to L3 and/or L1 suggesting higher epidemiological success. L4 demonstrated higher THD and clustering (OR 5.1 (95% CI 2.3-12.3) in multivariate models controlling for host factors and DR. CONCLUSION L2 shows a higher frequency of DR and both L2 and L4 demonstrate evidence of higher epidemiological success than L3 or L1 in the study setting. Our findings highlight the need for contact tracing around TB cases, and heightened surveillance of TB DR in India.
Collapse
Affiliation(s)
- Avika Dixit
- Division of Infectious Diseases, Boston Children's Hospital, Boston MA, USA
- Department of Biomedical Informatics, Harvard Medical School, Boston MA, USA
| | - Yasha Ektefaie
- Department of Biomedical Informatics, Harvard Medical School, Boston MA, USA
| | - Anju Kagal
- Byramjee-Jeejeebhoy Government Medical College, Pune, India
| | - Luca Freschi
- Department of Biomedical Informatics, Harvard Medical School, Boston MA, USA
| | | | - Rahul Lokhande
- Byramjee-Jeejeebhoy Government Medical College, Pune, India
| | - Matthias Groschel
- Department of Biomedical Informatics, Harvard Medical School, Boston MA, USA
| | - Jeffrey A Tornheim
- Center for Clinical Global Health Education, Division of Infectious Diseases, Department of Medicine, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Center for Tuberculosis Research, Department of Medicine, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Nikhil Gupte
- Center for Clinical Global Health Education, Division of Infectious Diseases, Department of Medicine, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Byramjee-Jeejeebhoy Medical College-Johns Hopkins University Clinical Research Site, Pune, India
- Johns Hopkins India, Pune, India
| | - Neeta N Pradhan
- Byramjee-Jeejeebhoy Medical College-Johns Hopkins University Clinical Research Site, Pune, India
- Johns Hopkins India, Pune, India
| | - Mandar S Paradkar
- Byramjee-Jeejeebhoy Medical College-Johns Hopkins University Clinical Research Site, Pune, India
- Johns Hopkins India, Pune, India
| | - Sona Deshmukh
- Byramjee-Jeejeebhoy Medical College-Johns Hopkins University Clinical Research Site, Pune, India
- Johns Hopkins India, Pune, India
| | - Dileep Kadam
- Byramjee-Jeejeebhoy Government Medical College, Pune, India
| | | | | | - Amita Gupta
- Center for Clinical Global Health Education, Division of Infectious Diseases, Department of Medicine, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Department of International Health, Johns Hopkins Bloomberg School of Public Heath, Baltimore, MD, USA
| | - Jonathan Golub
- Center for Tuberculosis Research, Department of Medicine, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Vidya Mave
- Center for Clinical Global Health Education, Division of Infectious Diseases, Department of Medicine, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Byramjee-Jeejeebhoy Medical College-Johns Hopkins University Clinical Research Site, Pune, India
- Johns Hopkins India, Pune, India
| | - Maha Farhat
- Department of Biomedical Informatics, Harvard Medical School, Boston MA, USA
- Division of Pulmonary and Critical Care Medicine, Massachusetts General Hospital, Boston, MA, USA
| |
Collapse
|
10
|
Sheinman M, Arndt PF, Massip F. Modeling the mosaic structure of bacterial genomes to infer their evolutionary history. Proc Natl Acad Sci U S A 2024; 121:e2313367121. [PMID: 38517978 PMCID: PMC10990148 DOI: 10.1073/pnas.2313367121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2023] [Accepted: 01/30/2024] [Indexed: 03/24/2024] Open
Abstract
The chronology and phylogeny of bacterial evolution are difficult to reconstruct due to a scarce fossil record. The analysis of bacterial genomes remains challenging because of large sequence divergence, the plasticity of bacterial genomes due to frequent gene loss, horizontal gene transfer, and differences in selective pressure from one locus to another. Therefore, taking advantage of the rich and rapidly accumulating genomic data requires accurate modeling of genome evolution. An important technical consideration is that loci with high effective mutation rates may diverge beyond the detection limit of the alignment algorithms used, biasing the genome-wide divergence estimates toward smaller divergences. In this article, we propose a novel method to gain insight into bacterial evolution based on statistical properties of genome comparisons. We find that the length distribution of sequence matches is shaped by the effective mutation rates of different loci, by the horizontal transfers, and by the aligner sensitivity. Based on these inputs, we build a model and show that it accounts for the empirically observed distributions, taking the Enterobacteriaceae family as an example. Our method allows to distinguish segments of vertical and horizontal origins and to estimate the time divergence and exchange rate between any pair of taxa from genome-wide alignments. Based on the estimated time divergences, we construct a time-calibrated phylogenetic tree to demonstrate the accuracy of the method.
Collapse
Affiliation(s)
- Michael Sheinman
- Institute for Advanced Studies, Sevastopol State University, Sevastopol299053, Crimea
| | - Peter F. Arndt
- Department of Computational Molecular Biology, Max Planck Institute for Molecular Genetics, Berlin12163, Germany
| | - Florian Massip
- Department U900, Centre for Computational Biology, Mines Paris, PSL University, Paris75006, France
- Department U900, Institut Curie, Université Paris Sciences et Lettres, Paris75005, France
- INSERM, U900, Paris75005, France
| |
Collapse
|
11
|
Nimmo C, Ortiz AT, Tan CCS, Pang J, Acman M, Millard J, Padayatchi N, Grant AD, O'Donnell M, Pym A, Brynildsrud OB, Eldholm V, Grandjean L, Didelot X, Balloux F, van Dorp L. Detection of a historic reservoir of bedaquiline/clofazimine resistance-associated variants in Mycobacterium tuberculosis. Genome Med 2024; 16:34. [PMID: 38374151 PMCID: PMC10877763 DOI: 10.1186/s13073-024-01289-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2023] [Accepted: 01/19/2024] [Indexed: 02/21/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Camus Nimmo
- UCL Genetics Institute, University College London, Darwin Building, Gower Street, London, UK.
- Division of Infection and Immunity, University College London, London, UK.
- Africa Health Research Institute, Durban, South Africa.
| | - Arturo Torres Ortiz
- UCL Genetics Institute, University College London, Darwin Building, Gower Street, London, UK
- Department of Medicine, Imperial College, London, UK
| | - Cedric C S Tan
- UCL Genetics Institute, University College London, Darwin Building, Gower Street, London, UK
| | - Juanita Pang
- UCL Genetics Institute, University College London, Darwin Building, Gower Street, London, UK
- Division of Infection and Immunity, University College London, London, UK
| | - Mislav Acman
- UCL Genetics Institute, University College London, Darwin Building, Gower Street, London, UK
| | - James Millard
- Africa Health Research Institute, Durban, South Africa
- Wellcome Trust Liverpool Glasgow Centre for Global Health Research, Liverpool, UK
- Institute of Infection and Global Health, University of Liverpool, Liverpool, UK
| | - Nesri Padayatchi
- CAPRISA MRC-HIV-TB Pathogenesis and Treatment Research Unit, Durban, South Africa
| | - Alison D Grant
- Africa Health Research Institute, Durban, South Africa
- TB Centre, London School of Hygiene & Tropical Medicine, London, UK
| | - Max O'Donnell
- CAPRISA MRC-HIV-TB Pathogenesis and Treatment Research Unit, Durban, South Africa
- Department of Medicine & Epidemiology, Columbia University Irving Medical Center, New York, NY, USA
| | - Alex Pym
- Africa Health Research Institute, Durban, South Africa
| | - Ola B Brynildsrud
- Division of Infectious Diseases and Environmental Health, Norwegian Institute of Public Health, Oslo, Norway
| | - Vegard Eldholm
- Division of Infectious Diseases and Environmental Health, Norwegian Institute of Public Health, Oslo, Norway
| | - Louis Grandjean
- Division of Infection and Immunity, University College London, London, UK
- Laboratorio de Investigacion y Enfermedades Infecciosas, Universidad Peruana Cayetano Heredia, Lima, Peru
- Department of Infection, Immunity and Inflammation, Institute of Child Health, University College London, London, UK
| | - Xavier Didelot
- School of Life Sciences and Department of Statistics, University of Warwick, Coventry, UK
| | - François Balloux
- UCL Genetics Institute, University College London, Darwin Building, Gower Street, London, UK.
| | - Lucy van Dorp
- UCL Genetics Institute, University College London, Darwin Building, Gower Street, London, UK.
| |
Collapse
|
12
|
Torres Ortiz A, Grandjean L. Phylogenetic Survival Analysis. Methods Mol Biol 2024; 2833:121-128. [PMID: 38949706 DOI: 10.1007/978-1-0716-3981-8_12] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/02/2024]
Abstract
Going back in time through a phylogenetic tree makes it possible to evaluate ancestral genomes and assess their potential to acquire key polymorphisms of interest over evolutionary time. Knowledge of this kind may allow for the emergence of key traits to be predicted and pre-empted from currently circulating strains in the future. Here, we present a novel genome-wide survival analysis and use the emergence of drug resistance in Mycobacterium tuberculosis as an example to demonstrate the potential and utility of the technique.
Collapse
Affiliation(s)
- Arturo Torres Ortiz
- Department of Infection, Immunity and Inflammation, Institute of Child Health, University College London, London, UK
| | - Louis Grandjean
- Department of Infection, Immunity and Inflammation, Institute of Child Health, University College London, London, UK.
| |
Collapse
|
13
|
Zein-Eddine R, Hak F, Le Meur A, Genestet C, Dumitrescu O, Guyeux C, Senelle G, Sola C, Refrégier G. The paradoxes of Mycobacterium tuberculosis molecular evolution and consequences for the inference of tuberculosis emergence date. Tuberculosis (Edinb) 2023; 143S:102378. [PMID: 38012921 DOI: 10.1016/j.tube.2023.102378] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2023] [Revised: 07/10/2023] [Accepted: 07/17/2023] [Indexed: 11/29/2023]
Abstract
The date of Mycobacterium tuberculosis complex emergence has been the subject of long debates. New studies joining archaeological efforts with sequencing methods raise high hopes for solving whether this emergence is closer to 70,000 or to 6000 years before present. Inferring the date of emergence of this pathogen based on sequence data requires a molecular clock. Several clocks inferred from different types of loci and/or different samples, using both sound reasoning and reliable data, are actually very different, which we refer to as the paradoxes of M. tuberculosis molecular evolution. After having presented these paradoxes, we will remind the limits of the molecular clocks used in the different studies such as the assumption of homogeneous substitution rate. We will then review recent results that shed new light on the characteristics of M. tuberculosis molecular evolution: traces of diverse selection pressures, the impact of host dynamics, etc. We provide some ideas on what to do next to get nearer to a reliable dating of Mycobacterium tuberculosis complex emergence. Among them, the collection of additional remains from ancient tuberculosis seems still essential.
Collapse
Affiliation(s)
- R Zein-Eddine
- Laboratoire d'Optique et Biosciences, Ecole Polytechnique, Institut National de la Santé et de la Recherche Médicale: U1182, Centre National de la Recherche Scientifique: UMR7645, France
| | - F Hak
- Université Paris-Saclay, CNRS, AgroParisTech, Ecologie Systématique et Evolution, 91190, Gif-sur-Yvette, France
| | - A Le Meur
- Université Paris-Saclay, CNRS, AgroParisTech, Ecologie Systématique et Evolution, 91190, Gif-sur-Yvette, France
| | - C Genestet
- CIRI - Centre International de Recherche en Infectiologie, Ecole Normale Supérieure de Lyon, Université Claude Bernard Lyon-1, Inserm U1111, CNRS UMR5308, Lyon, France; Hospices Civils de Lyon, Institut des Agents Infectieux, Laboratoire de bactériologie, Lyon, France
| | - O Dumitrescu
- CIRI - Centre International de Recherche en Infectiologie, Ecole Normale Supérieure de Lyon, Université Claude Bernard Lyon-1, Inserm U1111, CNRS UMR5308, Lyon, France; Hospices Civils de Lyon, Institut des Agents Infectieux, Laboratoire de bactériologie, Lyon, France
| | - C Guyeux
- DISC Computer Science Department, FEMTO-ST Institute, UMR 6174 CNRS, Univ. Bourgogne Franche-Comté (UBFC), 16 Route de Gray, 25000, Besançon, France
| | - G Senelle
- DISC Computer Science Department, FEMTO-ST Institute, UMR 6174 CNRS, Univ. Bourgogne Franche-Comté (UBFC), 16 Route de Gray, 25000, Besançon, France
| | - C Sola
- Université de Paris, IAME, UMR1137, INSERM, Paris, France; AP-HP, GHU Nord, Service de mycobactériologie spécialisée et de référence, Paris, France; Université Paris-Saclay, Saint-Aubin, France
| | - G Refrégier
- Université Paris-Saclay, CNRS, AgroParisTech, Ecologie Systématique et Evolution, 91190, Gif-sur-Yvette, France.
| |
Collapse
|
14
|
Pereira AC, Reis AC, Cunha MV. Genomic epidemiology sheds light on the emergence and spread of Mycobacterium bovis Eu2 Clonal Complex in Portugal. Emerg Microbes Infect 2023; 12:2253340. [PMID: 37640285 PMCID: PMC10484045 DOI: 10.1080/22221751.2023.2253340] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2023] [Revised: 08/21/2023] [Accepted: 08/24/2023] [Indexed: 08/31/2023]
Abstract
ABSTRACTAnimal tuberculosis (TB) remains a serious concern for animal and human health. Mycobacterium bovis circulates in multi-host systems, dominated by the European 2 clonal complex (Eu2) in Iberia. In this work, we use genomic epidemiology to infer the emergence, spread, and spatiotemporal patterns of Eu2 in the official epidemiological risk area of animal TB in Portugal. Phylogenetic analysis of 144 M. bovis whole-genome sequences from cattle, wild boar, and red deer, representing the 2002-2021 period, distinguished three Eu2 clades that evolved independently. The major Eu2 clade underwent phylodynamic inferences to estimate the time and location of outbreaks, host transitions, and spatial diffusion as well. The origin of this Eu2 clade was attributed to the red deer population in the Castelo Branco district, near the border with Spain. Most host transitions were intraspecific (80%), while interspecific transmissions between wildlife species (wild boar-red deer), and between wild boar and cattle, were highly supported. Phylogeographic reconstruction evidenced that most transitions (82%) occur within municipalities, highlighting local transmission corridors.Our study indicates that M. bovis continues to spread at the cattle-wildlife interface within the animal TB hotspot area, possibly driven by the foraging behaviour of wild boar near agricultural lands. Red deer seems to be an important driver of TB within wildlife hosts, while the wild boar links the multi-host wildlife community and livestock. This work highlights the value of combining genomic epidemiology with phylodynamic inference to resolve host jumps and spatial patterns of M. bovis, providing real-time clues about points of intervention.
Collapse
Affiliation(s)
- André C. Pereira
- Centre for Ecology, Evolution and Environmental Changes (cE3c) & CHANGE – Global Change and Sustainability Institute, Faculdade de Ciências, Universidade de Lisboa, Lisboa, Portugal
- Biosystems & Integrative Sciences Institute (BioISI), Faculdade de Ciências, Universidade de Lisboa, Lisboa, Portugal
| | - Ana C. Reis
- Centre for Ecology, Evolution and Environmental Changes (cE3c) & CHANGE – Global Change and Sustainability Institute, Faculdade de Ciências, Universidade de Lisboa, Lisboa, Portugal
- Biosystems & Integrative Sciences Institute (BioISI), Faculdade de Ciências, Universidade de Lisboa, Lisboa, Portugal
| | - Mónica V. Cunha
- Centre for Ecology, Evolution and Environmental Changes (cE3c) & CHANGE – Global Change and Sustainability Institute, Faculdade de Ciências, Universidade de Lisboa, Lisboa, Portugal
- Biosystems & Integrative Sciences Institute (BioISI), Faculdade de Ciências, Universidade de Lisboa, Lisboa, Portugal
| |
Collapse
|
15
|
Yenew B, Ghodousi A, Diriba G, Tesfaye E, Cabibbe AM, Amare M, Moga S, Alemu A, Dagne B, Sinshaw W, Mollalign H, Meaza A, Tadesse M, Gamtesa DF, Abebaw Y, Seid G, Zerihun B, Getu M, Chiacchiaretta M, Gaudin C, Marceau M, Didelot X, Tolera G, Abdella S, Kebede A, Getahun M, Mehammed Z, Supply P, Cirillo DM. A smooth tubercle bacillus from Ethiopia phylogenetically close to the Mycobacterium tuberculosis complex. Nat Commun 2023; 14:7519. [PMID: 37980337 PMCID: PMC10657438 DOI: 10.1038/s41467-023-42755-9] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2022] [Accepted: 10/18/2023] [Indexed: 11/20/2023] Open
Abstract
The Mycobacterium tuberculosis complex (MTBC) includes several human- and animal-adapted pathogens. It is thought to have originated in East Africa from a recombinogenic Mycobacterium canettii-like ancestral pool. Here, we describe the discovery of a clinical tuberculosis strain isolated in Ethiopia that shares archetypal phenotypic and genomic features of M. canettii strains, but represents a phylogenetic branch much closer to the MTBC clade than to the M. canettii strains. Analysis of genomic traces of horizontal gene transfer in this isolate and previously identified M. canettii strains indicates a persistent albeit decreased recombinogenic lifestyle near the emergence of the MTBC. Our findings support that the MTBC emergence from its putative free-living M. canettii-like progenitor is evolutionarily very recent, and suggest the existence of a continuum of further extant derivatives from ancestral stages, close to the root of the MTBC, along the Great Rift Valley.
Collapse
Affiliation(s)
- Bazezew Yenew
- Ethiopian Public Health Institute, Addis Ababa, Ethiopia
| | - Arash Ghodousi
- Vita-Salute San Raffaele University, Milan, Italy.
- IRCCS San Raffaele Scientific Institute, Milan, Italy.
| | - Getu Diriba
- Ethiopian Public Health Institute, Addis Ababa, Ethiopia
| | - Ephrem Tesfaye
- Ethiopian Public Health Institute, Addis Ababa, Ethiopia
| | | | - Misikir Amare
- Ethiopian Public Health Institute, Addis Ababa, Ethiopia
| | - Shewki Moga
- Ethiopian Public Health Institute, Addis Ababa, Ethiopia
| | - Ayinalem Alemu
- Ethiopian Public Health Institute, Addis Ababa, Ethiopia
| | - Binyam Dagne
- Ethiopian Public Health Institute, Addis Ababa, Ethiopia
| | | | | | - Abyot Meaza
- Ethiopian Public Health Institute, Addis Ababa, Ethiopia
| | | | | | | | - Getachew Seid
- Ethiopian Public Health Institute, Addis Ababa, Ethiopia
| | | | - Melak Getu
- Ethiopian Public Health Institute, Addis Ababa, Ethiopia
| | | | - Cyril Gaudin
- Univ. Lille, CNRS, Inserm, CHU Lille, Institut Pasteur de Lille, U1019 - UMR 9017 - CIIL - Center for Infection and Immunity of Lille, F-59000, Lille, France
| | - Michael Marceau
- Univ. Lille, CNRS, Inserm, CHU Lille, Institut Pasteur de Lille, U1019 - UMR 9017 - CIIL - Center for Infection and Immunity of Lille, F-59000, Lille, France
| | - Xavier Didelot
- School of Life Sciences and Department of Statistics, University of Warwick, CV4 7AL, Coventry, UK
| | | | - Saro Abdella
- Ethiopian Public Health Institute, Addis Ababa, Ethiopia
| | - Abebaw Kebede
- Ethiopian Public Health Institute, Addis Ababa, Ethiopia
| | | | | | - Philip Supply
- Univ. Lille, CNRS, Inserm, CHU Lille, Institut Pasteur de Lille, U1019 - UMR 9017 - CIIL - Center for Infection and Immunity of Lille, F-59000, Lille, France.
| | - Daniela Maria Cirillo
- Vita-Salute San Raffaele University, Milan, Italy.
- IRCCS San Raffaele Scientific Institute, Milan, Italy.
| |
Collapse
|
16
|
Rott P, Grinstead S, Dallot S, Foster ZSL, Daugrois JH, Fernandez E, Kaye CJ, Hendrickson L, Hu X, Adhikari B, Malapi M, Grünwald NJ, Roumagnac P, Mollov D. Genetic Diversity, Evolution, and Diagnosis of Sugarcane Yellow Leaf Virus from 19 Sugarcane-Producing Locations Worldwide. PLANT DISEASE 2023; 107:3437-3447. [PMID: 37079008 DOI: 10.1094/pdis-10-22-2405-re] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
Sugarcane yellow leaf virus (SCYLV), the causal agent of yellow leaf, has been reported in an increasing number of sugarcane-growing locations since its first report in the 1990s in Brazil, Florida, and Hawaii. In this study, the genetic diversity of SCYLV was investigated using the genome coding sequence (5,561 to 5,612 nt) of 109 virus isolates from 19 geographical locations, including 65 new isolates from 16 geographical regions worldwide. These isolates were distributed in three major phylogenetic lineages (BRA, CUB, and REU), except for one isolate from Guatemala. Twenty-two recombination events were identified among the 109 isolates of SCYLV, thus confirming that recombination was a significant driving force in the genetic diversity and evolution of this virus. No temporal signal was found in the genomic sequence dataset, most likely because of the short temporal window of the 109 SCYLV isolates (1998 to 2020). Among 27 primers reported in the literature for the detection of the virus by RT-PCR, none matched 100% with all 109 SCYLV sequences, suggesting that the use of some primer pairs may not result in the detection of all virus isolates. Primers YLS111/YLS462, which were the first primer pair used by numerous research organizations to detect the virus by RT-PCR, failed to detect isolates belonging to the CUB lineage. In contrast, primer pair ScYLVf1/ScYLVr1 efficiently detected isolates of all three lineages. Continuous pursuit of knowledge of SCYLV genetic variability is therefore critical for effective diagnosis of yellow leaf, especially in virus-infected and mainly asymptomatic sugarcane plants.
Collapse
Affiliation(s)
- Philippe Rott
- CIRAD, UMR PHIM, 34398 Montpellier, France
- PHIM Plant Health Institute, Univ Montpellier, CIRAD, INRAE, Institut Agro, IRD, Montpellier, France
| | - Sam Grinstead
- National Germplasm Resources Laboratory, USDA-ARS, Beltsville, MD 20705, U.S.A
| | - Sylvie Dallot
- CIRAD, UMR PHIM, 34398 Montpellier, France
- PHIM Plant Health Institute, Univ Montpellier, CIRAD, INRAE, Institut Agro, IRD, Montpellier, France
| | - Zachary S L Foster
- Horticultural Crops Disease and Pest Management Research Unit, USDA-ARS, Corvallis, OR 97330, U.S.A
| | - Jean H Daugrois
- CIRAD, UMR PHIM, 34398 Montpellier, France
- PHIM Plant Health Institute, Univ Montpellier, CIRAD, INRAE, Institut Agro, IRD, Montpellier, France
| | - Emmanuel Fernandez
- CIRAD, UMR PHIM, 34398 Montpellier, France
- PHIM Plant Health Institute, Univ Montpellier, CIRAD, INRAE, Institut Agro, IRD, Montpellier, France
| | | | | | - Xiaojun Hu
- Plant Germplasm Quarantine Program, USDA-APHIS, Beltsville, MD 20705, U.S.A
| | - Bishwo Adhikari
- Plant Germplasm Quarantine Program, USDA-APHIS, Beltsville, MD 20705, U.S.A
| | - Martha Malapi
- Plant Germplasm Quarantine Program, USDA-APHIS, Beltsville, MD 20705, U.S.A
| | - Niklaus J Grünwald
- Horticultural Crops Disease and Pest Management Research Unit, USDA-ARS, Corvallis, OR 97330, U.S.A
| | - Philippe Roumagnac
- CIRAD, UMR PHIM, 34398 Montpellier, France
- PHIM Plant Health Institute, Univ Montpellier, CIRAD, INRAE, Institut Agro, IRD, Montpellier, France
| | - Dimitre Mollov
- Horticultural Crops Disease and Pest Management Research Unit, USDA-ARS, Corvallis, OR 97330, U.S.A
| |
Collapse
|
17
|
Silcocks M, Dunstan SJ. Parallel signatures of Mycobacterium tuberculosis and human Y-chromosome phylogeography support the Two Layer model of East Asian population history. Commun Biol 2023; 6:1037. [PMID: 37833496 PMCID: PMC10575886 DOI: 10.1038/s42003-023-05388-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2023] [Accepted: 09/25/2023] [Indexed: 10/15/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Matthew Silcocks
- Department of Infectious Diseases, University of Melbourne at the Peter Doherty Institute for Infection and Immunity, Parkville, VIC, Australia.
| | - Sarah J Dunstan
- Department of Infectious Diseases, University of Melbourne at the Peter Doherty Institute for Infection and Immunity, Parkville, VIC, Australia
| |
Collapse
|
18
|
Lempens P, Van Deun A, Aung KJM, Hossain MA, Behruznia M, Decroo T, Rigouts L, de Jong BC, Meehan CJ. Borderline rpoB mutations transmit at the same rate as common rpoB mutations in a tuberculosis cohort in Bangladesh. Microb Genom 2023; 9:001109. [PMID: 37750750 PMCID: PMC10569737 DOI: 10.1099/mgen.0.001109] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2023] [Accepted: 09/18/2023] [Indexed: 09/27/2023] Open
Abstract
The spread of multidrug-resistant tuberculosis (MDR-TB) is a growing problem in many countries worldwide. Resistance to one of the primary first-line drugs, rifampicin, is caused by mutations in the Mycobacterium tuberculosis rpoB gene. So-called borderline rpoB mutations confer low-level resistance, in contrast to more common rpoB mutations which confer high-level resistance. While some borderline mutations show lower fitness in vitro than common mutations, their in vivo fitness is currently unknown. We used a dataset of 394 whole genome sequenced MDR-TB isolates from Bangladesh, representing around 44 % of notified MDR-TB cases over 6 years, to look at differences in transmission clustering between isolates with borderline rpoB mutations and those with common rpoB mutations. We found a relatively low percentage of transmission clustering in the dataset (34.8 %) but no difference in clustering between different types of rpoB mutations. Compensatory mutations in rpoA, rpoB, and rpoC were associated with higher levels of transmission clustering as were lineages two, three, and four relative to lineage one. Young people as well as patients with high sputum smear positive TB were more likely to be in a transmission cluster. Our findings show that although borderline rpoB mutations have lower in vitro growth potential this does not translate into lower transmission potential or in vivo fitness. Proper detection of these mutations is crucial to ensure they do not go unnoticed and spread MDR-TB within communities.
Collapse
Affiliation(s)
- Pauline Lempens
- Unit of Mycobacteriology, Department of Biomedical Sciences, Institute of Tropical Medicine, Antwerp, Belgium
- Department of Biomedical Sciences, University of Antwerp, Antwerp, Belgium
| | | | | | | | | | - Tom Decroo
- Unit of HIV and TB, Department of Clinical Sciences, Institute of Tropical Medicine, Antwerp, Belgium
| | - Leen Rigouts
- Unit of Mycobacteriology, Department of Biomedical Sciences, Institute of Tropical Medicine, Antwerp, Belgium
- Department of Biomedical Sciences, University of Antwerp, Antwerp, Belgium
| | - Bouke C. de Jong
- Unit of Mycobacteriology, Department of Biomedical Sciences, Institute of Tropical Medicine, Antwerp, Belgium
| | - Conor J. Meehan
- Unit of Mycobacteriology, Department of Biomedical Sciences, Institute of Tropical Medicine, Antwerp, Belgium
- Department of Biosciences, Nottingham Trent University, Nottingham, UK
| |
Collapse
|
19
|
Vargas R, Luna MJ, Freschi L, Marin M, Froom R, Murphy KC, Campbell EA, Ioerger TR, Sassetti CM, Farhat MR. Phase variation as a major mechanism of adaptation in Mycobacterium tuberculosis complex. Proc Natl Acad Sci U S A 2023; 120:e2301394120. [PMID: 37399390 PMCID: PMC10334774 DOI: 10.1073/pnas.2301394120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2023] [Accepted: 05/03/2023] [Indexed: 07/05/2023] Open
Abstract
Phase variation induced by insertions and deletions (INDELs) in genomic homopolymeric tracts (HT) can silence and regulate genes in pathogenic bacteria, but this process is not characterized in MTBC (Mycobacterium tuberculosis complex) adaptation. We leverage 31,428 diverse clinical isolates to identify genomic regions including phase-variants under positive selection. Of 87,651 INDEL events that emerge repeatedly across the phylogeny, 12.4% are phase-variants within HTs (0.02% of the genome by length). We estimated the in-vitro frameshift rate in a neutral HT at 100× the neutral substitution rate at [Formula: see text] frameshifts/HT/year. Using neutral evolution simulations, we identified 4,098 substitutions and 45 phase-variants to be putatively adaptive to MTBC (P < 0.002). We experimentally confirm that a putatively adaptive phase-variant alters the expression of espA, a critical mediator of ESX-1-dependent virulence. Our evidence supports the hypothesis that phase variation in the ESX-1 system of MTBC can act as a toggle between antigenicity and survival in the host.
Collapse
Affiliation(s)
- Roger Vargas
- Center for Computational Biomedicine, Harvard Medical School, Boston, MA02115
- Department of Biomedical Informatics, Harvard Medical School, Boston, MA02115
| | - Michael J. Luna
- Department of Microbiology and Physiological Systems, University of Massachusetts Chan Medical School, Worcester, MA01655
| | - Luca Freschi
- Department of Biomedical Informatics, Harvard Medical School, Boston, MA02115
| | - Maximillian Marin
- Department of Biomedical Informatics, Harvard Medical School, Boston, MA02115
| | - Ruby Froom
- Laboratory of Molecular Biophysics, The Rockefeller University, New York, NY10065
- Laboratory of Host-Pathogen Biology, The Rockefeller University, New York, NY10065
| | - Kenan C. Murphy
- Department of Microbiology and Physiological Systems, University of Massachusetts Chan Medical School, Worcester, MA01655
| | | | - Thomas R. Ioerger
- Department of Computer Science and Engineering, Texas A&M University, College Station, TX77843
| | - Christopher M. Sassetti
- Department of Microbiology and Physiological Systems, University of Massachusetts Chan Medical School, Worcester, MA01655
| | - Maha Reda Farhat
- Department of Biomedical Informatics, Harvard Medical School, Boston, MA02115
- Pulmonary and Critical Care Medicine, Massachusetts General Hospital, Boston, MA02114
| |
Collapse
|
20
|
Goig GA, Menardo F, Salaam-Dreyer Z, Dippenaar A, Streicher EM, Daniels J, Reuter A, Borrell S, Reinhard M, Doetsch A, Beisel C, Warren RM, Cox H, Gagneux S. Effect of compensatory evolution in the emergence and transmission of rifampicin-resistant Mycobacterium tuberculosis in Cape Town, South Africa: a genomic epidemiology study. THE LANCET. MICROBE 2023; 4:e506-e515. [PMID: 37295446 PMCID: PMC10319636 DOI: 10.1016/s2666-5247(23)00110-6] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/28/2022] [Revised: 01/26/2023] [Accepted: 03/01/2023] [Indexed: 06/12/2023]
Abstract
BACKGROUND Experimental data show that drug-resistance-conferring mutations are often associated with a decrease in the replicative fitness of bacteria in vitro, and that this fitness cost can be mitigated by compensatory mutations; however, the role of compensatory evolution in clinical settings is less clear. We assessed whether compensatory evolution was associated with increased transmission of rifampicin-resistant tuberculosis in Khayelitsha, Cape Town, South Africa. METHODS We did a genomic epidemiological study by analysing available M tuberculosis isolates and their associated clinical data from individuals routinely diagnosed with rifampicin-resistant tuberculosis in primary care and hospitals in Khayelitsha, Cape Town, South Africa. Isolates were collected as part of a previous study. All individuals diagnosed with rifampicin-resistant tuberculosis and with linked biobanked specimens were included in this study. We applied whole-genome sequencing, Bayesian reconstruction of transmission trees, and phylogenetic multivariable regression analysis to identify individual and bacterial factors associated with the transmission of rifampicin-resistant M tuberculosis strains. FINDINGS Between Jan 1, 2008, and Dec 31, 2017, 2161 individuals were diagnosed with multidrug-resistant or rifampicin-resistant tuberculosis in Khayelitsha, Cape Town, South Africa. Whole-genome sequences were available for 1168 (54%) unique individual M tuberculosis isolates. Compensatory evolution was associated with smear-positive pulmonary disease (adjusted odds ratio 1·49, 95% CI 1·08-2·06) and a higher number of drug-resistance-conferring mutations (incidence rate ratio 1·38, 95% CI 1·28-1·48). Compensatory evolution was also associated with increased transmission of rifampicin-resistant disease between individuals (adjusted odds ratio 1·55; 95% CI 1·13-2·12), independent of other patient and bacterial factors. INTERPRETATION Our findings suggest that compensatory evolution enhances the in vivo fitness of drug-resistant M tuberculosis genotypes, both within and between patients, and that the in vitro replicative fitness of rifampicin-resistant M tuberculosis measured in the laboratory correlates with the bacterial fitness measured in clinical settings. These results emphasise the importance of enhancing surveillance and monitoring efforts to prevent the emergence of highly transmissible clones capable of rapidly accumulating new drug resistance mutations. This concern becomes especially crucial at present, because treatment regimens incorporating novel drugs are being implemented. FUNDING Funding for this study was provided by a Swiss and South Africa joint research award (grant numbers 310030_188888, CRSII5_177163, and IZLSZ3_170834), the European Research Council (grant number 883582), and a Wellcome Trust fellowship (to HC; reference number 099818/Z/12/Z). ZS-D was funded through a PhD scholarship from the South African National Research Foundation and RMW was funded through the South African Medical Research Council.
Collapse
Affiliation(s)
- Galo A Goig
- Medical Parasitology and Infection Biology, Swiss Tropical and Public Health Institute, Allschwil, Switzerland; University of Basel, Basel, Switzerland.
| | - Fabrizio Menardo
- Department of Plant and Microbial Biology, University of Zürich, Zürich, Switzerland
| | - Zubeida Salaam-Dreyer
- Division of Medical Microbiology, Department of Pathology, University of Cape Town, Cape Town, South Africa
| | - Anzaan Dippenaar
- Tuberculosis Omics Research Consortium, Family Medicine and Population Health, Institute of Global Health, Faculty of Medicine and Health Sciences, University of Antwerp, Antwerp, Belgium
| | - Elizabeth M Streicher
- Department of Science and Innovation - National Research Foundation Centre of Excellence for Biomedical Tuberculosis Research, Stellenbosch University, Stellenbosch, South Africa; South African Medical Research Council Centre for Tuberculosis Research, Division of Molecular Biology and Human Genetics, Faculty of Medicine and Health Sciences, Stellenbosch University, Stellenbosch, South Africa
| | - Johnny Daniels
- Médecins Sans Frontières, Khayelitsha, Cape Town, South Africa
| | - Anja Reuter
- Médecins Sans Frontières, Khayelitsha, Cape Town, South Africa
| | - Sonia Borrell
- Medical Parasitology and Infection Biology, Swiss Tropical and Public Health Institute, Allschwil, Switzerland; University of Basel, Basel, Switzerland
| | - Miriam Reinhard
- Medical Parasitology and Infection Biology, Swiss Tropical and Public Health Institute, Allschwil, Switzerland; University of Basel, Basel, Switzerland
| | - Anna Doetsch
- Medical Parasitology and Infection Biology, Swiss Tropical and Public Health Institute, Allschwil, Switzerland; University of Basel, Basel, Switzerland
| | - Christian Beisel
- Department of Biosystems Science and Engineering, Eidgenössische Technische Hochschule Zürich, Zürich, Swizterland
| | - Robin M Warren
- Department of Science and Innovation - National Research Foundation Centre of Excellence for Biomedical Tuberculosis Research, Stellenbosch University, Stellenbosch, South Africa; South African Medical Research Council Centre for Tuberculosis Research, Division of Molecular Biology and Human Genetics, Faculty of Medicine and Health Sciences, Stellenbosch University, Stellenbosch, South Africa
| | - Helen Cox
- Division of Medical Microbiology, Department of Pathology, University of Cape Town, Cape Town, South Africa; Institute of Infectious Disease and Molecular Medicine and Wellcome Centre for Infectious Disease Research, University of Cape Town, Cape Town, South Africa
| | - Sebastien Gagneux
- Medical Parasitology and Infection Biology, Swiss Tropical and Public Health Institute, Allschwil, Switzerland; University of Basel, Basel, Switzerland
| |
Collapse
|
21
|
Abstract
Mycobacteria are responsible for several human and animal diseases. NOD2 is a pattern recognition receptor that has an important role in mycobacterial recognition. However, the mechanisms by which mutations in NOD2 alter the course of mycobacterial infection remain unclear. Herein, we aimed to review the totality of studies directly addressing the relationship between NOD2 and mycobacteria as a foundation for moving the field forward. NOD2 was linked to mycobacterial infection at 3 levels: (1) genetic, through association with mycobacterial diseases of humans; (2) chemical, through the distinct NOD2 ligand in the mycobacterial cell wall; and (3) immunologic, through heightened NOD2 signaling caused by the unique modification of the NOD2 ligand. The immune response to mycobacteria is shaped by NOD2 signaling, responsible for NF-κB and MAPK activation, and the production of various immune effectors like cytokines and nitric oxide, with some evidence linking this to bacteriologic control. Absence of NOD2 during mycobacterial infection of mice can be detrimental, but the mechanism remains unknown. Conversely, the success of immunization with mycobacteria has been linked to NOD2 signaling and NOD2 has been targeted as an avenue of immunotherapy for diseases even beyond mycobacteria. The mycobacteria-NOD2 interaction remains an important area of study, which may shed light on immune mechanisms in disease.
Collapse
Affiliation(s)
- Jean-Yves Dubé
- Department of Microbiology and Immunology, McGill University, Montréal, Canada
| | - Marcel A. Behr
- Department of Medicine, McGill University Health Centre, Montréal, Canada
| |
Collapse
|
22
|
Sanders JG, Sprockett DD, Li Y, Mjungu D, Lonsdorf EV, Ndjango JBN, Georgiev AV, Hart JA, Sanz CM, Morgan DB, Peeters M, Hahn BH, Moeller AH. Widespread extinctions of co-diversified primate gut bacterial symbionts from humans. Nat Microbiol 2023; 8:1039-1050. [PMID: 37169918 PMCID: PMC10860671 DOI: 10.1038/s41564-023-01388-w] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2022] [Accepted: 04/19/2023] [Indexed: 05/13/2023]
Abstract
Humans and other primates harbour complex gut bacterial communities that influence health and disease, but the evolutionary histories of these symbioses remain unclear. This is partly due to limited information about the microbiota of ancestral primates. Here, using phylogenetic analyses of metagenome-assembled genomes (MAGs), we show that hundreds of gut bacterial clades diversified in parallel (that is, co-diversified) with primate species over millions of years, but that humans have experienced widespread losses of these ancestral symbionts. Analyses of 9,460 human and non-human primate MAGs, including newly generated MAGs from chimpanzees and bonobos, revealed significant co-diversification within ten gut bacterial phyla, including Firmicutes, Actinobacteriota and Bacteroidota. Strikingly, ~44% of the co-diversifying clades detected in African apes were absent from available metagenomic data from humans and ~54% were absent from industrialized human populations. In contrast, only ~3% of non-co-diversifying clades detected in African apes were absent from humans. Co-diversifying clades present in both humans and chimpanzees displayed consistent genomic signatures of natural selection between the two host species but differed in functional content from co-diversifying clades lost from humans, consistent with selection against certain functions. This study discovers host-species-specific bacterial symbionts that predate hominid diversification, many of which have undergone accelerated extinctions from human populations.
Collapse
Affiliation(s)
- Jon G Sanders
- Department of Ecology and Evolutionary Biology, Cornell University, Ithaca, NY, USA
| | - Daniel D Sprockett
- Department of Ecology and Evolutionary Biology, Cornell University, Ithaca, NY, USA
| | - Yingying Li
- Departments of Medicine and Microbiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Deus Mjungu
- Gombe Stream Research Center, Kigoma, Tanzania
| | - Elizabeth V Lonsdorf
- Department of Psychology and Biological Foundations of Behavior Program, Franklin and Marshall College, Lancaster, PA, USA
- Department of Anthropology, Emory University, Atlanta, GA, USA
| | - Jean-Bosco N Ndjango
- Department of Ecology and Management of Plant and Animal Resources, Faculty of Sciences, University of Kisangani, Kisangani, Democratic Republic of the Congo
| | - Alexander V Georgiev
- Department of Human Evolutionary Biology, Harvard University, Cambridge, MA, USA
- School of Natural Sciences, Bangor University, Bangor, UK
| | - John A Hart
- Lukuru Wildlife Research Foundation, Tshuapa-Lomami-Lualaba Project, Kinshasa, Democratic Republic of the Congo
| | - Crickette M Sanz
- Department of Anthropology, Washington University in St Louis, Saint Louis, MO, USA
- Wildlife Conservation Society, Congo Program, Brazzaville, Republic of Congo
| | - David B Morgan
- Lester E. Fisher Center for the Study and Conservation of Apes, Lincoln Park Zoo, Chicago, IL, USA
| | - Martine Peeters
- Recherche Translationnelle Appliquée au VIH et aux Maladies Infectieuses, Institut de Recherche pour le Développement, University of Montpellier, INSERM, Montpellier, France
| | - Beatrice H Hahn
- Departments of Medicine and Microbiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Andrew H Moeller
- Department of Ecology and Evolutionary Biology, Cornell University, Ithaca, NY, USA.
| |
Collapse
|
23
|
Cramer N, Klockgether J, Tümmler B. Microevolution of Pseudomonas aeruginosa in the airways of people with cystic fibrosis. Curr Opin Immunol 2023; 83:102328. [PMID: 37116385 DOI: 10.1016/j.coi.2023.102328] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2022] [Revised: 03/26/2023] [Accepted: 03/27/2023] [Indexed: 04/30/2023]
Abstract
The chronic infections of cystic fibrosis (CF) airways with Pseudomonas aeruginosa are a paradigm of how environmental bacteria can conquer, adapt, and persist in an atypical habitat and successfully evade defense mechanisms and chemotherapy in a susceptible host. The within-host evolution of intraclonal diversity has been examined by whole-genome sequencing, phenotyping, and competitive fitness experiments of serial P. aeruginosa isolates collected from CF airways since onset of colonization for a period of up to 40 years. The spectrum of de novo mutations and the adaptation of phenotype and fitness of the bacterial progeny were more influenced by the living conditions in the CF lung than by the clone type of their ancestor and its genetic repertoire.
Collapse
Affiliation(s)
- Nina Cramer
- Department of Pediatric Pneumology, Allergology and Neonatology, Hannover Medical School, D-30625 Hannover, Germany
| | - Jens Klockgether
- Department of Pediatric Pneumology, Allergology and Neonatology, Hannover Medical School, D-30625 Hannover, Germany
| | - Burkhard Tümmler
- Department of Pediatric Pneumology, Allergology and Neonatology, Hannover Medical School, D-30625 Hannover, Germany; Biomedical Research in Endstage and Obstructive Lung Disease Hannover (BREATH), Member of the German Center for Lung Research (DZL), Hannover, Germany.
| |
Collapse
|
24
|
Żukowska L, Zygała-Pytlos D, Struś K, Zabost A, Kozińska M, Augustynowicz-Kopeć E, Dziadek J, Minias A. An overview of tuberculosis outbreaks reported in the years 2011-2020. BMC Infect Dis 2023; 23:253. [PMID: 37081448 PMCID: PMC10116450 DOI: 10.1186/s12879-023-08197-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2022] [Accepted: 03/24/2023] [Indexed: 04/22/2023] Open
Abstract
BACKGROUND In many countries tuberculosis (TB) remains a highly prevalent disease and a major contributor to infectious disease mortality. The fight against TB requires surveillance of the population of strains circulating worldwide and the analysis of the prevalence of certain strains in populations. Nowadays, whole genome sequencing (WGS) allows for accurate tracking of TB transmission. Currently, there is a lack of a comprehensive summary of the characteristics of TB outbreaks. METHODS We systematically analyzed studies reporting TB outbreaks worldwide, monitored through WGS of Mycobacterium tuberculosis. We 1) mapped the reported outbreaks from 2011- 2020, 2) estimated the average size of the outbreaks, 3) indicated genetic lineages causing the outbreaks, and 4) determined drug-resistance patterns of M. tuberculosis strains involved in the outbreaks. RESULTS Most data originated from Europe, Asia, and North America. We found that TB outbreaks were reported throughout the globe, on all continents, and in countries with both high and low incidences. The detected outbreaks contained a median of five M. tuberculosis isolates. Most strains causing the outbreaks belonged to lineage four, more rarely to lineage two. Reported outbreak isolates were often drug resistant. CONCLUSIONS We conclude that more WGS surveillance of M. tuberculosis outbreaks is needed. Globally standardized procedures might improve the control of M. tuberculosis infections.
Collapse
Affiliation(s)
- Lidia Żukowska
- Institute of Medical Biology, Polish Academy of Sciences, Lodz, Poland
- The Bio-Med-Chem Doctoral School of the University of Lodz and Lodz Institutes of the Polish Academy of Sciences, Lodz, Poland
| | - Daria Zygała-Pytlos
- Institute of Medical Biology, Polish Academy of Sciences, Lodz, Poland
- The Bio-Med-Chem Doctoral School of the University of Lodz and Lodz Institutes of the Polish Academy of Sciences, Lodz, Poland
| | - Katarzyna Struś
- Institute of Medical Biology, Polish Academy of Sciences, Lodz, Poland
- Institute of Biology and Biotechnology, College of Natural Sciences, University of Rzeszów, Rzeszów, Poland
| | - Anna Zabost
- Department of Microbiology, National Tuberculosis and Lung Diseases Research Institute, Warsaw, Poland
| | - Monika Kozińska
- Department of Microbiology, National Tuberculosis and Lung Diseases Research Institute, Warsaw, Poland
| | - Ewa Augustynowicz-Kopeć
- Department of Microbiology, National Tuberculosis and Lung Diseases Research Institute, Warsaw, Poland
| | - Jarosław Dziadek
- Institute of Medical Biology, Polish Academy of Sciences, Lodz, Poland
| | - Alina Minias
- Institute of Medical Biology, Polish Academy of Sciences, Lodz, Poland.
| |
Collapse
|
25
|
Jiang Q, Liu HC, Liu QY, Phelan JE, Tao FX, Zhao XQ, Wang J, Glynn JR, Takiff HE, Clark TG, Wan KL, Gao Q. The Evolution and Transmission Dynamics of Multidrug-Resistant Tuberculosis in an Isolated High-Plateau Population of Tibet, China. Microbiol Spectr 2023; 11:e0399122. [PMID: 36912683 PMCID: PMC10101056 DOI: 10.1128/spectrum.03991-22] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2022] [Accepted: 02/15/2023] [Indexed: 03/14/2023] Open
Abstract
On the Tibetan Plateau, most tuberculosis is caused by indigenous Mycobacterium tuberculosis strains with a monophyletic structure and high-level drug resistance. This study investigated the emergence, evolution, and transmission dynamics of multidrug-resistant tuberculosis (MDR-TB) in Tibet. The whole-genome sequences of 576 clinical strains from Tibet were analyzed with the TB-profiler tool to identify drug-resistance mutations. The evolution of the drug resistance was then inferred based on maximum-likelihood phylogeny and dated trees that traced the serial acquisition of mutations conferring resistance to different drugs. Among the 576 clinical M. tuberculosis strains, 346 (60.1%) carried at least 1 resistance-conferring mutation and 231 (40.1%) were MDR-TB. Using a pairwise distance of 50 single nucleotide polymorphisms (SNPs), most strains (89.9%, 518/576) were phylogenetically separated into 50 long-term transmission clusters. Eleven large drug-resistant clusters contained 76.1% (176/231) of the local multidrug-resistant strains. A total of 85.2% of the isoniazid-resistant strains were highly transmitted with an average of 6.6 cases per cluster, of which most shared the mutation KatG Ser315Thr. A lower proportion (71.6%) of multidrug-resistant strains were transmitted, with an average cluster size of 2.9 cases. The isoniazid-resistant clusters appear to have undergone substantial bacterial population growth in the 1970s to 1990s and then subsequently accumulated multiple rifampicin-resistance mutations and caused the current local MDR-TB burden. These findings highlight the importance of detecting and curing isoniazid-resistant strains to prevent the emergence of endemic MDR-TB. IMPORTANCE Emerging isoniazid resistance in the 1970s allowed M. tuberculosis strains to spread and form into large multidrug-resistant tuberculosis clusters in the isolated plateau of Tibet, China. The epidemic was driven by the high risk of transmission as well as the potential of acquiring further drug resistance from isoniazid-resistant strains. Eleven large drug-resistant clusters consisted of the majority of local multidrug-resistant cases. Among the clusters, isoniazid resistance overwhelmingly evolved before all the other resistance types. A large bacterial population growth of isoniazid-resistant clusters occurred between 1970s and 1990s, which subsequently accumulated rifampicin-resistance-conferring mutations in parallel and accounted for the local multidrug-resistant tuberculosis burden. The results of our study indicate that it may be possible to restrict MDR-TB evolution and dissemination by prioritizing screening for isoniazid (INH)-resistant TB strains before they become MDR-TB and by adopting measures that can limit their transmission.
Collapse
Affiliation(s)
- Qi Jiang
- Department of Epidemiology and Biostatistics, School of Public Health, Wuhan University, Wuhan, China
| | - Hai-Can Liu
- State Key Laboratory for Infectious Disease Prevention and Control and National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing, China
| | - Qing-Yun Liu
- Department of Immunology and Infectious Diseases, Harvard T. H. Chan School of Public Health, Boston, Massachusetts, USA
| | - Jody E. Phelan
- Department of Infection Biology, London School of Hygiene & Tropical Medicine, London, United Kingdom
| | - Feng-Xi Tao
- Department of Epidemiology and Biostatistics, School of Public Health, Wuhan University, Wuhan, China
| | - Xiu-Qin Zhao
- State Key Laboratory for Infectious Disease Prevention and Control and National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing, China
| | - Jian Wang
- Tibet Center for Disease Control and Prevention, Lhasa, Tibet Autonomous Region, China
| | - Judith R. Glynn
- Department of Infectious Disease Epidemiology, London School of Hygiene & Tropical Medicine, London, United Kingdom
| | - Howard E. Takiff
- Laboratorio de Genética Molecular, CMBC, IVIC, Caracas, Venezuela
| | - Taane G. Clark
- Department of Infection Biology, London School of Hygiene & Tropical Medicine, London, United Kingdom
- Department of Infectious Disease Epidemiology, London School of Hygiene & Tropical Medicine, London, United Kingdom
| | - Kang-Lin Wan
- State Key Laboratory for Infectious Disease Prevention and Control and National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing, China
| | - Qian Gao
- National Clinical Research Center for Infectious Diseases, The Third People’s Hospital of Shenzhen, Shenzhen, Guangdong, China
| |
Collapse
|
26
|
López MG, Campos-Herrero MI, Torres-Puente M, Cañas F, Comín J, Copado R, Wintringer P, Iqbal Z, Lagarejos E, Moreno-Molina M, Pérez-Lago L, Pino B, Sante L, García de Viedma D, Samper S, Comas I. Deciphering the Tangible Spatio-Temporal Spread of a 25-Year Tuberculosis Outbreak Boosted by Social Determinants. Microbiol Spectr 2023; 11:e0282622. [PMID: 36786614 PMCID: PMC10100973 DOI: 10.1128/spectrum.02826-22] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2022] [Accepted: 01/18/2023] [Indexed: 02/15/2023] Open
Abstract
Outbreak strains of Mycobacterium tuberculosis are promising candidates as targets in the search for intrinsic determinants of transmissibility, as they are responsible for many cases with sustained transmission; however, the use of low-resolution typing methods and restricted geographical investigations represent flaws in assessing the success of long-lived outbreak strains. We can now address the nature of outbreak strains by combining large genomic data sets and phylodynamic approaches. We retrospectively sequenced the whole genome of representative samples assigned to an outbreak circulating in the Canary Islands (the GC strain) since 1993, which accounts for ~20% of local tuberculosis cases. We selected a panel of specific single nucleotide polymorphism (SNP) markers for an in-silico search for additional outbreak-related sequences within publicly available tuberculosis genomic data. Using this information, we inferred the origin, spread, and epidemiological parameters of the GC strain. Our approach allowed us to accurately trace the historical and more recent dispersion of the GC strain. We provide evidence of a highly successful nature within the Canarian archipelago but limited expansion abroad. Estimation of epidemiological parameters from genomic data disagree with a distinctive biology of the GC strain. With the increasing availability of genomic data allowing for the accurate inference of strain spread and critical epidemiological parameters, we can now revisit the link between Mycobacterium tuberculosis genotypes and transmission, as is routinely carried out for SARS-CoV-2 variants of concern. We demonstrate that social determinants rather than intrinsically higher bacterial transmissibility better explain the success of the GC strain. Importantly, our approach can be used to trace and characterize strains of interest worldwide. IMPORTANCE Infectious disease outbreaks represent a significant problem for public health. Tracing outbreak expansion and understanding the main factors behind emergence and persistence remain critical to effective disease control. Our study allows researchers and public health authorities to use Whole-Genome Sequencing-based methods to trace outbreaks, and shows how available epidemiological information helps to evaluate the factors underpinning outbreak persistence. Taking advantage of all the freely available information placed in public repositories, researchers can accurately establish the expansion of an outbreak beyond original boundaries, and determine the potential risk of a strain to inform health authorities which, in turn, can define target strategies to mitigate expansion and persistence. Finally, we show the need to evaluate strain transmissibility in different geographic contexts to unequivocally associate spread to local or pathogenic factors, an important lesson taken from genomic surveillance of SARS-CoV-2.
Collapse
Affiliation(s)
- Mariana G. López
- Tuberculosis Genomics Unit, Instituto de Biomedicina de Valencia (IBV), CSIC, Valencia, Spain
| | - Ma Isolina Campos-Herrero
- Servicio de Microbiología, Hospital Universitario de Gran Canaria Dr. Negrín, Las Palmas de Gran Canaria, Spain
| | - Manuela Torres-Puente
- Tuberculosis Genomics Unit, Instituto de Biomedicina de Valencia (IBV), CSIC, Valencia, Spain
| | - Fernando Cañas
- Hospital Universitario Insular de Gran Canaria, Las Palmas de Gran Canaria, Spain
| | - Jessica Comín
- Instituto Aragonés de Ciencias de la Salud, Fundación IIS Aragón, Zaragoza, Spain
| | - Rodolfo Copado
- Hospital José Molina Orosa, Las Palmas de Gran Canaria, Spain
| | - Penelope Wintringer
- European Molecular Biology Laboratory – European Bioinformatics Institute, Hinxton, UK
| | - Zamin Iqbal
- European Molecular Biology Laboratory – European Bioinformatics Institute, Hinxton, UK
| | - Eduardo Lagarejos
- Servicio de Microbiología, Hospital Universitario de Gran Canaria Dr. Negrín, Las Palmas de Gran Canaria, Spain
| | - Miguel Moreno-Molina
- Tuberculosis Genomics Unit, Instituto de Biomedicina de Valencia (IBV), CSIC, Valencia, Spain
| | - Laura Pérez-Lago
- Servicio Microbiología Clínica y Enfermedades Infecciosas, Hospital General Universitario Gregorio Marañón, Instituto de Investigación Sanitaria Gregorio Marañón, Madrid, Spain
| | - Berta Pino
- Hospital Nuestra Señora de la Candelaria, Santa Cruz de Tenerife, Spain
| | - Laura Sante
- Hospital Universitario de Canarias, Santa Cruz de Tenerife, Spain
| | - Darío García de Viedma
- Servicio Microbiología Clínica y Enfermedades Infecciosas, Hospital General Universitario Gregorio Marañón, Instituto de Investigación Sanitaria Gregorio Marañón, Madrid, Spain
- CIBER Enfermedades Respiratorias, Instituto de Salud Carlos III, Madrid, Spain
| | - Sofía Samper
- Instituto Aragonés de Ciencias de la Salud, Fundación IIS Aragón, Zaragoza, Spain
- CIBER Enfermedades Respiratorias, Instituto de Salud Carlos III, Madrid, Spain
| | - Iñaki Comas
- Tuberculosis Genomics Unit, Instituto de Biomedicina de Valencia (IBV), CSIC, Valencia, Spain
- CIBER Epidemiología y Salud Pública, Instituto de Salud Carlos III, Madrid, Spain
| |
Collapse
|
27
|
Ashton PM, Cha J, Anscombe C, Thuong NTT, Thwaites GE, Walker TM. Distribution and origins of Mycobacterium tuberculosis L4 in Southeast Asia. Microb Genom 2023; 9:mgen000955. [PMID: 36729036 PMCID: PMC9997747 DOI: 10.1099/mgen.0.000955] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2022] [Accepted: 12/21/2022] [Indexed: 02/03/2023] Open
Abstract
Molecular and genomic studies have revealed that Mycobacterium tuberculosis Lineage 4 (L4, Euro-American lineage) emerged in Europe before becoming distributed around the globe by trade routes, colonial migration and other historical connections. Although L4 accounts for tens or hundreds of thousands of tuberculosis (TB) cases in multiple Southeast Asian countries, phylogeographical studies have either focused on a single country or just included Southeast Asia as part of a global analysis. Therefore, we interrogated public genomic data to investigate the historical patterns underlying the distribution of L4 in Southeast Asia and surrounding countries. We downloaded 6037 genomes associated with 29 published studies, focusing on global analyses of L4 and Asian studies of M. tuberculosis. We identified 2256 L4 genomes including 968 from Asia. We show that 81 % of L4 in Thailand, 51 % of L4 in Vietnam and 9 % of L4 in Indonesia belong to sub-lineages of L4 that are rarely seen outside East and Southeast Asia (L4.2.2, L4.4.2 and L4.5). These sub-lineages have spread between East and Southeast Asian countries, with no recent European ancestor. Although there is considerable uncertainty about the exact direction and order of intra-Asian M. tuberculosis dispersal, due to differing sampling frames between countries, our analysis suggests that China may be the intermediate location between Europe and Southeast Asia for two of the three predominantly East and Southeast Asian L4 sub-lineages (L4.2.2 and L4.5). This new perspective on L4 in Southeast Asia raises the possibility of investigating host population-specific evolution and highlights the need for more structured sampling from Southeast Asian countries to provide more certainty of the historical and current routes of dispersal.
Collapse
Affiliation(s)
- Philip M. Ashton
- Oxford University Clinical Research Unit, Ho Chi Minh City, Vietnam
- Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, UK
- Institute of Infection, Veterinary and Ecological Sciences, University of Liverpool, Liverpool, UK
| | - Jaeyoon Cha
- Department of Molecular Biology, Princeton University, Princeton, New Jersey, USA
| | - Catherine Anscombe
- Oxford University Clinical Research Unit, Ho Chi Minh City, Vietnam
- Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | - Nguyen T. T. Thuong
- Oxford University Clinical Research Unit, Ho Chi Minh City, Vietnam
- Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | - Guy E. Thwaites
- Oxford University Clinical Research Unit, Ho Chi Minh City, Vietnam
- Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | - Timothy M. Walker
- Oxford University Clinical Research Unit, Ho Chi Minh City, Vietnam
- Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| |
Collapse
|
28
|
Sanabria GE, Sequera G, Aguirre S, Méndez J, Dos Santos PCP, Gustafson NW, Godoy M, Ortiz A, Cespedes C, Martínez G, García-Basteiro AL, Andrews JR, Croda J, Walter KS. Phylogeography and transmission of Mycobacterium tuberculosis spanning prisons and surrounding communities in Paraguay. Nat Commun 2023; 14:303. [PMID: 36658111 PMCID: PMC9849832 DOI: 10.1038/s41467-023-35813-9] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2022] [Accepted: 01/04/2023] [Indexed: 01/20/2023] Open
Abstract
Recent rises in incident tuberculosis (TB) cases in Paraguay and the increasing concentration of TB within prisons highlight the urgency of targeting strategies to interrupt transmission and prevent new infections. However, whether specific cities or carceral institutions play a disproportionate role in transmission remains unknown. We conducted prospective genomic surveillance, sequencing 471 Mycobacterium tuberculosis complex genomes, from inside and outside prisons in Paraguay's two largest urban areas, Asunción and Ciudad del Este, from 2016 to 2021. We found genomic evidence of frequent recent transmission within prisons and transmission linkages spanning prisons and surrounding populations. We identified a signal of frequent M. tuberculosis spread between urban areas and marked recent population size expansion of the three largest genomic transmission clusters. Together, our findings highlight the urgency of strengthening TB control programs to reduce transmission risk within prisons in Paraguay, where incidence was 70 times that outside prisons in 2021.
Collapse
Affiliation(s)
| | - Guillermo Sequera
- Instituto de Salud Global de Barcelona (ISGLOBAL), Barcelona, Spain
- Programa Nacional de Control de la Tuberculosis, Ministerio de Salud Pública y Bienestar Social (MSPyBS), Asunción, Paraguay
| | - Sarita Aguirre
- Programa Nacional de Control de la Tuberculosis, Ministerio de Salud Pública y Bienestar Social (MSPyBS), Asunción, Paraguay
| | - Julieta Méndez
- Instituto Regional de Investigación en Salud, Caaguazú, Paraguay
| | - Paulo César Pereira Dos Santos
- Postgraduate Program in Infectious and Parasitic Diseases, Federal University of Mato Grosso do Sul, Mato Grosso do Sul, Brazil
| | - Natalie Weiler Gustafson
- Laboratorio Central de Salud Pública (LCSP), Ministerio de Salud Publica y Bienestar Social (MSPyBS), Asunción, Paraguay
| | - Margarita Godoy
- Laboratorio Central de Salud Pública (LCSP), Ministerio de Salud Publica y Bienestar Social (MSPyBS), Asunción, Paraguay
| | - Analía Ortiz
- Instituto Regional de Investigación en Salud, Caaguazú, Paraguay
| | - Cynthia Cespedes
- Programa Nacional de Control de la Tuberculosis, Ministerio de Salud Pública y Bienestar Social (MSPyBS), Asunción, Paraguay
| | - Gloria Martínez
- Instituto Regional de Investigación en Salud, Caaguazú, Paraguay
| | - Alberto L García-Basteiro
- Instituto de Salud Global de Barcelona (ISGLOBAL), Barcelona, Spain
- Centro de Investigação em Saude de Manhiça (CISM), Maputo, Mozambique
- Centro de Investigación Biomédica en Red de Enfermedades Infecciosas (CIBERINFEC), Barcelona, Spain
| | - Jason R Andrews
- Division of Infectious Diseases and Geographic Medicine, Stanford University School of Medicine, Stanford, CA, USA
| | - Julio Croda
- Federal University of Mato Grosso do Sul - UFMS, Campo Grande, MS, Brazil
- Oswaldo Cruz Foundation Mato Grosso do Sul, Mato Grosso do Sul, Brazil
- Department of Epidemiology of Microbial Diseases, Yale School of Public Health, New Haven, USA
| | - Katharine S Walter
- Division of Epidemiology, University of Utah, Salt Lake City, UT, 84105, USA.
| |
Collapse
|
29
|
Pan J, Li X, Zhang M, Lu Y, Zhu Y, Wu K, Wu Y, Wang W, Chen B, Liu Z, Wang X, Gao J. TransFlow: a Snakemake workflow for transmission analysis of Mycobacterium tuberculosis whole-genome sequencing data. Bioinformatics 2022; 39:6873737. [PMID: 36469333 PMCID: PMC9825751 DOI: 10.1093/bioinformatics/btac785] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2022] [Revised: 10/26/2022] [Accepted: 12/02/2022] [Indexed: 12/11/2022] Open
Abstract
MOTIVATION Whole-genome sequencing (WGS) is increasingly used to aid the understanding of Mycobacterium tuberculosis (MTB) transmission. The epidemiological analysis of tuberculosis based on the WGS technique requires a diverse collection of bioinformatics tools. Effectively using these analysis tools in a scalable and reproducible way can be challenging, especially for non-experts. RESULTS Here, we present TransFlow (Transmission Workflow), a user-friendly, fast, efficient and comprehensive WGS-based transmission analysis pipeline. TransFlow combines some state-of-the-art tools to take transmission analysis from raw sequencing data, through quality control, sequence alignment and variant calling, into downstream transmission clustering, transmission network reconstruction and transmission risk factor inference, together with summary statistics and data visualization in a summary report. TransFlow relies on Snakemake and Conda to resolve dependencies among consecutive processing steps and can be easily adapted to any computation environment. AVAILABILITY AND IMPLEMENTATION TransFlow is free available at https://github.com/cvn001/transflow. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
| | | | - Mingwu Zhang
- The Institute of TB Control, Zhejiang Provincial Center for Disease Control and Prevention, Hangzhou, Zhejiang 310051, China
| | - Yewei Lu
- Key Laboratory of Precision Medicine in Diagnosis and Monitoring Research of Zhejiang Province, Hangzhou, Zhejiang 310020, China
| | - Yelei Zhu
- The Institute of TB Control, Zhejiang Provincial Center for Disease Control and Prevention, Hangzhou, Zhejiang 310051, China
| | - Kunyang Wu
- The Institute of TB Control, Zhejiang Provincial Center for Disease Control and Prevention, Hangzhou, Zhejiang 310051, China
| | - Yiwen Wu
- Department of Medical Oncology, Zhejiang Chinese Medical University, Hangzhou, Zhejiang 310053, China
| | - Weixin Wang
- Key Laboratory of Precision Medicine in Diagnosis and Monitoring Research of Zhejiang Province, Hangzhou, Zhejiang 310020, China
| | - Bin Chen
- The Institute of TB Control, Zhejiang Provincial Center for Disease Control and Prevention, Hangzhou, Zhejiang 310051, China
| | - Zhengwei Liu
- To whom correspondence should be addressed. or or
| | | | - Junshun Gao
- To whom correspondence should be addressed. or or
| |
Collapse
|
30
|
Goldstein IH, Bayer D, Barilar I, Kizito B, Matsiri O, Modongo C, Zetola NM, Niemann S, Minin VM, Shin SS. Using genetic data to identify transmission risk factors: Statistical assessment and application to tuberculosis transmission. PLoS Comput Biol 2022; 18:e1010696. [PMID: 36469509 PMCID: PMC9754595 DOI: 10.1371/journal.pcbi.1010696] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2022] [Revised: 12/15/2022] [Accepted: 10/31/2022] [Indexed: 12/12/2022] Open
Abstract
Identifying host factors that influence infectious disease transmission is an important step toward developing interventions to reduce disease incidence. Recent advances in methods for reconstructing infectious disease transmission events using pathogen genomic and epidemiological data open the door for investigation of host factors that affect onward transmission. While most transmission reconstruction methods are designed to work with densely sampled outbreaks, these methods are making their way into surveillance studies, where the fraction of sampled cases with sequenced pathogens could be relatively low. Surveillance studies that use transmission event reconstruction then use the reconstructed events as response variables (i.e., infection source status of each sampled case) and use host characteristics as predictors (e.g., presence of HIV infection) in regression models. We use simulations to study estimation of the effect of a host factor on probability of being an infection source via this multi-step inferential procedure. Using TransPhylo-a widely-used method for Bayesian estimation of infectious disease transmission events-and logistic regression, we find that low sensitivity of identifying infection sources leads to dilution of the signal, biasing logistic regression coefficients toward zero. We show that increasing the proportion of sampled cases improves sensitivity and some, but not all properties of the logistic regression inference. Application of these approaches to real world data from a population-based TB study in Botswana fails to detect an association between HIV infection and probability of being a TB infection source. We conclude that application of a pipeline, where one first uses TransPhylo and sparsely sampled surveillance data to infer transmission events and then estimates effects of host characteristics on probabilities of these events, should be accompanied by a realistic simulation study to better understand biases stemming from imprecise transmission event inference.
Collapse
Affiliation(s)
- Isaac H. Goldstein
- Department of Statistics, University of California, Irvine, California, United States of America
| | - Damon Bayer
- Department of Statistics, University of California, Irvine, California, United States of America
| | - Ivan Barilar
- German Center for Infection Research, Research Center Borstel, Borstel, Germany
| | | | | | | | | | - Stefan Niemann
- German Center for Infection Research, Research Center Borstel, Borstel, Germany
| | - Volodymyr M. Minin
- Department of Statistics, University of California, Irvine, California, United States of America
| | - Sanghyuk S. Shin
- Sue & Bill Gross School of Nursing, University of California, Irvine, California, United States of America
| |
Collapse
|
31
|
van Doren TP. Biocultural perspectives of infectious diseases and demographic evolution: Tuberculosis and its comorbidities through history. Evol Anthropol 2022; 32:100-117. [PMID: 36436141 DOI: 10.1002/evan.21970] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2021] [Revised: 06/09/2022] [Accepted: 11/11/2022] [Indexed: 11/28/2022]
Abstract
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.
Collapse
|
32
|
Guellil M, van Dorp L, Inskip SA, Dittmar JM, Saag L, Tambets K, Hui R, Rose A, D’Atanasio E, Kriiska A, Varul L, Koekkelkoren AMHC, Goldina RD, Cessford C, Solnik A, Metspalu M, Krause J, Herbig A, Robb JE, Houldcroft CJ, Scheib CL. Ancient herpes simplex 1 genomes reveal recent viral structure in Eurasia. SCIENCE ADVANCES 2022; 8:eabo4435. [PMID: 35895820 PMCID: PMC9328674 DOI: 10.1126/sciadv.abo4435] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/03/2022] [Accepted: 06/10/2022] [Indexed: 05/05/2023]
Abstract
Human herpes simplex virus 1 (HSV-1), a life-long infection spread by oral contact, infects a majority of adults globally. Phylogeographic clustering of sampled diversity into European, pan-Eurasian, and African groups has suggested the virus codiverged with human migrations out of Africa, although a much younger origin has also been proposed. We present three full ancient European HSV-1 genomes and one partial genome, dating from the 3rd to 17th century CE, sequenced to up to 9.5× with paired human genomes up to 10.16×. Considering a dataset of modern and ancient genomes, we apply phylogenetic methods to estimate the age of sampled modern Eurasian HSV-1 diversity to 4.68 (3.87 to 5.65) ka. Extrapolation of estimated rates to a global dataset points to the age of extant sampled HSV-1 as 5.29 (4.60 to 6.12) ka, suggesting HSV-1 lineage replacement coinciding with the late Neolithic period and following Bronze Age migrations.
Collapse
Affiliation(s)
- Meriam Guellil
- Estonian Biocentre, Institute of Genomics, University of Tartu, Riia 23B, Tartu 51010, Estonia
| | - Lucy van Dorp
- UCL Genetics Institute, Department of Genetics, Evolution, and Environment, University College London, London WC1E 6BT, UK
| | - Sarah A. Inskip
- McDonald Institute for Archaeological Research, University of Cambridge, Cambridge, UK
- Department of Archaeology and Ancient History, University of Leicester, University Road, Leicester, LE1 7RH, UK
| | - Jenna M. Dittmar
- McDonald Institute for Archaeological Research, University of Cambridge, Cambridge, UK
- Department of Archaeology, University of Aberdeen, UK
| | - Lehti Saag
- Estonian Biocentre, Institute of Genomics, University of Tartu, Riia 23B, Tartu 51010, Estonia
- UCL Genetics Institute, Department of Genetics, Evolution, and Environment, University College London, London WC1E 6BT, UK
| | - Kristiina Tambets
- Estonian Biocentre, Institute of Genomics, University of Tartu, Riia 23B, Tartu 51010, Estonia
| | - Ruoyun Hui
- McDonald Institute for Archaeological Research, University of Cambridge, Cambridge, UK
- Alan Turing Institute, 2QR, John Dodson House, 96 Euston Rd., London NW1 2DB, UK
| | - Alice Rose
- McDonald Institute for Archaeological Research, University of Cambridge, Cambridge, UK
| | | | - Aivar Kriiska
- Department of Archaeology, Institute of History and Archaeology, University of Tartu, Tartu 51014, Estonia
| | - Liivi Varul
- Archaeological Research Collection, School of Humanities, Tallinn University, Tallinn 10130, Estonia
| | | | - Rimma D. Goldina
- Department History of Udmurtia, Archaeology and Ethnology, Udmurt State University, 1, Universitetskaya St. 1, 426034 Izhevsk, Russia
| | - Craig Cessford
- Cambridge Archaeological Unit, Department of Archaeology, University of Cambridge, Cambridge, UK
| | - Anu Solnik
- Core Facility, Institute of Genomics, University of Tartu, Riia 23B, Tartu 51010 Estonia
| | - Mait Metspalu
- Estonian Biocentre, Institute of Genomics, University of Tartu, Riia 23B, Tartu 51010, Estonia
| | - Johannes Krause
- Department of Archaeogenetics, Max Planck Institute for the Science of Human History, Jena, Germany
- Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany
| | - Alexander Herbig
- Department of Archaeogenetics, Max Planck Institute for the Science of Human History, Jena, Germany
- Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany
| | - John E. Robb
- Department of Archaeology, University of Cambridge, Cambridge, UK
| | | | - Christiana L. Scheib
- Estonian Biocentre, Institute of Genomics, University of Tartu, Riia 23B, Tartu 51010, Estonia
- St. John’s College, University of Cambridge, Cambridge, CB2 1TP, UK
| |
Collapse
|
33
|
Saavedra Cervera B, López MG, Chiner-Oms Á, García AM, Cancino-Muñoz I, Torres-Puente M, Villamayor L, Madrazo-Moya C, Mambuque E, Sequera GV, Respeito D, Blanco S, Augusto O, López-Varela E, García-Basteiro AL, Comas I. Fine-grain population structure and transmission patterns of Mycobacterium tuberculosis in southern Mozambique, a high TB/HIV burden area. Microb Genom 2022; 8. [PMID: 35787782 PMCID: PMC9455694 DOI: 10.1099/mgen.0.000844] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
Genomic studies of the Mycobacterium tuberculosis complex (MTBC) might shed light on the dynamics of its transmission, especially in high-burden settings, where recent outbreaks are embedded in the complex natural history of the disease. To this end, we conducted a 1 year prospective surveillance-based study in Mozambique. We applied whole-genome sequencing (WGS) to 295 positive cultures. We fully characterized MTBC isolates by phylogenetics and dating evaluation, and carried out a molecular epidemiology analysis to investigate further associations with pre-defined transmission risk factors. The majority of strains (49.5%, 136/275) belonged to lineage (L) 4; 57.8 % of them (159/275) were in genomic transmission clusters (cut-off 5 SNPs), and a strikingly high proportion (45.5%) shared an identical genotype (0 SNP pairwise distance). We found two ‘likely endemic’ clades, comprising 67 strains, belonging to L1.2, which dated back to the late 19th century and were associated with recent spread among people living with human immunodeficiency virus (PLHIV). We describe for the first time the population structure of MTBC in our region, a high tuberculosis (TB)/HIV burden area. Clustering analysis revealed an unforeseen pattern of spread and high rates of progression to active TB, suggesting weaknesses in TB control activities. The long-term presence of local strains in Mozambique, which were responsible for large transmission among HIV/TB-coinfected patients, calls into question the role of HIV in TB transmission.
Collapse
Affiliation(s)
- Belén Saavedra Cervera
- PhD Programin Medicine and Translational Research, Universitat de Barcelona, Barcelona, Spain.,Centro de Investigação em Saúde de Manhiça (CISM), Maputo, Mozambique.,ISGlobal, Hospital Clínic - Universitat de Barcelona, Barcelona, Spain
| | - Mariana G López
- Instituto de Biomedicina de Valencia (IBV), CSIC, Valencia, Spain
| | | | - Ana María García
- Instituto de Biomedicina de Valencia (IBV), CSIC, Valencia, Spain.,Universidad de Valencia, Valencia, Spain
| | | | | | | | | | - Edson Mambuque
- Centro de Investigação em Saúde de Manhiça (CISM), Maputo, Mozambique
| | | | - Durval Respeito
- Centro de Investigação em Saúde de Manhiça (CISM), Maputo, Mozambique
| | - Silvia Blanco
- Centro de Investigação em Saúde de Manhiça (CISM), Maputo, Mozambique
| | - Orvalho Augusto
- Centro de Investigação em Saúde de Manhiça (CISM), Maputo, Mozambique
| | - Elisa López-Varela
- Centro de Investigação em Saúde de Manhiça (CISM), Maputo, Mozambique.,ISGlobal, Hospital Clínic - Universitat de Barcelona, Barcelona, Spain
| | - Alberto L García-Basteiro
- Centro de Investigação em Saúde de Manhiça (CISM), Maputo, Mozambique.,ISGlobal, Hospital Clínic - Universitat de Barcelona, Barcelona, Spain.,Centro de Investigación Biomédica en Red de Enfermedades Infecciosas (CIBERINFEC), Barcelona, Spain
| | - Iñaki Comas
- Instituto de Biomedicina de Valencia (IBV), CSIC, Valencia, Spain.,CIBER in Epidemiology and Public Health, Madrid, Spain
| |
Collapse
|
34
|
Menardo F. Understanding drivers of phylogenetic clustering and terminal branch lengths distribution in epidemics of Mycobacterium tuberculosis. eLife 2022; 11:76780. [PMID: 35762734 PMCID: PMC9239681 DOI: 10.7554/elife.76780] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2022] [Accepted: 06/15/2022] [Indexed: 11/13/2022] Open
Abstract
Detecting factors associated with transmission is important to understand disease epidemics, and to design effective public health measures. Clustering and terminal branch lengths (TBL) analyses are commonly applied to genomic data sets of Mycobacterium tuberculosis (MTB) to identify sub-populations with increased transmission. Here, I used a simulation-based approach to investigate what epidemiological processes influence the results of clustering and TBL analyses, and whether differences in transmission can be detected with these methods. I simulated MTB epidemics with different dynamics (latency, infectious period, transmission rate, basic reproductive number R0, sampling proportion, sampling period, and molecular clock), and found that all considered factors, except for the length of the infectious period, affect the results of clustering and TBL distributions. I show that standard interpretations of this type of analyses ignore two main caveats: (1) clustering results and TBL depend on many factors that have nothing to do with transmission, (2) clustering results and TBL do not tell anything about whether the epidemic is stable, growing, or shrinking, unless all the additional parameters that influence these metrics are known, or assumed identical between sub-populations. An important consequence is that the optimal SNP threshold for clustering depends on the epidemiological conditions, and that sub-populations with different epidemiological characteristics should not be analyzed with the same threshold. Finally, these results suggest that different clustering rates and TBL distributions, that are found consistently between different MTB lineages, are probably due to intrinsic bacterial factors, and do not indicate necessarily differences in transmission or evolutionary success.
Collapse
Affiliation(s)
- Fabrizio Menardo
- Department of Plant and Microbial Biology, University of Zurich, Zurich, Switzerland
| |
Collapse
|
35
|
Swargam S, Kumari I, Kumar A, Pradhan D, Alam A, Singh H, Jain A, Devi KR, Trivedi V, Sarma J, Hanif M, Narain K, Ehtesham NZ, Hasnain SE, Ahmad S. MycoVarP: Mycobacterium Variant and Drug Resistance Prediction Pipeline for Whole-Genome Sequence Data Analysis. FRONTIERS IN BIOINFORMATICS 2022; 1:805338. [PMID: 36303799 PMCID: PMC9580932 DOI: 10.3389/fbinf.2021.805338] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2021] [Accepted: 12/13/2021] [Indexed: 11/13/2022] Open
Abstract
Whole-genome sequencing (WGS) provides a comprehensive tool to analyze the bacterial genomes for genotype–phenotype correlations, diversity of single-nucleotide variant (SNV), and their evolution and transmission. Several online pipelines and standalone tools are available for WGS analysis of Mycobacterium tuberculosis (Mtb) complex (MTBC). While they facilitate the processing of WGS data with minimal user expertise, they are either too general, providing little insights into bacterium-specific issues such as gene variations, INDEL/synonymous/PE-PPE (IDP family), and drug resistance from sample data, or are limited to specific objectives, such as drug resistance. It is understood that drug resistance and lineage-specific issues require an elaborate prioritization of identified variants to choose the best target for subsequent therapeutic intervention. Mycobacterium variant pipeline (MycoVarP) addresses these specific issues with a flexible battery of user-defined and default filters. It provides an end-to-end solution for WGS analysis of Mtb variants from the raw reads and performs two quality checks, viz, before trimming and after alignments of reads to the reference genome. MycoVarP maps the annotated variants to the drug-susceptible (DS) database and removes the false-positive variants, provides lineage identification, and predicts potential drug resistance. We have re-analyzed the WGS data reported by Advani et al. (2019) using MycoVarP and identified some additional variants not reported so far. We conclude that MycoVarP will help in identifying nonsynonymous, true-positive, drug resistance–associated variants more effectively and comprehensively, including those within the IDP of the PE-PPE/PGRS family, than possible from the currently available pipelines.
Collapse
Affiliation(s)
- Sandeep Swargam
- Department of Biochemical Engineering and Biotechnology, Indian Institute of Technology, Hauz Khas, New Delhi, India
- Department of Molecular Medicine, School of Interdisciplinary Sciences, Jamia Hamdard, New Delhi, India
| | - Indu Kumari
- Inflammation Biology and Cell Signalling Lab, Safdarjung Hospital Campus, ICMR National Institute of Pathology, New Delhi, India
| | - Amit Kumar
- ICMR Computational Genomics Centre, Informatics Systems and Research Management (ISRM) Division, Indian Council of Medical Research (ICMR), New Delhi, India
| | - Dibyabhaba Pradhan
- ICMR Computational Genomics Centre, Informatics Systems and Research Management (ISRM) Division, Indian Council of Medical Research (ICMR), New Delhi, India
| | - Anwar Alam
- Inflammation Biology and Cell Signalling Lab, Safdarjung Hospital Campus, ICMR National Institute of Pathology, New Delhi, India
| | - Harpreet Singh
- ICMR Computational Genomics Centre, Informatics Systems and Research Management (ISRM) Division, Indian Council of Medical Research (ICMR), New Delhi, India
| | - Anuja Jain
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| | | | - Vishal Trivedi
- Department of Biosciences and Bioengineering, Indian Institute of Technology-Guwahati, Guwahati, India
| | - Jogesh Sarma
- Department of Pulmonary Medicine, Guwahati, India
| | | | - Kanwar Narain
- ICMR-Regional Medical Research Centre, Dibrugarh, India
| | - Nasreen Zafar Ehtesham
- Inflammation Biology and Cell Signalling Lab, Safdarjung Hospital Campus, ICMR National Institute of Pathology, New Delhi, India
| | - Seyed Ehtesham Hasnain
- Department of Biochemical Engineering and Biotechnology, Indian Institute of Technology, Hauz Khas, New Delhi, India
- Department of Life Sciences, Sharda University, Greater NOIDA, India
| | - Shandar Ahmad
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| |
Collapse
|
36
|
Ceres KM, Stanhope MJ, Gröhn YT. A critical evaluation of Mycobacterium bovis pangenomics, with reference to its utility in outbreak investigation. Microb Genom 2022; 8:mgen000839. [PMID: 35763423 PMCID: PMC9455707 DOI: 10.1099/mgen.0.000839] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2021] [Accepted: 04/29/2022] [Indexed: 11/18/2022] Open
Abstract
The increased accessibility of next generation sequencing has allowed enough genomes from a given bacterial species to be sequenced to describe the distribution of genes in the pangenome, without limiting analyses to genes present in reference strains. Although some taxa have thousands of whole genome sequences available on public databases, most genomes were sequenced with short read technology, resulting in incomplete assemblies. Studying pangenomes could lead to important insights into adaptation, pathogenicity, or molecular epidemiology, however given the known information loss inherent in analyzing contig-level assemblies, these inferences may be biased or inaccurate. In this study we describe the pangenome of a clonally evolving pathogen, Mycobacterium bovis , and examine the utility of gene content variation in M. bovis outbreak investigation. We constructed the M. bovis pangenome using 1463 de novo assembled genomes. We tested the assumption of strict clonal evolution by studying evidence of recombination in core genes and analyzing the distribution of accessory genes among core monophyletic groups. To determine if gene content variation could be utilized in outbreak investigation, we carefully examined accessory genes detected in a well described M. bovis outbreak in Minnesota. We found significant errors in accessory gene classification. After accounting for these errors, we show that M. bovis has a much smaller accessory genome than previously described and provide evidence supporting ongoing clonal evolution and a closed pangenome, with little gene content variation generated over outbreaks. We also identified frameshift mutations in multiple genes, including a mutation in glpK , which has recently been associated with antibiotic tolerance in Mycobacterium tuberculosis . A pangenomic approach enables a more comprehensive analysis of genome dynamics than is possible with reference-based approaches; however, without critical evaluation of accessory gene content, inferences of transmission patterns employing these loci could be misguided.
Collapse
Affiliation(s)
- Kristina M. Ceres
- Department of Population Medicine and Diagnostic Sciences, College of Veterinary Medicine, Cornell University, Ithaca, New York, USA
- Population and Ecosystem Health, College of Veterinary Medicine, Cornell University, Ithaca, NY, USA
| | - Michael J Stanhope
- Department of Population Medicine and Diagnostic Sciences, College of Veterinary Medicine, Cornell University, Ithaca, New York, USA
- Population and Ecosystem Health, College of Veterinary Medicine, Cornell University, Ithaca, NY, USA
| | - Yrjö T. Gröhn
- Department of Population Medicine and Diagnostic Sciences, College of Veterinary Medicine, Cornell University, Ithaca, New York, USA
- Population and Ecosystem Health, College of Veterinary Medicine, Cornell University, Ithaca, NY, USA
| |
Collapse
|
37
|
Bainomugisa A, Lavu E, Pandey S, Majumdar S, Banamu J, Coulter C, Marais B, Coin L, Graham SM, du Cros P. Evolution and spread of a highly drug resistant strain of Mycobacterium tuberculosis in Papua New Guinea. BMC Infect Dis 2022; 22:437. [PMID: 35524232 PMCID: PMC9077924 DOI: 10.1186/s12879-022-07414-2] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2021] [Accepted: 04/12/2022] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Molecular mechanisms determining the transmission and prevalence of drug resistant tuberculosis (DR-TB) in Papua New Guinea (PNG) are poorly understood. We used genomic and drug susceptibility data to explore the evolutionary history, temporal acquisition of resistance and transmission dynamics of DR-TB across PNG. METHODS We performed whole genome sequencing on isolates from Central Public Health Laboratory, PNG, collected 2017-2019. Data analysis was done on a composite dataset that also included 100 genomes previously sequenced from Daru, PNG (2012-2015). RESULTS Sampled isolates represented 14 of the 22 PNG provinces, the majority (66/94; 70%) came from the National Capital District (NCD). In the composite dataset, 91% of strains were Beijing 2.2.1.1, identified in 13 provinces. Phylogenetic tree of Beijing strains revealed two clades, Daru dominant clade (A) and NCD dominant clade (B). Multi-drug resistance (MDR) was repeatedly and independently acquired, with the first MDR cases in both clades noted to have emerged in the early 1990s, while fluoroquinolone resistance emerged in 2009 (95% highest posterior density 2000-2016). We identified the presence of a frameshift mutation within Rv0678 (p.Asp47fs) which has been suggested to confer resistance to bedaquiline, despite no known exposure to the drug. Overall genomic clustering was significantly associated with rpoC compensatory and inhA promoter mutations (p < 0.001), with high percentage of most genomic clusters (12/14) identified in NCD, reflecting its role as a potential national amplifier. CONCLUSIONS The acquisition and evolution of drug resistance among the major clades of Beijing strain threaten the success of DR-TB treatment in PNG. With continued transmission of this strain in PNG, genotypic drug resistance surveillance using whole genome sequencing is essential for improved public health response to outbreaks. With occurrence of resistance to newer drugs such as bedaquiline, knowledge of full drug resistance profiles will be important for optimal treatment selection.
Collapse
Affiliation(s)
| | - Evelyn Lavu
- University of Papua New Guinea, Port Moresby, Papua New Guinea.,Central Public Health Laboratory, Port Moresby, Papua New Guinea
| | - Sushil Pandey
- Queensland Mycobacteria Reference Laboratory, Brisbane, QLD, Australia
| | - Suman Majumdar
- Burnet Institute, 85 Commercial Road, Melbourne, VIC, 3004, Australia.,University of Melbourne Department of Paediatrics and Murdoch Childrens Research Institute, Royal Children's Hospital, Melbourne, VIC, Australia
| | - Jennifer Banamu
- Central Public Health Laboratory, Port Moresby, Papua New Guinea
| | - Chris Coulter
- Queensland Mycobacteria Reference Laboratory, Brisbane, QLD, Australia
| | - Ben Marais
- University of Sydney, Sydney, NSW, Australia
| | - Lachlan Coin
- Peter Doherty Institute, Melbourne, VIC, Australia
| | - Stephen M Graham
- Burnet Institute, 85 Commercial Road, Melbourne, VIC, 3004, Australia.,University of Melbourne Department of Paediatrics and Murdoch Childrens Research Institute, Royal Children's Hospital, Melbourne, VIC, Australia
| | - Philipp du Cros
- Burnet Institute, 85 Commercial Road, Melbourne, VIC, 3004, Australia.
| |
Collapse
|
38
|
Walter KS, Dos Santos PCP, Gonçalves TO, da Silva BO, da Silva Santos A, de Cássia Leite A, da Silva AM, Figueira Moreira FM, de Oliveira RD, Lemos EF, Cunha E, Liu YE, Ko AI, Colijn C, Cohen T, Mathema B, Croda J, Andrews JR. The role of prisons in disseminating tuberculosis in Brazil: A genomic epidemiology study. LANCET REGIONAL HEALTH. AMERICAS 2022; 9. [PMID: 35647574 PMCID: PMC9140320 DOI: 10.1016/j.lana.2022.100186] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
Background Globally, prisons are high-incidence settings for tuberculosis. Yet the role of prisons as reservoirs of M. tuberculosis, propagating epidemics through spillover to surrounding communities, has been difficult to measure directly. Methods To quantify the role of prisons in driving wider community M. tuberculosis transmission, we conducted prospective genomic surveillance in Central West Brazil from 2014 to 2019. We whole genome sequenced 1152 M. tuberculosis isolates collected during active and passive surveillance inside and outside prisons and linked genomes to detailed incarceration histories. We applied multiple phylogenetic and genomic clustering approaches and inferred timed transmission trees. Findings M. tuberculosis sequences from incarcerated and non-incarcerated people were closely related in a maximum likelihood phylogeny. The majority (70.8%; 46/65) of genomic clusters including people with no incarceration history also included individuals with a recent history of incarceration. Among cases in individuals with no incarceration history, 50.6% (162/320) were in clusters that included individuals with recent incarceration history, suggesting that transmission chains often span prisons and communities. We identified a minimum of 18 highly probable spillover events, M. tuberculosis transmission from people with a recent incarceration history to people with no prior history of incarceration, occurring in the state’s four largest cities and across sampling years. We additionally found that frequent transfers of people between the state’s prisons creates a highly connected prison network that likely disseminates M. tuberculosis across the state. Interpretation We developed a framework for measuring spillover from high-incidence environments to surrounding communities by integrating genomic and spatial information. Our findings indicate that, in this setting, prisons serve not only as disease reservoirs, but also disseminate M. tuberculosis across highly connected prison networks, both amplifying and propagating M. tuberculosis risk in surrounding communities. Funding Brazil’s National Council for Scientific and Technological Development and US National Institutes of Health.
Collapse
Affiliation(s)
- Katharine S Walter
- Division of Infectious Diseases and Geographic Medicine, Stanford University School of Medicine, Stanford, CA 94305, United States
| | | | | | - Bruna Oliveira da Silva
- Health Sciences Research Laboratory, Federal University of Grande Dourados, Dourados, Brazil
| | - Andrea da Silva Santos
- Health Sciences Research Laboratory, Federal University of Grande Dourados, Dourados, Brazil
| | | | - Alessandra Moura da Silva
- School of Medicine, Federal University of Mato Grosso do Sul, School of Medicine, Campo Grande, Brazil
| | | | | | - Everton Ferreira Lemos
- School of Medicine, Federal University of Mato Grosso do Sul, School of Medicine, Campo Grande, Brazil
| | - Eunice Cunha
- Laboratory of Bacteriology, Central Laboratory of Mato Grosso do Sul, Campo Grande, Brazil
| | - Yiran E Liu
- Division of Infectious Diseases and Geographic Medicine, Stanford University School of Medicine, Stanford, CA 94305, United States.,Cancer Biology Graduate Program, Stanford University School of Medicine, Stanford, United States
| | - Albert I Ko
- Department of Epidemiology of Microbial Diseases, Yale School of Public Health, New Haven, United States.,Instituto Gonçalo¸ Moniz, Fundação Oswaldo Cruz, Salvador, BA, Brazil
| | - Caroline Colijn
- Department of Mathematics, Simon Fraser University, Burnaby, Canada
| | - Ted Cohen
- Department of Epidemiology of Microbial Diseases, Yale School of Public Health, New Haven, United States
| | - Barun Mathema
- Department of Epidemiology, Columbia University Mailman School of Public Health, New York, United States
| | - Julio Croda
- School of Medicine, Federal University of Mato Grosso do Sul, School of Medicine, Campo Grande, Brazil.,Department of Epidemiology of Microbial Diseases, Yale School of Public Health, New Haven, United States.,Mato Grosso do Sul Office, Oswaldo Cruz Foundation, Campo Grande, Brazil
| | - Jason R Andrews
- Division of Infectious Diseases and Geographic Medicine, Stanford University School of Medicine, Stanford, CA 94305, United States
| |
Collapse
|
39
|
Evidence for local and international spread of Mycobacterium avium subspecies paratuberculosis through whole genome sequencing of isolates from the island of Ireland. Vet Microbiol 2022; 268:109416. [DOI: 10.1016/j.vetmic.2022.109416] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2021] [Revised: 01/14/2022] [Accepted: 04/01/2022] [Indexed: 12/18/2022]
|
40
|
Chiner-Oms Á, López MG, Moreno-Molina M, Furió V, Comas I. Gene evolutionary trajectories in Mycobacterium tuberculosis reveal temporal signs of selection. Proc Natl Acad Sci U S A 2022; 119:e2113600119. [PMID: 35452305 PMCID: PMC9173582 DOI: 10.1073/pnas.2113600119] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2021] [Accepted: 02/17/2022] [Indexed: 12/20/2022] Open
Abstract
Genetic differences between different Mycobacterium tuberculosis complex (MTBC) strains determine their ability to transmit within different host populations, their latency times, and their drug resistance profiles. Said differences usually emerge through de novo mutations and are maintained or discarded by the balance of evolutionary forces. Using a dataset of ∼5,000 strains representing global MTBC diversity, we determined the past and present selective forces that have shaped the current variability observed in the pathogen population. We identified regions that have evolved under changing types of selection since the time of the MTBC common ancestor. Our approach highlighted striking differences in the genome regions relevant for host–pathogen interaction and, in particular, suggested an adaptive role for the sensor protein of two-component systems. In addition, we applied our approach to successfully identify potential determinants of resistance to drugs administered as second-line tuberculosis treatments.
Collapse
Affiliation(s)
- Álvaro Chiner-Oms
- Instituto de Biomedicina de Valencia (IBV-CSIC), Valencia, 46010, Spain
| | - Mariana G. López
- Instituto de Biomedicina de Valencia (IBV-CSIC), Valencia, 46010, Spain
| | | | - Victoria Furió
- Instituto de Biomedicina de Valencia (IBV-CSIC), Valencia, 46010, Spain
| | - Iñaki Comas
- Instituto de Biomedicina de Valencia (IBV-CSIC), Valencia, 46010, Spain
- CIBER en Epidemiología y Salud Pública, Valencia, Spain
| |
Collapse
|
41
|
Cardona PJ, Català M, Prats C. The Origin and Maintenance of Tuberculosis Is Explained by the Induction of Smear-Negative Disease in the Paleolithic. Pathogens 2022; 11:pathogens11030366. [PMID: 35335692 PMCID: PMC8955375 DOI: 10.3390/pathogens11030366] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2022] [Revised: 03/05/2022] [Accepted: 03/15/2022] [Indexed: 12/05/2022] Open
Abstract
Is it possible that the origin of Mycobacterium tuberculosis (Mtb) infection was around 70,000 years before the common era? At that time Homo sapiens was just another primate species with discrete growth and a very low-density geographic occupation. Therefore, it is difficult to understand the origin of a highly virulent obligate human pathogen. We have designed a new SEIR model (TBSpectr) that allows the differentiation of smear-positive and -negative tuberculosis. The model reconciles currently accepted growth rates for the Middle Paleolithic (0.003%/year) and Neolithic (0.1%/year). The obtained data link the origin of Mtb infection in the Middle Paleolithic to the induction of smear-negative TB, and reveal that its persistence required interrelations among hunter–gatherer groups, while the risk of human extinction was negligible. It also highlights the number of people infected per case and the fast progression to disease for Mtb infection maintenance, as well as the link between poor health in the Neolithic with the increased incidence of more severe forms of TB (smear-positive). In conclusion, our data support the origin of TB as a well-tolerated, highly persistent disease, even in low-density populations, showing the difficulty of its eradication and highlighting the necessity for providing better health conditions to humans to reduce its severity.
Collapse
Affiliation(s)
- Pere-Joan Cardona
- Unitat de Tuberculosi Experimental, Germans Trias i Pujol Research Institute (IGTP), 08916 Badalona, Spain
- Microbiology Department, North Metropolitan Clinical Laboratory, ‘Germans Trias i Pujol’ University Hospital, 08916 Badalona, Spain
- Genetics and Microbiology Department, Universitat Autònoma de Barcelona, 08916 Barcelona, Spain
- Centro de Investigación Biomédica en Red de Enfermedades Respiratorias (CIBERES), 28029 Madrid, Spain
- Correspondence:
| | - Martí Català
- Comparative Medicine and Bioimage Centre of Catalonia (CMCiB), Germans Trias i Pujol Research Institute (IGTP), 08916 Badalona, Spain; (M.C.); (C.P.)
- Departament de Física, Escola d’Enginyeria Agroalimentària i de Biosistemes de Barcelona, Universitat Politècnica de Catalunya (UPC)-BarcelonaTech, 08916 Badalona, Spain
| | - Clara Prats
- Comparative Medicine and Bioimage Centre of Catalonia (CMCiB), Germans Trias i Pujol Research Institute (IGTP), 08916 Badalona, Spain; (M.C.); (C.P.)
- Departament de Física, Escola d’Enginyeria Agroalimentària i de Biosistemes de Barcelona, Universitat Politècnica de Catalunya (UPC)-BarcelonaTech, 08916 Badalona, Spain
| |
Collapse
|
42
|
Wollenberg KR, Jeffrey BM, Harris MA, Gabrielian A, Hurt DE, Rosenthal A. Patterns of genomic interrelatedness of publicly available samples in the TB portals database. Tuberculosis (Edinb) 2022; 133:102171. [PMID: 35101846 PMCID: PMC8997244 DOI: 10.1016/j.tube.2022.102171] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2021] [Revised: 01/18/2022] [Accepted: 01/23/2022] [Indexed: 10/19/2022]
Abstract
The TB Portals program is an international collaboration for the collection and dissemination of tuberculosis data from patient cases focused on drug resistance. The central database is a patient-oriented resource containing both patient and pathogen clinical and genomic information. Herein we provide a summary of the pathogen genomic data available through the TB Portals and show one potential application by examining patterns of genomic pairwise distances. Distributions of pairwise distances highlight overall patterns of genome variability within and between Mycobacterium tuberculosis phylogenomic lineages. Closely related isolates (based on whole-genome pairwise distances and time between sample collection dates) from different countries were identified as potential evidence of international transmission of drug-resistant tuberculosis. These high-level views of genomic relatedness provide information that can stimulate hypotheses for further and more detailed research.
Collapse
Affiliation(s)
- Kurt R. Wollenberg
- Office of Cyber Infrastructure and Computational Biology, National Institute of Allergy and Infectious Disease, National Institutes of Health, Bethesda, Maryland, USA
| | - Brendan M. Jeffrey
- Office of Cyber Infrastructure and Computational Biology, National Institute of Allergy and Infectious Disease, National Institutes of Health, Bethesda, Maryland, USA
| | - Michael A. Harris
- Office of Cyber Infrastructure and Computational Biology, National Institute of Allergy and Infectious Disease, National Institutes of Health, Bethesda, Maryland, USA,To whom correspondence should be addressed: . Telephone: 301-761-6746
| | - Andrei Gabrielian
- Office of Cyber Infrastructure and Computational Biology, National Institute of Allergy and Infectious Disease, National Institutes of Health, Bethesda, MD, USA.
| | - Darrell E. Hurt
- Office of Cyber Infrastructure and Computational Biology, National Institute of Allergy and Infectious Disease, National Institutes of Health, Bethesda, Maryland, USA
| | - Alex Rosenthal
- Office of Cyber Infrastructure and Computational Biology, National Institute of Allergy and Infectious Disease, National Institutes of Health, Bethesda, MD, USA.
| |
Collapse
|
43
|
Malyarchuk AB, Andreeva TV, Kuznetsova IL, Kunizheva SS, Protasova MS, Uralsky LI, Tyazhelova TV, Gusev FE, Manakhov AD, Rogaev EI. Genomics of Ancient Pathogens: First Advances and Prospects. BIOCHEMISTRY. BIOKHIMIIA 2022; 87:242-258. [PMID: 35526849 PMCID: PMC8916790 DOI: 10.1134/s0006297922030051] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/14/2021] [Revised: 01/08/2022] [Accepted: 01/21/2022] [Indexed: 11/23/2022]
Abstract
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: Yersinia pestis (plague), Variola virus (smallpox), Vibrio cholerae (cholera), HBV (hepatitis B virus), as well as the equally important endemic human infectious agents: Mycobacterium tuberculosis (tuberculosis), Mycobacterium leprae (leprosy), and Treponema pallidum (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.
Collapse
Affiliation(s)
- Alexandra B Malyarchuk
- Center for Genetics and Genetic Technologies, Faculty of Biology, Lomonosov Moscow State University, Moscow, 119234, Russia.
| | - Tatiana V Andreeva
- Center for Genetics and Genetic Technologies, Faculty of Biology, Lomonosov Moscow State University, Moscow, 119234, Russia
- Vavilov Institute of General Genetics, Russian Academy of Sciences, Moscow, 119333, Russia
| | - Irina L Kuznetsova
- Vavilov Institute of General Genetics, Russian Academy of Sciences, Moscow, 119333, Russia
- Center for Genetics and Life Science, Sirius University of Science and Technology, Sochi, 354340, Russia
| | - Svetlana S Kunizheva
- Vavilov Institute of General Genetics, Russian Academy of Sciences, Moscow, 119333, Russia
- Center for Genetics and Life Science, Sirius University of Science and Technology, Sochi, 354340, Russia
| | - Maria S Protasova
- Vavilov Institute of General Genetics, Russian Academy of Sciences, Moscow, 119333, Russia
| | - Lev I Uralsky
- Vavilov Institute of General Genetics, Russian Academy of Sciences, Moscow, 119333, Russia
- Center for Genetics and Life Science, Sirius University of Science and Technology, Sochi, 354340, Russia
| | - Tatiana V Tyazhelova
- Vavilov Institute of General Genetics, Russian Academy of Sciences, Moscow, 119333, Russia
| | - Fedor E Gusev
- Vavilov Institute of General Genetics, Russian Academy of Sciences, Moscow, 119333, Russia
| | - Andrey D Manakhov
- Vavilov Institute of General Genetics, Russian Academy of Sciences, Moscow, 119333, Russia
- Center for Genetics and Life Science, Sirius University of Science and Technology, Sochi, 354340, Russia
| | - Evgeny I Rogaev
- Vavilov Institute of General Genetics, Russian Academy of Sciences, Moscow, 119333, Russia.
- Center for Genetics and Life Science, Sirius University of Science and Technology, Sochi, 354340, Russia
- Department of Psychiatry, UMass Chan Medical School, Shrewsbury, MA 01545, USA
| |
Collapse
|
44
|
Yang C, Sobkowiak B, Naidu V, Codreanu A, Ciobanu N, Gunasekera KS, Chitwood MH, Alexandru S, Bivol S, Russi M, Havumaki J, Cudahy P, Fosburgh H, Allender CJ, Centner H, Engelthaler DM, Menzies NA, Warren JL, Crudu V, Colijn C, Cohen T. Phylogeography and transmission of M. tuberculosis in Moldova: A prospective genomic analysis. PLoS Med 2022; 19:e1003933. [PMID: 35192619 PMCID: PMC8903246 DOI: 10.1371/journal.pmed.1003933] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/17/2021] [Revised: 03/08/2022] [Accepted: 01/31/2022] [Indexed: 11/19/2022] Open
Abstract
BACKGROUND The incidence of multidrug-resistant tuberculosis (MDR-TB) remains critically high in countries of the former Soviet Union, where >20% of new cases and >50% of previously treated cases have resistance to rifampin and isoniazid. Transmission of resistant strains, as opposed to resistance selected through inadequate treatment of drug-susceptible tuberculosis (TB), is the main driver of incident MDR-TB in these countries. METHODS AND FINDINGS We conducted a prospective, genomic analysis of all culture-positive TB cases diagnosed in 2018 and 2019 in the Republic of Moldova. We used phylogenetic methods to identify putative transmission clusters; spatial and demographic data were analyzed to further describe local transmission of Mycobacterium tuberculosis. Of 2,236 participants, 779 (36%) had MDR-TB, of whom 386 (50%) had never been treated previously for TB. Moreover, 92% of multidrug-resistant M. tuberculosis strains belonged to putative transmission clusters. Phylogenetic reconstruction identified 3 large clades that were comprised nearly uniformly of MDR-TB: 2 of these clades were of Beijing lineage, and 1 of Ural lineage, and each had additional distinct clade-specific second-line drug resistance mutations and geographic distributions. Spatial and temporal proximity between pairs of cases within a cluster was associated with greater genomic similarity. Our study lasted for only 2 years, a relatively short duration compared with the natural history of TB, and, thus, the ability to infer the full extent of transmission is limited. CONCLUSIONS The MDR-TB epidemic in Moldova is associated with the local transmission of multiple M. tuberculosis strains, including distinct clades of highly drug-resistant M. tuberculosis with varying geographic distributions and drug resistance profiles. This study demonstrates the role of comprehensive genomic surveillance for understanding the transmission of M. tuberculosis and highlights the urgency of interventions to interrupt transmission of highly drug-resistant M. tuberculosis.
Collapse
Affiliation(s)
- Chongguang Yang
- School of Public Health (Shenzhen), Shenzhen Campus of Sun Yat-sen University, Shenzhen, China
- Department of Epidemiology of Microbial Diseases, Yale School of Public Health, New Haven, Connecticut, United States of America
| | | | - Vijay Naidu
- Department of Mathematics, Simon Fraser University, Burnaby, Canada
| | | | - Nelly Ciobanu
- Phthisiopneumology Institute, Chisinau, Republic of Moldova
| | - Kenneth S. Gunasekera
- Department of Epidemiology of Microbial Diseases, Yale School of Public Health, New Haven, Connecticut, United States of America
| | - Melanie H. Chitwood
- Department of Epidemiology of Microbial Diseases, Yale School of Public Health, New Haven, Connecticut, United States of America
| | | | - Stela Bivol
- Center for Health Policies and Studies, Chisinau, Republic of Moldova
| | - Marcus Russi
- Department of Epidemiology of Microbial Diseases, Yale School of Public Health, New Haven, Connecticut, United States of America
| | - Joshua Havumaki
- Department of Epidemiology of Microbial Diseases, Yale School of Public Health, New Haven, Connecticut, United States of America
| | - Patrick Cudahy
- Department of Epidemiology of Microbial Diseases, Yale School of Public Health, New Haven, Connecticut, United States of America
| | - Heather Fosburgh
- Department of Epidemiology of Microbial Diseases, Yale School of Public Health, New Haven, Connecticut, United States of America
| | | | - Heather Centner
- Translational Genomics Research Institute, Flagstaff, Arizona, United States of America
| | - David M. Engelthaler
- Translational Genomics Research Institute, Flagstaff, Arizona, United States of America
| | - Nicolas A. Menzies
- Department of Global Health and Population, and Center for Health Decision Science, Harvard TH Chan School of Public Health, Boston, Massachusetts, United States of America
| | - Joshua L. Warren
- Department of Biostatistics, Yale School of Public Health, New Haven, Connecticut, United States of America
| | - Valeriu Crudu
- Phthisiopneumology Institute, Chisinau, Republic of Moldova
| | - Caroline Colijn
- Department of Mathematics, Simon Fraser University, Burnaby, Canada
| | - Ted Cohen
- Department of Epidemiology of Microbial Diseases, Yale School of Public Health, New Haven, Connecticut, United States of America
| |
Collapse
|
45
|
Guyeux C, Senelle G, Refrégier G, Bretelle-Establet F, Cambau E, Sola C. Connection between two historical tuberculosis outbreak sites in Japan, Honshu, by a new ancestral Mycobacterium tuberculosis L2 sublineage. Epidemiol Infect 2022; 150:1-25. [PMID: 35042579 PMCID: PMC8931808 DOI: 10.1017/s0950268822000048] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2021] [Revised: 12/24/2021] [Accepted: 01/03/2022] [Indexed: 11/07/2022] Open
Abstract
By gathering 680 publicly available Sequence Read Archives from isolates of Mycobacterium tuberculosis complex (MTBC) including 190 belonging to the lineage 2 Beijing , and using an in-house bioinformatical pipeline, the TB-Annotator , that analyses more than 50 000 characters, we describe herein a new L2 sublineage from 20 isolates found in the Tochigi province, (Japan), that we designate as asia ancestral 5 (AAnc5). These isolates harbour a number of specific criteria (42 SNPs) and their intra-cluster pairwise distance suggests historical and not epidemiological transmission. These isolates harbour a mutation in rpoC , and do not fulfil, any of the modern Beijing lineage criteria, nor any of the other ancestral Beijing lineages described so far. Asia ancestral 5 isolates do not possess mutT2 58 and ogt 12 characteristics of modern Beijing , but possess ancestral Beijing SNPs characteristics. By looking into the literature, we found a reference isolate ID381, described in Kobe and Osaka belonging to the ‘G3’ group, sharing 36 out of the 42 specific SNPs found in AAnc5. We also assessed the intermediate position of the asia ancestral 4 (AAnc4) sublineage recently described in Thailand and propose an improved classification of the L2 that now includes AAnc4 and AAnc5. By increasing the recruitment into TB-Annotator to around 3000 genomes (including 642 belonging to L2), we confirmed our results and discovered additional historical ancestral L2 branches that remain to be investigated in more detail. We also present, in addition, some anthropological and historical data from Chinese and Japan history of tuberculosis, as well as from Korea, that could support our results on L2 evolution. This study shows that the reconstruction of the early history of tuberculosis in Asia is likely to reveal complex patterns since its emergence.
Collapse
Affiliation(s)
- Christophe Guyeux
- DISC Computer Science Department, FEMTO-ST Institute, UMR 6174 CNRS, Univ. Bourgogne Franche-Comté (UBFC), 16 Route de Gray, 25000Besançon, France
| | - Gaetan Senelle
- DISC Computer Science Department, FEMTO-ST Institute, UMR 6174 CNRS, Univ. Bourgogne Franche-Comté (UBFC), 16 Route de Gray, 25000Besançon, France
| | - Guislaine Refrégier
- Université Paris-Saclay, Saint-Aubin, France
- Université Paris-Saclay, CNRS, AgroParisTech, UMR ESE, 91405, Orsay, France
| | | | - Emmanuelle Cambau
- Université de Paris, IAME, UMR1137, INSERM, Paris, France
- AP-HP, GHU Nord, service de mycobactériologie spécialisée et de référence, Laboratoire associé du Centre National de Référence des mycobactéries et résistance des mycobactéries aux antituberculeux (CNR-MyRMA), Paris, France
| | - Christophe Sola
- Université Paris-Saclay, Saint-Aubin, France
- Université de Paris, IAME, UMR1137, INSERM, Paris, France
| |
Collapse
|
46
|
Torres Ortiz A, Coronel J, Vidal JR, Bonilla C, Moore DAJ, Gilman RH, Balloux F, Kon OM, Didelot X, Grandjean L. Genomic signatures of pre-resistance in Mycobacterium tuberculosis. Nat Commun 2021; 12:7312. [PMID: 34911948 PMCID: PMC8674244 DOI: 10.1038/s41467-021-27616-7] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2021] [Accepted: 11/29/2021] [Indexed: 11/29/2022] Open
Abstract
Recent advances in bacterial whole-genome sequencing have resulted in a comprehensive catalog of antibiotic resistance genomic signatures in Mycobacterium tuberculosis. With a view to pre-empt the emergence of resistance, we hypothesized that pre-existing polymorphisms in susceptible genotypes (pre-resistance mutations) could increase the risk of becoming resistant in the future. We sequenced whole genomes from 3135 isolates sampled over a 17-year period. After reconstructing ancestral genomes on time-calibrated phylogenetic trees, we developed and applied a genome-wide survival analysis to determine the hazard of resistance acquisition. We demonstrate that M. tuberculosis lineage 2 has a higher risk of acquiring resistance than lineage 4, and estimate a higher hazard of rifampicin resistance evolution following isoniazid mono-resistance. Furthermore, we describe loci and genomic polymorphisms associated with a higher risk of resistance acquisition. Identifying markers of future antibiotic resistance could enable targeted therapy to prevent resistance emergence in M. tuberculosis and other pathogens.
Collapse
Affiliation(s)
- Arturo Torres Ortiz
- grid.7445.20000 0001 2113 8111Imperial College London, Department of Infectious Diseases, London, UK
| | - Jorge Coronel
- grid.11100.310000 0001 0673 9488Universidad Peruana Cayetano Heredia, Lima, Perú
| | - Julia Rios Vidal
- grid.419858.90000 0004 0371 3700Unidad Técnica de Tuberculosis MDR, Ministerio de Salud, Lima, Perú
| | - Cesar Bonilla
- grid.419858.90000 0004 0371 3700Unidad Técnica de Tuberculosis MDR, Ministerio de Salud, Lima, Perú ,grid.441740.20000 0004 0542 2122Universidad Privada San Juan Bautista, Lima, Perú
| | - David A. J. Moore
- grid.8991.90000 0004 0425 469XLondon School of Hygiene and Tropical Medicine, London, UK
| | - Robert H. Gilman
- grid.21107.350000 0001 2171 9311Johns Hopkins Bloomberg School of Public Health, Baltimore, MD USA
| | | | - Onn Min Kon
- grid.7445.20000 0001 2113 8111Respiratory Medicine, National Heart and Lung Institute, Imperial College London, London, UK
| | - Xavier Didelot
- grid.7372.10000 0000 8809 1613University of Warwick, School of Life Sciences and Department of Statistics, Warwick, UK
| | - Louis Grandjean
- Imperial College London, Department of Infectious Diseases, London, UK. .,UCL Department of Infection, Institute of Child Health, London, UK.
| |
Collapse
|
47
|
van Tonder AJ, Thornton MJ, Conlan AJK, Jolley KA, Goolding L, Mitchell AP, Dale J, Palkopoulou E, Hogarth PJ, Hewinson RG, Wood JLN, Parkhill J. Inferring Mycobacterium bovis transmission between cattle and badgers using isolates from the Randomised Badger Culling Trial. PLoS Pathog 2021; 17:e1010075. [PMID: 34843579 PMCID: PMC8659364 DOI: 10.1371/journal.ppat.1010075] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2021] [Revised: 12/09/2021] [Accepted: 10/29/2021] [Indexed: 11/18/2022] Open
Abstract
Mycobacterium bovis (M. bovis) is a causative agent of bovine tuberculosis, a significant source of morbidity and mortality in the global cattle industry. The Randomised Badger Culling Trial was a field experiment carried out between 1998 and 2005 in the South West of England. As part of this trial, M. bovis isolates were collected from contemporaneous and overlapping populations of badgers and cattle within ten defined trial areas. We combined whole genome sequences from 1,442 isolates with location and cattle movement data, identifying transmission clusters and inferred rates and routes of transmission of M. bovis. Most trial areas contained a single transmission cluster that had been established shortly before sampling, often contemporaneous with the expansion of bovine tuberculosis in the 1980s. The estimated rate of transmission from badger to cattle was approximately two times higher than from cattle to badger, and the rate of within-species transmission considerably exceeded these for both species. We identified long distance transmission events linked to cattle movement, recurrence of herd breakdown by infection within the same transmission clusters and superspreader events driven by cattle but not badgers. Overall, our data suggests that the transmission clusters in different parts of South West England that are still evident today were established by long-distance seeding events involving cattle movement, not by recrudescence from a long-established wildlife reservoir. Clusters are maintained primarily by within-species transmission, with less frequent spill-over both from badger to cattle and cattle to badger.
Collapse
Affiliation(s)
- Andries J. van Tonder
- Department of Veterinary Medicine, University of Cambridge, Cambridge, United Kingdom
| | - Mark J. Thornton
- Department of Veterinary Medicine, University of Cambridge, Cambridge, United Kingdom
| | - Andrew J. K. Conlan
- Department of Veterinary Medicine, University of Cambridge, Cambridge, United Kingdom
| | - Keith A. Jolley
- Department of Zoology, University of Oxford, Oxford, United Kingdom
| | - Lee Goolding
- Animal and Plant Health Agency, New Haw, United Kingdom
| | | | - James Dale
- Animal and Plant Health Agency, New Haw, United Kingdom
| | | | | | | | - James L. N. Wood
- Department of Veterinary Medicine, University of Cambridge, Cambridge, United Kingdom
| | - Julian Parkhill
- Department of Veterinary Medicine, University of Cambridge, Cambridge, United Kingdom
| |
Collapse
|
48
|
Ruis C, Bryant JM, Bell SC, Thomson R, Davidson RM, Hasan NA, van Ingen J, Strong M, Floto RA, Parkhill J. Dissemination of Mycobacterium abscessus via global transmission networks. Nat Microbiol 2021; 6:1279-1288. [PMID: 34545208 PMCID: PMC8478660 DOI: 10.1038/s41564-021-00963-3] [Citation(s) in RCA: 43] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2021] [Accepted: 08/18/2021] [Indexed: 02/07/2023]
Abstract
Mycobacterium abscessus, a multidrug-resistant nontuberculous mycobacterium, has emerged as a major pathogen affecting people with cystic fibrosis (CF). Although originally thought to be acquired independently from the environment, most individuals are infected with one of several dominant circulating clones (DCCs), indicating the presence of global transmission networks of M. abscessus. How and when these clones emerged and spread globally is unclear. Here, we use evolutionary analyses of isolates from individuals both with and without CF to reconstruct the population history, spatiotemporal spread and recent transmission networks of the DCCs. We demonstrate synchronous expansion of six unrelated DCCs in the 1960s, a period associated with major changes in CF care and survival. Each of these clones has spread globally as a result of rare intercontinental transmission events. We show that the DCCs, but not environmentally acquired isolates, exhibit a specific smoking-associated mutational signature and that current transmission networks include individuals both with and without CF. We therefore propose that the DCCs initially emerged in non-CF populations but were then amplified and spread through the CF community. While individuals with CF are probably the most permissive host, non-CF individuals continue to play a key role in transmission networks and may facilitate long-distance transmission.
Collapse
Affiliation(s)
- Christopher Ruis
- Molecular Immunity Unit, University of Cambridge Department of Medicine, MRC-Laboratory of Molecular Biology, Cambridge, UK
- Department of Veterinary Medicine, University of Cambridge, Cambridge, UK
| | - Josephine M Bryant
- Molecular Immunity Unit, University of Cambridge Department of Medicine, MRC-Laboratory of Molecular Biology, Cambridge, UK
| | - Scott C Bell
- Children's Health Research Institute, The University of Queensland, Brisbane, Australia
- Thoracic Medicine, The Prince Charles Hospital, Brisbane, Australia
| | - Rachel Thomson
- Gallipoli Medical Research Institute, The University of Queensland, Brisbane, Australia
| | - Rebecca M Davidson
- Center for Genes, Environment and Health, National Jewish Health, Denver, CO, USA
| | - Nabeeh A Hasan
- Center for Genes, Environment and Health, National Jewish Health, Denver, CO, USA
| | - Jakko van Ingen
- Center for Infectious Diseases, Department of Medical Microbiology, Radboud University Medical Center, Nijmegen, Netherlands
| | - Michael Strong
- Center for Genes, Environment and Health, National Jewish Health, Denver, CO, USA
| | - R Andres Floto
- Molecular Immunity Unit, University of Cambridge Department of Medicine, MRC-Laboratory of Molecular Biology, Cambridge, UK.
- Cambridge Centre for Lung Infection, Papworth Hospital, Cambridge, UK.
- Cambridge Centre for AI in Medicine, University of Cambridge, Cambridge, UK.
| | - Julian Parkhill
- Department of Veterinary Medicine, University of Cambridge, Cambridge, UK.
| |
Collapse
|
49
|
Boom WH, Schaible UE, Achkar JM. The knowns and unknowns of latent Mycobacterium tuberculosis infection. J Clin Invest 2021; 131:136222. [PMID: 33529162 DOI: 10.1172/jci136222] [Citation(s) in RCA: 56] [Impact Index Per Article: 18.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Humans have been infected with Mycobacterium tuberculosis (Mtb) for thousands of years. While tuberculosis (TB), one of the deadliest infectious diseases, is caused by uncontrolled Mtb infection, over 90% of presumed infected individuals remain asymptomatic and contain Mtb in a latent TB infection (LTBI) without ever developing disease, and some may clear the infection. A small number of heavily Mtb-exposed individuals appear to resist developing traditional LTBI. Because Mtb has mechanisms for intracellular survival and immune evasion, successful control involves all of the arms of the immune system. Here, we focus on immune responses to Mtb in humans and nonhuman primates and discuss new concepts and outline major knowledge gaps in our understanding of LTBI, ranging from the earliest events of exposure and infection to success or failure of Mtb control.
Collapse
Affiliation(s)
- W Henry Boom
- Department of Medicine.,Department of Pathology, and.,Department of Molecular Biology and Microbiology, Case Western Reserve University and University Hospitals Cleveland Medical Center, Cleveland, Ohio, USA
| | - Ulrich E Schaible
- Division of Cellular Microbiology, Research Center Borstel-Leibniz Lung Center, Borstel, Germany.,German Center for Infection Research, partner site Hamburg-Lübeck-Borstel-Riems, Germany
| | - Jacqueline M Achkar
- Department of Medicine and.,Department of Microbiology and Immunology, Albert Einstein College of Medicine, Bronx, New York, USA
| |
Collapse
|
50
|
Calland JK, Pascoe B, Bayliss SC, Mourkas E, Berthenet E, Thorpe HA, Hitchings MD, Feil EJ, Corander J, Blaser MJ, Falush D, Sheppard SK. Quantifying bacterial evolution in the wild: A birthday problem for Campylobacter lineages. PLoS Genet 2021; 17:e1009829. [PMID: 34582435 PMCID: PMC8500405 DOI: 10.1371/journal.pgen.1009829] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2021] [Revised: 10/08/2021] [Accepted: 09/20/2021] [Indexed: 11/20/2022] Open
Abstract
Measuring molecular evolution in bacteria typically requires estimation of the rate at which nucleotide changes accumulate in strains sampled at different times that share a common ancestor. This approach has been useful for dating ecological and evolutionary events that coincide with the emergence of important lineages, such as outbreak strains and obligate human pathogens. However, in multi-host (niche) transmission scenarios, where the pathogen is essentially an opportunistic environmental organism, sampling is often sporadic and rarely reflects the overall population, particularly when concentrated on clinical isolates. This means that approaches that assume recent common ancestry are not applicable. Here we present a new approach to estimate the molecular clock rate in Campylobacter that draws on the popular probability conundrum known as the 'birthday problem'. Using large genomic datasets and comparative genomic approaches, we use isolate pairs that share recent common ancestry to estimate the rate of nucleotide change for the population. Identifying synonymous and non-synonymous nucleotide changes, both within and outside of recombined regions of the genome, we quantify clock-like diversification to estimate synonymous rates of nucleotide change for the common pathogenic bacteria Campylobacter coli (2.4 x 10-6 s/s/y) and Campylobacter jejuni (3.4 x 10-6 s/s/y). Finally, using estimated total rates of nucleotide change, we infer the number of effective lineages within the sample time frame-analogous to a shared birthday-and assess the rate of turnover of lineages in our sample set over short evolutionary timescales. This provides a generalizable approach to calibrating rates in populations of environmental bacteria and shows that multiple lineages are maintained, implying that large-scale clonal sweeps may take hundreds of years or more in these species.
Collapse
Affiliation(s)
- Jessica K. Calland
- The Milner Centre for Evolution, University of Bath, Bath, United Kingdom
| | - Ben Pascoe
- The Milner Centre for Evolution, University of Bath, Bath, United Kingdom
| | - Sion C. Bayliss
- The Milner Centre for Evolution, University of Bath, Bath, United Kingdom
| | - Evangelos Mourkas
- The Milner Centre for Evolution, University of Bath, Bath, United Kingdom
| | - Elvire Berthenet
- French National Reference Center for Campylobacters and Helicobacters, University of Bordeaux, Bordeaux, France
- Institute of Life Sciences, Swansea University Medical School, Swansea University, Singleton Park, Swansea, United Kingdom
| | - Harry A. Thorpe
- The Milner Centre for Evolution, University of Bath, Bath, United Kingdom
- Department of Biostatistics, University of Oslo, Oslo, Norway
| | - Matthew D. Hitchings
- Institute of Life Sciences, Swansea University Medical School, Swansea University, Singleton Park, Swansea, United Kingdom
| | - Edward J. Feil
- The Milner Centre for Evolution, University of Bath, Bath, United Kingdom
| | - Jukka Corander
- Department of Biostatistics, University of Oslo, Oslo, Norway
- Department of Mathematics and Statistics, Helsinki Institute for Information Technology, University of Helsinki, Helsinki, Finland
- Parasites and Microbes, Wellcome Sanger Institute, Cambridge, United Kingdom
| | - Martin J. Blaser
- Center for Advanced Biotechnology and Medicine, Rutgers University, New Brunswick, New Jersey, United States of America
| | - Daniel Falush
- Centre for Microbes, Development and Health, Institute Pasteur of Shanghai, Shanghai, China
- * E-mail: (DF); (SKS)
| | - Samuel K. Sheppard
- The Milner Centre for Evolution, University of Bath, Bath, United Kingdom
- Department of Zoology, University of Oxford, Oxford, United Kingdom
- * E-mail: (DF); (SKS)
| |
Collapse
|