<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Med.</journal-id>
<journal-title>Frontiers in Medicine</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Med.</abbrev-journal-title>
<issn pub-type="epub">2296-858X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmed.2023.1292665</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Medicine</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Prolonged survival of a patient with active MDR-TB HIV co-morbidity: insights from a <italic>Mycobacterium tuberculosis</italic> strain with a unique genomic deletion</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" equal-contrib="yes">
<name>
<surname>Rubinstein</surname>
<given-names>Mor</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="author-notes" rid="fn0004"><sup>&#x2020;</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/951907/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author" equal-contrib="yes">
<name>
<surname>Makhon</surname>
<given-names>Andrei</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="author-notes" rid="fn0004"><sup>&#x2020;</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Losev</surname>
<given-names>Yelena</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/2354894/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Valenci</surname>
<given-names>Gal Zizelski</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Gatt</surname>
<given-names>Yair E.</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Margalit</surname>
<given-names>Hanah</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1045120/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Fass</surname>
<given-names>Ephraim</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Kutikov</surname>
<given-names>Ina</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Murik</surname>
<given-names>Omer</given-names>
</name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zeevi</surname>
<given-names>David A.</given-names>
</name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Savyon</surname>
<given-names>Michal</given-names>
</name>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Tau</surname>
<given-names>Luba</given-names>
</name>
<xref ref-type="aff" rid="aff5"><sup>5</sup></xref>
<xref rid="aff6" ref-type="aff"><sup>6</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Kaidar Shwartz</surname>
<given-names>Hasia</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Dveyrin</surname>
<given-names>Zeev</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1485239/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Rorman</surname>
<given-names>Efrat</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Nissan</surname>
<given-names>Israel</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x002A;</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/663948/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>National Public Health Laboratory, Public Health Directorate, Ministry of Health</institution>, <addr-line>Tel Aviv</addr-line>, <country>Israel</country></aff>
<aff id="aff2"><sup>2</sup><institution>Department of Microbiology and Molecular Genetics, Institute for Medical Research Israel-Canada, Faculty of Medicine, The Hebrew University of Jerusalem</institution>, <addr-line>Jerusalem</addr-line>, <country>Israel</country></aff>
<aff id="aff3"><sup>3</sup><institution>Translational Genomics Laboratory, Medical Genetics Institute, Shaare Zedek Medical Center</institution>, <addr-line>Jerusalem</addr-line>, <country>Israel</country></aff>
<aff id="aff4"><sup>4</sup><institution>Tel Aviv District Health Office, Ministry of Health</institution>, <addr-line>Tel Aviv</addr-line>, <country>Israel</country></aff>
<aff id="aff5"><sup>5</sup><institution>Tel Aviv Sourasky Medical Center</institution>, <addr-line>Tel Aviv</addr-line>, <country>Israel</country></aff>
<aff id="aff6"><sup>6</sup><institution>Faculty of Medicine, Tel Aviv University</institution>, <addr-line>Tel Aviv</addr-line>, <country>Israel</country></aff>
<author-notes>
<fn fn-type="edited-by" id="fn0005"><p>Edited by: Farzam Vaziri, University of California, Davis, United States</p></fn>
<fn fn-type="edited-by" id="fn0006"><p>Reviewed by: Shima Hadifar, Pasteur Institute of Iran (PII), Iran; John S. Spencer, Colorado State University, United States</p></fn>
<corresp id="c001">&#x002A;Correspondence: Israel Nissan, <email>Israel.nissan@gmail.com</email></corresp>
<fn fn-type="equal" id="fn0004"><p><sup>&#x2020;</sup>These authors have contributed equally to this work</p></fn>
</author-notes>
<pub-date pub-type="epub">
<day>06</day>
<month>11</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>10</volume>
<elocation-id>1292665</elocation-id>
<history>
<date date-type="received">
<day>11</day>
<month>09</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>16</day>
<month>10</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2023 Rubinstein, Makhon, Losev, Valenci, Gatt, Margalit, Fass, Kutikov, Murik, Zeevi, Savyon, Tau, Kaidar Shwartz, Dveyrin, Rorman and Nissan.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Rubinstein, Makhon, Losev, Valenci, Gatt, Margalit, Fass, Kutikov, Murik, Zeevi, Savyon, Tau, Kaidar Shwartz, Dveyrin, Rorman and Nissan</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>Coinfection of HIV and multidrug-resistant tuberculosis (MDR-TB) presents significant challenges in terms of the treatment and prognosis of tuberculosis, leading to complexities in managing the disease and impacting the overall outcome for TB patients. This study presents a remarkable case of a patient with MDR-TB and HIV coinfection who survived for over 8&#x2009;years, despite poor treatment adherence and comorbidities. Whole genome sequencing (WGS) of the infecting <italic>Mycobacterium tuberculosis</italic> (<italic>Mtb</italic>) strain revealed a unique genomic deletion, spanning 18 genes, including key genes involved in hypoxia response, intracellular survival, immunodominant antigens, and dormancy. This deletion, that we have called &#x201C;Del-X,&#x201D; potentially exerts a profound influence on the bacterial physiology and its virulence. Only few similar deletions were detected in other non-related <italic>Mtb</italic> genomes worldwide. <italic>In vivo</italic> evolution analysis identified drug resistance and metabolic adaptation mutations and their temporal dynamics during the patient&#x2019;s treatment course.</p>
</abstract>
<kwd-group>
<kwd><italic>Mycobacterium tuberculosis</italic> (MTB)</kwd>
<kwd>acquired immunodeficiency syndrome (AIDS)</kwd>
<kwd>human immunodeficiency virus (HIV)</kwd>
<kwd>multidrug-resistant tuberculosis (MDR-TB)</kwd>
<kwd>HspX</kwd>
<kwd>whole genome sequencing (WGS)</kwd>
<kwd><italic>in vivo</italic> evolution</kwd>
<kwd>dormancy</kwd>
</kwd-group>
<counts>
<fig-count count="5"/>
<table-count count="0"/>
<equation-count count="0"/>
<ref-count count="85"/>
<page-count count="12"/>
<word-count count="9660"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Infectious Diseases: Pathogenesis and Therapy</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="sec1">
<title>Introduction</title>
<p>Tuberculosis (TB) is a major cause of morbidity and one of the leading causes of death worldwide (<xref ref-type="bibr" rid="ref1">1</xref>). TB/HIV coinfection and drug resistant TB pose challenges to the TB elimination efforts. Effective treatment of TB, particularly multidrug-resistant TB (MDR-TB) and extensively drug resistant TB (XDR-TB) is long, expensive and toxic. The long treatment duration [6&#x2013;22&#x2009;months (<xref ref-type="bibr" rid="ref2">2</xref>)], and severe adverse effect of anti MDR-TB medications often lead to poor patient adherence to the treatment, putting not only these patients, but also the people in their close proximity, at risk (<xref ref-type="bibr" rid="ref3">3</xref>, <xref ref-type="bibr" rid="ref4">4</xref>). Discontinuous treatment and partial treatment are two leading causes for drug resistance (<xref ref-type="bibr" rid="ref5">5</xref>).</p>
<p>Coinfection of TB and HIV exacerbates the symptoms of both diseases in a synergetic manner, as both pathologies impair the immune system&#x2019;s ability to defend against these pathogens (<xref ref-type="bibr" rid="ref6">6</xref>). Coinfection with HIV significantly increases mortality among TB patients (<xref ref-type="bibr" rid="ref7">7</xref>). For example, a study conducted in Eastern Europe reported a median survival time of 5.9&#x2009;years following diagnosis with MDR-TB or XDR-TB, which decreases to 1.9&#x2009;years for patients also diagnosed with HIV (<xref ref-type="bibr" rid="ref8">8</xref>). Additional risk factors for death among TB patients include alcoholism, male gender and old age, among others (<xref ref-type="bibr" rid="ref7">7</xref>, <xref ref-type="bibr" rid="ref9">9</xref>).</p>
<p><italic>Mycobacterium tuberculosis (Mtb)</italic>, the major causative agent of TB, can exist inside the human host either in an active, replicating state, or in a dormant state. Dormancy is a programed response which enables <italic>Mtb</italic> to survive inside macrophage phagosomes under conditions of starvation, hypoxia, low pH, and oxidative and nitrosative stress. During dormancy, bacteria cease replication, downregulate central metabolism, and switch to anaerobic metabolism (<xref ref-type="bibr" rid="ref10">10</xref>). The induction of dormancy is regulated by the dormancy survival DosR regulon (<xref ref-type="bibr" rid="ref11">11</xref>, <xref ref-type="bibr" rid="ref12">12</xref>). This dormant state allows latent infection that can persist for many years or even a lifetime without clinical symptoms. However, if a critical number of bacteria manage to evade the host immune system and resume an active, replicating state, the patient develops an active, symptomatic TB disease (<xref ref-type="bibr" rid="ref13">13</xref>, <xref ref-type="bibr" rid="ref14">14</xref>).</p>
<p>The chronic nature of infection with <italic>Mtb</italic> allows for within-host evolution of the infecting bacteria over time. This process involves the accumulation of mutations, which are believed to contribute to the pathogen&#x2019;s adaptation to the environmental conditions within the host. Indirect evidence supporting the selective advantage of these mutations is the observation that they often occur in the same genes across different patients, suggesting homoplastic evolution. A major mycobacterial function influenced by these mutations is drug resistance, frequently accompanied by compensatory mutations that mitigate the fitness cost of drug resistance mutations. Genes linked to virulence, lipid and cell wall metabolism and glycerol-3-phosphate metabolism have also been demonstrated to accumulate mutations during host infection (<xref ref-type="bibr" rid="ref15">15</xref>, <xref ref-type="bibr" rid="ref16">16</xref>).</p>
<p>Compared to other bacteria, the genome of <italic>Mtb</italic> is highly conserved. Modern <italic>Mtb</italic> strains do not typically acquire new genomic features through horizontal gene transfer (<xref ref-type="bibr" rid="ref17">17</xref>, <xref ref-type="bibr" rid="ref18">18</xref>). However, <italic>Mtb</italic> strains often lose fragments of their genomes. Baena et al. (<xref ref-type="bibr" rid="ref19">19</xref>) found 4,111 genomic deletions longer than 1,000&#x2009;bp in 522 <italic>Mtb</italic> isolates of L4 lineage. One of those deletions, that is 6,479&#x2009;bp long and results in the loss of 10 genes, was observed in 249 genomes from four different countries. Other large deletions were also detected in multiple strains. Most of these strains were obtained from sputum, which is produced exclusively during active respiratory tuberculosis. Assuming each unique deletion originated from a distinct event, it is therefore evident that some large genomic deletions do not impede pathogeny and transmission.</p>
<p>Israel has maintained ongoing surveillance of all TB cases since 1997. Patients with suspected pulmonary TB in Israel are requested to provide sputum or other biological specimens for culture. All isolates are sent by primary laboratories to the National Mycobacterium Reference Laboratory (NMRL) in Tel Aviv, Israel. NMRL identifies the strain, performs drug susceptibility tests, and maintains a bank of strains as part of routine work. Since 2019, all new TB strains are fully sequenced (WGS) for epidemiological purposes and for drug susceptibility prediction. Between 2010 and 2019, NMRL received 35&#x2009;TB samples from a patient diagnosed with both TB and HIV. These samples were found to belong to a single clone of MDR-TB characterized by a unique genomic deletion spanning 18 genes. In spite of the patient&#x2019;s long progressive AIDS and active TB, incomplete treatment, and multiple risk factors, this patient survived for almost 9&#x2009;years. Here we describe this exceptional case, and discuss the distinct features of the infecting pathogen that may have contributed to this survival.</p>
</sec>
<sec sec-type="results" id="sec2">
<title>Results</title>
<sec id="sec3">
<title>Case description</title>
<p>A man in his late 50&#x2019;s was diagnosed with multidrug resistant tuberculosis (MDR-TB) in 2010, 5&#x2009;years after being diagnosed with HIV. Over the next eight and a half years, he suffered from this severe and notorious co-morbidity, and eventually passed away in 2019 at his late 60s (<xref ref-type="fig" rid="fig1">Figure 1</xref>). The recorded causes of deaths included TB, alcoholism, AIDS, and acute respiratory failure.</p>
<fig position="float" id="fig1">
<label>Figure 1</label>
<caption>
<p>Timeline medical history of Mr. X. Colored dots represent susceptibility to antibiotics of the isolates. Yellow rectangles show the duration of antimicrobial treatment with a specific antibiotic. The ticks on the right side indicate the sequenced samples.</p>
</caption>
<graphic xlink:href="fmed-10-1292665-g001.tif"/>
</fig>
<p>The patient, who will be referred to as &#x201C;Mr. X&#x201D; from this point on, was born in Russia and immigrated to Israel. Shortly after his immigration he was diagnosed with tuberculosis. His compliance with medical treatment was poor and sporadic, therefore his medical records were sparse, difficult to track and incomplete.</p>
<p>Throughout the years, Mr. X started treatment multiple times, either in community settings or during hospitalization. Importantly, he never completed any of his therapies and often disappeared from the medical authorities&#x2019; radar. The 35&#x2009;TB positive sputum samples collected in each of his therapy were sent to NMRL and tested for anti-TB drugs susceptibility (<xref ref-type="fig" rid="fig1">Figure 1</xref>; <xref rid="SM1" ref-type="supplementary-material">Supplementary Table S1</xref>). Eleven samples were subjected to Whole Genome Sequencing (WGS), which confirmed that the infection comprised a single clone of <italic>Mtb</italic> Beijing lineage throughout the entire illness duration.</p>
<p>Initially, Mr. X was treated using standard therapy for sensitive TB, i.e., a combination of rifampicin (RIF), isoniazid (INH), ethambutol (EMB), and pyrazinamide (PZA; <xref ref-type="fig" rid="fig1">Figure 1</xref>). However, after the laboratory reported that the infecting strain is resistant to RIF, INH and PZA, these drugs were replaced, and within 2&#x2009;months of the first TB diagnosis, his therapy regimens included primarily second line drugs. Overall, he was treated with 13 different anti-tuberculosis antibiotics (<xref ref-type="fig" rid="fig1">Figure 1</xref>).</p>
<p>Despite his HIV infection, Mr. X firmly refused antiretroviral therapy (ART). In 2017, his CD4 count was around 80/&#x03BC;L, corresponding to 10% of lymphocytes, with CD4/CD8 ratio lower than 0.2 and a viral load near 350,000 per mL. In the following year, his tests revealed a CD4 count around 80 per &#x03BC;L (10%), CD4/CD8 ratio lower than 0.2, and viral load around 900,000 per mL. These tests indicate an active AIDS, with a severely compromised immune system.</p>
<p>In light of his HIV and MDR-TB comorbidity, his poor compliance to treatment, his age, male gender and alcoholism, the survival of this patient around 9&#x2009;years is very unusual.</p>
</sec>
<sec id="sec4">
<title>Drug resistance mutations</title>
<p>Drug Susceptibility Testing (DST) results are shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>; <xref rid="SM1" ref-type="supplementary-material">Supplementary Table S1</xref>. Antibiotic resistance was also predicted based on whole genome sequencing (WGS), by mapping short reads of each TB isolate to the reference genome, and comparing the variations identified to a database of known resistance mutations (<xref ref-type="bibr" rid="ref20">20</xref>). The genomic resistance prediction generally agreed with the phenotypic resistance determined by DST. We identified mutations that are known to confer resistance to INH, RIF, EMB, and STM (<xref rid="SM1" ref-type="supplementary-material">Supplementary Table S2</xref>). The situation was different for PZA resistance. The gene <italic>pncA</italic> encodes an enzyme required for the transformation of PZA into its active form, pyrazinoic acid. Genetic resistance to PZA is usually caused by mutations in this gene (<xref ref-type="bibr" rid="ref20">20</xref>). However, while DST results showed PZA resistance, no mutation was found that could account for it. Upon careful inspection, we noticed there was no coverage of the genomic region around <italic>pncA</italic> by the strain&#x2019;s short reads. This uncovered region proved to be a genomic deletion of size 15.7&#x2009;kb (Del-X, see next section, <xref ref-type="fig" rid="fig2">Figure 2</xref>). The lack of functional <italic>pncA</italic> gene accounted for the PZA resistance observed by DST. The acquisition of an A90V mutation in <italic>gyrA</italic>, explains the transition to fluoroquinolones (FQs) resistance, shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>. This mutation, which is known to confer FQ resistance, was not present in any of the isolates until May 2016. It first appeared in August 2016 and was present in the 2017 isolates, but disappeared in the genome of the last, FQs sensitive, sample taken in 2019. A frameshift mutation in the gene <italic>ethA</italic>, classified by the WHO as associated with &#x201C;ethionamide resistance&#x2014;interim,&#x201D; was found in the genomes of 8 out of the 11 sequenced samples. A careful inspection revealed a significant coverage of this mutation in some reads of an additional sequenced sample. The significance of <italic>ethA</italic> to ethionamide (ETO) resistance is similar to that of <italic>pncA</italic> to PZA resistance: it alters this pro-drug to its active form. Nevertheless, this mutation only provided intermediate resistance (<xref ref-type="fig" rid="fig1">Figure 1</xref>). Two of the sequenced genomes did not have this mutation, but had another <italic>ethA</italic> mutation of uncertain significance, which was not shared by the other samples (<xref rid="SM1" ref-type="supplementary-material">Supplementary Table S2</xref>). In contrast to genomic prediction, no differences in ETO DST were detected among the different samples. No clofazimine (CFZ) resistance mutation was identified to explain the intermediate resistance observed in most isolates. We identified one mutation with uncertain significance regarding capreomycin (CAP) and amikacin (AMK) resistance in all samples. We were unable to explain the differences among AMK DST results using genomic data. Genomic-based susceptibility states were not predicted for cycloserine (CYC) and clarithromycin (CLR).</p>
<fig position="float" id="fig2">
<label>Figure 2</label>
<caption>
<p>Del-X genes and coverage. On the top: a genomic map of the Del-X region according to the Mtb H37Rv reference genome. Genes in red are immunogenic. Symbols represent gene functions of interest. On the bottom: coverage of the Del-X region with short Illumina reads of sample &#x201C;Mr.X 10&#x2013;10&#x201D;; combined coverage of all Nanopore long reads from Mr. X&#x2019;s samples except for &#x201C;Mr.X 02&#x2013;19&#x201D;; coverage of the isolate that is most closely related to Mr. X&#x2019;s samples, according to Pathogen Detection. Dashed lines represent deletion borders, and left labels on ticks show coverage depth. BWA-MEM (<xref ref-type="bibr" rid="ref21">21</xref>) was used to map short reads, and minimap2 v2.24-r1122 (<xref ref-type="bibr" rid="ref22">22</xref>) was used to map long reads.</p>
</caption>
<graphic xlink:href="fmed-10-1292665-g002.tif"/>
</fig>
</sec>
<sec id="sec5">
<title>Del-X</title>
<p>As mentioned above, WGS of the <italic>Mtb</italic> strain from Mr. X revealed a genomic deletion of size 15.7&#x2009;kb, spanning bases 2,274,648&#x2013;2,290,386 when compared to the reference <italic>Mycobacterium tuberculosis</italic> strain H37Rv (RefSeq id NC_000962.3, <xref ref-type="fig" rid="fig2">Figure 2</xref>). We have designated this deletion as Del-X. The precise boundaries of Del-X were verified using both short and long read sequencing technologies (<xref ref-type="fig" rid="fig2">Figure 2</xref>). The reference genome contains 18 genes in this region. These genes play important roles in the bacterial response to hypoxia, antibiotics susceptibility, interaction with the host, intracellular survival, and dormancy. Many of them are considered immunogenic (<xref rid="SM1" ref-type="supplementary-material">Supplementary Table S3</xref>). One of these genes is <italic>pncA</italic>, the absence of which explains the pathogen&#x2019;s resistance to PZA, as discussed previously. Another key gene in Del-X is <italic>hspX</italic>, which encodes an alpha-crystallin&#x2014;a small heat-shock protein, which is essential for long-term viability during latent, asymptomatic infections (<xref ref-type="bibr" rid="ref23">23</xref>). The deletion of these 18 genes is thus expected to affect the pathogen&#x2019;s physiology, in particular to impose a pressure against intra-phagosome existence, as further discussed in the discussion section.</p>
</sec>
<sec id="sec6">
<title>Closely related genomes</title>
<p>To find the origin of the <italic>Mtb</italic> strain that infected Mr. X, we utilized the NCBI Pathogen Detection platform (<xref ref-type="bibr" rid="ref24">24</xref>) to identify closely related sequenced isolates. NCBI Pathogen Detection Project clusters related pathogen genome sequences in order to identify transmission networks and outbreaks. Based on genomic similarity, 13 <italic>Mtb</italic> isolates closely related to the strain from Mr. X were identified (<xref ref-type="fig" rid="fig3">Figure 3</xref>; <xref rid="SM1" ref-type="supplementary-material">Supplementary Table S4</xref>). These isolates were collected between 2007 and 2015 in Russia, Sweden, Moldova, and Iran. The genome most closely related to the clone from Mr. X was collected in 2009 in Russia, differing from it by only four SNPs (<xref ref-type="fig" rid="fig3">Figure 3</xref>). None of these genomes contains Del-X, suggesting that the deletion event might have occurred in Mr. X immediately following his infection, or in another unknown host prior to his infection. Mr. X immigrated to Israel from Russia the same year in which he was diagnosed with active tuberculosis. Thus, it is reasonable to assume that he was infected in Russia, where a significant number of the isolates, which belong to the same cluster, were collected.</p>
<fig position="float" id="fig3">
<label>Figure 3</label>
<caption>
<p>Minimum spanning tree of Mr. X isolates and other closely related genomes: Minimum spanning tree of 11 Mr. X&#x2019;s samples and 13 closest samples from the Pathogen detection cluster PDS000010075.9. Numbers on branches represent the distance between nodes in SNPs, while nodes present a distance of zero SNPs. Nodes are colored according to their countries of origin. For visual purposes, the SAMN23441947-SAMEA1403736 branch has been shortened.</p>
</caption>
<graphic xlink:href="fmed-10-1292665-g003.tif"/>
</fig>
</sec>
<sec id="sec7">
<title>Other genomes with deletions in Del-X region</title>
<p>In order to estimate the uniqueness of the Del-X deletion, we searched for deletions overlapping the Del-X region in all non-redundant, well-annotated, publicly available <italic>Mtb</italic> genomes. Among the 7,028 RefSeq <italic>Mtb</italic> genomes found, 252 had deletions of various lengths overlapping Del-X (<xref rid="SM1" ref-type="supplementary-material">Supplementary File S5</xref>). Of these, 51 had a deletion longer than 1&#x2009;kb, 33 had a deletion longer than 3&#x2009;kb, and 22 had a deletion longer than 5&#x2009;kb (<xref ref-type="fig" rid="fig4">Figure 4</xref>; <xref rid="SM1" ref-type="supplementary-material">Supplementary File S5</xref>). These genomes were isolated from various parts of the world, and belonged to different <italic>Mtb</italic> lineages, primarily the Beijing and the LAM lineages. It is important to notice that not all these 22 genomes represent unique deletion events. In fact, nine of them belong to three different clusters of nearly identical genomes. One such pair is particularly interesting: two clonal genomes with a 23.2&#x2009;kb genomic deletion spanning almost the entire Del-X region, collected from a parent&#x2013;child pair in Lima, Peru (<xref ref-type="bibr" rid="ref25">25</xref>) (BioSamples SAMN10697670 and SAMN10697510, rows 4 and 5 in <xref ref-type="fig" rid="fig4">Figure 4</xref>). This observation suggests that despite its potential impact on dormancy, hypoxia response and intracellular survival, Del-X does not prevent host-to-host transmission. This pair of clonal isolates belong to a larger cluster in NCBI Pathogen Detection (<xref ref-type="bibr" rid="ref24">24</xref>), comprising 223 isolates (cluster PDT000456105.1). Several of these 223 isolates have smaller deletions in Del-X region, such as a 194&#x2009;bp deletion inside the Rv2028c gene, a deletion of 195&#x2009;bp inside the <italic>acg</italic> gene (Rv2032), and five clonal isolates with a deletion spanning seven genes inside Del-X, including <italic>pncA</italic>. However, the larger 23.2&#x2009;kb genomic deletion is not likely to have emerged from any of these smaller deletions, as this region is intact in the genome, which is closest to the pair with this large deletion (BioSample SAMEA2683069). Similarly, BioSample SAMEA2682984, a strain that is clonal (0 SNPs in NCBI Pathogen Detection) to the cluster of strains with large deletions appearing in rows 12&#x2013;16 of <xref ref-type="fig" rid="fig4">Figure 4</xref>, has no deletion in Del-X region.</p>
<fig position="float" id="fig4">
<label>Figure 4</label>
<caption>
<p>Deletions overlapping Del-X in other <italic>Mtb</italic> genomes: On the top: a genomic map of the Del-X region according to the H37Rv reference genome. Genes in red are immunogenic. Symbols represent gene functions of interest. On the bottom: Samples from the NCBI RefSeq containing a genomic deletion of at least 5 kb that overlaps the Del-X region. The deletions are shown as orange rectangles. Ticks on the right display country of origin and lineage. The left ticks indicate the BioSample accession number, and the asterisk denotes that the deletion was confirmed using the sequencing data from the SRA database. Bars to the left of BioSample IDs indicate identical deletions. Pathogen detection was used to find the number of SNP between samples with identical deletions. The turquoise bar shows the parent&#x2013;child pair with one SNP difference. The samples marked with the red bar have zero distance between them. The brown bar shows two samples with 16 SNPs difference.</p>
</caption>
<graphic xlink:href="fmed-10-1292665-g004.tif"/>
</fig>
</sec>
<sec id="sec8">
<title><italic>In vivo</italic> evolution</title>
<p>All <italic>Mtb</italic> sequenced isolates from Mr. X were clonal. Nevertheless, due to the large deletion, the lack of adherence to therapy, exposure to multiple drugs and the comorbidities, the <italic>Mtb</italic> strain infecting Mr. X was subjected to unique selective pressures during the infection. We therefore utilized the pipeline developed by Gatt and Margalit (<xref ref-type="bibr" rid="ref15">15</xref>), and identified 26 <italic>Mtb</italic> genes which underwent non-synonymous mutations during the eight and a half years of active tuberculosis (<xref ref-type="fig" rid="fig5">Figure 5</xref>). Nine of these genes have been previously identified by Gatt and Margalit as &#x201C;adaptive genes,&#x201D; that tend to undergo changes at a relatively high frequency during in-host infection of <italic>Mtb</italic> (<xref ref-type="bibr" rid="ref15">15</xref>), many of them belong to the PE-PGRS family. Among the 26 genes, seven underwent mutations predicted to lead to loss of function such as large deletions or frameshift mutations, and the remaining 19 genes underwent other non-synonymous mutations. The biological processes these genes are involved in are diverse, and include antibiotic resistance, virulence, fatty acid metabolism, central metabolism and more, reflecting the dynamic nature of bacterial adaptation to the <italic>in vivo</italic> environment.</p>
<fig position="float" id="fig5">
<label>Figure 5</label>
<caption>
<p><italic>In vivo</italic> evolution: Genes determined to undergo non-synonymous mutations during the infection of Mr. The left column includes genes determined as adaptive in <italic>Mtb</italic> by Gatt and Margalit (<xref ref-type="bibr" rid="ref15">15</xref>) and the right column includes additional genes. Functional categories of the are denoted next to the table. Genes undergoing mutations in mechanisms predicted to lead to loss of function are marked in red, and genes undergoing other non-synonymous mutations are marked in gray.</p>
</caption>
<graphic xlink:href="fmed-10-1292665-g005.tif"/>
</fig>
</sec>
</sec>
<sec sec-type="discussion" id="sec9">
<title>Discussion</title>
<p>In this study, we present a case of exceptionally long survival in a man over 60&#x2009;years of age who concurrently had progressive AIDS, active tuberculosis caused by MDR-TB and alcoholism. This remarkable extended survival is particularly notable, given the patient&#x2019;s failure to complete his anti-TB treatment regimens, and his refusal of anti-HIV treatment as reflected by his very low CD4 cells counts and high viremic state. A previous study on survival rates in patients with drug resistance TB found a median survival of 5.9&#x2009;years in MDR-TB and XDR-TB patients and 1.9&#x2009;years in patients co-infected with HIV (<xref ref-type="bibr" rid="ref8">8</xref>). The advanced age, male gender, and alcoholism of our patient makes this case even more exceptional (<xref ref-type="bibr" rid="ref8">8</xref>). The severe consequences and poor prognosis for patients co-infected with HIV/MDR-TB, particularly those that do not take ART is well documented in the literature (<xref ref-type="bibr" rid="ref26 ref27 ref28 ref29">26&#x2013;29</xref>). Furthermore, prior HIV infection, like in our case, significantly increase the chance of acquiring MDR-TB (1.42 times higher) (<xref ref-type="bibr" rid="ref30">30</xref>).</p>
<p>Over the years, since his first TB diagnosis, Mr. X started treatment several times, but never completed it. He was treated with 13 different antibiotics in total. His treatment ultimately failed to cure his TB, presumably due to the combination of his co-morbidities, the pathogen&#x2019;s multiple drug resistance and his low compliance (<xref ref-type="fig" rid="fig1">Figure 1</xref>; <xref rid="SM1" ref-type="supplementary-material">Supplementary Table S1</xref>). Initially the standard combination of RIF, INH, EMB, and PZA was administered, but when DST results from the NMRL indicated that the strain was resistant to RIF, INH, and PZA, these drugs were replaced with second line drugs. His treatment later included EMB again, despite most of his <italic>Mtb</italic> isolates exhibiting intermediate resistance to this anti-mycobacterial agent, and two isolate being fully resistant.</p>
<p>The isolation and genome sequencing of the <italic>Mtb</italic> isolates from this patient revealed an unusual genomic deletion, termed Del-X. This deletion spans from bases 2,274,648 to 2,290,386 in the H37Rv reference genome (RefSeq id NC_000962.3), causing the loss of 18 genes (Rv2028c to Rv2045c in the standard <italic>Mtb</italic> loci numbering). What effect could Del-X have on the bacterium&#x2019;s physiology, and the nature of the disease? Many of the genes inside Del-X are involved in dormancy (<xref ref-type="fig" rid="fig2">Figure 2</xref>; <xref rid="SM1" ref-type="supplementary-material">Supplementary Table S3</xref>). The <italic>hspX</italic> gene (Rv2031c, also known as <italic>acr</italic>) encodes alpha-Crystallin-like heat shock protein (Hsp16.3), which plays a pivotal role in <italic>Mtb</italic> dormancy response (<xref ref-type="bibr" rid="ref23">23</xref>, <xref ref-type="bibr" rid="ref31">31</xref>, <xref ref-type="bibr" rid="ref32">32</xref>). DosR, a transcription factor that activates the transcription of <italic>Mtb</italic> dormancy regulon, was shown to induce the expression of <italic>hspX</italic> and <italic>acg</italic> (<xref ref-type="bibr" rid="ref33">33</xref>). In turn, Rv2028 and Rv2034 positively regulate the expression of the <italic>dosR</italic> gene (<xref ref-type="bibr" rid="ref32">32</xref>, <xref ref-type="bibr" rid="ref34">34</xref>). Many of the genes within Del-X are induced by the conditions typical of macrophages phagosomes, such as hypoxia, starvation, low levels of nitric oxide, and acidic pH (<xref ref-type="bibr" rid="ref32">32</xref>, <xref ref-type="bibr" rid="ref35 ref36 ref37 ref38 ref39 ref40 ref41">35&#x2013;41</xref>). <italic>hspX</italic> actively contributes to slowing the growth of <italic>Mtb</italic> immediately following infection and is essential for its resistance to low oxygen stress (<xref ref-type="bibr" rid="ref31">31</xref>, <xref ref-type="bibr" rid="ref32">32</xref>). Therefore, it is reasonable to speculate that the strain infecting Mr. X was unable to enter a dormant state. Studies have shown that a knockout strain of <italic>Mtb</italic> lacking <italic>hspX</italic> exhibits reduced growth in macrophages (<xref ref-type="bibr" rid="ref23">23</xref>), and that depletion of HspX protein deteriorates <italic>Mtb</italic> tolerance to anaerobic conditions (<xref ref-type="bibr" rid="ref32">32</xref>). Rv2037c too was shown to enhance <italic>in vivo</italic> and <italic>ex vivo</italic> intracellular survival (<xref ref-type="bibr" rid="ref39">39</xref>). Thus, we suggest that this genomic deletion not only prevented the pathogen from entering a dormant state, but might have also abolished its ability to survive inside macrophages&#x2019; phagosomes, imposing an escape from the phagosomes into cells&#x2019; cytosol (<xref ref-type="bibr" rid="ref42">42</xref>).</p>
<p>Several of the genes in Del-X region play an important role in the host immune response to <italic>Mtb</italic>. <italic>hspX</italic>, Rv2029c, Rv2034, Rv2041c, and possibly Rv2028c encode strong immunogenic antigens (<xref ref-type="bibr" rid="ref43">43</xref>, <xref ref-type="bibr" rid="ref44">44</xref>) (<xref ref-type="fig" rid="fig2">Figure 2</xref>; <xref rid="SM1" ref-type="supplementary-material">Supplementary Table S3</xref>). HspX in particular plays a major role in activation of dendritic cells and induces the production of TNF-&#x03B1; in lymphocytes (<xref ref-type="bibr" rid="ref44">44</xref>, <xref ref-type="bibr" rid="ref45">45</xref>). <italic>acg</italic> is important for development of the typical tuberculosis granulomas (<xref ref-type="bibr" rid="ref46">46</xref>). In a study on MDR-TB patients, Rv2037c was shown to induce the production of antibodies and of pro-inflammatory cytokines (<xref ref-type="bibr" rid="ref39">39</xref>). Rv2041c increases the production of cytokines, which are important for macrophage maturation (M1 polarization) (<xref ref-type="bibr" rid="ref47">47</xref>). Therefore, Del-X is predicted to weaken the immune response against the pathogen, compared to the response against wild type <italic>Mtb</italic>. The patient&#x2019;s immune system was severely compromised by his HIV infection, as evidenced by the low CD4 counts (about 80 per &#x03BC;L), low CD4/CD8 ratio (10%), and high viral load (near 350,000 and 900,000 per mL) measured in 2017 and in 2018, respectively. The combined effect of impaired immune system and lack of <italic>Mtb</italic> determinants important for immunogenic recognition and activation, likely diminished the host ability to neutralize the infection, resulting in a constantly active TB.</p>
<p>On the other hand, Del-X may have severely affected the pathogen&#x2019;s proliferation inside the host. Two of the deleted genes, Rv2028 and Rv2040c, have been previously identified in a comprehensive screen for genes that are important for in-host survival (<xref ref-type="bibr" rid="ref48">48</xref>). <italic>acg</italic> has also been shown to be an essential virulence factor, as its deletion attenuated the bacterial growth in acute and persistent mice infections (<xref ref-type="bibr" rid="ref46">46</xref>). The Rv2038c-Rv2041c operon encodes a sugar importer of the ABC transporter family, which is thought to play an important role in virulence or in-host survival (<xref ref-type="bibr" rid="ref49">49</xref>), and its absence may limit the bacterial energy source. Furthermore, Del-X may impair lipid metabolism, due to the deletion of Rv2034 and Rv2037c (<xref ref-type="bibr" rid="ref34">34</xref>, <xref ref-type="bibr" rid="ref39">39</xref>), and nucleotide metabolism, due to the deletion of Rv2030c (<xref ref-type="bibr" rid="ref32">32</xref>). These deficiencies could further attenuate the growth of this strain. We hypothesize that Del-X genomic deletion resulted in an attenuated <italic>Mtb</italic> strain, unable to enter dormancy, yet with poor proliferation rate. The reduced virulence of this strain may have been a major factor contributing to the long survival of the patient. However, in contrast to this hypothesis, it is worth noting that some evidences have shown that deletion of <italic>hspX</italic> increases bacterial growth <italic>in vivo</italic> (<xref ref-type="bibr" rid="ref31">31</xref>). Notably, the <italic>Mtb</italic> mutant used in the paper of Hu et al. (<xref ref-type="bibr" rid="ref31">31</xref>), was deleted only in <italic>hspX</italic> gene, on the other hand, the del-X mutants in our study, contain an unusually large genomic deletion that covers 18 genes involved in multiple cellular functions, as described. Therefore, these strains cannot be compared. An obvious effect of Del-X was PZA resistance, due to the deletion of <italic>pncA</italic>. Treatment with PZA likely selected for the strain with this deletion, allowing it to take over the population, despite its defects. Indeed, although deletions in this genomic region are rare, almost all known <italic>Mtb</italic> genomes with a deletion of 5&#x2009;kb or more in this region also exhibit a deletion of <italic>pncA</italic> (<xref ref-type="fig" rid="fig4">Figure 4</xref>; <xref rid="SM1" ref-type="supplementary-material">Supplementary File S5</xref>). In line with this, the work of Godfroid et al. (<xref ref-type="bibr" rid="ref50">50</xref>), focusing on analysis of multidrug resistant MTB outbreak isolates, discovered that genomic insertions and deletions are significantly enriched in genes conferring antibiotic resistance.</p>
<p>The discovery of Del-X underscores the importance of detecting large genomic deletions in routine analysis of <italic>Mtb</italic> whole genome sequencing. Standard pipelines based on mapping to reference genomes and variant call may miss such deletions. Indeed, following this case, our laboratory has incorporated a step of identifying deletions based on uncovered genomic regions to its routine analysis.</p>
<p>We found a few <italic>Mtb</italic> genomes with genomic deletions in Del-X, from different parts of the world, mainly belonging to the Beijing or LAM lineages. As mentioned, the isolates from Mr. X also belong to the Beijing lineage that is known to be hyper-mutagenic, virulent, that rapidly acquires mutations for resistance to anti-tuberculosis drugs (<xref ref-type="bibr" rid="ref51 ref52 ref53 ref54">51&#x2013;54</xref>).</p>
<p>Interestingly, two of these genomes were isolated from a parent and their child, and are genetically clonal. Hence, one of these family members probably infected the other (<xref ref-type="bibr" rid="ref25">25</xref>). This finding suggests that the Del-X genomic deletion does not prevent transmission from one patient to another and may have already been present in the strain when Mr. X was infected. However, none of this strain&#x2019;s known closely related isolates had a deletion in the Del-X region. Thus, we cannot exclude the possibility that the deletion arose within this patient (<xref ref-type="fig" rid="fig3">Figure 3</xref>; <xref rid="SM1" ref-type="supplementary-material">Supplementary Table S4</xref>).</p>
<p>The isolates obtained from Mr. X at different time points during his nearly 9&#x2009;years of infection belong to the same clone, yet there are several variations differing between them. These variations may originate from co-existing sub populations randomly sampled, or reflect changes that accumulated over time during the infection. Our phylogenetic analysis also attempts to separate these cases and include only mutations likely to have emerged during the infection. The later possibility is generally the case in <italic>Mtb</italic> infections (<xref ref-type="bibr" rid="ref17">17</xref>). Many of the genes mutated during the infection were previously identified as &#x201C;hot spots&#x201D; for <italic>in vivo</italic> evolution during TB infection (<xref ref-type="fig" rid="fig5">Figure 5</xref>) (<xref ref-type="bibr" rid="ref15">15</xref>, <xref ref-type="bibr" rid="ref16">16</xref>). For example, the <italic>gyrA</italic> A90V mutation discussed above, which explains the emergence of FQ resistance in isolates obtained between August 2016 and November 2017, and was probably selected for by the treatment with MFX and LFX (<xref ref-type="bibr" rid="ref55">55</xref>). Furthermore, the fact that the patient did not receive continuous treatment is known to be a major risk factor for the appearance of resistance (<xref ref-type="bibr" rid="ref5">5</xref>). After the exposure to these drugs was removed, an isolate taken in February 2019 lacked this mutation, probably representing a sub-population that evaded FQs and became dominant once they were no longer used. Another interesting variation is the insertion of C in the homopolymeric tract of 7 cytosine in the <italic>glpK</italic> gene, relative to the TB reference H37Rv genome. This mutation is common to all isolates except for the one from August 2016 and the last isolate from February 2019. The <italic>glpK</italic> gene encodes a glycerol kinase involved in glycerol metabolism. This mutation is known to serve as a mechanism for drug tolerance and growth slow-down, which occurs and reverses at high frequency (<xref ref-type="bibr" rid="ref56">56</xref>). Notably, among the 26 genes mutated in-host there were six PE-PGRS genes, the functions of which are largely unknown (<xref ref-type="fig" rid="fig5">Figure 5</xref>).</p>
<p>The medical records of Mr. X indicate that he started TB treatment at least six different times, probably reflecting worsening of the symptoms (<xref ref-type="fig" rid="fig1">Figure 1</xref>). In all of these episodes, the patient stopped the treatment before completing it, at times by simply disappearing. This demonstrates the challenges faced by the health authorities in treating TB, which often involve social and personal difficulties. The case of Mr. X underscores the importance of patient adherence to appropriate medical treatment and highlights the significant challenges in achieving this goal due to personal and social issues. It is probable that with a combined therapy, including ART and suitable anti-TB drugs, this patient could have successfully overcome the infection with the del-X mutant. Treatment of low compliant and unmotivated patients is often unsuccessful and requires a tailored plan based on establishing a long-term trust relationship. Directly observed therapy (DOT) is an important part of such a plan (<xref ref-type="bibr" rid="ref57">57</xref>). Mr. X&#x2019;s case raises complex ethical, legal, and personal dilemmas. Israeli law grants district doctors the authority to forcibly hospitalize individuals suffering from diseases that jeopardize public health, provided they have obtained judicial approval. The issue of forced medical treatment is even more intricate, and, on rare occasions, the law allows for the treatment of such patients. However, it is important to note that the successful treatment of drug-resistant tuberculosis hinges on the patient&#x2019;s cooperation, which, unfortunately, was unattainable in this particular case.</p>
<p>Taken together, we hypothesize that the unusual prolonged TB illness of Mr. X was a result of a unique intricate interplay between the host and the pathogen. The potentially low immunogenic signature of the pathogen together with the severe immune deficiency of the host prevented the elimination of the pathogen. The Del-X genomic deletion prevented latency, resulting in constant active TB. At the same time, Del-X probably resulted in an attenuated pathogen, allowing the long-term survival of the patient.</p>
<p>The importance of this report stems from the unusual long survival of Mr. X, considering the circumstances, as well as the uniqueness of the genomic deletion described here. This case raises important questions about the influence of <italic>Mtb</italic> genetic factors on disease progression and the potential implications for prognosis and treatment. It is interesting to find out whether in additional cases of patients who survive active tuberculosis for many years without treatment, mutated <italic>Mtb</italic> with genomic mutations in important genes are involved. Further investigations may elucidate the precise mechanisms by which the Del-X deletion affects the pathogenicity and virulence of <italic>Mtb</italic>, as well as the host immune response to the infection. In addition, this special isolate, can be used as a model bacterium for a better understanding of host-pathogen fundamental interactions, such as <italic>Mtb</italic> escape from the phagosome to the cytosol, dormancy and more.</p>
</sec>
<sec sec-type="materials|methods" id="sec10">
<title>Materials and methods</title>
<sec id="sec12">
<title>Data classification and identification of <italic>Mycobacterium tuberculosis complex</italic></title>
<p>All cultures processed in Israel are sent to the National Mycobacterium Reference Laboratory (NMRL), in Tel Aviv, Israel. At the NMRL, smears from samples and TB cultures are stained using Ziehl-Neelsen staining and cultured on Lowenstein-Jensen (LJ) media using standard methods (<xref ref-type="bibr" rid="ref58">58</xref>), and in addition cultured on BACTECe MGITe 960 (BD, Sparks, MD, United States) system. The species are identified using conventional biochemical methods (<xref ref-type="bibr" rid="ref59">59</xref>, <xref ref-type="bibr" rid="ref60">60</xref>) and a commercially available strip DNA probe test (Hain Lifesciences, Nehren, Germany).</p>
</sec>
<sec id="sec13">
<title>Drug susceptibility testing for MTBC strains</title>
<p>Drug susceptibility testing for all the drugs was done using resistance ratio method (<xref ref-type="bibr" rid="ref61">61</xref>). Pyrazinamide (PZA) was tested using Mark&#x2019;s stepped pH method (<xref ref-type="bibr" rid="ref62">62</xref>).</p>
<p>Drugs threshold levels used to determine resistance are described in <xref rid="SM1" ref-type="supplementary-material">Supplementary Table S6</xref>. Importantly, for the purpose of the article, the categories Borderline and RR4 were redefined as Intermediate.</p>
</sec>
<sec id="sec14">
<title>Whole genome sequencing and bioinformatics</title>
<p>Genomic DNA from inactivated bacterial cultures was extracted according to the manufacturer&#x2019;s protocol using Maxwell RSC cultured cells DNA kit and Maxwell&#x00AE; 16 System (Promega) (<xref ref-type="bibr" rid="ref63">63</xref>). DNA Paired-end libraries were prepared using the Illumina Nextera XT DNA Library Preparation Kit according to Illumina protocols. For sequencing, we utilized the Illumina MiSeq platform using a MiSeq Reagent Kit v2 (500-cycles; catalog MS-102-2003) or a MiSeq Reagent Kit v3 (600-cycle), (catalog MS-102-3003).</p>
<p>Raw reads fastq files were quality analyzed by FastQC (<xref ref-type="bibr" rid="ref64">64</xref>), Kraken2 taxonomic classification (<xref ref-type="bibr" rid="ref65">65</xref>) was used to identify mixed cultures. MTBseq, a comprehensive pipeline for whole genome sequence analysis of MTBC isolates, was used as the main tool for lineage annotation, identification of genomic variations, and clusters analysis (<xref ref-type="bibr" rid="ref66">66</xref>). Lineages were identified as part of the MTBseq pipeline, based on the SNP scheme developed by Coll et al. (<xref ref-type="bibr" rid="ref67">67</xref>). Genomic variations were compared to the WHO 2021 Catalog of drug resistance associated mutations in MTBC and used to identify resistant strains to antibiotics (<xref ref-type="bibr" rid="ref20">20</xref>).</p>
</sec>
<sec id="sec15">
<title>Searching for deletions among assemblies from the NCBI</title>
<p>A search for deletions was carried out on 7,298 <italic>Mtb</italic> (taxid1773) assemblies available in the NCBI RefSeq database (<xref ref-type="bibr" rid="ref68">68</xref>) (accessed March 15 2023). After excluding replaced and anomalous assemblies, 7,028 valid assemblies remained. Assembly metadata were retrieved using the Entrez Direct,<xref ref-type="fn" rid="fn0001"><sup>1</sup></xref> and sequences were downloaded using the ncbi-genome-download v0.3.1 tool.<xref ref-type="fn" rid="fn0002"><sup>2</sup></xref> All 7,028 sequences were then aligned to the reference sequence of <italic>Mtb</italic> H37Rv (NC_000962.3) using the NUCmer program from the MUMmer v4.0.0rc1 package (<xref ref-type="bibr" rid="ref69">69</xref>) with the default options. Assembly containing uncovered regions of any length that overlap with H37Rv region 2,274,648&#x2013;2,290,368 was considered to have a deletion (<xref rid="SM1" ref-type="supplementary-material">Supplementary File S5</xref>).</p>
</sec>
<sec id="sec16">
<title>Deletions verification with sequence read archive data</title>
<p>Samples with assembly-predicted deletions longer than 1&#x2009;kb were additionally evaluated using sequencing data from the SRA database (<xref ref-type="bibr" rid="ref70">70</xref>). 32 samples out of 51 with total deletion sizes of more than 1&#x2009;kb have sequencing data in the SRA. The GRIDSS v2.13.2 software for the detection of genomic rearrangements (<xref ref-type="bibr" rid="ref71">71</xref>, <xref ref-type="bibr" rid="ref72">72</xref>) was used to examine these samples. Additionally, Samtools v1.16.1 (<xref ref-type="bibr" rid="ref73">73</xref>) was used to calculate per position read depth for all samples. The deletion is considered to be correctly verified if: (1) The deletion predicted by GRIDSS overlaps with the deletion predicted from the assembly by more than 99% of the length and vice versa; (2) More than 90% of the bases in the given region have a read depth of less than 5% of the mean genome depth. Assembly GCF_000679915.1 had four SRA runs, three of them passed, therefore we considered this assembly as verified. After applying these filters, seven samples were discarded. Samples that lack sequencing data in the SRA database were used in the further analysis without verification (<xref rid="SM1" ref-type="supplementary-material">Supplementary File S5</xref>).</p>
</sec>
<sec id="sec17">
<title>Long reads</title>
<p>For Oxford Nanopore sequencing, DNA from Mr. X 10 isolates (except 02&#x2013;19) were converted into a sequencing library using Rapid Barcoding Kit 96 (SQK-RBK110.96) chemistry according to the manufacturer&#x2019;s protocol. Sequencing was performed on a MinION sequencer with an R9.4.1 (FLO-MIN106D) flow cell, and basecalling was carried out by Guppy version 6.0.1. Sequences of all 11 isolates were used to create a single assembly.</p>
</sec>
<sec id="sec18">
<title>MST construction</title>
<p>13 samples from the Pathogen detection cluster PDS000010075.9 and 11 Mr. X&#x2019;s isolates were used to construct the minimum spanning tree. There are three samples from Pathogen detection (BioSample accession numbers SAMN23441947, SAMN15963001, and SAMN12325235) for which there are no sequencing data in the SRA database (<xref ref-type="bibr" rid="ref70">70</xref>). Assembled sequences of these samples were artificially converted to paired-end reads using a custom python script. The forward and reverse read lengths were both set to 150&#x2009;bp with a minimal coverage of 50. Prefetch and fasterq-dump from the SRA Toolkit v3.0.1<xref ref-type="fn" rid="fn0003"><sup>3</sup></xref> were used to obtain FASTQ files from the remaining 10 samples. BWA-MEM v0.7.17-r1188 (<xref ref-type="bibr" rid="ref21">21</xref>) was used to map the reads to the <italic>Mtb</italic> H37Rv (NC_000962.3) reference sequence. Viewing, sorting, and duplication removal in BAM files were done using Samtools v1.16.1 (<xref ref-type="bibr" rid="ref73">73</xref>). FreeBayes v1.3.6 (<xref ref-type="bibr" rid="ref74">74</xref>) was used to call SNPs with the &#x2212;p 1&#x2212;q three parameters. SNPs were filtered using the following criteria: Site quality&#x2009;&#x2265;&#x2009;30, variant read depth&#x2009;&#x2265;&#x2009;5, distance to nearest SNP&#x2009;&#x2265;&#x2009;3&#x2009;bp, &#x2265;5% read mapping up-or downstream, &#x2265;5% read mapping to least-covered strand, variant allele frequency&#x2009;&#x2265;&#x2009;90%. In addition, SNPs in genes associated with antibiotic resistance as well as those located in PE/PPE were excluded in the subsequent analysis (<xref ref-type="bibr" rid="ref75">75</xref>). Sites with more than 5% of missing variants were filtered out as the final stage of filtering. SNP filtering was carried out using VCFtools v 0.1.16 (<xref ref-type="bibr" rid="ref76">76</xref>), bcftools v1.16 (<xref ref-type="bibr" rid="ref73">73</xref>), vcffilter from vcflib v 1.0.3 (<xref ref-type="bibr" rid="ref77">77</xref>), and bedtools v2.30.0 (<xref ref-type="bibr" rid="ref78">78</xref>). The GrapeTree v 1.5.0 (<xref ref-type="bibr" rid="ref79">79</xref>) with MSTree method was used to create the minimum spanning tree.</p>
</sec>
<sec id="sec19">
<title>Lineage assignment</title>
<p>Genome assemblies were mapped to the reference <italic>Mtb</italic> H37Rv genome (RefSeq id NC_000962.3), and variants were detected as previously described (<xref ref-type="bibr" rid="ref80">80</xref>). These variants were compared to a SNP barcode of <italic>Mtb</italic> lineages to assign a lineage (<xref ref-type="bibr" rid="ref67">67</xref>).</p>
</sec>
<sec id="sec20">
<title>Determination of within-host evolution</title>
<p>Detection of mutations was performed as in Gatt and Margalit (<xref ref-type="bibr" rid="ref15">15</xref>). Briefly, assemblies were constructed for all paired-end libraries using kSNP v.3.1 (<xref ref-type="bibr" rid="ref81">81</xref>, <xref ref-type="bibr" rid="ref82">82</xref>), and phylogenetic structure was determined using the TRACE algorithm with high confidence progenitor-progeny isolate pairs ascertained from the different isolates of Mr. X. For each progenitor-progeny pair, <italic>breseq</italic> v.0.32 (<xref ref-type="bibr" rid="ref83">83</xref>, <xref ref-type="bibr" rid="ref84">84</xref>) was then used to determine differential variations from the reference genome <italic>Mtb</italic> H37Rv (RefSeq id NC_000962.3) appearing in the progeny isolate and not in the progenitor isolate. SnpEff v.4.1 (<xref ref-type="bibr" rid="ref85">85</xref>) was used to predict the effect of the different mutations, with high impact mutations predicted to lead to loss of function.</p>
</sec>
</sec>
<sec sec-type="data-availability" id="sec21">
<title>Data availability statement</title>
<p>The raw genomic sequences from this work have been deposited to NCBI and are available in BioProject PRJNA834625.</p>
</sec>
<sec sec-type="ethics-statement" id="sec22">
<title>Ethics statement</title>
<p>This study was approved by the National Helsinki Committee for Medical Experiments on Humans, of the Israeli Ministry of Health (# MOH-081-2021). The studies were conducted in accordance with the local legislation and institutional requirements. The human samples used in this study were acquired from a by-product of routine care or industry. Written informed consent for participation was not required from the participant or the participant&#x2019;s legal guardians/next of kin in accordance with the national legislation and institutional requirements. Written informed consent was not obtained from the individual for the publication of any potentially identifiable images or data included in this article because the patient this article is about died several years before publication. Therefore, a written informed consent was obtained from the participant&#x2019;s next of kin, for the publication of this case report, and any potentially identifiable images or data included in this article.</p>
</sec>
<sec sec-type="author-contributions" id="sec23">
<title>Author contributions</title>
<p>MR: Conceptualization, Data curation, Formal analysis, Writing &#x2013; original draft, Writing &#x2013; review &#x0026; editing. AM: Data curation, Formal analysis, Visualization, Writing &#x2013; original draft, Writing &#x2013; review &#x0026; editing. YL: Formal analysis, Writing &#x2013; original draft, Writing &#x2013; review &#x0026; editing. GV: Data curation, Formal analysis, Writing &#x2013; review &#x0026; editing. YG: Data curation, Formal analysis, Writing &#x2013; review &#x0026; editing. HM: Writing &#x2013; review &#x0026; editing. EF: Formal analysis, Writing &#x2013; review &#x0026; editing. IK: Formal analysis, Writing &#x2013; review &#x0026; editing. OM: Formal analysis, Writing &#x2013; review &#x0026; editing. DZ: Formal analysis, Writing &#x2013; review &#x0026; editing. MS: Investigation, Writing &#x2013; review &#x0026; editing. LT: Investigation, Writing &#x2013; review &#x0026; editing. HK: Writing &#x2013; review &#x0026; editing. ZD: Writing &#x2013; review &#x0026; editing. ER: Project administration, Supervision, Writing &#x2013; review &#x0026; editing. IN: Conceptualization, Data curation, Formal analysis, Methodology, Writing &#x2013; original draft, Writing &#x2013; review &#x0026; editing.</p>
</sec>
</body>
<back>
<ack>
<p>We would like to thank Adv. Nomi Bibi, Tel Aviv Health District, for the legal consultation.</p>
</ack>
<sec sec-type="funding-information" id="sec24">
<title>Funding</title>
<p>This work was supported by Israeli Ministry of Health.</p>
</sec>
<sec sec-type="COI-statement" id="sec25">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="sec100" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec sec-type="supplementary-material" id="sec26">
<title>Supplementary material</title>
<p>The Supplementary material for this article can be found online at: <ext-link xlink:href="https://www.frontiersin.org/articles/10.3389/fmed.2023.1292665/full#supplementary-material" ext-link-type="uri">https://www.frontiersin.org/articles/10.3389/fmed.2023.1292665/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Table_1.XLSX" id="SM1" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Table_2.XLSX" id="SM2" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Table_3.XLSX" id="SM3" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Table_4.XLSX" id="SM4" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Table_5.XLSX" id="SM5" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Table_6.XLSX" id="SM6" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<fn-group>
<fn id="fn0001"><p><sup>1</sup><ext-link xlink:href="https://www.ncbi.nlm.nih.gov/books/NBK179288/" ext-link-type="uri">https://www.ncbi.nlm.nih.gov/books/NBK179288/</ext-link></p></fn>
<fn id="fn0002"><p><sup>2</sup><ext-link xlink:href="https://github.com/kblin/ncbi-genome-download" ext-link-type="uri">https://github.com/kblin/ncbi-genome-download</ext-link></p></fn>
<fn id="fn0003"><p><sup>3</sup><ext-link xlink:href="https://github.com/ncbi/sra-tools/" ext-link-type="uri">https://github.com/ncbi/sra-tools/</ext-link></p></fn>
</fn-group>
<ref-list>
<title>References</title>
<ref id="ref1"><label>1.</label><citation citation-type="other"><person-group person-group-type="author"><collab id="coll1">World Health Organization</collab></person-group> (<year>2022</year>). Global tuberculosis report 2022. Licence: CC BY-NC-SA 3.0 IGO. 2022. Geneva: World Health Organization.</citation></ref>
<ref id="ref2"><label>2.</label><citation citation-type="book"><person-group person-group-type="author"><collab id="coll2">World Health Organization</collab></person-group>. <source>WHO Operational Handbook on Tuberculosis. Module 4: Treatment&#x2014;Drug-Resistant Tuberculosis Treatment, 2022 update</source>. <publisher-loc>Geneva</publisher-loc>: <publisher-name>World Health Organization</publisher-name> (<year>2022</year>).</citation></ref>
<ref id="ref3"><label>3.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sotgiu</surname><given-names>G</given-names></name> <name><surname>Sulis</surname><given-names>G</given-names></name> <name><surname>Matteelli</surname><given-names>A</given-names></name></person-group>. <article-title>Tuberculosis-a World Health Organization perspective</article-title>. <source>Microbiol Spectr</source>. (<year>2017</year>) <volume>5</volume>. doi: <pub-id pub-id-type="doi">10.1128/microbiolspec.TNMI7-0036-2016</pub-id>, PMID: <pub-id pub-id-type="pmid">28185618</pub-id></citation></ref>
<ref id="ref4"><label>4.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Linh</surname><given-names>NN</given-names></name> <name><surname>Viney</surname><given-names>K</given-names></name> <name><surname>Gegia</surname><given-names>M</given-names></name> <name><surname>Falzon</surname><given-names>D</given-names></name> <name><surname>Glaziou</surname><given-names>P</given-names></name> <name><surname>Floyd</surname><given-names>K</given-names></name> <etal/></person-group>. <article-title>World Health Organization treatment outcome definitions for tuberculosis: 2021 update</article-title>. <source>Eur Respir J</source>. (<year>2021</year>) <volume>58</volume>:<fpage>2100804</fpage>. doi: <pub-id pub-id-type="doi">10.1183/13993003.00804-2021</pub-id>, PMID: <pub-id pub-id-type="pmid">34413124</pub-id></citation></ref>
<ref id="ref5"><label>5.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gillespie</surname><given-names>SH</given-names></name></person-group>. <article-title>Evolution of drug resistance in <italic>Mycobacterium tuberculosis</italic>: clinical and molecular perspective</article-title>. <source>Antimicrob Agents Chemother</source>. (<year>2002</year>) <volume>46</volume>:<fpage>267</fpage>&#x2013;<lpage>74</lpage>. doi: <pub-id pub-id-type="doi">10.1128/AAC.46.2.267-274.2002</pub-id>, PMID: <pub-id pub-id-type="pmid">11796329</pub-id></citation></ref>
<ref id="ref6"><label>6.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bell</surname><given-names>LCK</given-names></name> <name><surname>Noursadeghi</surname><given-names>M</given-names></name></person-group>. <article-title>Pathogenesis of HIV-1 and <italic>Mycobacterium tuberculosis</italic> co-infection</article-title>. <source>Nat Rev Microbiol</source>. (<year>2018</year>) <volume>16</volume>:<fpage>80</fpage>&#x2013;<lpage>90</lpage>. doi: <pub-id pub-id-type="doi">10.1038/nrmicro.2017.128</pub-id></citation></ref>
<ref id="ref7"><label>7.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nordholm</surname><given-names>AC</given-names></name> <name><surname>Andersen</surname><given-names>AB</given-names></name> <name><surname>Wejse</surname><given-names>C</given-names></name> <name><surname>Norman</surname><given-names>A</given-names></name> <name><surname>Ekstrom</surname><given-names>CT</given-names></name> <name><surname>Andersen</surname><given-names>PH</given-names></name> <etal/></person-group>. <article-title>Mortality, risk factors, and causes of death among people with tuberculosis in Denmark, 1990-2018</article-title>. <source>Int J Infect Dis</source>. (<year>2023</year>) <volume>130</volume>:<fpage>76</fpage>&#x2013;<lpage>82</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.ijid.2023.02.024</pub-id>, PMID: <pub-id pub-id-type="pmid">36907546</pub-id></citation></ref>
<ref id="ref8"><label>8.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Balabanova</surname><given-names>Y</given-names></name> <name><surname>Ignatyeva</surname><given-names>O</given-names></name> <name><surname>Fiebig</surname><given-names>L</given-names></name> <name><surname>Riekstina</surname><given-names>V</given-names></name> <name><surname>Danilovits</surname><given-names>M</given-names></name> <name><surname>Jaama</surname><given-names>K</given-names></name> <etal/></person-group>. <article-title>Survival of patients with multidrug-resistant TB in Eastern Europe: what makes a difference?</article-title> <source>Thorax</source>. (<year>2016</year>) <volume>71</volume>:<fpage>854</fpage>&#x2013;<lpage>61</lpage>. doi: <pub-id pub-id-type="doi">10.1136/thoraxjnl-2015-207638</pub-id>, PMID: <pub-id pub-id-type="pmid">27012887</pub-id></citation></ref>
<ref id="ref9"><label>9.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wigger</surname><given-names>GW</given-names></name> <name><surname>Bouton</surname><given-names>TC</given-names></name> <name><surname>Jacobson</surname><given-names>KR</given-names></name> <name><surname>Auld</surname><given-names>SC</given-names></name> <name><surname>Yeligar</surname><given-names>SM</given-names></name> <name><surname>Staitieh</surname><given-names>BS</given-names></name></person-group>. <article-title>The impact of alcohol use disorder on tuberculosis: a review of the epidemiology and potential immunologic mechanisms</article-title>. <source>Front Immunol</source>. (<year>2022</year>) <volume>13</volume>:<fpage>864817</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fimmu.2022.864817</pub-id></citation></ref>
<ref id="ref10"><label>10.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Singh</surname><given-names>PR</given-names></name> <name><surname>Vijjamarri</surname><given-names>AK</given-names></name> <name><surname>Sarkar</surname><given-names>D</given-names></name></person-group>. <article-title>Metabolic switching of <italic>Mycobacterium tuberculosis</italic> during hypoxia is controlled by the virulence regulator PhoP</article-title>. <source>J Bacteriol</source>. (<year>2020</year>) <volume>202</volume>. doi: <pub-id pub-id-type="doi">10.1128/JB.00705-19</pub-id>, PMID: <pub-id pub-id-type="pmid">31932312</pub-id></citation></ref>
<ref id="ref11"><label>11.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Leistikow</surname><given-names>RL</given-names></name> <name><surname>Morton</surname><given-names>RA</given-names></name> <name><surname>Bartek</surname><given-names>IL</given-names></name> <name><surname>Frimpong</surname><given-names>I</given-names></name> <name><surname>Wagner</surname><given-names>K</given-names></name> <name><surname>Voskuil</surname><given-names>MI</given-names></name></person-group>. <article-title>The <italic>Mycobacterium tuberculosis</italic> DosR regulon assists in metabolic homeostasis and enables rapid recovery from nonrespiring dormancy</article-title>. <source>J Bacteriol</source>. (<year>2010</year>) <volume>192</volume>:<fpage>1662</fpage>&#x2013;<lpage>70</lpage>. doi: <pub-id pub-id-type="doi">10.1128/JB.00926-09</pub-id></citation></ref>
<ref id="ref12"><label>12.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Peddireddy</surname><given-names>V</given-names></name> <name><surname>Doddam</surname><given-names>SN</given-names></name> <name><surname>Ahmed</surname><given-names>N</given-names></name></person-group>. <article-title>Mycobacterial dormancy systems and host responses in tuberculosis</article-title>. <source>Front Immunol</source>. (<year>2017</year>) <volume>8</volume>:<fpage>84</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fimmu.2017.00084</pub-id>, PMID: <pub-id pub-id-type="pmid">28261197</pub-id></citation></ref>
<ref id="ref13"><label>13.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Delogu</surname><given-names>G</given-names></name> <name><surname>Sali</surname><given-names>M</given-names></name> <name><surname>Fadda</surname><given-names>G</given-names></name></person-group>. <article-title>The biology of <italic>mycobacterium tuberculosis</italic> infection</article-title>. <source>Mediterr J Hematol Infect Dis</source>. (<year>2013</year>) <volume>5</volume>:<fpage>e2013070</fpage>. doi: <pub-id pub-id-type="doi">10.4084/MJHID.2013.070</pub-id>, PMID: <pub-id pub-id-type="pmid">24363885</pub-id></citation></ref>
<ref id="ref14"><label>14.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Verma</surname><given-names>A</given-names></name> <name><surname>Ghoshal</surname><given-names>A</given-names></name> <name><surname>Dwivedi</surname><given-names>VP</given-names></name> <name><surname>Bhaskar</surname><given-names>A</given-names></name></person-group>. <article-title>Tuberculosis: the success tale of less explored dormant <italic>Mycobacterium tuberculosis</italic></article-title>. <source>Front Cell Infect Microbiol</source>. (<year>2022</year>) <volume>12</volume>:<fpage>1079569</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fcimb.2022.1079569</pub-id>, PMID: <pub-id pub-id-type="pmid">36619761</pub-id></citation></ref>
<ref id="ref15"><label>15.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gatt</surname><given-names>YE</given-names></name> <name><surname>Margalit</surname><given-names>H</given-names></name></person-group>. <article-title>Common adaptive strategies underlie within-host evolution of bacterial pathogens</article-title>. <source>Mol Biol Evol</source>. (<year>2021</year>) <volume>38</volume>:<fpage>1101</fpage>&#x2013;<lpage>21</lpage>. doi: <pub-id pub-id-type="doi">10.1093/molbev/msaa278</pub-id></citation></ref>
<ref id="ref16"><label>16.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ley</surname><given-names>SD</given-names></name> <name><surname>de Vos</surname><given-names>M</given-names></name> <name><surname>Van Rie</surname><given-names>A</given-names></name> <name><surname>Warren</surname><given-names>RM</given-names></name></person-group>. <article-title>Deciphering within-host microevolution of <italic>Mycobacterium tuberculosis</italic> through whole-genome sequencing: the phenotypic impact and way forward</article-title>. <source>Microbiol Mol Biol Rev</source>. (<year>2019</year>) <volume>83</volume>. doi: <pub-id pub-id-type="doi">10.1128/MMBR.00062-18</pub-id>, PMID: <pub-id pub-id-type="pmid">30918049</pub-id></citation></ref>
<ref id="ref17"><label>17.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gagneux</surname><given-names>S</given-names></name></person-group>. <article-title>Ecology and evolution of <italic>Mycobacterium tuberculosis</italic></article-title>. <source>Nat Rev Microbiol</source>. (<year>2018</year>) <volume>16</volume>:<fpage>202</fpage>&#x2013;<lpage>13</lpage>. doi: <pub-id pub-id-type="doi">10.1038/nrmicro.2018.8</pub-id>, PMID: <pub-id pub-id-type="pmid">29456241</pub-id></citation></ref>
<ref id="ref18"><label>18.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dos Vultos</surname><given-names>T</given-names></name> <name><surname>Mestre</surname><given-names>O</given-names></name> <name><surname>Rauzier</surname><given-names>J</given-names></name> <name><surname>Golec</surname><given-names>M</given-names></name> <name><surname>Rastogi</surname><given-names>N</given-names></name> <name><surname>Rasolofo</surname><given-names>V</given-names></name> <etal/></person-group>. <article-title>Evolution and diversity of clonal bacteria: the paradigm of <italic>Mycobacterium tuberculosis</italic></article-title>. <source>PLoS One</source>. (<year>2008</year>) <volume>3</volume>:<fpage>e1538</fpage>. doi: <pub-id pub-id-type="doi">10.1371/journal.pone.0001538</pub-id></citation></ref>
<ref id="ref19"><label>19.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Baena</surname><given-names>A</given-names></name> <name><surname>Cabarcas</surname><given-names>F</given-names></name> <name><surname>Ocampo</surname><given-names>JC</given-names></name> <name><surname>Barrera</surname><given-names>LF</given-names></name> <name><surname>Alzate</surname><given-names>JF</given-names></name></person-group>. <article-title>Large genomic deletions delineate <italic>Mycobacterium tuberculosis</italic> L4 sublineages in south American countries</article-title>. <source>PLoS One</source>. (<year>2023</year>) <volume>18</volume>:<fpage>e0285417</fpage>. doi: <pub-id pub-id-type="doi">10.1371/journal.pone.0285417</pub-id>, PMID: <pub-id pub-id-type="pmid">37205685</pub-id></citation></ref>
<ref id="ref20"><label>20.</label><citation citation-type="other"><person-group person-group-type="author"><collab id="coll3">World Health Organization</collab></person-group> (<year>2013</year>). Catalogue of mutations in <italic>Mycobacterium tuberculosis</italic> complex and their association with drug resistance. Licence: CC BY-NC-SA 3.0 IGO. Geneva: World Health Organization.</citation></ref>
<ref id="ref21"><label>21.</label><citation citation-type="other"><person-group person-group-type="author"><name><surname>Li</surname><given-names>H.</given-names></name></person-group> (<year>2013</year>) Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv [Preprint]. doi: <pub-id pub-id-type="doi">10.48550/arXiv.1303.3997</pub-id></citation></ref>
<ref id="ref22"><label>22.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Li</surname><given-names>H</given-names></name></person-group>. <article-title>Minimap2: pairwise alignment for nucleotide sequences</article-title>. <source>Bioinformatics</source>. (<year>2018</year>) <volume>34</volume>:<fpage>3094</fpage>&#x2013;<lpage>100</lpage>. doi: <pub-id pub-id-type="doi">10.1093/bioinformatics/bty191</pub-id>, PMID: <pub-id pub-id-type="pmid">29750242</pub-id></citation></ref>
<ref id="ref23"><label>23.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yuan</surname><given-names>Y</given-names></name> <name><surname>Crane</surname><given-names>DD</given-names></name> <name><surname>Simpson</surname><given-names>RM</given-names></name> <name><surname>Zhu</surname><given-names>YQ</given-names></name> <name><surname>Hickey</surname><given-names>MJ</given-names></name> <name><surname>Sherman</surname><given-names>DR</given-names></name> <etal/></person-group>. <article-title>The 16-kDa alpha-crystallin (Acr) protein of <italic>Mycobacterium tuberculosis</italic> is required for growth in macrophages</article-title>. <source>Proc Natl Acad Sci U S A</source>. (<year>1998</year>) <volume>95</volume>:<fpage>9578</fpage>&#x2013;<lpage>83</lpage>. doi: <pub-id pub-id-type="doi">10.1073/pnas.95.16.9578</pub-id>, PMID: <pub-id pub-id-type="pmid">9689123</pub-id></citation></ref>
<ref id="ref24"><label>24.</label><citation citation-type="other"><person-group person-group-type="author"><collab id="coll4">NCBI</collab></person-group> (<year>2016</year>). The NCBI pathogen detection project [internet]. Bethesda (MD): National Library of Medicine (US), National Center for Biotechnology Information. Available at: <ext-link xlink:href="https://www.ncbi.nlm.nih.gov/pathogens/" ext-link-type="uri">https://www.ncbi.nlm.nih.gov/pathogens/</ext-link></citation></ref>
<ref id="ref25"><label>25.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dixit</surname><given-names>A</given-names></name> <name><surname>Freschi</surname><given-names>L</given-names></name> <name><surname>Vargas</surname><given-names>R</given-names></name> <name><surname>Calderon</surname><given-names>R</given-names></name> <name><surname>Sacchettini</surname><given-names>J</given-names></name> <name><surname>Drobniewski</surname><given-names>F</given-names></name> <etal/></person-group>. <article-title>Whole genome sequencing identifies bacterial factors affecting transmission of multidrug-resistant tuberculosis in a high-prevalence setting</article-title>. <source>Sci Rep</source>. (<year>2019</year>) <volume>9</volume>:<fpage>5602</fpage>. doi: <pub-id pub-id-type="doi">10.1038/s41598-019-41967-8</pub-id></citation></ref>
<ref id="ref26"><label>26.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bayowa</surname><given-names>JR</given-names></name> <name><surname>Kalyango</surname><given-names>JN</given-names></name> <name><surname>Baluku</surname><given-names>JB</given-names></name> <name><surname>Katuramu</surname><given-names>R</given-names></name> <name><surname>Ssendikwanawa</surname><given-names>E</given-names></name> <name><surname>Zalwango</surname><given-names>JF</given-names></name> <etal/></person-group>. <article-title>Mortality rate and associated factors among patients co-infected with drug resistant tuberculosis/HIV at Mulago National Referral Hospital, Uganda, a retrospective cohort study</article-title>. <source>PLOS Glob Public Health</source>. (<year>2023</year>) <volume>3</volume>:<fpage>e0001020</fpage>. doi: <pub-id pub-id-type="doi">10.1371/journal.pgph.0001020</pub-id>, PMID: <pub-id pub-id-type="pmid">37410761</pub-id></citation></ref>
<ref id="ref27"><label>27.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chem</surname><given-names>ED</given-names></name> <name><surname>Van Hout</surname><given-names>MC</given-names></name> <name><surname>Hope</surname><given-names>V</given-names></name></person-group>. <article-title>Treatment outcomes and antiretroviral uptake in multidrug-resistant tuberculosis and HIV co-infected patients in sub Saharan Africa: a systematic review and meta-analysis</article-title>. <source>BMC Infect Dis</source>. (<year>2019</year>) <volume>19</volume>:<fpage>723</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s12879-019-4317-4</pub-id></citation></ref>
<ref id="ref28"><label>28.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Isaakidis</surname><given-names>P</given-names></name> <name><surname>Casas</surname><given-names>EC</given-names></name> <name><surname>Das</surname><given-names>M</given-names></name> <name><surname>Tseretopoulou</surname><given-names>X</given-names></name> <name><surname>Ntzani</surname><given-names>EE</given-names></name> <name><surname>Ford</surname><given-names>N</given-names></name></person-group>. <article-title>Treatment outcomes for HIV and MDR-TB co-infected adults and children: systematic review and meta-analysis</article-title>. <source>Int J Tuberc Lung Dis</source>. (<year>2015</year>) <volume>19</volume>:<fpage>969</fpage>&#x2013;<lpage>78</lpage>. doi: <pub-id pub-id-type="doi">10.5588/ijtld.15.0123</pub-id>, PMID: <pub-id pub-id-type="pmid">26162364</pub-id></citation></ref>
<ref id="ref29"><label>29.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lelisho</surname><given-names>ME</given-names></name> <name><surname>Wotale</surname><given-names>TW</given-names></name> <name><surname>Tareke</surname><given-names>SA</given-names></name> <name><surname>Alemu</surname><given-names>BD</given-names></name> <name><surname>Hassen</surname><given-names>SS</given-names></name> <name><surname>Yemane</surname><given-names>DM</given-names></name> <etal/></person-group>. <article-title>Survival rate and predictors of mortality among TB/HIV co-infected adult patients: retrospective cohort study</article-title>. <source>Sci Rep</source>. (<year>2022</year>) <volume>12</volume>:<fpage>18360</fpage>. doi: <pub-id pub-id-type="doi">10.1038/s41598-022-23316-4</pub-id>, PMID: <pub-id pub-id-type="pmid">36319734</pub-id></citation></ref>
<ref id="ref30"><label>30.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sultana</surname><given-names>ZZ</given-names></name> <name><surname>Hoque</surname><given-names>FU</given-names></name> <name><surname>Beyene</surname><given-names>J</given-names></name> <name><surname>Akhlak-Ul-Islam</surname><given-names>M</given-names></name> <name><surname>Khan</surname><given-names>MHR</given-names></name> <name><surname>Ahmed</surname><given-names>S</given-names></name> <etal/></person-group>. <article-title>HIV infection and multidrug resistant tuberculosis: a systematic review and meta-analysis</article-title>. <source>BMC Infect Dis</source>. (<year>2021</year>) <volume>21</volume>:<fpage>51</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s12879-020-05749-2</pub-id>, PMID: <pub-id pub-id-type="pmid">33430786</pub-id></citation></ref>
<ref id="ref31"><label>31.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hu</surname><given-names>Y</given-names></name> <name><surname>Movahedzadeh</surname><given-names>F</given-names></name> <name><surname>Stoker</surname><given-names>NG</given-names></name> <name><surname>Coates</surname><given-names>AR</given-names></name></person-group>. <article-title>Deletion of the <italic>Mycobacterium tuberculosis</italic> alpha-crystallin-like hspX gene causes increased bacterial growth in vivo</article-title>. <source>Infect Immun</source>. (<year>2006</year>) <volume>74</volume>:<fpage>861</fpage>&#x2013;<lpage>8</lpage>. doi: <pub-id pub-id-type="doi">10.1128/IAI.74.2.861-868.2006</pub-id>, PMID: <pub-id pub-id-type="pmid">16428728</pub-id></citation></ref>
<ref id="ref32"><label>32.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mushtaq</surname><given-names>K</given-names></name> <name><surname>Sheikh</surname><given-names>JA</given-names></name> <name><surname>Amir</surname><given-names>M</given-names></name> <name><surname>Khan</surname><given-names>N</given-names></name> <name><surname>Singh</surname><given-names>B</given-names></name> <name><surname>Agrewala</surname><given-names>JN</given-names></name></person-group>. <article-title>Rv2031c of <italic>Mycobacterium tuberculosis</italic>: a master regulator of Rv2028-Rv2031 (HspX) operon</article-title>. <source>Front Microbiol</source>. (<year>2015</year>) <volume>6</volume>:<fpage>351</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fmicb.2015.00351</pub-id>, PMID: <pub-id pub-id-type="pmid">25964780</pub-id></citation></ref>
<ref id="ref33"><label>33.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Park</surname><given-names>HD</given-names></name> <name><surname>Guinn</surname><given-names>KM</given-names></name> <name><surname>Harrell</surname><given-names>MI</given-names></name> <name><surname>Liao</surname><given-names>R</given-names></name> <name><surname>Voskuil</surname><given-names>MI</given-names></name> <name><surname>Tompa</surname><given-names>M</given-names></name> <etal/></person-group>. <article-title>Rv3133c/dosR is a transcription factor that mediates the hypoxic response of <italic>Mycobacterium tuberculosis</italic></article-title>. <source>Mol Microbiol</source>. (<year>2003</year>) <volume>48</volume>:<fpage>833</fpage>&#x2013;<lpage>43</lpage>. doi: <pub-id pub-id-type="doi">10.1046/j.1365-2958.2003.03474.x</pub-id>, PMID: <pub-id pub-id-type="pmid">12694625</pub-id></citation></ref>
<ref id="ref34"><label>34.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gao</surname><given-names>CH</given-names></name> <name><surname>Yang</surname><given-names>M</given-names></name> <name><surname>He</surname><given-names>ZG</given-names></name></person-group>. <article-title>Characterization of a novel ArsR-like regulator encoded by Rv2034 in <italic>Mycobacterium tuberculosis</italic></article-title>. <source>PLoS One</source>. (<year>2012</year>) <volume>7</volume>:<fpage>e36255</fpage>. doi: <pub-id pub-id-type="doi">10.1371/journal.pone.0036255</pub-id>, PMID: <pub-id pub-id-type="pmid">22558408</pub-id></citation></ref>
<ref id="ref35"><label>35.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Betts</surname><given-names>JC</given-names></name> <name><surname>Lukey</surname><given-names>PT</given-names></name> <name><surname>Robb</surname><given-names>LC</given-names></name> <name><surname>McAdam</surname><given-names>RA</given-names></name> <name><surname>Duncan</surname><given-names>K</given-names></name></person-group>. <article-title>Evaluation of a nutrient starvation model of <italic>Mycobacterium tuberculosis</italic> persistence by gene and protein expression profiling</article-title>. <source>Mol Microbiol</source>. (<year>2002</year>) <volume>43</volume>:<fpage>717</fpage>&#x2013;<lpage>31</lpage>. doi: <pub-id pub-id-type="doi">10.1046/j.1365-2958.2002.02779.x</pub-id>, PMID: <pub-id pub-id-type="pmid">11929527</pub-id></citation></ref>
<ref id="ref36"><label>36.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Barreto</surname><given-names>C</given-names></name> <name><surname>de Oliveira</surname><given-names>M</given-names></name> <name><surname>Balan</surname><given-names>A</given-names></name></person-group>. <article-title>The ATP-binding cassette (ABC) transport Systems in <italic>Mycobacterium tuberculosis</italic>: structure, function, and possible targets for therapeutics</article-title>. <source>Biology</source>. (<year>2020</year>) <volume>9</volume>:<fpage>443</fpage>. doi: <pub-id pub-id-type="doi">10.3390/biology9120443</pub-id>, PMID: <pub-id pub-id-type="pmid">33291531</pub-id></citation></ref>
<ref id="ref37"><label>37.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Garbe</surname><given-names>TR</given-names></name> <name><surname>Hibler</surname><given-names>NS</given-names></name> <name><surname>Deretic</surname><given-names>V</given-names></name></person-group>. <article-title>Response to reactive nitrogen intermediates in <italic>Mycobacterium tuberculosis</italic>: induction of the 16-kilodalton alpha-crystallin homolog by exposure to nitric oxide donors</article-title>. <source>Infect Immun</source>. (<year>1999</year>) <volume>67</volume>:<fpage>460</fpage>&#x2013;<lpage>5</lpage>. doi: <pub-id pub-id-type="doi">10.1128/IAI.67.1.460-465.1999</pub-id>, PMID: <pub-id pub-id-type="pmid">9864257</pub-id></citation></ref>
<ref id="ref38"><label>38.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kim</surname><given-names>SY</given-names></name> <name><surname>Lee</surname><given-names>BS</given-names></name> <name><surname>Shin</surname><given-names>SJ</given-names></name> <name><surname>Kim</surname><given-names>HJ</given-names></name> <name><surname>Park</surname><given-names>JK</given-names></name></person-group>. <article-title>Differentially expressed genes in <italic>Mycobacterium tuberculosis</italic> H37Rv under mild acidic and hypoxic conditions</article-title>. <source>J Med Microbiol</source>. (<year>2008</year>) <volume>57</volume>:<fpage>1473</fpage>&#x2013;<lpage>80</lpage>. doi: <pub-id pub-id-type="doi">10.1099/jmm.0.2008/001545-0</pub-id>, PMID: <pub-id pub-id-type="pmid">19018016</pub-id></citation></ref>
<ref id="ref39"><label>39.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kumari</surname><given-names>B</given-names></name> <name><surname>Saini</surname><given-names>V</given-names></name> <name><surname>Kaur</surname><given-names>J</given-names></name> <name><surname>Kaur</surname><given-names>J</given-names></name></person-group>. <article-title>Rv2037c, a stress induced conserved hypothetical protein of <italic>Mycobacterium tuberculosis</italic>, is a phospholipase: role in cell wall modulation and intracellular survival</article-title>. <source>Int J Biol Macromol</source>. (<year>2020</year>) <volume>153</volume>:<fpage>817</fpage>&#x2013;<lpage>35</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.ijbiomac.2020.03.037</pub-id>, PMID: <pub-id pub-id-type="pmid">32165202</pub-id></citation></ref>
<ref id="ref40"><label>40.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Saxena</surname><given-names>A</given-names></name> <name><surname>Srivastava</surname><given-names>V</given-names></name> <name><surname>Srivastava</surname><given-names>R</given-names></name> <name><surname>Srivastava</surname><given-names>BS</given-names></name></person-group>. <article-title>Identification of genes of <italic>Mycobacterium tuberculosis</italic> upregulated during anaerobic persistence by fluorescence and kanamycin resistance selection</article-title>. <source>Tuberculosis</source>. (<year>2008</year>) <volume>88</volume>:<fpage>518</fpage>&#x2013;<lpage>25</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.tube.2008.01.003</pub-id>, PMID: <pub-id pub-id-type="pmid">18434250</pub-id></citation></ref>
<ref id="ref41"><label>41.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Snasel</surname><given-names>J</given-names></name> <name><surname>Machova</surname><given-names>I</given-names></name> <name><surname>Solinova</surname><given-names>V</given-names></name> <name><surname>Kasicka</surname><given-names>V</given-names></name> <name><surname>Krecmerova</surname><given-names>M</given-names></name> <name><surname>Pichova</surname><given-names>I</given-names></name></person-group>. <article-title>Phosphofructokinases a and B from <italic>Mycobacterium tuberculosis</italic> display different catalytic properties and allosteric regulation</article-title>. <source>Int J Mol Sci</source>. (<year>2021</year>) <volume>22</volume>. doi: <pub-id pub-id-type="doi">10.3390/ijms22031483</pub-id>, PMID: <pub-id pub-id-type="pmid">33540748</pub-id></citation></ref>
<ref id="ref42"><label>42.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jamwal</surname><given-names>SV</given-names></name> <name><surname>Mehrotra</surname><given-names>P</given-names></name> <name><surname>Singh</surname><given-names>A</given-names></name> <name><surname>Siddiqui</surname><given-names>Z</given-names></name> <name><surname>Basu</surname><given-names>A</given-names></name> <name><surname>Rao</surname><given-names>KV</given-names></name></person-group>. <article-title>Mycobacterial escape from macrophage phagosomes to the cytoplasm represents an alternate adaptation mechanism</article-title>. <source>Sci Rep</source>. (<year>2016</year>) <volume>6</volume>:<fpage>23089</fpage>. doi: <pub-id pub-id-type="doi">10.1038/srep23089</pub-id>, PMID: <pub-id pub-id-type="pmid">26980157</pub-id></citation></ref>
<ref id="ref43"><label>43.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Legesse</surname><given-names>M</given-names></name> <name><surname>Ameni</surname><given-names>G</given-names></name> <name><surname>Medhin</surname><given-names>G</given-names></name> <name><surname>Mamo</surname><given-names>G</given-names></name> <name><surname>Franken</surname><given-names>KL</given-names></name> <name><surname>Ottenhoff</surname><given-names>TH</given-names></name> <etal/></person-group>. <article-title>IgA response to ESAT-6/CFP-10 and Rv2031 antigens varies in patients with culture-confirmed pulmonary tuberculosis, healthy <italic>Mycobacterium tuberculosis</italic>-infected and non-infected individuals in a tuberculosis endemic setting</article-title>. <source>Ethiopia Scand J Immunol</source>. (<year>2013</year>) <volume>78</volume>:<fpage>266</fpage>&#x2013;<lpage>74</lpage>. doi: <pub-id pub-id-type="doi">10.1111/sji.12080</pub-id>, PMID: <pub-id pub-id-type="pmid">23713613</pub-id></citation></ref>
<ref id="ref44"><label>44.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Meier</surname><given-names>NR</given-names></name> <name><surname>Battegay</surname><given-names>M</given-names></name> <name><surname>Ottenhoff</surname><given-names>THM</given-names></name> <name><surname>Furrer</surname><given-names>H</given-names></name> <name><surname>Nemeth</surname><given-names>J</given-names></name> <name><surname>Ritz</surname><given-names>N</given-names></name></person-group>. <article-title>HIV-infected patients developing tuberculosis disease show early changes in the immune response to novel <italic>Mycobacterium tuberculosis</italic> antigens</article-title>. <source>Front Immunol</source>. (<year>2021</year>) <volume>12</volume>:<fpage>620622</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fimmu.2021.620622</pub-id>, PMID: <pub-id pub-id-type="pmid">33777000</pub-id></citation></ref>
<ref id="ref45"><label>45.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Amir</surname><given-names>M</given-names></name> <name><surname>Aqdas</surname><given-names>M</given-names></name> <name><surname>Nadeem</surname><given-names>S</given-names></name> <name><surname>Siddiqui</surname><given-names>KF</given-names></name> <name><surname>Khan</surname><given-names>N</given-names></name> <name><surname>Sheikh</surname><given-names>JA</given-names></name> <etal/></person-group>. <article-title>Diametric role of the latency-associated protein Acr1 of <italic>Mycobacterium tuberculosis</italic> in modulating the functionality of pre-and post-maturational stages of dendritic cells</article-title>. <source>Front Immunol</source>. (<year>2017</year>) <volume>8</volume>:<fpage>624</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fimmu.2017.00624</pub-id>, PMID: <pub-id pub-id-type="pmid">28611779</pub-id></citation></ref>
<ref id="ref46"><label>46.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hu</surname><given-names>Y</given-names></name> <name><surname>Coates</surname><given-names>AR</given-names></name></person-group>. <article-title><italic>Mycobacterium tuberculosis</italic> acg gene is required for growth and virulence in vivo</article-title>. <source>PLoS One</source>. (<year>2011</year>) <volume>6</volume>:<fpage>e20958</fpage>. doi: <pub-id pub-id-type="doi">10.1371/journal.pone.0020958</pub-id>, PMID: <pub-id pub-id-type="pmid">21687631</pub-id></citation></ref>
<ref id="ref47"><label>47.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kim</surname><given-names>S-Y</given-names></name></person-group>. <article-title>A-rum Shin, Byung-Soo Lee, Hwa-Jung Kim2 Bo young Jeon, sang-Nae Cho, Jeong-Kyu Park, sung Jae Shin characterization of immune responses to <italic>Mycobacterium tuberculosis</italic> Rv2041c protein</article-title>. <source>J Bacteriol Virol</source>. (<year>2009</year>) <volume>39</volume>:<fpage>183</fpage>&#x2013;<lpage>93</lpage>. doi: <pub-id pub-id-type="doi">10.4167/jbv.2009.39.3.183</pub-id></citation></ref>
<ref id="ref48"><label>48.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sassetti</surname><given-names>CM</given-names></name> <name><surname>Boyd</surname><given-names>DH</given-names></name> <name><surname>Rubin</surname><given-names>EJ</given-names></name></person-group>. <article-title>Genes required for mycobacterial growth defined by high density mutagenesis</article-title>. <source>Mol Microbiol</source>. (<year>2003</year>) <volume>48</volume>:<fpage>77</fpage>&#x2013;<lpage>84</lpage>. doi: <pub-id pub-id-type="doi">10.1046/j.1365-2958.2003.03425.x</pub-id>, PMID: <pub-id pub-id-type="pmid">12657046</pub-id></citation></ref>
<ref id="ref49"><label>49.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>De la Torre</surname><given-names>LI</given-names></name> <name><surname>Vergara Meza</surname><given-names>JG</given-names></name> <name><surname>Cabarca</surname><given-names>S</given-names></name> <name><surname>Costa-Martins</surname><given-names>AG</given-names></name> <name><surname>Balan</surname><given-names>A</given-names></name></person-group>. <article-title>Comparison of carbohydrate ABC importers from <italic>Mycobacterium tuberculosis</italic></article-title>. <source>BMC Genomics</source>. (<year>2021</year>) <volume>22</volume>:<fpage>841</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s12864-021-07972-w</pub-id>, PMID: <pub-id pub-id-type="pmid">34798821</pub-id></citation></ref>
<ref id="ref50"><label>50.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Godfroid</surname><given-names>M</given-names></name> <name><surname>Dagan</surname><given-names>T</given-names></name> <name><surname>Merker</surname><given-names>M</given-names></name> <name><surname>Kohl</surname><given-names>TA</given-names></name> <name><surname>Diel</surname><given-names>R</given-names></name> <name><surname>Maurer</surname><given-names>FP</given-names></name> <etal/></person-group>. <article-title>Insertion and deletion evolution reflects antibiotics selection pressure in a <italic>Mycobacterium tuberculosis</italic> outbreak</article-title>. <source>PLoS Pathog</source>. (<year>2020</year>) <volume>16</volume>:<fpage>e1008357</fpage>. doi: <pub-id pub-id-type="doi">10.1371/journal.ppat.1008357</pub-id>, PMID: <pub-id pub-id-type="pmid">32997707</pub-id></citation></ref>
<ref id="ref51"><label>51.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ford</surname><given-names>CB</given-names></name> <name><surname>Shah</surname><given-names>RR</given-names></name> <name><surname>Maeda</surname><given-names>MK</given-names></name> <name><surname>Gagneux</surname><given-names>S</given-names></name> <name><surname>Murray</surname><given-names>MB</given-names></name> <name><surname>Cohen</surname><given-names>T</given-names></name> <etal/></person-group>. <article-title><italic>Mycobacterium tuberculosis</italic> mutation rate estimates from different lineages predict substantial differences in the emergence of drug-resistant tuberculosis</article-title>. <source>Nat Genet</source>. (<year>2013</year>) <volume>45</volume>:<fpage>784</fpage>&#x2013;<lpage>90</lpage>. doi: <pub-id pub-id-type="doi">10.1038/ng.2656</pub-id>, PMID: <pub-id pub-id-type="pmid">23749189</pub-id></citation></ref>
<ref id="ref52"><label>52.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hakamata</surname><given-names>M</given-names></name> <name><surname>Takihara</surname><given-names>H</given-names></name> <name><surname>Iwamoto</surname><given-names>T</given-names></name> <name><surname>Tamaru</surname><given-names>A</given-names></name> <name><surname>Hashimoto</surname><given-names>A</given-names></name> <name><surname>Tanaka</surname><given-names>T</given-names></name> <etal/></person-group>. <article-title>Higher genome mutation rates of Beijing lineage of <italic>Mycobacterium tuberculosis</italic> during human infection</article-title>. <source>Sci Rep</source>. (<year>2020</year>) <volume>10</volume>:<fpage>17997</fpage>. doi: <pub-id pub-id-type="doi">10.1038/s41598-020-75028-2</pub-id>, PMID: <pub-id pub-id-type="pmid">33093577</pub-id></citation></ref>
<ref id="ref53"><label>53.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname><given-names>Q</given-names></name> <name><surname>Wang</surname><given-names>D</given-names></name> <name><surname>Martinez</surname><given-names>L</given-names></name> <name><surname>Lu</surname><given-names>P</given-names></name> <name><surname>Zhu</surname><given-names>L</given-names></name> <name><surname>Lu</surname><given-names>W</given-names></name> <etal/></person-group>. <article-title><italic>Mycobacterium tuberculosis</italic> Beijing genotype strains and unfavourable treatment outcomes: a systematic review and meta-analysis</article-title>. <source>Clin Microbiol Infect</source>. (<year>2020</year>) <volume>26</volume>:<fpage>180</fpage>&#x2013;<lpage>8</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.cmi.2019.07.016</pub-id>, PMID: <pub-id pub-id-type="pmid">31336202</pub-id></citation></ref>
<ref id="ref54"><label>54.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ribeiro</surname><given-names>SC</given-names></name> <name><surname>Gomes</surname><given-names>LL</given-names></name> <name><surname>Amaral</surname><given-names>EP</given-names></name> <name><surname>Andrade</surname><given-names>MR</given-names></name> <name><surname>Almeida</surname><given-names>FM</given-names></name> <name><surname>Rezende</surname><given-names>AL</given-names></name> <etal/></person-group>. <article-title><italic>Mycobacterium tuberculosis</italic> strains of the modern sublineage of the Beijing family are more likely to display increased virulence than strains of the ancient sublineage</article-title>. <source>J Clin Microbiol</source>. (<year>2014</year>) <volume>52</volume>:<fpage>2615</fpage>&#x2013;<lpage>24</lpage>. doi: <pub-id pub-id-type="doi">10.1128/JCM.00498-14</pub-id>, PMID: <pub-id pub-id-type="pmid">24829250</pub-id></citation></ref>
<ref id="ref55"><label>55.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Singh</surname><given-names>R</given-names></name> <name><surname>Dwivedi</surname><given-names>SP</given-names></name> <name><surname>Gaharwar</surname><given-names>US</given-names></name> <name><surname>Meena</surname><given-names>R</given-names></name> <name><surname>Rajamani</surname><given-names>P</given-names></name> <name><surname>Prasad</surname><given-names>T</given-names></name></person-group>. <article-title>Recent updates on drug resistance in <italic>Mycobacterium tuberculosis</italic></article-title>. <source>J Appl Microbiol</source>. (<year>2020</year>) <volume>128</volume>:<fpage>1547</fpage>&#x2013;<lpage>67</lpage>. doi: <pub-id pub-id-type="doi">10.1111/jam.14478</pub-id>, PMID: <pub-id pub-id-type="pmid">31595643</pub-id></citation></ref>
<ref id="ref56"><label>56.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Safi</surname><given-names>H</given-names></name> <name><surname>Gopal</surname><given-names>P</given-names></name> <name><surname>Lingaraju</surname><given-names>S</given-names></name> <name><surname>Ma</surname><given-names>S</given-names></name> <name><surname>Levine</surname><given-names>C</given-names></name> <name><surname>Dartois</surname><given-names>V</given-names></name> <etal/></person-group>. <article-title>Phase variation in <italic>Mycobacterium tuberculosis</italic> glpK produces transiently heritable drug tolerance</article-title>. <source>Proc Natl Acad Sci U S A</source>. (<year>2019</year>) <volume>116</volume>:<fpage>19665</fpage>&#x2013;<lpage>74</lpage>. doi: <pub-id pub-id-type="doi">10.1073/pnas.1907631116</pub-id>, PMID: <pub-id pub-id-type="pmid">31488707</pub-id></citation></ref>
<ref id="ref57"><label>57.</label><citation citation-type="other"><person-group person-group-type="author"><collab id="coll5">World Health Organization. Communicable Diseases Cluster</collab></person-group> (<year>1999</year>). What is DOTS?: A guide to understanding the WHO-recommended TB control strategy known as DOTS. World Health Organization. Available at: <ext-link xlink:href="https://apps.who.int/iris/handle/10665/65979" ext-link-type="uri">https://apps.who.int/iris/handle/10665/65979</ext-link></citation></ref>
<ref id="ref58"><label>58.</label><citation citation-type="book"><person-group person-group-type="author"><name><surname>Kent</surname><given-names>PT</given-names></name> <name><surname>Kubica</surname><given-names>GP</given-names></name></person-group>. <source>Public Health Mycobacteriology: A Guide for the Level III Laboratory</source>. <publisher-loc>Atlanta</publisher-loc>: <publisher-name>Centers for Disease Control</publisher-name> (<year>1985</year>).</citation></ref>
<ref id="ref59"><label>59.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Marks</surname><given-names>J</given-names></name></person-group>. <article-title>A system for the examination of tubercle bacilli and other mycobacteria</article-title>. <source>Tubercle</source>. (<year>1976</year>) <volume>57</volume>:<fpage>207</fpage>&#x2013;<lpage>25</lpage>. doi: <pub-id pub-id-type="doi">10.1016/0041-3879(76)90030-1</pub-id>, PMID: <pub-id pub-id-type="pmid">62439</pub-id></citation></ref>
<ref id="ref60"><label>60.</label><citation citation-type="other"><person-group person-group-type="author"><name><surname>Collins</surname><given-names>C. H.</given-names></name> <name><surname>Grange</surname><given-names>J. M.</given-names></name> <name><surname>Yates</surname><given-names>M.D.</given-names></name></person-group> (<year>1985</year>). Organization and Practice in Tuberculosis Bacteriology.</citation></ref>
<ref id="ref61"><label>61.</label><citation citation-type="book"><person-group person-group-type="author"><collab id="coll6">European Centre for Disease Prevention and Control</collab></person-group>. <source>Handbook on Tuberculosis Laboratory Diagnostic Methods in the European Union&#x2014;Updated 2018</source>. <publisher-loc>Stockholm</publisher-loc>: <publisher-name>ECDC</publisher-name> (<year>2018</year>).</citation></ref>
<ref id="ref62"><label>62.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Marks</surname><given-names>J</given-names></name></person-group>. <article-title>A 'Stepped Ph' technique for the estimation of pyrazinamide sensitivity</article-title>. <source>Tubercle</source>. (<year>1964</year>) <volume>45</volume>:<fpage>47</fpage>&#x2013;<lpage>50</lpage>. doi: <pub-id pub-id-type="doi">10.1016/s0041-3879(64)80087-8</pub-id>, PMID: <pub-id pub-id-type="pmid">14129199</pub-id></citation></ref>
<ref id="ref63"><label>63.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cabibbe</surname><given-names>AM</given-names></name> <name><surname>Walker</surname><given-names>TM</given-names></name> <name><surname>Niemann</surname><given-names>S</given-names></name> <name><surname>Cirillo</surname><given-names>DM</given-names></name></person-group>. <article-title>Whole genome sequencing of <italic>Mycobacterium tuberculosis</italic></article-title>. <source>Eur Respir J</source>. (<year>2018</year>) <volume>52</volume>:<fpage>1801163</fpage>. doi: <pub-id pub-id-type="doi">10.1183/13993003.01163-2018</pub-id>, PMID: <pub-id pub-id-type="pmid">30209198</pub-id></citation></ref>
<ref id="ref64"><label>64.</label><citation citation-type="other"><person-group person-group-type="author"><name><surname>Andrews</surname><given-names>S.</given-names></name></person-group> (<year>2010</year>). FASTQC. A quality control tool for high throughput sequence data.</citation></ref>
<ref id="ref65"><label>65.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wood</surname><given-names>DE</given-names></name> <name><surname>Lu</surname><given-names>J</given-names></name> <name><surname>Langmead</surname><given-names>B</given-names></name></person-group>. <article-title>Improved metagenomic analysis with kraken 2</article-title>. <source>Genome Biol</source>. (<year>2019</year>) <volume>20</volume>:<fpage>257</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s13059-019-1891-0</pub-id>, PMID: <pub-id pub-id-type="pmid">31779668</pub-id></citation></ref>
<ref id="ref66"><label>66.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kohl</surname><given-names>TA</given-names></name> <name><surname>Utpatel</surname><given-names>C</given-names></name> <name><surname>Schleusener</surname><given-names>V</given-names></name> <name><surname>De Filippo</surname><given-names>MR</given-names></name> <name><surname>Beckert</surname><given-names>P</given-names></name> <name><surname>Cirillo</surname><given-names>DM</given-names></name> <etal/></person-group>. <article-title>MTBseq: a comprehensive pipeline for whole genome sequence analysis of <italic>Mycobacterium tuberculosis</italic> complex isolates</article-title>. <source>PeerJ</source>. (<year>2018</year>) <volume>6</volume>:<fpage>e5895</fpage>. doi: <pub-id pub-id-type="doi">10.7717/peerj.5895</pub-id>, PMID: <pub-id pub-id-type="pmid">30479891</pub-id></citation></ref>
<ref id="ref67"><label>67.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Coll</surname><given-names>F</given-names></name> <name><surname>McNerney</surname><given-names>R</given-names></name> <name><surname>Guerra-Assuncao</surname><given-names>JA</given-names></name> <name><surname>Glynn</surname><given-names>JR</given-names></name> <name><surname>Perdigao</surname><given-names>J</given-names></name> <name><surname>Viveiros</surname><given-names>M</given-names></name> <etal/></person-group>. <article-title>A robust SNP barcode for typing <italic>Mycobacterium tuberculosis</italic> complex strains</article-title>. <source>Nat Commun</source>. (<year>2014</year>) <volume>5</volume>:<fpage>4812</fpage>. doi: <pub-id pub-id-type="doi">10.1038/ncomms5812</pub-id>, PMID: <pub-id pub-id-type="pmid">25176035</pub-id></citation></ref>
<ref id="ref68"><label>68.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>O'Leary</surname><given-names>NA</given-names></name> <name><surname>Wright</surname><given-names>MW</given-names></name> <name><surname>Brister</surname><given-names>JR</given-names></name> <name><surname>Ciufo</surname><given-names>S</given-names></name> <name><surname>Haddad</surname><given-names>D</given-names></name> <name><surname>McVeigh</surname><given-names>R</given-names></name> <etal/></person-group>. <article-title>Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation</article-title>. <source>Nucleic Acids Res</source>. (<year>2016</year>) <volume>44</volume>:<fpage>D733</fpage>&#x2013;<lpage>45</lpage>. doi: <pub-id pub-id-type="doi">10.1093/nar/gkv1189</pub-id>, PMID: <pub-id pub-id-type="pmid">26553804</pub-id></citation></ref>
<ref id="ref69"><label>69.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Marcais</surname><given-names>G</given-names></name> <name><surname>Delcher</surname><given-names>AL</given-names></name> <name><surname>Phillippy</surname><given-names>AM</given-names></name> <name><surname>Coston</surname><given-names>R</given-names></name> <name><surname>Salzberg</surname><given-names>SL</given-names></name> <name><surname>Zimin</surname><given-names>A</given-names></name></person-group>. <article-title>MUMmer4: A fast and versatile genome alignment system</article-title>. <source>PLoS Comput Biol</source>. (<year>2018</year>) <volume>14</volume>:<fpage>e1005944</fpage>. doi: <pub-id pub-id-type="doi">10.1371/journal.pcbi.1005944</pub-id>, PMID: <pub-id pub-id-type="pmid">29373581</pub-id></citation></ref>
<ref id="ref70"><label>70.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Leinonen</surname><given-names>R</given-names></name> <name><surname>Sugawara</surname><given-names>H</given-names></name> <name><surname>Shumway</surname><given-names>M</given-names></name></person-group>. <article-title>International nucleotide sequence database, collaboration the sequence read archive</article-title>. <source>Nucleic Acids Res</source>. (<year>2011</year>) <volume>39</volume>:<fpage>D19</fpage>&#x2013;<lpage>21</lpage>. doi: <pub-id pub-id-type="doi">10.1093/nar/gkq1019</pub-id>, PMID: <pub-id pub-id-type="pmid">21062823</pub-id></citation></ref>
<ref id="ref71"><label>71.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cameron</surname><given-names>DL</given-names></name> <name><surname>Baber</surname><given-names>J</given-names></name> <name><surname>Shale</surname><given-names>C</given-names></name> <name><surname>Valle-Inclan</surname><given-names>JE</given-names></name> <name><surname>Besselink</surname><given-names>N</given-names></name> <name><surname>van Hoeck</surname><given-names>A</given-names></name> <etal/></person-group>. <article-title>GRIDSS2: comprehensive characterisation of somatic structural variation using single breakend variants and structural variant phasing</article-title>. <source>Genome Biol</source>. (<year>2021</year>) <volume>22</volume>:<fpage>202</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s13059-021-02423-x</pub-id>, PMID: <pub-id pub-id-type="pmid">34253237</pub-id></citation></ref>
<ref id="ref72"><label>72.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cameron</surname><given-names>DL</given-names></name> <name><surname>Schroder</surname><given-names>J</given-names></name> <name><surname>Penington</surname><given-names>JS</given-names></name> <name><surname>Do</surname><given-names>H</given-names></name> <name><surname>Molania</surname><given-names>R</given-names></name> <name><surname>Dobrovic</surname><given-names>A</given-names></name> <etal/></person-group>. <article-title>GRIDSS: sensitive and specific genomic rearrangement detection using positional de Bruijn graph assembly</article-title>. <source>Genome Res</source>. (<year>2017</year>) <volume>27</volume>:<fpage>2050</fpage>&#x2013;<lpage>60</lpage>. doi: <pub-id pub-id-type="doi">10.1101/gr.222109.117</pub-id>, PMID: <pub-id pub-id-type="pmid">29097403</pub-id></citation></ref>
<ref id="ref73"><label>73.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Danecek</surname><given-names>P</given-names></name> <name><surname>Bonfield</surname><given-names>JK</given-names></name> <name><surname>Liddle</surname><given-names>J</given-names></name> <name><surname>Marshall</surname><given-names>J</given-names></name> <name><surname>Ohan</surname><given-names>V</given-names></name> <name><surname>Pollard</surname><given-names>MO</given-names></name> <etal/></person-group>. <article-title>Twelve years of SAMtools and BCFtools</article-title>. <source>Gigascience</source>. (<year>2021</year>) <volume>10</volume>. doi: <pub-id pub-id-type="doi">10.1093/gigascience/giab008</pub-id></citation></ref>
<ref id="ref74"><label>74.</label><citation citation-type="other"><person-group person-group-type="author"><name><surname>Garrison</surname><given-names>E.</given-names></name> <name><surname>Marth</surname><given-names>G. T.</given-names></name></person-group> (<year>2012</year>). Haplotype-based variant detection from short-read sequencing. arXiv [Preprint]. doi: <pub-id pub-id-type="doi">10.48550/arXiv.1207.3907</pub-id></citation></ref>
<ref id="ref75"><label>75.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kohl</surname><given-names>TA</given-names></name> <name><surname>Diel</surname><given-names>R</given-names></name> <name><surname>Harmsen</surname><given-names>D</given-names></name> <name><surname>Rothganger</surname><given-names>J</given-names></name> <name><surname>Walter</surname><given-names>KM</given-names></name> <name><surname>Merker</surname><given-names>M</given-names></name> <etal/></person-group>. <article-title>Whole-genome-based <italic>Mycobacterium tuberculosis</italic> surveillance: a standardized, portable, and expandable approach</article-title>. <source>J Clin Microbiol</source>. (<year>2014</year>) <volume>52</volume>:<fpage>2479</fpage>&#x2013;<lpage>86</lpage>. doi: <pub-id pub-id-type="doi">10.1128/JCM.00567-14</pub-id>, PMID: <pub-id pub-id-type="pmid">24789177</pub-id></citation></ref>
<ref id="ref76"><label>76.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Danecek</surname><given-names>P</given-names></name> <name><surname>Auton</surname><given-names>A</given-names></name> <name><surname>Abecasis</surname><given-names>G</given-names></name> <name><surname>Albers</surname><given-names>CA</given-names></name> <name><surname>Banks</surname><given-names>E</given-names></name> <name><surname>DePristo</surname><given-names>MA</given-names></name> <etal/></person-group>. <article-title>Genomes project analysis, group the variant call format and VCFtools</article-title>. <source>Bioinformatics</source>. (<year>2011</year>) <volume>27</volume>:<fpage>2156</fpage>&#x2013;<lpage>8</lpage>. doi: <pub-id pub-id-type="doi">10.1093/bioinformatics/btr330</pub-id>, PMID: <pub-id pub-id-type="pmid">21653522</pub-id></citation></ref>
<ref id="ref77"><label>77.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Garrison</surname><given-names>E</given-names></name> <name><surname>Kronenberg</surname><given-names>ZN</given-names></name> <name><surname>Dawson</surname><given-names>ET</given-names></name> <name><surname>Pedersen</surname><given-names>BS</given-names></name> <name><surname>Prins</surname><given-names>P</given-names></name></person-group>. <article-title>A spectrum of free software tools for processing the VCF variant call format: vcflib, bio-VCF, cyvcf2, hts-nim and slivar</article-title>. <source>PLoS Comput Biol</source>. (<year>2022</year>) <volume>18</volume>:<fpage>e1009123</fpage>. doi: <pub-id pub-id-type="doi">10.1371/journal.pcbi.1009123</pub-id>, PMID: <pub-id pub-id-type="pmid">35639788</pub-id></citation></ref>
<ref id="ref78"><label>78.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Quinlan</surname><given-names>AR</given-names></name> <name><surname>Hall</surname><given-names>IM</given-names></name></person-group>. <article-title>BEDTools: a flexible suite of utilities for comparing genomic features</article-title>. <source>Bioinformatics</source>. (<year>2010</year>) <volume>26</volume>:<fpage>841</fpage>&#x2013;<lpage>2</lpage>. doi: <pub-id pub-id-type="doi">10.1093/bioinformatics/btq033</pub-id>, PMID: <pub-id pub-id-type="pmid">20110278</pub-id></citation></ref>
<ref id="ref79"><label>79.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhou</surname><given-names>Z</given-names></name> <name><surname>Alikhan</surname><given-names>NF</given-names></name> <name><surname>Sergeant</surname><given-names>MJ</given-names></name> <name><surname>Luhmann</surname><given-names>N</given-names></name> <name><surname>Vaz</surname><given-names>C</given-names></name> <name><surname>Francisco</surname><given-names>AP</given-names></name> <etal/></person-group>. <article-title>GrapeTree: visualization of core genomic relationships among 100,000 bacterial pathogens</article-title>. <source>Genome Res</source>. (<year>2018</year>) <volume>28</volume>:<fpage>1395</fpage>&#x2013;<lpage>404</lpage>. doi: <pub-id pub-id-type="doi">10.1101/gr.232397.117</pub-id>, PMID: <pub-id pub-id-type="pmid">30049790</pub-id></citation></ref>
<ref id="ref80"><label>80.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rubinstein</surname><given-names>M</given-names></name> <name><surname>Grossman</surname><given-names>R</given-names></name> <name><surname>Nissan</surname><given-names>I</given-names></name> <name><surname>Schwaber</surname><given-names>MJ</given-names></name> <name><surname>Carmeli</surname><given-names>Y</given-names></name> <name><surname>Kaidar-Shwartz</surname><given-names>H</given-names></name> <etal/></person-group>. <article-title><italic>Mycobacterium intracellulare</italic> subsp. chimaera from cardio surgery heating-cooling units and from clinical samples in Israel are genetically unrelated</article-title>. <source>Pathogens</source>. (<year>2021</year>) <volume>10</volume>:<fpage>1392</fpage>. doi: <pub-id pub-id-type="doi">10.3390/pathogens10111392</pub-id></citation></ref>
<ref id="ref81"><label>81.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gardner</surname><given-names>SN</given-names></name> <name><surname>Hall</surname><given-names>BG</given-names></name></person-group>. <article-title>When whole-genome alignments just won't work: kSNP v2 software for alignment-free SNP discovery and phylogenetics of hundreds of microbial genomes</article-title>. <source>PLoS One</source>. (<year>2013</year>) <volume>8</volume>:<fpage>e81760</fpage>. doi: <pub-id pub-id-type="doi">10.1371/journal.pone.0081760</pub-id></citation></ref>
<ref id="ref82"><label>82.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gardner</surname><given-names>SN</given-names></name> <name><surname>Slezak</surname><given-names>T</given-names></name> <name><surname>Hall</surname><given-names>BG</given-names></name></person-group>. <article-title>kSNP3.0: SNP detection and phylogenetic analysis of genomes without genome alignment or reference genome</article-title>. <source>Bioinformatics</source>. (<year>2015</year>) <volume>31</volume>:<fpage>2877</fpage>&#x2013;<lpage>8</lpage>. doi: <pub-id pub-id-type="doi">10.1093/bioinformatics/btv271</pub-id></citation></ref>
<ref id="ref83"><label>83.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Barrick</surname><given-names>JE</given-names></name> <name><surname>Colburn</surname><given-names>G</given-names></name> <name><surname>Deatherage</surname><given-names>DE</given-names></name> <name><surname>Traverse</surname><given-names>CC</given-names></name> <name><surname>Strand</surname><given-names>MD</given-names></name> <name><surname>Borges</surname><given-names>JJ</given-names></name> <etal/></person-group>. <article-title>Identifying structural variation in haploid microbial genomes from short-read resequencing data using breseq</article-title>. <source>BMC Genomics</source>. (<year>2014</year>) <volume>15</volume>:<fpage>1039</fpage>. doi: <pub-id pub-id-type="doi">10.1186/1471-2164-15-1039</pub-id>, PMID: <pub-id pub-id-type="pmid">25432719</pub-id></citation></ref>
<ref id="ref84"><label>84.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Deatherage</surname><given-names>DE</given-names></name> <name><surname>Barrick</surname><given-names>JE</given-names></name></person-group>. <article-title>Identification of mutations in laboratory-evolved microbes from next-generation sequencing data using breseq</article-title>. <source>Methods Mol Biol</source>. (<year>2014</year>) <volume>1151</volume>:<fpage>165</fpage>&#x2013;<lpage>88</lpage>. doi: <pub-id pub-id-type="doi">10.1007/978-1-4939-0554-6_12</pub-id>, PMID: <pub-id pub-id-type="pmid">24838886</pub-id></citation></ref>
<ref id="ref85"><label>85.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cingolani</surname><given-names>P</given-names></name> <name><surname>Platts</surname><given-names>A</given-names></name> <name><surname>Wang le</surname><given-names>L</given-names></name> <name><surname>Coon</surname><given-names>M</given-names></name> <name><surname>Nguyen</surname><given-names>T</given-names></name> <name><surname>Wang</surname><given-names>L</given-names></name> <etal/></person-group>. <article-title>A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of <italic>Drosophila melanogaster</italic> strain w1118; iso-2; iso-3</article-title>. <source>Fly</source>. (<year>2012</year>) <volume>6</volume>:<fpage>80</fpage>&#x2013;<lpage>92</lpage>. doi: <pub-id pub-id-type="doi">10.4161/fly.19695</pub-id>, PMID: <pub-id pub-id-type="pmid">22728672</pub-id></citation></ref>
</ref-list>
</back>
</article>