51
|
Pérez-Martínez DE, Bermúdez-Hernández GA, Madrazo-Moya CF, Cancino-Muñoz I, Montero H, Licona-Cassani C, Muñiz-Salazar R, Comas I, Zenteno-Cuevas R. SNPs in Genes Related to DNA Damage Repair in Mycobacterium Tuberculosis: Their Association with Type 2 Diabetes Mellitus and Drug Resistance. Genes (Basel) 2022; 13:genes13040609. [PMID: 35456415 PMCID: PMC9029044 DOI: 10.3390/genes13040609] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2022] [Revised: 03/19/2022] [Accepted: 03/24/2022] [Indexed: 12/04/2022] Open
Abstract
Genes related to DNA damage repair in Mycobacterium tuberculosis are critical for survival and genomic diversification. The aim of this study is to compare the presence of SNPs in genes related to DNA damage repair in sensitive and drug-resistant M. tuberculosis genomes isolated from patients with and without type 2 diabetes mellitus (T2DM). We collected 399 M. tuberculosis L4 genomes from several public repositories; 224 genomes belonging to hosts without T2DM, of which 123 (54.9%) had drug sensitive tuberculosis (TB) and 101 (45.1%) had drug resistance (DR)-TB; and 175 genomes from individuals with T2DM, of which 100 (57.1%) had drug sensitive TB and 75 (42.9%) had DR-TB. The presence of SNPs in the coding regions of 65 genes related to DNA damage repair was analyzed and compared with the resistance profile and the presence/absence of T2DM in the host. The results show the phylogenetic relationships of some SNPS and L4 sub-lineages, as well as differences in the distribution of SNPs present in DNA damage repair-related genes related to the resistance profile of the infecting strain and the presence of T2DM in the host. Given these differences, it was possible to generate two discriminant functions to distinguish between drug sensitive and drug resistant genomes, as well as patients with or without T2DM.
Collapse
Affiliation(s)
- Damián E. Pérez-Martínez
- Programa de Doctorado en Ciencias de la Salud, Instituto de Ciencias de la Salud, Universidad Veracruzana, Av. Luis, Dr. Castelazo Ayala s/n, Col. Industrial Animas, Xalapa 91190, Mexico; (D.E.P.-M.); (G.A.B.-H.)
| | - Gustavo A. Bermúdez-Hernández
- Programa de Doctorado en Ciencias de la Salud, Instituto de Ciencias de la Salud, Universidad Veracruzana, Av. Luis, Dr. Castelazo Ayala s/n, Col. Industrial Animas, Xalapa 91190, Mexico; (D.E.P.-M.); (G.A.B.-H.)
| | - Carlos F. Madrazo-Moya
- Biomedical Institute of Valencia IBV-CSIC, C. de Jaume Roig, 11, 46010 Valencia, Spain; (C.F.M.-M.); (I.C.-M.); (I.C.)
| | - Irving Cancino-Muñoz
- Biomedical Institute of Valencia IBV-CSIC, C. de Jaume Roig, 11, 46010 Valencia, Spain; (C.F.M.-M.); (I.C.-M.); (I.C.)
- CIBER of Epidemiology and Public Health, 08908 Madrid, Spain
| | - Hilda Montero
- Instituto de Salud Pública, Universidad Veracruzana, Av. Luis Castelazo Ayala s/n, A.P. 57, Col. Industrial Animas, Xalapa 91190, Mexico;
| | - Cuauhtemoc Licona-Cassani
- Escuela de Ingeniería y Ciencias, Tecnológico de Monterrey, Ave. Eugenio Garza Sada 2501 Sur, Monterrey 64849, Mexico;
- Red Multidisciplinaria de Investigación en Tuberculosis, Mexico City 14080, Mexico;
- Division of Integrative Biology, The Institute for Obesity Research, Tecnológico de Monterrey, Monterrey 64849, Mexico
| | - Raquel Muñiz-Salazar
- Red Multidisciplinaria de Investigación en Tuberculosis, Mexico City 14080, Mexico;
- Laboratorio de Epidemiología y Ecología Molecular, Escuela de Ciencias de la Salud, Universidad Autónoma de Baja California, Ensenada 22890, Mexico
| | - Iñaki Comas
- Biomedical Institute of Valencia IBV-CSIC, C. de Jaume Roig, 11, 46010 Valencia, Spain; (C.F.M.-M.); (I.C.-M.); (I.C.)
- CIBER of Epidemiology and Public Health, 08908 Madrid, Spain
| | - Roberto Zenteno-Cuevas
- Instituto de Salud Pública, Universidad Veracruzana, Av. Luis Castelazo Ayala s/n, A.P. 57, Col. Industrial Animas, Xalapa 91190, Mexico;
- Red Multidisciplinaria de Investigación en Tuberculosis, Mexico City 14080, Mexico;
- Correspondence:
| |
Collapse
|
52
|
Yang T, Gan M, Liu Q, Liang W, Tang Q, Luo G, Zuo T, Guo Y, Hong C, Li Q, Tan W, Gao Q. SAM-TB: a whole genome sequencing data analysis website for detection of Mycobacterium tuberculosis drug resistance and transmission. Brief Bioinform 2022; 23:6535677. [PMID: 35211720 PMCID: PMC8921607 DOI: 10.1093/bib/bbac030] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2021] [Revised: 12/28/2021] [Accepted: 01/25/2022] [Indexed: 01/28/2023] Open
Abstract
Whole genome sequencing (WGS) can provide insight into drug-resistance, transmission chains and the identification of outbreaks, but data analysis remains an obstacle to its routine clinical use. Although several drug-resistance prediction tools have appeared, until now no website integrates drug-resistance prediction with strain genetic relationships and species identification of nontuberculous mycobacteria (NTM). We have established a free, function-rich, user-friendly online platform for MTB WGS data analysis (SAM-TB, http://samtb.szmbzx.com) that integrates drug-resistance prediction for 17 antituberculosis drugs, detection of variants, analysis of genetic relationships and NTM species identification. The accuracy of SAM-TB in predicting drug-resistance was assessed using 3177 sequenced clinical isolates with results of phenotypic drug-susceptibility tests (pDST). Compared to pDST, the sensitivity of SAM-TB for detecting multidrug-resistant tuberculosis was 93.9% [95% confidence interval (CI) 92.6–95.1%] with specificity of 96.2% (95% CI 95.2–97.1%). SAM-TB also analyzes the genetic relationships between multiple strains by reconstructing phylogenetic trees and calculating pairwise single nucleotide polymorphism (SNP) distances to identify genomic clusters. The incorporated mlstverse software identifies NTM species with an accuracy of 98.2% and Kraken2 software can detect mixed MTB and NTM samples. SAM-TB also has the capacity to share both sequence data and analysis between users. SAM-TB is a multifunctional integrated website that uses WGS raw data to accurately predict antituberculosis drug-resistance profiles, analyze genetic relationships between multiple strains and identify NTM species and mixed samples containing both NTM and MTB. SAM-TB is a useful tool for guiding both treatment and epidemiological investigation.
Collapse
Affiliation(s)
- Tingting Yang
- Key Laboratory of Medical Molecular Virology (MOE/NHC/CAMS), School of Basic Medical Sciences, Shanghai Medical College, Fudan University, Shanghai, China
| | | | | | | | | | - Geyang Luo
- Key Laboratory of Medical Molecular Virology (MOE/NHC/CAMS), School of Basic Medical Sciences, Shanghai Medical College, Shanghai Institute of Infectious Disease and Biosecurity and Shanghai Public Health Clinical Center, Fudan University, Shanghai, China
| | - Tianyu Zuo
- Fudan University, Shanghai, China. Presently, he is a data analyst in patsnap, Shanghai, China
| | - Yongchao Guo
- Shenzhen Uni-medica Technology Co., Ltd, Shenzhen, China
| | - Chuangyue Hong
- Shenzhen Center for Chronic Disease Control, Shenzhen, China
| | | | - Weiguo Tan
- Shenzhen Center for Chronic Disease Control, Shenzhen, China
| | - Qian Gao
- Key Laboratory of Medical Molecular Virology (MOE/NHC/CAMS), School of Basic Medical Sciences, Shanghai Medical College, Shanghai Institute of Infectious Disease and Biosecurity and Shanghai Public Health Clinical Center, Fudan University, Shanghai, China
| |
Collapse
|
53
|
Netikul T, Thawornwattana Y, Mahasirimongkol S, Yanai H, Maung HMW, Chongsuvivatwong V, Palittapongarnpim P. Whole-genome single nucleotide variant phylogenetic analysis of Mycobacterium tuberculosis Lineage 1 in endemic regions of Asia and Africa. Sci Rep 2022; 12:1565. [PMID: 35091638 PMCID: PMC8799649 DOI: 10.1038/s41598-022-05524-0] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2021] [Accepted: 01/11/2022] [Indexed: 12/14/2022] Open
Abstract
Mycobacterium tuberculosis (Mtb) lineage 1 (L1) contributes considerably to the disease morbidity. While whole genome sequencing (WGS) is increasingly used for studying Mtb, our understanding of genetic diversity of L1 remains limited. Using phylogenetic analysis of WGS data from endemic range in Asia and Africa, we provide an improved genotyping scheme for L1. Mapping deletion patterns of the 68 direct variable repeats (DVRs) in the CRISPR region of the genome onto the phylogeny provided supporting evidence that the CRISPR region evolves primarily by deletion, and hinted at a possible Southeast Asian origin of L1. Both phylogeny and DVR patterns clarified some relationships between different spoligotypes, and highlighted the limited resolution of spoligotyping. We identified a diverse repertoire of drug resistance mutations. Altogether, this study demonstrates the usefulness of WGS data for understanding the genetic diversity of L1, with implications for public health surveillance and TB control. It also highlights the need for more WGS studies in high-burden but underexplored regions.
Collapse
Affiliation(s)
- Thidarat Netikul
- Faculty of Medicine, Siam University, Phet Kasem Road, Bangkok, Thailand.,Pornchai Matangkasombut Center for Microbial Genomics, Department of Microbiology, Faculty of Science, Mahidol University, Rama 6 Road, Bangkok, Thailand
| | - Yuttapong Thawornwattana
- Pornchai Matangkasombut Center for Microbial Genomics, Department of Microbiology, Faculty of Science, Mahidol University, Rama 6 Road, Bangkok, Thailand.,Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA, USA
| | | | - Hideki Yanai
- Fukujuji Hospital and Research Institute of Tuberculosis, Japan Anti-Tuberculosis Association, Kiyose, 204-8533, Japan
| | - Htet Myat Win Maung
- National TB Control Programme, Department of Public Health, Ministry of Health and Sports, Naypyitaw, 15011, Myanmar.,Epidemiology Unit, Faculty of Medicine, Prince of Songkla University, Had Yai, 90110, Thailand
| | | | - Prasit Palittapongarnpim
- Pornchai Matangkasombut Center for Microbial Genomics, Department of Microbiology, Faculty of Science, Mahidol University, Rama 6 Road, Bangkok, Thailand. .,National Science and Technology Development Agency, Pathumthani, Thailand.
| |
Collapse
|
54
|
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
|
55
|
Deelder W, Napier G, Campino S, Palla L, Phelan J, Clark TG. A modified decision tree approach to improve the prediction and mutation discovery for drug resistance in Mycobacterium tuberculosis. BMC Genomics 2022; 23:46. [PMID: 35016609 PMCID: PMC8753810 DOI: 10.1186/s12864-022-08291-4] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2021] [Accepted: 01/03/2022] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Drug resistant Mycobacterium tuberculosis is complicating the effective treatment and control of tuberculosis disease (TB). With the adoption of whole genome sequencing as a diagnostic tool, machine learning approaches are being employed to predict M. tuberculosis resistance and identify underlying genetic mutations. However, machine learning approaches can overfit and fail to identify causal mutations if they are applied out of the box and not adapted to the disease-specific context. We introduce a machine learning approach that is customized to the TB setting, which extracts a library of genomic variants re-occurring across individual studies to improve genotypic profiling. RESULTS We developed a customized decision tree approach, called Treesist-TB, that performs TB drug resistance prediction by extracting and evaluating genomic variants across multiple studies. The application of Treesist-TB to rifampicin (RIF), isoniazid (INH) and ethambutol (EMB) drugs, for which resistance mutations are known, demonstrated a level of predictive accuracy similar to the widely used TB-Profiler tool (Treesist-TB vs. TB-Profiler tool: RIF 97.5% vs. 97.6%; INH 96.8% vs. 96.5%; EMB 96.8% vs. 95.8%). Application of Treesist-TB to less understood second-line drugs of interest, ethionamide (ETH), cycloserine (CYS) and para-aminosalisylic acid (PAS), led to the identification of new variants (52, 6 and 11, respectively), with a high number absent from the TB-Profiler library (45, 4, and 6, respectively). Thereby, Treesist-TB had improved predictive sensitivity (Treesist-TB vs. TB-Profiler tool: PAS 64.3% vs. 38.8%; CYS 45.3% vs. 30.7%; ETH 72.1% vs. 71.1%). CONCLUSION Our work reinforces the utility of machine learning for drug resistance prediction, while highlighting the need to customize approaches to the disease-specific context. Through applying a modified decision learning approach (Treesist-TB) across a range of anti-TB drugs, we identified plausible resistance-encoding genomic variants with high predictive ability, whilst potentially overcoming the overfitting challenges that can affect standard machine learning applications.
Collapse
Affiliation(s)
- Wouter Deelder
- London School of Hygiene & Tropical Medicine, Keppel Street, London, WC1E 7HT, UK
- Dalberg Advisors, 7 Rue de Chantepoulet, CH-1201, Geneva, Switzerland
| | - Gary Napier
- London School of Hygiene & Tropical Medicine, Keppel Street, London, WC1E 7HT, UK
| | - Susana Campino
- London School of Hygiene & Tropical Medicine, Keppel Street, London, WC1E 7HT, UK
| | - Luigi Palla
- London School of Hygiene & Tropical Medicine, Keppel Street, London, WC1E 7HT, UK
- Department of Public Health and Infectious Diseases, University of Rome La Sapienza, Rome, Italy
| | - Jody Phelan
- London School of Hygiene & Tropical Medicine, Keppel Street, London, WC1E 7HT, UK
| | - Taane G Clark
- London School of Hygiene & Tropical Medicine, Keppel Street, London, WC1E 7HT, UK.
- Department of Infection Biology, Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, London, UK.
| |
Collapse
|
56
|
Thawornwattana Y, Mahasirimongkol S, Yanai H, Maung HMW, Cui Z, Chongsuvivatwong V, Palittapongarnpim P. Revised nomenclature and SNP barcode for Mycobacterium tuberculosis lineage 2. Microb Genom 2021; 7. [PMID: 34787541 PMCID: PMC8743535 DOI: 10.1099/mgen.0.000697] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022] Open
Abstract
Mycobacterium tuberculosis (Mtb) lineage 2 (L2) strains are present globally, contributing to a widespread tuberculosis (TB) burden, particularly in Asia where both prevalence of TB and numbers of drug resistant TB are highest. The increasing availability of whole-genome sequencing (WGS) data worldwide provides an opportunity to improve our understanding of the global genetic diversity of Mtb L2 and its association with the disease epidemiology and pathogenesis. However, existing L2 sublineage classification schemes leave >20 % of the Modern Beijing isolates unclassified. Here, we present a revised SNP-based classification scheme of L2 in a genomic framework based on phylogenetic analysis of >4000 L2 isolates from 34 countries in Asia, Eastern Europe, Oceania and Africa. Our scheme consists of over 30 genotypes, many of which have not been described before. In particular, we propose six main genotypes of Modern Beijing strains, denoted L2.2.M1–L2.2.M6. We also provide SNP markers for genotyping L2 strains from WGS data. This fine-scale genotyping scheme, which can classify >98 % of the studied isolates, serves as a basis for more effective monitoring and reporting of transmission and outbreaks, as well as improving genotype-phenotype associations such as disease severity and drug resistance. This article contains data hosted by Microreact.
Collapse
Affiliation(s)
- Yuttapong Thawornwattana
- Pornchai Matangkasombut Center for Microbial Genomics, Department of Microbiology, Faculty of Science, Mahidol University, Bangkok 10400, Thailand.,Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA
| | | | - Hideki Yanai
- Fukujuji Hospital and Research Institute of Tuberculosis, Japan Anti-Tuberculosis Association, Kiyose 204-8533, Japan
| | - Htet Myat Win Maung
- National TB Control Programme, Department of Public Health, Ministry of Health and Sports, Naypyitaw 15011, Myanmar.,Epidemiology Unit, Faculty of Medicine, Prince of Songkla University, Had Yai 90110, Thailand
| | - Zhezhe Cui
- Epidemiology Unit, Faculty of Medicine, Prince of Songkla University, Had Yai 90110, Thailand.,Department of Tuberculosis Control, Guangxi Zhuang Autonomous Region Center for Disease Control and Prevention, Nanning, Guangxi, 530028, PR China
| | | | - Prasit Palittapongarnpim
- Pornchai Matangkasombut Center for Microbial Genomics, Department of Microbiology, Faculty of Science, Mahidol University, Bangkok 10400, Thailand.,National Science and Technology Development Agency, Pathumthani 12120, Thailand
| |
Collapse
|
57
|
Heupink TH, Verboven L, Warren RM, Van Rie A. Comprehensive and accurate genetic variant identification from contaminated and low-coverage Mycobacterium tuberculosis whole genome sequencing data. Microb Genom 2021; 7:000689. [PMID: 34793294 PMCID: PMC8743552 DOI: 10.1099/mgen.0.000689] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2021] [Accepted: 09/09/2021] [Indexed: 12/30/2022] Open
Abstract
Improved understanding of the genomic variants that allow Mycobacterium tuberculosis (Mtb ) to acquire drug resistance, or tolerance, and increase its virulence are important factors in controlling the current tuberculosis epidemic. Current approaches to Mtb sequencing, however, cannot reveal Mtb ’s full genomic diversity due to the strict requirements of low contamination levels, high Mtb sequence coverage and elimination of complex regions. We have developed the XBS (compleX Bacterial Samples) bioinformatics pipeline, which implements joint calling and machine-learning-based variant filtering tools to specifically improve variant detection in the important Mtb samples that do not meet these criteria, such as those from unbiased sputum samples. Using novel simulated datasets, which permit exact accuracy verification, XBS was compared to the UVP and MTBseq pipelines. Accuracy statistics showed that all three pipelines performed equally well for sequence data that resemble those obtained from culture isolates of high depth of coverage and low-level contamination. In the complex genomic regions, however, XBS accurately identified 9.0 % more SNPs and 8.1 % more single nucleotide insertions and deletions than the WHO-endorsed unified analysis variant pipeline. XBS also had superior accuracy for sequence data that resemble those obtained directly from sputum samples, where depth of coverage is typically very low and contamination levels are high. XBS was the only pipeline not affected by low depth of coverage (5–10×), type of contamination and excessive contamination levels (>50 %). Simulation results were confirmed using whole genome sequencing (WGS) data from clinical samples, confirming the superior performance of XBS with a higher sensitivity (98.8%) when analysing culture isolates and identification of 13.9 % more variable sites in WGS data from sputum samples as compared to MTBseq, without evidence for false positive variants when rRNA regions were excluded. The XBS pipeline facilitates sequencing of less-than-perfect Mtb samples. These advances will benefit future clinical applications of Mtb sequencing, especially WGS directly from clinical specimens, thereby avoiding in vitro biases and making many more samples available for drug resistance and other genomic analyses. The additional genetic resolution and increased sample success rate will improve genome-wide association studies and sequence-based transmission studies.
Collapse
Affiliation(s)
- Tim H. Heupink
- Family Medicine and Population Health (FAMPOP), Faculty of Medicine and Health Sciences, University of Antwerp, Antwerp, Belgium
| | - Lennert Verboven
- Family Medicine and Population Health (FAMPOP), Faculty of Medicine and Health Sciences, University of Antwerp, Antwerp, Belgium
| | - Robin M. Warren
- South African Medical Research Council Centre for Tuberculosis Research and DST/NRF Centre of Excellence for Biomedical Tuberculosis Research, Division of Molecular Biology and Human Genetics, Stellenbosch University, Stellenbosch, South Africa
| | - Annelies Van Rie
- Family Medicine and Population Health (FAMPOP), Faculty of Medicine and Health Sciences, University of Antwerp, Antwerp, Belgium
| |
Collapse
|
58
|
Practical approach to detection and surveillance of emerging highly resistant Mycobacterium tuberculosis Beijing 1071-32-cluster. Sci Rep 2021; 11:21392. [PMID: 34725411 PMCID: PMC8560753 DOI: 10.1038/s41598-021-00890-7] [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: 07/30/2021] [Accepted: 10/19/2021] [Indexed: 11/24/2022] Open
Abstract
Ancient sublineage of the Mycobacterium tuberculosis Beijing genotype is endemic and prevalent in East Asia and rare in other world regions. While these strains are mainly drug susceptible, we recently identified a novel clonal group Beijing 1071-32 within this sublineage emerging in Siberia, Russia and present in other Russian regions. This cluster included only multi/extensive drug resistant (MDR/XDR) isolates. Based on the phylogenetic analysis of the available WGS data, we identified three synonymous SNPs in the genes Rv0144, Rv0373c, and Rv0334 that were specific for the Beijing 1071-32-cluster and developed a real-time PCR assay for their detection. Analysis of the 2375 genetically diverse M. tuberculosis isolates collected between 1996 and 2020 in different locations (European and Asian parts of Russia, former Soviet Union countries, Albania, Greece, China, Vietnam, Japan and Brazil), confirmed 100% specificity and sensitivity of this real-time PCR assay. Moreover, the epidemiological importance of this strain and the newly developed screening assay is further stressed by the fact that all identified Beijing 1071-32 isolates were found to exhibit MDR genotypic profiles with concomitant resistance to additional first-line drugs due to a characteristic signature of six mutations in rpoB450, rpoC485, katG315, katG335, rpsL43 and embB497. In conclusion, this study provides a set of three concordant SNPs for the detection and screening of Beijing 1071-32 isolates along with a validated real-time PCR assay easily deployable across multiple settings for the epidemiological tracking of this significant MDR cluster.
Collapse
|
59
|
Meumann EM, Horan K, Ralph AP, Farmer B, Globan M, Stephenson E, Popple T, Boyd R, Kaestli M, Seemann T, Vandelannoote K, Lowbridge C, Baird RW, Stinear TP, Williamson DA, Currie BJ, Krause VL. Tuberculosis in Australia's tropical north: a population-based genomic epidemiological study. THE LANCET REGIONAL HEALTH. WESTERN PACIFIC 2021; 15:100229. [PMID: 34528010 PMCID: PMC8350059 DOI: 10.1016/j.lanwpc.2021.100229] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/27/2021] [Revised: 07/03/2021] [Accepted: 07/09/2021] [Indexed: 11/17/2022]
Abstract
BACKGROUND The Northern Territory (NT) has the highest tuberculosis (TB) rate of all Australian jurisdictions. We combined TB public health surveillance data with genomic sequencing of Mycobacterium tuberculosis isolates in the tropical 'Top End' of the NT to investigate trends in TB incidence and transmission. METHODS This retrospective observational study included all 741 culture-confirmed cases of TB in the Top End over three decades from 1989-2020. All 497 available M. tuberculosis isolates were sequenced. We used contact tracing data to define a threshold pairwise SNP distance for hierarchical single linkage clustering, and examined putative transmission clusters in the context of epidemiologic information. FINDINGS There were 359 (48%) cases born overseas, 329 (44%) cases among Australian First Nations peoples, and 52 (7%) cases were Australian-born and non-Indigenous. The annual incidence in First Nations peoples from 1989-2019 fell from average 50.4 to 11.0 per 100,000 (P<0·001). First Nations cases were more likely to die from TB (41/329, 12·5%) than overseas-born cases (11/359, 3·1%; P<0·001). Using a threshold of ≤12 SNPs, 28 clusters of between 2-64 individuals were identified, totalling 250 cases; 214 (86%) were First Nations cases and 189 (76%) were from a remote region. The time between cases and past epidemiologically- and genomically-linked contacts ranged from 4·5 months to 24 years. INTERPRETATION Our findings support prioritisation of timely case detection, contact tracing augmented by genomic sequencing, and latent TB treatment to break transmission chains in Top End remote hotspot regions.
Collapse
Affiliation(s)
- Ella M Meumann
- Global and Tropical Health Division, Menzies School of Health Research, Charles Darwin University, Darwin, Australia
- Department of Infectious Diseases, Division of Medicine, Royal Darwin Hospital, Darwin, Australia
- Territory Pathology, Royal Darwin Hospital, Darwin, Australia
- Nothern Territory Centre for Disease Control, Northern Territory Government, Darwin, Australia
| | - Kristy Horan
- Microbiological Diagnostic Unit Public Health Laboratory, Department of Microbiology and Immunology, The University of Melbourne at The Peter Doherty Institute for Infection and Immunity, Melbourne, Australia
| | - Anna P Ralph
- Global and Tropical Health Division, Menzies School of Health Research, Charles Darwin University, Darwin, Australia
- Department of Infectious Diseases, Division of Medicine, Royal Darwin Hospital, Darwin, Australia
- Nothern Territory Centre for Disease Control, Northern Territory Government, Darwin, Australia
| | - Belinda Farmer
- Nothern Territory Centre for Disease Control, Northern Territory Government, Darwin, Australia
| | - Maria Globan
- Mycobacterium Reference Laboratory, Victorian Infectious Diseases Reference Laboratory, Royal Melbourne Hospital at The Peter Doherty Institute for Infection and Immunity, Melbourne, Australia
| | - Elizabeth Stephenson
- Nothern Territory Centre for Disease Control, Northern Territory Government, Darwin, Australia
| | - Tracy Popple
- Nothern Territory Centre for Disease Control, Northern Territory Government, Darwin, Australia
| | - Rowena Boyd
- Nothern Territory Centre for Disease Control, Northern Territory Government, Darwin, Australia
| | - Mirjam Kaestli
- Global and Tropical Health Division, Menzies School of Health Research, Charles Darwin University, Darwin, Australia
| | - Torsten Seemann
- Microbiological Diagnostic Unit Public Health Laboratory, Department of Microbiology and Immunology, The University of Melbourne at The Peter Doherty Institute for Infection and Immunity, Melbourne, Australia
| | - Koen Vandelannoote
- Department of Microbiology and Immunology, The University of Melbourne at The Peter Doherty Institute for Infection and Immunity, Melbourne, Australia
| | - Christopher Lowbridge
- Global and Tropical Health Division, Menzies School of Health Research, Charles Darwin University, Darwin, Australia
- Nothern Territory Centre for Disease Control, Northern Territory Government, Darwin, Australia
| | - Robert W. Baird
- Department of Infectious Diseases, Division of Medicine, Royal Darwin Hospital, Darwin, Australia
- Territory Pathology, Royal Darwin Hospital, Darwin, Australia
| | - Timothy P. Stinear
- Department of Microbiology and Immunology, The University of Melbourne at The Peter Doherty Institute for Infection and Immunity, Melbourne, Australia
| | - Deborah A. Williamson
- Microbiological Diagnostic Unit Public Health Laboratory, Department of Microbiology and Immunology, The University of Melbourne at The Peter Doherty Institute for Infection and Immunity, Melbourne, Australia
| | - Bart J. Currie
- Global and Tropical Health Division, Menzies School of Health Research, Charles Darwin University, Darwin, Australia
- Department of Infectious Diseases, Division of Medicine, Royal Darwin Hospital, Darwin, Australia
| | - Vicki L. Krause
- Nothern Territory Centre for Disease Control, Northern Territory Government, Darwin, Australia
| |
Collapse
|
60
|
Genetic diversity of candidate loci linked to Mycobacterium tuberculosis resistance to bedaquiline, delamanid and pretomanid. Sci Rep 2021; 11:19431. [PMID: 34593898 PMCID: PMC8484543 DOI: 10.1038/s41598-021-98862-4] [Citation(s) in RCA: 34] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2021] [Accepted: 09/16/2021] [Indexed: 02/08/2023] Open
Abstract
Tuberculosis (TB), caused by Mycobacterium tuberculosis, is one of the deadliest infectious diseases worldwide. Multidrug and extensively drug-resistant strains are making disease control difficult, and exhausting treatment options. New anti-TB drugs bedaquiline (BDQ), delamanid (DLM) and pretomanid (PTM) have been approved for the treatment of multi-drug resistant TB, but there is increasing resistance to them. Nine genetic loci strongly linked to resistance have been identified (mmpR5, atpE, and pepQ for BDQ; ddn, fgd1, fbiA, fbiB, fbiC, and fbiD for DLM/PTM). Here we investigated the genetic diversity of these loci across >33,000 M. tuberculosis isolates. In addition, epistatic mutations in mmpL5-mmpS5 as well as variants in ndh, implicated for DLM/PTM resistance in M. smegmatis, were explored. Our analysis revealed 1,227 variants across the nine genes, with the majority (78%) present in isolates collected prior to the roll-out of BDQ and DLM/PTM. We identified phylogenetically-related mutations, which are unlikely to be resistance associated, but also high-impact variants such as frameshifts (e.g. in mmpR5, ddn) with likely functional effects, as well as non-synonymous mutations predominantly in MDR-/XDR-TB strains with predicted protein destabilising effects. Overall, our work provides a comprehensive mutational catalogue for BDQ and DLM/PTM associated genes, which will assist with establishing associations with phenotypic resistance; thereby, improving the understanding of the causative mechanisms of resistance for these drugs, leading to better treatment outcomes.
Collapse
|
61
|
Tunstall T, Phelan J, Eccleston C, Clark TG, Furnham N. Structural and Genomic Insights Into Pyrazinamide Resistance in Mycobacterium tuberculosis Underlie Differences Between Ancient and Modern Lineages. Front Mol Biosci 2021; 8:619403. [PMID: 34422898 PMCID: PMC8372558 DOI: 10.3389/fmolb.2021.619403] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2020] [Accepted: 04/14/2021] [Indexed: 11/30/2022] Open
Abstract
Resistance to drugs used to treat tuberculosis disease (TB) continues to remain a public health burden, with missense point mutations in the underlying Mycobacterium tuberculosis bacteria described for nearly all anti-TB drugs. The post-genomics era along with advances in computational and structural biology provide opportunities to understand the interrelationships between the genetic basis and the structural consequences of M. tuberculosis mutations linked to drug resistance. Pyrazinamide (PZA) is a crucial first line antibiotic currently used in TB treatment regimens. The mutational promiscuity exhibited by the pncA gene (target for PZA) necessitates computational approaches to investigate the genetic and structural basis for PZA resistance development. We analysed 424 missense point mutations linked to PZA resistance derived from ∼35K M. tuberculosis clinical isolates sourced globally, which comprised the four main M. tuberculosis lineages (Lineage 1-4). Mutations were annotated to reflect their association with PZA resistance. Genomic measures (minor allele frequency and odds ratio), structural features (surface area, residue depth and hydrophobicity) and biophysical effects (change in stability and ligand affinity) of point mutations on pncA protein stability and ligand affinity were assessed. Missense point mutations within pncA were distributed throughout the gene, with the majority (>80%) of mutations with a destabilising effect on protomer stability and on ligand affinity. Active site residues involved in PZA binding were associated with multiple point mutations highlighting mutational diversity due to selection pressures at these functionally important sites. There were weak associations between genomic measures and biophysical effect of mutations. However, mutations associated with PZA resistance showed statistically significant differences between structural features (surface area and residue depth), but not hydrophobicity score for mutational sites. Most interestingly M. tuberculosis lineage 1 (ancient lineage) exhibited a distinct protein stability profile for mutations associated with PZA resistance, compared to modern lineages.
Collapse
Affiliation(s)
- Tanushree Tunstall
- Department of Infection Biology, London School of Hygiene and Tropical Medicine, London, United Kingdom
| | - Jody Phelan
- Department of Infection Biology, London School of Hygiene and Tropical Medicine, London, United Kingdom
| | - Charlotte Eccleston
- Department of Infection Biology, London School of Hygiene and Tropical Medicine, London, United Kingdom
| | - Taane G. Clark
- Department of Infection Biology, London School of Hygiene and Tropical Medicine, London, United Kingdom
- Department of Infectious Disease Epidemiology, London School of Hygiene and Tropical Medicine, London, United Kingdom
| | - Nicholas Furnham
- Department of Infection Biology, London School of Hygiene and Tropical Medicine, London, United Kingdom
| |
Collapse
|
62
|
Characterization of rifampicin-resistant Mycobacterium tuberculosis in Khyber Pakhtunkhwa, Pakistan. Sci Rep 2021; 11:14194. [PMID: 34244539 PMCID: PMC8270973 DOI: 10.1038/s41598-021-93501-4] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2021] [Accepted: 06/24/2021] [Indexed: 11/08/2022] Open
Abstract
Tuberculosis (TB), caused by Mycobacterium tuberculosis, is endemic in Pakistan. Resistance to both firstline rifampicin and isoniazid drugs (multidrug-resistant TB; MDR-TB) is hampering disease control. Rifampicin resistance is attributed to rpoB gene mutations, but rpoA and rpoC loci may also be involved. To characterise underlying rifampicin resistance mutations in the TB endemic province of Khyber Pakhtunkhwa, we sequenced 51 M. tuberculosis isolates collected between 2016 and 2019; predominantly, MDR-TB (n = 44; 86.3%) and lineage 3 (n = 30, 58.8%) strains. We found that known mutations in rpoB (e.g. S405L), katG (e.g. S315T), or inhA promoter loci explain the MDR-TB. There were 24 unique mutations in rpoA, rpoB, and rpoC genes, including four previously unreported. Five instances of within-host resistance diversity were observed, where two were a mixture of MDR-TB strains containing mutations in rpoB, katG, and the inhA promoter region, as well as compensatory mutations in rpoC. Heteroresistance was observed in two isolates with a single lineage. Such complexity may reflect the high transmission nature of the Khyber Pakhtunkhwa setting. Our study reinforces the need to apply sequencing approaches to capture the full-extent of MDR-TB genetic diversity, to understand transmission, and to inform TB control activities in the highly endemic setting of Pakistan.
Collapse
|
63
|
Le Hang NT, Hijikata M, Maeda S, Miyabayashi A, Wakabayashi K, Seto S, Diem NTK, Yen NTT, Van Duc L, Thuong PH, Van Huan H, Hoang NP, Mitarai S, Keicho N, Kato S. Phenotypic and genotypic features of the Mycobacterium tuberculosis lineage 1 subgroup in central Vietnam. Sci Rep 2021; 11:13609. [PMID: 34193941 PMCID: PMC8245516 DOI: 10.1038/s41598-021-92984-5] [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: 01/04/2021] [Accepted: 06/15/2021] [Indexed: 11/09/2022] Open
Abstract
Mycobacterium tuberculosis (Mtb) has different features depending on different geographic areas. We collected Mtb strains from patients with smear-positive pulmonary tuberculosis in Da Nang, central Vietnam. Using a whole genome sequencing platform, including genome assembly complemented by long-read-sequencing data, genomic characteristics were studied. Of 181 Mtb isolates, predominant Vietnamese EAI4_VNM and EAI4-like spoligotypes (31.5%), ZERO strains (5.0%), and part of EAI5 (11.1%) were included in a lineage-1 (L1) sublineage, i.e., L1.1.1.1. These strains were found less often in younger people, and they genetically clustered less frequently than other modern strains. Patients infected with ZERO strains demonstrated less lung infiltration. A region in RD2bcg spanning six loci, i.e., PE_PGRS35, cfp21, Rv1985c, Rv1986, Rv1987, and erm(37), was deleted in EAI4_VNM, EAI4-like, and ZERO strains, whereas another 118 bp deletion in furA was specific only to ZERO strains. L1.1.1.1-sublineage-specific deletions in PE_PGRS4 and PE_PGRS22 were also identified. RD900, seen in ancestral lineages, was present in majority of the L1 members. All strains without IS6110 (5.0%) had the ZERO spoligo-pattern. Distinctive features of the ancestral L1 strains provide a basis for investigation of the modern versus ancestral Mtb lineages and allow consideration of countermeasures against this heterogeneous pathogen.
Collapse
Affiliation(s)
| | - Minako Hijikata
- Department of Pathophysiology and Host Defense, The Research Institute of Tuberculosis, JATA, Tokyo, Japan
| | - Shinji Maeda
- Faculty of Pharmaceutical Sciences, Hokkaido University of Science, Hokkaido, Japan
| | - Akiko Miyabayashi
- Department of Pathophysiology and Host Defense, The Research Institute of Tuberculosis, JATA, Tokyo, Japan
| | - Keiko Wakabayashi
- Department of Pathophysiology and Host Defense, The Research Institute of Tuberculosis, JATA, Tokyo, Japan
| | - Shintaro Seto
- Department of Pathophysiology and Host Defense, The Research Institute of Tuberculosis, JATA, Tokyo, Japan
| | | | | | - Le Van Duc
- Da Nang General Hospital, Da Nang, Vietnam
| | | | | | | | - Satoshi Mitarai
- Department of Mycobacterium Reference and Research, The Research Institute of Tuberculosis, JATA, Tokyo, Japan
| | - Naoto Keicho
- The Research Institute of Tuberculosis, JATA, Japan Anti-Tuberculosis Association, 3-1-24 Matsuyama, Kiyose, Tokyo, 204-8533, Japan. .,National Center for Global Health and Medicine, Tokyo, Japan.
| | - Seiya Kato
- The Research Institute of Tuberculosis, JATA, Japan Anti-Tuberculosis Association, 3-1-24 Matsuyama, Kiyose, Tokyo, 204-8533, Japan
| |
Collapse
|
64
|
Dubé JY, Fava VM, Schurr E, Behr MA. Underwhelming or Misunderstood? Genetic Variability of Pattern Recognition Receptors in Immune Responses and Resistance to Mycobacterium tuberculosis. Front Immunol 2021; 12:714808. [PMID: 34276708 PMCID: PMC8278570 DOI: 10.3389/fimmu.2021.714808] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2021] [Accepted: 06/17/2021] [Indexed: 12/23/2022] Open
Abstract
Human genetic control is thought to affect a considerable part of the outcome of infection with Mycobacterium tuberculosis (Mtb). Most of us deal with the pathogen by containment (associated with clinical "latency") or sterilization, but tragically millions each year do not. After decades of studies on host genetic susceptibility to Mtb infection, genetic variation has been discovered to play a role in tuberculous immunoreactivity and tuberculosis (TB) disease. Genes encoding pattern recognition receptors (PRRs) enable a consistent, molecularly direct interaction between humans and Mtb which suggests the potential for co-evolution. In this review, we explore the roles ascribed to PRRs during Mtb infection and ask whether such a longstanding and intimate interface between our immune system and this pathogen plays a critical role in determining the outcome of Mtb infection. The scientific evidence to date suggests that PRR variation is clearly implicated in altered immunity to Mtb but has a more subtle role in limiting the pathogen and pathogenesis. In contrast to 'effectors' like IFN-γ, IL-12, Nitric Oxide and TNF that are critical for Mtb control, 'sensors' like PRRs are less critical for the outcome of Mtb infection. This is potentially due to redundancy of the numerous PRRs in the innate arsenal, such that Mtb rarely goes unnoticed. Genetic association studies investigating PRRs during Mtb infection should therefore be designed to investigate endophenotypes of infection - such as immunological or clinical variation - rather than just TB disease, if we hope to understand the molecular interface between innate immunity and Mtb.
Collapse
Affiliation(s)
- Jean-Yves Dubé
- Department of Microbiology and Immunology, McGill University, Montreal, QC, Canada
- Program in Infectious Diseases and Immunity in Global Health, The Research Institute of the McGill University Health Centre, Montreal, QC, Canada
- McGill International TB Centre, McGill University, Montreal, QC, Canada
| | - Vinicius M. Fava
- Program in Infectious Diseases and Immunity in Global Health, The Research Institute of the McGill University Health Centre, Montreal, QC, Canada
- McGill International TB Centre, McGill University, Montreal, QC, Canada
| | - Erwin Schurr
- Department of Microbiology and Immunology, McGill University, Montreal, QC, Canada
- Program in Infectious Diseases and Immunity in Global Health, The Research Institute of the McGill University Health Centre, Montreal, QC, Canada
- McGill International TB Centre, McGill University, Montreal, QC, Canada
- Department of Human Genetics, Faculty of Medicine, McGill University, Montreal, QC, Canada
- Department of Medicine, Faculty of Medicine, McGill University, Montreal, QC, Canada
| | - Marcel A. Behr
- Department of Microbiology and Immunology, McGill University, Montreal, QC, Canada
- Program in Infectious Diseases and Immunity in Global Health, The Research Institute of the McGill University Health Centre, Montreal, QC, Canada
- McGill International TB Centre, McGill University, Montreal, QC, Canada
- Department of Medicine, Faculty of Medicine, McGill University, Montreal, QC, Canada
| |
Collapse
|
65
|
Khan MT, Ali S, Khan AS, Ali A, Khan A, Kaushik AC, Irfan M, Chinnasamy S, Zhang S, Zhang YJ, Cui Z, Wei AJ, Wang Y, Zhao M, Liu K, Wang H, Zeb MT, Wei DQ. Insight into the drug resistance whole genome of Mycobacterium tuberculosis isolates from Khyber Pakhtunkhwa, Pakistan. INFECTION GENETICS AND EVOLUTION 2021; 92:104861. [PMID: 33862292 DOI: 10.1016/j.meegid.2021.104861] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/27/2020] [Revised: 03/31/2021] [Accepted: 04/08/2021] [Indexed: 12/17/2022]
Abstract
Whole genome sequencing (WGS) is one of the most reliable methods for detection of drug resistance, genetic diversity in other virulence factor and also evolutionary dynamics of Mycobacterium tuberculosis complex (MTBC). First-line anti-tuberculosis drugs are the major weapons against Mycobacterium tuberculosis (MTB). However, the emergence of drug resistance remained a major obstacle towards global tuberculosis (TB) control program 2030, especially in high burden countries including Pakistan. To overcome the resistance and design potent drugs, genomic variations in drugs targets as well as in the virulence and evolutionary factors might be useful for better understanding and designing potential inhibitors. Here we aimed to find genomic variations in the first-line drugs targets, along with other virulence and evolutionary factors among the circulating isolates in Khyber Pakhtunkhwa, Pakistan. Samples were collected and drug susceptibility testing (DST) was performed as per WHO standard. The resistance samples were subjected to WGS. Among the five whole genome sequences, three samples (NCBI BioProject Accession: PRJNA629298, PRJNA629388) harbored 1997, 1162, and 2053 mutations. Some novel mutations have been detected in drugs targets. Similarly, numerous novel variants have also been detected in virulency and evolutionary factors, PE, PPE, and secretory system of MTB isolates. Exploring the genomic variations among the circulating isolates in geographical specific locations might be useful for future drug designing. To the best of our knowledge, this is the first study that provides useful data regarding the insight genomic variations in virulency, evolutionary factors including ESX and PE/PPE as well as drug targets, for better understanding and management of TB in a WHO declared high burden country.
Collapse
Affiliation(s)
- Muhammad Tahir Khan
- Institute of Molecular Biology and Biotechnology, The University of Lahore, Lahore, Pakistan
| | - Sajid Ali
- Department of Microbiology, Quaid-i-Azam University, Islamabad, Pakistan.
| | - Anwar Sheed Khan
- Department of Microbiology, Kohat University of Science and Technology and Provincial Tuberculosis Reference Laboratory, Peshawar, Pakistan.
| | - Arif Ali
- State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, and Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education, Shanghai Jiao Tong University, 800 Dongchuan Road, Minhang District, Shanghai 200240, China; Peng Cheng Laboratory, Vanke Cloud City Phase I Building 8, Xili Street, Nanshan District, Shenzhen, Guangdong 518055, China.
| | - Abbas Khan
- State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, and Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education, Shanghai Jiao Tong University, 800 Dongchuan Road, Minhang District, Shanghai 200240, China; Peng Cheng Laboratory, Vanke Cloud City Phase I Building 8, Xili Street, Nanshan District, Shenzhen, Guangdong 518055, China.
| | | | - Muhammad Irfan
- Department of Oral Biology, College of Dentistry, University of Florida, Gainesville, FL, USA.
| | - Sathishkumar Chinnasamy
- State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, and Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education, Shanghai Jiao Tong University, 800 Dongchuan Road, Minhang District, Shanghai 200240, China; Peng Cheng Laboratory, Vanke Cloud City Phase I Building 8, Xili Street, Nanshan District, Shenzhen, Guangdong 518055, China
| | - Shulin Zhang
- Department of Immunology and Microbiology, Shanghai Jiao Tong University, School of Medicine, Shanghai, China.
| | - Yu-Juan Zhang
- College of Life Sciences, Chongqing Normal University, Chongqing, China.
| | - Zhilei Cui
- Zhilei Cui, Department of Respiratory Medicine, XinHua Hospital Affiliated to Shanghai Jiao Tong University School of Medicine, Shanghai, China
| | - Amie Jinghua Wei
- State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, and Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education, Shanghai Jiao Tong University, 800 Dongchuan Road, Minhang District, Shanghai 200240, China; Peng Cheng Laboratory, Vanke Cloud City Phase I Building 8, Xili Street, Nanshan District, Shenzhen, Guangdong 518055, China
| | - Yanjie Wang
- State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, and Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education, Shanghai Jiao Tong University, 800 Dongchuan Road, Minhang District, Shanghai 200240, China; Peng Cheng Laboratory, Vanke Cloud City Phase I Building 8, Xili Street, Nanshan District, Shenzhen, Guangdong 518055, China
| | - Mingzhu Zhao
- State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, and Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education, Shanghai Jiao Tong University, 800 Dongchuan Road, Minhang District, Shanghai 200240, China; Peng Cheng Laboratory, Vanke Cloud City Phase I Building 8, Xili Street, Nanshan District, Shenzhen, Guangdong 518055, China.
| | - Kejia Liu
- State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, and Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education, Shanghai Jiao Tong University, 800 Dongchuan Road, Minhang District, Shanghai 200240, China; Peng Cheng Laboratory, Vanke Cloud City Phase I Building 8, Xili Street, Nanshan District, Shenzhen, Guangdong 518055, China.
| | - Heng Wang
- State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, and Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education, Shanghai Jiao Tong University, 800 Dongchuan Road, Minhang District, Shanghai 200240, China; Peng Cheng Laboratory, Vanke Cloud City Phase I Building 8, Xili Street, Nanshan District, Shenzhen, Guangdong 518055, China.
| | - Muhammad Tariq Zeb
- Khyber Medical University and Senior Research Officer, In-charge Genomic Laboratory, Veterinary Research Institute, Peshawar 25000, Pakistan
| | - Dong Qing Wei
- State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, and Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education, Shanghai Jiao Tong University, 800 Dongchuan Road, Minhang District, Shanghai 200240, China; Peng Cheng Laboratory, Vanke Cloud City Phase I Building 8, Xili Street, Nanshan District, Shenzhen, Guangdong 518055, China.
| |
Collapse
|