Mapping and manipulating the Mycobacterium tuberculosis transcriptome using a transcription factor overexpression-derived regulatory network

Background

Mycobacterium tuberculosis is characterized by a low mutation rate and a lack of genetic recombination. Yet, the rise of extensively resistant strains paints a picture of a microbe with an impressive adaptive potential. Here we describe the first documented case of extensively drug-resistant tuberculosis evolved from a susceptible ancestor within a single patient.

Results

Genome sequences of nine serial M. tuberculosis isolates from the same patient uncovered a dramatic turnover of competing lineages driven by the emergence, and subsequent fixation or loss of single nucleotide polymorphisms. For most drugs, resistance arose through independent emergence of mutations in more than one clone, of which only one ultimately prevailed as the clone carrying it expanded, displacing the other clones in the process. The vast majority of mutations identified over 3.5 years were either involved in drug resistance or hitchhiking in the genetic background of these. Additionally, RNA-sequencing of isolates grown in the absence of drug challenge revealed that the efflux-associated iniBAC operon was up-regulated over time, whereas down-regulated genes include those involved in mycolic acid synthesis.

Conclusions

We observed both rapid acquisitions of resistance to antimicrobial compounds mediated by individual mutations as well as a gradual increase in fitness in the presence of antibiotics, likely driven by stable gene expression reprogramming. The rapid turnover of resistance mutations and hitchhiking neutral mutations has major implications for inferring tuberculosis transmission events in situations where drug resistance evolves within transmission chains.     

Mutations in dnaA and a cryptic interaction site increase drug resistance in Mycobacterium tuberculosis

Genomic dissection of antibiotic resistance in bacterial pathogens has largely focused on genetic changes conferring growth above a single critical concentration of drug. However, reduced susceptibility to antibiotics—even below this breakpoint—is associated with poor treatment outcomes in the clinic, including in tuberculosis. Clinical strains of Mycobacterium tuberculosis exhibit extensive quantitative variation in antibiotic susceptibility but the genetic basis behind this spectrum of drug susceptibility remains ill-defined. Through a genome wide association study, we show that non-synonymous mutations in dnaA, which encodes an essential and highly conserved regulator of DNA replication, are associated with drug resistance in clinical M. tuberculosis strains. We demonstrate that these dnaA mutations specifically enhance M. tuberculosis survival during isoniazid treatment via reduced expression of katG, the activator of isoniazid. To identify DnaA interactors relevant to this phenotype, we perform the first genome-wide biochemical mapping of DnaA binding sites in mycobacteria which reveals a DnaA interaction site that is the target of recurrent mutation in clinical strains. Reconstructing clinically prevalent mutations in this DnaA interaction site reproduces the phenotypes of dnaA mutants, suggesting that clinical strains of M. tuberculosis have evolved mutations in a previously uncharacterized DnaA pathway that quantitatively increases resistance to the key first-line antibiotic isoniazid. Discovering genetic mechanisms that reduce drug susceptibility and support the evolution of high-level drug resistance will guide development of biomarkers capable of prospectively identifying patients at risk of treatment failure in the clinic.     

Identification of drug resistance mutations among Mycobacterium bovis lineages in the Americas

Identifying the Mycobacterium tuberculosis resistance mutation patterns is of the utmost importance to assure proper patient’s management and devising of control programs aimed to limit spread of disease. Zoonotic Mycobacterium bovis infection still represents a threat to human health, particularly in dairy production regions. Routinary, molecular characterization of M. bovis is performed primarily by spoligotyping and mycobacterial interspersed repetitive units (MIRU) while next generation sequencing (NGS) approaches are often performed by reference laboratories. However, spoligotyping and MIRU methodologies lack the resolution required for the fine characterization of tuberculosis isolates, particularly in outbreak settings. In conjunction with sophisticated bioinformatic algorithms, whole genome sequencing (WGS) analysis is becoming the method of choice for advanced genetic characterization of tuberculosis isolates. WGS provides valuable information on drug resistance and compensatory mutations that other technologies cannot assess. Here, we performed an analysis of the most frequently identified mutations associated with tuberculosis drug resistance and their genetic relationship among 2,074 Mycobacterium bovis WGS recovered primarily from non-human hosts. Full-length gene sequences harboring drug resistant associated mutations and their phylogenetic relationships were analyzed. The results showed that M. bovis isolates harbor mutations conferring resistance to both first- and second-line antibiotics. Mutations conferring resistance for isoniazid, fluoroquinolones, streptomycin, and aminoglycosides were identified among animal strains. Our findings highlight the importance of molecular surveillance to monitor the emergence of mutations associated with multi and extensive drug resistance in livestock and other non-human mammals.     

