1
|
Filardi ETM, Pucca MB, Araujo Junior JP, da Costa PI. Pandemic paradox: the impact of the COVID-19 on the global and Brazilian tuberculosis epidemics. Front Public Health 2024; 12:1399860. [PMID: 39131577 PMCID: PMC11312383 DOI: 10.3389/fpubh.2024.1399860] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2024] [Accepted: 07/10/2024] [Indexed: 08/13/2024] Open
Affiliation(s)
- Eloise T. M. Filardi
- Department of Clinical Analysis, School of Pharmaceutical Sciences, São Paulo State University (UNESP), Araraquara, São Paulo, Brazil
| | - Manuela B. Pucca
- Department of Clinical Analysis, School of Pharmaceutical Sciences, São Paulo State University (UNESP), Araraquara, São Paulo, Brazil
| | - João Pessoa Araujo Junior
- Department of Chemical and Biological Sciences, Institute of Biosciences, São Paulo State University (UNESP), Botucatu, São Paulo, Brazil
| | - Paulo I. da Costa
- Department of Clinical Analysis, School of Pharmaceutical Sciences, São Paulo State University (UNESP), Araraquara, São Paulo, Brazil
| |
Collapse
|
2
|
Morey-León G, Mejía-Ponce PM, Granda Pardo JC, Muñoz-Mawyin K, Fernández-Cadena JC, García-Moreira E, Andrade-Molina D, Licona-Cassani C, Berná L. A precision overview of genomic resistance screening in Ecuadorian isolates of Mycobacterium tuberculosis using web-based bioinformatics tools. PLoS One 2023; 18:e0294670. [PMID: 38051742 DOI: 10.1371/journal.pone.0294670] [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: 03/18/2023] [Accepted: 11/07/2023] [Indexed: 12/07/2023] Open
Abstract
INTRODUCTION Tuberculosis (TB) is among the deadliest diseases worldwide, and its impact is mainly due to the continuous emergence of resistant isolates during treatment due to the laborious process of resistance diagnosis, nonadherence to treatment and circulation of previously resistant isolates of Mycobacterium tuberculosis. In this study, we evaluated the performance and functionalities of web-based tools, including Mykrobe, TB-profiler, PhyResSE, KvarQ, and SAM-TB, for detecting resistance in 88 Ecuadorian isolates of Mycobacterium tuberculosis drug susceptibility tested previously. Statistical analysis was used to determine the correlation between genomic and phenotypic analysis. Our results showed that with the exception of KvarQ, all tools had the highest correlation with the conventional drug susceptibility test (DST) for global resistance detection (98% agreement and 0.941 Cohen's kappa), while SAM-TB, PhyResSE, TB-profiler and Mykrobe had better correlations with DST for first-line drug analysis individually. We also identified that in our study, only 50% of mutations characterized by the web-based tools in the rpoB, katG, embB, pncA, gyrA and rrs regions were canonical and included in the World Health Organization (WHO) catalogue. Our findings suggest that SAM-TB, PhyResSE, TB-profiler and Mykrobe were efficient in determining canonical resistance-related mutations, but more analysis is needed to improve second-line detection. Improving surveillance programs using whole-genome sequencing tools for first-line drugs, MDR-TB and XDR-TB is essential to understand the molecular epidemiology of TB in Ecuador. IMPORTANCE Tuberculosis, an infectious disease caused by Mycobacterium tuberculosis, most commonly affects the lungs and is often spread through the air when infected people cough, sneeze, or spit. However, despite the existence of effective drug treatment, patient adherence, long duration of treatment, and late diagnosis have reduced the effectiveness of therapy and increased drug resistance. The increase in resistant cases, added to the impact of the COVID-19 pandemic, has highlighted the importance of implementing efficient and timely diagnostic methodologies worldwide. The significance of our research is in evaluating and identifying a more efficient and user-friendly web-based tool to characterize resistance in Mycobacterium tuberculosis by whole-genome sequencing, which will allow more routine application to improve TB strain surveillance programs locally.
Collapse
Affiliation(s)
- Gabriel Morey-León
- Facultad de Ciencias de la Salud, Universidad Espíritu Santo, Samborondón, Ecuador
- Universidad de la República, Montevideo, Uruguay
- University of Guayaquil, Guayaquil, Ecuador
| | - Paulina M Mejía-Ponce
- Escuela de Ingeniería y Ciencias, Tecnológico de Monterrey, Monterrey, Nuevo León, México
| | - Juan Carlos Granda Pardo
- Centro de Referencia Nacional de Micobacterias, Instituto Nacional de Investigación en Salud Pública Dr Leopoldo Izquieta Perez, INSPI-LIP, Guayaquil, Ecuador
| | - Karen Muñoz-Mawyin
- Laboratorio de Ciencias Ómicas, Universidad Espíritu Santo, Samborondón, Ecuador
| | | | | | - Derly Andrade-Molina
- Facultad de Ciencias de la Salud, Universidad Espíritu Santo, Samborondón, Ecuador
- Laboratorio de Ciencias Ómicas, Universidad Espíritu Santo, Samborondón, Ecuador
| | | | - Luisa Berná
- Laboratorio de Interacciones Hospedero-Patógeno, Unidad de Biología Molecular, Institut Pasteur de Montevideo, Montevideo, Uruguay
- Unidad de Genómica Evolutiva, Facultad de Ciencias, Universidad de la República, Montevideo, Uruguay
| |
Collapse
|
3
|
Trisakul K, Nonghanphithak D, Chaiyachat P, Kaewprasert O, Sakmongkoljit K, Reechaipichitkul W, Chaiprasert A, Blair D, Clark TG, Faksri K. High clustering rate and genotypic drug-susceptibility screening for the newly recommended anti-tuberculosis drugs among global extensively drug-resistant Mycobacterium tuberculosis isolates. EMERGING MICROBES & INFECTIONS 2022; 11:1857-1866. [PMID: 35792049 PMCID: PMC9336503 DOI: 10.1080/22221751.2022.2099304] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
Multidrug-resistant tuberculosis (MDR-TB) and extensively drug-resistant tuberculosis (XDR-TB) make TB difficult to control. Global susceptibility data for six newly recommended anti-TB drugs against M/XDR-TB are still limited. Using publicly available whole-genome sequences, we determined the proportion of 513 phenotypically XDR-TB isolates that carried mutations associated with resistance against these drugs (bedaquiline, clofazimine, linezolid, delamanid, pretomanid and cycloserine). Mutations of Rv0678 and Rv1979c were detected in 69/513 isolates (13.5%) for bedaquiline resistance and 79/513 isolates (15.4%) for clofazimine resistance with additional mmpL5 mutations. Mutations conferring resistance to delamanid were detected in fbiB and ddn genes for 11/513 isolates (2.1%). For pretomanid, a mutation was detected in the ddn gene for 3/513 isolates (0.6%). Nineteen mutations of pykA, cycA, ald, and alr genes, conferring resistance to cycloserine, were found in 153/513 isolates (29.8%). No known mutations associated with linezolid resistance were detected. Cluster analysis showed that 408/513 isolates fell within 99 clusters and that 354 of these isolates were possible primary drug-resistant TB (292 XDR-TB, 57 pre-XDR-TB and 5 MDR-TB). Clonal transmission of primary XDR isolates might contribute significantly to the high prevalence of DR-TB globally.
Collapse
Affiliation(s)
- Kanwara Trisakul
- Department of Microbiology, Faculty of Medicine, Khon Kaen University, Khon Kaen, Thailand
- Research and Diagnostic Center for Emerging Infectious Diseases (RCEID), Khon Kaen University, Khon Kaen, Thailand
| | - Ditthawat Nonghanphithak
- Department of Microbiology, Faculty of Medicine, Khon Kaen University, Khon Kaen, Thailand
- Research and Diagnostic Center for Emerging Infectious Diseases (RCEID), Khon Kaen University, Khon Kaen, Thailand
| | - Pratchakan Chaiyachat
- Department of Microbiology, Faculty of Medicine, Khon Kaen University, Khon Kaen, Thailand
- Research and Diagnostic Center for Emerging Infectious Diseases (RCEID), Khon Kaen University, Khon Kaen, Thailand
| | - Orawee Kaewprasert
- Department of Microbiology, Faculty of Medicine, Khon Kaen University, Khon Kaen, Thailand
- Research and Diagnostic Center for Emerging Infectious Diseases (RCEID), Khon Kaen University, Khon Kaen, Thailand
| | - Kankanon Sakmongkoljit
- Department of Geotechnology, Faculty of Technology, Khon Kaen University, Khon Kaen, Thailand
| | - Wipa Reechaipichitkul
- Department of Microbiology, Faculty of Medicine, Khon Kaen University, Khon Kaen, Thailand
- Research and Diagnostic Center for Emerging Infectious Diseases (RCEID), Khon Kaen University, Khon Kaen, Thailand
| | - Angkana Chaiprasert
- Office for Research and Development, Faculty of Medicine Siriraj Hospital, Mahidol University, Bangkok, Thailand
| | - David Blair
- James Cook University, Townsville, Queensland, Australia
| | - Taane G. Clark
- London School of Hygiene and Tropical Medicine, London, United Kingdom
| | - Kiatichai Faksri
- Department of Microbiology, Faculty of Medicine, Khon Kaen University, Khon Kaen, Thailand
- Research and Diagnostic Center for Emerging Infectious Diseases (RCEID), Khon Kaen University, Khon Kaen, Thailand
| |
Collapse
|
4
|
Conceição EC, Salvato RS, Gomes KM, Guimarães AEDS, da Conceição ML, Souza e Guimarães RJDP, Sharma A, Furlaneto IP, Barcellos RB, Bollela VR, Anselmo LMP, Sisco MC, Niero CV, Ferrazoli L, Refrégier G, Lourenço MCDS, Gomes HM, de Brito AC, Catanho M, Duarte RS, Suffys PN, Lima KVB. Molecular epidemiology of Mycobacterium tuberculosis in Brazil before the whole genome sequencing era: a literature review. Mem Inst Oswaldo Cruz 2021; 116:e200517. [PMID: 33729319 PMCID: PMC7976556 DOI: 10.1590/0074-02760200517] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2020] [Accepted: 02/11/2021] [Indexed: 11/22/2022] Open
Abstract
Molecular-typing can help in unraveling epidemiological scenarios and improvement for disease control strategies. A literature review of Mycobacterium tuberculosis transmission in Brazil through genotyping on 56 studies published from 1996-2019 was performed. The clustering rate for mycobacterial interspersed repetitive units - variable tandem repeats (MIRU-VNTR) of 1,613 isolates were: 73%, 33% and 28% based on 12, 15 and 24-loci, respectively; while for RFLP-IS6110 were: 84% among prison population in Rio de Janeiro, 69% among multidrug-resistant isolates in Rio Grande do Sul, and 56.2% in general population in São Paulo. These findings could improve tuberculosis (TB) surveillance and set up a solid basis to build a database of Mycobacterium genomes.
Collapse
Affiliation(s)
- Emilyn Costa Conceição
- Fundação Oswaldo Cruz-Fiocruz, Instituto Nacional de Infectologia
Evandro Chagas, Programa de Pós-Graduação em Pesquisa Clínica e Doenças Infecciosas,
Rio de Janeiro, RJ, Brasil
- Fundação Oswaldo Cruz-Fiocruz, Instituto Nacional de Infectologia
Evandro Chagas, Laboratório de Bacteriologia e Bioensaios em Micobactérias, Rio de
Janeiro, RJ, Brasil
- Fundação Oswaldo Cruz-Fiocruz, Instituto Oswaldo Cruz, Laboratório
de Biologia Molecular Aplicada a Micobactérias, Rio de Janeiro, RJ, Brasil
| | - Richard Steiner Salvato
- Universidade Federal do Rio Grande do Sul, Programa de Pós-Graduação
em Biologia Celular e Molecular, Porto Alegre, RS, Brasil
- Secretaria Estadual de Saúde do Rio Grande do Sul, Centro Estadual
de Vigilância em Saúde, Centro de Desenvolvimento Científico e Tecnológico, Porto
Alegre, RS, Brasil
| | - Karen Machado Gomes
- Fundação Oswaldo Cruz-Fiocruz, Escola Nacional de Saúde Pública
Sergio Arouca, Centro de Referência Professor Hélio Fraga, Laboratório de Referência
Nacional para Tuberculose e outras Micobacterioses, Rio de Janeiro, RJ, Brasil
| | - Arthur Emil dos Santos Guimarães
- Universidade do Estado do Pará, Instituto de Ciências Biológicas e
da Saúde, Pós-Graduação Biologia Parasitária na Amazônia, Belém, PA, Brasil
- Instituto Evandro Chagas, Seção de Bacteriologia e Micologia,
Ananindeua, PA, Brasil
| | - Marília Lima da Conceição
- Universidade do Estado do Pará, Instituto de Ciências Biológicas e
da Saúde, Pós-Graduação Biologia Parasitária na Amazônia, Belém, PA, Brasil
- Instituto Evandro Chagas, Seção de Bacteriologia e Micologia,
Ananindeua, PA, Brasil
| | | | - Abhinav Sharma
- International Institute of Information Technology, Department of
Data Science, Bangalore, India
| | | | - Regina Bones Barcellos
- Secretaria Estadual de Saúde do Rio Grande do Sul, Centro Estadual
de Vigilância em Saúde, Centro de Desenvolvimento Científico e Tecnológico, Porto
Alegre, RS, Brasil
| | - Valdes Roberto Bollela
- Universidade de São Paulo, Departamento de Clínica Médica da
Faculdade de Medicina de Ribeirão Preto, Ribeirão Preto, SP, Brasil
| | - Lívia Maria Pala Anselmo
- Universidade de São Paulo, Departamento de Clínica Médica da
Faculdade de Medicina de Ribeirão Preto, Ribeirão Preto, SP, Brasil
| | - Maria Carolina Sisco
- Fundação Oswaldo Cruz-Fiocruz, Instituto Oswaldo Cruz, Laboratório
de Biologia Molecular Aplicada a Micobactérias, Rio de Janeiro, RJ, Brasil
- Universidade Federal do Rio de Janeiro, Instituto de Microbiologia
Paulo de Góes, Laboratório de Micobactérias, Rio de Janeiro, RJ, Brasil
| | - Cristina Viana Niero
- Universidade Federal de São Paulo, Departamento de Microbiologia,
Imunologia e Parasitologia, São Paulo, SP, Brasil
| | - Lucilaine Ferrazoli
- Instituto Adolfo Lutz, Centro de Bacteriologia, Núcleo de
Tuberculose e Micobacterioses, São Paulo, SP, Brasil
| | - Guislaine Refrégier
- Universit e Paris-Saclay, Ecologie Systematique Evolution, Centre
National de la Recherche Scientifique, AgroParisTech, Orsay, France
| | - Maria Cristina da Silva Lourenço
- Fundação Oswaldo Cruz-Fiocruz, Instituto Nacional de Infectologia
Evandro Chagas, Laboratório de Bacteriologia e Bioensaios em Micobactérias, Rio de
Janeiro, RJ, Brasil
| | - Harrison Magdinier Gomes
- Fundação Oswaldo Cruz-Fiocruz, Instituto Oswaldo Cruz, Laboratório
de Biologia Molecular Aplicada a Micobactérias, Rio de Janeiro, RJ, Brasil
| | - Artemir Coelho de Brito
- Coordenação Geral de Vigilância das Doenças de Transmissão
Respiratória de Condições Crônicas, Brasília, DF, Brasil
| | - Marcos Catanho
- Fundação Oswaldo Cruz-Fiocruz, Instituto Oswaldo Cruz, Laboratório
de Genética Molecular de Microrganismos, Rio de Janeiro, RJ, Brasil
| | - Rafael Silva Duarte
- Universidade Federal do Rio de Janeiro, Instituto de Microbiologia
Paulo de Góes, Laboratório de Micobactérias, Rio de Janeiro, RJ, Brasil
| | - Philip Noel Suffys
- Fundação Oswaldo Cruz-Fiocruz, Instituto Oswaldo Cruz, Laboratório
de Biologia Molecular Aplicada a Micobactérias, Rio de Janeiro, RJ, Brasil
| | - Karla Valéria Batista Lima
- Universidade do Estado do Pará, Instituto de Ciências Biológicas e
da Saúde, Pós-Graduação Biologia Parasitária na Amazônia, Belém, PA, Brasil
- Instituto Evandro Chagas, Seção de Bacteriologia e Micologia,
Ananindeua, PA, Brasil
| |
Collapse
|
5
|
Genomic Variations in Drug Resistant Mycobacterium tuberculosis Strains Collected from Patients with Different Localization of Infection. Antibiotics (Basel) 2020; 10:antibiotics10010027. [PMID: 33396320 PMCID: PMC7824472 DOI: 10.3390/antibiotics10010027] [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: 11/23/2020] [Revised: 12/23/2020] [Accepted: 12/25/2020] [Indexed: 11/28/2022] Open
Abstract
Mycobacterium tuberculosis is a highly studied pathogen due to public health importance. Despite this, problems like early drug resistance, diagnostics and treatment success prediction are still not fully resolved. Here, we analyze the incidence of point mutations widely used for drug resistance detection in laboratory practice and conduct comparative analysis of whole-genome sequence (WGS) for clinical M. tuberculosis strains collected from patients with pulmonary tuberculosis (PTB) and extra-pulmonary tuberculosis (XPTB) localization. A total of 72 pulmonary and 73 extrapulmonary microbiologically characterized M. tuberculosis isolates were collected from patients from 2007 to 2014 in Russia. Genomic DNA was used for WGS and obtained data allowed identifying major mutations known to be associated with drug resistance to first-line and second-line antituberculous drugs. In some cases previously described mutations were not identified. Using genome-based phylogenetic analysis we identified M. tuberculosis substrains associated with distinctions in the occurrence in PTB vs. XPTB cases. Phylogenetic analyses did reveal M. tuberculosis genetic substrains associated with TB localization. XPTB was associated with Beijing sublineages Central Asia (Beijing CAO), Central Asia Clade A (Beijing A) and 4.8 groups, while PTB localization was associated with group LAM (4.3). Further, the XPTB strain in some cases showed elevated drug resistance patterns relative to PTB isolates. HIV was significantly associated with the development of XPTB in the Beijing B0/W148 group and among unclustered Beijing isolates.
Collapse
|
6
|
Hu D, Liu B, Wang L, Reeves PR. Living Trees: High-Quality Reproducible and Reusable Construction of Bacterial Phylogenetic Trees. Mol Biol Evol 2020; 37:563-575. [PMID: 31633785 DOI: 10.1093/molbev/msz241] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
An ideal bacterial phylogenetic tree accurately retraces evolutionary history and accurately incorporates mutational, recombination and other events on the appropriate branches. Current strain-level bacterial phylogenetic analysis based on large numbers of genomes lacks reliability and resolution, and is hard to be replicated, confirmed and reused, because of the highly divergent nature of microbial genomes. We present SNPs and Recombination Events Tree (SaRTree), a pipeline using six "living trees" modules that addresses problems arising from the high numbers and variable quality of bacterial genome sequences. It provides for reuse of the tree and offers a major step toward global standardization of phylogenetic analysis by generating deposit files including all steps involved in phylogenetic inference. The tree itself is a "living tree" that can be extended by addition of more sequences, or the deposit can be used to vary the programs or parameters used, to assess the effect of such changes. This approach will allow phylogeny papers to meet the traditional responsibility of providing data and analysis that can be repeated and critically evaluated by others. We used the Acinetobacter baumannii global clone I to illustrate use of SaRTree to optimize tree resolution. An Escherichia coli tree was built from 351 sequences selected from 11,162 genome sequences, with the others added back onto well-defined branches, to show how this facility can greatly improve the outcomes from genome sequencing. SaRTree is designed for prokaryote strain-level analysis but could be adapted for other usage.
Collapse
Affiliation(s)
- Dalong Hu
- School of Life and Environmental Sciences, University of Sydney, Sydney, NSW, Australia
| | - Bin Liu
- TEDA Institute of Biological Sciences and Biotechnology, Nankai University, Tianjin Economic-Technological Development Area, Tianjin, People's Republic of China.,Tianjin Research Center for Functional Genomics and Biochip, Tianjin, People's Republic of China
| | - Lei Wang
- TEDA Institute of Biological Sciences and Biotechnology, Nankai University, Tianjin Economic-Technological Development Area, Tianjin, People's Republic of China.,Ministry of Education, The Key Laboratory of Molecular Microbiology and Technology, Tianjin, People's Republic of China.,Tianjin Key Laboratory of Microbial Functional Genomics, Tianjin, People's Republic of China.,State Key Laboratory of Medicinal Chemical Biology, Nankai University, Tianjin, People's Republic of China
| | - Peter R Reeves
- School of Life and Environmental Sciences, University of Sydney, Sydney, NSW, Australia
| |
Collapse
|
7
|
Zehady AK, Fordham BG, Ogg JG. Integrated species-phenon trees: visualizing infraspecific diversity within lineages. Sci Rep 2019; 9:18968. [PMID: 31831804 PMCID: PMC6908663 DOI: 10.1038/s41598-019-55435-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2019] [Accepted: 10/25/2019] [Indexed: 12/03/2022] Open
Abstract
The unprecedented detail with which contemporary molecular phylogenetics are visualizing infraspecific relationships within living species and species complexes cannot as yet be reliably extended into deep time. Yet paleontological systematics has routinely dealt in (mainly) morphotaxa envisaged in various ways to have been components of past species lineages. Bridging these perspectives can only enrich both. We present a visualization tool that digitally depicts infraspecific diversity within species through deep time. Our integrated species-phenon tree merges ancestor-descendant trees for fossil morphotaxa (phena) into reconstructed phylogenies of lineages (species) by expanding the latter into "species boxes" and placing the phenon trees inside. A key programming strategy to overcome the lack of a simple overall parent-child hierarchy in the integrated tree has been the progressive population of a species-phenon relationship map which then provides the graphical footprint for the overarching species boxes. Our initial case has been limited to planktonic foraminfera via Aze & others' important macroevolutionary dataset. The tool could potentially be appropriated for other organisms, to detail other kinds of infraspecific granularity within lineages, or more generally to visualize two nested but loosely coupled trees.
Collapse
Affiliation(s)
- Abdullah Khan Zehady
- Department of Earth, Atmospheric, and Planetary Sciences, Purdue University, West Lafayette, IN, USA
| | - Barry G Fordham
- Research School of Earth Sciences, Australian National University, Canberra, ACT, Australia.
| | - James G Ogg
- Department of Earth, Atmospheric, and Planetary Sciences, Purdue University, West Lafayette, IN, USA
- State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Chengdu University of Technology, Chengdu, Sichuan, China
| |
Collapse
|
8
|
Jabbar A, Phelan JE, de Sessions PF, Khan TA, Rahman H, Khan SN, Cantillon DM, Wildner LM, Ali S, Campino S, Waddell SJ, Clark TG. Whole genome sequencing of drug resistant Mycobacterium tuberculosis isolates from a high burden tuberculosis region of North West Pakistan. Sci Rep 2019; 9:14996. [PMID: 31628383 PMCID: PMC6802378 DOI: 10.1038/s41598-019-51562-6] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2019] [Accepted: 10/03/2019] [Indexed: 12/22/2022] Open
Abstract
Tuberculosis (TB), caused by Mycobacterium tuberculosis bacteria, is a leading infectious cause of mortality worldwide, including in Pakistan. Drug resistant M. tuberculosis is an emerging threat for TB control, making it important to detect the underlying genetic mutations, and thereby inform treatment decision making and prevent transmission. Whole genome sequencing has emerged as the new diagnostic to reliably predict drug resistance within a clinically relevant time frame, and its deployment will have the greatest impact on TB control in highly endemic regions. To evaluate the mutations leading to drug resistance and to assess for evidence of the transmission of resistant strains, 81 M. tuberculosis samples from Khyber Pakhtunkhwa province (North West Pakistan) were subjected to whole genome sequencing and standard drug susceptibility testing for eleven anti-TB drugs. We found the majority of M. tuberculosis isolates were the CAS/Delhi strain-type (lineage 3; n = 57; 70.4%) and multi-drug resistant (MDR; n = 62; 76.5%). The most frequent resistance mutations were observed in the katG and rpoB genes, conferring resistance to isoniazid and rifampicin respectively. Mutations were also observed in genes conferring resistance to other first and second-line drugs, including in pncA (pyrazinamide), embB (ethambutol), gyrA (fluoroquinolones), rrs (aminoglycosides), rpsL, rrs and giB (streptomycin) loci. Whilst the majority of mutations have been reported in global datasets, we describe unreported putative resistance markers in katG, ethA (ethionamide), gyrA and gyrB (fluoroquinolones), and pncA. Analysis of the mutations revealed that acquisition of rifampicin resistance often preceded isoniazid in our isolates. We also observed a high proportion (17.6%) of pre-MDR isolates with fluoroquinolone resistance markers, potentially due to unregulated anti-TB drug use. Our isolates were compared to previously sequenced strains from Pakistan in a combined phylogenetic tree analysis. The presence of lineage 2 was only observed in our isolates. Using a cut-off of less than ten genome-wide mutation differences between isolates, a transmission analysis revealed 18 M. tuberculosis isolates clustering within eight networks, thereby providing evidence of drug-resistant TB transmission in the Khyber Pakhtunkhwa province. Overall, we have demonstrated that drug-resistant TB isolates are circulating and transmitted in North West Pakistan. Further, we have shown the usefulness of whole genome sequencing as a diagnostic tool for characterizing M. tuberculosis isolates, which will assist future epidemiological studies and disease control activities in Pakistan.
Collapse
Affiliation(s)
- Abdul Jabbar
- Department of Medical Lab Technology, University of Haripur, Haripur, Pakistan.
- Department of Microbiology, Kohat University of Science and Technology, Kohat, Pakistan.
| | - Jody E Phelan
- Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, Keppel Street, London, WC1E 7HT, UK
| | | | - Taj Ali Khan
- Department of Microbiology, Kohat University of Science and Technology, Kohat, Pakistan
| | - Hazir Rahman
- Department of Microbiology, Abdul Wali Khan University, Mardan, Pakistan
| | - Sadiq Noor Khan
- Department of Medical Lab Technology, University of Haripur, Haripur, Pakistan
| | - Daire M Cantillon
- Department of Global Health and Infection, Brighton and Sussex Medical School, University of Sussex, Brighton, BN1 9PX, UK
| | - Leticia Muraro Wildner
- Department of Global Health and Infection, Brighton and Sussex Medical School, University of Sussex, Brighton, BN1 9PX, UK
| | - Sajid Ali
- Provincial Tuberculosis Reference Laboratory, Hayatabad Medical Complex Peshawar, Khyber Pakhtunkhwa, Pakistan
| | - Susana Campino
- Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, Keppel Street, London, WC1E 7HT, UK
| | - Simon J Waddell
- Department of Global Health and Infection, Brighton and Sussex Medical School, University of Sussex, Brighton, BN1 9PX, UK
| | - Taane G Clark
- Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, Keppel Street, London, WC1E 7HT, UK.
- Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, Keppel Street, London, WC1E 7HT, UK.
| |
Collapse
|
9
|
Yang Y, Niehaus KE, Walker TM, Iqbal Z, Walker AS, Wilson DJ, Peto TEA, Crook DW, Smith EG, Zhu T, Clifton DA. Machine learning for classifying tuberculosis drug-resistance from DNA sequencing data. Bioinformatics 2019; 34:1666-1671. [PMID: 29240876 PMCID: PMC5946815 DOI: 10.1093/bioinformatics/btx801] [Citation(s) in RCA: 66] [Impact Index Per Article: 13.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2017] [Accepted: 12/10/2017] [Indexed: 01/07/2023] Open
Abstract
Motivation Correct and rapid determination of Mycobacterium tuberculosis (MTB) resistance against available tuberculosis (TB) drugs is essential for the control and management of TB. Conventional molecular diagnostic test assumes that the presence of any well-studied single nucleotide polymorphisms is sufficient to cause resistance, which yields low sensitivity for resistance classification. Summary Given the availability of DNA sequencing data from MTB, we developed machine learning models for a cohort of 1839 UK bacterial isolates to classify MTB resistance against eight anti-TB drugs (isoniazid, rifampicin, ethambutol, pyrazinamide, ciprofloxacin, moxifloxacin, ofloxacin, streptomycin) and to classify multi-drug resistance. Results Compared to previous rules-based approach, the sensitivities from the best-performing models increased by 2-4% for isoniazid, rifampicin and ethambutol to 97% (P < 0.01), respectively; for ciprofloxacin and multi-drug resistant TB, they increased to 96%. For moxifloxacin and ofloxacin, sensitivities increased by 12 and 15% from 83 and 81% based on existing known resistance alleles to 95% and 96% (P < 0.01), respectively. Particularly, our models improved sensitivities compared to the previous rules-based approach by 15 and 24% to 84 and 87% for pyrazinamide and streptomycin (P < 0.01), respectively. The best-performing models increase the area-under-the-ROC curve by 10% for pyrazinamide and streptomycin (P < 0.01), and 4–8% for other drugs (P < 0.01). Availability and implementation The details of source code are provided at http://www.robots.ox.ac.uk/~davidc/code.php. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Yang Yang
- Department of Engineering Science, Institute of Biomedical Engineering, University of Oxford, Oxford, OX3 7DQ, UK
| | - Katherine E Niehaus
- Department of Engineering Science, Institute of Biomedical Engineering, University of Oxford, Oxford, OX3 7DQ, UK
| | - Timothy M Walker
- Nuffield Department of Medicine, University of Oxford, Oxford, OX3 7BN, UK
| | - Zamin Iqbal
- Nuffield Department of Medicine, University of Oxford, Oxford, OX3 7BN, UK
| | - A Sarah Walker
- Nuffield Department of Medicine, University of Oxford, Oxford, OX3 7BN, UK.,NIHR Oxford Biomedical Research Centre, Oxford, UK
| | - Daniel J Wilson
- Nuffield Department of Medicine, University of Oxford, Oxford, OX3 7BN, UK
| | - Tim E A Peto
- Nuffield Department of Medicine, University of Oxford, Oxford, OX3 7BN, UK.,National Infection Service, Public Health England, Colindale, London, UK
| | - Derrick W Crook
- Nuffield Department of Medicine, University of Oxford, Oxford, OX3 7BN, UK.,NIHR Oxford Biomedical Research Centre, Oxford, UK.,National Infection Service, Public Health England, Colindale, London, UK
| | | | - Tingting Zhu
- Department of Engineering Science, Institute of Biomedical Engineering, University of Oxford, Oxford, OX3 7DQ, UK
| | - David A Clifton
- Department of Engineering Science, Institute of Biomedical Engineering, University of Oxford, Oxford, OX3 7DQ, UK
| |
Collapse
|
10
|
Oppong YEA, Phelan J, Perdigão J, Machado D, Miranda A, Portugal I, Viveiros M, Clark TG, Hibberd ML. Genome-wide analysis of Mycobacterium tuberculosis polymorphisms reveals lineage-specific associations with drug resistance. BMC Genomics 2019; 20:252. [PMID: 30922221 PMCID: PMC6440112 DOI: 10.1186/s12864-019-5615-3] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2018] [Accepted: 03/15/2019] [Indexed: 12/30/2022] Open
Abstract
Background Continuing evolution of the Mycobacterium tuberculosis (Mtb) complex genomes associated with resistance to anti-tuberculosis drugs is threatening tuberculosis disease control efforts. Both multi- and extensively drug resistant Mtb (MDR and XDR, respectively) are increasing in prevalence, but the full set of Mtb genes involved are not known. There is a need for increased sensitivity of genome-wide approaches in order to elucidate the genetic basis of anti-microbial drug resistance and gain a more detailed understanding of Mtb genome evolution in a context of widespread antimicrobial therapy. Population structure within the Mtb complex, due to clonal expansion, lack of lateral gene transfer and low levels of recombination between lineages, may be reducing statistical power to detect drug resistance associated variants. Results To investigate the effect of lineage-specific effects on the identification of drug resistance associations, we applied genome-wide association study (GWAS) and convergence-based (PhyC) methods to multiple drug resistance phenotypes of a global dataset of Mtb lineages 2 and 4, using both lineage-wise and combined approaches. We identify both well-established drug resistance variants and novel associations; uniquely identifying associations for both lineage-specific and -combined GWAS analyses. We report 17 potential novel associations between antimicrobial resistance phenotypes and Mtb genomic variants. Conclusions For GWAS, both lineage-specific and -combined analyses are useful, whereas PhyC may perform better in contexts of greater diversity. Unique associations with XDR in lineage-specific analyses provide evidence of diverging evolutionary trajectories between lineages 2 and 4 in response to antimicrobial drug therapy. Electronic supplementary material The online version of this article (10.1186/s12864-019-5615-3) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Yaa E A Oppong
- Pathogen Molecular Biology Department, Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, Keppel Street, London, WC1E 7HT, UK.
| | - Jody Phelan
- Pathogen Molecular Biology Department, Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, Keppel Street, London, WC1E 7HT, UK
| | - João Perdigão
- iMed.ULisboa - Research Institute for Medicines, Faculdade de Farmácia, Universidade de Lisboa, Lisbon, Portugal
| | - Diana Machado
- Global Health and Tropical Medicine, Instituto de Higiene e Medicina Tropical, Universidade Nova de Lisboa, UNL, Lisbon, Portugal
| | - Anabela Miranda
- National Mycobacterium Reference Laboratory, Porto, Portugal
| | - Isabel Portugal
- iMed.ULisboa - Research Institute for Medicines, Faculdade de Farmácia, Universidade de Lisboa, Lisbon, Portugal
| | - Miguel Viveiros
- Global Health and Tropical Medicine, Instituto de Higiene e Medicina Tropical, Universidade Nova de Lisboa, UNL, Lisbon, Portugal
| | - Taane G Clark
- Pathogen Molecular Biology Department, Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, Keppel Street, London, WC1E 7HT, UK.,Faculty of Epidemiology and Population Health, LSHTM, London, UK
| | - Martin L Hibberd
- Pathogen Molecular Biology Department, Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, Keppel Street, London, WC1E 7HT, UK
| |
Collapse
|
11
|
An integrated whole genome analysis of Mycobacterium tuberculosis reveals insights into relationship between its genome, transcriptome and methylome. Sci Rep 2019; 9:5204. [PMID: 30914757 PMCID: PMC6435705 DOI: 10.1038/s41598-019-41692-2] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2018] [Accepted: 03/14/2019] [Indexed: 01/01/2023] Open
Abstract
Human tuberculosis disease (TB), caused by Mycobacterium tuberculosis (Mtb), is a complex disease, with a spectrum of outcomes. Genomic, transcriptomic and methylation studies have revealed differences between Mtb lineages, likely to impact on transmission, virulence and drug resistance. However, so far no studies have integrated sequence-based genomic, transcriptomic and methylation characterisation across a common set of samples, which is critical to understand how DNA sequence and methylation affect RNA expression and, ultimately, Mtb pathogenesis. Here we perform such an integrated analysis across 22 M. tuberculosis clinical isolates, representing ancient (lineage 1) and modern (lineages 2 and 4) strains. The results confirm the presence of lineage-specific differential gene expression, linked to specific SNP-based expression quantitative trait loci: with 10 eQTLs involving SNPs in promoter regions or transcriptional start sites; and 12 involving potential functional impairment of transcriptional regulators. Methylation status was also found to have a role in transcription, with evidence of differential expression in 50 genes across lineage 4 samples. Lack of methylation was associated with three novel variants in mamA, likely to cause loss of function of this enzyme. Overall, our work shows the relationship of DNA sequence and methylation to RNA expression, and differences between ancient and modern lineages. Further studies are needed to verify the functional consequences of the identified mechanisms of gene expression regulation.
Collapse
|
12
|
Chernyaeva E, Rotkevich M, Krasheninnikova K, Yurchenko A, Vyazovaya A, Mokrousov I, Solovieva N, Zhuravlev V, Yablonsky P, O'Brien SJ. Whole-Genome Analysis of Mycobacterium tuberculosis from Patients with Tuberculous Spondylitis, Russia. Emerg Infect Dis 2019; 24:579-583. [PMID: 29460750 PMCID: PMC5823328 DOI: 10.3201/eid2403.170151] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022] Open
Abstract
Whole-genome analysis of Mycobacterium tuberculosis isolates collected in Russia (N = 71) from patients with tuberculous spondylitis supports a detailed characterization of pathogen strain distributions and drug resistance phenotype, plus distinguished occurrence and association of known resistance mutations. We identify known and novel genome determinants related to bacterial virulence, pathogenicity, and drug resistance.
Collapse
|
13
|
Phelan J, O’Sullivan DM, Machado D, Ramos J, Whale AS, O’Grady J, Dheda K, Campino S, McNerney R, Viveiros M, Huggett JF, Clark TG. The variability and reproducibility of whole genome sequencing technology for detecting resistance to anti-tuberculous drugs. Genome Med 2016; 8:132. [PMID: 28003022 PMCID: PMC5178084 DOI: 10.1186/s13073-016-0385-x] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2016] [Accepted: 11/30/2016] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND The emergence of resistance to anti-tuberculosis drugs is a serious and growing threat to public health. Next-generation sequencing is rapidly gaining traction as a diagnostic tool for investigating drug resistance in Mycobacterium tuberculosis to aid treatment decisions. However, there are few little data regarding the precision of such sequencing for assigning resistance profiles. METHODS We investigated two sequencing platforms (Illumina MiSeq, Ion Torrent PGM™) and two rapid analytic pipelines (TBProfiler, Mykrobe predictor) using a well characterised reference strain (H37Rv) and clinical isolates from patients with tuberculosis resistant to up to 13 drugs. Results were compared to phenotypic drug susceptibility testing. To assess analytical robustness individual DNA samples were subjected to repeated sequencing. RESULTS The MiSeq and Ion PGM systems accurately predicted drug-resistance profiles and there was high reproducibility between biological and technical sample replicates. Estimated variant error rates were low (MiSeq 1 per 77 kbp, Ion PGM 1 per 41 kbp) and genomic coverage high (MiSeq 51-fold, Ion PGM 53-fold). MiSeq provided superior coverage in GC-rich regions, which translated into incremental detection of putative genotypic drug-specific resistance, including for resistance to para-aminosalicylic acid and pyrazinamide. The TBProfiler bioinformatics pipeline was concordant with reported phenotypic susceptibility for all drugs tested except pyrazinamide and para-aminosalicylic acid, with an overall concordance of 95.3%. When using the Mykrobe predictor concordance with phenotypic testing was 73.6%. CONCLUSIONS We have demonstrated high comparative reproducibility of two sequencing platforms, and high predictive ability of the TBProfiler mutation library and analytical pipeline, when profiling resistance to first- and second-line anti-tuberculosis drugs. However, platform-specific variability in coverage of some genome regions may have implications for predicting resistance to specific drugs. These findings may have implications for future clinical practice and thus deserve further scrutiny, set within larger studies and using updated mutation libraries.
Collapse
Affiliation(s)
- Jody Phelan
- Department of Pathogen Molecular Biology, Faculty of Infectious and Tropical Diseases, London School of Hygiene & Tropical Medicine, Keppel Street, WC1E 7HT London, UK
| | | | - Diana Machado
- Unidade de Microbiologia Médica, Global Health and Tropical Medicine, GHTM, Instituto de Higiene e Medicina Tropical, IHMT, Universidade NOVA de Lisboa, UNL, Lisbon, Portugal
| | - Jorge Ramos
- Unidade de Microbiologia Médica, Global Health and Tropical Medicine, GHTM, Instituto de Higiene e Medicina Tropical, IHMT, Universidade NOVA de Lisboa, UNL, Lisbon, Portugal
| | - Alexandra S. Whale
- Molecular Biology, LGC Ltd, Queens Road, Teddington, Middlesex TW11 0LY UK
| | - Justin O’Grady
- Norwich Medical School, University of East Anglia, Norwich Research Park, Norwich, NR4 7TJ UK
| | - Keertan Dheda
- Division of Pulmonary Medicine and UCT Lung Institute, Lung Infection and Immunity Unit, University of Cape Town, Groote Schuur Hospital, Observatory, 7925, Cape Town, South Africa
| | - Susana Campino
- Department of Pathogen Molecular Biology, Faculty of Infectious and Tropical Diseases, London School of Hygiene & Tropical Medicine, Keppel Street, WC1E 7HT London, UK
| | - Ruth McNerney
- Division of Pulmonary Medicine and UCT Lung Institute, Lung Infection and Immunity Unit, University of Cape Town, Groote Schuur Hospital, Observatory, 7925, Cape Town, South Africa
| | - Miguel Viveiros
- Unidade de Microbiologia Médica, Global Health and Tropical Medicine, GHTM, Instituto de Higiene e Medicina Tropical, IHMT, Universidade NOVA de Lisboa, UNL, Lisbon, Portugal
| | - Jim F. Huggett
- Molecular Biology, LGC Ltd, Queens Road, Teddington, Middlesex TW11 0LY UK
- School of Biosciences & Medicine, Faculty of Health & Medical Science, University of Surrey, Guildford, GU2 7XH UK
| | - Taane G. Clark
- Department of Pathogen Molecular Biology, Faculty of Infectious and Tropical Diseases, London School of Hygiene & Tropical Medicine, Keppel Street, WC1E 7HT London, UK
- Faculty of Epidemiology and Population Health, London School of Hygiene & Tropical Medicine, WC1E 7HT London, UK
| |
Collapse
|
14
|
Brown TS, Narechania A, Walker JR, Planet PJ, Bifani PJ, Kolokotronis SO, Kreiswirth BN, Mathema B. Genomic epidemiology of Lineage 4 Mycobacterium tuberculosis subpopulations in New York City and New Jersey, 1999-2009. BMC Genomics 2016; 17:947. [PMID: 27871225 PMCID: PMC5117616 DOI: 10.1186/s12864-016-3298-6] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2016] [Accepted: 11/15/2016] [Indexed: 12/31/2022] Open
Abstract
Background Whole genome sequencing (WGS) has rapidly become an important research tool in tuberculosis epidemiology and is likely to replace many existing methods in public health microbiology in the near future. WGS-based methods may be particularly useful in areas with less diverse Mycobacterium tuberculosis populations, such as New York City, where conventional genotyping is often uninformative and field epidemiology often difficult. This study applies four candidate strategies for WGS-based identification of emerging M. tuberculosis subpopulations, employing both phylogenomic and population genetics methods. Results M. tuberculosis subpopulations in New York City and New Jersey can be distinguished via phylogenomic reconstruction, evidence of demographic expansion and subpopulation-specific signatures of selection, and by determination of subgroup-defining nucleotide substitutions. These methods identified known historical outbreak clusters and previously unidentified subpopulations within relatively monomorphic M. tuberculosis endemic clone groups. Neutrality statistics based on the site frequency spectrum were less useful for identifying M. tuberculosis subpopulations, likely due to the low levels of informative genetic variation in recently diverged isolate groups. In addition, we observed that isolates from New York City endemic clone groups have acquired multiple non-synonymous SNPs in virulence- and growth-associated pathways, and relatively few mutations in drug resistance-associated genes, suggesting that overall pathoadaptive fitness, rather than the acquisition of drug resistance mutations, has played a central role in the evolutionary history and epidemiology of M. tuberculosis subpopulations in New York City. Conclusions Our results demonstrate that some but not all WGS-based methods are useful for detection of emerging M. tuberculosis clone groups, and support the use of phylogenomic reconstruction in routine tuberculosis laboratory surveillance, particularly in areas with relatively less diverse M. tuberculosis populations. Our study also supports the use of wider-reaching phylogenomic and population genomic methods in tuberculosis public health practice, which can support tuberculosis control activities by identifying genetic polymorphisms contributing to epidemiological success in local M. tuberculosis populations and possibly explain why certain isolate groups are apparently more successful in specific host populations. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-3298-6) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Tyler S Brown
- Department of Medicine, Columbia University, New York, NY, USA
| | - Apurva Narechania
- Sackler Institute for Comparative Genomics, American Museum of Natural History, New York, NY, USA
| | - John R Walker
- The Genomic Institute of the Novartis Research Foundation, San Diego, CA, USA
| | - Paul J Planet
- Department of Pediatrics, Division of Infectious Diseases, The Children's Hospital of Philadelphia, Philadelphia, PA, USA
| | - Pablo J Bifani
- Novartis Institute for Tropical Diseases, Singapore, Singapore
| | - Sergios-Orestis Kolokotronis
- Department of Epidemiology and Biostatistics, School of Public Health, SUNY Downstate Medical Center, Brooklyn, NY, USA
| | | | - Barun Mathema
- Department of Epidemiology, Mailman School of Public Health, Columbia University, New York, USA.
| |
Collapse
|
15
|
Refrégier G, Abadia E, Matsumoto T, Ano H, Takashima T, Tsuyuguchi I, Aktas E, Cömert F, Gomgnimbou MK, Panaiotov S, Phelan J, Coll F, McNerney R, Pain A, Clark TG, Sola C. Turkish and Japanese Mycobacterium tuberculosis sublineages share a remote common ancestor. INFECTION GENETICS AND EVOLUTION 2016; 45:461-473. [PMID: 27746295 DOI: 10.1016/j.meegid.2016.10.009] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/14/2016] [Revised: 09/23/2016] [Accepted: 10/12/2016] [Indexed: 12/19/2022]
Abstract
Two geographically distant M. tuberculosis sublineages, Tur from Turkey and T3-Osaka from Japan, exhibit partially identical genotypic signatures (identical 12-loci MIRU-VNTR profiles, distinct spoligotyping patterns). We investigated T3-Osaka and Tur sublineages characteristics and potential genetic relatedness, first using MIRU-VNTR locus analysis on 21 and 25 samples of each sublineage respectively, and second comparing Whole Genome Sequences of 8 new samples to public data from 45 samples uncovering human tuberculosis diversity. We then tried to date their Most Recent Common Ancestor (MRCA) using three calibrations of SNP accumulation rate (long-term=0.03SNP/genome/year, derived from a tuberculosis ancestor of around 70,000years old; intermediate=0.2SNP/genome/year derived from a Peruvian mummy; short-term=0.5SNP/genome/year). To disentangle between these scenarios, we confronted the corresponding divergence times with major human history events and knowledge on human genetic divergence. We identified relatively high intrasublineage diversity for both T3-Osaka and Tur. We definitively proved their monophyly; the corresponding super-sublineage (referred to as "T3-Osa-Tur") shares a common ancestor with T3-Ethiopia and Ural sublineages but is only remotely related to other Euro-American sublineages such as X, LAM, Haarlem and S. The evolutionary scenario based on long-term evolution rate being valid until T3-Osa-Tur MRCA was not supported by Japanese fossil data. The evolutionary scenario relying on short-term evolution rate since T3-Osa-Tur MRCA was contradicted by human history and potential traces of past epidemics. T3-Osaka and Tur sublineages were found likely to have diverged between 800y and 2000years ago, potentially at the time of Mongol Empire. Altogether, this study definitively proves a strong genetic link between Turkish and Japanese tuberculosis. It provides a first hypothesis for calibrating TB Euro-American lineage molecular clock; additional studies are needed to reliably date events corresponding to intermediate depths in tuberculosis phylogeny.
Collapse
Affiliation(s)
- Guislaine Refrégier
- Institute for Integrative Biology of the Cell (I2BC), CEA, CNRS, Univ Paris-Sud, Université Paris-Saclay, 91198, Gif-sur-Yvette cedex, France.
| | - Edgar Abadia
- Instituto Venezolano de Investigaciones Científicas (IVIC), Caracas, Venezuela
| | - Tomoshige Matsumoto
- Department of Clinical Research and Development, Osaka Prefectural Hospital Organization, Osaka Prefectural Medical Center for Respiratory and Allergic Diseases, Habikino-city, Japan
| | - Hiromi Ano
- Department of Clinical Research and Development, Osaka Prefectural Hospital Organization, Osaka Prefectural Medical Center for Respiratory and Allergic Diseases, Habikino-city, Japan
| | - Tetsuya Takashima
- Department of Clinical Research and Development, Osaka Prefectural Hospital Organization, Osaka Prefectural Medical Center for Respiratory and Allergic Diseases, Habikino-city, Japan
| | - Izuo Tsuyuguchi
- Department of Clinical Research and Development, Osaka Prefectural Hospital Organization, Osaka Prefectural Medical Center for Respiratory and Allergic Diseases, Habikino-city, Japan
| | - Elif Aktas
- Şişli Etfal Research and Training Hopital, Istanbul, Turkey
| | - Füsun Cömert
- Faculty of Medicine, Bülent Ecevit University, Zonguldak, Turkey
| | - Michel Kireopori Gomgnimbou
- Institute for Integrative Biology of the Cell (I2BC), CEA, CNRS, Univ Paris-Sud, Université Paris-Saclay, 91198, Gif-sur-Yvette cedex, France
| | - Stefan Panaiotov
- National Center of Parasitic and Infectious Diseases, Sofia, Bulgaria
| | - Jody Phelan
- Faculty of Infectious and Tropical Diseases, London School of Hygiene & Tropical Medicine, London, UK
| | - Francesc Coll
- Faculty of Infectious and Tropical Diseases, London School of Hygiene & Tropical Medicine, London, UK
| | - Ruth McNerney
- Faculty of Infectious and Tropical Diseases, London School of Hygiene & Tropical Medicine, London, UK; Infection and Immunity Unit, UCT Lung Institute, University of Cape Town, Old Main Building, Groote Schuur Hospital, Cape Town,South Africa
| | - Arnab Pain
- Pathogen Genomics Group, Biological, Environmental Sciences and Engineering Division, King Abdullah University of Science and Technology, Thuwal, Saudi Arabia
| | - Taane G Clark
- Faculty of Infectious and Tropical Diseases, London School of Hygiene & Tropical Medicine, London, UK
| | - Christophe Sola
- Institute for Integrative Biology of the Cell (I2BC), CEA, CNRS, Univ Paris-Sud, Université Paris-Saclay, 91198, Gif-sur-Yvette cedex, France
| |
Collapse
|
16
|
Bioinformatics tools and databases for whole genome sequence analysis of Mycobacterium tuberculosis. INFECTION GENETICS AND EVOLUTION 2016; 45:359-368. [PMID: 27637931 DOI: 10.1016/j.meegid.2016.09.013] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/17/2016] [Revised: 09/12/2016] [Accepted: 09/12/2016] [Indexed: 11/24/2022]
Abstract
Tuberculosis (TB) is an infectious disease of global public health importance caused by Mycobacterium tuberculosis complex (MTC) in which M. tuberculosis (Mtb) is the major causative agent. Recent advancements in genomic technologies such as next generation sequencing have enabled high throughput cost-effective generation of whole genome sequence information from Mtb clinical isolates, providing new insights into the evolution, genomic diversity and transmission of the Mtb bacteria, including molecular mechanisms of antibiotic resistance. The large volume of sequencing data generated however necessitated effective and efficient management, storage, analysis and visualization of the data and results through development of novel and customized bioinformatics software tools and databases. In this review, we aim to provide a comprehensive survey of the current freely available bioinformatics software tools and publicly accessible databases for genomic analysis of Mtb for identifying disease transmission in molecular epidemiology and in rapid determination of the antibiotic profiles of clinical isolates for prompt and optimal patient treatment.
Collapse
|
17
|
Faksri K, Tan JH, Disratthakit A, Xia E, Prammananan T, Suriyaphol P, Khor CC, Teo YY, Ong RTH, Chaiprasert A. Whole-Genome Sequencing Analysis of Serially Isolated Multi-Drug and Extensively Drug Resistant Mycobacterium tuberculosis from Thai Patients. PLoS One 2016; 11:e0160992. [PMID: 27518818 PMCID: PMC4982626 DOI: 10.1371/journal.pone.0160992] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2016] [Accepted: 07/28/2016] [Indexed: 01/08/2023] Open
Abstract
Multi-drug and extensively drug-resistant tuberculosis (MDR and XDR-TB) are problems that threaten public health worldwide. Only some genetic markers associated with drug-resistant TB are known. Whole-genome sequencing (WGS) is a promising tool for distinguishing between re-infection and persistent infection in isolates taken at different times from a single patient, but has not yet been applied in MDR and XDR-TB. We aim to detect genetic markers associated with drug resistance and distinguish between reinfection and persistent infection from MDR and XDR-TB patients based on WGS analysis. Samples of Mycobacterium tuberculosis (n = 7), serially isolated from 2 MDR cases and 1 XDR-TB case, were retrieved from Siriraj Hospital, Bangkok. The WGS analysis used an Illumina Miseq sequencer. In cases of persistent infection, MDR-TB isolates differed at an average of 2 SNPs across the span of 2–9 months whereas in the case of reinfection, isolates differed at 61 SNPs across 2 years. Known genetic markers associated with resistance were detected from strains susceptible to streptomycin (2/7 isolates), p-aminosalicylic acid (3/7 isolates) and fluoroquinolone drugs. Among fluoroquinolone drugs, ofloxacin had the highest phenotype-genotype concordance (6/7 isolates), whereas gatifloxcain had the lowest (3/7 isolates). A putative candidate SNP in Rv2477c associated with kanamycin and amikacin resistance was suggested for further validation. WGS provided comprehensive results regarding molecular epidemiology, distinguishing between persistent infection and reinfection in M/XDR-TB and potentially can be used for detection of novel mutations associated with drug resistance.
Collapse
Affiliation(s)
- Kiatichai Faksri
- Department of Microbiology Faculty of Medicine, Khon Kaen University, Khon Kaen, Thailand
- Research and Diagnostic Center for Emerging Infectious Diseases (RCEID), Khon Kaen University, Khon Kaen, Thailand
- * E-mail: (KF); (AC)
| | - Jun Hao Tan
- Saw Swee Hock School of Public Health, National University of Singapore, Singapore
| | - Areeya Disratthakit
- Department of Microbiology, Faculty of Medicine Siriraj Hospital, Mahidol University, Bangkok, Thailand
| | - Eryu Xia
- NUS Graduate School for Integrative Sciences and Engineering, National University of Singapore, Singapore
| | - Therdsak Prammananan
- National Center for Genetic Engineering and Biotechnology, National Science and Technology Development Agency, Ministry of Science and Technology, Pathum Thani, Thailand
| | - Prapat Suriyaphol
- Bioinformatics and Data Management for Research Unit, Office for Research and Development, Faculty of Medicine Siriraj Hospital, Mahidol University, Bangkok, Thailand
| | | | - Yik-Ying Teo
- Saw Swee Hock School of Public Health, National University of Singapore, Singapore
- NUS Graduate School for Integrative Sciences and Engineering, National University of Singapore, Singapore
- Genome Institute of Singapore, Singapore
- Department of Statistics and Applied Probability, National University of Singapore, Singapore
- Life Sciences Institute, National University of Singapore, Singapore
| | - Rick Twee-Hee Ong
- Saw Swee Hock School of Public Health, National University of Singapore, Singapore
| | - Angkana Chaiprasert
- Department of Microbiology, Faculty of Medicine Siriraj Hospital, Mahidol University, Bangkok, Thailand
- * E-mail: (KF); (AC)
| |
Collapse
|
18
|
Coker OO, Chaiprasert A, Ngamphiw C, Tongsima S, Regmi SM, Clark TG, Ong RTH, Teo YY, Prammananan T, Palittapongarnpim P. Genetic signatures of Mycobacterium tuberculosis Nonthaburi genotype revealed by whole genome analysis of isolates from tuberculous meningitis patients in Thailand. PeerJ 2016; 4:e1905. [PMID: 27114869 PMCID: PMC4841212 DOI: 10.7717/peerj.1905] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2015] [Accepted: 03/16/2016] [Indexed: 01/27/2023] Open
Abstract
Genome sequencing plays a key role in understanding the genetic diversity of Mycobacterium tuberculosis (M.tb). The genotype-specific character of M. tb contributes to tuberculosis severity and emergence of drug resistance. Strains of M. tb complex can be classified into seven lineages. The Nonthaburi (NB) genotype, belonging to the Indo-Oceanic lineage (lineage 1), has a unique spoligotype and IS6110-RFLP pattern but has not previously undergone a detailed whole genome analysis. In addition, there is not much information available on the whole genome analysis of M. tb isolates from tuberculous meningitis (TBM) patients in public databases. Isolates CSF3053, 46-5069 and 43-13838 of NB genotype were obtained from the cerebrospinal fluids of TBM Thai patients in Siriraj Hospital, Bangkok. The whole genomes were subjected to high throughput sequencing. The sequence data of each isolate were assembled into draft genome. The sequences were also aligned to reference genome, to determine genomic variations. Single nucleotide polymorphisms (SNPs) were obtained and grouped according to the functions of the genes containing them. They were compared with SNPs from 1,601 genomes, representing the seven lineages of M. tb complex, to determine the uniqueness of NB genotype. Susceptibility to first-line, second-line and other antituberculosis drugs were determined and related to the SNPs previously reported in drug-resistant related genes. The assembled genomes have an average size of 4,364,461 bp, 4,154 genes, 48 RNAs and 64 pseudogenes. A 500 base pairs deletion, which includes ppe50, was found in all isolates. RD239, specific for members of Indo Oceanic lineage, and RD147c were identified. A total of 2,202 SNPs were common to the isolates and used to classify the NB strains as members of sublineage 1.2.1. Compared with 1,601 genomes from the seven lineages of M. tb complex, mutation G2342203C was found novel to the isolates in this study. Three mutations (T28910C, C1180580T and C152178T) were found only in Thai NB isolates, including isolates from previous study. Although drug susceptibility tests indicated pan-susceptibility, non-synonymous SNPs previously reported to be associated with resistance to anti-tuberculous drugs; isoniazid, ethambutol, and ethionamide were identified in all the isolates. Non-synonymous SNPs were found in virulence genes such as the genes playing roles in apoptosis inhibition and phagosome arrest. We also report polymorphisms in essential genes, efflux pumps associated genes and genes with known epitopes. The analysis of the TBM isolates and the availability of the variations obtained will provide additional resources for global comparison of isolates from pulmonary tuberculosis and TBM. It will also contribute to the richness of genomic databases towards the prediction of antibiotic resistance, level of virulence and of origin of infection.
Collapse
Affiliation(s)
- Olabisi Oluwabukola Coker
- Department of Microbiology, Faculty of Medicine Siriraj Hospital, Mahidol University , Bangkok , Thailand
| | - Angkana Chaiprasert
- Department of Microbiology, Faculty of Medicine Siriraj Hospital, Mahidol University , Bangkok , Thailand
| | - Chumpol Ngamphiw
- National Center for Genetic Engineering and Biotechnology, National Science and Technology Development Agency , Pathum Thani , Thailand
| | - Sissades Tongsima
- National Center for Genetic Engineering and Biotechnology, National Science and Technology Development Agency , Pathum Thani , Thailand
| | - Sanjib Mani Regmi
- Department of Microbiology, Gandaki Medical College , Pokhara Kaski , Nepal
| | - Taane G Clark
- Faculty of Epidemiology and Population Health, Faculty of Infectious and Tropical Diseases, London School of Hygiene & Tropical Medicine, University of London , London , United Kingdom
| | - Rick Twee Hee Ong
- Saw Swee Hock School of Public Health, National University of Singapore , Singapore
| | - Yik-Ying Teo
- Saw Swee Hock School of Public Health, National University of Singapore , Singapore
| | - Therdsak Prammananan
- Tuberculosis Research Laboratory, Medical Molecular Biology Research Unit, National Center for Genetic Engineering and Biotechnology, National Science and Technology Development Agency , Pathum Thani , Thailand
| | | |
Collapse
|
19
|
Draft Genome Sequence of Amikacin- and Kanamycin-Resistant Mycobacterium tuberculosis MT433 without rrs and eis Mutations. GENOME ANNOUNCEMENTS 2015; 3:3/6/e01363-15. [PMID: 26586896 PMCID: PMC4653798 DOI: 10.1128/genomea.01363-15] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
Abstract
We announce the draft genome sequence of amikacin- and kanamycin-resistant Mycobacterium tuberculosis MT433, which has been previously described as the strain carrying an unknown resistance mechanism.
Collapse
|
20
|
Sekizuka T, Yamashita A, Murase Y, Iwamoto T, Mitarai S, Kato S, Kuroda M. TGS-TB: Total Genotyping Solution for Mycobacterium tuberculosis Using Short-Read Whole-Genome Sequencing. PLoS One 2015; 10:e0142951. [PMID: 26565975 PMCID: PMC4643978 DOI: 10.1371/journal.pone.0142951] [Citation(s) in RCA: 68] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2015] [Accepted: 10/28/2015] [Indexed: 11/18/2022] Open
Abstract
Whole-genome sequencing (WGS) with next-generation DNA sequencing (NGS) is an increasingly accessible and affordable method for genotyping hundreds of Mycobacterium tuberculosis (Mtb) isolates, leading to more effective epidemiological studies involving single nucleotide variations (SNVs) in core genomic sequences based on molecular evolution. We developed an all-in-one web-based tool for genotyping Mtb, referred to as the Total Genotyping Solution for TB (TGS-TB), to facilitate multiple genotyping platforms using NGS for spoligotyping and the detection of phylogenies with core genomic SNVs, IS6110 insertion sites, and 43 customized loci for variable number tandem repeat (VNTR) through a user-friendly, simple click interface. This methodology is implemented with a KvarQ script to predict MTBC lineages/sublineages and potential antimicrobial resistance. Seven Mtb isolates (JP01 to JP07) in this study showing the same VNTR profile were accurately discriminated through median-joining network analysis using SNVs unique to those isolates. An additional IS6110 insertion was detected in one of those isolates as supportive genetic information in addition to core genomic SNVs. The results of in silico analyses using TGS-TB are consistent with those obtained using conventional molecular genotyping methods, suggesting that NGS short reads could provide multiple genotypes to discriminate multiple strains of Mtb, although longer NGS reads (≥300-mer) will be required for full genotyping on the TGS-TB web site. Most available short reads (~100-mer) can be utilized to discriminate the isolates based on the core genome phylogeny. TGS-TB provides a more accurate and discriminative strain typing for clinical and epidemiological investigations; NGS strain typing offers a total genotyping solution for Mtb outbreak and surveillance. TGS-TB web site: https://gph.niid.go.jp/tgs-tb/.
Collapse
Affiliation(s)
- Tsuyoshi Sekizuka
- Pathogen Genomics Center, National Institute of Infectious Diseases, Shinjyuku-ku, Tokyo, Japan
| | - Akifumi Yamashita
- Pathogen Genomics Center, National Institute of Infectious Diseases, Shinjyuku-ku, Tokyo, Japan
| | - Yoshiro Murase
- Molecular Epidemiology Division, The Research Institute of Tuberculosis/Japan Anti-Tuberculosis Association, Kiyose, Tokyo, Japan
| | - Tomotada Iwamoto
- Department of Infectious Diseases, Kobe Institute of Health, Chuo-ku, Kobe, Japan
| | - Satoshi Mitarai
- Bacteriology Division, Department of Mycobacterium Reference and Research, Research Institute of Tuberculosis, Japan Anti-Tuberculosis Association, Kiyose, Tokyo, Japan
| | - Seiya Kato
- Research Institute of Tuberculosis, Japan Anti-Tuberculosis Association, Kiyose, Tokyo, Japan
| | - Makoto Kuroda
- Pathogen Genomics Center, National Institute of Infectious Diseases, Shinjyuku-ku, Tokyo, Japan
- * E-mail:
| |
Collapse
|