Multidrug-resistant tuberculosis in panama is driven by clonal expansion of a multidrug-resistant Mycobacterium tuberculosis strain related to the KZN extensively drug-resistant M. tuberculosis strain from South Africa

Heterogeneity of Mycobacterium tuberculosis Strains Circulating in Panama's Western Region

Tuberculosis remains a challenge in both rural and urban areas. Although a majority of countries display a higher burden in urban areas compared with rural areas, Panama continues to report the highest mortality rate in Central America. Urban areas, such as Panama City, report a high tuberculosis burden, whereas Panama's western region, including the provinces of Chiriquí, Bocas del Toro (both semiurban) and Ngäbe-Bugle (rural), show a lower burden. We aimed to identify highly transmitted Mycobacterium tuberculosis strains within rural and semiurban settings of Panama's western region during a 3-year period (2017, 2019, 2021). We randomly selected 87 M. tuberculosis isolates from a biobank from Panama's western region and analyzed them using allele-specific oligonucleotide polymerase chain reaction and 24-mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR). Our results show only 11.7% (10/85) of M. tuberculosis strains identified as prevalent A-Beijing, B-Haarlem, or C-LAM Strains. We found a low prevalence of A, B, and C M. tuberculosis strains in both rural and semirural settings compared with isolates collected from the Eastern Colon Province. MIRU-VNTR genotyping revealed a high degree of diversity with no clusters with single loci variation of ≥ 2 loci. These results support the notion that tuberculosis prevalence in the rural and semiurban western region of Panama are not due to previously described highly transmitted strains but is influenced instead by other health determinants, including poor health system access and a lack of systematic transmission chain monitoring. For remote rural and semiurban settings, we recommend allocating resources to reinforce efforts to prevent tuberculosis spread.

Investigating the intensity of social contacts associated with tuberculosis: a weighted networks model

BACKGROUND

Tuberculosis is known as one of the principal health problems, especially in developing countries. This study aimed to visualize, statistically model, and describe the weighted networks to investigate the intensity of social contacts associated with tuberculosis.

METHODS

In this case-control study, we applied weighted network analysis to assess the network of person-time spent in stores, workplaces, restaurants, mosques, Police bases, homes, hospitals, colleges, hairdressers, schools, contact homes, health centers, cinemas, parks, and markets. Modules will be determined based on the similarities between the variables in a topology overlap matrix. The most important variables will be found considering the association between each variable and module eigenvalues.

RESULTS

The result shows the extracted modules of locations based on the connectivity followed by the person-time at each place. The correlation (p-value) between TB and the turquoise, blue, and brown modules was 0.058 (0.351), 0.004 (0.943), and 0.117 (0.039), respectively. The brown module is the most important one, demonstrating a significant connection between homes, contact homes, health centers, and hospitals. Therefore, an association was found between person-time in four places and the occurrence of TB.

CONCLUSION

The finding of this study showed that most transmission of tuberculosis infection occurs in homes, contact homes, health centers, and hospitals. These place evaluations allow the identification of people with more contact and in need of screening, so critically leading to the identification of more patients with active TB.

Evolutionary history and spread of the Mycobacterium tuberculosis Latin American and Mediterranean (L4.3/LAM) sublineage, Tunisia

BACKGROUND

To infer the origin and spread of the Mycobacterium tuberculosis Latin American and Mediterranean (L4.3/LAM) sublineage in a Mediterranean country, Tunisia, where it predominates.

METHODS

We combined Bayesian (STRUCTURE) and maximum likelihood (MIGRAINE) estimation approaches based on a global 24-loci mycobacterial interspersed repetitive units-variable numbers of tandem repeats (MIRU-VNTR24) genotyping dataset consisting of 1573 L4.3/LAM clinical strains from four continents, including 252 isolates originating from Tunisia.

RESULTS

Phylogenetic analyses coupled with Bayesian estimations suggested that the most predominant L4.3/LAM subpopulation in Tunisia (65.07%), which is dominated by a single clonal complex, TUN4.3_CC1 (94.51%), has evolved from an ancestral pool that is restricted to Europe and Africa, contrasting with the remaining L4.3/LAM subpopulations whose ancestry was traced all over the word. Maximum likelihood analysis revealed that TUN4.3_CC1 has been undergoing a demographic expansion since 131 years ago (CI95%: 90.7-205), thus explaining its preponderance relative to the second most predominant CC, TUN4.3_CC2, whose population was found under contraction.

CONCLUSIONS

The preponderance of L4.3/LAM in Tunisia stems from a 130-year expansion process of a locally evolved clone.

Genomic Analysis of Multidrug-Resistant Mycobacterium tuberculosis Strains From Patients in Kazakhstan

Tuberculosis (TB) is an infectious disease that remains an essential public health problem in many countries. Despite decreasing numbers of new cases worldwide, the incidence of antibiotic-resistant forms (multidrug resistant and extensively drug-resistant) of TB is increasing. Next-generation sequencing technologies provide a high-throughput approach to identify known and novel potential genetic variants that are associated with drug resistance in Mycobacterium tuberculosis (Mtb). There are limited reports and data related to whole-genome characteristics of drug-resistant Mtb strains circulating in Kazakhstan. Here, we report whole-genome sequencing and analysis results of eight multidrug-resistant strains collected from TB patients in Kazakhstan. Genotyping and validation of all strains by MIRU-VNTR and spoligotyping methodologies revealed that these strains belong to the Beijing family. The spectrum of specific and potentially novel genomic variants (single-nucleotide polymorphisms, insertions, and deletions) related to drug resistance was identified and annotated. ResFinder, CARD, and CASTB antibiotic resistance databases were used for the characterization of genetic variants in genes associated with drug resistance. Our results provide reference data and genomic profiles of multidrug-resistant isolates for further comparative studies and investigations of genetic patterns in drug-resistant Mtb strains.

Characterization of Mutations Associated with Streptomycin Resistance in Multidrug-Resistant Mycobacterium tuberculosis in Zambia