Comparison of Gene Expression Profiles Between Pansensitive and Multidrug-Resistant Strains of Mycobacterium tuberculosis

Mycobacterium tuberculosis has developed resistance to anti-tuberculosis first-line drugs. Multidrug-resistant strains complicate the control of tuberculosis and have converted it into a worldwide public health problem. Mutational studies of target genes have tried to envisage the resistance in clinical isolates; however, detection of these mutations in some cases is not sufficient to identify drug resistance, suggesting that other mechanisms are involved. Therefore, the identification of new markers of susceptibility or resistance to first-line drugs could contribute (1) to specifically diagnose the type of M. tuberculosis strain and prescribe an appropriate therapy, and (2) to elucidate the mechanisms of resistance in multidrug-resistant strains. In order to identify specific genes related to resistance in M. tuberculosis, we compared the gene expression profiles between the pansensitive H37Rv strain and a clinical CIBIN:UMF:15:99 multidrug-resistant isolate using microarray analysis. Quantitative real-time PCR confirmed that in the clinical multidrug-resistant isolate, the esxG, esxH, rpsA, esxI, and rpmI genes were upregulated, while the lipF, groES, and narG genes were downregulated. The modified genes could be involved in the mechanisms of resistance to first-line drugs in M. tuberculosis and could contribute to increased efficiency in molecular diagnosis approaches of infections with drug-resistant strains.     

A genome-wide analysis of multidrug-resistant and extensively drug-resistant strains of Mycobacterium tuberculosis Beijing genotype

The Beijing genotype of Mycobacterium tuberculosis (MTB) is one of the most successful MTB lineages that has disseminated in the world. In China, the rate of multidrug-resistant (MDR) tuberculosis is significantly higher than the global average rate, and the Beijing genotype strains take the largest share of MDR strains. To study the genetic basis of the epidemiological findings that Beijing genotype has often been associated with tuberculosis outbreaks and drug resistance, we determined the genome sequences of four clinical isolates: two extensively drug resistant (XDR1219, XDR1221) and two multidrug resistant (WX1, WX3), using whole-genome sequencing. A large number of individual and shared SNPs of the four Beijing strains were identified. Our isolates harbored almost all classic drug resistance-associated mutations. The mutations responsible for drug resistance in the two XDR strains were consistent with the clinical quantitative drug resistance levels. COG analysis revealed that Beijing strains have significantly higher abundances of the mutations responsible for cell wall/membrane/envelope biogenesis (COG M), secondary metabolites biosynthesis, transport and catabolism (COG Q), lipid transport and metabolism (COG I) and defense mechanisms (COG V). The shared mutated genes of the four studied Beijing strains were significantly overrepresented in three DNA repair pathways. Our analyses promote the understanding of the genome polymorphism of the Beijing family strains and provide the molecular genetic basis for their wide dissemination capacity and drug resistance.     

Comparative Genomic and Proteomic Anatomy of Mycobacterium Ubiquitous Esx Family Proteins: Implications in Pathogenicity and Virulence

Secreted proteins are among the most important molecules involved in host—pathogen interaction of Mycobacterium tuberculosis, the etiological agent of human tuberculosis (TB). M. tuberculosis encodes five types of VII secretion systems (ESX-1 to ESX-5) responsible for the exportation of many proteins. This system mediated substrates including members of the Esx family implicated in tuberculosis pathogenesis and survival within host cells. However, the distribution and evolution of this family remain elusive. To explore the evolution and distribution of Esx family proteins, we analyzed all available Mycobacteria genomes. Interestingly, amino mutations among M. tuberculosis esx family proteins may relate to their functions. We further analyzed the differences between pathogenic Mycobacteria, the attenuated Mycobacteria and non-pathogenic Mycobacteria. The stability, the globular domains and the phosphorylation of serine/threonine residues of M. tuberculosis esx proteins with their homologies among other Mycoabcteria were analyzed. Our comparative genomic and proteomic analysis found that the change of stability, gain or loss of globular domains and phosphorylation of serine/threonine might be responsible for the difference between the pathogenesis and virulence of the esx proteins and its homolog widespread among Mycobacteria and related species, which may provide clues for novel anti-tuberculosis drug targets.     

