Using genomics to understand the origin and dispersion of multidrug and extensively drug resistant tuberculosis in Portugal

Delineating the evolutionary pathway to multidrug-resistant outbreaks of a Mycobacterium tuberculosis L4.1.2.1/Haarlem sublineage

OBJECTIVES

We sought to capture the evolutionary itinerary of the Mycobacterium tuberculosis L4.1.2.1/Haarlem sublineage in northern Tunisia, where it caused a major multidrug-resistant (MDR) tuberculosis outbreak in a context strictly negative for HIV infection.

METHODS

We combined whole genome sequencing and Bayesian approaches using a representative collection of drug-susceptible and drug-resistant L4.1.2.1/Haarlem clinical strains (n = 121) recovered from the outbreak region over 16 years.

RESULTS

In the absence of drug resistance, the L4.1.2.1/Haarlem sublineage showed a propensity for rapid transmission as witnessed by the high clustering (44.6%) and recent transmission rates (25%), as well as the reduced mean distance between genome pairs. The entire pool of L4.1.2.1/Haarlem MDR strains was found to be linked to either the aforementioned major outbreak (68 individuals, 2001-2016) or to a minor, newly uncovered outbreak (six cases, 2001-2011). Strikingly, the two outbreaks descended independently from a common ancestor that can be dated back to 1886.

CONCLUSIONS

Our data point to the intrinsic propensity for rapid transmission of the M. tuberculosis L4.1.2.1/Haarlem sublineage in northern Tunisia, linking the overall MDR tuberculosis epidemic to a single ancestor. These findings bring out the important role of the bacillus' genetic background in the emergence of successful MDR M. tuberculosis clones.

Genetic Characterization and Population Structure of Drug-Resistant Mycobacterium tuberculosis Isolated from Brazilian Patients Using Whole-Genome Sequencing

The present study aimed to determine the genetic diversity of isolates of Mycobacterium tuberculosis (Mtb) from presumed drug-resistant tuberculosis patients from several states of Brazil. The isolates had been submitted to conventional drug susceptibility testing for first- and second-line drugs. Multidrug-resistant (MDR-TB) (54.8%) was the most frequent phenotypic resistance profile, in addition to an important high frequency of pre-extensive resistance (p-XDR-TB) (9.2%). Using whole-genome sequencing (WGS), we characterized 298 Mtb isolates from Brazil. Besides the analysis of genotype distribution and possible correlations between molecular and clinical data, we determined the performance of an in-house WGS pipeline with other online pipelines for Mtb lineages and drug resistance profile definitions. Sub-lineage 4.3 (52%) was the most frequent genotype, and the genomic approach revealed a p-XDR-TB level of 22.5%. We detected twenty novel mutations in three resistance genes, and six of these were observed in eight phenotypically resistant isolates. A cluster analysis of 170 isolates showed that 43.5% of the TB patients belonged to 24 genomic clusters, suggesting considerable ongoing transmission of DR-TB, including two interstate transmissions. The in-house WGS pipeline showed the best overall performance in drug resistance prediction, presenting the best accuracy values for five of the nine drugs tested. Significant associations were observed between suffering from fatal disease and genotypic p-XDR-TB (p = 0.03) and either phenotypic (p = 0.006) or genotypic (p = 0.0007) ethambutol resistance. The use of WGS analysis improved our understanding of the population structure of MTBC in Brazil and the genetic and clinical data correlations and demonstrated its utility for surveillance efforts regarding the spread of DR-TB, hopefully helping to avoid the emergence of even more resistant strains and to reduce TB incidence and mortality rates.

Multidrug-resistant tuberculosis

Tuberculosis (TB) remains the foremost cause of death by an infectious disease globally. Multidrug-resistant or rifampicin-resistant TB (MDR/RR-TB; resistance to rifampicin and isoniazid, or rifampicin alone) is a burgeoning public health challenge in several parts of the world, and especially Eastern Europe, Russia, Asia and sub-Saharan Africa. Pre-extensively drug-resistant TB (pre-XDR-TB) refers to MDR/RR-TB that is also resistant to a fluoroquinolone, and extensively drug-resistant TB (XDR-TB) isolates are additionally resistant to other key drugs such as bedaquiline and/or linezolid. Collectively, these subgroups are referred to as drug-resistant TB (DR-TB). All forms of DR-TB can be as transmissible as rifampicin-susceptible TB; however, it is more difficult to diagnose, is associated with higher mortality and morbidity, and higher rates of post-TB lung damage. The various forms of DR-TB often consume >50% of national TB budgets despite comprising <5-10% of the total TB case-load. The past decade has seen a dramatic change in the DR-TB treatment landscape with the introduction of new diagnostics and therapeutic agents. However, there is limited guidance on understanding and managing various aspects of this complex entity, including the pathogenesis, transmission, diagnosis, management and prevention of MDR-TB and XDR-TB, especially at the primary care physician level.