Streptomycin (STR) is recommended for the management of multidrug-resistant tuberculosis (MDR-TB). Streptomycin resistance-conferring mutation types and frequency are shown to be influenced by genotypes of circulating strains in a population. This study aimed to characterize the mutations in MDR-TB isolates and examine their relationship with the genotypes in Zambia. A total of 138 MDR-TB isolates stored at the University Teaching Hospital Tuberculosis Reference Laboratory in Zambia were analyzed using spoligotyping and sequencing of STR resistance-associated genes. Streptomycin resistance was observed in 65.9% (91/138) of MDR-TB isolates. Mutations in rpsL, rrs, and gidB accounted for 33%, 12.1%, and 49.5%, respectively. Amino acid substitution K43R in rpsL was strongly associated with the CAS1_Kili genotype (p < 0.0001). The combination of three genes could predict 91.2% of STR resistance. Clustering of isolates based on resistance-conferring mutations and spoligotyping was observed. The clustering of isolates suggests that the increase in STR-resistant MDR-TB in Zambia is largely due to the spread of resistant strains from inadequate treatment. Therefore, rapid detection of STR resistance genetically is recommended before its use in MDR-TB treatment in Zambia.

Highly transmitted M. tuberculosis strains are more likely to evolve MDR/XDR and cause outbreaks, but what makes them highly transmitted?

Multi-Drug-Resistant strains of Mycobacterium tuberculosis (MDR-TB) are a serious obstacle to global TB eradication. While most MDR-TB strains are infrequently transmitted, a few cause large transmission clusters that contribute substantially to local MDR-TB burdens. Here we examine whether the known mutations in these strains can explain their success. Drug resistance mutations differ in fitness costs and strains can also acquire compensatory mutations (CM) to restore fitness, but some highly transmitted MDR strains have no CM. The acquisition of resistance mutations that maintain high transmissibility seems to occur by chance and are more likely in strains that are intrinsically highly transmitted and cause many cases. Modern Beijing lineage strains have caused several large outbreaks, but MDR outbreaks are also caused by ancient Beijing and lineage 4 strains, suggesting the lineage is less important than the characteristics of the individual strain. The development of fluoroquinolone resistance appears to represent another level of selection, in which strains must surmount unknown fitness costs of gyrA mutations. The genetic determinants of high transmission are poorly defined but may involve genes encoding proteins involved in molybdenum acquisition and the Esx systems. In addition, strains eliciting lower cytokine responses and producing more caseating granulomas may have advantages for transmission. Successful MDR/XDR strains generally evolve from highly transmitted drug sensitive parent strains due to selection pressures from deficiencies in local TB control programs. Until TB incidence is considerably reduced, there will likely be highly transmitted strains that develop resistance to any new antibiotic.

Low-Level Rifampin Resistance and rpoB Mutations in Mycobacterium tuberculosis: an Analysis of Whole-Genome Sequencing and Drug Susceptibility Test Data in New York

Rapid and reliable detection of rifampin (RIF) resistance is critical for the diagnosis and treatment of drug-resistant and multidrug-resistant (MDR) tuberculosis. Discordant RIF phenotype/genotype susceptibility results remain a challenge due to the presence of rpoB mutations that do not confer high levels of RIF resistance, as have been exhibited in strains with mutations such as Ser450Leu. These strains, termed low-level RIF resistant, exhibit elevated RIF MICs compared to fully susceptible strains but remain phenotypically susceptible by mycobacterial growth indicator tube (MGIT) testing and have been associated with poor patient outcomes. Here, we assess RIF resistance prediction by whole-genome sequencing (WGS) among a set of 1,779 prospectively tested strains by both prevalence of rpoB gene mutation and phenotype as part of routine clinical testing during a 2.5-year period. During this time, 139 strains were found to have nonsynonymous rpoB mutations, 53 of which were associated with RIF resistance, including both low-level and high-level resistance. Resistance to RIF (1.0 μg/ml in MGIT) was identified in 43 (81.1%) isolates. The remaining 10 (18.9%) strains were susceptible by MGIT but were confirmed to be low-level RIF resistant by MIC testing. Full rpoB gene sequencing overcame the limitations of critical concentration phenotyping, probe-based genotyping, and partial gene sequencing methods. Universal clinical WGS with concurrent phenotypic testing provided a more complete understanding of the prevalence and type of rpoB mutations and their association with RIF resistance in New York.

Reporting practices for genomic epidemiology of tuberculosis: a systematic review of the literature using STROME-ID guidelines as a benchmark

BACKGROUND

Pathogen genomics have become increasingly important in infectious disease epidemiology and public health. The Strengthening the Reporting of Molecular Epidemiology for Infectious Diseases (STROME-ID) guidelines were developed to outline a minimum set of criteria that should be reported in genomic epidemiology studies to facilitate assessment of study quality. We evaluate such reporting practices, using tuberculosis as an example.

METHODS

For this systematic review, we initially searched MEDLINE, Embase Classic, and Embase on May 3, 2017, using the search terms "tuberculosis" and "genom* sequencing". We updated this initial search on April 23, 2019, and also included a search of bioRxiv at this time. We included studies in English, French, or Spanish that recruited patients with microbiologically confirmed tuberculosis and used whole genome sequencing for typing of strains. Non-human studies, conference abstracts, and literature reviews were excluded. For each included study, the number and proportion of fulfilled STROME-ID criteria were recorded by two reviewers. A comparison of the mean proportion of fulfilled STROME-ID criteria before and after publication of the STROME-ID guidelines (in 2014) was done using a two-tailed t test. Quasi-Poisson regression and tobit regression were used to examine associations between study characteristics and the number and proportion of fulfilled STROME-ID criteria. This study was registered with PROSPERO, CRD42017064395.

FINDINGS

976 titles and abstracts were identified by our primary search, with an additional 16 studies identified in bioRxiv. 114 full texts (published between 2009 and 2019) were eligible for inclusion. The mean proportion of STROME-ID criteria fulfilled was 50% (SD 12; range 16-75). The proportion of criteria fulfilled was similar before and after STROME-ID publication (51% [SD 11] vs 46% [14], p=0·26). The number of criteria reported (among those applicable to all studies) was not associated with impact factor, h-index, country of affiliation of senior author, or sample size of isolates. Similarly, the proportion of criteria fulfilled was not associated with these characteristics, with the exception of a sample size of isolates of 277 or more (the highest quartile). In terms of reproducibility, 100 (88%) studies reported which bioinformatic tools were used, but only 33 (33%) reported corresponding version numbers. Sequencing data were available for 86 (75%) studies.

INTERPRETATION

