Armed conflict and population displacement as drivers of the evolution and dispersal of Mycobacterium tuberculosis

Significance We used population genomic analyses to reconstruct the recent history and dispersal of a major clade of Mycobacterium tuberculosis in central Asia and beyond. Our results indicate that the fall of the Soviet Union and the ensuing collapse of public health systems led to a rise in M. tuberculosis drug resistance. We also show that armed conflict and population displacement is likely to have aided the export of this clade from central Asia to war-torn Afghanistan and beyond. The “Beijing” Mycobacterium tuberculosis (Mtb) lineage 2 (L2) is spreading globally and has been associated with accelerated disease progression and increased antibiotic resistance. Here we performed a phylodynamic reconstruction of one of the L2 sublineages, the central Asian clade (CAC), which has recently spread to western Europe. We find that recent historical events have contributed to the evolution and dispersal of the CAC. Our timing estimates indicate that the clade was likely introduced to Afghanistan during the 1979–1989 Soviet–Afghan war and spread further after population displacement in the wake of the American invasion in 2001. We also find that drug resistance mutations accumulated on a massive scale in Mtb isolates from former Soviet republics after the fall of the Soviet Union, a pattern that was not observed in CAC isolates from Afghanistan. Our results underscore the detrimental effects of political instability and population displacement on tuberculosis control and demonstrate the power of phylodynamic methods in exploring bacterial evolution in space and time.

[1]  Y. Balabanova,et al.  Drug-resistant tuberculosis, clinical virulence, and the dominance of the Beijing strain family in Russia. , 2005, JAMA.

[2]  Xavier Didelot,et al.  Bayesian Inference of Infectious Disease Transmission from Whole-Genome Sequence Data , 2014, Molecular biology and evolution.

[3]  Gang Sun,et al.  Southern East Asian origin and coexpansion of Mycobacterium tuberculosis Beijing family with Han Chinese , 2015, Proceedings of the National Academy of Sciences.

[4]  F. Balloux,et al.  Evolution of extensively drug-resistant Mycobacterium tuberculosis from a susceptible ancestor in a single patient , 2014, Genome Biology.

[5]  Marc A Suchard,et al.  Using Time-Structured Data to Estimate Evolutionary Rates of Double-Stranded DNA Viruses , 2010, Molecular biology and evolution.

[6]  Gemma L. Kay,et al.  Eighteenth-century genomes show that mixed infections were common at time of peak tuberculosis in Europe , 2015, Nature Communications.

[7]  Stefan Niemann,et al.  Whole-genome sequencing of rifampicin-resistant Mycobacterium tuberculosis strains identifies compensatory mutations in RNA polymerase genes , 2011 .

[8]  Peer Bork,et al.  Interactive Tree Of Life v2: online annotation and display of phylogenetic trees made easy , 2011, Nucleic Acids Res..

[9]  F. Balloux,et al.  Antimicrobial Resistance in Mycobacterium tuberculosis: The Odd One Out. , 2016, Trends in microbiology.

[10]  M. Chase,et al.  The mutation rate of mycobacterial repetitive unit loci in strains of M. tuberculosis from cynomolgus macaque infection , 2013, BMC Genomics.

[11]  P. Small,et al.  Progression to active tuberculosis, but not transmission, varies by Mycobacterium tuberculosis lineage in The Gambia. , 2008, The Journal of infectious diseases.

[12]  Omar E. Cornejo,et al.  The Role of Selection in Shaping Diversity of Natural M. tuberculosis Populations , 2013, PLoS pathogens.

[13]  D. Minnikin,et al.  Human tuberculosis predates domestication in ancient Syria. , 2015, Tuberculosis.

[14]  Jukka Corander,et al.  Evolution and transmission of drug resistant tuberculosis in a Russian population , 2014, Nature Genetics.

[15]  Larry D. Martin,et al.  Mycobacterium tuberculosis Complex Lipid Virulence Factors Preserved in the 17,000-Year-Old Skeleton of an Extinct Bison, Bison antiquus , 2012, PloS one.

[16]  Genome-scale rates of evolutionary change in bacteria , 2016 .

[17]  Michael P. Cummings,et al.  BEAST (Bayesian Evolutionary Analysis by Sampling Trees) , 2014 .

[18]  Daniel J. Wilson,et al.  Whole-genome sequencing to delineate Mycobacterium tuberculosis outbreaks: a retrospective observational study , 2013, The Lancet. Infectious diseases.