Rapid Molecular Detection of Multidrug-Resistant Tuberculosis by PCR-Nucleic Acid Lateral Flow Immunoassay

Several existing molecular tests for multidrug-resistant tuberculosis (MDR-TB) are limited by complexity and cost, hindering their widespread application. The objective of this proof of concept study was to develop a simple Nucleic Acid Lateral Flow (NALF) immunoassay as a potential diagnostic alternative, to complement conventional PCR, for the rapid molecular detection of MDR-TB. The NALF device was designed using antibodies for the indirect detection of labeled PCR amplification products. Multiplex PCR was optimized to permit the simultaneous detection of the drug resistant determining mutations in the 81-bp hot spot region of the rpoB gene (rifampicin resistance), while semi-nested PCR was optimized for the S315T mutation detection in the katG gene (isoniazid resistance). The amplification process additionally targeted a conserved region of the genes as Mycobacterium tuberculosis (Mtb) DNA control. The optimized conditions were validated with the H37Rv wild-type (WT) Mtb isolate and Mtb isolates with known mutations (MT) within the rpoB and katG genes. Results indicate the correct identification of WT (drug susceptible) and MT (drug resistant) Mtb isolates, with the least limit of detection (LOD) being 10⁴ genomic copies per PCR reaction. NALF is a simple, rapid and low-cost device suitable for low resource settings where conventional PCR is already employed on a regular basis. Moreover, the use of antibody-based NALF to target primer-labels, without the requirement for DNA hybridization, renders the device generic, which could easily be adapted for the molecular diagnosis of other infectious and non-infectious diseases requiring nucleic acid detection.     

Efflux Pump Gene Expression in Multidrug-Resistant Mycobacterium tuberculosis Clinical Isolates

Isoniazid (INH) and rifampicin (RIF) are the two most effective drugs in tuberculosis therapy. Understanding the molecular mechanisms of resistance to these two drugs is essential to quickly diagnose multidrug-resistant (MDR) tuberculosis and extensive drug-resistant tuberculosis. Nine clinical Mycobacterium tuberculosis isolates resistant to only INH and RIF and 10 clinical pan-sensitive isolates were included to evaluate the expression of 20 putative drug efflux pump genes and sequence mutations in rpoB (RIF), katG (INH), the inhA promoter (INH), and oxyR-ahpC (INH). Nine and three MDR isolates were induced to overexpress efflux pump genes by INH and RIF, respectively. Eight and two efflux pump genes were induced to overexpress by INH and RIF in MDR isolates, respectively. drrA, drrB, efpA, jefA (Rv2459), mmr, Rv0849, Rv1634, and Rv1250 were overexpressed under INH or RIF stress. Most efflux pump genes were overexpressed under INH stress in a MDR isolates that carried the wild-type katG, inhA, and oxyR-ahpC associated with INH resistance than in those that carried mutations. The expression levels of 11 genes (efpA, Rv0849, Rv1250, P55 (Rv1410c), Rv1634, Rv2994, stp, Rv2459, pstB, drrA, and drrB) without drug inducement were significantly higher (P < 0.05) in nine MDR isolates than in 10 pan-sensitive isolates. In conclusion, efflux pumps may play an important role in INH acquired resistance in MDR M. tuberculosis, especially in those strains having no mutations in genes associated with INH resistance; basal expression levels of some efflux pump genes are higher in MDR isolates than in pan-sensitive isolates and the basal expressional differences may be helpful to diagnose and treat resistant tuberculosis.     

Detection of Rifampicin Resistance in Mycobacterium tuberculosis by Padlock Probes and Magnetic Nanobead-Based Readout

Control of the global epidemic tuberculosis is severely hampered by the emergence of drug-resistant Mycobacterium tuberculosis strains. Molecular methods offer a more rapid means of characterizing resistant strains than phenotypic drug susceptibility testing. We have developed a molecular method for detection of rifampicin-resistant M. tuberculosis based on padlock probes and magnetic nanobeads. Padlock probes were designed to target the most common mutations associated with rifampicin resistance in M. tuberculosis, i.e. at codons 516, 526 and 531 in the gene rpoB. For detection of the wild type sequence at all three codons simultaneously, a padlock probe and two gap-fill oligonucleotides were used in a novel assay configuration, requiring three ligation events for circularization. The assay also includes a probe for identification of the M. tuberculosis complex. Circularized probes were amplified by rolling circle amplification. Amplification products were coupled to oligonucleotide-conjugated magnetic nanobeads and detected by measuring the frequency-dependent magnetic response of the beads using a portable AC susceptometer.     

Genome Analysis of Multi- and Extensively-Drug-Resistant Tuberculosis from KwaZulu-Natal,South Africa

The KZN strain family of Mycobacterium tuberculosis is a highly virulent strain endemic to the KwaZulu-Natal region of South Africa, which has recently experienced an outbreak of extensively-drug resistant tuberculosis. To investigate the causes and evolution of drug-resistance, we determined the DNA sequences of several clinical isolates - one drug-susceptible, one multi-drug resistant, and nine extensively drug-resistant - using whole-genome sequencing. Analysis of polymorphisms among the strains is consistent with the drug-susceptibility profiles, in that well-known mutations are observed that are correlated with resistance to isoniazid, rifampicin, kanamycin, ofloxacin, ethambutol, and pyrazinamide. However, the mutations responsible for rifampicin resistance in rpoB and pyrazinamide in pncA are in different nucleotide positions in the multi-drug-resistant and extensively drug-resistant strains, clearly showing that they acquired these mutations independently, and that the XDR strain could not have evolved directly from the MDR strain (though it could have arisen from another similar MDR strain). Sequencing of eight additional XDR strains from other areas of KwaZulu-Natal shows that they have identical drug resistant mutations to the first one sequenced, including the same polymorphisms at sites associated with drug resistance, supporting the theory that this represents a case of clonal expansion.     

Resistome analysis of Mycobacterium tuberculosis: Identification of aminoglycoside 2'-Nacetyltransferase (AAC) as co-target for drug desigining

The emergence of multidrug resistant tuberculosis (MDRTB) highlights the urgent need to understand the mechanisms of resistance to the drugs and to develop a new arena of therapeutics to treat the disease. Ethambutol, isonazid, pyrazinamide, rifampicin are first line of drugs against TB, whereas aminoglycoside, polypeptides, fluoroquinolone, ethionamide are important second line of bactericidal drugs used to treat MDRTB, and resistance to one or both of these drugs are defining characteristic of extensively drug resistant TB. We retrieved 1,221 resistant genes from Antibiotic Resistance Gene Database (ARDB), which are responsible for resistance against first and second line antibiotics used in treatment of Mycobacterium tuberculosis infection. From network analysis of these resistance genes, 53 genes were found to be common. Phylogenetic analysis shows that more than 60% of these genes code for acetyltransferase. Acetyltransferases detoxify antibiotics by acetylation, this mechanism plays central role in antibiotic resistance. Seven acetyltransferase (AT-1 to AT-7) were selected from phylogenetic analysis. Structural alignment shows that these acetyltransferases share common ancestral core, which can be used as a template for structure based drug designing. From STRING analysis it is found that acetyltransferase interact with 10 different proteins and it shows that, all these interaction were specific to M. tuberculosis. These results have important implications in designing new therapeutic strategies with acetyltransferase as lead co-target to combat against MDR as well as Extreme drug resistant (XDR) tuberculosis.

Abbreviations

AA - amino acid, AT - Acetyltransferase, AAC - Aminoglycoside 2''-N-acetyltransferase, XDR - Extreme drug-resistant, MDR - Multidrug-resistant, Mtb - Mycobacterium tuberculosis, TB - Tuberculosis.     

Directed mutagenesis of Mycobacterium smegmatis 16S rRNA to reconstruct the in vivo evolution of aminoglycoside resistance in Mycobacterium tuberculosis

Drug resistance in Mycobacterium tuberculosis is a global problem, with major consequences for treatment and public health systems. As the emergence and spread of drug‐resistant tuberculosis epidemics is largely influenced by the impact of the resistance mechanism on bacterial fitness, we wished to investigate whether compensatory evolution occurs in drug‐resistant clinical isolates of M. tuberculosis. By combining information from molecular epidemiology studies of drug‐resistant clinical M. tuberculosis isolates with genetic reconstructions and measurements of aminoglycoside susceptibility and fitness in Mycobacterium smegmatis, we have reconstructed a plausible pathway for how aminoglycoside resistance develops in clinical isolates of M. tuberculosis. Thus, we show by reconstruction experiments that base changes in the highly conserved A‐site of 16S rRNA that: (i) cause aminoglycoside resistance, (ii) confer a high fitness cost and (iii) destabilize a stem‐loop structure, are associated with a particular compensatory point mutation that restores rRNA secondary structure and bacterial fitness, while maintaining to a large extent the drug‐resistant phenotype. The same types of resistance and associated mutations can be found in M. tuberculosis in clinical isolates, suggesting that compensatory evolution contributes to the spread of drug‐resistant tuberculosis disease.     