The reporting of STROME-ID criteria in genomic epidemiology studies of tuberculosis between 2009 and 2019 was low, with implications for assessment of study quality. The considerable proportion of studies without bioinformatics version numbers or sequencing data available highlights a key concern for reproducibility.

Genetic Diversity in Mycobacterium tuberculosis Clinical Isolates and Resulting Outcomes of Tuberculosis Infection and Disease

Tuberculosis claims more human lives than any other bacterial infectious disease and represents a clear and present danger to global health as new tools for vaccination, treatment, and interruption of transmission have been slow to emerge. Additionally, tuberculosis presents with notable clinical heterogeneity, which complicates diagnosis, treatment, and the establishment of nonrelapsing cure. How this heterogeneity is driven by the diversity ofclinical isolates of the causative agent, Mycobacterium tuberculosis, has recently garnered attention. Herein, we review advances in the understanding of how naturally occurring variation in clinical isolates affects transmissibility, pathogenesis, immune modulation, and drug resistance. We also summarize how specific changes in transcriptional responses can modulate infection or disease outcome, together with strain-specific effects on gene essentiality. Further understanding of how this diversity of M. tuberculosis isolates affects disease and treatment outcomes will enable the development of more effective therapeutic options and vaccines for this dreaded disease.

Circulation of M. tuberculosis Beijing genotype in Latin America and the Caribbean

Lineage 2 (East Asian), which includes the Beijing genotype, is one of the most prevalent lineages of Mycobacterium tuberculosis (Mtb) throughout the world. The Beijing family is associated to hypervirulence and drug-resistant tuberculosis. The study of this genotype's circulation in Latin America is crucial for achieving total control of TB, the goal established by the World Health Organization, for the American sub-continent, before 2035. In this sense, the present work presents an overview of the status of the Beijing genotype for this region, with a bibliographical review, and data analysis of MIRU-VNTRs for available Beijing isolates. Certain countries present a prevalent trend of <5%, suggesting low transmissibility for the region, with the exception of Cuba (17.2%), Perú (16%) and Colombia (5%). Minimum Spanning Tree analysis, obtained from MIRU-VNTR data, shows distribution of specific clonal complex strains in each country. From this data, in most countries, we found that molecular epidemiology has not been a tool used for the control of TB, suggesting that the Beijing genotype may be underestimated in Latin America. It is recommended that countries with the highest incidence of the Beijing genotype use effective control strategies and increased care, as a requirement for public health systems.

Macro-geographical specificities of the prevailing tuberculosis epidemic as seen through SITVIT2, an updated version of the Mycobacterium tuberculosis genotyping database

In order to provide a global overview of genotypic, epidemiologic, demographic, phylogeographical, and drug resistance characteristics related to the prevailing tuberculosis (TB) epidemic, we hereby report an update of the 6th version of the international genotyping database SITVIT2. We also make all the available information accessible through a dedicated website (available at http://www.pasteur-guadeloupe.fr:8081/SITVIT2). Thanks to the public release of SITVIT2 which is currently the largest international multimarker genotyping database with a compilation of 111,635 clinical isolates from 169 countries of patient origin (131 countries of isolation, representing 1032 cities), our major aim is to highlight macro- and micro-geographical cleavages and phylogeographical specificities of circulating Mycobacterium tuberculosis complex (MTBC) clones worldwide. For this purpose, we retained strains typed by the most commonly used PCR-based methodology for TB genotyping, i.e., spoligotyping based on the polymorphism of the direct repeat (DR) locus, 5-loci Exact Tandem Repeats (ETRs), and MIRU-VNTR minisatellites used in 12-, 15-, or 24-loci formats. We describe the SITVIT2 database and integrated online applications that permit to interrogate the database using easy drop-down menus to draw maps, graphics and tables versus a long list of parameters and variables available for individual clinical isolates (year and place of isolation, origin, sex, and age of patient, drug-resistance, etc.). Available tools further allow to generate phylogenetical snapshot of circulating strains as Lineage-specific WebLogos, as well as minimum spanning trees of their genotypes in conjunction with their geographical distribution, drug-resistance, demographic, and epidemiologic characteristics instantaneously; whereas online statistical analyses let a user to pinpoint phylogeographical specificities of circulating MTBC lineages and conclude on actual demographic trends. Available associated information on gender (n = 18,944), age (n = 16,968), drug resistance (n = 19,606), and HIV serology (n = 2673), allowed to draw some important conclusions on TB geo-epidemiology; e.g. a positive correlation exists between certain Mycobacterium tuberculosis lineages (such as CAS and Beijing) and drug resistance (p-value<.001), while other lineages (such as LAM, X, and BOV) are more frequently associated with HIV-positive serology (p-value<.001). Besides, availability of information on the year of isolation of strains (range 1759-2012), also allowed to make tentative correlations between drug resistance information and lineages - portraying probable evolution trends over time and space. To conclude, the present approach of geographical mapping of predominant clinical isolates of tubercle bacilli causing the bulk of the disease both at country and regional level in conjunction with epidemiologic and demographic characteristics allows to shed new light on TB geo-epidemiology in relation with the continued waves of peopling and human migration.

Global expansion of Mycobacterium tuberculosis lineage 4 shaped by colonial migration and local adaptation

On the basis of population genomic and phylogeographic analyses of 1669 Mycobacterium tuberculosis lineage 4 (L4) genomes, we find that dispersal of L4 has been completely dominated by historical migrations out of Europe. We demonstrate an intimate temporal relationship between European colonial expansion into Africa and the Americas and the spread of L4 tuberculosis (TB). Markedly, in the age of antibiotics, mutations conferring antimicrobial resistance overwhelmingly emerged locally (at the level of nations), with minimal cross-border transmission of resistance. The latter finding was found to reflect the relatively recent emergence of these mutations, as a similar degree of local restriction was observed for susceptible variants emerging on comparable time scales. The restricted international transmission of drug-resistant TB suggests that containment efforts at the level of individual countries could be successful.

Importance of differential identification of Mycobacterium tuberculosis strains for understanding differences in their prevalence, treatment efficacy, and vaccine development

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb), remains a serious global health problem in the 21^st century because of its high mortality. Mtb is an extremely successful human-adapted pathogen that displays a multifactorial ability to control the host immune response and to evade killing by drugs, resulting in the breakdown of BCG vaccine-conferred anti-TB immunity and development of multidrug-resistant (MDR) and extensively drug-resistant (XDR) Mtb. Although genetic components of the genomes of the Mtb complex strains are highly conserved, showing over 99% similarity to other bacterial genera, recently accumulated evidence suggests that the genetic diversity of the Mtb complex strains has implications for treatment outcomes, development of MDR/XDR Mtb, BCG vaccine efficacy, transmissibility, and epidemiological outbreaks. Thus, new insights into the pathophysiological features of the Mtb complex strains are required for development of novel vaccines and for control of MDR/XDR Mtb infection, eventually leading to refinement of treatment regimens and the health care system. Many studies have focused on the differential identification of Mtb complex strains belonging to different lineages because of differences in their virulence and geographical dominance. In this review, we discuss the impact of differing genetic characteristics among Mtb complex strains on vaccine efficacy, treatment outcome, development of MDR/XDR Mtb strains, and epidemiological outbreaks by focusing on the best-adapted human Mtb lineages. We further explore the rationale for differential identification of Mtb strains for more effective control of TB in clinical and laboratory settings by scrutinizing current diagnostic methods.