simpiTB - a pipeline designed to extract meaningful information from whole genome sequencing data of Mycobacterium tuberculosis complex, allows to combine genomic, phylogenetic and clustering analyses in existing SITVIT databases

Data obtained from new sequencing technologies are evolving rapidly, leading to the development of specific bioinformatic tools, pipelines and softwares. Several algorithms and tools are today available allowing a better identification and description of Mycobacterium tuberculosis complex (MTBC) isolates worldwide. Our approach consists in applying existing methods to analyze DNA sequencing data (from FASTA or FASTQ files), and tentatively extract meaningful information that would facilitate identification as well as a better understanding and management of MTBC isolates (taking into account whole genome sequencing and classical genotyping data). The aim of this study is to propose a pipeline analysis allowing to potentially simplify MTBC data analysis by providing different ways to interpret genomic or genotyping information based on existing tools. Furthermore, we propose a "reconciledTB" list making a link with results directly obtained from whole genome sequencing (WGS) data and results obtained from classical genotyping analysis (data inferred from SpoTyping and MIRUReader). Data visualization graphics and trees generated provide additional elements to better understand and confer associations among information overlap analyses. Additionally, comparison between data entered in an international genotyping database (SITVITEXTEND) and ensuing data obtained from the pipeline not only provide meaningful information, but further suggest that simpiTB could also be suitable for new data integration in specific TB genotyping databases.

Drug resistance prediction for Mycobacterium tuberculosis with reference graphs

Tuberculosis is a global pandemic disease with a rising burden of antimicrobial resistance. As a result, the World Health Organization (WHO) has a goal of enabling universal access to drug susceptibility testing (DST). Given the slowness of and infrastructure requirements for phenotypic DST, whole-genome sequencing, followed by genotype-based prediction of DST, now provides a route to achieving this. Since a central component of genotypic DST is to detect the presence of any known resistance-causing mutations, a natural approach is to use a reference graph that allows encoding of known variation. We have developed DrPRG (Drug resistance Prediction with Reference Graphs) using the bacterial reference graph method Pandora. First, we outline the construction of a Mycobacterium tuberculosis drug resistance reference graph. The graph is built from a global dataset of isolates with varying drug susceptibility profiles, thus capturing common and rare resistance- and susceptible-associated haplotypes. We benchmark DrPRG against the existing graph-based tool Mykrobe and the haplotype-based approach of TBProfiler using 44 709 and 138 publicly available Illumina and Nanopore samples with associated phenotypes. We find that DrPRG has significantly improved sensitivity and specificity for some drugs compared to these tools, with no significant decreases. It uses significantly less computational memory than both tools, and provides significantly faster runtimes, except when runtime is compared to Mykrobe with Nanopore data. We discover and discuss novel insights into resistance-conferring variation for M. tuberculosis - including deletion of genes katG and pncA - and suggest mutations that may warrant reclassification as associated with resistance.

Transcriptional regulation and drug resistance in Mycobacterium tuberculosis

Bacterial drug resistance is one of the major challenges to present and future human health, as the continuous selection of multidrug resistant bacteria poses at serious risk the possibility to treat infectious diseases in the near future. One of the infection at higher risk to become incurable is tuberculosis, due to the few drugs available in the market against Mycobacterium tuberculosis. Drug resistance in this species is usually due to point mutations in the drug target or in proteins required to activate prodrugs. However, another interesting and underexplored aspect of bacterial physiology with important impact on drug susceptibility is represented by the changes in transcriptional regulation following drug exposure. The main regulators involved in this phenomenon in M. tuberculosis are the sigma factors, and regulators belonging to the WhiB, GntR, XRE, Mar and TetR families. Better understanding the impact of these regulators in survival to drug treatment might contribute to identify new drug targets and/or to design new strategies of intervention.