Locked Nucleic Acid Probe-Based Real-Time PCR Assay for the Rapid Detection of Rifampin-Resistant Mycobacterium tuberculosis

Drug-resistant Mycobacterium tuberculosis can be rapidly diagnosed through nucleic acid amplification techniques by analyzing the variations in the associated gene sequences. In the present study, a locked nucleic acid (LNA) probe-based real-time PCR assay was developed to identify the mutations in the rpoB gene associated with rifampin (RFP) resistance in M. tuberculosis. Six LNA probes with the discrimination capability of one-base mismatch were designed to monitor the 23 most frequent rpoB mutations. The target mutations were identified using the probes in a “probe dropout” manner (quantification cycle = 0); thus, the proposed technique exhibited superiority in mutation detection. The LNA probe-based real-time PCR assay was developed in a two-tube format with three LNA probes and one internal amplification control probe in each tube. The assay showed excellent specificity to M. tuberculosis with or without RFP resistance by evaluating 12 strains of common non-tuberculosis mycobacteria. The limit of detection of M. tuberculosis was 10 genomic equivalents (GE)/reaction by further introducing a nested PCR method. In a blind validation of 154 clinical mycobacterium isolates, 142/142 (100%) were correctly detected through the assay. Of these isolates, 88/88 (100%) were determined as RFP susceptible and 52/54 (96.3%) were characterized as RFP resistant. Two unrecognized RFP-resistant strains were sequenced and were found to contain mutations outside the range of the 23 mutation targets. In conclusion, this study established a sensitive, accurate, and low-cost LNA probe-based assay suitable for a four-multiplexing real-time PCR instrument. The proposed method can be used to diagnose RFP-resistant tuberculosis in clinical laboratories.     

Molecular characterisation of Mycobacterium tuberculosis isolates in the First National Survey of Anti-tuberculosis Drug Resistance from Venezuela

Background

Molecular typing of Mycobacterium tuberculosis strains has become a valuable tool in the epidemiology of tuberculosis (TB) by allowing detection of outbreaks, tracking of epidemics, identification of genotypes and transmission events among patients who would have remained undetected by conventional contact investigation. This is the first genetic biodiversity study of M. tuberculosis in Venezuela. Thus, we investigated the genetic patterns of strains isolated in the first survey of anti-tuberculosis drug-resistance realised as part of the Global Project of Anti-tuberculosis Drug Resistance Surveillance (WHO/IUATLD).

Results

Clinical isolates (670/873) were genotyped by spoligotyping. The results were compared with the international spoligotyping database (SpolDB4). Multidrug resistant (MDR) strains (14/18) were also analysed by IS6110-RFLP assays, and resistance to isoniazid and rifampicin was characterised. Spoligotyping grouped 82% (548/670) of the strains into 59 clusters. Twenty new spoligotypes (SITs) specific to Venezuela were identified. Eight new inter-regional clusters were created. The Beijing genotype was not found. The genetic network shows that the Latin American and Mediterranean family constitutes the backbone of the genetic TB population-structure in Venezuela, responsible of >60% of total TB cases studied. MDR was 0.5% in never treated patients and 13.5% in previously treated patients. Mutations in rpoB gene and katG genes were detected in 64% and 43% of the MDR strains, respectively. Two clusters were found to be identical by the four different analysis methods, presumably representing cases of recent transmission of MDR tuberculosis.

Conclusion

This study gives a first overview of the M. tuberculosis strains circulating in Venezuela during the first survey of anti-tuberculosis drug-resistance. It may aid in the creation of a national database that will be a valuable support for further studies.     

Multidrug Resistant Mycobacterium tuberculosis: A Retrospective katG and rpoB Mutation Profile Analysis in Isolates from a Reference Center in Brazil

Background

Multidrug resistance is a critical factor in tuberculosis control. To gain better understanding of multidrug resistant tuberculosis in Brazil, a retrospective study was performed to compare genotypic diversity and drug resistance associated mutations in Mycobacterium tuberculosis isolates from a national reference center.

Methods and Findings

Ninety-nine multidrug resistant isolates from 12 Brazilian states were studied. Drug-resistance patterns were determined and the rpoB and katG genes were screened for mutations. Genotypic diversity was investigated by IS6110-RFLP and Luminex 47 spoligotyping. Mutations in rpoB and katG were seen in 91% and 93% of the isolates, respectively. Codon 315 katG mutations occurred in 82.8% of the isolates with a predominance of the Ser315Thr substitution. Twenty-five isolates were clustered in 11 groups with identical IS6110-RFLP patterns while 74 showed unique patterns with no association between mutation frequencies or susceptibility profiles. The most prevalent spoligotyping lineages were LAM (47%), T (17%) and Haarlen (12%). The Haarlen lineage showed a higher frequency of codon 516 rpoB mutations while codon 531 mutations prevailed in the other isolates.

Conclusions

Our data suggest that there were no major multidrug resistant M. tuberculosis strains transmitted among patients referred to the reference center, indicating an independent acquisition of resistance. In addition, drug resistance associated mutation profiles were well established among the main spoligotyping lineages found in these Brazilian multidrug resistant isolates, providing useful data for patient management and treatment.     

Pre-existing isoniazid resistance, but not the genotype of Mycobacterium tuberculosis drives rifampicin resistance codon preference in vitro

Both the probability of a mutation occurring and the ability of the mutant to persist will influence the distribution of mutants that arise in a population. We studied the interaction of these factors for the in vitro selection of rifampicin (RIF)-resistant mutants of Mycobacterium tuberculosis. We characterised two series of spontaneous RIF-resistant in vitro mutants from isoniazid (INH)-sensitive and -resistant laboratory strains and clinical isolates, representing various M. tuberculosis genotypes. The first series were selected from multiple parallel 1 ml cultures and the second from single 10 ml cultures. RIF-resistant mutants were screened by Multiplex Ligation-dependent Probe Amplification (MLPA) or by sequencing the rpoB gene. For all strains the mutation rate for RIF resistance was determined with a fluctuation assay. The most striking observation was a shift towards rpoB-S531L (TCG→TTG) mutations in a panel of laboratory-generated INH-resistant mutants selected from the 10-ml cultures (p<0.001). All tested strains showed similar mutation rates (1.33×10⁻⁸ to 2.49×10⁻⁷) except one of the laboratory-generated INH mutants with a mutation rate measured at 5.71×10⁻⁷, more than 10 times higher than that of the INH susceptible parental strain (5.46–7.44×10⁻⁸). No significant, systematic difference in the spectrum of rpoB-mutations between strains of different genotypes was observed. The dramatic shift towards rpoB-S531L in our INH-resistant laboratory mutants suggests that the relative fitness of resistant mutants can dramatically impact the distribution of (subsequent) mutations that accumulate in a M. tuberculosis population, at least in vitro. We conclude that, against specific genetic backgrounds, certain resistance mutations are particularly likely to spread. Molecular screening for these (combinations of) mutations in clinical isolates could rapidly identify these particular pathogenic strains. We therefore recommend that isolates are screened for the distribution of resistance mutations, especially in regions that are highly endemic for (multi)drug resistant tuberculosis.     

A molecular platform for the diagnosis of multidrug-resistant and pre-extensively drug-resistant tuberculosis based on single nucleotide polymorphism mutations present in Colombian isolates of Mycobacterium tuberculosis

Developing a fast, inexpensive, and specific test that reflects the mutations present in Mycobacterium tuberculosis isolates according to geographic region is the main challenge for drug-resistant tuberculosis (TB) control. The objective of this study was to develop a molecular platform to make a rapid diagnosis of multidrug-resistant (MDR) and extensively drug-resistant TB based on single nucleotide polymorphism (SNP) mutations present in therpoB, katG, inhA,ahpC, and gyrA genes from Colombian M. tuberculosis isolates. The amplification and sequencing of each target gene was performed. Capture oligonucleotides, which were tested before being used with isolates to assess the performance, were designed for wild type and mutated codons, and the platform was standardised based on the reverse hybridisation principle. This method was tested on DNA samples extracted from clinical isolates from 160 Colombian patients who were previously phenotypically and genotypically characterised as having susceptible or MDR M. tuberculosis. For our method, the kappa index of the sequencing results was 0,966, 0,825, 0,766, 0,740, and 0,625 forrpoB, katG, inhA,ahpC, and gyrA, respectively. Sensitivity and specificity were ranked between 90-100% compared with those of phenotypic drug susceptibility testing. Our assay helps to pave the way for implementation locally and for specifically adapted methods that can simultaneously detect drug resistance mutations to first and second-line drugs within a few hours.     