RDRioMycobacterium tuberculosis lineage in the Brazil/Paraguay/Argentina triple border

The high tuberculosis (TB) incidence rates, the closeness of the cities and the high migration flux on the Brazil/Paraguay/Argentina border deserves an in-depth study, using Mycobacterial Interspersed Repetitive Unit (MIRU) and Spoligotyping genetic markers to explore the impact of the Mycobacterium tuberculosis RD^Rio lineage on disease transmission and resistance to anti-TB drugs in this setting. Although without the totality of M. tuberculosis isolates causing TB in this studied setting, a number of 97 isolates obtained from sputa samples culture of patients with confirmed TB, from 2013 to 2015, were submitted to 24 loci MIRU, Spoligotyping, detection of RD^Rio lineage and detection of mutation related to isoniazid and rifampicin resistance by MTBDRplus/DNA STRIP. In this sample, it was observed high clonal variability of circulating M. tuberculosis isolates causing TB in Brazilian cities bordering Paraguay and Argentina. The percentage of RD^Rio lineage causing TB in this setting was 15.46%, and lower than the detected in different areas of Brazil. According to 24 loci MIRU, the major MIRU International Type (MIT) related with RD^Rio lineage were MIT 26, MIT 738, MIT 601 with four, two and one isolates, respectively. Eight isolates with RD^Rio marker were classified as orphans. The mainly Spoligofamily related with RD^Rio lineage was LAM1 and LAM9 and no relationship between RD^Rio lineage and resistance in M. tuberculosis isolates circulating in this setting could be established. This work is pioneer in studying the dynamics of RD^Rio lineage transmission on the Brazil/Paraguay/Argentina border and deserves further studies to analyze the real contribution of the RD^Rio lineage in outbreaks and the risk of significant development of MDR-TB in the setting studied.

Antimicrobial resistance in Mycobacterium tuberculosis: mechanistic and evolutionary perspectives

Antibiotic-resistant Mycobacterium tuberculosis strains are threatening progress in containing the global tuberculosis epidemic. Mycobacterium tuberculosis is intrinsically resistant to many antibiotics, limiting the number of compounds available for treatment. This intrinsic resistance is due to a number of mechanisms including a thick, waxy, hydrophobic cell envelope and the presence of drug degrading and modifying enzymes. Resistance to the drugs which are active against M. tuberculosis is, in the absence of horizontally transferred resistance determinants, conferred by chromosomal mutations. These chromosomal mutations may confer drug resistance via modification or overexpression of the drug target, as well as by prevention of prodrug activation. Drug resistance mutations may have pleiotropic effects leading to a reduction in the bacterium's fitness, quantifiable e.g. by a reduction in the in vitro growth rate. Secondary so-called compensatory mutations, not involved in conferring resistance, can ameliorate the fitness cost by interacting epistatically with the resistance mutation. Although the genetic diversity of M. tuberculosis is low compared to other pathogenic bacteria, the strain genetic background has been demonstrated to influence multiple aspects in the evolution of drug resistance. The rate of resistance evolution and the fitness costs of drug resistance mutations may vary as a function of the genetic background.

A Nonsynonymous SNP Catalog of Mycobacterium tuberculosis Virulence Genes and Its Use for Detecting New Potentially Virulent Sublineages

Mycobacterium tuberculosis is divided into several distinct lineages, and various genetic markers such as IS-elements, VNTR, and SNPs are used for lineage identification. We propose an M. tuberculosis classification approach based on functional polymorphisms in virulence genes. An M. tuberculosis virulence genes catalog has been established, including 319 genes from various protein groups, such as proteases, cell wall proteins, fatty acid and lipid metabolism proteins, sigma factors, toxin–antitoxin systems. Another catalog of 1,573 M. tuberculosis isolates of different lineages has been developed. The developed SNP-calling program has identified 3,563 nonsynonymous SNPs. The constructed SNP-based phylogeny reflected the evolutionary relationship between lineages and detected new sublineages. SNP analysis of sublineage F15/LAM4/KZN revealed four lineage-specific mutations in cyp125, mce3B, vapC25, and vapB34. The Ural lineage has been divided into two geographical clusters based on different SNPs in virulence genes. A new sublineage, B0/N-90, was detected inside the Beijing-B0/W-148 by SNPs in irtB, mce3F and vapC46. We have found 27 members of B0/N-90 among the 227 available genomes of the Beijing-B0/W-148 sublineage. Whole-genome sequencing of strain B9741, isolated from an HIV-positive patient, was demonstrated to belong to the new B0/N-90 group. A primer set for PCR detection of B0/N-90 lineage-specific mutations has been developed. The prospective use of mce3 mutant genes as genetically engineered vaccine is discussed.

Characterization of EccA3, a CbbX family ATPase from the ESX-3 secretion pathway of M. tuberculosis

EccA family proteins are conserved components of ESX secretion pathways in M. tuberculosis H37Rv. Here, we report the characterization of EccA3 (Rv0282), a CbbX family AAA (ATPases Associated with diverse cellular Activities) protein from the ESX-3 pathway that is required for in vitro growth of mycobacteria, secretion of virulence factors, and acquisition of iron and zinc. EccA3 is a thermostable ATPase with a molecular weight of ~68kDa. It exists as a dodecamer in the apo form and associates as a hexamer in the presence of ATP. Its C-terminal region consists of a CbbX-like AAA-domain while the N-terminal region contains a tetratricopeptide repeat (TPR) domain with lower homology to other EccA-type proteins. Further, the C-terminal domain functions as the oligomerization domain and also exhibits ATPase activity. Mutational analysis, steady state kinetics and molecular docking studies identify R573 as the important 'sensor arginine' and R505 as an 'arginine finger' in EccA3. Dynamic fluorescence quenching experiments suggest that the N-terminal domain moves closer to the C-terminal domain upon ATP-binding. The ATP-dependent 'open-close' relative movements of the two domains might help EccA3 interaction and secretion of essential virulence factors.

Strain Diversity and the Evolution of Antibiotic Resistance

Drug resistance is best thought of as an ongoing biological process. Resistant bacteria must emerge, become established and ultimately transmit in order to be relevant to human health. In this context, genetic diversity can influence the rate and likelihood of resistance emerging; it can also modulate the net physiological impact of resistance and the propensity of an organism to improve any defects that arise from it. Combined, these effects can have an impact on a larger scale, with highly transmissible drug-resistant bacterial strains posing a formidable threat to global health. These considerations are pertinent to the future of tuberculosis control as well. In this chapter, we review our current understanding of the impact of genetic diversity in the broadest sense on the evolution of drug-resistant members of the Mycobacterium tuberculosis complex.

Mycobacterium tuberculosis lineage 4 comprises globally distributed and geographically restricted sublineages

Generalist and specialist species differ in the breadth of their ecological niches. Little is known about the niche width of obligate human pathogens. Here we analyzed a global collection of Mycobacterium tuberculosis lineage 4 clinical isolates, the most geographically widespread cause of human tuberculosis. We show that lineage 4 comprises globally distributed and geographically restricted sublineages, suggesting a distinction between generalists and specialists. Population genomic analyses showed that, whereas the majority of human T cell epitopes were conserved in all sublineages, the proportion of variable epitopes was higher in generalists. Our data further support a European origin for the most common generalist sublineage. Hence, the global success of lineage 4 reflects distinct strategies adopted by different sublineages and the influence of human migration.

Whole genome sequencing for M/XDR tuberculosis surveillance and for resistance testing

Whole genome sequencing (WGS) can help to relate Mycobacterium tuberculosis genomes to one another to assess genetic relatedness and infer the likelihood of transmission between cases. The same sequence data are now increasingly being used to predict drug resistance and susceptibility. Controlling the spread of tuberculosis and providing patients with the correct treatment are central to the World Health Organization's target to 'End TB' by 2035, for which the global prevalence of drug-resistant tuberculosis remains one of the main obstacles to success. So far, WGS has been applied largely to drug-susceptible strains for the purposes of understanding transmission, leaving a number of analytical considerations before transferring what has been learnt from drug-susceptible disease to drug-resistant tuberculosis. We discuss these potential problems here, alongside some of the challenges to characterizing the Mycobacterium tuberculosis 'resistome'-the optimal knowledge-base required for WGS-based assays to successfully direct individualized treatment regimens through the prediction of drug resistance and susceptibility in the future.

Whole genome sequencing of Mycobacterium tuberculosis for detection of drug resistance: a systematic review

OBJECTIVES

We conducted a systematic review to determine the diagnostic accuracy of whole genome sequencing (WGS) of Mycobacterium tuberculosis for the detection of resistance to first- and second-line anti-tuberculosis (TB) drugs.

METHODS

The study was conducted according to the criteria of the Preferred Reporting Items for Systematic Reviews group. A total of 20 publications were included. The sensitivity, specificity, positive-predictive value and negative-predictive value of WGS using phenotypic drug susceptibility testing methods as a reference standard were determined.

RESULTS

Anti-TB agents tested included all first-line drugs, a variety of reserve drugs, as well as new drugs. Polymorphisms in a total of 53 genes were tested for associations with drug resistance. Pooled sensitivity and specificity values for detection of resistance to selected first-line drugs were 0.98 (95% CI 0.93-0.98) and 0.98 (95% CI 0.98-1.00) for rifampicin and 0.97 (95% CI 0.94-0.99) and 0.93 (95% CI 0.91-0.96) for isoniazid, respectively. Due to high heterogeneity in study designs, lack of data, knowledge of resistance mechanisms and clarity on exclusion of phylogenetic markers, there was a significant variation in analytical performance of WGS for the remaining first-line, reserved drugs and new drugs.

CONCLUSIONS

Whole genome sequencing could be considered a promising alternative to existing phenotypic and molecular drug susceptibility testing methods for rifampicin and isoniazid pending standardization of analytical pipelines. To ensure clinical relevance of WGS for detection of M. tuberculosis complex drug resistance, future studies should include information on clinical outcomes.

First Evaluation of GenoType MTBDRplus 2.0 Performed Directly on Respiratory Specimens in Central America

The turnaround times for conventional methods used to detect Mycobacterium tuberculosis in sputum samples and to obtain drug susceptibility information are long in many developing countries, including Panama, leading to delays in appropriate treatment initiation and continued transmission in the community. We evaluated the performance of a molecular line probe assay, the Genotype MTBDRplus version 2.0 assay, in detecting M. tuberculosis complex directly in respiratory specimens from smear-positive tuberculosis cases from four different regions in Panama, as well as the most frequent mutations in genes conferring resistance to isoniazid (katG and inhA) and rifampin (rpoB). Our results were confirmed with the nitrate reductase assay and genomic sequencing. M. tuberculosis complex was detected by the Genotype MTBDRplus 2.0 assay with 100% sensitivity and specificity. The sensitivity and specificity for rifampin resistance were 100% and 100%, respectively, and those for isoniazid resistance were 90.7% and 100%. Isoniazid monoresistance was detected in 5.2% of new cases. Genotype MTBDRplus 2.0 is highly accurate in detecting M. tuberculosis complex in respiratory specimens and is able to discriminate isoniazid-monoresistant cases from multidrug-resistant cases within 2 days.

The Distribution of Fitness Costs of Resistance-Conferring Mutations Is a Key Determinant for the Future Burden of Drug-Resistant Tuberculosis: A Model-Based Analysis

BACKGROUND

Drug resistance poses a serious challenge for the control of tuberculosis in many settings. It is well established that the expected future trend in resistance depends on the reproductive fitness of drug-resistant Mycobacterium tuberculosis. However, the variability in fitness between strains with different resistance-conferring mutations has been largely ignored when making these predictions.

