Dataset Information

TransFlow: a Snakemake workflow for transmission analysis of Mycobacterium tuberculosis whole-genome sequencing data.

ABSTRACT:

Motivation

Whole-genome sequencing (WGS) is increasingly used to aid the understanding of Mycobacterium tuberculosis (MTB) transmission. The epidemiological analysis of tuberculosis based on the WGS technique requires a diverse collection of bioinformatics tools. Effectively using these analysis tools in a scalable and reproducible way can be challenging, especially for non-experts.

Results

Here, we present TransFlow (Transmission Workflow), a user-friendly, fast, efficient and comprehensive WGS-based transmission analysis pipeline. TransFlow combines some state-of-the-art tools to take transmission analysis from raw sequencing data, through quality control, sequence alignment and variant calling, into downstream transmission clustering, transmission network reconstruction and transmission risk factor inference, together with summary statistics and data visualization in a summary report. TransFlow relies on Snakemake and Conda to resolve dependencies among consecutive processing steps and can be easily adapted to any computation environment.

Availability and implementation

TransFlow is free available at https://github.com/cvn001/transflow.

Supplementary information

Supplementary data are available at Bioinformatics online.

SUBMITTER: Pan J

PROVIDER: S-EPMC9825751 | biostudies-literature | 2023 Jan

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

TransFlow: a Snakemake workflow for transmission analysis of Mycobacterium tuberculosis whole-genome sequencing data.

Pan Junhang J Li Xiangchen X Zhang Mingwu M Lu Yewei Y Zhu Yelei Y Wu Kunyang K Wu Yiwen Y Wang Weixin W Chen Bin B Liu Zhengwei Z Wang Xiaomeng X Gao Junshun J

Bioinformatics (Oxford, England) 20230101 1

<h4>Motivation</h4>Whole-genome sequencing (WGS) is increasingly used to aid the understanding of Mycobacterium tuberculosis (MTB) transmission. The epidemiological analysis of tuberculosis based on the WGS technique requires a diverse collection of bioinformatics tools. Effectively using these analysis tools in a scalable and reproducible way can be challenging, especially for non-experts.<h4>Results</h4>Here, we present TransFlow (Transmission Workflow), a user-friendly, fast, efficient and co ...[more]

PMID: 36469333

Similar Datasets

Project description:Namibia is among 30 countries with a high burden of tuberculosis (TB), with an estimated incidence of 460 per 100,000 population and around 800 new multidrug-resistant (MDR) TB cases per year. Still, data on the transmission and evolution of drug-resistant Mycobacterium tuberculosis complex (Mtbc) strains are not available. Whole-genome sequencing data of 136 rifampicin-resistant (RIFr) Mtbc strains obtained from 2016 to 2018 were used for phylogenetic classification, resistance prediction, and cluster analysis and linked with phenotypic drug susceptibility testing (pDST) data. Roughly 50% of the strains investigated were resistant to all first-line drugs. Furthermore, 13% of the MDR Mtbc strains were already pre-extensively drug resistant (pre-XDR). The cluster rates were high, at 74.6% among MDR and 85% among pre-XDR strains. A significant proportion of strains had borderline resistance-conferring mutations, e.g., inhA promoter mutations or rpoB L430P. Accordingly, 25% of the RIFr strains tested susceptible by pDST. Finally, we determined a potentially new bedaquiline resistance mutation (Rv0678 D88G) occurring in two independent clusters. High rates of resistance to first-line drugs in line with emerging pre-XDR and likely bedaquiline resistance linked with the ongoing recent transmission of MDR Mtbc clones underline the urgent need for the implementation of interventions that allow rapid diagnostics to break MDR TB transmission chains in the country. A borderline RIFr mutation in the dominant outbreak strain causing discrepancies between phenotypic and genotypic resistance testing results may require breakpoint adjustments but also may allow individualized regimens with high-dose treatment. IMPORTANCE The transmission of drug-resistant tuberculosis (TB) is a major problem for global TB control. Using genome sequencing, we showed that 13% of the multidrug-resistant (MDR) M. tuberculosis complex strains from Namibia are already pre-extensively drug resistant (pre-XDR), which is substantial in an African setting. Our data also indicate that the ongoing transmission of MDR and pre-XDR strains contributes significantly to the problem. In contrast to other settings with higher rates of drug resistance, we found a high proportion of strains having so-called borderline low-level resistance mutations, e.g., inhA promoter mutations or rpoB L430P. This led to the misclassification of 25% of the rifampicin-resistant strains as susceptible by phenotypic drug susceptibility testing. This observation potentially allows individualized regimens with high-dose treatment as a potential option for patients with few treatment options. We also found a potentially new bedaquiline resistance mutation in rv0678.