Uncovering Beta-Lactam Susceptibility Patterns in Clinical Isolates of Mycobacterium tuberculosis through Whole-Genome Sequencing

The increasing threat of drug resistance and a stagnated pipeline of novel therapeutics endanger the eradication of tuberculosis. Beta-lactams constitute promising additions to the current therapeutic arsenal and two carbapenems are included in group C of medicines recommended by the WHO for use in longer multidrug-resistant tuberculosis regimens. However, the determinants underlining diverse Mycobacterium tuberculosis phenotypes to beta-lactams remain largely undefined. To decipher these, we present a proof-of-concept study based on a large-scale beta-lactam susceptibility screening for 172 M. tuberculosis clinical isolates from Portugal, including 72 antimycobacterial drug-resistant strains. MICs were determined for multiple beta-lactams and strains were subjected to whole-genome sequencing to identify core-genome single-nucleotide variant-based profiles. Global and cell wall-targeted approaches were then followed to detect putative drivers of beta-lactam response. We found that drug-resistant strains were more susceptible to beta-lactams, but significant differences were not observed between distinct drug-resistance profiles. Sublineage 4.3.4.2 strains were significantly more susceptible to beta-lactams, while the contrary was observed for Beijing and 4.1.2.1 sublineages. While mutations in beta-lactamase or cell wall biosynthesis genes were uncommon, a rise in beta-lactam MICs was detected in parallel with the accumulation of mutations in peptidoglycan cross-linking or cell division genes. Finally, we exposed that putative beta-lactam resistance markers occurred in genes for which relevant roles in cell wall processes have been ascribed, such as rpfC or pknA. Genetic studies to validate the relevance of the identified mutations for beta-lactam susceptibility and further improvement of the phenotype-genotype associations are needed in the future.

IMPORTANCE Associations between differential M. tuberculosis beta-lactam phenotypes and preexisting antimycobacterial drug resistance, strain sublineage, or specific mutational patterns were established. Importantly, we reveal that highly drug-resistant isolates of sublineage 4.3.4.2 have an increased susceptibility to beta-lactams compared with other strains. Thus, directing beta-lactams to treat infections by specific M. tuberculosis strains and refraining its use from others emerges as a potentially important strategy to avoid resistance development. Individual mutations in blaC or genes encoding canonical beta-lactam targets, such as peptidoglycan transpeptidases, are infrequent and do not greatly impact the MICs of potent carbapenem plus clavulanic acid combinations. An improved understanding of the global effect of cumulative mutations in relevant gene sets for peptidoglycan and cell division processes on beta-lactam susceptibility is also provided.

Transmission, distribution and drug resistance-conferring mutations of extensively drug-resistant tuberculosis in the Western Cape Province, South Africa

Extensively drug-resistant tuberculosis (XDR-TB), defined as resistance to at least isoniazid (INH), rifampicin (RIF), a fluoroquinolone (FQ) and a second-line injectable drug (SLID), is difficult to treat and poses a major threat to TB control. The transmission dynamics and distribution of XDR Mycobacterium tuberculosis (Mtb) strains have not been thoroughly investigated. Using whole genome sequencing data on 461 XDR-Mtb strains, we aimed to investigate the geographical distribution of XDR-Mtb strains in the Western Cape Province of South Africa over a 10 year period (2006–2017) and assess the association between Mtb sub-lineage, age, gender, geographical patient location and membership or size of XDR-TB clusters. First, we identified transmission clusters by excluding drug resistance-conferring mutations and using the 5 SNP cutoff, followed by merging clusters based on their most recent common ancestor. We then consecutively included variants conferring resistance to INH, RIF, ethambutol (EMB), pyrazinamide (PZA), SLIDs and FQs in the cluster definition. Cluster sizes were classified as small (2–4 isolates), medium (5–20 isolates), large (21–100 isolates) or very large (>100 isolates) to reflect the success of individual strains. We found that most XDR-TB strains were clustered and that including variants conferring resistance to INH, RIF, EMB, PZA and SLIDs in the cluster definition did not significantly reduce the proportion of clustered isolates (85.5–82.2 %) but increased the number of patients belonging to small clusters (4.3–12.4 %, P=0.56). Inclusion of FQ resistance-conferring variants had the greatest effect, with 11 clustered isolates reclassified as unique while the number of clusters increased from 17 to 37. Lineage 2 strains (lineage 2.2.1 typical Beijing or lineage 2.2.2 atypical Beijing) showed the large clusters which were spread across all health districts of the Western Cape Province. We identified a significant association between residence in the Cape Town metropole and cluster membership (P=0.016) but no association between gender, age and cluster membership or cluster size (P=0.39). Our data suggest that the XDR-TB epidemic in South Africa probably has its origin in the endemic spread of MDR Mtb and pre-XDR Mtb strains followed by acquisition of FQ resistance, with more limited transmission of XDR Mtb strains. This only became apparent with the inclusion of drug resistance-conferring variants in the definition of a cluster. In addition to the prevention of amplification of resistance, rapid diagnosis of MDR, pre-XDR and XDR-TB and timely initiation of appropriate treatment is needed to reduce transmission of difficult-to-treat TB.

Practical approach to detection and surveillance of emerging highly resistant Mycobacterium tuberculosis Beijing 1071-32-cluster

Ancient sublineage of the Mycobacterium tuberculosis Beijing genotype is endemic and prevalent in East Asia and rare in other world regions. While these strains are mainly drug susceptible, we recently identified a novel clonal group Beijing 1071-32 within this sublineage emerging in Siberia, Russia and present in other Russian regions. This cluster included only multi/extensive drug resistant (MDR/XDR) isolates. Based on the phylogenetic analysis of the available WGS data, we identified three synonymous SNPs in the genes Rv0144, Rv0373c, and Rv0334 that were specific for the Beijing 1071-32-cluster and developed a real-time PCR assay for their detection. Analysis of the 2375 genetically diverse M. tuberculosis isolates collected between 1996 and 2020 in different locations (European and Asian parts of Russia, former Soviet Union countries, Albania, Greece, China, Vietnam, Japan and Brazil), confirmed 100% specificity and sensitivity of this real-time PCR assay. Moreover, the epidemiological importance of this strain and the newly developed screening assay is further stressed by the fact that all identified Beijing 1071-32 isolates were found to exhibit MDR genotypic profiles with concomitant resistance to additional first-line drugs due to a characteristic signature of six mutations in rpoB450, rpoC485, katG315, katG335, rpsL43 and embB497. In conclusion, this study provides a set of three concordant SNPs for the detection and screening of Beijing 1071-32 isolates along with a validated real-time PCR assay easily deployable across multiple settings for the epidemiological tracking of this significant MDR cluster.

Whole-genome sequencing as a tool for studying the microevolution of drug-resistant serial Mycobacterium tuberculosis isolates

Treatment of drug-resistant tuberculosis requires extended use of more toxic and less effective drugs and may result in retreatment cases due to failure, abandonment or disease recurrence. It is therefore important to understand the evolutionary process of drug resistance in Mycobacterium tuberculosis. We here in describe the microevolution of drug resistance in serial isolates from six previously treated patients. Drug resistance was initially investigated through phenotypic methods, followed by genotypic approaches. The use of whole-genome sequencing allowed the identification of mutations in the katG, rpsL and rpoB genes associated with drug resistance, including the detection of rare mutations in katG and mixed populations of strains. Molecular docking simulation studies of the impact of observed mutations on isoniazid binding were also performed. Whole-genome sequencing detected 266 single nucleotide polymorphisms between two isolates obtained from one patient, suggesting a case of exogenous reinfection. In conclusion, sequencing technologies can detect rare mutations related to drug resistance, identify subpopulations of resistant strains, and identify diverse populations of strains due to exogenous reinfection, thus improving tuberculosis control by guiding early implementation of appropriate clinical and therapeutic interventions.

Genetic diversity of candidate loci linked to Mycobacterium tuberculosis resistance to bedaquiline, delamanid and pretomanid

Tuberculosis (TB), caused by Mycobacterium tuberculosis, is one of the deadliest infectious diseases worldwide. Multidrug and extensively drug-resistant strains are making disease control difficult, and exhausting treatment options. New anti-TB drugs bedaquiline (BDQ), delamanid (DLM) and pretomanid (PTM) have been approved for the treatment of multi-drug resistant TB, but there is increasing resistance to them. Nine genetic loci strongly linked to resistance have been identified (mmpR5, atpE, and pepQ for BDQ; ddn, fgd1, fbiA, fbiB, fbiC, and fbiD for DLM/PTM). Here we investigated the genetic diversity of these loci across >33,000 M. tuberculosis isolates. In addition, epistatic mutations in mmpL5-mmpS5 as well as variants in ndh, implicated for DLM/PTM resistance in M. smegmatis, were explored. Our analysis revealed 1,227 variants across the nine genes, with the majority (78%) present in isolates collected prior to the roll-out of BDQ and DLM/PTM. We identified phylogenetically-related mutations, which are unlikely to be resistance associated, but also high-impact variants such as frameshifts (e.g. in mmpR5, ddn) with likely functional effects, as well as non-synonymous mutations predominantly in MDR-/XDR-TB strains with predicted protein destabilising effects. Overall, our work provides a comprehensive mutational catalogue for BDQ and DLM/PTM associated genes, which will assist with establishing associations with phenotypic resistance; thereby, improving the understanding of the causative mechanisms of resistance for these drugs, leading to better treatment outcomes.

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.

Exploratory development of PCR-fluorescent probes in rapid detection of mutations associated with extensively drug-resistant tuberculosis

This study aims to evaluate the clinical value of PCR-fluorescent probes for detecting the mutation gene associated with extensively drug-resistant tuberculosis (XDR-TB). The molecular species identification of 900 sputum specimens was performed using polymerase chain reaction (PCR)-fluorescent probe. The mutations of the drug resistance genes rpoB, katG, inhA, embB, rpsL, rrs, and gyrA were detected. The conventional drug susceptibility testing (DST) and PCR-directed sequencing (PCR-DS) were carried out as control. DST demonstrated that there were 501 strains of rifampicin resistance, 451 strains of isoniazid resistance, 293 strains of quinolone resistance, 425 strains of streptomycin resistance, 235 strains of ethambutol resistance, and 204 strains of amikacin resistance. Furthermore, 427 (47.44%) or 146 (16.22%) strains were MDR-TB or XDR-TB, respectively. The mutations of the rpoB, katG, inhA, embB, rpsL, rrs, and gyrA genes were detected in 751 of 900 TB patients by PCR-fluorescent probe method, and the rate of drug resistance was 751/900 (83.44%). No mutant genes were detected in the other 149 patients. Compared with DST, the mutant rates of rpoB, katG/inhA, rpsL, rrs, embB, and gyrA of six drugs were higher than 88%; five of six drugs were higher than 90% except for SM (88.11%). The MDR and XDR mutant gene types were found in 398 (42.22%) and 137 (15.22%) samples. PCR-DS was also employed and confirmed the PCR-fluorescent probe method with the accordance rate of 100%. The PCR-fluorescent probe method is rapid and straightforward in detecting XDR-TB genotypes and is worthy of being applied in hospitals.

Ancient Antibiotics, Ancient Resistance

As the spread of antibiotic resistance threatens our ability to treat infections, avoiding the return of a preantibiotic era requires the discovery of new drugs. While therapeutic use of antibiotics followed by the inevitable selection of resistance is a modern phenomenon, these molecules and the genetic determinants of resistance were in use by environmental microbes long before humans discovered them. In this review, we discuss evidence that antibiotics and resistance were present in the environment before anthropogenic use, describing techniques including direct sampling of ancient DNA and phylogenetic analyses that are used to reconstruct the past. We also pay special attention to the ecological and evolutionary forces that have shaped the natural history of antibiotic biosynthesis, including a discussion of competitive versus signaling roles for antibiotics, proto-resistance, and substrate promiscuity of biosynthetic and resistance enzymes. Finally, by applying an evolutionary lens, we describe concepts governing the origins and evolution of biosynthetic gene clusters and cluster-associated resistance determinants. These insights into microbes' use of antibiotics in nature, a game they have been playing for millennia, can provide inspiration for discovery technologies and management strategies to combat the growing resistance crisis.

The 'thanato-resistome' - The funeral industry as a potential reservoir of antibiotic resistance: Early insights and perspectives

The funeral industry is a potential reservoir of antibiotic resistance. The occurrence, human exposure and health risks of antibiotic resistance in the funeral industry were examined. The funeral industry harbours antibiotic resistance to multiple common and last-resort antibiotics, hence constitutes the 'thanato-resistome'. Hydrological processes, air-borne particulates and vectors disseminate antibiotic resistance, while horizontal gene transfer circulates antibiotic resistance among resistomes, forming a complex network. Ingestion, inhalation of air-borne particulates, dermal intake and clothes of workers contribute to human exposure. Human health risks include; development of drug resistance in previously susceptible pathogens, and increased morbidity and mortality caused by increased pathogenicity and outbreaks of multi-drug resistant infections. Ecological risks include the proliferation of resistant organisms at the expense of susceptible ones, thereby disrupting ecosystem structure and function, including biogeochemical cycles. Barring inferential data, quantitative evidence linking antibiotic resistance to human infections is weak. This reflects the lack of systematic quantitative studies, rather than the absence of such health risks. Quantitative risk assessment is constrained by lack of quantitative data on antibiotic resistance in various reservoirs and exposure routes. A framework for risk assessment and mitigation is proposed. Finally, ten hypotheses and emerging tools such as genomics, in silico techniques and big data analytics are highlighted.

Genomic epidemiology of Mycobacterium tuberculosis in Santa Catarina, Southern Brazil

Mycobacterium tuberculosis (M.tb), the pathogen responsible for tuberculosis (TB) poses as the major cause of death among infectious diseases. The knowledge about the molecular diversity of M.tb enables the implementation of more effective surveillance and control measures and, nowadays, Whole Genome Sequencing (WGS) holds the potential to produce high-resolution epidemiological data in a high-throughput manner. Florianópolis, the state capital of Santa Catarina (SC) in south Brazil, shows a high TB incidence (46.0/100,000). Here we carried out a WGS-based evaluation of the M.tb strain diversity, drug-resistance and ongoing transmission in the capital metropolitan region. Resistance to isoniazid, rifampicin, streptomycin was identified respectively in 4.0% (n = 6), 2.0% (n = 3) and 1.3% (n = 2) of the 151 studied strains by WGS. Besides, resistance to pyrazinamide and ethambutol was detected in 0.7% (n = 1) and reistance to ethionamide and fluoroquinolone (FQ) in 1.3% (n = 2), while a single (0.7%) multidrug-resistant (MDR) strain was identified. SNP-based typing classified all isolates into M.tb Lineage 4, with high proportion of sublineages LAM (60.3%), T (16.4%) and Haarlem (7.9%). The average core-genome distance between isolates was 420.3 SNPs, with 43.7% of all isolates grouped across 22 genomic clusters thereby showing the presence of important ongoing TB transmission events. Most clusters were geographically distributed across the study setting which highlights the need for an urgent interruption of these large transmission chains. The data conveyed by this study shows the presence of important and uncontrolled TB transmission in the metropolitan area and provides precise data to support TB control measures in this region.

Genomic Analysis Identifies Mutations Concerning Drug-Resistance and Beijing Genotype in Multidrug-Resistant Mycobacterium tuberculosis Isolated From China

Development of modern genomics provides us an effective method to understand the molecular mechanism of drug resistance and diagnose drug-resistant Mycobacterium tuberculosis. In this study, mutations in 18 genes or intergenic regions acquired by whole-genome sequencing (WGS) of 183 clinical M. tuberculosis strains, including 137 multidrug-resistant and 46 pan-susceptible isolates from China, were identified and used to analyze their associations with resistance of isoniazid, rifampin, ethambutol, and streptomycin. Using the proportional method as the gold standard method, the accuracy values of WGS to predict resistance were calculated. The association between synonymous or lineage definition mutations with different genotypes were also analyzed. The results show that, compared to the phenotypic proportional method, the sensitivity and specificity of WGS for resistance detection were 94.2 and 100.0% for rifampicin (based on mutations in rpoB), 90.5 and 97.8% for isoniazid (katG), 83.0 and 97.8% for streptomycin (rpsL combined with rrs 530 loop and 912 loop), and 90.9 and 65.1% for ethambutol (embB), respectively. WGS data also showed that mutations in the inhA promoter increased only 2.2% sensitivity for INH based on mutations in katG. Synonymous mutation rpoB A1075A was confirmed to be associated with the Beijing genotype. This study confirmed that mutations in rpoB, katG, rrs 530 loop and 912 loop, and rpsL were excellent biomarkers for predicting rifampicin, isoniazid, and streptomycin resistance, respectively, and provided clues in clarifying the drug-resistance mechanism of M. tuberculosis isolates from China.