METHODS

We developed a novel approach for incorporating the variable fitness costs of drug resistance-conferring mutations and for tracking this distribution of fitness costs over time within a transmission model. We used this approach to describe the effects of realistic fitness cost distributions on the future prevalence of drug-resistant tuberculosis.

RESULTS

The shape of the distribution of fitness costs was a strong predictor of the long-term prevalence of resistance. While, as expected, lower average fitness costs of drug resistance-conferring mutations were associated with more severe epidemics of drug-resistant tuberculosis, fitness distributions with greater variance also led to higher levels of drug resistance. For example, compared to simulations in which the fitness cost of resistance was fixed, introducing a realistic amount of variance resulted in a 40% increase in prevalence of drug-resistant tuberculosis after 20 years.

CONCLUSIONS

The differences in the fitness costs associated with drug resistance-conferring mutations are a key determinant of the future burden of drug-resistant tuberculosis. Future studies that can better establish the range of fitness costs associated with drug resistance-conferring mutations will improve projections and thus facilitate better public health planning efforts.

Interpreting whole genome sequencing for investigating tuberculosis transmission: a systematic review

BACKGROUND

Whole genome sequencing (WGS) is becoming an important part of epidemiological investigations of infectious diseases due to greater resolution and cost reductions compared to traditional typing approaches. Many public health and clinical teams will increasingly use WGS to investigate clusters of potential pathogen transmission, making it crucial to understand the benefits and assumptions of the analytical methods for investigating the data. We aimed to understand how different approaches affect inferences of transmission dynamics and outline limitations of the methods.

METHODS

We comprehensively searched electronic databases for studies that presented methods used to interpret WGS data for investigating tuberculosis (TB) transmission. Two authors independently selected studies for inclusion and extracted data. Due to considerable methodological heterogeneity between studies, we present summary data with accompanying narrative synthesis rather than pooled analyses.

RESULTS

Twenty-five studies met our inclusion criteria. Despite the range of interpretation tools, the usefulness of WGS data in understanding TB transmission often depends on the amount of genetic diversity in the setting. Where diversity is small, distinguishing re-infections from relapses may be impossible; interpretation may be aided by the use of epidemiological data, examining minor variants and deep sequencing. Conversely, when within-host diversity is large, due to genetic hitchhiking or co-infection of two dissimilar strains, it is critical to understand how it arose. Greater understanding of microevolution and mixed infection will enhance interpretation of WGS data.

CONCLUSIONS

As sequencing studies have sampled more intensely and integrated multiple sources of information, the understanding of TB transmission and diversity has grown, but there is still much to be learnt about the origins of diversity that will affect inferences from these data. Public health teams and researchers should combine epidemiological, clinical and WGS data to strengthen investigations of transmission.

Latin-American-Mediterranean lineage of Mycobacterium tuberculosis: Human traces across pathogen's phylogeography

Currently, Mycobacterium tuberculosis isolates of Latin-American Mediterranean (LAM) family may be detected far beyond the geographic areas that coined its name 15years ago. Here, we established the framework phylogeny of this geographically intriguing and pathobiologically important mycobacterial lineage and hypothesized how human demographics and migration influenced its phylogeography. Phylogenetic analysis of LAM isolates from all continents based on 24 variable number of tandem repeats (VNTR) loci and other markers identified three global sublineages with certain geographic affinities and defined by large deletions RD115, RD174, and by spoligotype SIT33. One minor sublineage (spoligotype SIT388) appears endemic in Japan. One-locus VNTR signatures were established for sublineages and served for their search in published literature and geographic mapping. We suggest that the LAM family originated in the Western Mediterranean region. The most widespread RD115 sublineage seems the most ancient and encompasses genetically and geographically distant branches, including extremely drug resistant KZN in South Africa and LAM-RUS recently widespread across Northern Eurasia. The RD174 sublineage likely started its active spread in Brazil; its earlier branch is relatively dominated by isolates from South America and the derived one is dominated by Portuguese and South/Southeastern African isolates. The relatively most recent SIT33-sublineage is marked with enigmatic gaps and peaks across the Americas and includes South African clade F11/RD761, which likely emerged within the SIT33 subpopulation after its arrival to Africa. In addition to SIT388-sublineage, other deeply rooted, endemic LAM sublineages may exist that remain to be discovered. As a general conclusion, human mass migration appears to be the major factor that shaped the M. tuberculosis phylogeography over large time-spans.

Whole genome sequencing to complement tuberculosis drug resistance surveys in Uganda

Understanding the circulating Mycobacterium tuberculosis resistance mutations is vital for better TB control strategies, especially to inform a new MDR-TB treatment programme. We complemented the phenotypic drug susceptibility testing (DST) based drug resistance surveys (DRSs) conducted in Uganda between 2008 and 2011 with Whole Genome Sequencing (WGS) of 90 Mycobacterium tuberculosis isolates phenotypically resistant to rifampicin and/or isoniazid to better understand the extent of drug resistance. A total of 31 (34.4 %) patients had MDR-TB, 5 (5.6 %) mono-rifampicin resistance and 54 (60.0 %) mono-isoniazid resistance by phenotypic DST. Pyrazinamide resistance mutations were identified in 32.3% of the MDR-TB patients. Resistance to injectable agents was detected in 4/90 (4.4%), and none to fluoroquinolones or novel drugs. Compensatory mutations in rpoC were identified in two patients. The sensitivity and specificity of drug resistance mutations compared to phenotypic DST were for rpoB 88.6% and 98.1%, katG 60.0% and 100%, fabG1 16.5% and 100%, katG and/or fabG1 71.8% and 100%, embCAB 63.0% and 82.5%, rrs 11.4% and 100%, rpsL 20.5% and 95.7% and rrs and/or rpsL 31.8% and 95.7%. Phylogenetic analysis showed dispersed MDR-TB isolate, with only one cluster of three Beijing family from South West Uganda. Among tuberculosis patients in Uganda, resistance beyond first-line drugs as well as compensatory mutations remain low, and MDR-TB isolates did not arise from a dominant clone. Our findings show the potential use of sequencing for complementing DRSs or surveillance in this setting, with good specificity compared to phenotypic DST. The reported high confidence mutations can be included in molecular assays, and population-based studies can track transmission of MDR-TB including the Beijing family strains in the South West of the country.