Project description:Tuberculosis (TB) originating from expatriates that hail from high TB-burden countries is hypothesized to play a role in continued TB transmission in Oman. Here, we used whole-genome sequencing (WGS) to assess national TB transmission dynamics. The annual incidence per 100,000 population per year was calculated for nationals and expatriates. A convenience sample of Mycobacterium tuberculosis (MTB) isolates from 2018 to 2019 was sequenced and analyzed with publicly available TB sequences from Bangladesh, Tanzania, the Philippines, India, and Pakistan. Relatedness was assessed by generating core-genome single nucleotide polymorphism (SNP) distances. The incidence of TB was five cases per 100,000 persons in 2018 and seven cases per 100,000 persons in 2020 (R2 = 0.34, P = 0.60). Incidence among nationals was 3.9 per 100,000 persons in 2018 and 3.5 per 100,000 persons in 2020 (R2 = 0.20, P = 0.70), and incidence among expatriates was 7.2 per 100,000 persons in 2018 and 12.7 per 100,000 persons in 2020 (R2 = 0.74, P = 0.34). Sixty-eight local MTB isolates were sequenced and analyzed with 393 global isolates. Isolates belonged to nine distinct spoligotypes. Two isolates, originating from an expatriate and an Omani national, were grouped into a WGS-based cluster (SNP distance < 12), which was corroborated by an epidemiological investigation. Relatedness of local and global isolates (SNP distance < 100) was also seen. The relatedness between MTB strains in Oman and those in expatriate countries of origin can aid inform TB control policy. Our results provide evidence that WGS can complement epidemiological analysis to achieve the End TB strategy goal in Oman. IMPORTANCE Tuberculosis (TB) incidence in Oman remains above national program control targets. TB transmission originating from expatriates from high TB-burden countries has been hypothesized to play a role. We used whole-genome sequencing (WGS) to assess TB transmission dynamics between expatriates and Omani nationals to inform TB control efforts. Available Mycobacterium tuberculosis isolates from 2018 to 2019 underwent WGS and analysis with publicly available TB sequences from Bangladesh, the Philippines, India, and Pakistan to assess for genetic relatedness. Our analysis revealed evidence of previously unrecognized transmission between an expatriate and an Omani national, which was corroborated by epidemiological investigation. Analysis of local and global isolates revealed evidence of distant relatedness between local and global isolates. Our results provide evidence that WGS can complement classic public health surveillance to inform targeted interventions to achieve the End TB strategy goal in Oman.

Project description:Cambodia has one of the highest tuberculosis (TB) incidence rates in the WHO Western Pacific region. Remarkably though, the prevalence of multidrug-resistant TB (MDR-TB) remains low. We explored the genetic diversity of Mycobacterium tuberculosis (MTB) circulating in this unique setting using whole-genome sequencing (WGS). From October 2017 until January 2018, we collected one hundred sputum specimens from consenting adults older than 21 years of age, newly diagnosed with bacteriologically confirmed TB in 3 districts of Phnom Penh and Takeo provinces of Cambodia before they commence on their TB treatment, where eighty MTB isolates were successfully cultured and sequenced. Majority of the isolates belonged to Lineage 1 (Indo-Oceanic) (69/80, 86.25%), followed by Lineage 2 (East Asian) (10/80, 12.5%) and Lineage 4 (Euro-American) (1/80, 1.25%). Phenotypic resistance to both streptomycin and isoniazid was found in 3 isolates (3/80, 3.75%), while mono-resistance to streptomycin and isoniazid was identical at 2.5% (N = 2 each). None of the isolates tested was resistant to either rifampicin or ethambutol. The specificities of genotypic prediction for resistance to all drugs tested were 100%, while the sensitivities of genotypic resistance predictions to isoniazid and streptomycin were lower at 40% (2/5) and 80% (4/5) respectively. We identified 8 clusters each comprising of two to five individuals all residing in the Takeo province, making up half (28/56, 50%) of all individuals sampled in the province, indicating the presence of multiple ongoing transmission events. All clustered isolates were of Lineage 1 and none are resistant to any of the drugs tested. This study while demonstrating the relevance and utility of WGS in predicting drug resistance and inference of disease transmission, highlights the need to increase the representation of genotype-phenotype TB data from low and middle income countries in Asia and Africa to improve the accuracies for prediction of drug resistance.

Project description:Background: Tuberculosis disease, caused by Mycobacterium tuberculosis, is a major public health problem. The emergence of M. tuberculosis strains resistant to existing treatments threatens to derail control efforts. Resistance is mainly conferred by mutations in genes coding for drug targets or converting enzymes, but our knowledge of these mutations is incomplete. Whole genome sequencing (WGS) is an increasingly common approach to rapidly characterize isolates and identify mutations predicting antimicrobial resistance and thereby providing a diagnostic tool to assist clinical decision making. Methods: We applied machine learning approaches to 16,688 M. tuberculosis isolates that have undergone WGS and laboratory drug-susceptibility testing (DST) across 14 antituberculosis drugs, with 22.5% of samples being multidrug resistant and 2.1% being extensively drug resistant. We used non-parametric classification-tree and gradient-boosted-tree models to predict drug resistance and uncover any associated novel putative mutations. We fitted separate models for each drug, with and without "co-occurrent resistance" markers known to be causing resistance to drugs other than the one of interest. Predictive performance was measured using sensitivity, specificity, and the area under the receiver operating characteristic curve, assuming DST results as the gold standard. Results: The predictive performance was highest for resistance to first-line drugs, amikacin, kanamycin, ciprofloxacin, moxifloxacin, and multidrug-resistant tuberculosis (area under the receiver operating characteristic curve above 96%), and lowest for third-line drugs such as D-cycloserine and Para-aminosalisylic acid (area under the curve below 85%). The inclusion of co-occurrent resistance markers led to improved performance for some drugs and superior results when compared to similar models in other large-scale studies, which had smaller sample sizes. Overall, the gradient-boosted-tree models performed better than the classification-tree models. The mutation-rank analysis detected no new single nucleotide polymorphisms linked to drug resistance. Discordance between DST and genotypically inferred resistance may be explained by DST errors, novel rare mutations, hetero-resistance, and nongenomic drivers such as efflux-pump upregulation. Conclusion: Our work demonstrates the utility of machine learning as a flexible approach to drug resistance prediction that is able to accommodate a much larger number of predictors and to summarize their predictive ability, thus assisting clinical decision making and single nucleotide polymorphism detection in an era of increasing WGS data generation.

Dataset Information

TransFlow: a Snakemake workflow for transmission analysis of Mycobacterium tuberculosis whole-genome sequencing data.

Motivation

Results

Availability and implementation

Supplementary information

Publications

TransFlow: a Snakemake workflow for transmission analysis of Mycobacterium tuberculosis whole-genome sequencing data.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets