Computational protocol: Draft Genome Sequences of Two Extensively Drug Resistant Strains of Mycobacterium tuberculosis Belonging to the Euro American S Lineage

Similar protocols

Protocol publication

[…] Drug-resistant tuberculosis (TB) caused by Mycobacterium tuberculosis is a global threat and a major public health problem in several countries (). In South Africa, circulating M. tuberculosis strains are diverse (), with three spoligotypes being most common: the Beijing spoligotype predominant in Western Cape, and the Euro-American LAM4 and S spoligotypes predominant in Gauteng and KwaZulu-Natal (KZN) provinces (). The Euro-American S spoligotype (ST34) is prevalent in TB patients within KZN and Gauteng provinces ().We describe the draft genome sequences of two extensively drug-resistant (XDR) TB clinical strains of M. tuberculosis belonging to ST34. Permission to use these strains was granted by the University of Pretoria, Faculty of Health Sciences, Research Ethics Committee (206/2012). Both strains, RSA184 and RSA114, were isolated from patients from Swaziland. Spoligotyping and drug susceptibility testing (DST) were performed per standard protocols (, ).DNA was extracted from heat-killed M. tuberculosis grown on slants, using a previously described chemical method (). Illumina sequencing libraries were prepared as previously described () and sequenced using the Illumina HiSeq platform at the Broad Institute (Cambridge, MA, USA). Reads from RSA114 and RSA184 were assembled into draft genomes using ALLPATHS-LG with Pilon (). The genome assemblies of 4,416,700 bp and 4,389,272 bp for RSA114 and RSA 184, respectively, were annotated by aligning each assembly to the H37Rv genome (CP003248.2) using Nucmer (). For those genes not cleanly mapping to H37Rv, the protein-coding genes of 4,020 and 4,019 for RSA 114 and RSA 184, respectively, were predicted with Prodigal (). Both strains had 45 tRNAs identified by tRNAscan-SE () and 3 rRNA genes predicted using RNAmmer (). We also confirmed the experimental spoligotype predictions using a previously described computational spoligotyping approach ().Sequence reads were aligned to the M. tuberculosis H37Rv reference genome using BWA (), and Pilon was used to identify variants. We detected a total of 797 and 734 nonsynonymous changes relative to H37Rv for RSA114 and RSA184, respectively. We also detected nonsynonymous changes in rpoB (S450L, I491F), katG (S315T), and gyrA (D94G), previously implicated in drug resistance. Interestingly, RSA114, which lacked known resistance-conferring gyrA mutations, had 14, 7, and 4 nonsynonymous changes in genes encoding efflux pumps (EPs), phthiocerol dimycocerosates (PDIMs) and type VII secretion systems (ESXs), respectively. Drug resistance in M. tuberculosis can be acquired through mutations in EPs that increase their activity to expel a broad spectrum of antibiotics (), and ofloxacin-resistant strains, lacking DNA gyrase mutations, were found to overexpress EPs (). In RSA114, we identified mutations within the EP-encoding genes Rv0987, Rv2039c, and Rv0402c that are predicted by PROVEAN (http://provean.jcvi.org/index.php), to impact efflux activity. ESX export enzymes are involved in the synthesis of PDIM proteins, which have a role in virulence (), and are overexpressed in XDR TB strains, suggesting a contribution to this XDR-level drug resistance (). Future functional studies are needed to determine the impact of these mutations on drug resistance. […]

Pipeline specifications

Software tools ALLPATHS-LG, Pilon, MUMmer, Prodigal, tRNAscan-SE, RNAmmer, BWA, PROVEAN
Applications Genome annotation, Nucleotide sequence alignment
Organisms Osphranter robustus
Diseases Tuberculosis, Tuberculosis, Multidrug-Resistant