The complex evolution of antibiotic resistance in Mycobacterium tuberculosis

Multidrug-resistant and extensively drug-resistant tuberculosis (TB) represent a major threat to the control of the disease worldwide. The mechanisms and pathways that result in the emergence and subsequent fixation of resistant strains of Mycobacterium tuberculosis are not fully understood and recent studies suggest that they are much more complex than initially thought. In this review, we highlight the exciting new areas of research within TB resistance that are beginning to fill these gaps in our understanding, whilst also raising new questions and providing future directions.

Burden of transmitted multidrug resistance in epidemics of tuberculosis: a transmission modelling analysis

Background

Multidrug-resistant tuberculosis (MDR-TB) can be acquired through de novo mutation during TB treatment or through transmission from other individuals with active MDR-TB. Understanding the balance between these two mechanisms is essential when allocating resources for MDR-TB.

Methods

We constructed a dynamic transmission model of an MDR-TB epidemic, allowing for both treatment-related acquisition and person-to-person transmission of resistance. We used national TB notification data to inform Bayesian estimates of the fraction of each country’s 2013 MDR-TB incidence that resulted from MDR transmission rather than treatment-related MDR acquisition.

Findings

Global estimates of 3·5% MDR-TB prevalence among new TB notifications and 20·5% among retreatment notifications translate into an estimate that resistance transmission rather than acquisition accounts for a median 96% (95% UR: 68–100%) of all incident MDR-TB, and 61% (16–95%) of incident MDR-TB in previously-treated individuals. The estimated percentage of MDR-TB resulting from transmission varied substantially with different countries’ notification data; for example, we estimated this percentage at 48% (30–75%) of MDR-TB in Bangladesh, versus 99% (91–100%) in Uzbekistan. Estimates were most sensitive to estimates of the transmissibility of MDR strains, the probability of acquiring MDR during tuberculosis treatment, and the responsiveness of MDR TB to first-line treatment.

Interpretation

Notifications of MDR prevalence from most high-burden settings are most consistent with the vast majority of incident MDR-TB resulting from transmission rather than new treatment-related acquisition of resistance. Merely improving the treatment of drug-susceptible TB is unlikely to greatly reduce future MDR-TB incidence. Improved diagnosis and treatment of MDR-TB – including new tests and drug regimens – should be highly prioritized.

Clinical value of whole-genome sequencing of Mycobacterium tuberculosis

Whole-genome sequencing (WGS) is now common as a result of new technologies that can rapidly sequence a complete bacterial genome for US$500 or less. Many studies have addressed questions about tuberculosis with WGS, and knowing the sequence of the entire genome, rather than only a few fragments, has greatly increased the precision of molecular epidemiology and contact tracing. Additionally, topics such as the mutation rate, drug resistance, the target of new drugs, and the phylogeny and evolution of the Mycobacterium tuberculosis complex bacteria have been elucidated by WGS. Nonetheless, WGS has not explained differences in transmissibility between strains, or why some strains are more virulent than others or more prone to development of multidrug resistance. With advances in technology, WGS of clinical specimens could become routine in high-income countries; however, its relevance will probably depend on easy to use software to efficiently process the sequences produced and accessible genomic databases that can be mined in future studies.

Genome Analysis of the First Extensively Drug-Resistant (XDR) Mycobacterium tuberculosis in Malaysia Provides Insights into the Genetic Basis of Its Biology and Drug Resistance

The outbreak of extensively drug-resistant tuberculosis (XDR-TB) has become an increasing problem in many TB-burdened countries. The underlying drug resistance mechanisms, including the genetic variation favored by selective pressure in the resistant population, are partially understood. Recently, the first case of XDR-TB was reported in Malaysia. However, the detailed genotype family and mechanisms of the formation of multiple drugs resistance are unknown. We sequenced the whole genome of the UM 1072388579 strain with a 2-kb insert-size library and combined with that from previously sequenced 500-bp-insert paired-end reads to produce an improved sequence with maximal sequencing coverage across the genome. In silico spoligotyping and phylogenetic analyses demonstrated that UM 1072388579 strain belongs to an ancestral-like, non-Beijing clade of East Asia lineage. This is supported by the presence of a number of lineage-specific markers, including fadD28, embA, nuoD and pks7. Polymorphism analysis showed that the drug-susceptibility profile is correlated with the pattern of resistance mutations. Mutations in drug-efflux pumps and the cell wall biogenesis pathway such as mmpL, pks and fadD genes may play an important role in survival and adaptation of this strain to its surrounding environment. In this work, fifty-seven putative promoter SNPs were identified. Among them, we identified a novel SNP located at -4 T allele of TetR/acrR promoter as an informative marker to recognize strains of East Asian lineage. Our work indicates that the UM 1072388579 harbors both classical and uncommon SNPs that allow it to escape from inhibition by many antibiotics. This study provides a strong foundation to dissect the biology and underlying resistance mechanisms of the first reported XDR M. tuberculosis in Malaysia.

Map the gap: missing children with drug-resistant tuberculosis

BACKGROUND

The lack of published information about children with multidrug-resistant tuberculosis (MDR-TB) is an obstacle to efforts to advocate for better diagnostics and treatment.

OBJECTIVE

To describe the lack of recognition in the published literature of MDR-TB and extensively drug-resistant TB (XDR-TB) in children.

DESIGN

We conducted a systematic search of the literature published in countries that reported any MDR- or XDR-TB case by 2012 to identify MDR- or XDR-TB cases in adults and in children.

RESULTS

Of 184 countries and territories that reported any case of MDR-TB during 2005-2012, we identified adult MDR-TB cases in the published literature in 143 (78%) countries and pediatric MDR-TB cases in 78 (42%) countries. Of the 92 countries that reported any case of XDR-TB, we identified adult XDR-TB cases in the published literature in 55 (60%) countries and pediatric XDR-TB cases for 9 (10%) countries.

CONCLUSION

The absence of publications documenting child MDR- and XDR-TB cases in settings where MDR- and XDR-TB in adults have been reported indicates both exclusion of childhood disease from the public discourse on drug-resistant TB and likely underdetection of sick children. Our results highlight a large-scale lack of awareness about children with MDR- and XDR-TB.

Culture-independent detection and characterisation of Mycobacterium tuberculosis and M. africanum in sputum samples using shotgun metagenomics on a benchtop sequencer

