Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Whole genome sequencing, analyses of drug resistance-conferring mutations, and correlation with transmission of <italic>Mycobacterium tuberculosis</italic>
carrying <italic>katG</italic>
-S315T in Hanoi, Vietnam</title>
<author><name sortKey="Hang, Nguyen Thi Le" sort="Hang, Nguyen Thi Le" uniqKey="Hang N" first="Nguyen Thi Le" last="Hang">Nguyen Thi Le Hang</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 4691 4377</institution-id>
<institution-id institution-id-type="GRID">grid.414163.5</institution-id>
<institution>NCGM-BMH Medical Collaboration Center,</institution>
</institution-wrap>
Hanoi, Vietnam</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Hijikata, Minako" sort="Hijikata, Minako" uniqKey="Hijikata M" first="Minako" last="Hijikata">Minako Hijikata</name>
<affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 1545 6914</institution-id>
<institution-id institution-id-type="GRID">grid.419151.9</institution-id>
<institution>Department of Pathophysiology and Host Defense,</institution>
<institution>The Research Institute of Tuberculosis, JATA,</institution>
</institution-wrap>
Tokyo, Japan</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Maeda, Shinji" sort="Maeda, Shinji" uniqKey="Maeda S" first="Shinji" last="Maeda">Shinji Maeda</name>
<affiliation><nlm:aff id="Aff3"><institution-wrap><institution-id institution-id-type="GRID">grid.444700.3</institution-id>
<institution>Faculty of Pharmaceutical Sciences,</institution>
<institution>Hokkaido University of Science,</institution>
</institution-wrap>
Hokkaido, Japan</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Thuong, Pham Huu" sort="Thuong, Pham Huu" uniqKey="Thuong P" first="Pham Huu" last="Thuong">Pham Huu Thuong</name>
<affiliation><nlm:aff id="Aff4"><institution-wrap><institution-id institution-id-type="GRID">grid.470059.f</institution-id>
<institution>Hanoi Lung Hospital,</institution>
</institution-wrap>
Hanoi, Vietnam</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Ohashi, Jun" sort="Ohashi, Jun" uniqKey="Ohashi J" first="Jun" last="Ohashi">Jun Ohashi</name>
<affiliation><nlm:aff id="Aff5"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 2151 536X</institution-id>
<institution-id institution-id-type="GRID">grid.26999.3d</institution-id>
<institution>Department of Biological Sciences,</institution>
<institution>Graduate School of Science, The University of Tokyo,</institution>
</institution-wrap>
Tokyo, Japan</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Van Huan, Hoang" sort="Van Huan, Hoang" uniqKey="Van Huan H" first="Hoang" last="Van Huan">Hoang Van Huan</name>
<affiliation><nlm:aff id="Aff4"><institution-wrap><institution-id institution-id-type="GRID">grid.470059.f</institution-id>
<institution>Hanoi Lung Hospital,</institution>
</institution-wrap>
Hanoi, Vietnam</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Hoang, Nguyen Phuong" sort="Hoang, Nguyen Phuong" uniqKey="Hoang N" first="Nguyen Phuong" last="Hoang">Nguyen Phuong Hoang</name>
<affiliation><nlm:aff id="Aff6"><institution-wrap><institution-id institution-id-type="GRID">grid.470059.f</institution-id>
<institution>Department of Microbiology,</institution>
<institution>Hanoi Lung Hospital,</institution>
</institution-wrap>
Hanoi, Vietnam</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Miyabayashi, Akiko" sort="Miyabayashi, Akiko" uniqKey="Miyabayashi A" first="Akiko" last="Miyabayashi">Akiko Miyabayashi</name>
<affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 1545 6914</institution-id>
<institution-id institution-id-type="GRID">grid.419151.9</institution-id>
<institution>Department of Pathophysiology and Host Defense,</institution>
<institution>The Research Institute of Tuberculosis, JATA,</institution>
</institution-wrap>
Tokyo, Japan</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Cuong, Vu Cao" sort="Cuong, Vu Cao" uniqKey="Cuong V" first="Vu Cao" last="Cuong">Vu Cao Cuong</name>
<affiliation><nlm:aff id="Aff7">Hanoi Department of Health, Hanoi, Vietnam</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Seto, Shintaro" sort="Seto, Shintaro" uniqKey="Seto S" first="Shintaro" last="Seto">Shintaro Seto</name>
<affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 1545 6914</institution-id>
<institution-id institution-id-type="GRID">grid.419151.9</institution-id>
<institution>Department of Pathophysiology and Host Defense,</institution>
<institution>The Research Institute of Tuberculosis, JATA,</institution>
</institution-wrap>
Tokyo, Japan</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Van Hung, Nguyen" sort="Van Hung, Nguyen" uniqKey="Van Hung N" first="Nguyen" last="Van Hung">Nguyen Van Hung</name>
<affiliation><nlm:aff id="Aff8">Department of Microbiology, National Lung Hospital, Hanoi, Vietnam</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Keicho, Naoto" sort="Keicho, Naoto" uniqKey="Keicho N" first="Naoto" last="Keicho">Naoto Keicho</name>
<affiliation><nlm:aff id="Aff9"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 1545 6914</institution-id>
<institution-id institution-id-type="GRID">grid.419151.9</institution-id>
<institution>The Research Institute of Tuberculosis JATA,</institution>
</institution-wrap>
Tokyo, Japan</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff10"><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 0489 0290</institution-id>
<institution-id institution-id-type="GRID">grid.45203.30</institution-id>
<institution>National Center for Global Health and Medicine,</institution>
</institution-wrap>
Tokyo, Japan</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">31653940</idno>
<idno type="pmc">6814805</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC6814805</idno>
<idno type="RBID">PMC:6814805</idno>
<idno type="doi">10.1038/s41598-019-51812-7</idno>
<date when="2019">2019</date>
<idno type="wicri:Area/Pmc/Corpus">000462</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000462</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Whole genome sequencing, analyses of drug resistance-conferring mutations, and correlation with transmission of <italic>Mycobacterium tuberculosis</italic>
carrying <italic>katG</italic>
-S315T in Hanoi, Vietnam</title>
<author><name sortKey="Hang, Nguyen Thi Le" sort="Hang, Nguyen Thi Le" uniqKey="Hang N" first="Nguyen Thi Le" last="Hang">Nguyen Thi Le Hang</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 4691 4377</institution-id>
<institution-id institution-id-type="GRID">grid.414163.5</institution-id>
<institution>NCGM-BMH Medical Collaboration Center,</institution>
</institution-wrap>
Hanoi, Vietnam</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Hijikata, Minako" sort="Hijikata, Minako" uniqKey="Hijikata M" first="Minako" last="Hijikata">Minako Hijikata</name>
<affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 1545 6914</institution-id>
<institution-id institution-id-type="GRID">grid.419151.9</institution-id>
<institution>Department of Pathophysiology and Host Defense,</institution>
<institution>The Research Institute of Tuberculosis, JATA,</institution>
</institution-wrap>
Tokyo, Japan</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Maeda, Shinji" sort="Maeda, Shinji" uniqKey="Maeda S" first="Shinji" last="Maeda">Shinji Maeda</name>
<affiliation><nlm:aff id="Aff3"><institution-wrap><institution-id institution-id-type="GRID">grid.444700.3</institution-id>
<institution>Faculty of Pharmaceutical Sciences,</institution>
<institution>Hokkaido University of Science,</institution>
</institution-wrap>
Hokkaido, Japan</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Thuong, Pham Huu" sort="Thuong, Pham Huu" uniqKey="Thuong P" first="Pham Huu" last="Thuong">Pham Huu Thuong</name>
<affiliation><nlm:aff id="Aff4"><institution-wrap><institution-id institution-id-type="GRID">grid.470059.f</institution-id>
<institution>Hanoi Lung Hospital,</institution>
</institution-wrap>
Hanoi, Vietnam</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Ohashi, Jun" sort="Ohashi, Jun" uniqKey="Ohashi J" first="Jun" last="Ohashi">Jun Ohashi</name>
<affiliation><nlm:aff id="Aff5"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 2151 536X</institution-id>
<institution-id institution-id-type="GRID">grid.26999.3d</institution-id>
<institution>Department of Biological Sciences,</institution>
<institution>Graduate School of Science, The University of Tokyo,</institution>
</institution-wrap>
Tokyo, Japan</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Van Huan, Hoang" sort="Van Huan, Hoang" uniqKey="Van Huan H" first="Hoang" last="Van Huan">Hoang Van Huan</name>
<affiliation><nlm:aff id="Aff4"><institution-wrap><institution-id institution-id-type="GRID">grid.470059.f</institution-id>
<institution>Hanoi Lung Hospital,</institution>
</institution-wrap>
Hanoi, Vietnam</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Hoang, Nguyen Phuong" sort="Hoang, Nguyen Phuong" uniqKey="Hoang N" first="Nguyen Phuong" last="Hoang">Nguyen Phuong Hoang</name>
<affiliation><nlm:aff id="Aff6"><institution-wrap><institution-id institution-id-type="GRID">grid.470059.f</institution-id>
<institution>Department of Microbiology,</institution>
<institution>Hanoi Lung Hospital,</institution>
</institution-wrap>
Hanoi, Vietnam</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Miyabayashi, Akiko" sort="Miyabayashi, Akiko" uniqKey="Miyabayashi A" first="Akiko" last="Miyabayashi">Akiko Miyabayashi</name>
<affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 1545 6914</institution-id>
<institution-id institution-id-type="GRID">grid.419151.9</institution-id>
<institution>Department of Pathophysiology and Host Defense,</institution>
<institution>The Research Institute of Tuberculosis, JATA,</institution>
</institution-wrap>
Tokyo, Japan</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Cuong, Vu Cao" sort="Cuong, Vu Cao" uniqKey="Cuong V" first="Vu Cao" last="Cuong">Vu Cao Cuong</name>
<affiliation><nlm:aff id="Aff7">Hanoi Department of Health, Hanoi, Vietnam</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Seto, Shintaro" sort="Seto, Shintaro" uniqKey="Seto S" first="Shintaro" last="Seto">Shintaro Seto</name>
<affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 1545 6914</institution-id>
<institution-id institution-id-type="GRID">grid.419151.9</institution-id>
<institution>Department of Pathophysiology and Host Defense,</institution>
<institution>The Research Institute of Tuberculosis, JATA,</institution>
</institution-wrap>
Tokyo, Japan</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Van Hung, Nguyen" sort="Van Hung, Nguyen" uniqKey="Van Hung N" first="Nguyen" last="Van Hung">Nguyen Van Hung</name>
<affiliation><nlm:aff id="Aff8">Department of Microbiology, National Lung Hospital, Hanoi, Vietnam</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Keicho, Naoto" sort="Keicho, Naoto" uniqKey="Keicho N" first="Naoto" last="Keicho">Naoto Keicho</name>
<affiliation><nlm:aff id="Aff9"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 1545 6914</institution-id>
<institution-id institution-id-type="GRID">grid.419151.9</institution-id>
<institution>The Research Institute of Tuberculosis JATA,</institution>
</institution-wrap>
Tokyo, Japan</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff10"><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 0489 0290</institution-id>
<institution-id institution-id-type="GRID">grid.45203.30</institution-id>
<institution>National Center for Global Health and Medicine,</institution>
</institution-wrap>
Tokyo, Japan</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">Scientific Reports</title>
<idno type="eISSN">2045-2322</idno>
<imprint><date when="2019">2019</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><p id="Par1">Drug-resistant tuberculosis (TB) is a serious global problem, and pathogen factors involved in the transmission of isoniazid (INH)-resistant TB have not been fully investigated. We performed whole genome sequencing of 332 clinical <italic>Mycobacterium tuberculosis</italic>
(Mtb) isolates collected from patients newly diagnosed with smear-positive pulmonary TB in Hanoi, Vietnam. Using a bacterial genome-wide approach based on linear mixed models, we investigated the associations between 31-bp k-mers and clustered strains harboring <italic>katG</italic>
-S315T, a major INH-resistance mutation in the present cohort and in the second panel previously published in South Africa. Five statistically significant genes, namely, <italic>PPE18/19</italic>
, <italic>gid</italic>
, <italic>emrB</italic>
, <italic>Rv1588c</italic>
, and <italic>pncA</italic>
, were shared by the two panels. We further identified variants of the genes responsible for these k-mers, which are relevant to the spread of INH-resistant strains. Phylogenetic convergence test showed that variants relevant to <italic>PPE46/47</italic>
-like chimeric genes were significantly associated with the same phenotype in Hanoi. The associations were further confirmed after adjustment for the confounders. These findings suggest that genomic variations of the pathogen facilitate the expansion of INH-resistance TB, at least in part, and our study provides a new insight into the mechanisms by which drug-resistant Mtb maintains fitness and spreads in Asia and Africa.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Bastos, Ml" uniqKey="Bastos M">ML Bastos</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pooran, A" uniqKey="Pooran A">A Pooran</name>
</author>
<author><name sortKey="Pieterson, E" uniqKey="Pieterson E">E Pieterson</name>
</author>
<author><name sortKey="Davids, M" uniqKey="Davids M">M Davids</name>
</author>
<author><name sortKey="Theron, G" uniqKey="Theron G">G Theron</name>
</author>
<author><name sortKey="Dheda, K" uniqKey="Dheda K">K Dheda</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Muller, B" uniqKey="Muller B">B Muller</name>
</author>
<author><name sortKey="Borrell, S" uniqKey="Borrell S">S Borrell</name>
</author>
<author><name sortKey="Rose, G" uniqKey="Rose G">G Rose</name>
</author>
<author><name sortKey="Gagneux, S" uniqKey="Gagneux S">S Gagneux</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Manson, Al" uniqKey="Manson A">AL Manson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Menzies, Dick" uniqKey="Menzies D">Dick Menzies</name>
</author>
<author><name sortKey="Benedetti, Andrea" uniqKey="Benedetti A">Andrea Benedetti</name>
</author>
<author><name sortKey="Paydar, Anita" uniqKey="Paydar A">Anita Paydar</name>
</author>
<author><name sortKey="Royce, Sarah" uniqKey="Royce S">Sarah Royce</name>
</author>
<author><name sortKey="Pai, Madhukar" uniqKey="Pai M">Madhukar Pai</name>
</author>
<author><name sortKey="Burman, William" uniqKey="Burman W">William Burman</name>
</author>
<author><name sortKey="Vernon, Andrew" uniqKey="Vernon A">Andrew Vernon</name>
</author>
<author><name sortKey="Lienhardt, Christian" uniqKey="Lienhardt C">Christian Lienhardt</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Stagg, Hr" uniqKey="Stagg H">HR Stagg</name>
</author>
<author><name sortKey="Lipman, Mc" uniqKey="Lipman M">MC Lipman</name>
</author>
<author><name sortKey="Mchugh, Td" uniqKey="Mchugh T">TD McHugh</name>
</author>
<author><name sortKey="Jenkins, He" uniqKey="Jenkins H">HE Jenkins</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nebenzahl Guimaraes, H" uniqKey="Nebenzahl Guimaraes H">H Nebenzahl-Guimaraes</name>
</author>
<author><name sortKey="Jacobson, Kr" uniqKey="Jacobson K">KR Jacobson</name>
</author>
<author><name sortKey="Farhat, Mr" uniqKey="Farhat M">MR Farhat</name>
</author>
<author><name sortKey="Murray, Mb" uniqKey="Murray M">MB Murray</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Coll, F" uniqKey="Coll F">F Coll</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Farhat, Mr" uniqKey="Farhat M">MR Farhat</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nhung, Nv" uniqKey="Nhung N">NV Nhung</name>
</author>
<author><name sortKey="Hoa, Nb" uniqKey="Hoa N">NB Hoa</name>
</author>
<author><name sortKey="Sy, Dn" uniqKey="Sy D">DN Sy</name>
</author>
<author><name sortKey="Hennig, Cm" uniqKey="Hennig C">CM Hennig</name>
</author>
<author><name sortKey="Dean, As" uniqKey="Dean A">AS Dean</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hang, Nt" uniqKey="Hang N">NT Hang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Minh, Nn" uniqKey="Minh N">NN Minh</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nhu, Nt" uniqKey="Nhu N">NT Nhu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nguyen, Hq" uniqKey="Nguyen H">HQ Nguyen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Papaventsis, D" uniqKey="Papaventsis D">D Papaventsis</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zhang, H" uniqKey="Zhang H">H Zhang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Casali, N" uniqKey="Casali N">N Casali</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cohen, Ka" uniqKey="Cohen K">KA Cohen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Earle, Sg" uniqKey="Earle S">SG Earle</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Power, Ra" uniqKey="Power R">RA Power</name>
</author>
<author><name sortKey="Parkhill, J" uniqKey="Parkhill J">J Parkhill</name>
</author>
<author><name sortKey="De Oliveira, T" uniqKey="De Oliveira T">T de Oliveira</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Coll, F" uniqKey="Coll F">F Coll</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gagneux, S" uniqKey="Gagneux S">S Gagneux</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Choi, Y" uniqKey="Choi Y">Y Choi</name>
</author>
<author><name sortKey="Chan, Ap" uniqKey="Chan A">AP Chan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Vaser, R" uniqKey="Vaser R">R Vaser</name>
</author>
<author><name sortKey="Adusumalli, S" uniqKey="Adusumalli S">S Adusumalli</name>
</author>
<author><name sortKey="Leng, Sn" uniqKey="Leng S">SN Leng</name>
</author>
<author><name sortKey="Sikic, M" uniqKey="Sikic M">M Sikic</name>
</author>
<author><name sortKey="Ng, Pc" uniqKey="Ng P">PC Ng</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hicks, Nd" uniqKey="Hicks N">ND Hicks</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Huyen, Mn" uniqKey="Huyen M">MN Huyen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Seifert, M" uniqKey="Seifert M">M Seifert</name>
</author>
<author><name sortKey="Catanzaro, D" uniqKey="Catanzaro D">D Catanzaro</name>
</author>
<author><name sortKey="Catanzaro, A" uniqKey="Catanzaro A">A Catanzaro</name>
</author>
<author><name sortKey="Rodwell, Tc" uniqKey="Rodwell T">TC Rodwell</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ameeruddin, Nu" uniqKey="Ameeruddin N">NU Ameeruddin</name>
</author>
<author><name sortKey="Luke Elizabeth, H" uniqKey="Luke Elizabeth H">H Luke Elizabeth</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ano, H" uniqKey="Ano H">H Ano</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Unissa, An" uniqKey="Unissa A">AN Unissa</name>
</author>
<author><name sortKey="Selvakumar, N" uniqKey="Selvakumar N">N Selvakumar</name>
</author>
<author><name sortKey="Narayanan, S" uniqKey="Narayanan S">S Narayanan</name>
</author>
<author><name sortKey="Suganthi, C" uniqKey="Suganthi C">C Suganthi</name>
</author>
<author><name sortKey="Hanna, Le" uniqKey="Hanna L">LE Hanna</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Walker, Tm" uniqKey="Walker T">TM Walker</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mcevoy, Cr" uniqKey="Mcevoy C">CR McEvoy</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Karboul, A" uniqKey="Karboul A">A Karboul</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Fishbein, S" uniqKey="Fishbein S">S Fishbein</name>
</author>
<author><name sortKey="Van Wyk, N" uniqKey="Van Wyk N">N van Wyk</name>
</author>
<author><name sortKey="Warren, Rm" uniqKey="Warren R">RM Warren</name>
</author>
<author><name sortKey="Sampson, Sl" uniqKey="Sampson S">SL Sampson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nair, S" uniqKey="Nair S">S Nair</name>
</author>
<author><name sortKey="Pandey, Ad" uniqKey="Pandey A">AD Pandey</name>
</author>
<author><name sortKey="Mukhopadhyay, S" uniqKey="Mukhopadhyay S">S Mukhopadhyay</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bhat, Kh" uniqKey="Bhat K">KH Bhat</name>
</author>
<author><name sortKey="Ahmed, A" uniqKey="Ahmed A">A Ahmed</name>
</author>
<author><name sortKey="Kumar, S" uniqKey="Kumar S">S Kumar</name>
</author>
<author><name sortKey="Sharma, P" uniqKey="Sharma P">P Sharma</name>
</author>
<author><name sortKey="Mukhopadhyay, S" uniqKey="Mukhopadhyay S">S Mukhopadhyay</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Dillon, Dc" uniqKey="Dillon D">DC Dillon</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lewinsohn, Da" uniqKey="Lewinsohn D">DA Lewinsohn</name>
</author>
<author><name sortKey="Lines, Ra" uniqKey="Lines R">RA Lines</name>
</author>
<author><name sortKey="Lewinsohn, Dm" uniqKey="Lewinsohn D">DM Lewinsohn</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mortier, Mc" uniqKey="Mortier M">MC Mortier</name>
</author>
<author><name sortKey="Jongert, E" uniqKey="Jongert E">E Jongert</name>
</author>
<author><name sortKey="Mettens, P" uniqKey="Mettens P">P Mettens</name>
</author>
<author><name sortKey="Ruelle, Jl" uniqKey="Ruelle J">JL Ruelle</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Reed, S" uniqKey="Reed S">S Reed</name>
</author>
<author><name sortKey="Lobet, Y" uniqKey="Lobet Y">Y Lobet</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rodo, Mj" uniqKey="Rodo M">MJ Rodo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gaseitsiwe, S" uniqKey="Gaseitsiwe S">S Gaseitsiwe</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gaseitsiwe, S" uniqKey="Gaseitsiwe S">S Gaseitsiwe</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hebert, Am" uniqKey="Hebert A">AM Hebert</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Okamoto, S" uniqKey="Okamoto S">S Okamoto</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wong, Sy" uniqKey="Wong S">SY Wong</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Buu, Tn" uniqKey="Buu T">TN Buu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rodrigues, L" uniqKey="Rodrigues L">L Rodrigues</name>
</author>
<author><name sortKey="Parish, T" uniqKey="Parish T">T Parish</name>
</author>
<author><name sortKey="Balganesh, M" uniqKey="Balganesh M">M Balganesh</name>
</author>
<author><name sortKey="Ainsa, Ja" uniqKey="Ainsa J">JA Ainsa</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Louw, Ge" uniqKey="Louw G">GE Louw</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kanji, A" uniqKey="Kanji A">A Kanji</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pang, Y" uniqKey="Pang Y">Y Pang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gun, L" uniqKey="Gun L">L Gun</name>
</author>
<author><name sortKey="Yumiao, R" uniqKey="Yumiao R">R Yumiao</name>
</author>
<author><name sortKey="Haixian, P" uniqKey="Haixian P">P Haixian</name>
</author>
<author><name sortKey="Liang, Z" uniqKey="Liang Z">Z Liang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chionh, Yh" uniqKey="Chionh Y">YH Chionh</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Machado, D" uniqKey="Machado D">D Machado</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Garima, K" uniqKey="Garima K">K Garima</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gordon, Sv" uniqKey="Gordon S">SV Gordon</name>
</author>
<author><name sortKey="Heym, B" uniqKey="Heym B">B Heym</name>
</author>
<author><name sortKey="Parkhill, J" uniqKey="Parkhill J">J Parkhill</name>
</author>
<author><name sortKey="Barrell, B" uniqKey="Barrell B">B Barrell</name>
</author>
<author><name sortKey="Cole, St" uniqKey="Cole S">ST Cole</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chang, Kc" uniqKey="Chang K">KC Chang</name>
</author>
<author><name sortKey="Yew, Ww" uniqKey="Yew W">WW Yew</name>
</author>
<author><name sortKey="Zhang, Y" uniqKey="Zhang Y">Y Zhang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Miotto, P" uniqKey="Miotto P">P Miotto</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Xia, Q" uniqKey="Xia Q">Q Xia</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pang, Y" uniqKey="Pang Y">Y Pang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Marmiesse, M" uniqKey="Marmiesse M">M Marmiesse</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Unissa, An" uniqKey="Unissa A">AN Unissa</name>
</author>
<author><name sortKey="Subbian, S" uniqKey="Subbian S">S Subbian</name>
</author>
<author><name sortKey="Hanna, Le" uniqKey="Hanna L">LE Hanna</name>
</author>
<author><name sortKey="Selvakumar, N" uniqKey="Selvakumar N">N Selvakumar</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ng, Vh" uniqKey="Ng V">VH Ng</name>
</author>
<author><name sortKey="Cox, Js" uniqKey="Cox J">JS Cox</name>
</author>
<author><name sortKey="Sousa, Ao" uniqKey="Sousa A">AO Sousa</name>
</author>
<author><name sortKey="Macmicking, Jd" uniqKey="Macmicking J">JD MacMicking</name>
</author>
<author><name sortKey="Mckinney, Jd" uniqKey="Mckinney J">JD McKinney</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cui, Ze Jia" uniqKey="Cui Z">Ze-Jia Cui</name>
</author>
<author><name sortKey="Yang, Qing Yong" uniqKey="Yang Q">Qing-Yong Yang</name>
</author>
<author><name sortKey="Zhang, Hong Yu" uniqKey="Zhang H">Hong-Yu Zhang</name>
</author>
<author><name sortKey="Zhu, Qiang" uniqKey="Zhu Q">Qiang Zhu</name>
</author>
<author><name sortKey="Zhang, Qing Ye" uniqKey="Zhang Q">Qing-Ye Zhang</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Hung, Nv" uniqKey="Hung N">NV Hung</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Maeda, S" uniqKey="Maeda S">S Maeda</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hang, Nt" uniqKey="Hang N">NT Hang</name>
</author>
<author><name sortKey="Maeda, S" uniqKey="Maeda S">S Maeda</name>
</author>
<author><name sortKey="Keicho, N" uniqKey="Keicho N">N Keicho</name>
</author>
<author><name sortKey="Thuong, Ph" uniqKey="Thuong P">PH Thuong</name>
</author>
<author><name sortKey="Endo, H" uniqKey="Endo H">H Endo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hijikata, M" uniqKey="Hijikata M">M Hijikata</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chuang, Pc" uniqKey="Chuang P">PC Chuang</name>
</author>
<author><name sortKey="Chen, Hy" uniqKey="Chen H">HY Chen</name>
</author>
<author><name sortKey="Jou, R" uniqKey="Jou R">R Jou</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Stamatakis, A" uniqKey="Stamatakis A">A Stamatakis</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jaillard, M" uniqKey="Jaillard M">M Jaillard</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bankevich, A" uniqKey="Bankevich A">A Bankevich</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kajitani, R" uniqKey="Kajitani R">R Kajitani</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">Sci Rep</journal-id>
<journal-id journal-id-type="iso-abbrev">Sci Rep</journal-id>
<journal-title-group><journal-title>Scientific Reports</journal-title>
</journal-title-group>
<issn pub-type="epub">2045-2322</issn>
<publisher><publisher-name>Nature Publishing Group UK</publisher-name>
<publisher-loc>London</publisher-loc>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">31653940</article-id>
<article-id pub-id-type="pmc">6814805</article-id>
<article-id pub-id-type="publisher-id">51812</article-id>
<article-id pub-id-type="doi">10.1038/s41598-019-51812-7</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Article</subject>
</subj-group>
</article-categories>
<title-group><article-title>Whole genome sequencing, analyses of drug resistance-conferring mutations, and correlation with transmission of <italic>Mycobacterium tuberculosis</italic>
carrying <italic>katG</italic>
-S315T in Hanoi, Vietnam</article-title>
</title-group>
<contrib-group><contrib contrib-type="author" equal-contrib="yes"><name><surname>Hang</surname>
<given-names>Nguyen Thi Le</given-names>
</name>
<xref ref-type="aff" rid="Aff1">1</xref>
</contrib>
<contrib contrib-type="author" equal-contrib="yes"><name><surname>Hijikata</surname>
<given-names>Minako</given-names>
</name>
<xref ref-type="aff" rid="Aff2">2</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Maeda</surname>
<given-names>Shinji</given-names>
</name>
<xref ref-type="aff" rid="Aff3">3</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Thuong</surname>
<given-names>Pham Huu</given-names>
</name>
<xref ref-type="aff" rid="Aff4">4</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Ohashi</surname>
<given-names>Jun</given-names>
</name>
<xref ref-type="aff" rid="Aff5">5</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Van Huan</surname>
<given-names>Hoang</given-names>
</name>
<xref ref-type="aff" rid="Aff4">4</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Hoang</surname>
<given-names>Nguyen Phuong</given-names>
</name>
<xref ref-type="aff" rid="Aff6">6</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Miyabayashi</surname>
<given-names>Akiko</given-names>
</name>
<xref ref-type="aff" rid="Aff2">2</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Cuong</surname>
<given-names>Vu Cao</given-names>
</name>
<xref ref-type="aff" rid="Aff7">7</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Seto</surname>
<given-names>Shintaro</given-names>
</name>
<xref ref-type="aff" rid="Aff2">2</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Van Hung</surname>
<given-names>Nguyen</given-names>
</name>
<xref ref-type="aff" rid="Aff8">8</xref>
</contrib>
<contrib contrib-type="author" corresp="yes"><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0001-8232-9448</contrib-id>
<name><surname>Keicho</surname>
<given-names>Naoto</given-names>
</name>
<address><email>nkeicho-tky@umin.ac.jp</email>
</address>
<xref ref-type="aff" rid="Aff9">9</xref>
<xref ref-type="aff" rid="Aff10">10</xref>
</contrib>
<aff id="Aff1"><label>1</label>
<institution-wrap><institution-id institution-id-type="ISNI">0000 0004 4691 4377</institution-id>
<institution-id institution-id-type="GRID">grid.414163.5</institution-id>
<institution>NCGM-BMH Medical Collaboration Center,</institution>
</institution-wrap>
Hanoi, Vietnam</aff>
<aff id="Aff2"><label>2</label>
<institution-wrap><institution-id institution-id-type="ISNI">0000 0001 1545 6914</institution-id>
<institution-id institution-id-type="GRID">grid.419151.9</institution-id>
<institution>Department of Pathophysiology and Host Defense,</institution>
<institution>The Research Institute of Tuberculosis, JATA,</institution>
</institution-wrap>
Tokyo, Japan</aff>
<aff id="Aff3"><label>3</label>
<institution-wrap><institution-id institution-id-type="GRID">grid.444700.3</institution-id>
<institution>Faculty of Pharmaceutical Sciences,</institution>
<institution>Hokkaido University of Science,</institution>
</institution-wrap>
Hokkaido, Japan</aff>
<aff id="Aff4"><label>4</label>
<institution-wrap><institution-id institution-id-type="GRID">grid.470059.f</institution-id>
<institution>Hanoi Lung Hospital,</institution>
</institution-wrap>
Hanoi, Vietnam</aff>
<aff id="Aff5"><label>5</label>
<institution-wrap><institution-id institution-id-type="ISNI">0000 0001 2151 536X</institution-id>
<institution-id institution-id-type="GRID">grid.26999.3d</institution-id>
<institution>Department of Biological Sciences,</institution>
<institution>Graduate School of Science, The University of Tokyo,</institution>
</institution-wrap>
Tokyo, Japan</aff>
<aff id="Aff6"><label>6</label>
<institution-wrap><institution-id institution-id-type="GRID">grid.470059.f</institution-id>
<institution>Department of Microbiology,</institution>
<institution>Hanoi Lung Hospital,</institution>
</institution-wrap>
Hanoi, Vietnam</aff>
<aff id="Aff7"><label>7</label>
Hanoi Department of Health, Hanoi, Vietnam</aff>
<aff id="Aff8"><label>8</label>
Department of Microbiology, National Lung Hospital, Hanoi, Vietnam</aff>
<aff id="Aff9"><label>9</label>
<institution-wrap><institution-id institution-id-type="ISNI">0000 0001 1545 6914</institution-id>
<institution-id institution-id-type="GRID">grid.419151.9</institution-id>
<institution>The Research Institute of Tuberculosis JATA,</institution>
</institution-wrap>
Tokyo, Japan</aff>
<aff id="Aff10"><label>10</label>
<institution-wrap><institution-id institution-id-type="ISNI">0000 0004 0489 0290</institution-id>
<institution-id institution-id-type="GRID">grid.45203.30</institution-id>
<institution>National Center for Global Health and Medicine,</institution>
</institution-wrap>
Tokyo, Japan</aff>
</contrib-group>
<pub-date pub-type="epub"><day>25</day>
<month>10</month>
<year>2019</year>
</pub-date>
<pub-date pub-type="pmc-release"><day>25</day>
<month>10</month>
<year>2019</year>
</pub-date>
<pub-date pub-type="collection"><year>2019</year>
</pub-date>
<volume>9</volume>
<elocation-id>15354</elocation-id>
<history><date date-type="received"><day>26</day>
<month>7</month>
<year>2019</year>
</date>
<date date-type="accepted"><day>8</day>
<month>10</month>
<year>2019</year>
</date>
</history>
<permissions><copyright-statement>© The Author(s) 2019</copyright-statement>
<license license-type="OpenAccess"><license-p><bold>Open Access</bold>
This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">http://creativecommons.org/licenses/by/4.0/</ext-link>
.</license-p>
</license>
</permissions>
<abstract id="Abs1"><p id="Par1">Drug-resistant tuberculosis (TB) is a serious global problem, and pathogen factors involved in the transmission of isoniazid (INH)-resistant TB have not been fully investigated. We performed whole genome sequencing of 332 clinical <italic>Mycobacterium tuberculosis</italic>
(Mtb) isolates collected from patients newly diagnosed with smear-positive pulmonary TB in Hanoi, Vietnam. Using a bacterial genome-wide approach based on linear mixed models, we investigated the associations between 31-bp k-mers and clustered strains harboring <italic>katG</italic>
-S315T, a major INH-resistance mutation in the present cohort and in the second panel previously published in South Africa. Five statistically significant genes, namely, <italic>PPE18/19</italic>
, <italic>gid</italic>
, <italic>emrB</italic>
, <italic>Rv1588c</italic>
, and <italic>pncA</italic>
, were shared by the two panels. We further identified variants of the genes responsible for these k-mers, which are relevant to the spread of INH-resistant strains. Phylogenetic convergence test showed that variants relevant to <italic>PPE46/47</italic>
-like chimeric genes were significantly associated with the same phenotype in Hanoi. The associations were further confirmed after adjustment for the confounders. These findings suggest that genomic variations of the pathogen facilitate the expansion of INH-resistance TB, at least in part, and our study provides a new insight into the mechanisms by which drug-resistant Mtb maintains fitness and spreads in Asia and Africa.</p>
</abstract>
<kwd-group kwd-group-type="npg-subject"><title>Subject terms</title>
<kwd>Bacterial genes</kwd>
<kwd>Tuberculosis</kwd>
</kwd-group>
<funding-group><award-group><funding-source><institution>This research was supported by AMED under Grant Number JP19fm0108001 (Japan Initiative for Global Research Network on Infectious Diseases; J-GRID).</institution>
</funding-source>
</award-group>
</funding-group>
<custom-meta-group><custom-meta><meta-name>issue-copyright-statement</meta-name>
<meta-value>© The Author(s) 2019</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body><sec id="Sec1" sec-type="introduction"><title>Introduction</title>
<p id="Par2">In 2017, 558,000 tuberculosis (TB) cases of multidrug resistance (MDR) or rifampicin (RMP) resistance were estimated globally<sup><xref ref-type="bibr" rid="CR1">1</xref>
</sup>
. Resistance to anti-TB drugs increases the burden of TB because the treatment of drug-resistant TB is generally prolonged and costly, while the outcome is relatively poor<sup><xref ref-type="bibr" rid="CR2">2</xref>
,<xref ref-type="bibr" rid="CR3">3</xref>
</sup>
. Isoniazid (INH) resistance without concurrent RMP resistance accounts for 7.1% of the new TB cases in the world<sup><xref ref-type="bibr" rid="CR1">1</xref>
</sup>
. The spread of INH resistance serves as a reservoir of more combined drug resistance<sup><xref ref-type="bibr" rid="CR4">4</xref>
,<xref ref-type="bibr" rid="CR5">5</xref>
</sup>
; empirical treatment after rapid genetic assessment of RMP resistance alone increases the future risk for developing MDR-TB or extensively drug-resistant (XDR)-TB<sup><xref ref-type="bibr" rid="CR6">6</xref>
,<xref ref-type="bibr" rid="CR7">7</xref>
</sup>
. INH-containing preventive therapy is also affected heavily when the INH-resistance rate is high in a target population.</p>
<p id="Par3">Extensive studies have been performed to identify the genetic mutations responsible for <italic>Mycobacterium tuberculosis</italic>
(Mtb) drug resistance, and comprehensive lists of drug-resistance-conferring mutations have been provided<sup><xref ref-type="bibr" rid="CR8">8</xref>
,<xref ref-type="bibr" rid="CR9">9</xref>
</sup>
. However, unidentified genetic variants, including compensatory mutations, may also assist or facilitate the transmission of drug-resistant TB without reducing the fitness of the bacilli<sup><xref ref-type="bibr" rid="CR10">10</xref>
</sup>
. Investigation of this mechanism would contribute to the effective control of drug-resistant TB, and the accumulation of molecular epidemiological data in many areas around the world would deepen the understanding of the dynamics.</p>
<p id="Par4">Vietnam is one of the 30 TB high-burden countries, with approximately 124,000 incidence cases reported in 2017<sup><xref ref-type="bibr" rid="CR1">1</xref>
</sup>
. According to the national drug resistance survey conducted in 2011, among the new cases, resistance to any anti-TB drug accounted for 32.7% of the cases, and the proportion of INH-resistant TB reached 18.9%<sup><xref ref-type="bibr" rid="CR11">11</xref>
</sup>
. In our previous study cohort of 489 newly diagnosed patients in a city area, INH resistance was observed in 28.2% of the patients<sup><xref ref-type="bibr" rid="CR12">12</xref>
</sup>
, while RMP resistance remained in 4.9% of the patients. The predominant genetic mutations <italic>katG-</italic>
S315T for INH and <italic>rpoB-</italic>
S450L for RMP that confer drug resistance in Vietnam are similar to those reported in other Asian countries<sup><xref ref-type="bibr" rid="CR13">13</xref>
–<xref ref-type="bibr" rid="CR15">15</xref>
</sup>
. However, pathogen factors correlated with the transmission of INH-resistant TB have not been fully investigated.</p>
<p id="Par5">Whole genome sequencing (WGS) has been recently used globally, offering new opportunities in the management of drug-resistant TB, since it can provide a huge amount of information, including genetic variants that are relevant to drug resistance throughout the genome<sup><xref ref-type="bibr" rid="CR16">16</xref>
</sup>
. WGS data also offer critical insights into the dynamics of TB endemics, transmission route, and the evolutionary patterns of genomic mutations<sup><xref ref-type="bibr" rid="CR10">10</xref>
,<xref ref-type="bibr" rid="CR17">17</xref>
–<xref ref-type="bibr" rid="CR19">19</xref>
</sup>
. Recently, bacterial genome-wide association studies (GWAS) controlling for population structure have also been performed for identifying the genes or genetic variants relevant to the TB phenotype, including drug resistance, by analyzing all single-nucleotide polymorphisms, small and large insertions/deletions (indels), or k-mers obtained from massive short-read data from next-generation sequencers (NGS)<sup><xref ref-type="bibr" rid="CR20">20</xref>
,<xref ref-type="bibr" rid="CR21">21</xref>
</sup>
.</p>
<p id="Par6">In this study, we investigated drug resistance-conferring mutations carried by the clinical Mtb isolates from patients newly diagnosed with smear-positive pulmonary TB in Hanoi, Vietnam, by using WGS with a bacterial GWAS approach incorporating linear mixed models (LMMs). We then identified the genetic variants that may be relevant to the success in an extensive spread of INH-resistant strains, in reference to a previously published cohort study in KwaZulu-Natal, South Africa, as the second panel<sup><xref ref-type="bibr" rid="CR19">19</xref>
</sup>
.</p>
</sec>
<sec id="Sec2" sec-type="results"><title>Results</title>
<sec id="Sec3"><title>Prevalence of known drug resistance-conferring mutations</title>
<p id="Par7">Among the 332 Hanoi samples analyzed with WGS, known mutations, which have been registered in the TBProfiler’s mutation database<sup><xref ref-type="bibr" rid="CR9">9</xref>
,<xref ref-type="bibr" rid="CR22">22</xref>
</sup>
, accounted for 80 (90.9%) of the 88 isolates with phenotypic INH resistance, 12 (92.3%) of the 13 isolates with RMP resistance, 63 (73.3%) of the 86 isolates with streptomycin (SM) resistance, 5 (100.0%) of the 5 isolates with ethambutol resistance, and 5 (55.6%) of the 9 isolates with pyrazinamide (PZA) resistance (Table <xref rid="Tab1" ref-type="table">1</xref>
). The most prevalent drug-resistance-conferring mutations were <italic>katG-</italic>
S315T (26.2%) to INH, <italic>rpsL-</italic>
K43R (13.3%) to SM, and <italic>rpoB-</italic>
S450L (2.7%) to RMP. Mutations conferring resistance to the second-line drugs were <italic>rrs-</italic>
A514C (4.5%) to amikacin and those in the <italic>fabG1</italic>
-promoter (2.7%) to ethionamide. Others were rare mutations such as <italic>ethR-</italic>
F110L to ethionamide, and <italic>gyrA-</italic>
A90V and <italic>gyrA-</italic>
D94G to fluoroquinolone (0.6%, 0.3%, and 0.6%, respectively). Among the 87 strains harboring <italic>katG-</italic>
S315T mutations, 54 (62.1%) had at least one other drug-resistance-conferring mutation. Co-occurrence of <italic>katG</italic>
-S315T and <italic>rpsL</italic>
-K43R with or without other known mutations was observed most frequently (31 of the 87 isolates; 35.6%).<table-wrap id="Tab1"><label>Table 1</label>
<caption><p>Frequencies of known mutations conferring resistance to first-line drugs among all isolates, phenotypically resistant isolates, and susceptible isolates in the Hanoi sample set (N = 332).</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th>Drug resistance-conferring mutations</th>
<th>Study population<break></break>
n (%)</th>
<th>Phenotypically resistant isolates<break></break>
n (%)</th>
<th>Phenotypically susceptible isolates<break></break>
n (%)</th>
</tr>
</thead>
<tbody><tr><td><bold>Isoniazid</bold>
</td>
<td>N = 332</td>
<td>N = 88</td>
<td>N = 244</td>
</tr>
<tr><td>(−)</td>
<td>233 (70.2)</td>
<td>8 (9.1)</td>
<td>225 (92.2)</td>
</tr>
<tr><td><italic>katG</italic>
(S315T)</td>
<td>87 (26.2)</td>
<td>75 (85.3)</td>
<td>12 (4.9)</td>
</tr>
<tr><td>Rv1482c<italic>-fabG1</italic>
(C-15T)</td>
<td>7 (2.1)</td>
<td>3 (3.4)</td>
<td>4 (1.6)</td>
</tr>
<tr><td>Rv1482c<italic>-fabG1</italic>
(C-15T), <italic>inhA</italic>
(S94A)</td>
<td>1 (0.3)</td>
<td>1 (1.1)</td>
<td>0 (0.0)</td>
</tr>
<tr><td>Rv1482c<italic>-fabG1</italic>
(C-15T), <italic>katG</italic>
(S315N)</td>
<td>1 (0.3)</td>
<td>1 (1.1)</td>
<td>0 (0.0)</td>
</tr>
<tr><td><italic>katG</italic>
(G234R)</td>
<td>1 (0.3)</td>
<td>0 (0.0)</td>
<td>1 (0.4)</td>
</tr>
<tr><td><italic>katG</italic>
(S315G)</td>
<td>1 (0.3)</td>
<td>0 (0.0)</td>
<td>1 (0.4)</td>
</tr>
<tr><td><italic>katG</italic>
(W191R)</td>
<td>1 (0.3)</td>
<td>0 (0.0)</td>
<td>1 (0.4)</td>
</tr>
<tr><td><bold>Streptomycin</bold>
</td>
<td>N = 332</td>
<td>N = 86</td>
<td>N = 246</td>
</tr>
<tr><td>(−)</td>
<td>259 (78.0)</td>
<td>23 (26.6)</td>
<td>236 (96.0)</td>
</tr>
<tr><td><italic>rpsL</italic>
(K43R)</td>
<td>44 (13.3)</td>
<td>42 (48.8)</td>
<td>2 (0.8)</td>
</tr>
<tr><td><italic>rrs</italic>
(A514C)</td>
<td>14 (4.2)</td>
<td>11 (12.8)</td>
<td>3 (1.2)</td>
</tr>
<tr><td><italic>rpsL</italic>
(K88R)</td>
<td>6 (1.8)</td>
<td>6 (7.0)</td>
<td>0 (0.0)</td>
</tr>
<tr><td><italic>gid</italic>
(L16R), <italic>rrs</italic>
(A514C)</td>
<td>1 (0.3)</td>
<td>1 (1.2)</td>
<td>0 (0.0)</td>
</tr>
<tr><td><italic>rpsL</italic>
(K43R), <italic>rrs</italic>
(A1401G)</td>
<td>1 (0.3)</td>
<td>1 (1.2)</td>
<td>0 (0.0)</td>
</tr>
<tr><td><italic>rrs</italic>
(A907T)</td>
<td>1 (0.3)</td>
<td>1 (1.2)</td>
<td>0 (0.0)</td>
</tr>
<tr><td><italic>rrs</italic>
(C517T)</td>
<td>2 (0.6)</td>
<td>1 (1.2)</td>
<td>1 (0.4)</td>
</tr>
<tr><td><italic>rrs</italic>
(T1239C)</td>
<td>1 (0.3)</td>
<td>0 (0.0)</td>
<td>1 (0.4)</td>
</tr>
<tr><td><italic>rpsL</italic>
(K88T)</td>
<td>1 (0.3)</td>
<td>0 (0.0)</td>
<td>1 (0.4)</td>
</tr>
<tr><td><italic>gid</italic>
(L16R)</td>
<td>2 (0.6)</td>
<td>0 (0.0)</td>
<td>2 (0.8)</td>
</tr>
<tr><td><bold>Rifampicin</bold>
</td>
<td>N = 332</td>
<td>N = 13</td>
<td>N = 319</td>
</tr>
<tr><td>(−)</td>
<td>310 (93.4)</td>
<td>1 (7.7)</td>
<td>309 (96.9)</td>
</tr>
<tr><td><italic>rpoB</italic>
(S450L)</td>
<td>9 (2.7)</td>
<td>4 (30.8)</td>
<td>5 (1.6)</td>
</tr>
<tr><td><italic>rpoB</italic>
(H445D)</td>
<td>4 (1.2)</td>
<td>4 (30.8)</td>
<td>0 (0.0)</td>
</tr>
<tr><td><italic>rpoB</italic>
(H445L)</td>
<td>2 (0.6)</td>
<td>2 (15.4)</td>
<td>0 (0.0)</td>
</tr>
<tr><td><italic>rpoB</italic>
(K446Q)</td>
<td>1 (0.3)</td>
<td>1 (7.7)</td>
<td>0 (0.0)</td>
</tr>
<tr><td><italic>rpoB</italic>
(Q432K)</td>
<td>1 (0.3)</td>
<td>1 (7.7)</td>
<td>0 (0.0)</td>
</tr>
<tr><td><italic>rpoB</italic>
(L452P)</td>
<td>1 (0.3)</td>
<td>0 (0.0)</td>
<td>1 (0.3)</td>
</tr>
<tr><td><italic>rposB</italic>
(S450L), <italic>rpoB</italic>
(E761D)</td>
<td>1 (0.3)</td>
<td>0 (0.0)</td>
<td>1 (0.3)</td>
</tr>
<tr><td><italic>rpoB</italic>
(L430P)</td>
<td>3 (0.9)</td>
<td>0 (0.0)</td>
<td>3 (0.9)</td>
</tr>
<tr><td><bold>Ethambutol</bold>
</td>
<td>N = 332</td>
<td>N = 5</td>
<td>N = 327</td>
</tr>
<tr><td>(−)</td>
<td>303 (91.3)</td>
<td>0 (0.0)</td>
<td>303 (92.7)</td>
</tr>
<tr><td><italic>embB</italic>
(M306V)</td>
<td>6 (1.8)</td>
<td>3 (60.0)</td>
<td>3 (0.9)</td>
</tr>
<tr><td><italic>embB</italic>
(L370R)</td>
<td>13 (3.9)</td>
<td>1 (20.0)</td>
<td>12 (3.7)</td>
</tr>
<tr><td><italic>embC-embA</italic>
(C-12T)</td>
<td>1 (0.3)</td>
<td>1 (20.0)</td>
<td>0 (0.0)</td>
</tr>
<tr><td><italic>embB</italic>
(D354A)</td>
<td>1 (0.3)</td>
<td>0 (0.0)</td>
<td>1 (0.3)</td>
</tr>
<tr><td><italic>embB</italic>
(G406A)</td>
<td>2 (0.6)</td>
<td>0 (0.0)</td>
<td>2 (0.6)</td>
</tr>
<tr><td><italic>embB</italic>
(G406D)</td>
<td>1 (0.3)</td>
<td>0 (0.0)</td>
<td>1 (0.3)</td>
</tr>
<tr><td><italic>embB</italic>
(M306I)</td>
<td>4 (1.2)</td>
<td>0 (0.0)</td>
<td>4 (1.2)</td>
</tr>
<tr><td><italic>embB</italic>
(Q497R)</td>
<td>1 (0.3)</td>
<td>0 (0.0)</td>
<td>1 (0.3)</td>
</tr>
<tr><td><bold>Pyrazinamide</bold>
</td>
<td>N = 332</td>
<td>N = 9</td>
<td>N = 323</td>
</tr>
<tr><td>(−)</td>
<td>319 (96.1)</td>
<td>4 (44.4)</td>
<td>315 (97.5)</td>
</tr>
<tr><td><italic>pncA</italic>
(Q141P)</td>
<td>3 (0.9)</td>
<td>2 (22.2)</td>
<td>1 (0.3)</td>
</tr>
<tr><td><italic>pncA</italic>
(L159R)</td>
<td>1 (0.3)</td>
<td>1 (11.1)</td>
<td>0 (0.0)</td>
</tr>
<tr><td><italic>pncA</italic>
(S59P)</td>
<td>1 (0.3)</td>
<td>1 (11.1)</td>
<td>0 (0.0)</td>
</tr>
<tr><td><italic>pncA</italic>
-Rv2044c (T-11C)</td>
<td>3 (0.9)</td>
<td>1 (11.1)</td>
<td>2 (0.6)</td>
</tr>
<tr><td><italic>pncA</italic>
(A46V)</td>
<td>1 (0.3)</td>
<td>0 (0.0)</td>
<td>1 (0.3)</td>
</tr>
<tr><td><italic>pncA</italic>
(I31S)</td>
<td>1 (0.3)</td>
<td>0 (0.0)</td>
<td>1 (0.3)</td>
</tr>
<tr><td><italic>pncA</italic>
(Q10P)</td>
<td>1 (0.3)</td>
<td>0 (0.0)</td>
<td>1 (0.3)</td>
</tr>
<tr><td><italic>pncA</italic>
(Q10R)</td>
<td>1 (0.3)</td>
<td>0 (0.0)</td>
<td>1 (0.3)</td>
</tr>
<tr><td><italic>pncA</italic>
(S104R)</td>
<td>1 (0.3)</td>
<td>0 (0.0)</td>
<td>1 (0.3)</td>
</tr>
</tbody>
</table>
</table-wrap>
</p>
<p id="Par8">Other than SNVs, in one isolate, a 353-bp deletion was found in the <italic>pncA</italic>
region after the screening of zero or low-depth areas, which also covered the <italic>pncA</italic>
promoter and a part of the Rv2044c nearby. One sample had an in-frame 6-bp deletion (M434-D435) in the <italic>rpoB</italic>
gene; this strain was not resistant to RMP. Eight isolates harboring 1-bp deletion in the <italic>gid</italic>
gene, and two other isolates carrying 1-bp deletion of the <italic>ethA</italic>
gene were also found (Supplementary Fig. <xref rid="MOESM1" ref-type="media">S1</xref>
).</p>
</sec>
<sec id="Sec4"><title>Mtb lineages/sublineages, drug-resistance-conferring mutations, and genetic clusters</title>
<p id="Par9">The lineage (L)2 East Asian—mostly the Beijing strains—possessed the mutations that conferred drug resistance most frequently (52.6%); L1 (Indo-Oceanic; 23.2%) carried less and then L4 (Euro-American; 18.3%) the least (P < 0.0001) in Hanoi. The proportion of the above drug resistance was the highest in ancient Beijing strains (59.8%), followed by modern Beijing strains (40.9%), and then by non-Beijing strains (20.8%) (P < 0.0001).</p>
<p id="Par10">Figure <xref rid="Fig1" ref-type="fig">1</xref>
shows the distribution of the strains harboring mutations of interest. Of note, <italic>katG</italic>
-S315T, known as a major INH-resistance mutation in the world<sup><xref ref-type="bibr" rid="CR23">23</xref>
</sup>
, was reported with a frequency of 26.2% in Hanoi, which accounted for 85.3% of the primary INH-resistance, and it was distributed unevenly among the strains like scattered islands. This S315T mutation, regardless of whether it occurred alone or in combination with other mutations, was more frequently observed in the L2 and L1 strains than in L4 (P < 0.0001) and also more in ancient Beijing than in modern Beijing and non-Beijing strains (59.0%, 40.9%, and 19.4%, respectively, P < 0.001). <italic>katG</italic>
-S315T was often carried by ancient Beijing strains, whereas clustered strains defined by pairwise SNV differences of <6 alone were observed more frequently among modern Beijing strains than among ancient and non-Beijing strains (24.2%, 18.9%, and 9.7%, respectively, P = 0.016).<fig id="Fig1"><label>Figure 1</label>
<caption><p>Distribution of drug resistant-relevant mutations and their presence in correlation with Mtb lineages and clustering. Clustered: defined by pairwise SNV difference <6 SNVs. Different colors indicate different mutations.</p>
</caption>
<graphic xlink:href="41598_2019_51812_Fig1_HTML" id="d29e1458"></graphic>
</fig>
</p>
</sec>
<sec id="Sec5"><title>Framework of bacterial GWAS analysis</title>
<p id="Par11">To investigate the pathogen factors involved in the wide spread of INH-resistant strains in Hanoi, we used a combination of a representative INH resistance-conferring mutation <italic>katG</italic>
-S315T (S315T[+]) with genetic cluster (cluster[+]) <6 SNV differences as a “phenotype” or objective variable for GWAS. <italic>katG</italic>
-S315T and SNV-based genetic clusters alone served as phenotypes for comparison. Initially 31-bp k-mers throughout the genome were set as the “genotype” or explanatory variables, because the presence or absence of such k-mers provides clues for identifying the SNVs, indels, or even structural variants without a standard reference sequence for mapping reads<sup><xref ref-type="bibr" rid="CR20">20</xref>
</sup>
. When candidate genes were obtained from the k-mer analysis, individual variants, including SNVs and indels, were further analyzed for confirmation using other variant-based platforms.</p>
</sec>
<sec id="Sec6"><title>Association between phenotypes and 31-bp k-mers throughout the genome</title>
<p id="Par12">When k-mer GWAS was performed for our Hanoi samples, the number of genes harboring significant k-mers associated with cluster[+] alone and S315T[+] alone were 7 and 48, respectively, whereas the combined phenotype, cluster[+]/S315T[+] extracted 403 statistically significant genes, controlling for the population structure using LMM (Table <xref rid="Tab2" ref-type="table">2</xref>
). For comparison, 337 samples from the South African KwaZulu-Natal study were analyzed in a similar way as for the combined phenotype (Table <xref rid="Tab2" ref-type="table">2</xref>
), and 14 genes were eventually extracted in common (Fig. <xref rid="Fig2" ref-type="fig">2a</xref>
). Of these, eight (<italic>Rv1148c</italic>
, <italic>mmpL6, PE_PGRS21</italic>
, <italic>PE_PGRS53</italic>
, <italic>pks12, rpoB</italic>
, <italic>Rv2090</italic>
, and <italic>wag22</italic>
) were excluded from further analysis because the presence or absence of k-mers appeared only in one sample of either case or control group. <italic>PE_PGRS10</italic>
was further excluded since the presence of k-mers was not further confirmed by BLAST search for read sequences from both Hanoi and KwaZulu-Natal samples.<table-wrap id="Tab2"><label>Table 2</label>
<caption><p>Genes harboring k-mers significantly associated with different phenotypes (from (1) to (6)) in Hanoi and KwaZulu-Natal studies.</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th colspan="6">Phenotype</th>
<th rowspan="2">Gene name</th>
<th rowspan="2">Number of genes</th>
</tr>
<tr><th>(1)</th>
<th>(2)</th>
<th>(3)</th>
<th>(4)</th>
<th>(5)</th>
<th>(6)</th>
</tr>
</thead>
<tbody><tr><td>No</td>
<td>Yes</td>
<td><bold>Yes</bold>
</td>
<td>Yes</td>
<td>Yes</td>
<td><bold>Yes</bold>
</td>
<td><bold><italic>gid</italic>
</bold>
</td>
<td><bold>1</bold>
</td>
</tr>
<tr><td>No</td>
<td>Yes</td>
<td>No</td>
<td>Yes</td>
<td>No</td>
<td>Yes</td>
<td><italic>rrs</italic>
</td>
<td>1</td>
</tr>
<tr><td>No</td>
<td>Yes</td>
<td>Yes</td>
<td>No</td>
<td>Yes</td>
<td>No</td>
<td><italic>katG</italic>
</td>
<td>1</td>
</tr>
<tr><td>No</td>
<td>No</td>
<td><bold>Yes</bold>
</td>
<td>Yes</td>
<td>Yes</td>
<td><bold>Yes</bold>
</td>
<td><bold><italic>PPE18, Rv1148c</italic>
</bold>
</td>
<td><bold>2</bold>
</td>
</tr>
<tr><td>Yes</td>
<td>No</td>
<td>Yes</td>
<td>No</td>
<td>No</td>
<td>No</td>
<td><italic>ephB, glpK, mce1A</italic>
</td>
<td>3</td>
</tr>
<tr><td>No</td>
<td>No</td>
<td><bold>Yes</bold>
</td>
<td>No</td>
<td>No</td>
<td><bold>Yes</bold>
</td>
<td><bold><italic>emrB, mmpL6, Rv1588c</italic>
</bold>
</td>
<td><bold>3</bold>
</td>
</tr>
<tr><td>No</td>
<td>No</td>
<td><bold>Yes</bold>
</td>
<td>Yes</td>
<td>No</td>
<td><bold>Yes</bold>
</td>
<td><bold><italic>PE_PGRS10, PE_PGRS21, PE_PGRS53, pks12, pncA, rpoB, Rv2090, wag22</italic>
</bold>
</td>
<td><bold>8</bold>
</td>
</tr>
<tr><td>No</td>
<td>No</td>
<td>No</td>
<td>Yes</td>
<td>No</td>
<td>Yes</td>
<td><italic>atsD, fabG1, gyrA, helZ, murF, PE_PGRS52, Rv2000, Rv2161c, Rv3471c, ubiA</italic>
</td>
<td>10</td>
</tr>
<tr><td>No</td>
<td>No</td>
<td>No</td>
<td>Yes</td>
<td>Yes</td>
<td>Yes</td>
<td><italic>ctpB, drrA, ethA, fhaA, PPE19, PPE60, Rv0094c, Rv1945, Rv2019, Rv2141c, Rv3921c</italic>
</td>
<td>11</td>
</tr>
<tr><td>No</td>
<td>Yes</td>
<td>Yes</td>
<td>No</td>
<td>No</td>
<td>No</td>
<td><italic>accD2, eccC3, eccC4, fadE7, fas, glfT1, icd1, lepB, lprN, mbtJ, mce1R, menA, mez, PE_PGRS42, PPE30, PPE4, PPE54, rbsK, Rv1106c, Rv1147, Rv1215c, Rv1747, Rv2005c, Rv2566, Rv3085, Rv3401, Rv3677c, yfiH, zwf1</italic>
</td>
<td>29</td>
</tr>
<tr><td>No</td>
<td>No</td>
<td>No</td>
<td>No</td>
<td>Yes</td>
<td>No</td>
<td><italic>fadD5, Rv0336, Rv0348</italic>
</td>
<td>3</td>
</tr>
<tr><td>Yes</td>
<td>No</td>
<td>No</td>
<td>No</td>
<td>No</td>
<td>No</td>
<td><italic>atsA, gltA2, PE_PGRS12, Rv3863</italic>
</td>
<td>4</td>
</tr>
<tr><td>No</td>
<td>Yes</td>
<td>No</td>
<td>No</td>
<td>No</td>
<td>No</td>
<td><italic>mycP2, rpsL, Rv0012, Rv0219, Rv0893c, Rv1217c, Rv1443c, Rv1515c, Rv1516c, Rv2033c, Rv2277c, Rv2672, Rv2825c, Rv3177, Rv3618, yjcE</italic>
</td>
<td>16</td>
</tr>
<tr><td>No</td>
<td>No</td>
<td>No</td>
<td>No</td>
<td>No</td>
<td>Yes</td>
<td><italic>cysK1, espI, glmS, hsaA, lprQ, PE23, pks15, pks7, pmmB, rhlE, Rv0121c, Rv1073, Rv1639c, Rv2219, Rv3217c, Rv3254, Rv3660c, Rv3707c, thiL, vapB34</italic>
</td>
<td>20</td>
</tr>
<tr><td>No</td>
<td>No</td>
<td>Yes</td>
<td>No</td>
<td>No</td>
<td>No</td>
<td><italic>35kd_ag, accA3, adhE1, aglA, aldA, alkB, alr, amiB2, amiC, aprA, aroA, aroG, arsA, arsB2, bfrA, bfrB, blaC, cdh, cobN, cut1, cyp139, cysA2, cysD, deoC, devS, dinF, dinG, dinX, dnaB, dnaE2, dnaG, dnaK, dxs2, eccA3, eccB2, eccC2, eccD5, efp, egtE, eno, erg3, espH, esxJ, esxK, esxM, etgB, fadB3, fadB4, fadD15, fadD18, fadD22, fadD24, fadE18, fadE25, fadE34, fbiC, fecB, folC, frdD, fusA1, gatB, gcvB, glgA, gnd2, groEL2, guaB1, guaB2, hemB, hemY, hisH, hupB, hycQ, iniC, ispD, kshA, lepA, leuA, lipD, lipP, lppA, lppH, lppZ, lpqH, lprA, ltp1, mbtE, mce3B, mce3C, mcr7, menC, metH, mmpL1, mmpL11, mmpL5, mmsA, mmuM, moeB1, murE, mutT3, mymT, narG, nlhH, nrp, oppB, oxyR, papA2, papA5, PE_PGRS1, PE_PGRS13, PE_PGRS16, PE_PGRS17, PE_PGRS19, PE_PGRS2, PE_PGRS22, PE_PGRS23, PE_PGRS24, PE_PGRS25, PE_PGRS26, PE_PGRS27, PE_PGRS30, PE_PGRS31, PE_PGRS32, PE_PGRS33, PE_PGRS38, PE_PGRS4, PE_PGRS43, PE_PGRS44, PE_PGRS45, PE_PGRS46, PE_PGRS47, PE_PGRS48, PE_PGRS50, PE_PGRS54, PE_PGRS55, PE_PGRS56, PE_PGRS57, PE_PGRS6, PE_PGRS62, PE_PGRS7, PE_PGRS9, PE1, pepC, pepD, pepR, php, pknB, pknJ, pknK, pknL, pks13, pks9, pncB2, pntAb, ppa, PPE11, PPE12, PPE15, PPE21, PPE24, PPE34, PPE42, PPE46, PPE55, PPE56, PPE57, PPE59, PPE6, PPE8, ppm1, ppnK, ppsA, ppsC, proV, prpD, psd, purH, purL, purT, recC, relE, rfe, rimJ, rpfA, rpsD, Rv0057, Rv0060, Rv0074, Rv0095c, Rv0111, Rv0133, Rv0149, Rv0178, Rv0180c, Rv0188, Rv0192, Rv0193c, Rv0203, Rv0218, Rv0238, Rv0257, Rv0278c, Rv0302, Rv0323c, Rv0339c, Rv0367c, Rv0378, Rv0386, Rv0412c, Rv0457c, Rv0474, Rv0477, Rv0492c, Rv0537c, Rv0538, Rv0575c, Rv0584, Rv0647c, Rv0650, Rv0658c, Rv0686, Rv0696, Rv0740, Rv0790c, Rv0842, Rv0845, Rv0874c, Rv0876c, Rv0921, Rv0939, Rv0996, Rv1004c, Rv1024, Rv1043c, Rv1045, Rv1047, Rv1048c, Rv1050, Rv1069c, Rv1085c, Rv1178, Rv1225c, Rv1273c, Rv1313c, Rv1318c, Rv1372, Rv1395, Rv1425, Rv1455, Rv1488, Rv1498A, Rv1499, Rv1526c, Rv1610, Rv1667c, Rv1683, Rv1702c, Rv1707, Rv1714, Rv1725c, Rv1726, Rv1728c, Rv1729c, Rv1739c, Rv1749c, Rv1815, Rv1825, Rv1828, Rv1831, Rv1835c, Rv1879, Rv1895, Rv1941, Rv1999c, Rv2017, Rv2023A, Rv2024c, Rv2025c, Rv2038c, Rv2074, Rv2134c, Rv2143, Rv2164c, Rv2177c, Rv2190c, Rv2219A, Rv2228c, Rv2254c, Rv2256a, Rv2286c, Rv2314c, Rv2319c, Rv2327, Rv2424c, Rv2456c, Rv2532c, Rv2556c, Rv2567, Rv2621c, Rv2630, Rv2638, Rv2656c, Rv2668, Rv2669, Rv2729c, Rv2733c, Rv2792c, Rv2798c, Rv2851c, Rv2864c, Rv2893, Rv2897c, Rv2951c, Rv2974c, Rv2994, Rv3027c, Rv3060c, Rv3064c, Rv3074, Rv3104c, Rv3169, Rv3193c, Rv3201c, Rv3208, Rv3216, Rv3218, Rv3230c, Rv3263, Rv3365c, Rv3400, Rv3439c, Rv3549c, Rv3668c, Rv3689, Rv3717, Rv3728, Rv3729, Rv3807c, Rv3856c, Rv3896c, Rv3898c, sirA, smc, tal, TB9.4, tgs2, thiD, tig, treX, trxB1, uspC, uvrA, vapB36, vapC29, whiA, yrbE3A, yrbE3B</italic>
</td>
<td>356</td>
</tr>
</tbody>
</table>
<table-wrap-foot><p>Phenotype (1): clustered strains (cluster[+]) (pairwise SNV difference between two isolates is no more than five SNVs) in Hanoi’s study cohort; (2): strains harboring <italic>katG</italic>
-S315 T mutation (S315 T[+]) in Hanoi’s cohort; (3): cluster[+]/S315T[+] in Hanoi’s cohort; (4): cluster[+] in KwaZulu-Natal’s cohort; (5) S315T[+] in KwaZulu-Natal’s cohort; (6): cluster[+]/S315T[+] in KwaZulu-Natal’s cohort; SNV: single nucleotide variant.</p>
<p>The genes in bold are associated with cluster[+]/S315T[+], and shared between Hanoi and KwaZulu-Natal study population.</p>
</table-wrap-foot>
</table-wrap>
<fig id="Fig2"><label>Figure 2</label>
<caption><p>Venn diagram showing the number of genes significantly associated with the given phenotypes and shared by the two study cohorts. (<bold>a</bold>
) results from k-mer-based GWAS; (<bold>b</bold>
) results from phyOverlap; cluster[+]: clustered (pairwise SNV difference between two isolates is no more than five SNVs); S315T[+]: harboring <italic>katG</italic>
-S315T mutation; Hanoi: Hanoi’s study cohort; KwaZulu-Natal: KwaZulu-Natal’s study cohort.</p>
</caption>
<graphic xlink:href="41598_2019_51812_Fig2_HTML" id="d29e2040"></graphic>
</fig>
</p>
<p id="Par13">The remaining five genes, namely, <italic>PPE18</italic>
, <italic>gid, emrB, Rv1588c</italic>
, and <italic>pncA</italic>
, were confirmed by BLAST search and were aligned to the H37Rv reference sequence (Supplementary Fig. <xref rid="MOESM1" ref-type="media">S2</xref>
) and nominated as real candidates. In the Hanoi study population, significant k-mers included 36 k-mers from the <italic>PPE18</italic>
gene (all had the same P value = 1.840E-09), 33 k-mers from <italic>gid</italic>
(the best P value = 5.213E-08), 31 from <italic>emrB</italic>
(8.579E-08), 4 from <italic>Rv1588c</italic>
(8.579E-08), and 31 from the <italic>pncA</italic>
gene (4.437E-07) (Supplementary Table <xref rid="MOESM1" ref-type="media">S1</xref>
), when 6.208E-06 was applied as the threshold of statistical significance after Bonferroni correction. For the KwaZulu-Natal study population, within the same set of genes, the best P values obtained from 5, 101, 31, 12, and 29 k-mers were 2.998E-08, 2.998E-08, 7.802E-07, 2.453E-06, and 4.802E-11, respectively, with the threshold of significance after Bonferroni correction as 4.591E-06 (Supplementary Table <xref rid="MOESM1" ref-type="media">S1</xref>
). Positive associations with the phenotype were mainly observed in L2 of Hanoi’s samples and L4 of KwaZulu-Natal’s samples (Fig. <xref rid="Fig3" ref-type="fig">3</xref>
).<fig id="Fig3"><label>Figure 3</label>
<caption><p>Distribution of k-mers derived from variant (or wild) types of five genes (<italic>PPE18/19, gid, emrB, Rv1588c</italic>
, and <italic>pncA</italic>
) that showed positive or negative associations with clustered strains carrying <italic>katG</italic>
-S315T in the phylogenetic trees of Hanoi (<bold>a</bold>
) and KwaZulu-Natal (<bold>b</bold>
) study population.</p>
</caption>
<graphic xlink:href="41598_2019_51812_Fig3_HTML" id="d29e2104"></graphic>
</fig>
</p>
<p id="Par14">Next, we attempted to confirm significant k-mers by using a DBGWAS approach to identify all the relevant sequences as “unitigs” (Table <xref rid="Tab3" ref-type="table">3</xref>
). Consequently, one significant unitig indicating merged k-mers was identified in <italic>PPE18</italic>
, three were identified <italic>in gid</italic>
, two in <italic>emrB</italic>
, and two in <italic>Rv1588c</italic>
(the best q value in each gene was 5.562E-07, 3.676E-05, 1.704E-05, and 1.704E-05, respectively) for the Hanoi samples (Supplementary Table <xref rid="MOESM1" ref-type="media">S2</xref>
, Fig. <xref rid="MOESM1" ref-type="media">S3</xref>
). For the KwaZulu-Natal samples, one significant unitig was found in <italic>PPE18</italic>
, seven in <italic>gid</italic>
, two in <italic>emrB</italic>
, and one in <italic>pncA</italic>
(the best q value in each gene was 4.856E-10, 7.943E-09, 7.943E-09, and 7.639E-12, respectively) (Supplementary Table <xref rid="MOESM1" ref-type="media">S2</xref>
).<table-wrap id="Tab3"><label>Table 3</label>
<caption><p>Genes with k-mers significantly associated with clustered strains harboring <italic>katG</italic>
-S315T mutation (cluster[+]/S315T[+]) in both Hanoi and KwaZulu-Natal study cohorts, further investigated by using the DBGWAS platform and variant/deletion-based search.</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th>Gene name</th>
<th>Sample panel</th>
<th>K-mer-based GWAS</th>
<th>Strongest P value</th>
<th>DBGWAS*</th>
<th>Strongest P value</th>
<th>Variant-based GWAS/deletion-based GWAS/BLAST search</th>
<th>Strongest P value</th>
<th>Conclusion on variant(s) corresponding to significant k-mer(s)</th>
</tr>
</thead>
<tbody><tr><td><italic>PPE18/19</italic>
</td>
<td>Hanoi</td>
<td>36 k-mers, SNV at 1339644 and 1339649</td>
<td>1.840E-09</td>
<td>1 unitig, SNVs at 1339644 and 1339649</td>
<td>5.56188E-07</td>
<td>SNVs at 1339644 (E99A) and 1339649 (A101T) found by variant-based GWAS</td>
<td>1.941E-09</td>
<td>Significant SNVs E99A and A101T</td>
</tr>
<tr><td><italic>PPE18/19</italic>
</td>
<td>KwaZulu-Natal</td>
<td>5 k-mers, SNV at 1340231</td>
<td>2.998E-08</td>
<td>1 unitig**, SNV at 1532777</td>
<td>4.8561E-10</td>
<td>SNV at 1532777 (Q286R) detected by BLAST search. By variant-based GWAS, SNV was not called properly</td>
<td>NA</td>
<td>SNV Q286R</td>
</tr>
<tr><td><italic>gid</italic>
</td>
<td>Hanoi</td>
<td>33 k-mers, SNV at 4407686</td>
<td>5.213E-08</td>
<td>3 unitgs, SNV at 4407686</td>
<td>3.676E-05</td>
<td>SNV at 4407686 (E173*) fount by variant-based GWAS</td>
<td>2.991E-06</td>
<td>Significant SNV E173*</td>
</tr>
<tr><td><italic>gid</italic>
</td>
<td>KwaZulu-Natal</td>
<td>101 k-mers</td>
<td>7.802E-07</td>
<td>7 unitigs, no SNV</td>
<td>7.943E-09</td>
<td>Large deletions (from 120 bp to 675 bp)</td>
<td>2.178E-04</td>
<td>Large deletions</td>
</tr>
<tr><td><italic>emrB</italic>
</td>
<td>Hanoi</td>
<td>31 k-mers, SNV at 876918</td>
<td>8.579E-08</td>
<td>2 unitigs, SNV at 876918</td>
<td>1.704E-05</td>
<td>SNV at 876918 (F508S) found by variant-based GWAS</td>
<td>1.507E-08</td>
<td>Significant SNV F508S</td>
</tr>
<tr><td><italic>emrB</italic>
</td>
<td>KwaZulu-Natal</td>
<td>31 k-mers, SNV at 877058</td>
<td>7.802E-07</td>
<td>2 unitigs**, SNV at 877058</td>
<td>7.943E-09</td>
<td>SNV at 877058 (I461I) found by variant-based GWAS</td>
<td>1.077E-06</td>
<td>Significant SNV I461I</td>
</tr>
<tr><td><italic>Rv1588c</italic>
</td>
<td>Hanoi</td>
<td>4 k-mers without SNV (neg.***)</td>
<td>8.579E-08</td>
<td>2 unitigs**, no SNV</td>
<td>1.704E-05</td>
<td>SNV at 1789735 (P34P) found by variant-based GWAS and BLAST search</td>
<td>1.511E-08</td>
<td>Significant SNV P34P</td>
</tr>
<tr><td><italic>Rv1588c</italic>
</td>
<td>KwaZulu-Natal</td>
<td>12 k-mers without SNV (neg.***)</td>
<td>2.453E-06</td>
<td>0 unitig</td>
<td>NA</td>
<td>No SNV inside the corresponding k-mers</td>
<td>NA</td>
<td>None</td>
</tr>
<tr><td><italic>pncA</italic>
</td>
<td>Hanoi</td>
<td>31 k-mers without SNV (neg.***)</td>
<td>4.437E-07</td>
<td>0 unitig (<italic>pncA</italic>
not found)</td>
<td>NA</td>
<td>SNV at 2288820 (Q141P) found by variant-based GWAS</td>
<td>4.919E-07</td>
<td>Significant SNV Q141P</td>
</tr>
<tr><td><italic>pncA</italic>
</td>
<td>KwaZulu-Natal</td>
<td>29 k-mers, SNV at 2288785 and 2288788</td>
<td>4.802E-11</td>
<td>1 unitig, G insertion at 2288785</td>
<td>7.639E-12</td>
<td>SNV at 2288785 (T153fs) found by variant-based GWAS</td>
<td>6.569E-11</td>
<td>Significant SNV T153fs</td>
</tr>
</tbody>
</table>
<table-wrap-foot><p>cluster[+]: clustered (pairwise SNV difference between two isolates is no more than five SNVs); S315T[+]: harboring <italic>katG</italic>
-S315T mutation; Hanoi: Hanoi’s study cohort; KwaZulu-Natal: KwaZulu-Natal’s study cohort; GWAS: genome wide association study; SNV: single nucleotide variant.</p>
<p>*Only significant unitigs and their SNVs are shown.</p>
<p>**Found by sequence search.</p>
<p>***Showing negative association with Cluster[+]/S315T[+] by logistic regression analyses</p>
</table-wrap-foot>
</table-wrap>
</p>
<p id="Par15">All k-mers annotated with <italic>PPE18</italic>
in both study sets were mapped to either <italic>PPE18, 19</italic>
, or 60. Nucleotide sequences within these three <italic>PPE</italic>
genes were hardly distinguishable from each other (Supplementary Fig. <xref rid="MOESM1" ref-type="media">S4</xref>
); BWA-MEM mapping or simple BLAST search to the reference could not specify the exact <italic>PPE</italic>
gene that each 31-bp k-mer belonged to. Nevertheless, the analysis of <italic>de novo</italic>
assembled contigs demonstrated that the Hanoi variant k-mers were derived from <italic>PPE18</italic>
and KwaZulu-Natal’s k-mers initially annotated with <italic>PPE18</italic>
were derived from <italic>PPE19</italic>
. Such a high degree of sequence similarity was also observed in significant unitigs identified in both study sets with the DBGWAS approach. Therefore, these k-mers and unitigs were designated as <italic>PPE18/19</italic>
in our study (Table <xref rid="Tab3" ref-type="table">3</xref>
).</p>
</sec>
<sec id="Sec7"><title>Individual variant analysis using a genome-wide approach</title>
<p id="Par16">When mapping the aforementioned significant k-mers of Hanoi samples on the H37Rv genome, we identified variants corresponding to the significant k-mers. Next, we attempted to confirm the variants by SNV- and small indel-based GWAS or BLAST search, and finally identified mutations encoding E99A and A101T, in <italic>PPE18</italic>
, E173* (Glu173Stop) in the <italic>gid</italic>
, F508S in <italic>emrB</italic>
, P34P in <italic>Rv1588c</italic>
, and Q141P in <italic>pncA</italic>
. They were significantly associated with cluster[+]/S315T[+] (P values were 1.941E-09 for both <italic>PPE18</italic>
SNVs, and the P values were 2.991E-06, 1.507E-08, 1.511E-08, and 4.919E-07 for <italic>gid, emrB</italic>
, <italic>Rv1588c</italic>
, and <italic>pncA</italic>
SNV, respectively) (Table <xref rid="Tab3" ref-type="table">3</xref>
). The PROVEAN web server<sup><xref ref-type="bibr" rid="CR24">24</xref>
</sup>
and SIFT_4G tool<sup><xref ref-type="bibr" rid="CR25">25</xref>
</sup>
predicted <italic>gid</italic>
-E173* and <italic>emrB</italic>
-F508S as deleterious or that affecting the protein function.</p>
<p id="Par17">For the KwaZulu-Natal samples, mapping results showed SNVs in <italic>emrB</italic>
(I461I) and <italic>pncA</italic>
(T153fs) (variant-based GWAS P values = 1.077E-06 and 6.596E-11, respectively). Initially, the SNV in the <italic>PPE18/19</italic>
was not clearly mapped to the reference for the aforementioned reason. After the in-depth search for <italic>de novo</italic>
assembled contigs, Q286R in <italic>PPE19</italic>
was identified. Large deletions were found in <italic>gid</italic>
in 101 of the 337 samples (from 120 to 675 bp). <italic>Rv1588c</italic>
did not show any SNV in the corresponding area of the k-mers (Table <xref rid="Tab3" ref-type="table">3</xref>
).</p>
<p id="Par18">SNVs and small indels from VCF files, significantly associated with the cluster[+]/S315T[+] of the Hanoi cohort obtained by GWAS, included 329 SNVs and indels other than the six shown above (Supplementary Table <xref rid="MOESM1" ref-type="media">S3</xref>
) (the best P = 1.754E-13, Supplementary Table <xref rid="MOESM1" ref-type="media">S4</xref>
). None of them were lineage-specific SNVs that have been reported elsewhere<sup><xref ref-type="bibr" rid="CR22">22</xref>
</sup>
. After further excluding variants in ambiguous <italic>PE</italic>
and <italic>PPE</italic>
genes, synonymous SNVs, a mutation causing S315T, and those with frequency no more than two in case or control samples, 144 SNVs and small indels were extracted in the samples from Hanoi. Of these, 77 SNVs correlated well with the principal component (PC)-7 (Supplementary Table <xref rid="MOESM1" ref-type="media">S4</xref>
), which corresponded to a phylogenetic branch of the ancient Beijing strains (Supplementary Fig. <xref rid="MOESM1" ref-type="media">S5</xref>
), with a high percentage (92.9%) of INH resistance.</p>
<p id="Par19">We further analyzed the structural variants by detecting zero or low-depth areas when H37Rv and four complete genomes in Hanoi (AL123456, and AP018033 to AP018036) were used as references to be mapped. We found two groups of gene deletions associated with this phenotype. The first group included 1-bp deletion in <italic>Rv1043c</italic>
, 52-bp deletion in <italic>Rv2286c</italic>
, 3-bp insertion in <italic>Rv0790c</italic>
, 1-bp deletion in <italic>Rv3230c</italic>
, and a big deletion in <italic>Rv2025c</italic>
, and these correlated well with PC-9 (Table <xref rid="Tab4" ref-type="table">4</xref>
), which corresponded to the same branch of ancient Beijing strains shown in SNV analysis. The second group included a 238-bp deletion in <italic>accD2</italic>
, 2-bp deletion in <italic>eis</italic>
, a 459-bp insertion in <italic>Rv3077</italic>
, a 12-bp deletion in <italic>Rv2690c</italic>
, and correlated with PC-12 (Table <xref rid="Tab4" ref-type="table">4</xref>
) corresponding to another branch of modern Beijing strains.<table-wrap id="Tab4"><label>Table 4</label>
<caption><p>Genes with deletions/insertions significantly associated with clustered strains harboring <italic>katG</italic>
-S315T mutation among Hanoi samples obtained by GWAS, and corresponding k-mers.</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th rowspan="2">Reference gene*</th>
<th rowspan="2">Locus in H37Rv (AL123456)</th>
<th>Deletion/insertion</th>
<th rowspan="2">P value</th>
<th rowspan="2">Principal component</th>
<th>K-mers</th>
<th rowspan="2">Strongest P value</th>
</tr>
<tr><th>Details</th>
<th>No of significant k-mers/No of k-mers</th>
</tr>
</thead>
<tbody><tr><td>AL123456.3H37Rv|00828|<italic>Rv0790c</italic>
</td>
<td>Rv0790c</td>
<td>3-bp insertion at 884090 to 884092</td>
<td>1.507E-08</td>
<td>9</td>
<td>42/46</td>
<td>1.43E-08</td>
</tr>
<tr><td>AL123456.3H37Rv|01092|<italic>Rv1043c</italic>
</td>
<td>Rv1043c</td>
<td>1-bp deletion at 1166499</td>
<td>1.941E-09</td>
<td>9</td>
<td>45/61</td>
<td>1.43E-08</td>
</tr>
<tr><td>AL123456.3H37Rv|02117|<italic>Rv2025c</italic>
</td>
<td>Rv2025c</td>
<td>Deletion of almost all Rv2025c gene, also Rv2024c (789-bp)</td>
<td>9.053E-08</td>
<td>9</td>
<td>819/833</td>
<td>1.43E-08</td>
</tr>
<tr><td>AL123456.3H37Rv|02393|<italic>Rv2286c</italic>
</td>
<td>Rv2286c</td>
<td>52-bp deletion from 2559452 to 2559503</td>
<td>1.507E-08</td>
<td>9</td>
<td>69/123</td>
<td>1.43E-08</td>
</tr>
<tr><td>AL123456.3H37Rv|03394|<italic>Rv3230c</italic>
</td>
<td>Rv3230c</td>
<td>1-bp deletion at 3607421</td>
<td>9.053E-08</td>
<td>9</td>
<td>59/63</td>
<td>1.84E-09</td>
</tr>
<tr><td>AL123456.3H37Rv|02535|<italic>eis</italic>
</td>
<td>Rv2416c</td>
<td>2-bp deletion at 2714527–2714528</td>
<td>3.219E-06</td>
<td>12</td>
<td>0/60</td>
<td>1.32E-05</td>
</tr>
<tr><td>AP018033.1HN024|01030|<italic>accD2</italic>
</td>
<td>Best hit with Rv0974c</td>
<td>238-bp deletion from 1086195 to 1086432</td>
<td>2.991E-06</td>
<td>12</td>
<td>—</td>
<td>—</td>
</tr>
<tr><td>AP018033.1HN024|02797|<italic>HN024_02796</italic>
</td>
<td>Best hit with Rv2690c</td>
<td>12-bp deletion</td>
<td>1.720E-05</td>
<td>12</td>
<td>—</td>
<td>—</td>
</tr>
<tr><td>AP018036.1HN506|03217|<italic>HN506_03216</italic>
</td>
<td>Best hit with Rv3077</td>
<td>459-bp insertion</td>
<td>1.720E-05</td>
<td>12</td>
<td>—</td>
<td>—</td>
</tr>
</tbody>
</table>
<table-wrap-foot><p>PC-9 correlated gene variants were observed in 14 samples (HN-042, HN-065, HN-120, HN-127, HN-152, HN-154, HN-170, HN-179, HN-186, HN-222, HN-253, HN-391, HN-441, HN-456) that belong to ancient Beijing sublineages.</p>
<p>PC-12 correlated gene variants were observed in 10 samples (HN-019, HN-075, HN-109, HN-137, HN-169, HN-324, HN-330, HN-435, HN-462, HN-506) that belong to modern Beijing sublineages.</p>
<p>*Genome sequences of clinical isolates in our Hanoi cohort, AP018033 to AP018036 as well as AL123456 (H37Rv genome) were used for reference.</p>
</table-wrap-foot>
</table-wrap>
</p>
</sec>
<sec id="Sec8"><title>Convergence-based phyOverlap analysis</title>
<p id="Par20">Seeking variants caused by convergent evolution is another alternative for detecting mutations supporting drug resistance, such as compensatory mutations<sup><xref ref-type="bibr" rid="CR26">26</xref>
</sup>
. We also tried to identify the phenotype-associated variants caused by convergence evolution. As expected, <italic>rpoB</italic>
and <italic>pncA</italic>
were also associated with <italic>katG</italic>
-S315T mutation alone. However, no genes and inter-genic regions were significantly associated with the presence of cluster[+]/S315T[+] in both KwaZulu-Natal and our Hanoi panels (Table <xref rid="Tab5" ref-type="table">5</xref>
, Supplementary Table <xref rid="MOESM1" ref-type="media">S5</xref>
, Fig. <xref rid="Fig2" ref-type="fig">2b</xref>
). One <italic>PPE</italic>
gene, two <italic>PE</italic>
genes, and one intergenic region, <italic>PPE47, PE_PGRS55, PE_PGRS20</italic>
and <italic>PE_PGRS3-PE_PGRS4</italic>
were extracted in the Hanoi samples only (Table <xref rid="Tab5" ref-type="table">5</xref>
).<table-wrap id="Tab5"><label>Table 5</label>
<caption><p>Genes detected by the phyOverlap method and their significant associations with the Hanoi strains harboring the <italic>katG</italic>
-S315T mutation, with and without clustering, and loci shared with the KwaZulu-Natal study population.</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th>Locus</th>
<th>Gene name</th>
<th>No of SNV</th>
<th>Parsimony score</th>
<th>P value</th>
<th>FDR (q value)</th>
<th>Remark</th>
</tr>
</thead>
<tbody><tr><td colspan="7"><bold><italic>Associated with cluster</italic>
</bold>
[+]<bold><italic>/S315T</italic>
</bold>
[+]</td>
</tr>
<tr><td>Rv3021c</td>
<td><italic>PPE47</italic>
</td>
<td>12</td>
<td>73</td>
<td>0.00000</td>
<td>0.00000</td>
<td>PE/PPE</td>
</tr>
<tr><td>Rv3511</td>
<td><italic>PE_PGRS55</italic>
</td>
<td>13</td>
<td>18</td>
<td>0.00000</td>
<td>0.00000</td>
<td>PE-PGRS family protein</td>
</tr>
<tr><td>Rv0278c-Rv0279c</td>
<td><italic>PE_PGRS3-PE_PGRS4</italic>
</td>
<td>2</td>
<td>37</td>
<td>0.00002</td>
<td>0.02434</td>
<td>PE-PGRS family protein</td>
</tr>
<tr><td>Rv1068c</td>
<td><italic>PE_PGRS20</italic>
</td>
<td>13</td>
<td>27</td>
<td>0.00002</td>
<td>0.02434</td>
<td>PE-PGRS family protein</td>
</tr>
<tr><td colspan="7"><bold><italic>Associated with S315T</italic>
</bold>
[+]</td>
</tr>
<tr><td><bold>Rv0667</bold>
</td>
<td><bold><italic>rpoB</italic>
</bold>
</td>
<td><bold>14</bold>
</td>
<td><bold>27</bold>
</td>
<td><bold>0.00000</bold>
</td>
<td><bold>0.00000</bold>
</td>
<td><bold>Rifampicin-resistance related gene</bold>
</td>
</tr>
<tr><td>Rv0682</td>
<td><italic>rpsL</italic>
</td>
<td>2</td>
<td>21</td>
<td>0.00000</td>
<td>0.00000</td>
<td>Streptomycin-resistance related gene</td>
</tr>
<tr><td>Rv0758</td>
<td><italic>phoR</italic>
</td>
<td>8</td>
<td>8</td>
<td>0.00002</td>
<td>0.01217</td>
<td>Sensor part of a two component regulatory system</td>
</tr>
<tr><td><bold>Rv2043c</bold>
</td>
<td><bold><italic>pncA</italic>
</bold>
</td>
<td><bold>10</bold>
</td>
<td><bold>11</bold>
</td>
<td><bold>0.00000</bold>
</td>
<td><bold>0.00000</bold>
</td>
<td><bold>Pyrazinamide-resistance related gene</bold>
</td>
</tr>
<tr><td>Rv3021c</td>
<td><italic>PPE47</italic>
</td>
<td>12</td>
<td>73</td>
<td>0.00000</td>
<td>0.00000</td>
<td>PE/PPE</td>
</tr>
<tr><td>Rv3169</td>
<td><italic>Rv3169</italic>
</td>
<td>5</td>
<td>5</td>
<td>0.00002</td>
<td>0.01217</td>
<td>Unknown</td>
</tr>
<tr><td>Rv3428c-Rv3429</td>
<td><italic>Rv3428c-PPE59</italic>
</td>
<td>28</td>
<td>247</td>
<td>0.00000</td>
<td>0.00000</td>
<td></td>
</tr>
<tr><td>Rv3647c-Rv3648c</td>
<td><italic>Rv3647c-cspA</italic>
</td>
<td>4</td>
<td>4</td>
<td>0.00006</td>
<td>0.03245</td>
<td></td>
</tr>
<tr><td>Rv3680-Rv3681c</td>
<td><italic>Rv3680-whiB4</italic>
</td>
<td>1</td>
<td>78</td>
<td>0.00000</td>
<td>0.00000</td>
<td></td>
</tr>
</tbody>
</table>
<table-wrap-foot><p>The genes in bold are shared between Hanoi and KwaZulu-Natal study population.</p>
</table-wrap-foot>
</table-wrap>
</p>
<p id="Par21"><italic>PPE46</italic>
and <italic>PPE47</italic>
share large portions of identical nucleotide sequences. By BLAST search for variants in the <italic>de novo</italic>
assembled contigs in addition to genome-wide screening of zero or low-depth areas (<15% of average depth) in the reference genome, we found that 37 (4 in L1, 19 in L2, and 14 in L4) strains had large deletions in <italic>PPE47</italic>
(Supplementary Fig. <xref rid="MOESM1" ref-type="media">S6</xref>
) and all had fusion with <italic>PPE46</italic>
, resulting in <italic>PPE46</italic>
-like chimeric genes in 35 isolates and <italic>PPE47</italic>
-like chimeric genes in 2 isolates (Supplementary Fig. <xref rid="MOESM1" ref-type="media">S7</xref>
). SNVs identified from 3379708 to 3379763 of AL123456.3 (H37Rv) in the <italic>PPE47</italic>
region were the main reason for the significant association in the phylogenetic convergence test. In other candidates from phyOverlap, specific variants were not validated, presumably owing to the difficulties in mapping short reads followed by ambiguous base calling within <italic>PE_PGRS</italic>
genes.</p>
</sec>
<sec id="Sec9"><title>Analyses using logistic regression models adjusted for host confounders in Hanoi samples</title>
<p id="Par22">By multivariate analyses using conventional logistic regression models after adjustment for patients’ gender, age, living area, as well as Mtb lineages, all k-mers from <italic>PPE18/19</italic>
, <italic>emrB</italic>
, and a part of k-mers from <italic>gid</italic>
showed positive associations with clustered strains carrying <italic>katG</italic>
-S315T mutations in the Hanoi study population (adjusted odds ratio [aOR] with 95% confidence interval [CI] = 13.20 [3.49–49.96], 11.98 [3.24–44.29], and 12.42 [2.81–54.90], respectively), whereas <italic>Rv1588c</italic>
and <italic>pncA</italic>
k-mers showed negative associations (aOR with 95% CI = 0.08 [0.02–0.31], and 0.01 [0.00–0.25], respectively). All variants corresponding to significant k-mers showed positive associations with cluster[+]/S315T[+]. The <italic>PPE46/47</italic>
-like chimeric gene also showed positive association with cluster[+]/S315T[+] (aOR with 95% CI = 6.81 [2.13–21.72]) (Supplementary Table <xref rid="MOESM1" ref-type="media">S6</xref>
).</p>
</sec>
</sec>
<sec id="Sec10" sec-type="discussion"><title>Discussion</title>
<p id="Par23">We identified a variety of drug resistance-conferring mutations prevailing in Hanoi, northern Vietnam, which appeared most frequently in the East Asian Mtb lineage L2 particularly in ancient Beijing sublineage, and less in L1, and then least in L4. Large deletions that were not detected by conventional variant calling from mapped short reads were also found in <italic>pncA</italic>
. Using the bacterial GWAS approach, we extracted candidate genes that were significantly associated with the clustered strains harboring the <italic>katG</italic>
-S315T mutation and that were common to the two independent data sets—our cohort panel and a previous South African study, the KwaZulu-Natal cohort by Cohen <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR19">19</xref>
</sup>
—suggesting that expansions of INH-resistance TB can be facilitated by pathogen factors, at least in part.</p>
<p id="Par24">A major drug-resistance conferring mutation, <italic>katG</italic>
-S315T, accounted for 85.3% of the INH-resistance strains in our study population, which was similar to that reported in other studies in Vietnam, e.g., 73.2% was reported by Minh <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR13">13</xref>
</sup>
, 73.6% was reported by Huyen <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR27">27</xref>
</sup>
, 81.3% by Nguyen <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR15">15</xref>
</sup>
and 78.4% was reported in Southeast Asian countries<sup><xref ref-type="bibr" rid="CR28">28</xref>
</sup>
. This <italic>katG</italic>
-S315T is known as a mutation with low-fitness cost, spreads to Beijing strains and others, and is more likely to be clustered<sup><xref ref-type="bibr" rid="CR23">23</xref>
,<xref ref-type="bibr" rid="CR29">29</xref>
,<xref ref-type="bibr" rid="CR30">30</xref>
</sup>
. The acquisition of <italic>katG</italic>
-S315T mutation precedes other resistance mutations, including RMP<sup><xref ref-type="bibr" rid="CR4">4</xref>
,<xref ref-type="bibr" rid="CR5">5</xref>
,<xref ref-type="bibr" rid="CR19">19</xref>
</sup>
. Indeed, this mutation occurs more frequently in MDR-TB strains than other mutations<sup><xref ref-type="bibr" rid="CR30">30</xref>
,<xref ref-type="bibr" rid="CR31">31</xref>
</sup>
, and it has an important implication in the transmission and control of MDR-TB.</p>
<p id="Par25">In this setting, we performed bacterial GWAS with phenotypic combination of clusters defined by <6 SNVs<sup><xref ref-type="bibr" rid="CR32">32</xref>
</sup>
and <italic>katG</italic>
-S315T mutation as a surrogate, and searched for pathogenic variants contributing to the spread of INH-resistance. We identified five genes, namely, <italic>PPE18/19</italic>
, <italic>gid, emrB, Rv1588c</italic>
, and <italic>pncA</italic>
, which were shared by two different sample panels in Asia and Africa with different major Mtb lineages, L2 and L4, respectively.</p>
<p id="Par26">Of the five genes extracted from the k-mer GWAS, except for drug-resistance-conferring genes, the best P values were obtained from <italic>PPE18/19</italic>
both in the Hanoi population and in the KwaZulu-Natal cohort. The association was mainly observed in L2 of Hanoi’s samples and L4 of KwaZulu-Natal’s samples. <italic>PPE18/19</italic>
genes are members of a multigene family and share high sequence similarity with another PPE gene, <italic>PPE60</italic>
<sup><xref ref-type="bibr" rid="CR33">33</xref>
</sup>
. Among these <italic>PPE</italic>
genes, homologous recombination events frequently occur and contribute to the sequence diversity<sup><xref ref-type="bibr" rid="CR34">34</xref>
,<xref ref-type="bibr" rid="CR35">35</xref>
</sup>
. Among nonsynonymous SNVs and small indels detected in the <italic>PPE18</italic>
gene of our study population leading to amino acid change in Mtb39A, two SNVs were significantly associated with cluster[+]/S315T[+], suggesting the spread of INH-resistant TB. Although their role in virulence is not fully understood, the PPE18 protein, also known as Mtb39A, was shown to downregulate the proinflammatory response and Th1-type immunity, interacting with host TLR2<sup><xref ref-type="bibr" rid="CR36">36</xref>
</sup>
, and facilitated survival and multiplication of Mtb bacilli in a mouse model<sup><xref ref-type="bibr" rid="CR37">37</xref>
</sup>
.</p>
<p id="Par27">The T-cell epitopes in Mtb39A<sup><xref ref-type="bibr" rid="CR38">38</xref>
–<xref ref-type="bibr" rid="CR40">40</xref>
</sup>
lead to strong T-cell proliferation and IFN-gamma production<sup><xref ref-type="bibr" rid="CR38">38</xref>
</sup>
. Therefore it is used as a subunit for the human TB vaccine candidate Mtb72F and its successor M72. It has been proved to be immunogenic and can stimulate both cellular and humoral immune response<sup><xref ref-type="bibr" rid="CR41">41</xref>
,<xref ref-type="bibr" rid="CR42">42</xref>
</sup>
. By comparing with the lists of the T-cell epitopes published elsewhere<sup><xref ref-type="bibr" rid="CR38">38</xref>
,<xref ref-type="bibr" rid="CR40">40</xref>
,<xref ref-type="bibr" rid="CR43">43</xref>
,<xref ref-type="bibr" rid="CR44">44</xref>
</sup>
, nonsynonymous mutations E99A and A101T in <italic>PPE18</italic>
identified in the Hanoi population and Q286R in <italic>PPE19</italic>
in the KwaZulu-Natal’s panel were both found to be included in the sequences acting as T-cell epitopes for cellular immunity. Despite the presence of <italic>PPE18</italic>
variants in the Hanoi strains, the original T-cell epitope sequences were mostly conserved in either <italic>PPE19</italic>
or <italic>PPE60</italic>
, when the BLAST search was applied. Variants of these <italic>PPE</italic>
genes associated with spread of INH-resistance strains may thus increase the antigenic diversity of the bacilli, which may help evade or exploit human immunity by unidentified mechanisms through the process of human-Mtb coevolution<sup><xref ref-type="bibr" rid="CR33">33</xref>
,<xref ref-type="bibr" rid="CR35">35</xref>
,<xref ref-type="bibr" rid="CR45">45</xref>
</sup>
.</p>
<p id="Par28">Mutations in the <italic>gid</italic>
gene have been reported to be associated with low-level resistance to SM<sup><xref ref-type="bibr" rid="CR46">46</xref>
</sup>
, which has been used as the first-line drug since a long time until recently. Because <italic>gid</italic>
encodes methyltransferase that is responsible for the methylation of 16 S rRNA involved in the translational fidelity<sup><xref ref-type="bibr" rid="CR47">47</xref>
</sup>
, it is thus conceivable that the <italic>gid</italic>
mutation may modulate the fitness of INH-resistance conferring mutations through the change in mRNA translational fidelity. Our study revealed that <italic>gid</italic>
k-mers with E173* mutant were significantly associated with cluster[+]/S315T[+] and S315T[+], but not with cluster[+] alone, indicating that a concurrence of <italic>gid</italic>
-E173* and <italic>katG</italic>
-S315T may facilitate transmission even after controlling for population structure. The concurrence of <italic>gid</italic>
130-bp deletion and <italic>katG</italic>
-S315T is the first step toward XDR-level drug resistance in Africa<sup><xref ref-type="bibr" rid="CR19">19</xref>
</sup>
. In our study population, all strains with <italic>gid</italic>
-E173* mutation had at least one genetic mutation conferring resistance to first-line drugs. This non-sense mutation <italic>gid</italic>
-E173* may facilitate the expansion of <italic>katG</italic>
-S315T mutant strains. Indeed, other groups have reported that SM-resistant strains seem to be more clustered in Vietnam<sup><xref ref-type="bibr" rid="CR48">48</xref>
</sup>
.</p>
<p id="Par29">Efflux pumps, including <italic>emrB</italic>
, have been reported to be associated with pathogenicity since the up-regulation of efflux gene expression is involved in the development of resistance to anti-TB drugs<sup><xref ref-type="bibr" rid="CR49">49</xref>
</sup>
and a wide array of physiologic processes such as the growth kinetics or transportation of a variety of compounds<sup><xref ref-type="bibr" rid="CR50">50</xref>
</sup>
. The combination of the <italic>katG</italic>
-S315T mutant with <italic>emrB</italic>
variant F508S in the Hanoi study may thus increase the drug efflux activity and facilitate Mtb survival and spread, by mitigating drug pressure. <italic>emrB</italic>
(<italic>Rv0783c</italic>
) belongs to the major facilitator superfamily (MFS) characteristically energized by the proton motive force (H<sup>+</sup>
or Na<sup>+</sup>
)<sup><xref ref-type="bibr" rid="CR50">50</xref>
,<xref ref-type="bibr" rid="CR51">51</xref>
</sup>
, and may confer low-level resistance to RMP<sup><xref ref-type="bibr" rid="CR52">52</xref>
</sup>
. Although the role of I461I in <italic>emrB</italic>
in KwaZulu-Natal’s population is unknown at present, this synonymous mutation (c.1383 C > A, ATC > ATA) is very rare in terms of codon usage of Mtb<sup><xref ref-type="bibr" rid="CR53">53</xref>
</sup>
, and the significance of codon usage bias and t-RNA modification should be taken into account, because rare codons are sometimes advantageous to the survival of Mtb under stress conditions<sup><xref ref-type="bibr" rid="CR54">54</xref>
</sup>
. When exposed to INH, various MFS efflux pump genes were reported to be overexpressed<sup><xref ref-type="bibr" rid="CR50">50</xref>
</sup>
, and these may induce sustained increased efflux activity with selection and stabilization of drug-resistant mutations<sup><xref ref-type="bibr" rid="CR55">55</xref>
</sup>
. This may also be relevant to the acquisition of additional drug resistance. Further studies are required for elucidating the function of the mutations in efflux pump genes. Indeed, SNVs in efflux pump genes are often found in XDR strains but not in drug-susceptible strains<sup><xref ref-type="bibr" rid="CR51">51</xref>
</sup>
, although it is often difficult to identify the extrusion of a drug to a specific gene<sup><xref ref-type="bibr" rid="CR50">50</xref>
,<xref ref-type="bibr" rid="CR56">56</xref>
</sup>
. The association pattern characterized by the presence or absence of the variant of <italic>emrB</italic>
with the phenotype was quite similar to that of <italic>PPE18/19</italic>
, whose phenomenon was shown beyond lineages, L2 in Hanoi, and L4 in KwaZulu-Natal. Although membrane localization of the <italic>PPE</italic>
genes may be functionally linked to efflux pump activities, it is currently unknown.</p>
<p id="Par30"><italic>Rv1588c</italic>
is a partial REP13E12 repeat protein<sup><xref ref-type="bibr" rid="CR57">57</xref>
</sup>
. Although k-mers carrying the reference sequences in <italic>Rv1588c</italic>
showed negative association with the clustered strains harboring <italic>katG</italic>
-S315T in the two panels, their functional significance remains unclear. As a variant, only a synonymous variant P34P was found in Hanoi, which was associated with ancient Beijing sublineage.</p>
<p id="Par31">The reference sequence (=wild type) k-mers in <italic>pncA</italic>
were also associated negatively with the phenotype in the Hanoi cohort, and variant-carrying k-mers showed a positive association in the KwaZulu-Natal cohort. As a variant found in Hanoi, Q141P in <italic>pncA</italic>
has been reported as a high-confidence mutation leading to PZA resistance<sup><xref ref-type="bibr" rid="CR58">58</xref>
,<xref ref-type="bibr" rid="CR59">59</xref>
</sup>
. PZA resistance is often observed among MDR-TB isolates<sup><xref ref-type="bibr" rid="CR60">60</xref>
,<xref ref-type="bibr" rid="CR61">61</xref>
</sup>
. Thus, the possibility of <italic>pncA</italic>
mutations facilitating the transmission of <italic>katG</italic>
-S315T mutant Mtb isolates may make TB management more challenging.</p>
<p id="Par32">Convergence-based phyOverlap analysis, which is a different approach, revealed that four different gene/intergenic regions only in the Hanoi study population that may have been caused by convergence evolution were significantly associated with the clustered strains carrying the <italic>katG</italic>
-S315T mutation; and these four were present in the <italic>PE or PPE</italic>
region. The impact of genetic variation on the function of PE or PPE proteins remains largely unknown<sup><xref ref-type="bibr" rid="CR35">35</xref>
</sup>
. However, at least large deletions between <italic>PPE46</italic>
and <italic>PPE47</italic>
genes were observed in all three lineages (L1, L2, and L4) in Hanoi, and these were positively associated with the spread of INH-resistant TB. The deletion between the identical sequences of the two <italic>PPE</italic>
genes leads to in-frame gene fusions through homologous recombination<sup><xref ref-type="bibr" rid="CR62">62</xref>
</sup>
, and a relatively high prevalence, indicating a clonal expansion of Haarlem strains (L4) in Tunisia<sup><xref ref-type="bibr" rid="CR34">34</xref>
</sup>
, suggested that the generation of the new chimeric genes may facilitate antigenic diversity and provide new determinants for pathogen’s virulence across the lineages.</p>
<p id="Par33">Further analyses using logistic regression models confirmed that all variants corresponding to significant k-mers of the five genes and even variants detected by the phylogenetic convergence test in Hanoi samples were positively associated with the spread of INH-resistant TB, even after adjustment for other possible confounders.</p>
<p id="Par34">Our study has some limitations. First, we were not able to trace the epidemiological link among the patients to corroborate the transmission chain. Hanoi is the capital city of Vietnam with on-going urbanization and this city consists of a floating population coming from many other provinces; thus, pursuing an epidemiological link is rather difficult. However, we have detailed information on the patients’ residential districts and we have added this information to the logistic regression analyses. Second, our samples were obtained in a population-based setting in an Asian city; but to generalize our finding, we analyzed another African cohort set available in the public by using the same methodology. Third, performing <italic>in vitro</italic>
experiments to elucidate the functional significance of each genetic variant was beyond our scope owing to resource limitation. Nevertheless, the extracted genes were associated with the spread of INH-resistant strains carrying <italic>katG</italic>
-S315T mutation and these reached statistically significant levels by using the bacterial GWAS approach based on LMM.</p>
<p id="Par35">Previous studies suggest that the <italic>katG</italic>
gene’s physiological function is not largely reduced by S315T substitution<sup><xref ref-type="bibr" rid="CR63">63</xref>
</sup>
. Its catalase-peroxidase-peroxynitritase activities may play a part to protect Mtb against reactive oxygen and nitrogen species derived from the phagocyte oxidative burst in human macrophages<sup><xref ref-type="bibr" rid="CR63">63</xref>
,<xref ref-type="bibr" rid="CR64">64</xref>
</sup>
. This may link KatG with other pathogen factors relevant to immune evasion or virulence such as PE/PPE<sup><xref ref-type="bibr" rid="CR65">65</xref>
,<xref ref-type="bibr" rid="CR66">66</xref>
</sup>
, although possible additive or synergistic effects on fitness should be further investigated. It is desirable to conduct validation studies in different populations. These Mtb genes are attractive candidates, presumably because of their relevance to the pathogen’s virulence, and they could be important sources to consider in <italic>in vitro</italic>
and in animal models.</p>
<p id="Par36">In conclusion, WGS data demonstrated the status of primary drug resistance at gene levels in the Hanoi city, and bacterial GWAS was performed to identify candidate genes that may facilitate the spreading of INH-resistant strains. Our findings provide new insights into the pathogenic mechanisms possibly mediated by the candidate genes including <italic>PE/PPE</italic>
, by which drug-resistant Mtb can maintain epidemiological fitness and spread in high-burden countries such as in Asia and Africa.</p>
</sec>
<sec id="Sec11"><title>Methods</title>
<sec id="Sec12"><title>Study sites, patient recruitment, and sample collection</title>
<p id="Par37">This was a part of our cohort study of patients who were over 16 years of age and who were newly diagnosed with smear-positive pulmonary TB without any treatment history in Hanoi, Vietnam during 2007–2009, in which basic data with clinical interpretation were published in a previous report<sup><xref ref-type="bibr" rid="CR12">12</xref>
,<xref ref-type="bibr" rid="CR67">67</xref>
–<xref ref-type="bibr" rid="CR69">69</xref>
</sup>
. In brief, we included 7 of the 14 districts in Hanoi as the catchment area, where more than half of new smear-positive TB patients in the city were diagnosed and treated in the area during the study period<sup><xref ref-type="bibr" rid="CR12">12</xref>
</sup>
. Sputum specimens were collected before starting the treatment, 92.7% of which revealed culture-positive, and drug susceptibility testing for first-line drugs was performed using the WHO standard proportional method<sup><xref ref-type="bibr" rid="CR12">12</xref>
</sup>
. The patients’ clinico-epidemiological information was also collected. Written informed consent was obtained from all the patients.</p>
</sec>
<sec id="Sec13"><title>Ethics statement</title>
<p id="Par38">This cohort study was approved by the Ethical Committees of the Ministry of Health, Vietnam, National Center for Global Health and Medicine, and the Research Institute of Tuberculosis, Japan Anti-Tuberculosis Association, Japan. All experiments were performed in accordance with relevant guidelines and regulations. In the case of minors, their parents provided written informed consent.</p>
</sec>
<sec id="Sec14"><title>WGS</title>
<p id="Par39">Mycobacterial DNA samples from Löwenstein-Jensen culture media were extracted using the Isoplant kit (Nippon Gene, Tokyo, Japan) and analyzed using Illumina HiSeq and MiSeq systems (Illumina, San Diego, CA, USA). These experiments were performed using a class II safety cabinet in a biosafety level 3 laboratory to prevent contamination. For Hiseq. 2500, a library of WGS was prepared using an automated sample preparation system (Agilent Technologies Inc.) with the TruSeq DNA PCR free sample prep kit (Illumina). For Miseq, a library was prepared from 200 ng of genomic DNA with the TruSeq Nano DNA LT Sample Preparation Kit (Illumina), following the manufacturer’s instructions. Paired-end (2 × 150 bp) sequencing was performed using Hiseq. 2500. For Miseq, the paired-end (2 × 250 bp or 2 × 300 bp) sequencing system was used. The sequence data are available in the DDBJ/EMBL/GenBank databases under the accession numbers DRA008666-7 and DRA008677.</p>
</sec>
<sec id="Sec15"><title>Extracting single nucleotide variants (SNVs) for Mtb lineages/sublineages and genetic clustering</title>
<p id="Par40">Briefly, after trimming and excluding severely contaminated samples, sequence reads were mapped to the H37Rv genome (AL123456.3) by using BWA-MEM 0.7.15 (<ext-link ext-link-type="uri" xlink:href="https://github.com/lh3/bwa">https://github.com/lh3/bwa</ext-link>
), followed by variant calling with the Genome analysis toolkit (GATK version 3.7)<sup><xref ref-type="bibr" rid="CR70">70</xref>
</sup>
. Only paired-end fastq files with average depth more than 25X were accepted for the analysis. The criteria set for identifying SNVs and small indels included Q30 minimum base call quality score and a minimum coverage depth of 10X. Drug resistance-conferring mutations, including small indels and lineage-specific variations, were extracted using the TB-Profiler version 0.3.7<sup><xref ref-type="bibr" rid="CR9">9</xref>
,<xref ref-type="bibr" rid="CR22">22</xref>
</sup>
. The Beijing genotype of lineage-2 (L2.2) was further classified into ancient and modern Beijing sublineages by detecting the SNVs at the nucleotide position 649,345, which is equivalent to the presence of IS6110 in the NTF region<sup><xref ref-type="bibr" rid="CR71">71</xref>
</sup>
.</p>
<p id="Par41">Large deletions were screened throughout the mapped reads by seeking zero or low-depth areas (<15% of the average depth) using an in-house python script and then visualized for confirmation with the Integrative Genomics Viewer (IGV) version 2.3.91. For this deletion screening, complete genome sequences of the clinical isolates in our Hanoi cohort, AP018033 to AP018036<sup><xref ref-type="bibr" rid="CR72">72</xref>
,<xref ref-type="bibr" rid="CR73">73</xref>
</sup>
as well as H37Rv genome, were used as reference sequences. After excluding ambiguous variants in categories of repetitive and insertion sequences and phages, genetic clusters were defined by the pairwise differences of no more than five SNVs<sup><xref ref-type="bibr" rid="CR74">74</xref>
</sup>
. A phylogenetic tree was constructed by the maximum likelihood method using RAxML version 8.2.8<sup><xref ref-type="bibr" rid="CR75">75</xref>
</sup>
and then visualized with plotTree (<ext-link ext-link-type="uri" xlink:href="https://github.com/katholt/plotTree">https://github.com/katholt/plotTree</ext-link>
) by using a lineage-7 strain ERR181435 as an out-group.</p>
</sec>
<sec id="Sec16"><title>Analyses of bacterial GWAS</title>
<p id="Par42">The associations between the phenotypes and the presence or absence of 31-bp short sequences, k-mers, in the genome were investigated using a genome-wide efficient mixed model association algorithm, the GEMMA software (<ext-link ext-link-type="uri" xlink:href="https://github.com/genetics-statistics/GEMMA">https://github.com/genetics-statistics/GEMMA</ext-link>
). At first, the DSK software (<ext-link ext-link-type="uri" xlink:href="https://github.com/GATB/dsk">https://github.com/GATB/dsk</ext-link>
) was used for listing all the unique 31-bp DNA k-mers, and then their presence or absence in all the samples was analyzed as mentioned above<sup><xref ref-type="bibr" rid="CR20">20</xref>
</sup>
.</p>
<p id="Par43">DBGWAS, an extended k-mer-based GWAS tool with compacted De Bruijn graph<sup><xref ref-type="bibr" rid="CR76">76</xref>
</sup>
, was further used to confirm the genetic variants associated with the phenotypes of interest. Sequence reads were assembled using SPAdes v3.13.0<sup><xref ref-type="bibr" rid="CR77">77</xref>
</sup>
and Platanus 1.2.4<sup><xref ref-type="bibr" rid="CR78">78</xref>
</sup>
when appropriate, and the generated contigs were used for BLAST search (ncbi-blast-2.8.1+) to identify the location of the phenotype-associated k-mers. Bonferroni correction was applied for multiple testing; the threshold of the significance after correction was calculated as 0.05 divided by the number of variants identified in the study samples.</p>
<p id="Par44">To investigate whether any variant, including SNVs or indels in the whole genome, has any possible association with the phenotypes of interest, we used bugwas R package<sup><xref ref-type="bibr" rid="CR20">20</xref>
</sup>
with built-in GEMMA. The bi-allelic SNVs were used to calculate the relatedness matrix of the samples for LMMs to control for the population structure. Likelihood ratio tests were used for assessing the significant associations.</p>
<p id="Par45">Phylogenetic convergence tests (phyOverlap)<sup><xref ref-type="bibr" rid="CR26">26</xref>
</sup>
were also performed to identify the convergent variants associated with the phenotypes. Herein, Benjamini-Hochberg adjustment at 0.05 false discovery rate level was applied as the q value for phyOverlap.</p>
<p id="Par46">To compare the findings obtained from the GWAS analysis of our 332 samples in Vietnam, another set of WGS data from 337 clinical isolates, which were collected in the KwaZulu-Natal province of South Africa from 2008 to 2013 in a study conducted by Cohen <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR19">19</xref>
</sup>
to investigate the emergence of drug-resistant TB (hereafter referred to as KwaZulu-Natal study), were retrieved from the public database, and analyzed in a similar way.</p>
</sec>
<sec id="Sec17"><title>Other statistical analyses</title>
<p id="Par47">Chi square and Fisher’s exact tests were performed to compare the differences in the proportions among the groups. Venn diagram (R version 3.4.4 VennDiagram package) was used to demonstrate the common gene(s) harboring variants, including k-mers, associated with different phenotypes. Possible associations between the given genetic variations and INH-resistant clusters, adjusted for Mtb lineages and patients’ age, gender, and living area were further studied using logistic regression models. These analyses were performed using STATA version 12 (StataCorp, College Station, TX, USA), and P values less than 0.05 were considered statistically significant.</p>
</sec>
</sec>
<sec sec-type="supplementary-material"><title>Supplementary information</title>
<sec id="Sec18"><p><supplementary-material content-type="local-data" id="MOESM1"><media xlink:href="41598_2019_51812_MOESM1_ESM.pdf"><caption><p>Supplementary information</p>
</caption>
</media>
</supplementary-material>
</p>
</sec>
</sec>
</body>
<back><fn-group><fn><p><bold>Publisher’s note</bold>
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.</p>
</fn>
<fn><p>These authors contributed equally: Nguyen Thi Le Hang and Minako Hijikata.</p>
</fn>
</fn-group>
<sec><title>Supplementary information</title>
<p>is available for this paper at 10.1038/s41598-019-51812-7.</p>
</sec>
<ack><title>Acknowledgements</title>
<p>The authors would like to thank Dr. Luu Thi Lien for her leadership and Dr. Shinsaku Sakurada for his monitoring of the on-site implementation of the original cohort study, and Ms. Ikumi Matsushita for her technical support and resource contribution. The authors would like to thank Enago (<ext-link ext-link-type="uri" xlink:href="http://www.enago.jp">www.enago.jp</ext-link>
) for the English language review. This research was supported by AMED under Grant Number JP19fm0108001 (Japan Initiative for Global Research Network on Infectious Diseases; J-GRID). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</p>
</ack>
<notes notes-type="author-contribution"><title>Author contributions</title>
<p>N.T.L.H. supervised the on-site implementation, performed the analyses, and wrote and finalized the manuscript. M.H. performed whole genome sequencing and analyzed the sequencing results. S.M. provided technical support and performed genotyping. J.O. provided statistical advice and support. A.M. and S.S. provided technical support for the microbiological experiments. P.H.T., H.V.H. and V.C.C. supervised on-site sample and data collection. N.P.H. and N.V.H. supervised on-site bacterial sample management. N.K. conceptualized the project, wrote the in-house scripts for data analyses, performed the analyses, and corrected and finalized the manuscript.</p>
</notes>
<notes notes-type="data-availability"><title>Data availability</title>
<p>All data pertaining to the manuscript have been provided in the forms of tables and figures. Supporting Information is available as Supplementary Tables <xref rid="MOESM1" ref-type="media">S1</xref>
–<xref rid="MOESM1" ref-type="media">S6</xref>
and Figs <xref rid="MOESM1" ref-type="media">S1</xref>
–<xref rid="MOESM1" ref-type="media">S7</xref>
. Datasets pertaining to the sequence searches described here are available from the corresponding author on request.</p>
</notes>
<notes notes-type="COI-statement"><title>Competing interests</title>
<p id="Par48">The authors declare no competing interests.</p>
</notes>
<ref-list id="Bib1"><title>References</title>
<ref id="CR1"><label>1.</label>
<mixed-citation publication-type="other">WHO. World Health Organization. Global Tuberculosis Report 2018 (2018).</mixed-citation>
</ref>
<ref id="CR2"><label>2.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bastos</surname>
<given-names>ML</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Treatment outcomes of patients with multidrug-resistant and extensively drug-resistant tuberculosis according to drug susceptibility testing to first- and second-line drugs: an individual patient data meta-analysis</article-title>
<source>Clin Infect Dis</source>
<year>2014</year>
<volume>59</volume>
<fpage>1364</fpage>
<lpage>1374</lpage>
<pub-id pub-id-type="doi">10.1093/cid/ciu619</pub-id>
<pub-id pub-id-type="pmid">25097082</pub-id>
</element-citation>
</ref>
<ref id="CR3"><label>3.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pooran</surname>
<given-names>A</given-names>
</name>
<name><surname>Pieterson</surname>
<given-names>E</given-names>
</name>
<name><surname>Davids</surname>
<given-names>M</given-names>
</name>
<name><surname>Theron</surname>
<given-names>G</given-names>
</name>
<name><surname>Dheda</surname>
<given-names>K</given-names>
</name>
</person-group>
<article-title>What is the cost of diagnosis and management of drug resistant tuberculosis in South Africa?</article-title>
<source>PLoS One</source>
<year>2013</year>
<volume>8</volume>
<fpage>e54587</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pone.0054587</pub-id>
<pub-id pub-id-type="pmid">23349933</pub-id>
</element-citation>
</ref>
<ref id="CR4"><label>4.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Muller</surname>
<given-names>B</given-names>
</name>
<name><surname>Borrell</surname>
<given-names>S</given-names>
</name>
<name><surname>Rose</surname>
<given-names>G</given-names>
</name>
<name><surname>Gagneux</surname>
<given-names>S</given-names>
</name>
</person-group>
<article-title>The heterogeneous evolution of multidrug-resistant Mycobacterium tuberculosis</article-title>
<source>Trends Genet</source>
<year>2013</year>
<volume>29</volume>
<fpage>160</fpage>
<lpage>169</lpage>
<pub-id pub-id-type="doi">10.1016/j.tig.2012.11.005</pub-id>
<pub-id pub-id-type="pmid">23245857</pub-id>
</element-citation>
</ref>
<ref id="CR5"><label>5.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Manson</surname>
<given-names>AL</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Genomic analysis of globally diverse Mycobacterium tuberculosis strains provides insights into the emergence and spread of multidrug resistance</article-title>
<source>Nat Genet</source>
<year>2017</year>
<volume>49</volume>
<fpage>395</fpage>
<lpage>402</lpage>
<pub-id pub-id-type="doi">10.1038/ng.3767</pub-id>
<pub-id pub-id-type="pmid">28092681</pub-id>
</element-citation>
</ref>
<ref id="CR6"><label>6.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Menzies</surname>
<given-names>Dick</given-names>
</name>
<name><surname>Benedetti</surname>
<given-names>Andrea</given-names>
</name>
<name><surname>Paydar</surname>
<given-names>Anita</given-names>
</name>
<name><surname>Royce</surname>
<given-names>Sarah</given-names>
</name>
<name><surname>Pai</surname>
<given-names>Madhukar</given-names>
</name>
<name><surname>Burman</surname>
<given-names>William</given-names>
</name>
<name><surname>Vernon</surname>
<given-names>Andrew</given-names>
</name>
<name><surname>Lienhardt</surname>
<given-names>Christian</given-names>
</name>
</person-group>
<article-title>Standardized Treatment of Active Tuberculosis in Patients with Previous Treatment and/or with Mono-resistance to Isoniazid: A Systematic Review and Meta-analysis</article-title>
<source>PLoS Medicine</source>
<year>2009</year>
<volume>6</volume>
<issue>9</issue>
<fpage>e1000150</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pmed.1000150</pub-id>
<pub-id pub-id-type="pmid">20101802</pub-id>
</element-citation>
</ref>
<ref id="CR7"><label>7.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Stagg</surname>
<given-names>HR</given-names>
</name>
<name><surname>Lipman</surname>
<given-names>MC</given-names>
</name>
<name><surname>McHugh</surname>
<given-names>TD</given-names>
</name>
<name><surname>Jenkins</surname>
<given-names>HE</given-names>
</name>
</person-group>
<article-title>Isoniazid-resistant tuberculosis: a cause for concern?</article-title>
<source>Int J Tuberc Lung Dis</source>
<year>2017</year>
<volume>21</volume>
<fpage>129</fpage>
<lpage>139</lpage>
<pub-id pub-id-type="doi">10.5588/ijtld.16.0716</pub-id>
<pub-id pub-id-type="pmid">28234075</pub-id>
</element-citation>
</ref>
<ref id="CR8"><label>8.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nebenzahl-Guimaraes</surname>
<given-names>H</given-names>
</name>
<name><surname>Jacobson</surname>
<given-names>KR</given-names>
</name>
<name><surname>Farhat</surname>
<given-names>MR</given-names>
</name>
<name><surname>Murray</surname>
<given-names>MB</given-names>
</name>
</person-group>
<article-title>Systematic review of allelic exchange experiments aimed at identifying mutations that confer drug resistance in Mycobacterium tuberculosis</article-title>
<source>J Antimicrob Chemother</source>
<year>2014</year>
<volume>69</volume>
<fpage>331</fpage>
<lpage>342</lpage>
<pub-id pub-id-type="doi">10.1093/jac/dkt358</pub-id>
<pub-id pub-id-type="pmid">24055765</pub-id>
</element-citation>
</ref>
<ref id="CR9"><label>9.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Coll</surname>
<given-names>F</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Rapid determination of anti-tuberculosis drug resistance from whole-genome sequences</article-title>
<source>Genome Med</source>
<year>2015</year>
<volume>7</volume>
<fpage>51</fpage>
<pub-id pub-id-type="doi">10.1186/s13073-015-0164-0</pub-id>
<pub-id pub-id-type="pmid">26019726</pub-id>
</element-citation>
</ref>
<ref id="CR10"><label>10.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Farhat</surname>
<given-names>MR</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Genomic analysis identifies targets of convergent positive selection in drug-resistant Mycobacterium tuberculosis</article-title>
<source>Nat Genet</source>
<year>2013</year>
<volume>45</volume>
<fpage>1183</fpage>
<lpage>1189</lpage>
<pub-id pub-id-type="doi">10.1038/ng.2747</pub-id>
<pub-id pub-id-type="pmid">23995135</pub-id>
</element-citation>
</ref>
<ref id="CR11"><label>11.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nhung</surname>
<given-names>NV</given-names>
</name>
<name><surname>Hoa</surname>
<given-names>NB</given-names>
</name>
<name><surname>Sy</surname>
<given-names>DN</given-names>
</name>
<name><surname>Hennig</surname>
<given-names>CM</given-names>
</name>
<name><surname>Dean</surname>
<given-names>AS</given-names>
</name>
</person-group>
<article-title>The fourth national anti-tuberculosis drug resistance survey in Viet Nam</article-title>
<source>Int J Tuberc Lung Dis</source>
<year>2015</year>
<volume>19</volume>
<fpage>670</fpage>
<lpage>675</lpage>
<pub-id pub-id-type="doi">10.5588/ijtld.14.0785</pub-id>
<pub-id pub-id-type="pmid">25946357</pub-id>
</element-citation>
</ref>
<ref id="CR12"><label>12.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hang</surname>
<given-names>NT</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Primary drug-resistant tuberculosis in Hanoi, Viet Nam: present status and risk factors</article-title>
<source>PLoS One</source>
<year>2013</year>
<volume>8</volume>
<fpage>e71867</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pone.0071867</pub-id>
<pub-id pub-id-type="pmid">23967255</pub-id>
</element-citation>
</ref>
<ref id="CR13"><label>13.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Minh</surname>
<given-names>NN</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Molecular characteristics of rifampin- and isoniazid-resistant mycobacterium tuberculosis strains isolated in Vietnam</article-title>
<source>J Clin Microbiol</source>
<year>2012</year>
<volume>50</volume>
<fpage>598</fpage>
<lpage>601</lpage>
<pub-id pub-id-type="doi">10.1128/JCM.05171-11</pub-id>
<pub-id pub-id-type="pmid">22170905</pub-id>
</element-citation>
</ref>
<ref id="CR14"><label>14.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nhu</surname>
<given-names>NT</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Association of streptomycin resistance mutations with level of drug resistance and Mycobacterium tuberculosis genotypes</article-title>
<source>Int J Tuberc Lung Dis</source>
<year>2012</year>
<volume>16</volume>
<fpage>527</fpage>
<lpage>531</lpage>
<pub-id pub-id-type="doi">10.5588/ijtld.11.0202</pub-id>
<pub-id pub-id-type="pmid">22640514</pub-id>
</element-citation>
</ref>
<ref id="CR15"><label>15.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nguyen</surname>
<given-names>HQ</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Quadruple-first line drug resistance in Mycobacterium tuberculosis in Vietnam: What can we learn from genes?</article-title>
<source>Infect Genet Evol</source>
<year>2017</year>
<volume>50</volume>
<fpage>55</fpage>
<lpage>61</lpage>
<pub-id pub-id-type="doi">10.1016/j.meegid.2017.02.012</pub-id>
<pub-id pub-id-type="pmid">28214557</pub-id>
</element-citation>
</ref>
<ref id="CR16"><label>16.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Papaventsis</surname>
<given-names>D</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Whole genome sequencing of Mycobacterium tuberculosis for detection of drug resistance: a systematic review</article-title>
<source>Clin Microbiol Infect</source>
<year>2017</year>
<volume>23</volume>
<fpage>61</fpage>
<lpage>68</lpage>
<pub-id pub-id-type="doi">10.1016/j.cmi.2016.09.008</pub-id>
<pub-id pub-id-type="pmid">27665704</pub-id>
</element-citation>
</ref>
<ref id="CR17"><label>17.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname>
<given-names>H</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Genome sequencing of 161 Mycobacterium tuberculosis isolates from China identifies genes and intergenic regions associated with drug resistance</article-title>
<source>Nat Genet</source>
<year>2013</year>
<volume>45</volume>
<fpage>1255</fpage>
<lpage>1260</lpage>
<pub-id pub-id-type="doi">10.1038/ng.2735</pub-id>
<pub-id pub-id-type="pmid">23995137</pub-id>
</element-citation>
</ref>
<ref id="CR18"><label>18.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Casali</surname>
<given-names>N</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Evolution and transmission of drug-resistant tuberculosis in a Russian population</article-title>
<source>Nat Genet</source>
<year>2014</year>
<volume>46</volume>
<fpage>279</fpage>
<lpage>286</lpage>
<pub-id pub-id-type="doi">10.1038/ng.2878</pub-id>
<pub-id pub-id-type="pmid">24464101</pub-id>
</element-citation>
</ref>
<ref id="CR19"><label>19.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cohen</surname>
<given-names>KA</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Evolution of Extensively Drug-Resistant Tuberculosis over Four Decades: Whole Genome Sequencing and Dating Analysis of Mycobacterium tuberculosis Isolates from KwaZulu-Natal</article-title>
<source>PLoS Med</source>
<year>2015</year>
<volume>12</volume>
<fpage>e1001880</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pmed.1001880</pub-id>
<pub-id pub-id-type="pmid">26418737</pub-id>
</element-citation>
</ref>
<ref id="CR20"><label>20.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Earle</surname>
<given-names>SG</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Identifying lineage effects when controlling for population structure improves power in bacterial association studies</article-title>
<source>Nat Microbiol</source>
<year>2016</year>
<volume>1</volume>
<fpage>16041</fpage>
<pub-id pub-id-type="doi">10.1038/nmicrobiol.2016.41</pub-id>
<pub-id pub-id-type="pmid">27572646</pub-id>
</element-citation>
</ref>
<ref id="CR21"><label>21.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Power</surname>
<given-names>RA</given-names>
</name>
<name><surname>Parkhill</surname>
<given-names>J</given-names>
</name>
<name><surname>de Oliveira</surname>
<given-names>T</given-names>
</name>
</person-group>
<article-title>Microbial genome-wide association studies: lessons from human GWAS</article-title>
<source>Nat Rev Genet</source>
<year>2017</year>
<volume>18</volume>
<fpage>41</fpage>
<lpage>50</lpage>
<pub-id pub-id-type="doi">10.1038/nrg.2016.132</pub-id>
<pub-id pub-id-type="pmid">27840430</pub-id>
</element-citation>
</ref>
<ref id="CR22"><label>22.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Coll</surname>
<given-names>F</given-names>
</name>
<etal></etal>
</person-group>
<article-title>A robust SNP barcode for typing Mycobacterium tuberculosis complex strains</article-title>
<source>Nat Commun</source>
<year>2014</year>
<volume>5</volume>
<fpage>4812</fpage>
<pub-id pub-id-type="doi">10.1038/ncomms5812</pub-id>
<pub-id pub-id-type="pmid">25176035</pub-id>
</element-citation>
</ref>
<ref id="CR23"><label>23.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gagneux</surname>
<given-names>S</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Impact of bacterial genetics on the transmission of isoniazid-resistant Mycobacterium tuberculosis</article-title>
<source>PLoS Pathog</source>
<year>2006</year>
<volume>2</volume>
<fpage>e61</fpage>
<pub-id pub-id-type="doi">10.1371/journal.ppat.0020061</pub-id>
<pub-id pub-id-type="pmid">16789833</pub-id>
</element-citation>
</ref>
<ref id="CR24"><label>24.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Choi</surname>
<given-names>Y</given-names>
</name>
<name><surname>Chan</surname>
<given-names>AP</given-names>
</name>
</person-group>
<article-title>PROVEAN web server: a tool to predict the functional effect of amino acid substitutions and indels</article-title>
<source>Bioinformatics</source>
<year>2015</year>
<volume>31</volume>
<fpage>2745</fpage>
<lpage>2747</lpage>
<pub-id pub-id-type="doi">10.1093/bioinformatics/btv195</pub-id>
<pub-id pub-id-type="pmid">25851949</pub-id>
</element-citation>
</ref>
<ref id="CR25"><label>25.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vaser</surname>
<given-names>R</given-names>
</name>
<name><surname>Adusumalli</surname>
<given-names>S</given-names>
</name>
<name><surname>Leng</surname>
<given-names>SN</given-names>
</name>
<name><surname>Sikic</surname>
<given-names>M</given-names>
</name>
<name><surname>Ng</surname>
<given-names>PC</given-names>
</name>
</person-group>
<article-title>SIFT missense predictions for genomes</article-title>
<source>Nat Protoc</source>
<year>2016</year>
<volume>11</volume>
<fpage>1</fpage>
<lpage>9</lpage>
<pub-id pub-id-type="doi">10.1038/nprot.2015.123</pub-id>
<pub-id pub-id-type="pmid">26633127</pub-id>
</element-citation>
</ref>
<ref id="CR26"><label>26.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hicks</surname>
<given-names>ND</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Clinically prevalent mutations in Mycobacterium tuberculosis alter propionate metabolism and mediate multidrug tolerance</article-title>
<source>Nat Microbiol</source>
<year>2018</year>
<volume>3</volume>
<fpage>1032</fpage>
<lpage>1042</lpage>
<pub-id pub-id-type="doi">10.1038/s41564-018-0218-3</pub-id>
<pub-id pub-id-type="pmid">30082724</pub-id>
</element-citation>
</ref>
<ref id="CR27"><label>27.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Huyen</surname>
<given-names>MN</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Epidemiology of isoniazid resistance mutations and their effect on tuberculosis treatment outcomes</article-title>
<source>Antimicrob Agents Chemother</source>
<year>2013</year>
<volume>57</volume>
<fpage>3620</fpage>
<lpage>3627</lpage>
<pub-id pub-id-type="doi">10.1128/AAC.00077-13</pub-id>
<pub-id pub-id-type="pmid">23689727</pub-id>
</element-citation>
</ref>
<ref id="CR28"><label>28.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Seifert</surname>
<given-names>M</given-names>
</name>
<name><surname>Catanzaro</surname>
<given-names>D</given-names>
</name>
<name><surname>Catanzaro</surname>
<given-names>A</given-names>
</name>
<name><surname>Rodwell</surname>
<given-names>TC</given-names>
</name>
</person-group>
<article-title>Genetic mutations associated with isoniazid resistance in Mycobacterium tuberculosis: a systematic review</article-title>
<source>PLoS One</source>
<year>2015</year>
<volume>10</volume>
<fpage>e0119628</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pone.0119628</pub-id>
<pub-id pub-id-type="pmid">25799046</pub-id>
</element-citation>
</ref>
<ref id="CR29"><label>29.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ameeruddin</surname>
<given-names>NU</given-names>
</name>
<name><surname>Luke Elizabeth</surname>
<given-names>H</given-names>
</name>
</person-group>
<article-title>Impact of isoniazid resistance on virulence of global and south Indian clinical isolates of Mycobacterium tuberculosis</article-title>
<source>Tuberculosis (Edinb)</source>
<year>2014</year>
<volume>94</volume>
<fpage>557</fpage>
<lpage>563</lpage>
<pub-id pub-id-type="doi">10.1016/j.tube.2014.08.011</pub-id>
<pub-id pub-id-type="pmid">25270728</pub-id>
</element-citation>
</ref>
<ref id="CR30"><label>30.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ano</surname>
<given-names>H</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Relationship between the isoniazid-resistant mutation katGS315T and the prevalence of MDR-/XDR-TB in Osaka, Japan</article-title>
<source>Int J Tuberc Lung Dis</source>
<year>2008</year>
<volume>12</volume>
<fpage>1300</fpage>
<lpage>1305</lpage>
<pub-id pub-id-type="pmid">18926041</pub-id>
</element-citation>
</ref>
<ref id="CR31"><label>31.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Unissa</surname>
<given-names>AN</given-names>
</name>
<name><surname>Selvakumar</surname>
<given-names>N</given-names>
</name>
<name><surname>Narayanan</surname>
<given-names>S</given-names>
</name>
<name><surname>Suganthi</surname>
<given-names>C</given-names>
</name>
<name><surname>Hanna</surname>
<given-names>LE</given-names>
</name>
</person-group>
<article-title>Investigation of Ser315 substitutions within katG gene in isoniazid-resistant clinical isolates of Mycobacterium tuberculosis from south India</article-title>
<source>Biomed Res Int</source>
<year>2015</year>
<volume>2015</volume>
<fpage>257983</fpage>
<pub-id pub-id-type="doi">10.1155/2015/257983</pub-id>
<pub-id pub-id-type="pmid">25699262</pub-id>
</element-citation>
</ref>
<ref id="CR32"><label>32.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Walker</surname>
<given-names>TM</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Whole-genome sequencing to delineate Mycobacterium tuberculosis outbreaks: a retrospective observational study</article-title>
<source>The Lancet Infectious Diseases</source>
<year>2013</year>
<volume>13</volume>
<fpage>137</fpage>
<lpage>146</lpage>
<pub-id pub-id-type="doi">10.1016/s1473-3099(12)70277-3</pub-id>
<pub-id pub-id-type="pmid">23158499</pub-id>
</element-citation>
</ref>
<ref id="CR33"><label>33.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>McEvoy</surname>
<given-names>CR</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Comparative analysis of Mycobacterium tuberculosis pe and ppe genes reveals high sequence variation and an apparent absence of selective constraints</article-title>
<source>PLoS One</source>
<year>2012</year>
<volume>7</volume>
<fpage>e30593</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pone.0030593</pub-id>
<pub-id pub-id-type="pmid">22496726</pub-id>
</element-citation>
</ref>
<ref id="CR34"><label>34.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Karboul</surname>
<given-names>A</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Frequent homologous recombination events in Mycobacterium tuberculosis PE/PPE multigene families: potential role in antigenic variability</article-title>
<source>J Bacteriol</source>
<year>2008</year>
<volume>190</volume>
<fpage>7838</fpage>
<lpage>7846</lpage>
<pub-id pub-id-type="doi">10.1128/JB.00827-08</pub-id>
<pub-id pub-id-type="pmid">18820012</pub-id>
</element-citation>
</ref>
<ref id="CR35"><label>35.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Fishbein</surname>
<given-names>S</given-names>
</name>
<name><surname>van Wyk</surname>
<given-names>N</given-names>
</name>
<name><surname>Warren</surname>
<given-names>RM</given-names>
</name>
<name><surname>Sampson</surname>
<given-names>SL</given-names>
</name>
</person-group>
<article-title>Phylogeny to function: PE/PPE protein evolution and impact on Mycobacterium tuberculosis pathogenicity</article-title>
<source>Mol Microbiol</source>
<year>2015</year>
<volume>96</volume>
<fpage>901</fpage>
<lpage>916</lpage>
<pub-id pub-id-type="doi">10.1111/mmi.12981</pub-id>
<pub-id pub-id-type="pmid">25727695</pub-id>
</element-citation>
</ref>
<ref id="CR36"><label>36.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nair</surname>
<given-names>S</given-names>
</name>
<name><surname>Pandey</surname>
<given-names>AD</given-names>
</name>
<name><surname>Mukhopadhyay</surname>
<given-names>S</given-names>
</name>
</person-group>
<article-title>The PPE18 protein of Mycobacterium tuberculosis inhibits NF-kappaB/rel-mediated proinflammatory cytokine production by upregulating and phosphorylating suppressor of cytokine signaling 3 protein</article-title>
<source>J Immunol</source>
<year>2011</year>
<volume>186</volume>
<fpage>5413</fpage>
<lpage>5424</lpage>
<pub-id pub-id-type="doi">10.4049/jimmunol.1000773</pub-id>
<pub-id pub-id-type="pmid">21451109</pub-id>
</element-citation>
</ref>
<ref id="CR37"><label>37.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bhat</surname>
<given-names>KH</given-names>
</name>
<name><surname>Ahmed</surname>
<given-names>A</given-names>
</name>
<name><surname>Kumar</surname>
<given-names>S</given-names>
</name>
<name><surname>Sharma</surname>
<given-names>P</given-names>
</name>
<name><surname>Mukhopadhyay</surname>
<given-names>S</given-names>
</name>
</person-group>
<article-title>Role of PPE18 protein in intracellular survival and pathogenicity of Mycobacterium tuberculosis in mice</article-title>
<source>PLoS One</source>
<year>2012</year>
<volume>7</volume>
<fpage>e52601</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pone.0052601</pub-id>
<pub-id pub-id-type="pmid">23300718</pub-id>
</element-citation>
</ref>
<ref id="CR38"><label>38.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dillon</surname>
<given-names>DC</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Molecular characterization and human T-cell responses to a member of a novel Mycobacterium tuberculosis mtb39 gene family</article-title>
<source>Infect Immun</source>
<year>1999</year>
<volume>67</volume>
<fpage>2941</fpage>
<lpage>2950</lpage>
<pub-id pub-id-type="pmid">10338503</pub-id>
</element-citation>
</ref>
<ref id="CR39"><label>39.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lewinsohn</surname>
<given-names>DA</given-names>
</name>
<name><surname>Lines</surname>
<given-names>RA</given-names>
</name>
<name><surname>Lewinsohn</surname>
<given-names>DM</given-names>
</name>
</person-group>
<article-title>Human dendritic cells presenting adenovirally expressed antigen elicit Mycobacterium tuberculosis–specific CD8+ T cells</article-title>
<source>Am J Respir Crit Care Med</source>
<year>2002</year>
<volume>166</volume>
<fpage>843</fpage>
<lpage>848</lpage>
<pub-id pub-id-type="doi">10.1164/rccm.2110094</pub-id>
<pub-id pub-id-type="pmid">12231495</pub-id>
</element-citation>
</ref>
<ref id="CR40"><label>40.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mortier</surname>
<given-names>MC</given-names>
</name>
<name><surname>Jongert</surname>
<given-names>E</given-names>
</name>
<name><surname>Mettens</surname>
<given-names>P</given-names>
</name>
<name><surname>Ruelle</surname>
<given-names>JL</given-names>
</name>
</person-group>
<article-title>Sequence conservation analysis and in silico human leukocyte antigen-peptide binding predictions for the Mtb72F and M72 tuberculosis candidate vaccine antigens</article-title>
<source>BMC Immunol</source>
<year>2015</year>
<volume>16</volume>
<fpage>63</fpage>
<pub-id pub-id-type="doi">10.1186/s12865-015-0119-7</pub-id>
<pub-id pub-id-type="pmid">26493839</pub-id>
</element-citation>
</ref>
<ref id="CR41"><label>41.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Reed</surname>
<given-names>S</given-names>
</name>
<name><surname>Lobet</surname>
<given-names>Y</given-names>
</name>
</person-group>
<article-title>Tuberculosis vaccine development; from mouse to man</article-title>
<source>Microbes Infect</source>
<year>2005</year>
<volume>7</volume>
<fpage>922</fpage>
<lpage>931</lpage>
<pub-id pub-id-type="doi">10.1016/j.micinf.2005.03.011</pub-id>
<pub-id pub-id-type="pmid">15935717</pub-id>
</element-citation>
</ref>
<ref id="CR42"><label>42.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rodo</surname>
<given-names>MJ</given-names>
</name>
<etal></etal>
</person-group>
<article-title>A comparison of antigen-specific T cell responses induced by six novel tuberculosis vaccine candidates</article-title>
<source>PLoS Pathog</source>
<year>2019</year>
<volume>15</volume>
<fpage>e1007643</fpage>
<pub-id pub-id-type="doi">10.1371/journal.ppat.1007643</pub-id>
<pub-id pub-id-type="pmid">30830940</pub-id>
</element-citation>
</ref>
<ref id="CR43"><label>43.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gaseitsiwe</surname>
<given-names>S</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Pattern recognition in pulmonary tuberculosis defined by high content peptide microarray chip analysis representing 61 proteins from M. tuberculosis</article-title>
<source>PLoS One</source>
<year>2008</year>
<volume>3</volume>
<fpage>e3840</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pone.0003840</pub-id>
<pub-id pub-id-type="pmid">19065269</pub-id>
</element-citation>
</ref>
<ref id="CR44"><label>44.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gaseitsiwe</surname>
<given-names>S</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Peptide microarray-based identification of Mycobacterium tuberculosis epitope binding to HLA-DRB1*0101, DRB1*1501, and DRB1*0401</article-title>
<source>Clin Vaccine Immunol</source>
<year>2010</year>
<volume>17</volume>
<fpage>168</fpage>
<lpage>175</lpage>
<pub-id pub-id-type="doi">10.1128/CVI.00208-09</pub-id>
<pub-id pub-id-type="pmid">19864486</pub-id>
</element-citation>
</ref>
<ref id="CR45"><label>45.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hebert</surname>
<given-names>AM</given-names>
</name>
<etal></etal>
</person-group>
<article-title>DNA polymorphisms in the pepA and PPE18 genes among clinical strains of Mycobacterium tuberculosis: implications for vaccine efficacy</article-title>
<source>Infect Immun</source>
<year>2007</year>
<volume>75</volume>
<fpage>5798</fpage>
<lpage>5805</lpage>
<pub-id pub-id-type="doi">10.1128/IAI.00335-07</pub-id>
<pub-id pub-id-type="pmid">17893137</pub-id>
</element-citation>
</ref>
<ref id="CR46"><label>46.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Okamoto</surname>
<given-names>S</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Loss of a conserved 7-methylguanosine modification in 16S rRNA confers low-level streptomycin resistance in bacteria</article-title>
<source>Mol Microbiol</source>
<year>2007</year>
<volume>63</volume>
<fpage>1096</fpage>
<lpage>1106</lpage>
<pub-id pub-id-type="doi">10.1111/j.1365-2958.2006.05585.x</pub-id>
<pub-id pub-id-type="pmid">17238915</pub-id>
</element-citation>
</ref>
<ref id="CR47"><label>47.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wong</surname>
<given-names>SY</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Functional role of methylation of G518 of the 16S rRNA 530 loop by GidB in Mycobacterium tuberculosis</article-title>
<source>Antimicrob Agents Chemother</source>
<year>2013</year>
<volume>57</volume>
<fpage>6311</fpage>
<lpage>6318</lpage>
<pub-id pub-id-type="doi">10.1128/AAC.00905-13</pub-id>
<pub-id pub-id-type="pmid">24100503</pub-id>
</element-citation>
</ref>
<ref id="CR48"><label>48.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Buu</surname>
<given-names>TN</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Increased transmission of Mycobacterium tuberculosis Beijing genotype strains associated with resistance to streptomycin: a population-based study</article-title>
<source>PLoS One</source>
<year>2012</year>
<volume>7</volume>
<fpage>e42323</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pone.0042323</pub-id>
<pub-id pub-id-type="pmid">22912700</pub-id>
</element-citation>
</ref>
<ref id="CR49"><label>49.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rodrigues</surname>
<given-names>L</given-names>
</name>
<name><surname>Parish</surname>
<given-names>T</given-names>
</name>
<name><surname>Balganesh</surname>
<given-names>M</given-names>
</name>
<name><surname>Ainsa</surname>
<given-names>JA</given-names>
</name>
</person-group>
<article-title>Antituberculosis drugs: reducing efflux=increasing activity</article-title>
<source>Drug Discov Today</source>
<year>2017</year>
<volume>22</volume>
<fpage>592</fpage>
<lpage>599</lpage>
<pub-id pub-id-type="doi">10.1016/j.drudis.2017.01.002</pub-id>
<pub-id pub-id-type="pmid">28089787</pub-id>
</element-citation>
</ref>
<ref id="CR50"><label>50.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Louw</surname>
<given-names>GE</given-names>
</name>
<etal></etal>
</person-group>
<article-title>A balancing act: efflux/influx in mycobacterial drug resistance</article-title>
<source>Antimicrob Agents Chemother</source>
<year>2009</year>
<volume>53</volume>
<fpage>3181</fpage>
<lpage>3189</lpage>
<pub-id pub-id-type="doi">10.1128/AAC.01577-08</pub-id>
<pub-id pub-id-type="pmid">19451293</pub-id>
</element-citation>
</ref>
<ref id="CR51"><label>51.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kanji</surname>
<given-names>A</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Single nucleotide polymorphisms in efflux pumps genes in extensively drug resistant Mycobacterium tuberculosis isolates from Pakistan</article-title>
<source>Tuberculosis (Edinb)</source>
<year>2017</year>
<volume>107</volume>
<fpage>20</fpage>
<lpage>30</lpage>
<pub-id pub-id-type="doi">10.1016/j.tube.2017.07.012</pub-id>
<pub-id pub-id-type="pmid">29050768</pub-id>
</element-citation>
</ref>
<ref id="CR52"><label>52.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pang</surname>
<given-names>Y</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Study of the rifampin monoresistance mechanism in Mycobacterium tuberculosis</article-title>
<source>Antimicrob Agents Chemother</source>
<year>2013</year>
<volume>57</volume>
<fpage>893</fpage>
<lpage>900</lpage>
<pub-id pub-id-type="doi">10.1128/AAC.01024-12</pub-id>
<pub-id pub-id-type="pmid">23208715</pub-id>
</element-citation>
</ref>
<ref id="CR53"><label>53.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gun</surname>
<given-names>L</given-names>
</name>
<name><surname>Yumiao</surname>
<given-names>R</given-names>
</name>
<name><surname>Haixian</surname>
<given-names>P</given-names>
</name>
<name><surname>Liang</surname>
<given-names>Z</given-names>
</name>
</person-group>
<article-title>Comprehensive Analysis and Comparison on the Codon Usage Pattern of Whole Mycobacterium tuberculosis Coding Genome from Different Area</article-title>
<source>Biomed Res Int</source>
<year>2018</year>
<volume>2018</volume>
<fpage>3574976</fpage>
<pub-id pub-id-type="doi">10.1155/2018/3574976</pub-id>
<pub-id pub-id-type="pmid">29854746</pub-id>
</element-citation>
</ref>
<ref id="CR54"><label>54.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chionh</surname>
<given-names>YH</given-names>
</name>
<etal></etal>
</person-group>
<article-title>tRNA-mediated codon-biased translation in mycobacterial hypoxic persistence</article-title>
<source>Nat Commun</source>
<year>2016</year>
<volume>7</volume>
<fpage>13302</fpage>
<pub-id pub-id-type="doi">10.1038/ncomms13302</pub-id>
<pub-id pub-id-type="pmid">27834374</pub-id>
</element-citation>
</ref>
<ref id="CR55"><label>55.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Machado</surname>
<given-names>D</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Contribution of efflux to the emergence of isoniazid and multidrug resistance in Mycobacterium tuberculosis</article-title>
<source>PLoS One</source>
<year>2012</year>
<volume>7</volume>
<fpage>e34538</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pone.0034538</pub-id>
<pub-id pub-id-type="pmid">22493700</pub-id>
</element-citation>
</ref>
<ref id="CR56"><label>56.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Garima</surname>
<given-names>K</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Differential expression of efflux pump genes of Mycobacterium tuberculosis in response to varied subinhibitory concentrations of antituberculosis agents</article-title>
<source>Tuberculosis (Edinb)</source>
<year>2015</year>
<volume>95</volume>
<fpage>155</fpage>
<lpage>161</lpage>
<pub-id pub-id-type="doi">10.1016/j.tube.2015.01.005</pub-id>
<pub-id pub-id-type="pmid">25680943</pub-id>
</element-citation>
</ref>
<ref id="CR57"><label>57.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gordon</surname>
<given-names>SV</given-names>
</name>
<name><surname>Heym</surname>
<given-names>B</given-names>
</name>
<name><surname>Parkhill</surname>
<given-names>J</given-names>
</name>
<name><surname>Barrell</surname>
<given-names>B</given-names>
</name>
<name><surname>Cole</surname>
<given-names>ST</given-names>
</name>
</person-group>
<article-title>New insertion sequences and a novel repeated sequence in the genome of Mycobacterium tuberculosis H37Rv</article-title>
<source>Microbiology</source>
<year>1999</year>
<volume>145</volume>
<issue>Pt 4</issue>
<fpage>881</fpage>
<lpage>892</lpage>
<pub-id pub-id-type="doi">10.1099/13500872-145-4-881</pub-id>
<pub-id pub-id-type="pmid">10220167</pub-id>
</element-citation>
</ref>
<ref id="CR58"><label>58.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chang</surname>
<given-names>KC</given-names>
</name>
<name><surname>Yew</surname>
<given-names>WW</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>Y</given-names>
</name>
</person-group>
<article-title>Pyrazinamide susceptibility testing in Mycobacterium tuberculosis: a systematic review with meta-analyses</article-title>
<source>Antimicrob Agents Chemother</source>
<year>2011</year>
<volume>55</volume>
<fpage>4499</fpage>
<lpage>4505</lpage>
<pub-id pub-id-type="doi">10.1128/AAC.00630-11</pub-id>
<pub-id pub-id-type="pmid">21768515</pub-id>
</element-citation>
</ref>
<ref id="CR59"><label>59.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Miotto</surname>
<given-names>P</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Mycobacterium tuberculosis pyrazinamide resistance determinants: a multicenter study</article-title>
<source>MBio</source>
<year>2014</year>
<volume>5</volume>
<fpage>e01819</fpage>
<lpage>01814</lpage>
<pub-id pub-id-type="doi">10.1128/mBio.01819-14</pub-id>
<pub-id pub-id-type="pmid">25336456</pub-id>
</element-citation>
</ref>
<ref id="CR60"><label>60.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Xia</surname>
<given-names>Q</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Phenotypic and genotypic characterization of pyrazinamide resistance among multidrug-resistant Mycobacterium tuberculosis isolates in Zhejiang, China</article-title>
<source>Antimicrob Agents Chemother</source>
<year>2015</year>
<volume>59</volume>
<fpage>1690</fpage>
<lpage>1695</lpage>
<pub-id pub-id-type="doi">10.1128/AAC.04541-14</pub-id>
<pub-id pub-id-type="pmid">25583712</pub-id>
</element-citation>
</ref>
<ref id="CR61"><label>61.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pang</surname>
<given-names>Y</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Prevalence and molecular characterization of pyrazinamide resistance among multidrug-resistant Mycobacterium tuberculosis isolates from Southern China</article-title>
<source>BMC Infect Dis</source>
<year>2017</year>
<volume>17</volume>
<fpage>711</fpage>
<pub-id pub-id-type="doi">10.1186/s12879-017-2761-6</pub-id>
<pub-id pub-id-type="pmid">29110640</pub-id>
</element-citation>
</ref>
<ref id="CR62"><label>62.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Marmiesse</surname>
<given-names>M</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Macro-array and bioinformatic analyses reveal mycobacterial ‘core’ genes, variation in the ESAT-6 gene family and new phylogenetic markers for the Mycobacterium tuberculosis complex</article-title>
<source>Microbiology</source>
<year>2004</year>
<volume>150</volume>
<fpage>483</fpage>
<lpage>496</lpage>
<pub-id pub-id-type="doi">10.1099/mic.0.26662-0</pub-id>
<pub-id pub-id-type="pmid">14766927</pub-id>
</element-citation>
</ref>
<ref id="CR63"><label>63.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Unissa</surname>
<given-names>AN</given-names>
</name>
<name><surname>Subbian</surname>
<given-names>S</given-names>
</name>
<name><surname>Hanna</surname>
<given-names>LE</given-names>
</name>
<name><surname>Selvakumar</surname>
<given-names>N</given-names>
</name>
</person-group>
<article-title>Overview on mechanisms of isoniazid action and resistance in Mycobacterium tuberculosis</article-title>
<source>Infect Genet Evol</source>
<year>2016</year>
<volume>45</volume>
<fpage>474</fpage>
<lpage>492</lpage>
<pub-id pub-id-type="doi">10.1016/j.meegid.2016.09.004</pub-id>
<pub-id pub-id-type="pmid">27612406</pub-id>
</element-citation>
</ref>
<ref id="CR64"><label>64.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ng</surname>
<given-names>VH</given-names>
</name>
<name><surname>Cox</surname>
<given-names>JS</given-names>
</name>
<name><surname>Sousa</surname>
<given-names>AO</given-names>
</name>
<name><surname>MacMicking</surname>
<given-names>JD</given-names>
</name>
<name><surname>McKinney</surname>
<given-names>JD</given-names>
</name>
</person-group>
<article-title>Role of KatG catalase-peroxidase in mycobacterial pathogenesis: countering the phagocyte oxidative burst</article-title>
<source>Mol Microbiol</source>
<year>2004</year>
<volume>52</volume>
<fpage>1291</fpage>
<lpage>1302</lpage>
<pub-id pub-id-type="doi">10.1111/j.1365-2958.2004.04078.x</pub-id>
<pub-id pub-id-type="pmid">15165233</pub-id>
</element-citation>
</ref>
<ref id="CR65"><label>65.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cui</surname>
<given-names>Ze-Jia</given-names>
</name>
<name><surname>Yang</surname>
<given-names>Qing-Yong</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>Hong-Yu</given-names>
</name>
<name><surname>Zhu</surname>
<given-names>Qiang</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>Qing-Ye</given-names>
</name>
</person-group>
<article-title>Bioinformatics Identification of Drug Resistance-Associated Gene Pairs in Mycobacterium tuberculosis</article-title>
<source>International Journal of Molecular Sciences</source>
<year>2016</year>
<volume>17</volume>
<issue>9</issue>
<fpage>1417</fpage>
<pub-id pub-id-type="doi">10.3390/ijms17091417</pub-id>
</element-citation>
</ref>
<ref id="CR66"><label>66.</label>
<mixed-citation publication-type="other">Bainomugisa, A. <italic>et al</italic>
. A complete high-quality MinION nanopore assembly of an extensively drug-resistant Mycobacterium tuberculosis Beijing lineage strain identifies novel variation in repetitive PE/PPE gene regions. <italic>Microb Genom</italic>
<bold>4</bold>
, 10.1099/mgen.0.000188 (2018).</mixed-citation>
</ref>
<ref id="CR67"><label>67.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hung</surname>
<given-names>NV</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Clonal expansion of Mycobacterium tuberculosis isolates and coexisting drug resistance in patients newly diagnosed with pulmonary tuberculosis in Hanoi, Vietnam</article-title>
<source>BMC research notes</source>
<year>2013</year>
<volume>6</volume>
<fpage>444</fpage>
<pub-id pub-id-type="doi">10.1186/1756-0500-6-444</pub-id>
<pub-id pub-id-type="pmid">24188178</pub-id>
</element-citation>
</ref>
<ref id="CR68"><label>68.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Maeda</surname>
<given-names>S</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Mycobacterium tuberculosis strains spreading in Hanoi, Vietnam: Beijing sublineages, genotypes, drug susceptibility patterns, and host factors</article-title>
<source>Tuberculosis (Edinb)</source>
<year>2014</year>
<volume>94</volume>
<fpage>649</fpage>
<lpage>656</lpage>
<pub-id pub-id-type="doi">10.1016/j.tube.2014.09.005</pub-id>
<pub-id pub-id-type="pmid">25459163</pub-id>
</element-citation>
</ref>
<ref id="CR69"><label>69.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hang</surname>
<given-names>NT</given-names>
</name>
<name><surname>Maeda</surname>
<given-names>S</given-names>
</name>
<name><surname>Keicho</surname>
<given-names>N</given-names>
</name>
<name><surname>Thuong</surname>
<given-names>PH</given-names>
</name>
<name><surname>Endo</surname>
<given-names>H</given-names>
</name>
</person-group>
<article-title>Sublineages of Mycobacterium tuberculosis Beijing genotype strains and unfavorable outcomes of anti-tuberculosis treatment</article-title>
<source>Tuberculosis (Edinb)</source>
<year>2015</year>
<volume>95</volume>
<fpage>336</fpage>
<lpage>342</lpage>
<pub-id pub-id-type="doi">10.1016/j.tube.2015.02.040</pub-id>
<pub-id pub-id-type="pmid">25732626</pub-id>
</element-citation>
</ref>
<ref id="CR70"><label>70.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hijikata</surname>
<given-names>M</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Spoligotyping and whole-genome sequencing analysis of lineage 1 strains of Mycobacterium tuberculosis in Da Nang, Vietnam</article-title>
<source>PLoS One</source>
<year>2017</year>
<volume>12</volume>
<fpage>e0186800</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pone.0186800</pub-id>
<pub-id pub-id-type="pmid">29049400</pub-id>
</element-citation>
</ref>
<ref id="CR71"><label>71.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chuang</surname>
<given-names>PC</given-names>
</name>
<name><surname>Chen</surname>
<given-names>HY</given-names>
</name>
<name><surname>Jou</surname>
<given-names>R</given-names>
</name>
</person-group>
<article-title>Single-nucleotide polymorphism in the fadD28 gene as a genetic marker for East Asia Lineage Mycobacterium tuberculosis</article-title>
<source>J Clin Microbiol</source>
<year>2010</year>
<volume>48</volume>
<fpage>4245</fpage>
<lpage>4247</lpage>
<pub-id pub-id-type="doi">10.1128/JCM.00970-10</pub-id>
<pub-id pub-id-type="pmid">20826639</pub-id>
</element-citation>
</ref>
<ref id="CR72"><label>72.</label>
<mixed-citation publication-type="other">Wada, T. <italic>et al</italic>
. Complete Genome Sequence of a Mycobacterium tuberculosis Strain Belonging to the East African-Indian Family in the Indo-Oceanic Lineage, Isolated in Hanoi, Vietnam. <italic>Genome Announc</italic>
<bold>5</bold>
, 10.1128/genomeA.00509-17 (2017).</mixed-citation>
</ref>
<ref id="CR73"><label>73.</label>
<mixed-citation publication-type="other">Wada, T. et al. Complete Genome Sequences of Three Representative Mycobacterium tuberculosis Beijing Family Strains Belonging to Distinct Genotype Clusters in Hanoi, Vietnam, during 2007 to 2009. <italic>Genome Announc</italic>
<bold>5</bold>
, 10.1128/genomeA.00510-17 (2017).</mixed-citation>
</ref>
<ref id="CR74"><label>74.</label>
<mixed-citation publication-type="other">Nikolayevskyy, V. <italic>et al</italic>
. Role and value of whole genome sequencing in studying tuberculosis transmission. <italic>Clin Microbiol Infect</italic>
, 10.1016/j.cmi.2019.03.022 (2019).</mixed-citation>
</ref>
<ref id="CR75"><label>75.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Stamatakis</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies</article-title>
<source>Bioinformatics</source>
<year>2014</year>
<volume>30</volume>
<fpage>1312</fpage>
<lpage>1313</lpage>
<pub-id pub-id-type="doi">10.1093/bioinformatics/btu033</pub-id>
<pub-id pub-id-type="pmid">24451623</pub-id>
</element-citation>
</ref>
<ref id="CR76"><label>76.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jaillard</surname>
<given-names>M</given-names>
</name>
<etal></etal>
</person-group>
<article-title>A fast and agnostic method for bacterial genome-wide association studies: Bridging the gap between k-mers and genetic events</article-title>
<source>PLoS Genet</source>
<year>2018</year>
<volume>14</volume>
<fpage>e1007758</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pgen.1007758</pub-id>
<pub-id pub-id-type="pmid">30419019</pub-id>
</element-citation>
</ref>
<ref id="CR77"><label>77.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bankevich</surname>
<given-names>A</given-names>
</name>
<etal></etal>
</person-group>
<article-title>SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing</article-title>
<source>J Comput Biol</source>
<year>2012</year>
<volume>19</volume>
<fpage>455</fpage>
<lpage>477</lpage>
<pub-id pub-id-type="doi">10.1089/cmb.2012.0021</pub-id>
<pub-id pub-id-type="pmid">22506599</pub-id>
</element-citation>
</ref>
<ref id="CR78"><label>78.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kajitani</surname>
<given-names>R</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Efficient de novo assembly of highly heterozygous genomes from whole-genome shotgun short reads</article-title>
<source>Genome Res</source>
<year>2014</year>
<volume>24</volume>
<fpage>1384</fpage>
<lpage>1395</lpage>
<pub-id pub-id-type="doi">10.1101/gr.170720.113</pub-id>
<pub-id pub-id-type="pmid">24755901</pub-id>
</element-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Sante/explor/MersV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000462 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000462 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Sante |area= MersV1 |flux= Pmc |étape= Corpus |type= RBID |clé= |texte= }}
![]() | This area was generated with Dilib version V0.6.33. | ![]() |