1
|
Khdhiri M, Thomas E, de Smet C, Chandar P, Chandrakumar I, Davidson JM, Anderson P, Chorlton SD. refMLST: reference-based multilocus sequence typing enables universal bacterial typing. BMC Bioinformatics 2024; 25:280. [PMID: 39192191 DOI: 10.1186/s12859-024-05913-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2023] [Accepted: 08/22/2024] [Indexed: 08/29/2024] Open
Abstract
BACKGROUND Commonly used approaches for genomic investigation of bacterial outbreaks, including SNP and gene-by-gene approaches, are limited by the requirement for background genomes and curated allele schemes, respectively. As a result, they only work on a select subset of known organisms, and fail on novel or less studied pathogens. We introduce refMLST, a gene-by-gene approach using the reference genome of a bacterium to form a scalable, reproducible and robust method to perform outbreak investigation. RESULTS When applied to multiple outbreak causing bacteria including 1263 Salmonella enterica, 331 Yersinia enterocolitica and 6526 Campylobacter jejuni genomes, refMLST enabled consistent clustering, improved resolution, and faster processing in comparison to commonly used tools like chewieSnake. CONCLUSIONS refMLST is a novel multilocus sequence typing approach that is applicable to any bacterial species with a public reference genome, does not require a curated scheme, and automatically accounts for genetic recombination. AVAILABILITY AND IMPLEMENTATION refMLST is freely available for academic use at https://bugseq.com/academic .
Collapse
Affiliation(s)
| | - Ella Thomas
- California Polytechnic State University, San Luis Obispo, CA, USA
| | - Chanel de Smet
- California Polytechnic State University, San Luis Obispo, CA, USA
| | - Priyanka Chandar
- California Polytechnic State University, San Luis Obispo, CA, USA
| | | | - Jean M Davidson
- California Polytechnic State University, San Luis Obispo, CA, USA
| | - Paul Anderson
- California Polytechnic State University, San Luis Obispo, CA, USA
| | | |
Collapse
|
2
|
Quan Z, Li M, Chen Y, Liang J, Takiff H, Gao Q. Performance evaluation of core genome multilocus sequence typing for genotyping of Mycobacterium tuberculosis strains in China: based on multicenter, population-based collection. Eur J Clin Microbiol Infect Dis 2024; 43:297-304. [PMID: 38041721 DOI: 10.1007/s10096-023-04720-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2023] [Accepted: 11/16/2023] [Indexed: 12/03/2023]
Abstract
PURPOSE To evaluate the performance of core genome multilocus sequence typing (cgMLST) for genotyping Mycobacterium tuberculosis (M.tuberculosis) Strains in regions where the lineage 2 strains predominate. METHODS We compared clustering by whole-genome SNP typing with cgMLST clustering in the analysis of WGS data of 6240 strains from five regions of China. Using both the receiver operating characteristic (ROC) curve and epidemiological investigation to determine the optimal threshold for defining genomic clustering by cgMLST. The performance of cgMLST was evaluated by quantifying the sensitivity, specificity and concordance of clustering between two methods. Logistic regression was used to gauge the impact of strain genetic diversity and lineage on cgMLST clustering. RESULTS The optimal threshold for cgMLST to define genomic clustering was determined to be ≤ 10 allelic differences between strains. The overall sensitivity and specificity of cgMLST averaged 99.6% and 96.3%, respectively; the concordance of clustering between two methods averaged 97.1%. Concordance was significantly correlated with strain genetic diversity and was 3.99 times (95% CI, 2.94-5.42) higher in regions with high genetic diversity (π > 1.55 × 10-4) compared to regions with low genetic diversity. The difference missed statistical significance, while concordance for lineage 2 strains (96.8%) was less than that for lineage 4 strains (98.3%). CONCLUSION : cgMLST showed a discriminatory power comparable to whole-genome SNP typing and could be used to genotype clinical M.tuberculosis strains in different regions of China. The discriminative power of cgMLST was significantly correlated with strain genetic diversity and was slightly lower with strains from regions with low genetic diversity.
Collapse
Affiliation(s)
- Zhuo Quan
- Shanghai Institute of Infectious Disease and Biosecurity, Key Laboratory of Medical Molecular Virology (MOE/NHC/CAMS), School of Basic Medical Science, Fudan University, 131 Dongan Road, Shanghai, 200032, China
| | - Meng Li
- Shanghai Institute of Infectious Disease and Biosecurity, Key Laboratory of Medical Molecular Virology (MOE/NHC/CAMS), School of Basic Medical Science, Fudan University, 131 Dongan Road, Shanghai, 200032, China
| | - Yiwang Chen
- Shanghai Institute of Infectious Disease and Biosecurity, Key Laboratory of Medical Molecular Virology (MOE/NHC/CAMS), School of Basic Medical Science, Fudan University, 131 Dongan Road, Shanghai, 200032, China
| | - Jialei Liang
- Shanghai Institute of Infectious Disease and Biosecurity, Key Laboratory of Medical Molecular Virology (MOE/NHC/CAMS), School of Basic Medical Science, Fudan University, 131 Dongan Road, Shanghai, 200032, China
| | - Howard Takiff
- Laboratorio de Genética Molecular, CMBC, Instituto Venezolano de Investigaciones Científicas, IVIC, Caracas, Venezuela
| | - Qian Gao
- Shanghai Institute of Infectious Disease and Biosecurity, Key Laboratory of Medical Molecular Virology (MOE/NHC/CAMS), School of Basic Medical Science, Fudan University, 131 Dongan Road, Shanghai, 200032, China.
| |
Collapse
|
3
|
Song Z, He W, Cao X, Ma A, He P, Zhao B, Wang S, Liu C, Zhao Y. The Recent Transmission and Associated Risk Factor of Mycobacterium tuberculosis in Golmud City, China. Infect Drug Resist 2024; 17:417-425. [PMID: 38318210 PMCID: PMC10840525 DOI: 10.2147/idr.s437026] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2023] [Accepted: 12/05/2023] [Indexed: 02/07/2024] Open
Abstract
Background Tuberculosis (TB) remains a severe public health problem globally, and it is essential to comprehend the transmission pattern to control tuberculosis. Herein, we evaluated the drug-resistant characteristics, recent transmission, and associated risk factors of TB in Golmud, Qinghai, China. Methods In this study, we performed a population-based study of patients diagnosed with TB in Golmud from 2013 to 2018. Drug-susceptibility testing and whole-genome sequencing were performed on 133 Mycobacterium tuberculosis strains. The genomic clustering rate was calculated to evaluate the level of recent transmission. Risk factors were identified by logistic regression analysis. Results Our results showed that 46.97% (62/132) of strains were phylogenetically clustered and formed into 23 transmission clusters, suggesting a high recent transmission of TB in the area. 12.78% (17/133) strains were multidrug-resistant/rifampicin tuberculosis (MDR/RR-TB), with a high drug-resistant burden. Based on drug resistance gene analysis, we found 23 strains belonging to genotype MDR/RR-TB, where some strains may have borderline mutations. Among these strains, 65.2% (15/23) were found within putative transmission clusters. Additionally, risk factor analysis showed that recent transmission of TB happened more in patients with Tibetan nationality or older age. Conclusion Overall our study indicates that the recent transmissions of MTB strains, especially genotypic MDR/RR strains, drive the tuberculosis epidemic in Golmud, which could contribute to developing effective TB prevention and control strategies.
Collapse
Affiliation(s)
- Zexuan Song
- National Key Laboratory of Intelligent Tracking and Forecasting for Infectious Diseases, National Center for Tuberculosis Control and Prevention, Chinese Center for Disease Control and Prevention, Changping, Beijing, 102206, People’s Republic of China
| | - Wencong He
- National Key Laboratory of Intelligent Tracking and Forecasting for Infectious Diseases, National Center for Tuberculosis Control and Prevention, Chinese Center for Disease Control and Prevention, Changping, Beijing, 102206, People’s Republic of China
| | - Xiaolong Cao
- National Key Laboratory of Intelligent Tracking and Forecasting for Infectious Diseases, National Center for Tuberculosis Control and Prevention, Chinese Center for Disease Control and Prevention, Changping, Beijing, 102206, People’s Republic of China
| | - Aijing Ma
- National Key Laboratory of Intelligent Tracking and Forecasting for Infectious Diseases, National Center for Tuberculosis Control and Prevention, Chinese Center for Disease Control and Prevention, Changping, Beijing, 102206, People’s Republic of China
| | - Ping He
- National Key Laboratory of Intelligent Tracking and Forecasting for Infectious Diseases, National Center for Tuberculosis Control and Prevention, Chinese Center for Disease Control and Prevention, Changping, Beijing, 102206, People’s Republic of China
| | - Bing Zhao
- National Key Laboratory of Intelligent Tracking and Forecasting for Infectious Diseases, National Center for Tuberculosis Control and Prevention, Chinese Center for Disease Control and Prevention, Changping, Beijing, 102206, People’s Republic of China
| | - Shengfen Wang
- National Key Laboratory of Intelligent Tracking and Forecasting for Infectious Diseases, National Center for Tuberculosis Control and Prevention, Chinese Center for Disease Control and Prevention, Changping, Beijing, 102206, People’s Republic of China
| | - Chunfa Liu
- National Key Laboratory of Intelligent Tracking and Forecasting for Infectious Diseases, National Center for Tuberculosis Control and Prevention, Chinese Center for Disease Control and Prevention, Changping, Beijing, 102206, People’s Republic of China
- Animal Science and Technology College, Beijing University of Agriculture, Huilongguan, Changping, Beijing, 102206, People’s Republic of China
| | - Yanlin Zhao
- National Key Laboratory of Intelligent Tracking and Forecasting for Infectious Diseases, National Center for Tuberculosis Control and Prevention, Chinese Center for Disease Control and Prevention, Changping, Beijing, 102206, People’s Republic of China
| |
Collapse
|
4
|
Wang TT, Hu YL, Li YF, Kong XL, Li YM, Sun PY, Wang DX, Li YY, Zhang YZ, Han QL, Zhu XH, An QQ, Liu LL, Liu Y, Li HC. Polyketide synthases mutation in tuberculosis transmission revealed by whole genomic sequence, China, 2011-2019. Front Genet 2024; 14:1217255. [PMID: 38259610 PMCID: PMC10800454 DOI: 10.3389/fgene.2023.1217255] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2023] [Accepted: 11/30/2023] [Indexed: 01/24/2024] Open
Abstract
Introduction: Tuberculosis (TB) is an infectious disease caused by a bacterium called Mycobacterium tuberculosis (Mtb). Previous studies have primarily focused on the transmissibility of multidrug-resistant (MDR) or extensively drug-resistant (XDR) Mtb. However, variations in virulence across Mtb lineages may also account for differences in transmissibility. In Mtb, polyketide synthase (PKS) genes encode large multifunctional proteins which have been shown to be major mycobacterial virulence factors. Therefore, this study aimed to identify the role of PKS mutations in TB transmission and assess its risk and characteristics. Methods: Whole genome sequences (WGSs) data from 3,204 Mtb isolates was collected from 2011 to 2019 in China. Whole genome single nucleotide polymorphism (SNP) profiles were used for phylogenetic tree analysis. Putative transmission clusters (≤10 SNPs) were identified. To identify the role of PKS mutations in TB transmission, we compared SNPs in the PKS gene region between "clustered isolates" and "non-clustered isolates" in different lineages. Results: Cluster-associated mutations in ppsA, pks12, and pks13 were identified among different lineage isolates. They were statistically significant among clustered strains, indicating that they may enhance the transmissibility of Mtb. Conclusion: Overall, this study provides new insights into the function of PKS and its localization in M. tuberculosis. The study found that ppsA, pks12, and pks13 may contribute to disease progression and higher transmission of certain strains. We also discussed the prospective use of mutant ppsA, pks12, and pks13 genes as drug targets.
Collapse
Affiliation(s)
- Ting-Ting Wang
- Shandong University of Traditional Chinese Medicine, Jinan, China
| | - Yuan-Long Hu
- Shandong University of Traditional Chinese Medicine, Jinan, China
| | - Yi-Fan Li
- Department of Pulmonary and Critical Care Medicine, The Third Affiliated Hospital of Shandong First Medical University (Affiliated Hospital of Shandong Academy of Medical Sciences), Jinan, China
| | - Xiang-Long Kong
- Shandong Artificial Intelligence Institute Qilu University of Technology (Shandong Academy of Sciences), Jinan, China
| | - Ya-Meng Li
- Shandong University of Traditional Chinese Medicine, Jinan, China
| | | | - Da-Xing Wang
- People’s Hospital of Huaiyin Jinan, Jinan, China
| | - Ying-Ying Li
- Shandong University of Traditional Chinese Medicine, Jinan, China
| | - Yu-Zhen Zhang
- Shandong First Medical University and Shandong Academy of Medical Sciences, Jinan, China
| | - Qi-Lin Han
- Shandong First Medical University and Shandong Academy of Medical Sciences, Jinan, China
| | - Xue-Han Zhu
- Shandong First Medical University and Shandong Academy of Medical Sciences, Jinan, China
| | - Qi-Qi An
- Department of Respiratory and Critical Care Medicine, Shandong Provincial Hospital Affiliated to 11 Shandong University, Shandong Provincial Hospital Affiliated to Shandong First Medical University, Jinan, Shandong, China
| | - Li-Li Liu
- People’s Hospital of Huaiyin Jinan, Jinan, China
| | - Yao Liu
- Department of Respiratory and Critical Care Medicine, Shandong Provincial Hospital Affiliated to 11 Shandong University, Shandong Provincial Hospital Affiliated to Shandong First Medical University, Jinan, Shandong, China
| | - Huai-Chen Li
- Shandong University of Traditional Chinese Medicine, Jinan, China
- Department of Respiratory and Critical Care Medicine, Shandong Provincial Hospital Affiliated to 11 Shandong University, Shandong Provincial Hospital Affiliated to Shandong First Medical University, Jinan, Shandong, China
| |
Collapse
|
5
|
Shibabaw A, Gelaw B, Ghanem M, Legall N, Schooley AM, Soehnlen MK, Salvador LCM, Gebreyes W, Wang SH, Tessema B. Molecular epidemiology and transmission dynamics of multi-drug resistant tuberculosis strains using whole genome sequencing in the Amhara region, Ethiopia. BMC Genomics 2023; 24:400. [PMID: 37460951 DOI: 10.1186/s12864-023-09502-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2023] [Accepted: 07/03/2023] [Indexed: 07/20/2023] Open
Abstract
BACKGROUND Drug resistant Mycobacterium tuberculosis prevention and care is a major challenge in Ethiopia. The World health organization has designated Ethiopia as one of the 30 high burden multi-drug resistant tuberculosis (MDR-TB) countries. There is limited information regarding genetic diversity and transmission dynamics of MDR-TB in Ethiopia. OBJECTIVE To investigate the molecular epidemiology and transmission dynamics of MDR-TB strains using whole genome sequence (WGS) in the Amhara region. METHODS Forty-five MDR-TB clinical isolates from Amhara region were collected between 2016 and 2018, and characterized using WGS and 24-loci Mycobacterium Interspersed Repetitive Units Variable Number of Tandem Repeats (MIRU-VNTR) typing. Clusters were defined based on the maximum distance of 12 single nucleotide polymorphisms (SNPs) or alleles as the upper threshold of genomic relatedness. Five or less SNPs or alleles distance or identical 24-loci VNTR typing is denoted as surrogate marker for recent transmission. RESULTS Forty-one of the 45 isolates were analyzed by WGS and 44% (18/41) of the isolates were distributed into 4 clusters. Of the 41 MDR-TB isolates, 58.5% were classified as lineage 4, 36.5% lineage 3 and 5% lineage 1. Overall, TUR genotype (54%) was the predominant in MDR-TB strains. 41% (17/41) of the isolates were clustered into four WGS groups and the remaining isolates were unique strains. The predominant cluster (Cluster 1) was composed of nine isolates belonging to lineage 4 and of these, four isolates were in the recent transmission links. CONCLUSIONS Majority of MDR-TB strain cluster and predominance of TUR lineage in the Amhara region give rise to concerns for possible ongoing transmission. Efforts to strengthen TB laboratory to advance diagnosis, intensified active case finding, and expanded contact tracing activities are needed in order to improve rapid diagnosis and initiate early treatment. This would lead to the interruption of the transmission chain and stop the spread of MDR-TB in the Amhara region.
Collapse
Affiliation(s)
- Agumas Shibabaw
- Department of Medical Laboratory Sciences, College of Medicine and Health Sciences, Wollo University, Dessie, Ethiopia.
- Global One Health Initiative (GOHi), The Ohio State University, Columbus, OH, USA.
- Department of Medical Microbiology, School of Medical Laboratory Sciences, College of Medicine and Health Sciences, University of Gondar, Gondar, Ethiopia.
- Michigan Department of Health and Human Services, Infectious disease, Lansing, MI, USA.
- Department of Veterinary Preventive Medicine, College of Veterinary Medicine, The Ohio State University, Columbus, OH, USA.
| | - Baye Gelaw
- Department of Medical Microbiology, School of Medical Laboratory Sciences, College of Medicine and Health Sciences, University of Gondar, Gondar, Ethiopia
| | - Mostafa Ghanem
- Department of Veterinary Medicine, Virginia-Maryland College of Veterinary Medicine, University of Maryland, College Park, MD, USA
| | - Noah Legall
- Institute of Bioinformatics, University of Georgia, Athens, GA, USA
| | - Angie M Schooley
- Michigan Department of Health and Human Services, Infectious disease, Lansing, MI, USA
| | - Marty K Soehnlen
- Michigan Department of Health and Human Services, Infectious disease, Lansing, MI, USA
| | - Liliana C M Salvador
- School of Animal and Comparative Biomedical Sciences, College of Agriculture and life sciences, University of Arizona, Tucson, AZ, USA
| | - Wondwossen Gebreyes
- Global One Health Initiative (GOHi), The Ohio State University, Columbus, OH, USA
- Department of Veterinary Preventive Medicine, College of Veterinary Medicine, The Ohio State University, Columbus, OH, USA
| | - Shu-Hua Wang
- Global One Health Initiative (GOHi), The Ohio State University, Columbus, OH, USA
- Department of Internal Medicine, Division of Infectious Diseases, College of Medicine, The Ohio State University, Columbus, OH, USA
| | - Belay Tessema
- Department of Medical Microbiology, School of Medical Laboratory Sciences, College of Medicine and Health Sciences, University of Gondar, Gondar, Ethiopia
| |
Collapse
|
6
|
Mekonnen D, Munshea A, Nibret E, Adnew B, Herrera-Leon S, Amor Aramendia A, Benito A, Abascal E, Jacqueline C, Aseffa A, Herrera-Leon L. Comparative whole-genome sequence analysis of Mycobacterium tuberculosis isolated from pulmonary tuberculosis and tuberculous lymphadenitis patients in Northwest Ethiopia. Front Microbiol 2023; 14:1211267. [PMID: 37455714 PMCID: PMC10348828 DOI: 10.3389/fmicb.2023.1211267] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2023] [Accepted: 05/30/2023] [Indexed: 07/18/2023] Open
Abstract
Background Tuberculosis (TB), caused by the Mycobacterium tuberculosis complex (MTBC), is a chronic infectious disease with both pulmonary and extrapulmonary forms. This study set out to investigate and compare the genomic diversity and transmission dynamics of Mycobacterium tuberculosis (Mtb) isolates obtained from tuberculous lymphadenitis (TBLN) and pulmonary TB (PTB) cases in Northwest Ethiopia. Methods A facility-based cross-sectional study was conducted using two groups of samples collected between February 2021 and June 2022 (Group 1) and between June 2020 and June 2022 (Group 2) in Northwest Ethiopia. Deoxyribonucleic acid (DNA) was extracted from 200 heat-inactivated Mtb isolates. Whole-genome sequencing (WGS) was performed from 161 isolates having ≥1 ng DNA/μl using Illumina NovaSeq 6000 technology. Results From the total 161 isolates sequenced, 146 Mtb isolates were successfully genotyped into three lineages (L) and 18 sub-lineages. The Euro-American (EA, L4) lineage was the prevailing (n = 100; 68.5%) followed by Central Asian (CAS, L3, n = 43; 25.3%) and then L7 (n = 3; 2.05%). The L4.2.2.ETH sub-lineage accounted for 19.9%, while Haarlem estimated at 13.7%. The phylogenetic tree revealed distinct Mtb clusters between PTB and TBLN isolates even though there was no difference at lineages and sub-lineages levels. The clustering rate (CR) and recent transmission index (RTI) for PTB were 30 and 15%, respectively. Similarly, the CR and RTI for TBLN were 31.1 and 18 %, respectively. Conclusion and recommendations PTB and TBLN isolates showed no Mtb lineages and sub-lineages difference. However, at the threshold of five allelic distances, Mtb isolates obtained from PTB and TBLN form distinct complexes in the phylogenetic tree, which indicates the presence of Mtb genomic variation among the two clinical forms. The high rate of clustering and RTI among TBLN implied that TBLN was likely the result of recent transmission and/or reactivation from short latency. Hence, the high incidence rate of TBLN in the Amhara region could be the result of Mtb genomic diversity and rapid clinical progression from primary infection and/or short latency. To validate this conclusion, a similar community-based study with a large sample size and better sampling technique is highly desirable. Additionally, analysis of genomic variants other than phylogenetic informative regions could give insightful information. Combined analysis of the host and the pathogen genome (GXG) together with environmental (GxGxE) factors could give comprehensive co-evolutionary information.
Collapse
Affiliation(s)
- Daniel Mekonnen
- Department of Medical Laboratory Sciences, School of Health Science, College of Medicine and Health Sciences, Bahir Dar University, Bahir Dar, Ethiopia
- Health Biotechnology Division, Institute of Biotechnology, Bahir Dar University, Bahir Dar, Ethiopia
- Amhara Public Health Institute, Bahir Dar, Ethiopia
| | - Abaineh Munshea
- Health Biotechnology Division, Institute of Biotechnology, Bahir Dar University, Bahir Dar, Ethiopia
- Department of Biology, Bahir Dar University, Bahir Dar, Ethiopia
| | - Endalkachew Nibret
- Health Biotechnology Division, Institute of Biotechnology, Bahir Dar University, Bahir Dar, Ethiopia
- Department of Biology, Bahir Dar University, Bahir Dar, Ethiopia
| | | | - Silvia Herrera-Leon
- National Centre for Microbiology, Instituto de Salud Carlos III, Madrid, Spain
| | | | - Agustín Benito
- National Center of Tropical Medicine, Institute of Health Carlos III, Centro de Investigación Biomédica en Red de Enfermedades Infecciosas, Madrid, Spain
| | - Estefanía Abascal
- National Centre for Microbiology, Instituto de Salud Carlos III, Madrid, Spain
| | - Camille Jacqueline
- National Centre for Microbiology, Instituto de Salud Carlos III, Madrid, Spain
- European Public Health Microbiology Training Programme, European Centre for Disease Prevention and Control, Stockholm, Sweden
| | - Abraham Aseffa
- Armauer Hansen Research Institute, Addis Ababa, Ethiopia
| | - Laura Herrera-Leon
- National Centre for Microbiology, Instituto de Salud Carlos III, Madrid, Spain
- CIBER Epidemiologia y Salud Publica, Madrid, Spain
| |
Collapse
|
7
|
Diel R, Meywald-Walter K, Schwarzbach C, Voss K, Dreyer V, Niemann S. Risk of tuberculosis transmission by children in Hamburg, Germany. Respir Med 2023; 209:107152. [PMID: 36781053 DOI: 10.1016/j.rmed.2023.107152] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/24/2022] [Revised: 02/06/2023] [Accepted: 02/08/2023] [Indexed: 02/13/2023]
Abstract
BACKGROUND Data from a prospective molecular-epidemiologic study (1997-2021) in Hamburg, Germany, were evaluated to assess the transmission risk of Mycobacterium tuberculosis complex (Mtbc) by children <15 years in a low-incidence setting. METHODS Isolates of Mtbc were genotyped whole genome sequencing, applying a core genome multilocus sequence typing scheme. Close contacts of culture-confirmed children were examined for latent Mtbc infections (LTBI) with particular focus on IGRA testing. RESULTS Out of 3154 culture-confirmed tuberculosis (TB) cases, 79 (2.5%) were children <15 years. Of those, 52 (58%) had pulmonary TB. Genotyping revealed that 35 of the 52 children (67%) were epidemiologically confirmed secondary cluster members; all of their source cases were adults. Six immigrant children presented without a presumed source case; their TB diagnoses came on average 48 weeks (interquartile range [IQR] 71) after their arrival in Germany. Three German-born children were determined to have been infected by adult relatives while visiting their parents' home country. Of the 317 children's close contacts tested with QuantiFERON-TB Gold-In Tube for LTBI, only 21 (6.6%) were positive. Absent a history of prior exposure or immigration from a high-incidence country, none of the contacts of younger (<10 years) TB-afflicted children was latently infected, whereas 2 older children infected 12 of their contacts, children and adults. During a mean observational period of 551 weeks (IQR 735) on average, no secondary TB cases appeared. CONCLUSIONS Children with pulmonary TB disease, especially those aged below 10 years, rarely transmit Mtbc to their close contacts in a low-incidence setting.
Collapse
Affiliation(s)
- Roland Diel
- Institute for Epidemiology, University Medical Hospital Schleswig-Holstein, Kiel, Germany; LungClinic Grosshansdorf, Airway Research Center North (ARCN), German Center for Lung Research (DZL), Germany.
| | | | | | - Klaas Voss
- Public Health Department Hamburg-Central, Hamburg, Germany
| | - Viola Dreyer
- Molecular and Experimental Mycobacteriology Group, National Reference Laboratory for Mycobacteria, Research Center Borstel, Borstel, Germany; German Center for Infection Research, Borstel, Germany
| | - Stefan Niemann
- Molecular and Experimental Mycobacteriology Group, National Reference Laboratory for Mycobacteria, Research Center Borstel, Borstel, Germany; German Center for Infection Research, Borstel, Germany
| |
Collapse
|
8
|
Transmission and Drug Resistance Genotype of Multidrug-Resistant or Rifampicin-Resistant Mycobacterium tuberculosis in Chongqing, China. Microbiol Spectr 2022; 10:e0240521. [PMID: 36214695 PMCID: PMC9604020 DOI: 10.1128/spectrum.02405-21] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Abstract
Multidrug-resistant or rifampicin-resistant tuberculosis (MDR/RR-TB) is a global barrier for the Stop TB plan. To identify risk factors for treatment outcome and cluster transmission of MDR/RR-TB, whole-genome sequencing (WGS) data of isolates from patients of the Chongqing Tuberculosis Control Institute were used for phylogenetic classifications, resistance predictions, and cluster analysis. A total of 223 MDR/RR-TB cases were recorded between 1 January 2018 and 31 December 2020. Elderly patients and those with lung cavitation are at increased risk of death due to MDR/RR-TB. A total of 187 MDR/RR strains were obtained from WGS data; 152 were classified as lineage 2 strains. Eighty (42.8%) strains differing by a distance of 12 or fewer single nucleotide polymorphisms were classified as 20 genomic clusters, indicating recent transmission. Patients infected with lineage 2 strains or those with occupations listed as "other" are significantly associated with a transmission cluster of MDR/RR-TB. Analysis of resistant mutations against first-line tuberculosis drugs found that 76 (95.0%) of all 80 strains had the same mutations within each cluster. A total of 55.0% (44 of 80) of the MDR/RR-TB strains accumulated additional drug resistance mutations along the transmission chain, especially against fluoroquinolones (63.6% [28 of 44]). Recent transmission of MDR/RR strains is driving the MDR/RR-TB epidemics, leading to the accumulation of more serious resistance along the transmission chains. IMPORTANCE The drug resistance molecular characteristics of MDR/RR-TB were elucidated by genome-wide analysis, and risk factors for death by MDR/RR-TB were identified in combination with patient information. Cluster characteristics of MDR/RR-TB in the region were analyzed by genome-wide analysis, and risk factors for cluster transmission (recent transmission) were analyzed. These analyses provide reference for the prevention and treatment of MDR/RR-TB in Chongqing.
Collapse
|
9
|
Delineating Mycobacterium abscessus population structure and transmission employing high-resolution core genome multilocus sequence typing. Nat Commun 2022; 13:4936. [PMID: 35999208 PMCID: PMC9399081 DOI: 10.1038/s41467-022-32122-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: 03/23/2022] [Accepted: 07/19/2022] [Indexed: 11/08/2022] Open
Abstract
Mycobacterium abscessus is an emerging multidrug-resistant non-tuberculous mycobacterium that causes a wide spectrum of infections and has caused several local outbreaks worldwide. To facilitate standardized prospective molecular surveillance, we established a novel core genome multilocus sequence typing (cgMLST) scheme. Whole genome sequencing data of 1991 isolates were employed to validate the scheme, re-analyze global population structure and set genetic distance thresholds for cluster detection and taxonomic identification. We confirmed and amended the nomenclature of the main dominant circulating clones and found that these also correlate well with traditional 7-loci MLST. Dominant circulating clones could be linked to a corresponding reference genome with less than 250 alleles while 99% of pairwise comparisons between epidemiologically linked isolates were below 25 alleles and 90% below 10 alleles. These thresholds can be used to guide further epidemiological investigations. Overall, the scheme will help to unravel the apparent global spread of certain clonal complexes and as yet undiscovered transmission routes.
Collapse
|
10
|
Dreyer V, Mandal A, Dev P, Merker M, Barilar I, Utpatel C, Nilgiriwala K, Rodrigues C, Crook DW, Crook DW, Peto TEA, Walker AS, Hoosdally SJ, Gibertoni Cruz AL, Carter J, Earle S, Kouchaki S, Yang Y, Walker TM, Fowler PW, Wilson D, Clifton DA, Iqbal Z, Hunt M, Knaggs J, Cirillo DM, Borroni E, Battaglia S, Ghodousi A, Spitaleri A, Cabibbe A, Tahseen S, Nilgiriwala K, Shah S, Rodrigues C, Kambli P, Surve U, Khot R, Niemann S, Kohl T, Merker M, Hoffmann H, Todt K, Plesnik S, Ismail N, Omar SV, Ngcamu LJD, Okozi N, Yao SY, Thwaites G, Thuong TNT, Ngoc NH, Srinivasan V, Moore D, Coronel J, Solano W, Gao GF, He G, Zhao Y, Ma A, Liu C, Zhu B, Laurenson I, Claxton P, Wilkinson RJ, Koch A, Lalvani A, Posey J, Gardy J, Werngren J, Paton N, Jou R, Wu MH, Xiao YX, Ferrazoli L, de Oliveira RS, Millard J, Warren R, Van Rie A, Lapierre SG, Rabodoarivelo MS, Rakotosamimanana N, Nimmo C, Musser K, Escuyer V, Cohen T, Rasigade JP, Wirth T, Mistry N, Niemann S. High fluoroquinolone resistance proportions among multidrug-resistant tuberculosis driven by dominant L2 Mycobacterium tuberculosis clones in the Mumbai Metropolitan Region. Genome Med 2022; 14:95. [PMID: 35989319 PMCID: PMC9394022 DOI: 10.1186/s13073-022-01076-0] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2022] [Accepted: 06/20/2022] [Indexed: 11/10/2022] Open
Abstract
Abstract
Background
Multidrug-resistant (MDR) Mycobacterium tuberculosis complex (MTBC) strains are a serious health problem in India, also contributing to one-fourth of the global MDR tuberculosis (TB) burden. About 36% of the MDR MTBC strains are reported fluoroquinolone (FQ) resistant leading to high pre-extensively drug-resistant (pre-XDR) and XDR-TB (further resistance against bedaquiline and/or linezolid) rates. Still, factors driving the MDR/pre-XDR epidemic in India are not well defined.
Methods
In a retrospective study, we analyzed 1852 consecutive MTBC strains obtained from patients from a tertiary care hospital laboratory in Mumbai by whole genome sequencing (WGS). Univariate and multivariate statistics was used to investigate factors associated with pre-XDR. Core genome multi locus sequence typing, time scaled haplotypic density (THD) method and homoplasy analysis were used to analyze epidemiological success, and positive selection in different strain groups, respectively.
Results
In total, 1016 MTBC strains were MDR, out of which 703 (69.2%) were pre-XDR and 45 (4.4%) were XDR. Cluster rates were high among MDR (57.8%) and pre-XDR/XDR (79%) strains with three dominant L2 (Beijing) strain clusters (Cl 1–3) representing half of the pre-XDR and 40% of the XDR-TB cases. L2 strains were associated with pre-XDR/XDR-TB (P < 0.001) and, particularly Cl 1–3 strains, had high first-line and FQ resistance rates (81.6–90.6%). Epidemic success analysis using THD showed that L2 strains outperformed L1, L3, and L4 strains in short- and long-term time scales. More importantly, L2 MDR and MDR + strains had higher THD success indices than their not-MDR counterparts. Overall, compensatory mutation rates were highest in L2 strains and positive selection was detected in genes of L2 strains associated with drug tolerance (prpB and ppsA) and virulence (Rv2828c). Compensatory mutations in L2 strains were associated with a threefold increase of THD indices, suggesting improved transmissibility.
Conclusions
Our data indicate a drastic increase of FQ resistance, as well as emerging bedaquiline resistance which endangers the success of newly endorsed MDR-TB treatment regimens. Rapid changes in treatment and control strategies are required to contain transmission of highly successful pre-XDR L2 strains in the Mumbai Metropolitan region but presumably also India-wide.
Collapse
|
11
|
Meumann EM, Krause VL, Baird R, Currie BJ. Using Genomics to Understand the Epidemiology of Infectious Diseases in the Northern Territory of Australia. Trop Med Infect Dis 2022; 7:tropicalmed7080181. [PMID: 36006273 PMCID: PMC9413455 DOI: 10.3390/tropicalmed7080181] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2022] [Revised: 08/09/2022] [Accepted: 08/11/2022] [Indexed: 11/16/2022] Open
Abstract
The Northern Territory (NT) is a geographically remote region of northern and central Australia. Approximately a third of the population are First Nations Australians, many of whom live in remote regions. Due to the physical environment and climate, and scale of social inequity, the rates of many infectious diseases are the highest nationally. Molecular typing and genomic sequencing in research and public health have provided considerable new knowledge on the epidemiology of infectious diseases in the NT. We review the applications of genomic sequencing technology for molecular typing, identification of transmission clusters, phylogenomics, antimicrobial resistance prediction, and pathogen detection. We provide examples where these methodologies have been applied to infectious diseases in the NT and discuss the next steps in public health implementation of this technology.
Collapse
Affiliation(s)
- Ella M. Meumann
- Global and Tropical Health Division, Menzies School of Health Research, Charles Darwin University, Darwin 0810, Australia
- Department of Infectious Diseases, Division of Medicine, Royal Darwin Hospital, Darwin 0810, Australia
- Correspondence:
| | - Vicki L. Krause
- Northern Territory Centre for Disease Control, Northern Territory Government, Darwin 0810, Australia
| | - Robert Baird
- Territory Pathology, Royal Darwin Hospital, Darwin 0810, Australia
| | - Bart J. Currie
- Global and Tropical Health Division, Menzies School of Health Research, Charles Darwin University, Darwin 0810, Australia
- Department of Infectious Diseases, Division of Medicine, Royal Darwin Hospital, Darwin 0810, Australia
| |
Collapse
|
12
|
Volk D, Yang-Turner F, Didelot X, Crook DW, Wyllie D. Catwalk: identifying closely related sequences in large microbial sequence databases. Microb Genom 2022; 8. [PMID: 35771206 PMCID: PMC9455716 DOI: 10.1099/mgen.0.000850] [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] [Indexed: 11/18/2022] Open
Abstract
There is a need to identify microbial sequences that may form part of transmission chains, or that may represent importations across national boundaries, amidst large numbers of SARS-CoV-2 and other bacterial or viral sequences. Reference-based compression is a sequence analysis technique that allows both a compact storage of sequence data and comparisons between sequences. Published implementations of the approach are being challenged by the large sample collections now being generated. Our aim was to develop a fast software detecting highly similar sequences in large collections of microbial genomes, including millions of SARS-CoV-2 genomes. To do so, we developed Catwalk, a tool that bypasses bottlenecks in the generation, comparison and in-memory storage of microbial genomes generated by reference mapping. It is a compiled solution, coded in Nim to increase performance. It can be accessed via command line, rest api or web server interfaces. We tested Catwalk using both SARS-CoV-2 and Mycobacterium tuberculosis genomes generated by prospective public-health sequencing programmes. Pairwise sequence comparisons, using clinically relevant similarity cut-offs, took about 0.39 and 0.66 μs, respectively; in 1 s, between 1 and 2 million sequences can be searched. Catwalk operates about 1700 times faster than, and uses about 8 % of the RAM of, a Python reference-based compression and comparison tool in current use for outbreak detection. Catwalk can rapidly identify close relatives of a SARS-CoV-2 or M. tuberculosis genome amidst millions of samples.
Collapse
Affiliation(s)
- Denis Volk
- Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | - Fan Yang-Turner
- Nuffield Department of Medicine, University of Oxford, Oxford, UK.,Present address: UKRI Science and Technologies Facilities Council, Harwell, UK
| | - Xavier Didelot
- School of Life Sciences, University of Warwick, Coventry CV4 7AL, UK.,Department of Statistics, University of Warwick, Coventry CV4 7AL, UK
| | - Derrick W Crook
- Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | - David Wyllie
- UK Health Security Agency, Forvie Site, Addenbrookes' Campus, Robinson Way, Cambridge CB2 0SR, UK
| |
Collapse
|
13
|
Mudliar SKR, Kulsum U, Rufai SB, Umpo M, Nyori M, Singh S. Snapshot of Mycobacterium tuberculosis Phylogenetics from an Indian State of Arunachal Pradesh Bordering China. Genes (Basel) 2022. [DOI: https://doi.org/10.3390/genes13020263] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Uncontrolled transmission of Mycobacterium tuberculosis (M. tuberculosis, MTB) drug resistant strains is a challenge to control efforts of the global tuberculosis program. Due to increasing multi-drug resistant (MDR) cases in Arunachal Pradesh, a northeastern state of India, the tracking and tracing of these resistant MTB strains is crucial for infection control and spread of drug resistance. This study aims to correlate the phenotypic DST, genomic DST (gDST) and phylogenetic analysis of MDR-MTB strains in the region. Of the total 200 samples 22 (11%) patients suspected of MDR-TB and 160 (80%) previously treated MDR-TB cases, 125 (62.5%) were identified as MTB. MGIT-960 SIRE DST detected 71/125 (56.8%) isolates as MDR/RR-MTB of which 22 (30.9%) were detected resistant to second-line drugs. Whole-genome sequencing of 65 isolates and their gDST found Ser315Thr mutation in katG (35/45; 77.8%) and Ser531Leu mutation in rpoB (21/41; 51.2%) associated with drug resistance. SNP barcoding categorized the dataset with Lineage2 (41; 63.1%) being predominant followed by Lineage3 (10; 15.4%), Lineage1 (8; 12.3%) and Lineage4 (6; 9.2%) respectively. Phylogenetic assignment by cgMLST gave insights of two Beijing sub-lineages viz; 2.2.1 (SNP difference < 19) and 2.2.1.2 (SNP difference < 9) associated with recent ongoing transmission in Arunachal Pradesh. This study provides insights in identifying two virulent Beijing sub-lineages (sub-lineage 2.2.1 and 2.2.1.2) with ongoing transmission of TB drug resistance in Arunachal Pradesh.
Collapse
|
14
|
Mudliar SKR, Kulsum U, Rufai SB, Umpo M, Nyori M, Singh S. Snapshot of Mycobacterium tuberculosis Phylogenetics from an Indian State of Arunachal Pradesh Bordering China. Genes (Basel) 2022. [DOI: https:/doi.org/10.3390/genes13020263] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Abstract
Uncontrolled transmission of Mycobacterium tuberculosis (M. tuberculosis, MTB) drug resistant strains is a challenge to control efforts of the global tuberculosis program. Due to increasing multi-drug resistant (MDR) cases in Arunachal Pradesh, a northeastern state of India, the tracking and tracing of these resistant MTB strains is crucial for infection control and spread of drug resistance. This study aims to correlate the phenotypic DST, genomic DST (gDST) and phylogenetic analysis of MDR-MTB strains in the region. Of the total 200 samples 22 (11%) patients suspected of MDR-TB and 160 (80%) previously treated MDR-TB cases, 125 (62.5%) were identified as MTB. MGIT-960 SIRE DST detected 71/125 (56.8%) isolates as MDR/RR-MTB of which 22 (30.9%) were detected resistant to second-line drugs. Whole-genome sequencing of 65 isolates and their gDST found Ser315Thr mutation in katG (35/45; 77.8%) and Ser531Leu mutation in rpoB (21/41; 51.2%) associated with drug resistance. SNP barcoding categorized the dataset with Lineage2 (41; 63.1%) being predominant followed by Lineage3 (10; 15.4%), Lineage1 (8; 12.3%) and Lineage4 (6; 9.2%) respectively. Phylogenetic assignment by cgMLST gave insights of two Beijing sub-lineages viz; 2.2.1 (SNP difference < 19) and 2.2.1.2 (SNP difference < 9) associated with recent ongoing transmission in Arunachal Pradesh. This study provides insights in identifying two virulent Beijing sub-lineages (sub-lineage 2.2.1 and 2.2.1.2) with ongoing transmission of TB drug resistance in Arunachal Pradesh.
Collapse
|
15
|
Mudliar SKR, Kulsum U, Rufai SB, Umpo M, Nyori M, Singh S. Snapshot of Mycobacterium tuberculosis Phylogenetics from an Indian State of Arunachal Pradesh Bordering China. Genes (Basel) 2022; 13:genes13020263. [PMID: 35205308 PMCID: PMC8872330 DOI: 10.3390/genes13020263] [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: 12/16/2021] [Revised: 01/21/2022] [Accepted: 01/25/2022] [Indexed: 02/07/2023] Open
Abstract
Uncontrolled transmission of Mycobacterium tuberculosis (M. tuberculosis, MTB) drug resistant strains is a challenge to control efforts of the global tuberculosis program. Due to increasing multi-drug resistant (MDR) cases in Arunachal Pradesh, a northeastern state of India, the tracking and tracing of these resistant MTB strains is crucial for infection control and spread of drug resistance. This study aims to correlate the phenotypic DST, genomic DST (gDST) and phylogenetic analysis of MDR-MTB strains in the region. Of the total 200 samples 22 (11%) patients suspected of MDR-TB and 160 (80%) previously treated MDR-TB cases, 125 (62.5%) were identified as MTB. MGIT-960 SIRE DST detected 71/125 (56.8%) isolates as MDR/RR-MTB of which 22 (30.9%) were detected resistant to second-line drugs. Whole-genome sequencing of 65 isolates and their gDST found Ser315Thr mutation in katG (35/45; 77.8%) and Ser531Leu mutation in rpoB (21/41; 51.2%) associated with drug resistance. SNP barcoding categorized the dataset with Lineage2 (41; 63.1%) being predominant followed by Lineage3 (10; 15.4%), Lineage1 (8; 12.3%) and Lineage4 (6; 9.2%) respectively. Phylogenetic assignment by cgMLST gave insights of two Beijing sub-lineages viz; 2.2.1 (SNP difference < 19) and 2.2.1.2 (SNP difference < 9) associated with recent ongoing transmission in Arunachal Pradesh. This study provides insights in identifying two virulent Beijing sub-lineages (sub-lineage 2.2.1 and 2.2.1.2) with ongoing transmission of TB drug resistance in Arunachal Pradesh.
Collapse
Affiliation(s)
- Shiv kumar Rashmi Mudliar
- Department of Microbiology, All India Institute of Medical Sciences, Bhopal 462020, Madhya Pradesh, India; (S.k.R.M.); (U.K.)
| | - Umay Kulsum
- Department of Microbiology, All India Institute of Medical Sciences, Bhopal 462020, Madhya Pradesh, India; (S.k.R.M.); (U.K.)
| | - Syed Beenish Rufai
- Infectious Diseases and Immunity in Global Health Program, Research Institute of the McGill University Health Center, Montreal, QC H4A 3J1, Canada;
- McGill International TB Center, Montreal, QC H4A 3J1, Canada
| | - Mika Umpo
- Tomo Riba Institute of Health & Medical Sciences, Naharlagun 791110, Arunachal Pradesh, India;
| | - Moi Nyori
- State TB Cell, Naharlagun 791110, Arunachal Pradesh, India;
| | - Sarman Singh
- Department of Microbiology, All India Institute of Medical Sciences, Bhopal 462020, Madhya Pradesh, India; (S.k.R.M.); (U.K.)
- Correspondence: ; Tel.: +91-9810813435
| |
Collapse
|
16
|
Peker N, Schuele L, Kok N, Terrazos M, Neuenschwander SM, de Beer J, Akkerman O, Peter S, Ramette A, Merker M, Niemann S, Couto N, Sinha B, Rossen JWA. Evaluation of whole-genome sequence data analysis approaches for short- and long-read sequencing of Mycobacterium tuberculosis. Microb Genom 2021; 7:000695. [PMID: 34825880 PMCID: PMC8743536 DOI: 10.1099/mgen.0.000695] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2021] [Accepted: 09/15/2021] [Indexed: 12/17/2022] Open
Abstract
Whole-genome sequencing (WGS) of Mycobacterium tuberculosis (MTB) isolates can be used to get an accurate diagnosis, to guide clinical decision making, to control tuberculosis (TB) and for outbreak investigations. We evaluated the performance of long-read (LR) and/or short-read (SR) sequencing for anti-TB drug-resistance prediction using the TBProfiler and Mykrobe tools, the fraction of genome recovery, assembly accuracies and the robustness of two typing approaches based on core-genome SNP (cgSNP) typing and core-genome multi-locus sequence typing (cgMLST). Most of the discrepancies between phenotypic drug-susceptibility testing (DST) and drug-resistance prediction were observed for the first-line drugs rifampicin, isoniazid, pyrazinamide and ethambutol, mainly with LR sequence data. Resistance prediction to second-line drugs made by both TBProfiler and Mykrobe tools with SR- and LR-sequence data were in complete agreement with phenotypic DST except for one isolate. The SR assemblies were more accurate than the LR assemblies, having significantly (P <0.05) fewer indels and mismatches per 100 kbp. However, the hybrid and LR assemblies had slightly higher genome fractions. For LR assemblies, Canu followed by Racon, and Medaka polishing was the most accurate approach. The cgSNP approach, based on either reads or assemblies, was more robust than the cgMLST approach, especially for LR sequence data. In conclusion, anti-TB drug-resistance prediction, particularly with only LR sequence data, remains challenging, especially for first-line drugs. In addition, SR assemblies appear more accurate than LR ones, and reproducible phylogeny can be achieved using cgSNP approaches.
Collapse
Affiliation(s)
- Nilay Peker
- University of Groningen, University Medical Center Groningen, Department of Medical Microbiology and Infection Prevention, Groningen, The Netherlands
| | - Leonard Schuele
- University of Groningen, University Medical Center Groningen, Department of Medical Microbiology and Infection Prevention, Groningen, The Netherlands
| | - Nienke Kok
- University of Groningen, University Medical Center Groningen, Department of Medical Microbiology and Infection Prevention, Groningen, The Netherlands
| | - Miguel Terrazos
- University of Bern, Institute for Infectious Diseases, Bern, Switzerland
| | | | - Jessica de Beer
- University of Groningen, University Medical Center Groningen, Department of Medical Microbiology and Infection Prevention, Groningen, The Netherlands
| | - Onno Akkerman
- University of Groningen, University Medical Center Groningen, Department of Pulmonary diseases and Tuberculosis, Groningen, The Netherlands
- University of Groningen, University Medical Center Groningen, TB Center Beatrixoord, Haren, The Netherlands
| | - Silke Peter
- University of Tübingen, Institute of Medical Microbiology and Hygiene, Tübingen, Germany
| | - Alban Ramette
- University of Bern, Institute for Infectious Diseases, Bern, Switzerland
| | - Matthias Merker
- Molecular and Experimental Mycobacteriology, Research Center Borstel, Borstel, Germany
| | - Stefan Niemann
- Molecular and Experimental Mycobacteriology, Research Center Borstel, Borstel, Germany
| | - Natacha Couto
- University of Groningen, University Medical Center Groningen, Department of Medical Microbiology and Infection Prevention, Groningen, The Netherlands
- The Milner Centre for Evolution, Department of Biology and Biochemistry, University of Bath, Bath, UK
| | - Bhanu Sinha
- University of Groningen, University Medical Center Groningen, Department of Medical Microbiology and Infection Prevention, Groningen, The Netherlands
| | - John WA Rossen
- University of Groningen, University Medical Center Groningen, Department of Medical Microbiology and Infection Prevention, Groningen, The Netherlands
- Department of Pathology, University of Utah School of Medicine, Salt Lake City, UT, USA
- IDbyDNA Inc., San Carlos, CA, USA
| |
Collapse
|
17
|
Merker M, Egbe NF, Ngangue YR, Vuchas C, Kohl TA, Dreyer V, Kuaban C, Noeske J, Niemann S, Sander MS. Transmission patterns of rifampicin resistant Mycobacterium tuberculosis complex strains in Cameroon: a genomic epidemiological study. BMC Infect Dis 2021; 21:891. [PMID: 34465301 PMCID: PMC8406724 DOI: 10.1186/s12879-021-06593-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2021] [Accepted: 08/17/2021] [Indexed: 11/29/2022] Open
Abstract
Background Determining factors affecting the transmission of rifampicin (RR) and multidrug-resistant (MDR) Mycobacterium tuberculosis complex strains under standardized tuberculosis (TB) treatment is key to control TB and prevent the evolution of drug resistance. Methods We combined bacterial whole genome sequencing (WGS) and epidemiological investigations for 37% (n = 195) of all RR/MDR-TB patients in Cameroon (2012–2015) to identify factors associated with recent transmission. Results Patients infected with a strain resistant to high-dose isoniazid, and ethambutol had 7.4 (95% CI 2.6–21.4), and 2.4 (95% CI 1.2–4.8) times increased odds of being in a WGS-cluster, a surrogate for recent transmission. Furthermore, age between 30 and 50 was positively correlated with recent transmission (adjusted OR 3.8, 95% CI 1.3–11.4). We found high drug-resistance proportions against three drugs used in the short standardized MDR-TB regimen in Cameroon, i.e. high-dose isoniazid (77.4%), ethambutol (56.9%), and pyrazinamide (43.1%). Virtually all strains were susceptible to fluoroquinolones, kanamycin, and clofazimine, and treatment outcomes were mostly favourable (87.5%). Conclusion Pre-existing resistance to high-dose isoniazid, and ethambutol is associated with recent transmission of RR/MDR strains in our study. A possible contributing factor for this observation is the absence of universal drug susceptibility testing in Cameroon, likely resulting in prolonged exposure of new RR/MDR-TB patients to sub-optimal or failing first-line drug regimens. Supplementary Information The online version contains supplementary material available at 10.1186/s12879-021-06593-8.
Collapse
Affiliation(s)
- Matthias Merker
- Molecular and Experimental Mycobacteriology, Research Center Borstel, Borstel, Germany. .,Evolution of the Resistome, Research Center Borstel, Borstel, Germany. .,German Center for Infection Research, Partner Site Hamburg-Lübeck-Borstel-Riems, Borstel, Germany.
| | - Nkongho F Egbe
- Tuberculosis Reference Laboratory Bamenda, Center for Health Promotion and Research, Bamenda, Cameroon.,School of Life Sciences, College of Science, University of Lincoln, Lincoln, England, UK
| | - Yannick R Ngangue
- Tuberculosis Reference Laboratory Bamenda, Center for Health Promotion and Research, Bamenda, Cameroon
| | - Comfort Vuchas
- Tuberculosis Reference Laboratory Bamenda, Center for Health Promotion and Research, Bamenda, Cameroon
| | - Thomas A Kohl
- Molecular and Experimental Mycobacteriology, Research Center Borstel, Borstel, Germany
| | - Viola Dreyer
- Molecular and Experimental Mycobacteriology, Research Center Borstel, Borstel, Germany
| | | | | | - Stefan Niemann
- Molecular and Experimental Mycobacteriology, Research Center Borstel, Borstel, Germany.,German Center for Infection Research, Partner Site Hamburg-Lübeck-Borstel-Riems, Borstel, Germany
| | - Melissa S Sander
- Tuberculosis Reference Laboratory Bamenda, Center for Health Promotion and Research, Bamenda, Cameroon
| |
Collapse
|
18
|
Lin D, Wang J, Cui Z, Ou J, Huang L, Wang Y. A genome epidemiological study of mycobacterium tuberculosis in subpopulations with high and low incidence rate in Guangxi, South China. BMC Infect Dis 2021; 21:840. [PMID: 34412585 PMCID: PMC8377953 DOI: 10.1186/s12879-021-06385-0] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2021] [Accepted: 06/18/2021] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Tuberculosis (TB) is caused by a bacterium called Mycobacterium tuberculosis (Mtb). China is the third in top 8 high TB burden countries and Guangxi is one of the high incidence areas in South China. Determine bacterial factors that affected TB incidence rate is a step toward Ending the TB epidemic. RESULTS Genomes of M. tuberculosis cultures from a relatively high and low incidence region in Guangxi have been sequenced. 347 of 358(96.9%) were identified as M. tuberculosis. All the strains belong to Lineage 2 and Lineage 4, except for one in Lineage 1. We found that the genetic structure of the M. tuberculosis population in each county varies enormously. Low incidence rate regions have a lower prevalence of Beijing genotypes than other regions. Four isolates which harbored mutT4-48 also had mutT2-58 mutations. It is suggested that strains from the ancestors of modern Beijing lineage is circulating in Guangxi. Strains of modern Beijing lineage (OR=2.04) were more likely to acquire drug resistances than Lineage 4. Most of the lineage differentiation SNPs are related to cell wall biosynthetic pathways. CONCLUSIONS These results provided a higher resolution to better understand the history of transmission of M. tuberculosis from/to South China. And the incidence rate of tuberculosis might be affected by bacterial population structure shaped by demographic history. Our findings also support the hypothesis that Modern Beijing lineage originated in South China.
Collapse
Affiliation(s)
- Dingwen Lin
- Department of Nutrition and School Health, Guangxi Zhuang Autonomous Region Center for Disease Control and Prevention, Nanning, China
| | - Junning Wang
- Zeta Biosciences(Shanghai) Co.,Ltd., Shanghai, China
| | - Zhezhe Cui
- Department of Tuberculosis Control, Guangxi Zhuang Autonomous Region Center for Disease Control and Prevention, Nanning, China
| | - Jing Ou
- Department of Tuberculosis Control, Guangxi Zhuang Autonomous Region Center for Disease Control and Prevention, Nanning, China
| | - Liwen Huang
- Department of Tuberculosis Control, Guangxi Zhuang Autonomous Region Center for Disease Control and Prevention, Nanning, China
| | - Ya Wang
- Zeta Biosciences(Shanghai) Co.,Ltd., Shanghai, China
| |
Collapse
|
19
|
A Whole-Genome-Based Gene-by-Gene Typing System for Standardized High-Resolution Strain Typing of Bacillus anthracis. J Clin Microbiol 2021; 59:e0288920. [PMID: 33827898 PMCID: PMC8218748 DOI: 10.1128/jcm.02889-20] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
Whole-genome sequencing (WGS) has been established for bacterial subtyping and is regularly used to study pathogen transmission, to investigate outbreaks, and to perform routine surveillance. Core-genome multilocus sequence typing (cgMLST) is a bacterial subtyping method that uses WGS data to provide a high-resolution strain characterization. This study aimed at developing a novel cgMLST scheme for Bacillus anthracis, a notorious pathogen that causes anthrax in livestock and humans worldwide. The scheme comprises 3,803 genes that were conserved in 57 B. anthracis genomes spanning the whole phylogeny. The scheme has been evaluated and applied to 584 genomes from 50 countries. On average, 99.5% of the cgMLST targets were detected. The cgMLST results confirmed the classical canonical single-nucleotide-polymorphism (SNP) grouping of B. anthracis into major clades and subclades. Genetic distances calculated based on cgMLST were comparable to distances from whole-genome-based SNP analysis with similar phylogenetic topology and comparable discriminatory power. Additionally, the application of the cgMLST scheme to anthrax outbreaks from Germany and Italy led to a definition of a cutoff threshold of five allele differences to trace epidemiologically linked strains for cluster typing and transmission analysis. Finally, the association of two clusters of B. anthracis with human cases of injectional anthrax in four European countries was confirmed using cgMLST. In summary, this study presents a novel cgMLST scheme that provides high-resolution strain genotyping for B. anthracis. This scheme can be used in parallel with SNP typing methods to facilitate rapid and harmonized interlaboratory comparisons, essential for global surveillance and outbreak analysis. The scheme is publicly available for application by users, including those with little bioinformatics knowledge.
Collapse
|
20
|
Sanchini A, Jandrasits C, Tembrockhaus J, Kohl TA, Utpatel C, Maurer FP, Niemann S, Haas W, Renard BY, Kröger S. Improving tuberculosis surveillance by detecting international transmission using publicly available whole genome sequencing data. ACTA ACUST UNITED AC 2021; 26. [PMID: 33446303 PMCID: PMC7809720 DOI: 10.2807/1560-7917.es.2021.26.2.1900677] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
IntroductionImproving the surveillance of tuberculosis (TB) is especially important for multidrug-resistant (MDR) and extensively drug-resistant (XDR) TB. The large amount of publicly available whole genome sequencing (WGS) data for TB gives us the chance to re-use data and to perform additional analyses at a large scale.AimWe assessed the usefulness of raw WGS data of global MDR/XDR Mycobacterium tuberculosis isolates available from public repositories to improve TB surveillance.MethodsWe extracted raw WGS data and the related metadata of M. tuberculosis isolates available from the Sequence Read Archive. We compared this public dataset with WGS data and metadata of 131 MDR- and XDR M. tuberculosis isolates from Germany in 2012 and 2013.ResultsWe aggregated a dataset that included 1,081 MDR and 250 XDR isolates among which we identified 133 molecular clusters. In 16 clusters, the isolates were from at least two different countries. For example, Cluster 2 included 56 MDR/XDR isolates from Moldova, Georgia and Germany. When comparing the WGS data from Germany with the public dataset, we found that 11 clusters contained at least one isolate from Germany and at least one isolate from another country. We could, therefore, connect TB cases despite missing epidemiological information.ConclusionWe demonstrated the added value of using WGS raw data from public repositories to contribute to TB surveillance. Comparing the German with the public dataset, we identified potential international transmission events. Thus, using this approach might support the interpretation of national surveillance results in an international context.
Collapse
Affiliation(s)
- Andrea Sanchini
- These authors contributed equally to this manuscript.,Respiratory Infections Unit (FG36), Department of Infectious Disease Epidemiology, Robert Koch Institute, Berlin, Germany
| | - Christine Jandrasits
- Bioinformatics Unit (MF1), Department of Methodology and Research Infrastructure, Robert Koch Institute, Berlin, Germany.,These authors contributed equally to this manuscript
| | - Julius Tembrockhaus
- Bioinformatics Unit (MF1), Department of Methodology and Research Infrastructure, Robert Koch Institute, Berlin, Germany
| | - Thomas Andreas Kohl
- German Center for Infection Research (DZIF), partner site Hamburg - Lübeck - Borstel - Riems, Germany.,Molecular and Experimental Mycobacteriology, Research Center Borstel, Borstel, Germany
| | - Christian Utpatel
- German Center for Infection Research (DZIF), partner site Hamburg - Lübeck - Borstel - Riems, Germany.,Molecular and Experimental Mycobacteriology, Research Center Borstel, Borstel, Germany
| | - Florian P Maurer
- National and WHO Supranational Reference Laboratory for Mycobacteria, Research Center Borstel, Borstel, Germany
| | - Stefan Niemann
- German Center for Infection Research (DZIF), partner site Hamburg - Lübeck - Borstel - Riems, Germany.,Molecular and Experimental Mycobacteriology, Research Center Borstel, Borstel, Germany
| | - Walter Haas
- Respiratory Infections Unit (FG36), Department of Infectious Disease Epidemiology, Robert Koch Institute, Berlin, Germany
| | - Bernhard Y Renard
- Hasso Plattner Institute, Faculty for Digital Engineering, University of Potsdam, Potsdam, Germany.,Bioinformatics Unit (MF1), Department of Methodology and Research Infrastructure, Robert Koch Institute, Berlin, Germany
| | - Stefan Kröger
- German Center for Infection Research (DZIF), partner site Hannover - Brunswick, Germany.,Respiratory Infections Unit (FG36), Department of Infectious Disease Epidemiology, Robert Koch Institute, Berlin, Germany
| |
Collapse
|
21
|
Establishment and Evaluation of a Core Genome Multilocus Sequence Typing Scheme for Whole-Genome Sequence-Based Typing of Pseudomonas aeruginosa. J Clin Microbiol 2021; 59:JCM.01987-20. [PMID: 33328175 DOI: 10.1128/jcm.01987-20] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2020] [Accepted: 12/07/2020] [Indexed: 01/04/2023] Open
Abstract
The environmental bacterium Pseudomonas aeruginosa, particularly multidrug-resistant clones, is often associated with nosocomial infections and outbreaks. Today, core genome multilocus sequence typing (cgMLST) is frequently applied to delineate sporadic cases from nosocomial transmissions. However, until recently, no cgMLST scheme for a standardized typing of P. aeruginosa was available. To establish a novel cgMLST scheme for P. aeruginosa, we initially determined the breadth of the P. aeruginosa population based on MLST data with a Bayesian approach (BAPS). Using genomic data of representative isolates for the whole population and all 12 serogroups, we extracted target genes and further refined them using a random data set of 1,000 P. aeruginosa genomes. Subsequently, we investigated reproducibility and discriminatory ability with repeatedly sequenced isolates and isolates from well-defined outbreak scenarios, respectively, and compared clustering applying two recently published cgMLST schemes. BAPS generated seven P. aeruginosa groups. To cover these and all serogroups, 15 reference strains were used to determine genes common in all strains. After refinement with the data set of 1,000 genomes, the cgMLST scheme consisted of 3,867 target genes, which are representative of the P. aeruginosa population and highly reproducible using biological replicates. We finally evaluated the scheme by reanalyzing two published outbreaks where the authors used single-nucleotide polymorphism (SNP) typing. In both cases, cgMLST was concordant with the previous SNP results and the results of the two other cgMLST schemes. In conclusion, the highly reproducible novel P. aeruginosa cgMLST scheme facilitates outbreak investigations due to the publicly available cgMLST nomenclature.
Collapse
|
22
|
Tagliani E, Anthony R, Kohl TA, de Neeling A, Nikolayevskyy V, Ködmön C, Maurer FP, Niemann S, van Soolingen D, van der Werf MJ, Cirillo DM. Use of a whole genome sequencing-based approach for Mycobacterium tuberculosis surveillance in Europe in 2017-2019: an ECDC pilot study. Eur Respir J 2021; 57:13993003.02272-2020. [PMID: 32732329 PMCID: PMC7784142 DOI: 10.1183/13993003.02272-2020] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2020] [Accepted: 07/17/2020] [Indexed: 01/31/2023]
Abstract
Whole genome sequencing (WGS) can be used for molecular typing and characterisation of Mycobacterium tuberculosis complex (MTBC) strains. We evaluated the systematic use of a WGS-based approach for MTBC surveillance involving all European Union/European Economic Area (EU/EEA) countries and highlight the challenges and lessons learnt to be considered for the future development of a WGS-based surveillance system. WGS and epidemiological data of patients with rifampicin-resistant (RR) and multidrug-resistant (MDR) tuberculosis (TB) were collected from EU/EEA countries between January 2017 and December 2019. WGS-based genetic relatedness analysis was performed using a standardised approach including both core genome multilocus sequence typing (cgMLST) and single nucleotide polymorphism (SNP)-based calculation of distances on all WGS data that fulfilled minimum quality criteria to ensure data comparability. A total of 2218 RR/MDR-MTBC isolates were collected from 25 countries. Among these, 56 cross-border clusters with increased likelihood of recent transmission (≤5 SNPs distance) comprising 316 RR/MDR-MTBC isolates were identified. The cross-border clusters included between two and 30 resistant isolates from two to six countries, demonstrating different RR/MDR-TB transmission patterns in Western and Eastern EU countries. This pilot study shows that a WGS-based surveillance system is not only feasible but can efficiently elucidate the dynamics of in-country and cross-border RR/MDR-TB transmission across EU/EEA countries. Lessons learnt from this study highlight that the establishment of an EU/EEA centralised WGS-based surveillance system for TB will require strengthening of national integrated systems performing prospective WGS surveillance and the development of clear procedures to facilitate international collaboration for the investigation of cross-border clusters. The implementation of a WGS-based surveillance system for monitoring the emergence of MDR-TB outbreaks in Europe is feasible and has the potential to provide supporting evidence to better elucidate cross-border transmission patternshttps://bit.ly/2ZTnPjk
Collapse
Affiliation(s)
- Elisa Tagliani
- Emerging Bacterial Pathogens Unit, Division of Immunology, Transplantation and Infectious Diseases, IRCCS San Raffaele Scientific Institute, Milan, Italy
| | - Richard Anthony
- Tuberculosis Reference Laboratory, Infectious Diseases Research, Diagnostics and Laboratory Surveillance (IDS), National Institute for Public Health and the Environment (RIVM), Bilthoven, The Netherlands.,Both authors contributed equally
| | - Thomas A Kohl
- Molecular and Experimental Mycobacteriology, National Reference Center for Mycobacteria, Research Center, Borstel, Germany.,German Center for Infection Research, partner site Borstel-Hamburg-Lübeck-Riems, Borstel, Germany.,Both authors contributed equally
| | - Albert de Neeling
- Tuberculosis Reference Laboratory, Infectious Diseases Research, Diagnostics and Laboratory Surveillance (IDS), National Institute for Public Health and the Environment (RIVM), Bilthoven, The Netherlands
| | | | - Csaba Ködmön
- European Centre for Disease Prevention and Control, Stockholm, Sweden
| | - Florian P Maurer
- Diagnostic Mycobacteriology, National Reference Center for Mycobacteria, Borstel, Germany.,Institute of Medical Microbiology, Virology and Hygiene, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Stefan Niemann
- Molecular and Experimental Mycobacteriology, National Reference Center for Mycobacteria, Research Center, Borstel, Germany.,German Center for Infection Research, partner site Borstel-Hamburg-Lübeck-Riems, Borstel, Germany
| | - Dick van Soolingen
- Tuberculosis Reference Laboratory, Infectious Diseases Research, Diagnostics and Laboratory Surveillance (IDS), National Institute for Public Health and the Environment (RIVM), Bilthoven, The Netherlands
| | | | | | | |
Collapse
|
23
|
Napier G, Campino S, Merid Y, Abebe M, Woldeamanuel Y, Aseffa A, Hibberd ML, Phelan J, Clark TG. Robust barcoding and identification of Mycobacterium tuberculosis lineages for epidemiological and clinical studies. Genome Med 2020; 12:114. [PMID: 33317631 PMCID: PMC7734807 DOI: 10.1186/s13073-020-00817-3] [Citation(s) in RCA: 58] [Impact Index Per Article: 14.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2020] [Accepted: 12/03/2020] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Tuberculosis, caused by bacteria in the Mycobacterium tuberculosis complex (MTBC), is a major global public health burden. Strain-specific genomic diversity in the known lineages of MTBC is an important factor in pathogenesis that may affect virulence, transmissibility, host response and emergence of drug resistance. Fast and accurate tracking of MTBC strains is therefore crucial for infection control, and our previous work developed a 62-single nucleotide polymorphism (SNP) barcode to inform on the phylogenetic identity of 7 human lineages and 64 sub-lineages. METHODS To update this barcode, we analysed whole genome sequencing data from 35,298 MTBC isolates (~ 1 million SNPs) covering 9 main lineages and 3 similar animal-related species (M. tuberculosis var. bovis, M. tuberculosis var. caprae and M. tuberculosis var. orygis). The data was partitioned into training (N = 17,903, 50.7%) and test (N = 17,395, 49.3%) sets and were analysed using an integrated phylogenetic tree and population differentiation (FST) statistical approach. RESULTS By constructing a phylogenetic tree on the training MTBC isolates, we characterised 90 lineages or sub-lineages or species, of which 30 are new, and identified 421 robust barcoding mutations, of which a minimal set of 90 was selected that included 20 markers from the 62-SNP barcode. The barcoding SNPs (90 and 421) discriminated perfectly the 86 MTBC isolate (sub-)lineages in the test set and could accurately reconstruct the clades across the combined 35k samples. CONCLUSIONS The validated 90 SNPs can be used for the rapid diagnosis and tracking of MTBC strains to assist public health surveillance and control. To facilitate this, the SNP markers have now been incorporated into the TB-Profiler informatics platform ( https://github.com/jodyphelan/TBProfiler ).
Collapse
Affiliation(s)
- Gary Napier
- Faculty of Infectious and Tropical Diseases, London School of Hygiene & Tropical Medicine, Keppel Street, London, WC1E 7HT UK
| | - Susana Campino
- Faculty of Infectious and Tropical Diseases, London School of Hygiene & Tropical Medicine, Keppel Street, London, WC1E 7HT UK
| | - Yared Merid
- Armauer Hansen Research Institute, Addis Ababa, Ethiopia
- Department of Microbiology, Immunology and Parasitology, College of Health Sciences, Addis Ababa University, Addis Ababa, Ethiopia
- Hawassa University College of Medicine and Health Sciences, Hawassa, Ethiopia
| | - Markos Abebe
- Armauer Hansen Research Institute, Addis Ababa, Ethiopia
| | - Yimtubezinash Woldeamanuel
- Department of Microbiology, Immunology and Parasitology, College of Health Sciences, Addis Ababa University, Addis Ababa, Ethiopia
| | - Abraham Aseffa
- Armauer Hansen Research Institute, Addis Ababa, Ethiopia
| | - Martin L. Hibberd
- Faculty of Infectious and Tropical Diseases, London School of Hygiene & Tropical Medicine, Keppel Street, London, WC1E 7HT UK
| | - Jody Phelan
- Faculty of Infectious and Tropical Diseases, London School of Hygiene & Tropical Medicine, Keppel Street, London, WC1E 7HT UK
| | - Taane G. Clark
- Faculty of Infectious and Tropical Diseases, London School of Hygiene & Tropical Medicine, Keppel Street, London, WC1E 7HT UK
- Faculty of Epidemiology and Population Health, London School of Hygiene & Tropical Medicine, London, WC1E 7HT UK
| |
Collapse
|
24
|
Jajou R, Kohl TA, Walker T, Norman A, Cirillo DM, Tagliani E, Niemann S, de Neeling A, Lillebaek T, Anthony RM, van Soolingen D. Towards standardisation: comparison of five whole genome sequencing (WGS) analysis pipelines for detection of epidemiologically linked tuberculosis cases. ACTA ACUST UNITED AC 2020; 24. [PMID: 31847944 PMCID: PMC6918587 DOI: 10.2807/1560-7917.es.2019.24.50.1900130] [Citation(s) in RCA: 32] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
Abstract
Background Whole genome sequencing (WGS) is a reliable tool for studying tuberculosis (TB) transmission. WGS data are usually processed by custom-built analysis pipelines with little standardisation between them. Aim To compare the impact of variability of several WGS analysis pipelines used internationally to detect epidemiologically linked TB cases. Methods From the Netherlands, 535 Mycobacterium tuberculosis complex (MTBC) strains from 2016 were included. Epidemiological information obtained from municipal health services was available for all mycobacterial interspersed repeat unit-variable number of tandem repeat (MIRU-VNTR) clustered cases. WGS data was analysed using five different pipelines: one core genome multilocus sequence typing (cgMLST) approach and four single nucleotide polymorphism (SNP)-based pipelines developed in Oxford, United Kingdom; Borstel, Germany; Bilthoven, the Netherlands and Copenhagen, Denmark. WGS clusters were defined using a maximum pairwise distance of 12 SNPs/alleles. Results The cgMLST approach and Oxford pipeline clustered all epidemiologically linked cases, however, in the other three SNP-based pipelines one epidemiological link was missed due to insufficient coverage. In general, the genetic distances varied between pipelines, reflecting different clustering rates: the cgMLST approach clustered 92 cases, followed by 84, 83, 83 and 82 cases in the SNP-based pipelines from Copenhagen, Oxford, Borstel and Bilthoven respectively. Conclusion Concordance in ruling out epidemiological links was high between pipelines, which is an important step in the international validation of WGS data analysis. To increase accuracy in identifying TB transmission clusters, standardisation of crucial WGS criteria and creation of a reference database of representative MTBC sequences would be advisable.
Collapse
Affiliation(s)
- Rana Jajou
- These authors contributed equally.,Center of Epidemiology and Surveillance of infectious diseases, National Institute for Public Health and the Environment (RIVM), Bilthoven, the Netherlands.,Tuberculosis Reference Laboratory, National Institute for Public Health and the Environment (RIVM), Bilthoven, the Netherlands
| | - Thomas A Kohl
- German Center for Infection Research, Borstel site, Borstel, Germany.,Molecular and Experimental Mycobacteriology, Forschungszentrum Borstel, Borstel, Germany.,These authors contributed equally
| | - Timothy Walker
- Nuffield Department of Medicine, University of Oxford, John Radcliffe Hospital, Oxford, United Kingdom
| | - Anders Norman
- International Reference Laboratory of Mycobacteriology, Statens Serum Institut, Copenhagen, Denmark
| | - Daniela Maria Cirillo
- Emerging Bacterial Pathogens Unit, Division of Immunology, Transplantation and Infectious Diseases, IRCCS San Raffaele Scientific Institute, Milan, Italy
| | - Elisa Tagliani
- Emerging Bacterial Pathogens Unit, Division of Immunology, Transplantation and Infectious Diseases, IRCCS San Raffaele Scientific Institute, Milan, Italy
| | - Stefan Niemann
- German Center for Infection Research (DZIF), Partner Site Hamburg-Lübeck-Borstel-Riems, Germany.,Molecular and Experimental Mycobacteriology, Forschungszentrum Borstel, Borstel, Germany
| | - Albert de Neeling
- Tuberculosis Reference Laboratory, National Institute for Public Health and the Environment (RIVM), Bilthoven, the Netherlands
| | - Troels Lillebaek
- Global Health Section, Department of Public Health, University of Copenhagen, Copenhagen, Denmark.,International Reference Laboratory of Mycobacteriology, Statens Serum Institut, Copenhagen, Denmark
| | - Richard M Anthony
- Tuberculosis Reference Laboratory, National Institute for Public Health and the Environment (RIVM), Bilthoven, the Netherlands
| | - Dick van Soolingen
- Tuberculosis Reference Laboratory, National Institute for Public Health and the Environment (RIVM), Bilthoven, the Netherlands
| |
Collapse
|
25
|
Guimaraes AMS, Zimpel CK. Mycobacterium bovis: From Genotyping to Genome Sequencing. Microorganisms 2020; 8:E667. [PMID: 32375210 PMCID: PMC7285088 DOI: 10.3390/microorganisms8050667] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2020] [Revised: 04/17/2020] [Accepted: 04/21/2020] [Indexed: 12/15/2022] Open
Abstract
Mycobacterium bovis is the main pathogen of bovine, zoonotic, and wildlife tuberculosis. Despite the existence of programs for bovine tuberculosis (bTB) control in many regions, the disease remains a challenge for the veterinary and public health sectors, especially in developing countries and in high-income nations with wildlife reservoirs. Current bTB control programs are mostly based on test-and-slaughter, movement restrictions, and post-mortem inspection measures. In certain settings, contact tracing and surveillance has benefited from M. bovis genotyping techniques. More recently, whole-genome sequencing (WGS) has become the preferential technique to inform outbreak response through contact tracing and source identification for many infectious diseases. As the cost per genome decreases, the application of WGS to bTB control programs is inevitable moving forward. However, there are technical challenges in data analyses and interpretation that hinder the implementation of M. bovis WGS as a molecular epidemiology tool. Therefore, the aim of this review is to describe M. bovis genotyping techniques and discuss current standards and challenges of the use of M. bovis WGS for transmission investigation, surveillance, and global lineages distribution. We compiled a series of associated research gaps to be explored with the ultimate goal of implementing M. bovis WGS in a standardized manner in bTB control programs.
Collapse
Affiliation(s)
- Ana M. S. Guimaraes
- Laboratory of Applied Research in Mycobacteria, Department of Microbiology, University of São Paulo, São Paulo 01246-904, Brazil;
| | - Cristina K. Zimpel
- Laboratory of Applied Research in Mycobacteria, Department of Microbiology, University of São Paulo, São Paulo 01246-904, Brazil;
- Department of Preventive Veterinary Medicine and Animal Health, University of São Paulo, São Paulo 01246-904, Brazil
| |
Collapse
|
26
|
Schwœbel V, Trébucq A, Kashongwe Z, Bakayoko AS, Kuaban C, Noeske J, Harouna SH, Souleymane MB, Piubello A, Ciza F, Fikouma V, Gasana M, Ouedraogo M, Gninafon M, Van Deun A, Tagliani E, Cirillo DM, Koura KG, Rieder HL. Outcomes of a nine-month regimen for rifampicin-resistant tuberculosis up to 24 months after treatment completion in nine African countries. EClinicalMedicine 2020; 20:100268. [PMID: 32300732 PMCID: PMC7152829 DOI: 10.1016/j.eclinm.2020.100268] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 09/27/2019] [Revised: 12/09/2019] [Accepted: 01/15/2020] [Indexed: 10/26/2022] Open
Abstract
BACKGROUND Treatment outcomes of the shorter regimen for rifampicin-resistant tuberculosis are not completely established. We report on these outcomes two years after treatment completion among patients enrolled in an observational cohort study in nine African countries. METHODS 1,006 patients treated with the nine-month regimen were followed every six months with sputum cultures up to 24 months after treatment completion. The risk of any unfavourable outcome, of failure and relapse, and of death during and after treatment was analysed according to patient's characteristics and initial drug susceptibility by Cox proportional hazard models. FINDINGS Respectively 67.8% and 57.2% patients had >=1 culture result six months and 12 months after treatment completion. Fourteen relapses were diagnosed. The probability of relapse-free success was 79.3% (95% confidence interval [CI] 76.6-82.0%) overall, 80.9% (95% CI 78.0-84.0%) among HIV-negative and 72.5% (95% CI 66.5-78.9%) among HIV-infected patients. Initial fluoroquinolone (adjusted hazard ratio [aHR] 6.7 [95% CI 3.4-13.1]) and isoniazid resistance (aHR 9.4 [95% CI 1.3-68.0]) were significantly associated with increased risk of failure/relapse and of any unfavourable outcome. INTERPRETATION The close to 80% relapse-free success indicates the good outcome of the regimen in low-and middle-income settings. Results confirm the lesser effectiveness of the regimen in patients with initial resistance to fluoroquinolones and support the use of high-dose isoniazid, but do not support exclusion of patients for resistance to drugs other than fluoroquinolones. FUNDING Expertise-France and Agence Française de Développement.
Collapse
Affiliation(s)
- Valérie Schwœbel
- International Union Against Tuberculosis and Lung Disease, 68 boulevard Saint-Michel, 75006 Paris, France
| | - Arnaud Trébucq
- International Union Against Tuberculosis and Lung Disease, 68 boulevard Saint-Michel, 75006 Paris, France
| | - Zacharie Kashongwe
- Kinshasa University School of Medicine, Kinshasa, Democratic Republic of Congo
| | | | | | | | | | | | - Alberto Piubello
- International Union Against Tuberculosis and Lung Disease, 68 boulevard Saint-Michel, 75006 Paris, France
- Action Damien, Niamey, Niger
| | | | - Valentin Fikouma
- Faculté des Sciences de la Santé, Bangui, Central African Republic
| | | | | | | | - Armand Van Deun
- Mycobacteriology Unit, Institute of Tropical Medicine, Antwerp, Belgium
| | - Elisa Tagliani
- Emerging Bacterial Pathogens, IRCCS San Raffaele Scientific Institute, Milan, Italy
| | - Daniela M. Cirillo
- Emerging Bacterial Pathogens, IRCCS San Raffaele Scientific Institute, Milan, Italy
| | - Kobto G. Koura
- International Union Against Tuberculosis and Lung Disease, 68 boulevard Saint-Michel, 75006 Paris, France
- Institut de Recherche pour le Développement, UMR216, Paris, France
- COMUE Sorbonne Paris Cité, Université Paris Descartes, Faculté des Sciences Pharmaceutiques et Biologiques, Paris, France
- École Nationale de Formation des Techniciens Supérieurs en Santé Publique et en Surveillance Epidémiologique, Université de Parakou, Parakou, Bénin
| | - Hans L. Rieder
- Epidemiology, Biostatistics and Prevention Institute, University of Zurich, Switzerland
- Tuberculosis Consultant Services, Kirchlindach, Switzerland
| |
Collapse
|
27
|
Jandrasits C, Kröger S, Haas W, Renard BY. Computational pan-genome mapping and pairwise SNP-distance improve detection of Mycobacterium tuberculosis transmission clusters. PLoS Comput Biol 2019; 15:e1007527. [PMID: 31815935 PMCID: PMC6922483 DOI: 10.1371/journal.pcbi.1007527] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2019] [Revised: 12/19/2019] [Accepted: 11/03/2019] [Indexed: 12/30/2022] Open
Abstract
Next-generation sequencing based base-by-base distance measures have become an integral complement to epidemiological investigation of infectious disease outbreaks. This study introduces PANPASCO, a computational pan-genome mapping based, pairwise distance method that is highly sensitive to differences between cases, even when located in regions of lineage specific reference genomes. We show that our approach is superior to previously published methods in several datasets and across different Mycobacterium tuberculosis lineages, as its characteristics allow the comparison of a high number of diverse samples in one analysis—a scenario that becomes more and more likely with the increased usage of whole-genome sequencing in transmission surveillance. Tuberculosis still is a threat to global health. It is essential to detect and interrupt transmissions to stop the spread of this infectious disease. With the rising use of next-generation sequencing methods, its application in the surveillance of Mycobacterium tuberculosis has become increasingly important in the last years. The main goal of molecular surveillance is the identification of patient-patient transmission and cluster detection. The mutation rate of M. tuberculosis is very low and stable. Therefore, many existing methods for comparative analysis of isolates provide inadequate results since their resolution is too limited. There is a need for a method that takes every detectable difference into account. We developed PANPASCO, a novel approach for comparing pairs of isolates using all genomic information available for each pair. We combine improved SNP-distance calculation with the use of a pan-genome incorporating more than 100 M. tuberculosis reference genomes representing lineages 1-4 for read mapping prior to variant detection. We thereby enable the collective analysis and comparison of similar and diverse isolates associated with different M. tuberculosis strains.
Collapse
Affiliation(s)
| | - Stefan Kröger
- Respiratory Infections Unit, Robert Koch Institute, Berlin, Germany
| | - Walter Haas
- Respiratory Infections Unit, Robert Koch Institute, Berlin, Germany
| | | |
Collapse
|
28
|
Al-Mutairi NM, Ahmad S, Mokaddas EM. Molecular characterization of multidrug-resistant Mycobacterium tuberculosis (MDR-TB) isolates identifies local transmission of infection in Kuwait, a country with a low incidence of TB and MDR-TB. Eur J Med Res 2019; 24:38. [PMID: 31806020 PMCID: PMC6894303 DOI: 10.1186/s40001-019-0397-2] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2019] [Accepted: 11/27/2019] [Indexed: 02/07/2023] Open
Abstract
Background Increasing incidence of multidrug-resistant Mycobacterium tuberculosis infections is hampering global tuberculosis control efforts. Kuwait is a low-tuberculosis-incidence country, and ~ 1% of M. tuberculosis strains are resistant to rifampicin and isoniazid (MDR-TB). This study detected mutations in seven genes predicting resistance to rifampicin, isoniazid, pyrazinamide, ethambutol and streptomycin in MDR-TB strains. Sequence data were combined with spoligotypes for detecting local transmission of MDR-TB in Kuwait. Methods Ninety-three MDR-TB strains isolated from 12 Kuwaiti and 81 expatriate patients and 50 pansusceptible strains were used. Phenotypic drug susceptibility was determined by MGIT 460 TB/960 system. Mutations conferring resistance to rifampicin, isoniazid, pyrazinamide, ethambutol and streptomycin were detected by genotype MTBDRplus assay and/or PCR sequencing of three rpoB regions, katG codon 315 (katG315) + inhA regulatory region, pncA, three embB regions and rpsL + rrs-500–900 regions. Spoligotyping kit was used, spoligotypes were identified by SITVIT2, and phylogenetic tree was constructed by using MIRU-VNTRplus software. Phylogenetic tree was also constructed from concatenated sequences by MEGA7 software. Additional PCR sequencing of gidB and rpsA was performed for cluster isolates. Results Pansusceptible isolates contained wild-type sequences. Mutations in rpoB and katG and/or inhA were detected in 93/93 and 92/93 MDR-TB strains, respectively. Mutations were also detected for pyrazinamide resistance, ethambutol resistance and streptomycin resistance in MDR-TB isolates in pncA, embB and rpsL + rrs, respectively. Spoligotyping identified 35 patterns with 18 isolates exhibiting unique patterns while 75 isolates grouped in 17 patterns. Beijing genotype was most common (32/93), and 11 isolates showed nine orphan patterns. Phylogenetic analysis of concatenated sequences showed unique patterns for 51 isolates while 42 isolates grouped in 16 clusters. Interestingly, 22 isolates in eight clusters by both methods were isolated from TB patients typically within a span of 2 years. Five of eight clusters were confirmed by additional gidB and rpsA sequence data. Conclusions Our study provides the first insight into molecular epidemiology of MDR-TB in Kuwait and identified several potential clusters of local transmission of MDR-TB involving 2–6 subjects which had escaped detection by routine surveillance studies. Prospective detection of resistance-conferring mutations can identify possible cases of local transmission of MDR-TB in low MDR-TB settings.
Collapse
Affiliation(s)
- Noura M Al-Mutairi
- Department of Microbiology, Faculty of Medicine, Kuwait University, P. O. Box 24923, 13110, Safat, Kuwait
| | - Suhail Ahmad
- Department of Microbiology, Faculty of Medicine, Kuwait University, P. O. Box 24923, 13110, Safat, Kuwait.
| | - Eiman M Mokaddas
- Department of Microbiology, Faculty of Medicine, Kuwait University, P. O. Box 24923, 13110, Safat, Kuwait
| |
Collapse
|
29
|
García de Viedma D. Pathways and strategies followed in the genomic epidemiology of Mycobacterium tuberculosis. INFECTION GENETICS AND EVOLUTION 2019; 72:4-9. [DOI: 10.1016/j.meegid.2019.01.027] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/29/2018] [Revised: 01/15/2019] [Accepted: 01/18/2019] [Indexed: 10/27/2022]
|
30
|
Barcellini L, Borroni E, Cimaglia C, Girardi E, Matteelli A, Marchese V, Stancanelli G, Abubakar I, Cirillo DM. App-based symptoms screening with Xpert MTB/RIF Ultra assay used for active tuberculosis detection in migrants at point of arrivals in Italy: The E-DETECT TB intervention analysis. PLoS One 2019; 14:e0218039. [PMID: 31260481 PMCID: PMC6602175 DOI: 10.1371/journal.pone.0218039] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2019] [Accepted: 05/23/2019] [Indexed: 11/26/2022] Open
Abstract
Background From 2014 to 2017, the number of migrants who came to Italy via the Mediterranean route has reached an unprecedented level. The majority of refugees and migrants were rescued in the Central Mediterranean and disembarked at ports in the Sicily region. Rapid on-spot active TB screening intervention at the point of arrival will cover most migrants arriving in EU and by detecting TB prevalent cases will limit further transmission of the disease. Material and methods Between November 2016 and December 2017 newly arrived migrants at point of arrivals in Sicily, were screened for active Tuberculosis using a smartphone application, followed in symptomatic individuals by fast molecular test, Xpert MTB/RIF Ultra, on collected sputum samples. Results In the study period 3787 migrants received a medical evaluation. Eight hundred and ninety-one (23.5%) reported at least one protocol-defined Tuberculosis symptom. Fifteen (2.7%) were positive to at least one microbiological test revealing a post-entry screening prevalence rate of 396 per 100.000 individuals screened (95% CI: 2.22–6.53). In logistic regression analysis, those with cough and at least one other symptom had an increased probability of testing positive compared to persons with symptoms other than cough. Whole-genome-sequencing demonstrate two separate cases of transmission. Discussion To our knowledge this study reports first-time results of an active TB case finding strategy based on on-spot symptom screening using a smartphone application, followed by fast molecular test on collected sputum samples. Our preliminary findings reveal a post-entry screening prevalence rate of 396 per 100.000 individuals screened (95% CI: 2.22–6.53).
Collapse
Affiliation(s)
- Lucia Barcellini
- Emerging Bacterial Pathogens Unit, Division of Immunology and Infectious Diseases IRCCS San Raffaele Scientific Institute, Milano, Italy
| | - Emanuele Borroni
- Emerging Bacterial Pathogens Unit, Division of Immunology and Infectious Diseases IRCCS San Raffaele Scientific Institute, Milano, Italy
| | - Claudia Cimaglia
- Clinical Epidemiology Unit, National Institute for Infectious Disease “Lazzaro Spallanzani”–IRCCS, Rome, Italy
| | - Enrico Girardi
- Clinical Epidemiology Unit, National Institute for Infectious Disease “Lazzaro Spallanzani”–IRCCS, Rome, Italy
| | - Alberto Matteelli
- Institute of Infectious and Tropical Diseases, Infectious Disease, University of Brescia, Brescia, Italy
| | - Valentina Marchese
- Institute of Infectious and Tropical Diseases, Infectious Disease, University of Brescia, Brescia, Italy
| | - Giovanna Stancanelli
- Emerging Bacterial Pathogens Unit, Division of Immunology and Infectious Diseases IRCCS San Raffaele Scientific Institute, Milano, Italy
| | - Ibrahim Abubakar
- Tuberculosis Section, Health Protection Agency Colindale and University College London, London, United Kingdom
| | | | - Daniela Maria Cirillo
- Emerging Bacterial Pathogens Unit, Division of Immunology and Infectious Diseases IRCCS San Raffaele Scientific Institute, Milano, Italy
- * E-mail:
| |
Collapse
|
31
|
Whole genome sequencing of Mycobacterium tuberculosis: current standards and open issues. Nat Rev Microbiol 2019; 17:533-545. [DOI: 10.1038/s41579-019-0214-5] [Citation(s) in RCA: 155] [Impact Index Per Article: 31.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
|
32
|
Phylogenetic Analysis of Mycobacterium tuberculosis Strains in Wales by Use of Core Genome Multilocus Sequence Typing To Analyze Whole-Genome Sequencing Data. J Clin Microbiol 2019; 57:JCM.02025-18. [PMID: 30944195 DOI: 10.1128/jcm.02025-18] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2018] [Accepted: 03/27/2019] [Indexed: 01/08/2023] Open
Abstract
An inability to standardize the bioinformatic data produced by whole-genome sequencing (WGS) has been a barrier to its widespread use in tuberculosis phylogenetics. The aim of this study was to carry out a phylogenetic analysis of tuberculosis in Wales, United Kingdom, using Ridom SeqSphere software for core genome multilocus sequence typing (cgMLST) analysis of whole-genome sequencing data. The phylogenetics of tuberculosis in Wales have not previously been studied. Sixty-six Mycobacterium tuberculosis isolates (including 42 outbreak-associated isolates) from south Wales were sequenced using an Illumina platform. Isolates were assigned to principal genetic groups, single nucleotide polymorphism (SNP) cluster groups, lineages, and sublineages using SNP-calling protocols. WGS data were submitted to the Ridom SeqSphere software for cgMLST analysis and analyzed alongside 179 previously lineage-defined isolates. The data set was dominated by the Euro-American lineage, with the sublineage composition being dominated by T, X, and Haarlem family strains. The cgMLST analysis successfully assigned 58 isolates to major lineages, and the results were consistent with those obtained by traditional SNP mapping methods. In addition, the cgMLST scheme was used to resolve an outbreak of tuberculosis occurring in the region. This study supports the use of a cgMLST method for standardized phylogenetic assignment of tuberculosis isolates and for outbreak resolution and provides the first insight into Welsh tuberculosis phylogenetics, identifying the presence of the Haarlem sublineage commonly associated with virulent traits.
Collapse
|
33
|
van der Werf MJ, Ködmön C. Whole-Genome Sequencing as Tool for Investigating International Tuberculosis Outbreaks: A Systematic Review. Front Public Health 2019; 7:87. [PMID: 31058125 PMCID: PMC6478655 DOI: 10.3389/fpubh.2019.00087] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2019] [Accepted: 04/01/2019] [Indexed: 12/31/2022] Open
Abstract
Background: Whole-genome sequencing (WGS) can support the investigation of tuberculosis (TB) outbreaks. The technique has been applied to estimate the timing and directionality of transmission and to exclude cases from an investigation. This review assesses how WGS was applied in international outbreak investigations and discusses the advantages and challenges of the application of WGS. Methods: Databases were searched for reports on international TB outbreak investigations. Information was extracted on: Why was WGS applied?; How was WGS applied?; Organizational issues; WGS methodology; What was learned/what were the implications of the WGS investigation?; and challenges and lessons learned. Results: Three studies reporting on international outbreak investigations were identified. Retrospective WGS sequencing was performed in all studies and prospective typing in two to study TB transmission. In one study, WGS data were produced centrally (i.e., in one laboratory) and analysis was done centrally. In two studies, WGS data production was done in a decentralized manner, and analysis was centralized in one laboratory. Three groups of professionals were involved in the international outbreak investigation: public health authorities, laboratory experts, and clinicians. The reported WGS methodology applied differed between the studies in some aspects, e.g., sequencing platform; quality measures, percentage of the reference genome covered, and the mean genomic coverage; analysis, use of a reference genome or de novo assembly; and software used for alignment and analysis. In all three studies, in-house scripts were used for variance calling, and the single nucleotide polymorphism (SNP) approach was used for analysis. All outbreak investigation reports stated that WGS refuted suspected transmission events and provided supporting evidence for epidemiological data. Several challenges were reported of which most were not related to WGS. The only challenge related to WGS was the timeframe of getting WGS data if WGS is not routinely performed. Conclusions: WGS was considered a useful addition in international TB outbreak investigations. Further standardization of the WGS methodology and good structures for international collaboration and coordination are needed to take full advantage of this new technology. Whether the use of WGS results in earlier detection of cases and thus limits transmission still needs to be determined.
Collapse
Affiliation(s)
| | - Csaba Ködmön
- European Centre for Disease Prevention and Control, Stockholm, Sweden
| |
Collapse
|
34
|
Macedo R, Pinto M, Borges V, Nunes A, Oliveira O, Portugal I, Duarte R, Gomes JP. Evaluation of a gene-by-gene approach for prospective whole-genome sequencing-based surveillance of multidrug resistant Mycobacterium tuberculosis. Tuberculosis (Edinb) 2019; 115:81-88. [PMID: 30948181 DOI: 10.1016/j.tube.2019.02.006] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2018] [Revised: 01/15/2019] [Accepted: 02/22/2019] [Indexed: 11/27/2022]
Abstract
Whole-genome sequencing (WGS) offers unprecedented resolution for tracking Mycobacterium tuberculosis transmission and antibiotic-resistance spread. Still, the establishment of standardized WGS-based pipelines and the definition of epidemiological clusters based on genetic relatedness are under discussion. We aimed to implement a dynamic gene-by-gene approach, fully relying on freely available software, for prospective WGS-based tuberculosis surveillance, demonstrating its application for detecting transmission chains by retrospectively analysing all M/XDR strains isolated in 2013-2017 in Portugal. We observed a good correlation between genetic relatedness and epidemiological links, with strongly epilinked clusters displaying mean pairwise allele differences (AD) always below 0.3% (ratio of mean AD over the total number of shared loci between same-cluster strains). This data parallels the genetic distances acquired by the core-SNV analysis, while providing higher resolution and epidemiological concordance than MIRU-VNTR genotyping. The dynamic analysis of strain sub-sets (i.e., increasing the number of shared loci within each sub-set) also strengthens the confidence in detecting epilinked clusters. This gene-by-gene strategy also offers several practical benefits (e.g., reliance on freely-available software, scalability and low computational requirements) that further consolidated its suitability for a timely and robust prospective WGS-based laboratory surveillance of M/XDR-TB cases.
Collapse
Affiliation(s)
- Rita Macedo
- National Reference Laboratory for Mycobacteria, Department of Infectious Diseases, National Institute of Health, 1649-016 Lisbon, Portugal.
| | - Miguel Pinto
- Bioinformatics Unit, Department of Infectious Diseases, National Institute of Health, 1649-016 Lisbon, Portugal.
| | - Vítor Borges
- Bioinformatics Unit, Department of Infectious Diseases, National Institute of Health, 1649-016 Lisbon, Portugal.
| | - Alexandra Nunes
- Bioinformatics Unit, Department of Infectious Diseases, National Institute of Health, 1649-016 Lisbon, Portugal.
| | - Olena Oliveira
- Life and Health Sciences Research Institute (ICVS), School of Medicine, University of Minho, 4710-057 Braga, Portugal; ICVS/3B, PT Government Associate Laboratory, 4710-057, Braga/4805-017, Guimarães, Portugal; EPI Unit, Institute of Public Health, University of Porto, 4050-600 Porto, Portugal.
| | - Isabel Portugal
- iMed.ULisboa-Research Institute for Medicines, University of Lisbon, Lisbon, Portugal.
| | - Raquel Duarte
- EPI Unit, Institute of Public Health, University of Porto, 4050-600 Porto, Portugal; Clinical Epidemiology, Predictive Medicine and Public Health Department, Faculty of Medicine, University of Porto, 4200-319 Porto, Portugal; Pulmonology Department, Centro Hospitalar de Vila Nova de Gaia/Espinho EPE, 4400-129 Vila Nova de Gaia, Portugal.
| | - João Paulo Gomes
- Bioinformatics Unit, Department of Infectious Diseases, National Institute of Health, 1649-016 Lisbon, Portugal.
| |
Collapse
|
35
|
Meehan CJ, Moris P, Kohl TA, Pečerska J, Akter S, Merker M, Utpatel C, Beckert P, Gehre F, Lempens P, Stadler T, Kaswa MK, Kühnert D, Niemann S, de Jong BC. The relationship between transmission time and clustering methods in Mycobacterium tuberculosis epidemiology. EBioMedicine 2018; 37:410-416. [PMID: 30341041 PMCID: PMC6284411 DOI: 10.1016/j.ebiom.2018.10.013] [Citation(s) in RCA: 86] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2018] [Revised: 09/17/2018] [Accepted: 10/03/2018] [Indexed: 11/20/2022] Open
Abstract
BACKGROUND Tracking recent transmission is a vital part of controlling widespread pathogens such as Mycobacterium tuberculosis. Multiple methods with specific performance characteristics exist for detecting recent transmission chains, usually by clustering strains based on genotype similarities. With such a large variety of methods available, informed selection of an appropriate approach for determining transmissions within a given setting/time period is difficult. METHODS This study combines whole genome sequence (WGS) data derived from 324 isolates collected 2005-2010 in Kinshasa, Democratic Republic of Congo (DRC), a high endemic setting, with phylodynamics to unveil the timing of transmission events posited by a variety of standard genotyping methods. Clustering data based on Spoligotyping, 24-loci MIRU-VNTR typing, WGS based SNP (Single Nucleotide Polymorphism) and core genome multi locus sequence typing (cgMLST) typing were evaluated. FINDINGS Our results suggest that clusters based on Spoligotyping could encompass transmission events that occurred almost 200 years prior to sampling while 24-loci-MIRU-VNTR often represented three decades of transmission. Instead, WGS based genotyping applying low SNP or cgMLST allele thresholds allows for determination of recent transmission events, e.g. in timespans of up to 10 years for a 5 SNP/allele cut-off. INTERPRETATION With the rapid uptake of WGS methods in surveillance and outbreak tracking, the findings obtained in this study can guide the selection of appropriate clustering methods for uncovering relevant transmission chains within a given time-period. For high resolution cluster analyses, WGS-SNP and cgMLST based analyses have similar clustering/timing characteristics even for data obtained from a high incidence setting.
Collapse
Affiliation(s)
- Conor J Meehan
- Unit of Mycobacteriology, Biomedical Sciences, Institute of Tropical Medicine, Antwerp 2000, Belgium.
| | - Pieter Moris
- Unit of Mycobacteriology, Biomedical Sciences, Institute of Tropical Medicine, Antwerp 2000, Belgium; Adrem Data Lab (Adrem), Department of Mathematics and Computer Science, University of Antwerp, Antwerp 2020, Belgium; Biomedical Informatics Research Network Antwerp (biomina), University of Antwerp, Antwerp 2020, Belgium
| | - Thomas A Kohl
- German Center for Infection Research, Partner Site Hamburg-Lübeck-Borstel-Riems, D-23845 Borstel, Germany; Molecular and Experimental Mycobacteriology, Priority Area Infections, Research Center Borstel, D-23845 Borstel, Germany
| | - Jūlija Pečerska
- Swiss Institute of Bioinformatics (SIB), 1015 Lausanne, Switzerland
| | - Suriya Akter
- Unit of Mycobacteriology, Biomedical Sciences, Institute of Tropical Medicine, Antwerp 2000, Belgium
| | - Matthias Merker
- German Center for Infection Research, Partner Site Hamburg-Lübeck-Borstel-Riems, D-23845 Borstel, Germany; Molecular and Experimental Mycobacteriology, Priority Area Infections, Research Center Borstel, D-23845 Borstel, Germany
| | - Christian Utpatel
- German Center for Infection Research, Partner Site Hamburg-Lübeck-Borstel-Riems, D-23845 Borstel, Germany; Molecular and Experimental Mycobacteriology, Priority Area Infections, Research Center Borstel, D-23845 Borstel, Germany
| | - Patrick Beckert
- German Center for Infection Research, Partner Site Hamburg-Lübeck-Borstel-Riems, D-23845 Borstel, Germany; Molecular and Experimental Mycobacteriology, Priority Area Infections, Research Center Borstel, D-23845 Borstel, Germany
| | - Florian Gehre
- Unit of Mycobacteriology, Biomedical Sciences, Institute of Tropical Medicine, Antwerp 2000, Belgium; Vaccines and Immunity Theme, Medical Research Council Unit The Gambia, Serekunda, Gambia; Department Infectious Diseases Epidemiology, Bernhard Nocht Institute for Tropical Medicine, Hamburg 20359, Germany
| | - Pauline Lempens
- Unit of Mycobacteriology, Biomedical Sciences, Institute of Tropical Medicine, Antwerp 2000, Belgium
| | - Tanja Stadler
- Swiss Institute of Bioinformatics (SIB), 1015 Lausanne, Switzerland
| | - Michel K Kaswa
- Unit of Mycobacteriology, Biomedical Sciences, Institute of Tropical Medicine, Antwerp 2000, Belgium; National Tuberculosis Program, Kinshasa, DR Congo
| | - Denise Kühnert
- Max Planck Institute for the Science of Human History, 07745 JENA, Germany
| | - Stefan Niemann
- German Center for Infection Research, Partner Site Hamburg-Lübeck-Borstel-Riems, D-23845 Borstel, Germany; Molecular and Experimental Mycobacteriology, Priority Area Infections, Research Center Borstel, D-23845 Borstel, Germany
| | - Bouke C de Jong
- Unit of Mycobacteriology, Biomedical Sciences, Institute of Tropical Medicine, Antwerp 2000, Belgium
| |
Collapse
|
36
|
O'Toole RF. Development of a New Genome-Wide MLST Scheme for High-Resolution Typing of Diverse Mycobacterium tuberculosis Complex Strains. EBioMedicine 2018; 34:6-7. [PMID: 30087031 PMCID: PMC6116476 DOI: 10.1016/j.ebiom.2018.07.038] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2018] [Accepted: 07/30/2018] [Indexed: 11/23/2022] Open
Abstract
In this issue of EBioMedicine, Kohl and colleagues describe the development of a new core genome MLST scheme (cgMLST) for Mycobacterium tuberculosis complex strains based on a set of 2891 genes. Here, the application of the scheme to a number of tuberculosis surveillance studies is examined.
Collapse
Affiliation(s)
- Ronan F O'Toole
- School of Medicine, College of Health and Medicine, University of Tasmania, Hobart, Australia; Department of Clinical Microbiology, Trinity College Dublin, Ireland.
| |
Collapse
|