[19]  Scott A. Sisson,et al.  A Model-Based Bayesian Estimation of the Rate of Evolution of VNTR Loci in Mycobacterium tuberculosis , 2012, PLoS Comput. Biol..

[20]  Ramón Doallo,et al.  CircadiOmics: integrating circadian genomics, transcriptomics, proteomics and metabolomics , 2012, Nature Methods.

[21]  M. Chase,et al.  Use of whole genome sequencing to estimate the mutation rate of Mycobacterium tuberculosis during latent infection , 2011, Nature Genetics.

[22]  A. Ruiz-Linares,et al.  Improved Calibration of the Human Mitochondrial Clock Using Ancient Genomes , 2014, Molecular biology and evolution.

[23]  M. Nei,et al.  MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. , 2011, Molecular biology and evolution.

[24]  L. Jelsbak,et al.  Substantial molecular evolution and mutation rates in prolonged latent Mycobacterium tuberculosis infection in humans. , 2016, International journal of medical microbiology : IJMM.

[25]  Nalin Rastogi,et al.  Evolutionary history and global spread of the Mycobacterium tuberculosis Beijing lineage , 2015, Nature Genetics.

[26]  David J. Edwards,et al.  Genome-scale rates of evolutionary change in bacteria , 2016, bioRxiv.

[27]  M. Suchard,et al.  Bayesian Phylogenetics with BEAUti and the BEAST 1.7 , 2012, Molecular biology and evolution.

[28]  Kelly M. Harkins,et al.  Pre-Columbian mycobacterial genomes reveal seals as a source of New World human tuberculosis , 2014, Nature.

[29]  A. Rambaut,et al.  BEAST: Bayesian evolutionary analysis by sampling trees , 2007, BMC Evolutionary Biology.

[30]  O. Pybus,et al.  Measurably evolving pathogens in the genomic era. , 2015, Trends in ecology & evolution.

[31]  S. Hoffner,et al.  Drug-Susceptible Mycobacterium tuberculosis Beijing Genotype Does Not Develop Mutation-Conferred Resistance to Rifampin at an Elevated Rate , 2003, Journal of Clinical Microbiology.

[32]  M. Borgdorff,et al.  Beijing/W Genotype Mycobacterium tuberculosis and Drug Resistance , 2006, Emerging infectious diseases.

[33]  Marc Lipsitch,et al.  Mycobacterium tuberculosis mutation rate estimates from different lineages predict substantial differences in the emergence of drug resistant tuberculosis , 2013, Nature Genetics.

[34]  F. Balloux,et al.  Four decades of transmission of a multidrug-resistant Mycobacterium tuberculosis outbreak strain , 2015, Nature Communications.

[35]  M. Suchard,et al.  Smooth skyride through a rough skyline: Bayesian coalescent-based inference of population dynamics. , 2008, Molecular biology and evolution.

[36]  C. Garzelli Beijing/W Genotype Mycobacterium tuberculosis and Drug Resistance , 2006 .

[37]  Stefan Niemann,et al.  Whole Genome Sequencing versus Traditional Genotyping for Investigation of a Mycobacterium tuberculosis Outbreak: A Longitudinal Molecular Epidemiological Study , 2013, PLoS medicine.

[38]  A. Vogler,et al.  Effect of Repeat Copy Number on Variable-Number Tandem Repeat Mutations in Escherichia coli O157:H7 , 2006, Journal of bacteriology.

[39]  K. Holt,et al.  Out-of-Africa migration and Neolithic co-expansion of Mycobacterium tuberculosis with modern humans , 2013, Nature Genetics.

[40]  Falk Hildebrand,et al.  Origin, Spread and Demography of the Mycobacterium tuberculosis Complex , 2008, PLoS pathogens.

[41]  Søren Brunak,et al.  Early Divergent Strains of Yersinia pestis in Eurasia 5,000 Years Ago , 2015, Cell.

[42]  O. Gascuel,et al.  SeaView version 4: A multiplatform graphical user interface for sequence alignment and phylogenetic tree building. , 2010, Molecular biology and evolution.

[43]  A. Rodrigo,et al.  Measurably evolving populations , 2003 .

[44]  Alexei J Drummond,et al.  Time dependency of molecular rate estimates and systematic overestimation of recent divergence times. , 2005, Molecular biology and evolution.