Tuberculosis (TB) remains one of the leading causes of morbidity and mortality worldwide. Extrapulmonary tuberculosis (EPTB) constitutes around 15–20% of TB cases in immunocompetent individuals. Extrapulmonary sites that are affected by TB include bones, lymph nodes, meningitis, pleura, and genitourinary tract. Whole genome sequencing has emerged as a powerful tool to map genetic diversity among Mycobacterium tuberculosis (MTB) isolates and identify the genomic signatures associated with drug resistance, pathogenesis, and disease transmission. Several pulmonary isolates of MTB have been sequenced over the years. However, availability of whole genome sequences of MTB isolates from extrapulmonary sites is limited. Some studies suggest that genetic variations in MTB might contribute to disease presentation in extrapulmonary sites. This can be addressed if whole genome sequence data from large number of extrapulmonary isolates becomes available. In this study, we have performed whole genome sequencing of five MTB clinical isolates derived from EPTB sites using next-generation sequencing platform. We identified 1434 nonsynonymous single nucleotide variations (SNVs), 143 insertions and 105 deletions. This includes 279 SNVs that were not reported before in publicly available datasets. We found several mutations that are known to confer resistance to drugs. All the five isolates belonged to East-African-Indian lineage (lineage 3). We identified 9 putative prophage DNA integrations and 14 predicted clustered regularly interspaced short palindromic repeats (CRISPR) in MTB genome. Our analysis indicates that more work is needed to map the genetic diversity of MTB. Whole genome sequencing in conjunction with comprehensive drug susceptibility testing can reveal clinically relevant mutations associated with drug resistance.
Get full access to this article
View all access options for this article.
References
1.
AkhtarP, SrivastavaS, SrivastavaA, SrivastavaM, SrivastavaBS, and SrivastavaR. (2006). Rv3303c of Mycobacterium tuberculosis protects tubercle bacilli against oxidative stress in vivo and contributes to virulence in mice. Microbes Infect, 8, 2855–2862.
2.
AliA, HasanZ, McNerneyR, et al. (2015). Whole genome sequencing based characterization of extensively drug-resistant Mycobacterium tuberculosis isolates from Pakistan. PLoS One, 10, e0117771.
3.
AnsariA, TalatN, JamilB, et al. (2009). Cytokine gene polymorphisms across tuberculosis clinical spectrum in Pakistani patients. PLoS One, 4, e4778.
4.
ArndtD, GrantJR, MarcuA, et al. (2016). PHASTER: A better, faster version of the PHAST phage search tool. Nucleic Acids Res, 44, W16–W21.
5.
ArnvigKB, ComasI, ThomsonNR, et al. (2011). Sequence-based analysis uncovers an abundance of non-coding RNA in the total transcriptome of Mycobacterium tuberculosis. PLoS Pathog, 7, e1002342.
6.
BanulsAL, SanouA, AnhNT, and GodreuilS. (2015). Mycobacterium tuberculosis: Ecology and evolution of a human bacterium. J Med Microbiol, 64, 1261–1269.
7.
BhuwanM, AroraN, SharmaA, et al. (2016). Interaction of Mycobacterium tuberculosis virulence factor RipA with chaperone MoxR1 is required for transport through the TAT secretion system. MBio, 7, e02259.
8.
BrandisG, and HughesD. (2013). Genetic characterization of compensatory evolution in strains carrying rpoB Ser531Leu, the rifampicin resistance mutation most frequently found in clinical isolates. J Antimicrob Chemother, 68, 2493–2497.
9.
BrandisG, WrandeM, LiljasL, and HughesD. (2012). Fitness-compensatory mutations in rifampicin-resistant RNA polymerase. Mol Microbiol, 85, 142–151.
10.
CadeCE, DlouhyAC, MedzihradszkyKF, Salas-CastilloSP, and GhiladiRA. (2010). Isoniazid-resistance conferring mutations in Mycobacterium tuberculosis KatG: Catalase, peroxidase, and INH-NADH adduct formation activities. Protein Sci, 19, 458–474.
11.
CasaliN, NikolayevskyyV, BalabanovaY, et al. (2014). Evolution and transmission of drug-resistant tuberculosis in a Russian population. Nat Genet, 46, 279–286.
12.
CawsM, ThwaitesG, DunstanS, et al. (2008). The influence of host and bacterial genotype on the development of disseminated disease with Mycobacterium tuberculosis. PLoS Pathog, 4, e1000034.
13.
ChernyaevaEN, ShulginaMV, RotkevichMS, et al. (2014). Genome-wide Mycobacterium tuberculosis variation (GMTV) database: A new tool for integrating sequence variations and epidemiology. BMC Genomics, 15, 308.
14.
ChooSW, WeeWY, NgeowYF, et al. (2014). Genomic reconnaissance of clinical isolates of emerging human pathogen Mycobacterium abscessus reveals high evolutionary potential. Sci Rep, 4, 4061.
15.
ColeST, BroschR, ParkhillJ, et al. (1998). Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence. Nature, 393, 537–544.
16.
CollF, McNerneyR, Guerra-AssuncaoJA, et al. (2014). A robust SNP barcode for typing Mycobacterium tuberculosis complex strains. Nat Commun, 5, 4812.
17.
ConverseSE, MougousJD, LeavellMD, LearyJA, BertozziCR, and CoxJS. (2003). MmpL8 is required for sulfolipid-1 biosynthesis and Mycobacterium tuberculosis virulence. Proc Natl Acad Sci U S A, 100, 6121–6126.
18.
CoscollaM, and GagneuxS. (2010). Does M. tuberculosis genomic diversity explain disease diversity?. Drug Discov Today Dis Mech, 7, e43–e59.
19.
CoscollaM, and GagneuxS. (2014). Consequences of genomic diversity in Mycobacterium tuberculosis. Semin Immunol, 26, 431–444.
20.
DasS, RoychowdhuryT, KumarP, et al. (2013). Genetic heterogeneity revealed by sequence analysis of Mycobacterium tuberculosis isolates from extra-pulmonary tuberculosis patients. BMC Genomics, 14, 404.
21.
DavidS, MateusAR, DuarteEL, et al. (2015). Determinants of the sympatric host-pathogen relationship in tuberculosis. PLoS One, 10, e0140625.
22.
de VosM, MullerB, BorrellS, et al. (2013). Putative compensatory mutations in the rpoC gene of rifampin-resistant Mycobacterium tuberculosis are associated with ongoing transmission. Antimicrob Agents Chemother, 57, 827–832.
23.
DrobniewskiF, BalabanovaY, RuddyM, et al. (2002). Rifampin- and multidrug-resistant tuberculosis in Russian civilians and prison inmates: Dominance of the Beijing strain family. Emerg Infect Dis, 8, 1320–1326.
24.
FaksriK, DrobniewskiF, NikolayevskyyV, et al. (2011). Epidemiological trends and clinical comparisons of Mycobacterium tuberculosis lineages in Thai TB meningitis. Tuberculosis (Edinb), 91, 594–600.
FernandoSL, SaundersBM, SluyterR, et al. (2007). A polymorphism in the P2X7 gene increases susceptibility to extrapulmonary tuberculosis. Am J Respir Crit Care Med, 175, 360–366.
27.
FilliolI, DriscollJR, van SoolingenD, et al. (2003). Snapshot of moving and expanding clones of Mycobacterium tuberculosis and their global distribution assessed by spoligotyping in an international study. J Clin Microbiol, 41, 1963–1970.
28.
FortierLC, and SekulovicO. (2013). Importance of prophages to evolution and virulence of bacterial pathogens. Virulence, 4, 354–365.
29.
FriedmanCR, QuinnGC, KreiswirthBN, et al. (1997). Widespread dissemination of a drug-susceptible strain of Mycobacterium tuberculosis. J Infect Dis, 176, 478–484.
30.
Garcia de ViedmaD, LorenzoG, CardonaPJ, et al. (2005). Association between the infectivity of Mycobacterium tuberculosis strains and their efficiency for extrarespiratory infection. J Infect Dis, 192, 2059–2065.
31.
GlickmanMS, CoxJS, and JacobsWRJr. (2000). A novel mycolic acid cyclopropane synthetase is required for cording, persistence, and virulence of Mycobacterium tuberculosis. Mol Cell, 5, 717–727.
32.
GoldenMP, and VikramHR. (2005). Extrapulmonary tuberculosis: An overview. Am Fam Physician, 72, 1761–1768.
33.
GorenMB, CernichM, and BroklO. (1978). Some observations of mycobacterial acid-fastness. Am Rev Respir Dis, 118, 151–154.
34.
GrissaI, VergnaudG, and PourcelC. (2007). CRISPRFinder: A web tool to identify clustered regularly interspaced short palindromic repeats. Nucleic Acids Res, 35, W52–W57.
35.
GuS, ChenJ, DobosKM, BradburyEM, BelisleJT, and ChenX. (2003). Comprehensive proteomic profiling of the membrane constituents of a Mycobacterium tuberculosis strain. Mol Cell Proteomics, 2, 1284–1296.
36.
HasanZ, ZaidiI, JamilB, KhanMA, KanjiA, and HussainR. (2005). Elevated ex vivo monocyte chemotactic protein-1 (CCL2) in pulmonary as compared with extra-pulmonary tuberculosis. BMC Immunol, 6, 14.
37.
HusonDH, and ScornavaccaC. (2012). Dendroscope 3: An interactive tool for rooted phylogenetic trees and networks. Syst Biol, 61, 1061–1067.
38.
IsemanMD. (1993). Treatment of multidrug-resistant tuberculosis. N Engl J Med, 329, 784–791.
39.
KansalRG, Gomez-FloresR, and MehtaRT. (1998). Change in colony morphology influences the virulence as well as the biochemical properties of the Mycobacterium avium complex. Microb Pathog, 25, 203–214.
40.
KolattukudyPE, FernandesND, AzadAK, FitzmauriceAM, and SirakovaTD. (1997). Biochemistry and molecular genetics of cell-wall lipid biosynthesis in mycobacteria. Mol Microbiol, 24, 263–270.
41.
LiH, and DurbinR. (2009). Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics, 25, 1754–1760.
42.
MartinezA, TorelloS, and KolterR. (1999). Sliding motility in mycobacteria. J Bacteriol, 181, 7331–7338.
43.
MassiMN, BiatkoKT, HandayaniI, et al. (2017). Evaluation of rapid GeneXpert MTB/RIF method using DNA tissue specimens of vertebral bones in patients with suspected spondylitis TB. J Orthop, 14, 189–191.
44.
McKennaA, HannaM, BanksE, et al. (2010). The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res, 20, 1297–1303.
45.
MeghdadiH, KhosraviAD, GhadiriAA, SinaAH, and AlamiA. (2015). Detection of Mycobacterium tuberculosis in extrapulmonary biopsy samples using PCR targeting IS6110, rpoB, and nested-rpoB PCR cloning. Front Microbiol, 6, 675.
46.
MiH, HuangX, MuruganujanA, et al. (2017). PANTHER version 11: Expanded annotation data from gene ontology and reactome pathways, and data analysis tool enhancements. Nucleic Acids Res, 45, D183–D189.
47.
MiottoP, FortiF, AmbrosiA, et al. (2012). Genome-wide discovery of small RNAs in Mycobacterium tuberculosis. PLoS One, 7, e51950.
48.
NathavitharanaRR, CudahyPG, SchumacherSG, SteingartKR, PaiM, and DenkingerCM. (2017). Accuracy of line probe assays for the diagnosis of pulmonary and multidrug-resistant tuberculosis: A systematic review and meta-analysis. Eur Respir J. 49.
49.
OkiNO, Motsinger-ReifAA, AntasPR, LevyS, HollandSM, and SterlingTR. (2011). Novel human genetic variants associated with extrapulmonary tuberculosis: A pilot genome wide association study. BMC Res Notes, 4, 28.
50.
PengL, PeiJ, PangH, GuoY, LinL, and HuangR. (2014). Whole genome sequencing reveals a novel CRISPR system in industrial Clostridium acetobutylicum. J Ind Microbiol Biotechnol, 41, 1677–1685.
51.
PriceMN, DehalPS, and ArkinAP. (2010). FastTree 2—approximately maximum-likelihood trees for large alignments. PLoS One, 5, e9490.
SaeedM, IramS, HussainS, AhmedA, AkbarM, and AslamM. (2017). GeneXpert: A new tool for the rapid detection of rifampicin resistance in Mycobacterium tuberculosis. J Pak Med Assoc, 67, 270–274.
54.
SampsonSL. (2011). Mycobacterial PE/PPE proteins at the host-pathogen interface. Clin Dev Immunol, 2011, 497203.
55.
ShemyakinIG, StepanshinaVN, IvanovIY, et al. (2004). Characterization of drug-resistant isolates of Mycobacterium tuberculosis derived from Russian inmates. Int J Tuberc Lung Dis, 8, 1194–1203.
56.
SongT, and WaiSN. (2009). A novel sRNA that modulates virulence and environmental fitness of Vibrio cholerae. RNA Biol, 6, 254–258.
57.
SrilohasinP, ChaiprasertA, TokunagaK, et al. (2014). Genetic diversity and dynamic distribution of Mycobacterium tuberculosis isolates causing pulmonary and extrapulmonary tuberculosis in Thailand. J Clin Microbiol, 52, 4267–4274.
58.
SteinerA, StuckiD, CoscollaM, BorrellS, and GagneuxS. (2014). KvarQ: Targeted and direct variant calling from FastQ reads of bacterial genomes. BMC Genomics, 15, 881.
59.
ToungoussovaOS, CaugantDA, SandvenP, MariandyshevAO, and BjuneG. (2004). Impact of drug resistance on fitness of Mycobacterium tuberculosis strains of the W-Beijing genotype. FEMS Immunol Med Microbiol, 42, 281–290.
60.
TsoHW, IpWK, ChongWP, TamCM, ChiangAK, and LauYL. (2005). Association of interferon gamma and interleukin 10 genes with tuberculosis in Hong Kong Chinese. Genes Immun, 6, 358–363.
61.
van SoolingenD, HermansPW, de HaasPE, SollDR, and van EmbdenJD. (1991). Occurrence and stability of insertion sequences in Mycobacterium tuberculosis complex strains: Evaluation of an insertion sequence-dependent DNA polymorphism as a tool in the epidemiology of tuberculosis. J Clin Microbiol, 29, 2578–2586.
62.
van SoolingenD, JajouR, MulderA, and de NeelingH. (2016). Whole genome sequencing as the ultimate tool to diagnose tuberculosis. Int J Mycobacteriol, 5,(Suppl 1),S60–S61.
63.
WilkinsonRJ, LlewelynM, ToossiZ, et al. (2000). Influence of vitamin D deficiency and vitamin D receptor polymorphisms on tuberculosis among Gujarati Asians in west London: A case-control study. Lancet, 355, 618–621.
64.
WongKC, LeongWM, LawHK, et al. (2007). Molecular characterization of clinical isolates of Mycobacterium tuberculosis and their association with phenotypic virulence in human macrophages. Clin Vaccine Immunol, 14, 1279–1284.
65.
YeK, SchulzMH, LongQ, ApweilerR, and NingZ. (2009). Pindel: A pattern growth approach to detect break points of large deletions and medium sized insertions from paired-end short reads. Bioinformatics, 25, 2865–2871.
66.
ZhangH, LiD, ZhaoL, et al. (2013). Genome sequencing of 161 Mycobacterium tuberculosis isolates from China identifies genes and intergenic regions associated with drug resistance. Nat Genet, 45, 1255–1260.
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