A single-step sequencing method for the identification of Mycobacterium tuberculosis complex species

Background

The Mycobacterium tuberculosis complex (MTC) comprises closely related species responsible for strictly human and zoonotic tuberculosis. Accurate species determination is useful for the identification of outbreaks and epidemiological links. Mycobacterium africanum and Mycobacterium canettii are typically restricted to Africa and M. bovis is a re-emerging pathogen. Identification of these species is difficult and expensive.

Methodology/Principal Findings

The Exact Tandem Repeat D (ETR-D; alias Mycobacterial Interspersed Repetitive Unit 4) was sequenced in MTC species type strains and 110 clinical isolates, in parallel to reference polyphasic identification based on phenotype profiling and sequencing of pncA, oxyR, hsp65, gyrB genes and the major polymorphism tandem repeat. Inclusion of M. tuberculosis isolates in the expanding, antibiotic-resistant Beijing clone was determined by Rv0927c gene sequencing. The ETR-D (780-bp) sequence unambiguously identified MTC species type strain except M. pinnipedii and M. microti thanks to six single nucleotide polymorphisms, variable numbers (1–7 copies) of the tandem repeat and two deletions/insertions. The ETR-D sequencing agreed with phenotypic identification in 107/110 clinical isolates and with reference polyphasic molecular identification in all isolates, comprising 98 M. tuberculosis, 5 M. bovis BCG type, 5 M. canettii, and 2 M. africanum. For M. tuberculosis isolates, the ETR-D sequence was not significantly associated with the Beijing clone.

Conclusions/Significance

ETR-D sequencing allowed accurate, single-step identification of the MTC at the species level. It circumvented the current expensive, time-consuming polyphasic approach. It could be used to depict epidemiology of zoonotic and human tuberculosis, especially in African countries where several MTC species are emerging.     

Phylogeography and transmission of M. tuberculosis in Moldova: A prospective genomic analysis

BackgroundThe incidence of multidrug-resistant tuberculosis (MDR-TB) remains critically high in countries of the former Soviet Union, where >20% of new cases and >50% of previously treated cases have resistance to rifampin and isoniazid. Transmission of resistant strains, as opposed to resistance selected through inadequate treatment of drug-susceptible tuberculosis (TB), is the main driver of incident MDR-TB in these countries.Methods and findingsWe conducted a prospective, genomic analysis of all culture-positive TB cases diagnosed in 2018 and 2019 in the Republic of Moldova. We used phylogenetic methods to identify putative transmission clusters; spatial and demographic data were analyzed to further describe local transmission of Mycobacterium tuberculosis. Of 2,236 participants, 779 (36%) had MDR-TB, of whom 386 (50%) had never been treated previously for TB. Moreover, 92% of multidrug-resistant M. tuberculosis strains belonged to putative transmission clusters. Phylogenetic reconstruction identified 3 large clades that were comprised nearly uniformly of MDR-TB: 2 of these clades were of Beijing lineage, and 1 of Ural lineage, and each had additional distinct clade-specific second-line drug resistance mutations and geographic distributions. Spatial and temporal proximity between pairs of cases within a cluster was associated with greater genomic similarity. Our study lasted for only 2 years, a relatively short duration compared with the natural history of TB, and, thus, the ability to infer the full extent of transmission is limited.ConclusionsThe MDR-TB epidemic in Moldova is associated with the local transmission of multiple M. tuberculosis strains, including distinct clades of highly drug-resistant M. tuberculosis with varying geographic distributions and drug resistance profiles. This study demonstrates the role of comprehensive genomic surveillance for understanding the transmission of M. tuberculosis and highlights the urgency of interventions to interrupt transmission of highly drug-resistant M. tuberculosis.

In a prospective genome surveillance study, Chongguang Yang and colleagues investigate the dynamics of multidrug-resistant tuberculosis transmission in Moldova.