Tuberculosis remains a major global health problem. Laboratory diagnostic methods that allow effective, early detection of cases are central to management of tuberculosis in the individual patient and in the community. Since the 1880s, laboratory diagnosis of tuberculosis has relied primarily on microscopy and culture. However, microscopy fails to provide species- or lineage-level identification and culture-based workflows for diagnosis of tuberculosis remain complex, expensive, slow, technically demanding and poorly able to handle mixed infections. We therefore explored the potential of shotgun metagenomics, sequencing of DNA from samples without culture or target-specific amplification or capture, to detect and characterise strains from the Mycobacterium tuberculosis complex in smear-positive sputum samples obtained from The Gambia in West Africa. Eight smear- and culture-positive sputum samples were investigated using a differential-lysis protocol followed by a kit-based DNA extraction method, with sequencing performed on a benchtop sequencing instrument, the Illumina MiSeq. The number of sequence reads in each sputum-derived metagenome ranged from 989,442 to 2,818,238. The proportion of reads in each metagenome mapping against the human genome ranged from 20% to 99%. We were able to detect sequences from the M. tuberculosis complex in all eight samples, with coverage of the H37Rv reference genome ranging from 0.002X to 0.7X. By analysing the distribution of large sequence polymorphisms (deletions and the locations of the insertion element IS6110) and single nucleotide polymorphisms (SNPs), we were able to assign seven of eight metagenome-derived genomes to a species and lineage within the M. tuberculosis complex. Two metagenome-derived mycobacterial genomes were assigned to M. africanum, a species largely confined to West Africa; the others that could be assigned belonged to lineages T, H or LAM within the clade of "modern" M. tuberculosis strains. We have provided proof of principle that shotgun metagenomics can be used to detect and characterise M. tuberculosis sequences from sputum samples without culture or target-specific amplification or capture, using an accessible benchtop-sequencing platform, the Illumina MiSeq, and relatively simple DNA extraction, sequencing and bioinformatics protocols. In our hands, sputum metagenomics does not yet deliver sufficient depth of coverage to allow sequence-based sensitivity testing; it remains to be determined whether improvements in DNA extraction protocols alone can deliver this or whether culture, capture or amplification steps will be required. Nonetheless, we can foresee a tipping point when a unified automated metagenomics-based workflow might start to compete with the plethora of methods currently in use in the diagnostic microbiology laboratory.

Evolution of drug resistance in tuberculosis: recent progress and implications for diagnosis and therapy

Drug-resistant tuberculosis is a growing threat to global public health. Recent efforts to understand the evolution of drug resistance have shown that changes in drug-target interactions are only the first step in a longer adaptive process. The emergence of transmissible drug-resistant Mycobacterium tuberculosis is the result of a multitude of additional genetic mutations, many of which interact, a phenomenon known as epistasis. The varied effects of these epistatic interactions include compensating for the reduction of the biological cost associated with the development of drug resistance, increasing the level of resistance, and possibly accommodating broader changes in the physiology of resistant bacteria. Knowledge of these processes and our ability to detect them as they happen informs the development of diagnostic tools and better control strategies. In particular, the use of whole genome sequencing combined with surveillance efforts in the field could provide a powerful instrument to prevent future epidemics of drug-resistant tuberculosis.

Mycobacterium tuberculosis isolates from single outpatient clinic in Panama City exhibit wide genetic diversity

Understanding Mycobacterium tuberculosis biodiversity and transmission is significant for tuberculosis control. This short report aimed to determine the genetic diversity of M. tuberculosis isolates from an outpatient clinic in Panama City. A total of 62 M. tuberculosis isolates were genotyped by 12 loci mycobacterial interspersed repetitive units-variable number of tandem repeats (MIRU-VNTR) and Spoligotyping. Forty-five (72.6%) of the isolates showed unique MIRU-VNTR genotypes, and 13 (21%) of the isolates were grouped into four clusters. Four isolates showed polyclonal MIRU-VNTR genotypes. The MIRU-VNTR Hunter-Gaston discriminatory index reached 0.988. The Spoligotyping analysis revealed 16 M. tuberculosis families, including Latin American-Mediterranean, Harlem, and Beijing. These findings suggest a wide genetic diversity of M. tuberculosis isolates at one outpatient clinic. A detailed molecular epidemiology survey is now warranted, especially following second massive immigration for local Panama Canal expansion activities.

Mycobacterium tuberculosis Latin American-Mediterranean family and its sublineages in the light of robust evolutionary markers

Mycobacterium tuberculosis has a clonal population structure, and the Latin American-Mediterranean (LAM) family is one of the largest and most widespread within this species, showing evidence for remarkable pathobiology and a confusing phylogeny. Here, we applied robust phylogenetic markers to study the evolution of the LAM family and its major sublineages circulating in Russia and neighboring countries. A total of 250 M. tuberculosis isolates were confirmed to belong to the LAM family based on the analysis of the LAM-specific single-nucleotide polymorphisms (SNPs) in the Rv3062 and Rv0129c genes. At this stage, the family status was rectified for 121 isolates misleadingly assigned by CRISPR spoligotyping to non-LAM families (T1- or T5-RUS1). Consequently, the reestimated LAM prevalence rate increased 2-fold in Russia and Kazakhstan and 4-fold in Belarus. The majority (91.8 to 98.7%) of the LAM isolates from all three countries belonged to the LAM-RUS sublineage. In contrast, the Ibero-American LAM RD-Rio sublineage was identified in only 7 Russian isolates. Taken together, our findings and further analyses suggest a monophyletic origin of LAM-RUS: at a historically distant time, in Russia, in a small founding bacterial/human population. Its dissemination pattern and high prevalence rate in Northern Eurasia may indicate a long-term coexistence of the LAM-RUS sublineage and local human populations hypothetically leading to coadaptation and reduced pathogenicity of the relatively more ancient clones, such as spoligotype international type 254 (SIT254), compared to the more recent SIT252 and SIT266 clones. In contrast, rare LAM RD-Rio isolates were likely brought to Russia through occasional human contact. The spread of RD-Rio strains is not as global as commonly claimed and is determined largely by human migration flows (rather than by pathobiological properties of these strains). Consequently, a host population factor appears to play a major role in shaping the in situ dissemination pattern of the imported strains in an autochthonous population.