Evolution of extensively drug-resistant Mycobacterium tuberculosis from a susceptible ancestor in a single patient

Ecological and evolutionary mechanisms driving within-patient emergence of antimicrobial resistance

The ecological and evolutionary mechanisms of antimicrobial resistance (AMR) emergence within patients and how these vary across bacterial infections are poorly understood. Increasingly widespread use of pathogen genome sequencing in the clinic enables a deeper understanding of these processes. In this Review, we explore the clinical evidence to support four major mechanisms of within-patient AMR emergence in bacteria: spontaneous resistance mutations; in situ horizontal gene transfer of resistance genes; selection of pre-existing resistance; and immigration of resistant lineages. Within-patient AMR emergence occurs across a wide range of host niches and bacterial species, but the importance of each mechanism varies between bacterial species and infection sites within the body. We identify potential drivers of such differences and discuss how ecological and evolutionary analysis could be embedded within clinical trials of antimicrobials, which are powerful but underused tools for understanding why these mechanisms vary between pathogens, infections and individuals. Ultimately, improving understanding of how host niche, bacterial species and antibiotic mode of action combine to govern the ecological and evolutionary mechanism of AMR emergence in patients will enable more predictive and personalized diagnosis and antimicrobial therapies.

Genetic mechanisms of co-emergence of INH-resistant Mycobacterium tuberculosis strains during the standard course of antituberculosis therapy

The incidence of isoniazid (INH) resistant Mycobacterium tuberculosis is increasing globally. This study aimed to identify the molecular mechanisms behind the development of INH resistance in M. tuberculosis strains collected from the same patients during the standard course of treatment. Three M. tuberculosis strains were collected from a patient before and during antituberculosis (anti-TB) therapy. The strains were characterized using phenotypic drug susceptibility tests, Mycobacterial Interspersed Repeated Unit-Variable-Number Tandem Repeats (MIRU-VNTR), and whole-genome sequencing (WGS) to identify mutations associated with INH resistance. To validate the role of the novel mutations in INH resistance, the mutated katG genes were electroporated into a KatG-deleted M. tuberculosis strain (GA03). Three-dimensional structures of mutated KatG were modeled to predict their impact on INH binding. The pre-treatment strain was susceptible to INH. However, two INH-resistant strains were isolated from the patient after anti-TB therapy. MIRU-VNTR and WGS revealed that the three strains were clonally identical. A missense mutation (P232L) and a nonsense mutation (Q461Stop) were identified in the katG of the two post-treatment strains, respectively. Transformation experiments showed that katG of the pre-treatment strain restored INH susceptibility in GA03, whereas the mutated katG genes from the post-treatment strains rendered negative catalase activity and INH resistance. The protein model indicated that P232L reduced INH-KatG binding affinity while Q461Stop truncated gene transcription. Our results showed that the two katG mutations, P232L and Q461Stop, accounted for the co-emergence of INH-resistant clones during anti-TB therapy. The inclusion of these mutations in the design of molecular assays could increase the diagnostic performance.IMPORTANCEThe evolution of drug-resistant strains of Mycobacterium tuberculosis within the lung lesions of a patient has a detrimental impact on treatment outcomes. This is particularly concerning for isoniazid (INH), which is the most potent first-line antimycobacterial drug. However, the precise genetic factors responsible for drug resistance in patients have not been fully elucidated, with approximately 15% of INH-resistant strains harboring unknown genetic factors. This raises concerns about the emergence of drug-resistant clones within patients, further contributing to the global epidemic of resistance. In this study, we revealed the presence of two novel katG mutations, which emerged independently due to the stress exerted by antituberculosis (anti-TB) treatment on a parental strain. Importantly, we experimentally demonstrated the functional significance of both mutations in conferring resistance to INH. Overall, this research sheds light on the genetic mechanisms underlying the evolution of INH resistance within patients and provides valuable insights for improving diagnostic performance by targeting specific mutations.

Mathematical models of drug-resistant tuberculosis lack bacterial heterogeneity: A systematic review

Drug-resistant tuberculosis (DR-TB) threatens progress in the control of TB. Mathematical models are increasingly being used to guide public health decisions on managing both antimicrobial resistance (AMR) and TB. It is important to consider bacterial heterogeneity in models as it can have consequences for predictions of resistance prevalence, which may affect decision-making. We conducted a systematic review of published mathematical models to determine the modelling landscape and to explore methods for including bacterial heterogeneity. Our first objective was to identify and analyse the general characteristics of mathematical models of DR-mycobacteria, including M. tuberculosis. The second objective was to analyse methods of including bacterial heterogeneity in these models. We had different definitions of heterogeneity depending on the model level. For between-host models of mycobacterium, heterogeneity was defined as any model where bacteria of the same resistance level were further differentiated. For bacterial population models, heterogeneity was defined as having multiple distinct resistant populations. The search was conducted following PRISMA guidelines in five databases, with studies included if they were mechanistic or simulation models of DR-mycobacteria. We identified 195 studies modelling DR-mycobacteria, with most being dynamic transmission models of non-treatment intervention impact in M. tuberculosis (n = 58). Studies were set in a limited number of specific countries, and 44% of models (n = 85) included only a single level of "multidrug-resistance (MDR)". Only 23 models (8 between-host) included any bacterial heterogeneity. Most of these also captured multiple antibiotic-resistant classes (n = 17), but six models included heterogeneity in bacterial populations resistant to a single antibiotic. Heterogeneity was usually represented by different fitness values for bacteria resistant to the same antibiotic (61%, n = 14). A large and growing body of mathematical models of DR-mycobacterium is being used to explore intervention impact to support policy as well as theoretical explorations of resistance dynamics. However, the majority lack bacterial heterogeneity, suggesting that important evolutionary effects may be missed.

Genotype-Phenotype Characterization of Serial Mycobacterium tuberculosis Isolates in Bedaquiline-Resistant Tuberculosis

BACKGROUND

Emerging resistance to bedaquiline (BDQ) threatens to undermine advances in the treatment of drug-resistant tuberculosis (DRTB). Characterizing serial Mycobacterium tuberculosis (Mtb) isolates collected during BDQ-based treatment can provide insights into the etiologies of BDQ resistance in this important group of DRTB patients.

METHODS

We measured mycobacteria growth indicator tube (MGIT)-based BDQ minimum inhibitory concentrations (MICs) of Mtb isolates collected from 195 individuals with no prior BDQ exposure who were receiving BDQ-based treatment for DRTB. We conducted whole-genome sequencing on serial Mtb isolates from all participants who had any isolate with a BDQ MIC >1 collected before or after starting treatment (95 total Mtb isolates from 24 participants).

RESULTS

Sixteen of 24 participants had BDQ-resistant TB (MGIT MIC ≥4 µg/mL) and 8 had BDQ-intermediate infections (MGIT MIC = 2 µg/mL). Participants with pre-existing resistance outnumbered those with resistance acquired during treatment, and 8 of 24 participants had polyclonal infections. BDQ resistance was observed across multiple Mtb strain types and involved a diverse catalog of mmpR5 (Rv0678) mutations, but no mutations in atpE or pepQ. Nine pairs of participants shared genetically similar isolates separated by <5 single nucleotide polymorphisms, concerning for potential transmitted BDQ resistance.

CONCLUSIONS

BDQ-resistant TB can arise via multiple, overlapping processes, including transmission of strains with pre-existing resistance. Capturing the within-host diversity of these infections could potentially improve clinical diagnosis, population-level surveillance, and molecular diagnostic test development.

The Role of Efflux Pumps transporter in Multi-drug Resistant Tuberculosis: Mycobacterial memberane protein(MmpL5)

BACKGROUND

The overexpression of efflux pumps (Eps) was reported to contribute to multidrug resistant tuberculosis (MDR-TB). Increases in Eps that expel structurally unrelated drugs contribute to reduced susceptibility by decreasing the intracellular concentration of antibiotics. In the present study, an association of mycobacterial membrane protein (MmpS5-MmpL5) Ep and its gene regulator (Rv0678) was investigated in MDR-tuberculosis isolates.

METHODS

MTB strains were isolated from patients at two different intervals, i.e., once when they had persistent symptoms despite 3-15 ≥ months of treatment and once when they had started new combination therapy ≥2-3 months. Sputum specimens were subjected to Xpert MTB/rifampicin test and then further susceptibility testing using proportional method and multiplex polymerase chain reaction (PCR) were performed on them. The isolates were characterized using both 16S-23S RNA and hsp65 genes spacer (PCR-restriction fragment length polymorphism). Whole-genome sequencing (WGS) was investigated on two isolates from culture-positive specimen per patient. The protein structure was simulated using the SWISS-MODEL. The input format used for this web server was FASTA (amino acid sequence). Protein structure was also analysis using Ramachandran plot.

RESULTS

WGS documented deletion, insertion, and substitution in transmembrane transport protein MmpL5 (Rv0676) of Eps. Majority of the studied isolates (n = 12; 92.3%) showed a unique deletion mutation at three positions: (a) from amino acid number 771 (isoleucine) to 776 (valine), (b) from amino acid number 785 (valine) to 793 (histidine), and (c) from amino acid number 798 (leucine) to 806 (glycine)." One isolate (7.6%) had no deletion mutation. In all isolates (n = 13; 100%), a large insertion mutation consisting of 94 amino acid was observed "from amino acid number 846 (isoleucine) to amino acid number 939 (leucine)". Thirty-eight substitutions in Rv0676 were detected, of which 92.3% were identical in the studied isolates. WGS of mycobacterial membrane proteins (MmpS5; Rv0677) and its gene regulator (Rv0678) documented no deletion, insertion, and substitution. No differences were observed between MmpS5-MmpL5 and its gene regulator in isolates that were collected at different intervals.

CONCLUSIONS

Significant genetic mutation like insertion, deletion, and substitution within transmembrane transport protein MmpL5 (Rv0676) can change the functional balance of Eps and cause a reduction in drug susceptibility. This is the first report documenting a unique amino acid mutation (insertion and deletion ≥4-94) in Rv0676 among drug-resistant MTB. We suggest the changes in Mmpl5 (Rv0676) might occurred due to in-vivo sub-therapeutic drug stress within the host cell. Changes in MmpL5 are stable and detected through subsequent culture-positive specimens.

Micro-nanoemulsion and nanoparticle-assisted drug delivery against drug-resistant tuberculosis: recent developments

Tuberculosis (TB) is a major global health problem and the second most prevalent infectious killer after COVID-19. It is caused by Mycobacterium tuberculosis (Mtb) and has become increasingly challenging to treat due to drug resistance. The World Health Organization declared TB a global health emergency in 1993. Drug resistance in TB is driven by mutations in the bacterial genome that can be influenced by prolonged drug exposure and poor patient adherence. The development of drug-resistant forms of TB, such as multidrug resistant, extensively drug resistant, and totally drug resistant, poses significant therapeutic challenges. Researchers are exploring new drugs and novel drug delivery systems, such as nanotechnology-based therapies, to combat drug resistance. Nanodrug delivery offers targeted and precise drug delivery, improves treatment efficacy, and reduces adverse effects. Along with nanoscale drug delivery, a new generation of antibiotics with potent therapeutic efficacy, drug repurposing, and new treatment regimens (combinations) that can tackle the problem of drug resistance in a shorter duration could be promising therapies in clinical settings. However, the clinical translation of nanomedicines faces challenges such as safety, large-scale production, regulatory frameworks, and intellectual property issues. In this review, we present the current status, most recent findings, challenges, and limiting barriers to the use of emulsions and nanoparticles against drug-resistant TB.

Native microbiome dominates over host factors in shaping the probiotic genetic evolution in the gut

Probiotics often acquire potentially adaptive mutations in vivo, gaining new functional traits through gut selection. While both the host and microbiome can contribute to probiotics' genetic evolution, separating the microbiome and the host's contribution to such selective pressures remains challenging. Here, we introduced germ-free (GF) and specific pathogen-free (SPF) mouse models to track how probiotic strains, i.e., Lactiplantibacillus plantarum HNU082 (Lp082) and Bifidobacterium animalis subsp. lactis V9 (BV9), genetically evolved under selection pressures derived from host factors alone and both host and microbial ecological factors. Notably, compared to the genome of a probiotic strain before consumption, the host only elicited <15 probiotic mutations in probiotic genomes that emerged in the luminal environment of GF mice, while a total of 840 mutations in Lp082 mutants and 21,579 mutations in BV9 were found in SPF mice, <0.25% of those derived from both factors that were never captured by other experimental evolution studies, indicating that keen microbial competitions exhibited the predominant evolutionary force in shaping probiotic genetic composition (>99.75%). For a given probiotic, functional genes occurring in potentially adaptive mutations induced by hosts (GF mice) were all shared with those found in mutants of SPF mice. Collectively, the native microbiome consistently drove a more rapid and divergent genetic evolution of probiotic strains in seven days of colonization than host factors did. Our study further laid a theoretical foundation for genetically engineering probiotics for better gut adaptation through in vitro artificial gut ecosystems without the selection pressures derived from host factors.

Overcoming Mycobacterium tuberculosis Drug Resistance: Novel Medications and Repositioning Strategies

Mycobacterium tuberculosis, the bacterium responsible for tuberculosis, is a global health concern, affecting millions worldwide. This bacterium has earned a reputation as a formidable adversary due to its multidrug-resistant nature, allowing it to withstand many antibiotics. The development of this drug resistance in Mycobacterium tuberculosis is attributed to innate and acquired mechanisms. In the past, rifampin was considered a potent medication for treating tuberculosis infections. However, the rapid development of resistance to this drug by the bacterium underscores the pressing need for new therapeutic agents. Fortunately, several other medications previously overlooked for tuberculosis treatment are already available in the market. Moreover, several innovative drugs are under clinical investigation, offering hope for more effective treatments. To enhance the effectiveness of these drugs, it is recommended that researchers concentrate on identifying unique target sites within the bacterium during the drug development process. This strategy could potentially circumvent the issues presented by Mycobacterium drug resistance. This review primarily focuses on the characteristics of novel drug resistance mechanisms in Mycobacterium tuberculosis. It also discusses potential medications being repositioned or sourced from novel origins. The ultimate objective of this review is to discover efficacious treatments for tuberculosis that can successfully tackle the hurdles posed by Mycobacterium drug resistance.

Structural dynamics insights into the M306L, M306V, and D1024N mutations in Mycobacterium tuberculosis inducing resistance to ethambutol

Resistance to anti-tuberculosis drugs, especially ethambutol (EMB), has been widely reported worldwide. EMB resistance is caused by mutations in the embB gene, which encodes the arabinosyl transferase enzyme. This study aimed to detect mutations in the embB gene of Mycobacterium tuberculosis from Papua and to evaluate their impact on the effectiveness of EMB. We analyzed 20 samples of M. tuberculosis culture that had undergone whole-genome sequencing, of which 19 samples were of sufficient quality for further bioinformatics analysis. Mutation analysis was performed using TBProfiler, which identified M306L, M306V, D1024N, and E378A mutations. In sample TB035, the M306L mutation was present along with E378A. The binding affinity of EMB to arabinosyl transferase was calculated using AutoDock Vina. The molecular docking results revealed that all mutants demonstrated an increased binding affinity to EMB compared to the native protein (-0.948 kcal/mol). The presence of the M306L mutation, when coexisting with E378A, resulted in a slight increase in binding affinity compared to the M306L mutation alone. The molecular dynamics simulation results indicated that the M306L, M306L + E378A, M306V, and E378A mutants decreased protein stability. Conversely, the D1024N mutant exhibited stability comparable to the native protein. In conclusion, this study suggests that the M306L, M306L + E378A, M306V, and E378A mutations may contribute to EMB resistance, while the D1024N mutation may be consistent with continued susceptibility to EMB.

Mixed strain pathogen populations accelerate the evolution of antibiotic resistance in patients

Antibiotic resistance poses a global health threat, but the within-host drivers of resistance remain poorly understood. Pathogen populations are often assumed to be clonal within hosts, and resistance is thought to emerge due to selection for de novo variants. Here we show that mixed strain populations are common in the opportunistic pathogen P. aeruginosa. Crucially, resistance evolves rapidly in patients colonized by multiple strains through selection for pre-existing resistant strains. In contrast, resistance evolves sporadically in patients colonized by single strains due to selection for novel resistance mutations. However, strong trade-offs between resistance and growth rate occur in mixed strain populations, suggesting that within-host diversity can also drive the loss of resistance in the absence of antibiotic treatment. In summary, we show that the within-host diversity of pathogen populations plays a key role in shaping the emergence of resistance in response to treatment.

Mycobacterium tuberculosis Intra-Host Evolution Among Drug-Resistant Tuberculosis Patients Failing Treatment

Background

Understanding Mycobacterium tuberculosis (Mtb) intra-host evolution of drug resistance is important for successful drug-resistant tuberculosis (DR-TB) treatment and control strategies. This study aimed to characterise the acquisition of genetic mutations and low-frequency variants associated with treatment-emergent Mtb drug resistance in longitudinally profiled clinical isolates from patients who experienced DR-TB treatment failure.

Patients and Methods

We performed deep Whole Genome Sequencing on 23 clinical isolates obtained longitudinally across nine timepoints from five patients who experienced DR-TB treatment failure enrolled in the CAPRISA 020 InDEX study. The minimum inhibitory concentrations (MICs) were established on the BACTEC™ MGIT 960™ instrument on 15/23 longitudinal clinical isolates for eight anti-TB drugs (rifampicin, isoniazid, ethambutol, levofloxacin, moxifloxacin, linezolid, clofazimine, bedaquiline).

Results

In total, 22 resistance associated mutations/variants were detected. We observed four treatment-emergent mutations in two out of the five patients. Emerging resistance to the fluoroquinolones was associated with 16- and 64-fold elevated levofloxacin (2-8 mg/L) and moxifloxacin (1-2 mg/L) MICs, respectively, resulting from the D94G/N and A90V variants in the gyrA gene. We identified two novel mutations associated with elevated bedaquiline MICs (>66-fold): an emerging frameshift variant (D165) on the Rv0678 gene and R409Q variant on the Rv1979c gene present from baseline.

Conclusion

Genotypic and phenotypic resistance to the fluoroquinolones and bedaquiline was acquired in two out of five patients who experienced DR-TB treatment failure. Deep sequencing of multiple longitudinal clinical isolates for resistance-associated mutations coupled with phenotypic MIC testing confirmed intra-host Mtb evolution.

Intermittent antibiotic treatment of bacterial biofilms favors the rapid evolution of resistance

Bacterial antibiotic resistance is a global health concern of increasing importance and intensive study. Although biofilms are a common source of infections in clinical settings, little is known about the development of antibiotic resistance within biofilms. Here, we use experimental evolution to compare selection of resistance mutations in planktonic and biofilm Escherichia coli populations exposed to clinically relevant cycles of lethal treatment with the aminoglycoside amikacin. Consistently, mutations in sbmA, encoding an inner membrane peptide transporter, and fusA, encoding the essential elongation factor G, are rapidly selected in biofilms, but not in planktonic cells. This is due to a combination of enhanced mutation rate, increased adhesion capacity and protective biofilm-associated tolerance. These results show that the biofilm environment favors rapid evolution of resistance and provide new insights into the dynamic evolution of antibiotic resistance in biofilms.

GWAS and functional studies suggest a role for altered DNA repair in the evolution of drug resistance in Mycobacterium tuberculosis

The emergence of drug resistance in Mycobacterium tuberculosis (Mtb) is alarming and demands in-depth knowledge for timely diagnosis. We performed genome-wide association analysis using 2237 clinical strains of Mtb to identify novel genetic factors that evoke drug resistance. In addition to the known direct targets, we identified for the first time, a strong association between mutations in DNA repair genes and the multidrug-resistant phenotype. To evaluate the impact of variants identified in the clinical samples in the evolution of drug resistance, we utilized knockouts and complemented strains in Mycobacterium smegmatis and Mtb. Results show that variant mutations compromised the functions of MutY and UvrB. MutY variant showed enhanced survival compared with wild-type (Rv) when the Mtb strains were subjected to multiple rounds of ex vivo antibiotic stress. In an in vivo guinea pig infection model, the MutY variant outcompeted the wild-type strain. We show that novel variant mutations in the DNA repair genes collectively compromise their functions and contribute to better survival under antibiotic/host stress conditions.

Macrophage activation highlight an important role for NER proteins in the survival, latency and multiplication of Mycobacterium tuberculosis

Nucleotide excision repair (NER) is one of the most extensively studied DNA repair processes in both prokaryotes and eukaryotes. The NER pathway is a highly conserved, ATP-dependent multi-step process involving several proteins/enzymes that function in a concerted manner to recognize and excise a wide spectrum of helix-distorting DNA lesions and bulky adducts by nuclease cleavage on either side of the damaged bases. As such, the NER pathway of Mycobacterium tuberculosis (Mtb) is essential for its survival within the hostile environment of macrophages and disease progression. This review focuses on present published knowledge about the crucial roles of Mtb NER proteins in the survival and multiplication of the pathogen within the macrophages and as potential targets for drug discovery.

Within-host evolution of bacterial pathogens during persistent infection of humans

Many bacterial pathogens can form persistent infections, providing an infectious reservoir, which allows for infection of new hosts. Currently, the molecular mechanisms and evolutionary dynamics driving persistence are still not well-understood. High-throughput sequencing methods have enabled the study of within-host evolution of persistent bacterial pathogens, revealing common trends among bacterial species in how they adapt to persist. We will focus on trends emerging from longitudinal human-cohort studies, including i) genome-size reduction, ii) metabolic adaptation to the host, iii) antimicrobial resistance, iv) changes in virulence and the bacterial cell surface, and v) hypermutation, and comment on where the field should focus going forward.

Evolution of Mycobacterium tuberculosis drug resistance in the genomic era

Mycobacterium tuberculosis has acquired drug resistance to all drugs that have been used against it, including those only recently introduced into clinical practice. Compared to other bacteria, it has a well conserved genome due to its role as an obligate human pathogen that has adapted to a niche over five to ten thousand years. These features facilitate reconstruction and dating of M. tuberculosis phylogenies, giving key insights into how resistance has been acquired and spread globally. Resistance to each new drug has occurred within five to ten years of clinical use and has occurred even more rapidly with recently introduced drugs. In most cases, resistance-conferring mutations come with a fitness cost, but this can be overcome by compensatory mutations which restore fitness to that of wild-type bacteria. It is likely that M. tuberculosis acquires drug resistance while maintaining limited genomic variability due the generation of low frequency within-host variation, combined with ongoing purifying selection causing loss of variants without a clear fitness advantage. However, variants that do confer an advantage, such as drug resistance, can increase in prevalence amongst all bacteria within a host and become the dominant clone. These resistant strains can then be transmitted leading to primary drug resistant infection in a new host. As many countries move towards genomic methods for diagnosis of M. tuberculosis infection and drug resistance, it is important to be aware of the implications for the evolution of resistance. Currently, understanding of resistance-conferring mutations is incomplete, and some targeted genetic diagnostics create their own selective pressures. We discuss an example where a rifampicin resistance-conferring mutation which was not routinely covered by standard testing became dominant. Finally, resistance to new drugs such as bedaquiline and delamanid is caused by individually rare mutations occurring across a large mutational genomic target that have been detected over a short time, and do not provide statistical power for genotype-phenotype correlation – in contrast to longer-established drugs that form the backbone of drug-sensitive antituberculosis therapy. Therefore, we need a different approach to identify resistance-conferring mutations of new drugs before their resistance becomes widespread, abrogating their usefulness.

Antibiotic resistance acquisition versus primary transmission in the presentation of extensively drug-resistant tuberculosis

Mycobacterium tuberculosis is the leading cause of mortality worldwide due to a single bacterial pathogen. Of concern is the negative impact that the COVID-19 pandemic has had on the control of tuberculosis (TB) including drug-resistant forms of the disease. Antimicrobial resistance increases the likelihood of worsened outcomes in TB patients including treatment failure and death. Multidrug-resistant (MDR) strains, resistant to first-line drugs isoniazid and rifampin, and extensively drug-resistant (XDR) strains with further resistance to second-line drugs (SLD), threaten control programs designed to lower TB incidence and end the disease as a public health challenge by 2030, in accordance with UN Sustainable Development Goals. Tackling TB requires an understanding of the pathways through which drug resistance emerges. Here, the roles of acquired resistance mutation, and primary transmission, are examined with regard to XDR-TB. It is apparent that XDR-TB can emerge from MDR-TB through a small number of additional resistance mutations that occur in patients undergoing drug treatment. Rapid detection of resistance, to first-line drugs and SLD, at the initiation of and during treatment, and prompt adjustment of regimens are required to ensure treatment success in these patients. Primary transmission is predicted to make an increasing contribution to the XDR-TB caseload in the future. Much work is required to improve the implementation of the World Health Organization-recommended infection control practices and block onward transmission of XDR-TB patients to contacts including health-care workers. Finally, limiting background resistance to fluoroquinolones in pre-XDR strains of M. tuberculosis will necessitate better antimicrobial stewardship in the broader use of this drug class.

An Outer Membrane Protein YiaD Contributes to Adaptive Resistance of Meropenem in Acinetobacter baumannii

Acinetobacter baumannii is an important nosocomial pathogen that can develop various resistance mechanisms to many antibiotics. However, little is known about how it evolves from an antibiotic sensitive to a resistant phenotype. In this study, we investigated the transition of outer membrane proteins (OMPs) under antibiotic stress and identified YiaD as an OMP marker involved in the development of adaptive resistance to meropenem (MEM) in A. baumannii. Following stimulation of a carbapenem-sensitive strain AB5116 with sub-MIC of MEM, yiaD showed significantly decreased expression, and this decrease continued with prolonged stimulation for 8 h. The downregulation of yiaD was not only observed in clinically sensitive strains but also in 45 carbapenem-resistant isolates that produced the β-lactamases TEM and OXA-23. However, the extent of the reduction of yiaD expression in resistant strains was less than that in sensitive strains. Lack of yiaD resulted in a 4-fold increase in the MIC of AB5116 to MEM. The same level of depressed susceptibility induced by yiaD deletion was observed in both a growth curve test and a survival rate assay. Moreover, the colony shape became enlarged and irregular after loss of yiaD, and the biofilm formation ability of A. baumannii was influenced by YiaD. These results suggest that YiaD could respond to the stimulus of MEM in A. baumannii with a downregulation trend that kept pace with the prolonged stimulation time, indicating that it participates in various routes to benefit MEM resistance evolution in both carbapenem-sensitive and -resistant A. baumannii strains.

IMPORTANCEAcinetobacter baumannii can develop various resistance mechanisms to carbapenems. However, the factors involved in the evolutionary process that leads from transition to the sensitive to resistant phenotype are not clear. The outer membrane protein YiaD of A. baumannii was downregulated under the stress of meropenem (MEM), and its expression level was continuously reduced with prolonged stimulation time. The downregulation of yiaD was not only observed in sensitive strains but also in carbapenem-resistant isolates producing the β-lactamases TEM and OXA-23. However, the extent of yiaD reduction was less in resistant strains than in sensitive strains. Lack of yiaD resulted in an increased MEM MIC, enlarged and irregular colonies, and decreased biofilm formation ability. These results suggest that YiaD responds to MEM stimulus in A. baumannii and participates in the adaptive resistance of MEM in both carbapenem-sensitive and -resistant strains.

Contribution and Future of High-Throughput Transcriptomics in Battling Tuberculosis

While Tuberculosis (TB) infection remains a serious challenge worldwide, big data and “omic” approaches have greatly contributed to the understanding of the disease. Transcriptomics have been used to tackle a wide variety of queries including diagnosis, treatment evolution, latency and reactivation, novel target discovery, vaccine response or biomarkers of protection. Although a powerful tool, the elevated cost and difficulties in data interpretation may hinder transcriptomics complete potential. Technology evolution and collaborative efforts among multidisciplinary groups might be key in its exploitation. Here, we discuss the main fields explored in TB using transcriptomics, and identify the challenges that need to be addressed for a real implementation in TB diagnosis, prevention and therapy.

Patterns of genomic interrelatedness of publicly available samples in the TB portals database

The TB Portals program is an international collaboration for the collection and dissemination of tuberculosis data from patient cases focused on drug resistance. The central database is a patient-oriented resource containing both patient and pathogen clinical and genomic information. Herein we provide a summary of the pathogen genomic data available through the TB Portals and show one potential application by examining patterns of genomic pairwise distances. Distributions of pairwise distances highlight overall patterns of genome variability within and between Mycobacterium tuberculosis phylogenomic lineages. Closely related isolates (based on whole-genome pairwise distances and time between sample collection dates) from different countries were identified as potential evidence of international transmission of drug-resistant tuberculosis. These high-level views of genomic relatedness provide information that can stimulate hypotheses for further and more detailed research.

Genome reorganization during emergence of host-associated Mycobacterium abscessus

Mycobacterium abscessus is a rapid growing, free-living species of bacterium that also causes lung infections in humans. Human infections are usually acquired from the environment; however, dominant circulating clones (DCCs) have emerged recently in both M. abscessus subsp. massiliense and subsp. abscessus that appear to be transmitted among humans and are now globally distributed. These recently emerged clones are potentially informative about the ecological and evolutionary mechanisms of pathogen emergence and host adaptation. The geographical distribution of DCCs has been reported, but the genomic processes underlying their transition from environmental bacterium to human pathogen are not well characterized. To address this knowledge gap, we delineated the structure of M. abscessus subspecies abscessus and massiliense using genomic data from 200 clinical isolates of M. abscessus from seven geographical regions. We identified differences in overall patterns of lateral gene transfer (LGT) and barriers to LGT between subspecies and between environmental and host-adapted bacteria. We further characterized genome reorganization that accompanied bacterial host adaptation, inferring selection pressures acting at both genic and intergenic loci. We found that both subspecies encode an expansive pangenome with many genes at rare frequencies. Recombination appears more frequent in M. abscessus subsp. massiliense than in subsp. abscessus, consistent with prior reports. We found evidence suggesting that phage are exchanged between subspecies, despite genetic barriers evident elsewhere throughout the genome. Patterns of LGT differed according to niche, with less LGT observed among host-adapted DCCs versus environmental bacteria. We also found evidence suggesting that DCCs are under distinct selection pressures at both genic and intergenic sites. Our results indicate that host adaptation of M. abscessus was accompanied by major changes in genome evolution, including shifts in the apparent frequency of LGT and impacts of selection. Differences were evident among the DCCs as well, which varied in the degree of gene content remodelling, suggesting they were placed differently along the evolutionary trajectory toward host adaptation. These results provide insight into the evolutionary forces that reshape bacterial genomes as they emerge into the pathogenic niche.

Advantages of laboratory natural selection in the applied sciences

In the past three decades, laboratory natural selection has become a widely used technique in biological research. Most studies which have utilized this technique are in the realm of basic science, often testing hypotheses related to mechanisms of evolutionary change or ecological dynamics. While laboratory natural selection is currently utilized heavily in this setting, there is a significant gap with its usage in applied studies, especially when compared to the other selection experiment methodologies like artificial selection and directed evolution. This is despite avenues of research in the applied sciences which seem well suited to laboratory natural selection. In this review, we place laboratory natural selection in context with other selection experiments, identify the characteristics which make it well suited for particular kinds of applied research and briefly cover key examples of the usefulness of selection experiments within applied science. Finally, we identify three promising areas of inquiry for laboratory natural selection in the applied sciences: bioremediation technology, identifying mechanisms of drug resistance and optimizing biofuel production. Although laboratory natural selection is currently less utilized in applied science when compared to basic research, the method has immense promise in the field moving forward.

Genetic Profile of Linezolid-Resistant M. tuberculosis Clinical Strains from Moscow

Background: Linezolid, bedaquiline, and newer fluoroquinolones are currently placed as priority Group A drugs for the treatment of drug-resistant tuberculosis. The number of reported linezolid-resistant clinical strains is still low, and the correlation of molecular determinants with phenotype is not perfect. Methods: We determined the linezolid MICs for clinical isolates from the Moscow region and identified mutations in rplC and rrl genes. Results: All 16 linezolid-resistant isolates had previously reported mutations in the rplC or rrl loci, and 13 of them bore a RplC C154R substitution. Detection of this substitution in a heteroresistant state was not successful, probably, due to the more stable DNA secondary structure of the mutated fragment, which precludes its amplification in mixes with the wild-type DNA. Strains with an rplC mutation had higher linezolid MIC compared to isolates with rrl mutations. Conclusions: Linezolid resistance mostly emerged during treatment with the latest regimen. Three primary cases with linezolid resistance question the possible transmission of totally drug-resistant tuberculosis in the Moscow region, which demands further investigation.

The within-host evolution of antimicrobial resistance in Mycobacterium tuberculosis

Tuberculosis (TB) has been responsible for the greatest number of human deaths due to an infectious disease in general, and due to antimicrobial resistance (AMR) in particular. The etiological agents of human TB are a closely-related group of human-adapted bacteria that belong to the Mycobacterium tuberculosis complex (MTBC). Understanding how MTBC populations evolve within-host may allow for improved TB treatment and control strategies. In this review, we highlight recent works that have shed light on how AMR evolves in MTBC populations within individual patients. We discuss the role of heteroresistance in AMR evolution, and review the bacterial, patient and environmental factors that likely modulate the magnitude of heteroresistance within-host. We further highlight recent works on the dynamics of MTBC genetic diversity within-host, and discuss how spatial substructures in patients' lungs, spatiotemporal heterogeneity in antimicrobial concentrations and phenotypic drug tolerance likely modulates the dynamics of MTBC genetic diversity in patients during treatment. We note the general characteristics that are shared between how the MTBC and other bacterial pathogens evolve in humans, and highlight the characteristics unique to the MTBC.

Common Adaptive Strategies Underlie Within-Host Evolution of Bacterial Pathogens

Within-host adaptation is a hallmark of chronic bacterial infections, involving substantial genomic changes. Recent large-scale genomic data from prolonged infections allow the examination of adaptive strategies employed by different pathogens and open the door to investigate whether they converge toward similar strategies. Here, we compiled extensive data of whole-genome sequences of bacterial isolates belonging to miscellaneous species sampled at sequential time points during clinical infections. Analysis of these data revealed that different species share some common adaptive strategies, achieved by mutating various genes. Although the same genes were often mutated in several strains within a species, different genes related to the same pathway, structure, or function were changed in other species utilizing the same adaptive strategy (e.g., mutating flagellar genes). Strategies exploited by various bacterial species were often predicted to be driven by the host immune system, a powerful selective pressure that is not species specific. Remarkably, we find adaptive strategies identified previously within single species to be ubiquitous. Two striking examples are shifts from siderophore-based to heme-based iron scavenging (previously shown for Pseudomonas aeruginosa) and changes in glycerol-phosphate metabolism (previously shown to decrease sensitivity to antibiotics in Mycobacterium tuberculosis). Virulence factors were often adaptively affected in different species, indicating shifts from acute to chronic virulence and virulence attenuation during infection. Our study presents a global view on common within-host adaptive strategies employed by different bacterial species and provides a rich resource for further studying these processes.

Microevolution of Mycobacterium tuberculosis Subpopulations and Heteroresistance in a Patient Receiving 27 Years of Tuberculosis Treatment in Germany

Preexisting and newly emerging resistant pathogen subpopulations (heteroresistance) are potential risk factors for treatment failure of multi/extensively drug resistant (MDR/XDR) tuberculosis (TB). Intrapatient evolutionary dynamics of Mycobacterium tuberculosis complex (Mtbc) strains and their implications on treatment outcomes are still not completely understood. To elucidate how Mtbc strains escape therapy, we analyzed 13 serial isolates from a German patient by whole-genome sequencing. Sequencing data were compared with phenotypic drug susceptibility profiles and the patient's collective 27-year treatment history to further elucidate factors fostering intrapatient resistance evolution. The patient endured five distinct TB episodes, ending in resistance to 16 drugs and a nearly untreatable XDR-TB infection. The first isolate obtained, during the patient's 5th TB episode, presented fixed resistance mutations to 7 anti-TB drugs, including isoniazid, rifampin, streptomycin, pyrazinamide, prothionamide, para-aminosalicylic acid, and cycloserine-terizidone. Over the next 13 years, a dynamic evolution with coexisting, heterogeneous subpopulations was observed in 6 out of 13 sequential bacterial isolates. The emergence of drug-resistant subpopulations coincided with frequent changes in treatment regimens, which often included two or fewer active compounds. This evolutionary arms race between competing subpopulations ultimately resulted in the fixation of a single XDR variant. Our data demonstrate the complex intrapatient microevolution of Mtbc subpopulations during failing MDR/XDR-TB treatment. Designing effective treatment regimens based on rapid detection of (hetero) resistance is key to avoid resistance development and treatment failure.

Import and transmission of Mycobacterium orygis and Mycobacterium africanum, Norway

Background

The aim of the current study was to improve our understanding of the origins and transmission of Mycobacterium africanum (MAF) in Norway.

Methods

Whole-genome sequences (WGS) were generated for all (n = 29) available clinical isolates received at the Norwegian National Reference Laboratory for Mycobacteria (NRL) and identified as MAF in Norway, in the period 2010–2020. Phylogenetic analyses were performed.

Results

The analyses indicated several imports of MAF lineage 6 from both East and West African countries, whereas MAF lineage 5 was restricted to patients with West African connections. We also find evidence for transmission of MAF in Norway. Finally, our analyses revealed that a group of isolates from patients originating in South Asia, identified as MAF by means of a commercial line-probe assay, in fact belonged to Mycobacterium orygis.

Conclusions

Most MAF cases in Norway are the result of import, but transmission is occurring within Norway.

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.

Collateral Sensitivity to β-Lactam Drugs in Drug-Resistant Tuberculosis Is Driven by the Transcriptional Wiring of BlaI Operon Genes

The evolution of resistance to one antimicrobial can result in enhanced sensitivity to another, known as "collateral sensitivity." This underexplored phenomenon opens new therapeutic possibilities for patients infected with pathogens unresponsive to classical treatments. Intrinsic resistance to β-lactams in Mycobacterium tuberculosis (the causative agent of tuberculosis) has traditionally curtailed the use of these low-cost and easy-to-administer drugs for tuberculosis treatment. Recently, β-lactam sensitivity has been reported in strains resistant to classical tuberculosis therapy, resurging the interest in β-lactams for tuberculosis. However, a lack of understanding of the molecular underpinnings of this sensitivity has delayed exploration in the clinic. We performed gene expression and network analyses and in silico knockout simulations of genes associated with β-lactam sensitivity and genes associated with resistance to classical tuberculosis drugs to investigate regulatory interactions and identify key gene mediators. We found activation of the key inhibitor of β-lactam resistance, blaI, following classical drug treatment as well as transcriptional links between genes associated with β-lactam sensitivity and those associated with resistance to classical treatment, suggesting that regulatory links might explain collateral sensitivity to β-lactams. Our results support M. tuberculosis β-lactam sensitivity as a collateral consequence of the evolution of resistance to classical tuberculosis drugs, mediated through changes to transcriptional regulation. These findings support continued exploration of β-lactams for the treatment of patients infected with tuberculosis strains resistant to classical therapies. IMPORTANCE Tuberculosis remains a significant cause of global mortality, with strains resistant to classical drug treatment considered a major health concern by the World Health Organization. Challenging treatment regimens and difficulty accessing drugs in low-income communities have led to a high prevalence of strains resistant to multiple drugs, making the development of alternative therapies a priority. Although Mycobacterium tuberculosis is naturally resistant to β-lactam drugs, previous studies have shown sensitivity in strains resistant to classical drug treatment, but we currently lack understanding of the molecular underpinnings behind this phenomenon. We found that genes involved in β-lactam susceptibility are activated after classical drug treatment resulting from tight regulatory links with genes involved in drug resistance. Our study supports the hypothesis that β-lactam susceptibility observed in drug-resistant strains results from the underlying regulatory network of M. tuberculosis, supporting further exploration of the use of β-lactams for tuberculosis treatment.

Control of Genome Stability by EndoMS/NucS-Mediated Non-Canonical Mismatch Repair

The DNA repair endonuclease EndoMS/NucS is highly conserved in Archaea and Actinobacteria. This enzyme is able to recognize and cleave dsDNA carrying a mismatched base pair, and its activity is enhanced by the interaction with the sliding clamp of the replisome. Today, EndoMS/NucS has been established as the key protein of a non-canonical mismatch repair (MMR) pathway, acting specifically in the repair of transitions and being essential for maintaining genome stability. Despite having some particularities, such as its lower activity on transversions and the inability to correct indels, EndoMS/NucS meets the main hallmarks of a MMR. Its absence leads to a hypermutator phenotype, a transition-biased mutational spectrum and an increase in homeologous recombination. Interestingly, polymorphic EndoMS/NucS variants with a possible effect in mutation rate have been detected in clinical isolates of the relevant actinobacterial pathogen Mycobacterium tuberculosis. Considering that MMR defects are often associated with the emergence of resistant bacteria, the existence of EndoMS/NucS-defective mutators could have an important role in the acquisition of antibiotic resistance in M. tuberculosis. Therefore, a further understanding of the EndoMS/NucS-mediated non-canonical MMR pathway may reveal new strategies to predict and fight drug resistance. This review is focused on the recent progress in NucS, with special emphasis on its effect on genome stability and evolvability in Actinobacteria.

Five-year microevolution of a multidrug-resistant Mycobacterium tuberculosis strain within a patient with inadequate compliance to treatment

BACKGROUND

Whole-genome sequencing has shown that the Mycobacterium tuberculosis infection process can be more heterogeneous than previously thought. Compartmentalized infections, exogenous reinfections, and microevolution are manifestations of this clonal complexity. The analysis of the mechanisms causing the microevolution -the genetic variability of M. tuberculosis at short time scales- of a parental strain into clonal variants with a patient is a relevant issue that has not been yet completely addressed. To our knowledge, a whole genome sequence microevolution analysis in a single patient with inadequate adherence to treatment has not been previously reported.

CASE PRESENTATION

In this work, we applied whole genome sequencing analysis for a more in-depth analysis of the microevolution of a parental Mycobacterium tuberculosis strain into clonal variants within a patient with poor treatment compliance in Argentina. We analyzed the whole-genome sequence of 8 consecutive Mycobacterium tuberculosis isolates obtained from a patient within 57-months of intermittent therapy. Nineteen mutations (9 short-term, 10 fixed variants) emerged, most of them associated with drug resistance. The first isolate was already resistant to isoniazid, rifampicin, and streptomycin, thereafter the strain developed resistance to fluoroquinolones and pyrazinamide. Surprisingly, isolates remained susceptible to the pro-drug ethionamide after acquiring a frameshift mutation in ethA, a gene required for its activation. We also found a novel variant, (T-54G), in the 5' untranslated region of whiB7 (T-54G), a region allegedly related to kanamycin resistance. Notably, discrepancies between canonical and phage-based susceptibility testing to kanamycin were previously found for the isolate harboring this mutation. In our patient, microevolution was mainly driven by drug selective pressure. Rare short-term mutations fixed together with resistance-conferring mutations during therapy.

CONCLUSIONS

This report highlights the relevance of whole-genome sequencing analysis in the clinic for characterization of pre-XDR and MDR resistance profile, particularly in patients with incomplete and/or intermittent treatment.

Whole genome sequencing based differentiation between re-infection and relapse in Indian patients with tuberculosis recurrence, with and without HIV co-infection

INTRODUCTION

Differentiation between relapse and reinfection in cases with tuberculosis (TB) recurrence has important implications for public health, especially in patients with human immunodeficiency virus (HIV) co-infection. We compared Mycobacterial Interspersed Repeat Unit (MIRU) typing and spoligotyping with whole genome sequencing (WGS) to differentiate between relapse and reinfection in patients (HIV-positive and HIV-negative) with TB recurrence. We also assessed the value of WGS to track acquired drug resistance in those with relapse after successful treatment.

METHOD

Forty-one paired M. tuberculosis isolates collected from 20 HIV-positive and 21 HIV-negative patients were subjected to WGS in addition to spoligotyping and MIRU typing. Phylogenetic and Single Nucleotide Substitution (SNP) clustering analyses were performed to determine whether recurrences were due to relapse or re-infection.

RESULTS

Comparison of M. tuberculosis genomes indicated that 95% of TB recurrences in the HIV-negative cohort were due to relapse, while the majority of TB recurrences (75%) in the HIV-positive cohort was due to reinfection (P = 0.0001). New drug resistance mutations were acquired in 5/24 cases (20.8%) that experienced relapse.

CONCLUSIONS

WGS provided increased resolution, but differentiation between relapse and reinfection was broadly consistent with MIRU and spoligotyping. The high contribution of reinfection among HIV infected patients experiencing TB recurrence warrants further study to explore risk factors for TB exposure.

Lesion Heterogeneity Coincides With Long-Term Heteroresistance in MDR-TB

Tuberculosis (TB) heteroresistance, in which only a fraction of the bacteria in a TB patient contains drug-resistant mutations, has been a rising concern. However, its origins and prevalence remain elusive. Here, whole-genome sequencing was performed on 83 serial isolates from 31 MDR-TB patients and heteroresistance was detected in isolates from 21 (67.74%) patients. Heteroresistance persisted in the host for long periods, spanning months to years, and was associated with having multiple tubercular lesions. Our findings indicate that heteroresistance is common and persistent in MDR-TB patients and may affect the success of their treatment regimens.

Nonsynonymous Polymorphism Counts in Bacterial Genomes: a Comparative Examination

Genomic data reveal single-nucleotide polymorphisms (SNPs) that may carry information about the evolutionary history of bacteria. However, it remains unclear what inferences about selection can be made from genomic SNP data. Bacterial species are often sampled during epidemic outbreaks or within hosts during the course of chronic infections. SNPs obtained from genomic analysis of these data are not necessarily fixed. Treating them as fixed during analysis by using measures such as the ratio of nonsynonymous to synonymous evolutionary changes (dN/dS) may lead to incorrect inferences about the strength and direction of selection. In this study, we consider data from a range of whole-genome sequencing studies of bacterial pathogens and explore patterns of nonsynonymous variation to assess whether evidence of selection can be identified by investigating SNP counts alone across multiple WGS studies. We visualize these SNP data in ways that highlight their relationship to neutral baseline expectations. These neutral expectations are based on a simple model of mutation, from which we simulate SNP accumulation to investigate how SNP counts are distributed under alternative assumptions about positive and negative selection. We compare these patterns with empirical SNP data and illustrate the general difficulty of detecting positive selection from SNP data. Finally, we consider whether SNP counts observed at the between-host population level differ from those observed at the within-host level and find some evidence that suggests that dynamics across these two scales are driven by different underlying processes.IMPORTANCE Identifying selection from SNP data obtained from whole-genome sequencing studies is challenging. Some current measures used to identify and quantify selection acting on genomes rely on fixed differences; thus, these are inappropriate for SNP data where variants are not fixed. With the increase in whole-genome sequencing studies, it is important to consider SNP data in the context of evolutionary processes. How SNPs are counted and analyzed can help in understanding mutation accumulation and trajectories of strains. We developed a tool for identifying possible evidence of selection and for comparative analysis with other SNP data. We propose a model that provides a rule-of-thumb guideline and two new visualization techniques that can be used to interpret and compare SNP data. We quantify the expected proportion of nonsynonymous SNPs in coding regions under neutrality and demonstrate its use in identifying evidence of positive and negative selection from simulations and empirical data.

Mechanisms of Drug-Induced Tolerance in Mycobacterium tuberculosis

Successful treatment of tuberculosis (TB) can be hampered by Mycobacterium tuberculosis populations that are temporarily able to survive antibiotic pressure in the absence of drug resistance-conferring mutations, a phenomenon termed drug tolerance. We summarize findings on M. tuberculosis tolerance published in the past 20 years. Key M. tuberculosis responses to drug pressure are reduced growth rates, metabolic shifting, and the promotion of efflux pump activity. Metabolic shifts upon drug pressure mainly occur in M. tuberculosis's lipid metabolism and redox homeostasis, with reduced tricarboxylic acid cycle activity in favor of lipid anabolism. Increased lipid anabolism plays a role in cell wall thickening, which reduces sensitivity to most TB drugs. In addition to these general mechanisms, drug-specific mechanisms have been described. Upon isoniazid exposure, M. tuberculosis reprograms several pathways associated with mycolic acid biosynthesis. Upon rifampicin exposure, M. tuberculosis upregulates the expression of its drug target rpoB Upon bedaquiline exposure, ATP synthesis is stimulated, and the transcription factors Rv0324 and Rv0880 are activated. A better understanding of M. tuberculosis's responses to drug pressure will be important for the development of novel agents that prevent the development of drug tolerance following treatment initiation. Such agents could then contribute to novel TB treatment-shortening strategies.

Dynamics of Extensive Drug Resistance Evolution of Mycobacterium tuberculosis in a Single Patient During 9 Years of Disease and Treatment

Mycobacterium tuberculosis is one of the hardest to treat bacterial pathogens with a high capacity to develop antibiotic resistance by mutations. Here we have performed whole-genome sequencing of consecutive M. tuberculosis isolates obtained during 9 years from a patient with pulmonary tuberculosis. The infecting strain was isoniazid resistant and during treatment it stepwise accumulated resistance mutations to 8 additional antibiotics. Heteroresistance was common and subpopulations with up to 3 different resistance mutations to the same drug coexisted. Sweeps of different resistant clones dominated the population at different time points, always coupled to resistance mutations coinciding with changes in the treatment regimens. Resistance mutations were predominant and no hitch-hiking, compensatory, or virulence-increasing mutations were detected, showing that the dominant selection pressure was antibiotic treatment. The results highlight the dynamic nature of M. tuberculosis infection, population structure, and resistance evolution and the importance of rapid antibiotic susceptibility tests to battle this pathogen.

Insertion and deletion evolution reflects antibiotics selection pressure in a Mycobacterium tuberculosis outbreak

In genome evolution, genetic variants are the source of diversity, which natural selection acts upon. Treatment of human tuberculosis (TB) induces a strong selection pressure for the emergence of antibiotic resistance-conferring variants in the infecting Mycobacterium tuberculosis (MTB) strains. MTB evolution in response to treatment has been intensively studied and mainly attributed to point substitutions. However, the frequency and contribution of insertions and deletions (indels) to MTB genome evolution remains poorly understood. Here, we analyzed a multi-drug resistant MTB outbreak for the presence of high-quality indels and substitutions. We find that indels are significantly enriched in genes conferring antibiotic resistance. Furthermore, we show that indels are inherited during the outbreak and follow a molecular clock with an evolutionary rate of 5.37e-9 indels/site/year, which is 23 times lower than the substitution rate. Inherited indels may co-occur with substitutions in genes along related biological pathways; examples are iron storage and resistance to second-line antibiotics. This suggests that epistatic interactions between indels and substitutions affect antibiotic resistance and compensatory evolution in MTB.

Factors favouring the evolution of multidrug resistance in bacteria

The evolution of multidrug antibiotic resistance in commensal bacteria is an important public health concern. Commensal bacteria such as Escherichia coli, Streptococcus pneumoniae or Staphylococcus aureus, are also opportunistic pathogens causing a large fraction of the community-acquired and hospital-acquired bacterial infections. Multidrug resistance (MDR) makes these infections harder to treat with antibiotics and may thus cause substantial additional morbidity and mortality. Here, we develop an evolutionary epidemiology model to identify the factors favouring the evolution of MDR in commensal bacteria. The model describes the evolution of antibiotic resistance in a commensal bacterial species evolving in a host population subjected to multiple antibiotic treatments. We combine statistical analysis of a large number of simulations and mathematical analysis to understand the model behaviour. We find that MDR evolves more readily when it is less costly than expected from the combinations of single resistances (positive epistasis). MDR frequently evolves when bacteria are in contact with multiple drugs prescribed in the host population, even if individual hosts are only treated with a single drug at a time. MDR is favoured when the host population is structured in different classes that vary in their rates of antibiotic treatment. However, under most circumstances, recombination between loci involved in resistance does not meaningfully affect the equilibrium frequency of MDR. Together, these results suggest that MDR is a frequent evolutionary outcome in commensal bacteria that encounter the variety of antibiotics prescribed in the host population. A better characterization of the variability in antibiotic use across the host population (e.g. across age classes or geographical location) would help predict which MDR genotypes will most readily evolve.

Microevolution of Campylobacter jejuni during long-term infection in an immunocompromised host

Campylobacteriosis typically manifests as a short-lived, self-limiting gastrointestinal infection in humans, however prolonged infection can be seen in cases with underlying immunodeficiency. Public Health England received 25 isolates of Campylobacter jejuni from an individual with combined variable immunodeficiency over a period of 15 years. All isolates were typed and archived at the time of receipt. Whole genome sequencing (WGS) and antimicrobial susceptibility testing were performed to examine the relatedness of the isolates and to investigate the changes in the genome that had taken place over the course of the infection. Genomic typing methods were compared to conventional phenotypic methods, and revealed that the infection was caused by a single, persistent strain of C. jejuni belonging to clonal complex ST-45, with evidence of adaptation and selection in the genome over the course of the infection. Genomic analysis of sequence variants associated with antimicrobial resistance identified the genetic background behind rRNA gene mutations causing variable levels of resistance to erythromycin. This application of WGS to examine a persistent case of campylobacteriosis provides insight into the mutations and selective pressures occurring over the course of an infection, some of which have important clinical relevance.

Mycobacterium tuberculosis clinical isolates carry mutational signatures of host immune environments

Mycobacterium tuberculosis (Mtb) infection results in a spectrum of clinical and histopathologic manifestations. It has been proposed that the environmental and immune pressures associated with different contexts of infection have different consequences for the associated bacterial populations, affecting drug susceptibility and the emergence of resistance. However, there is little concrete evidence for this model. We prospectively collected sputum samples from 18 newly diagnosed and treatment-naïve patients with tuberculosis and sequenced 795 colony-derived Mtb isolates. Mutant accumulation rates varied considerably between different bacilli isolated from the same individual, and where high rates of mutation were observed, the mutational spectrum was consistent with reactive oxygen species-induced mutagenesis. Elevated bacterial mutation rates were identified in isolates from HIV-negative but not HIV-positive individuals, suggesting that they were immune-driven. These results support the model that mutagenesis of Mtb in vivo is modulated by the host environment, which could drive the emergence of variants associated with drug resistance in a host-dependent manner.

Dynamics of within-host Mycobacterium tuberculosis diversity and heteroresistance during treatment

BACKGROUND

Studying within-host genetic diversity of Mycobacterium tuberculosis (Mtb) in patients during treatment may identify adaptations to antibiotic and immune pressure. Understanding the significance of genetic heteroresistance, and more specifically heterozygous resistance-associated variants (RAVs), is clinically important given increasing use of rapid molecular tests and whole genome sequencing (WGS).

METHODS

We analyse data from six studies in KwaZulu-Natal, South Africa. Most patients (>75%) had baseline rifampicin resistance. Sputum was collected for culture at baseline and at between two and nine intervals until month six. Positive cultures underwent WGS. Mixed infections and reinfections were excluded from analysis.

FINDINGS

Baseline Mtb overall genetic diversity (at treatment initiation or major change to regimen) was associated with cavitary disease, not taking antiretroviral therapy if HIV infected, infection with lineage 2 strains and absence of second-line drug resistance on univariate analyses. Baseline genetic diversity was not associated with six-month outcome. Genetic diversity increased from baseline to weeks one and two before returning to previous levels. Baseline genetic heteroresistance was most common for bedaquiline (6/10 [60%] of isolates with RAVs) and fluoroquinolones (9/62 [13%]). Most patients with heterozygous RAVs on WGS with sequential isolates available demonstrated RAV persistence or fixation (17/20, 85%). New RAVs emerged in 9/286 (3%) patients during treatment. We could detect low-frequency RAVs preceding emergent resistance in only one case, although validation of deep sequencing to detect rare variants is required.

INTERPRETATION

In this study of single-strain Mtb infections, baseline within-host bacterial genetic diversity did not predict outcome but may reveal adaptations to host and drug pressures. Predicting emergent resistance from low-frequency RAVs requires further work to separate transient from consequential mutations.

FUNDING

Wellcome Trust, NIH/NIAID.

Metabolic Changes of Mycobacterium tuberculosis during the Anti-Tuberculosis Therapy

Tuberculosis, caused by Mycobacterium tuberculosis complex bacteria, remains one of the most pressing health problems. Despite the general trend towards reduction of the disease incidence rate, the situation remains extremely tense due to the distribution of the resistant forms. Most often, these strains emerge through the intra-host microevolution of the pathogen during treatment failure. In the present study, the focus was on three serial clinical isolates of Mycobacterium tuberculosis Beijing B0/W148 cluster from one patient with pulmonary tuberculosis, to evaluate their changes in metabolism during anti-tuberculosis therapy. Using whole genome sequencing (WGS), 9 polymorphisms were determined, which occurred in a stepwise or transient manner during treatment and were linked to the resistance (GyrA D94A; inhA t-8a) or virulence. The effect of the inhA t-8a mutation was confirmed on both proteomic and transcriptomic levels. Additionally, the amount of RpsL protein, which is a target of anti-tuberculosis drugs, was reduced. At the systemic level, profound changes in metabolism, linked to the evolution of the pathogen in the host and the effects of therapy, were documented. An overabundance of the FAS-II system proteins (HtdX, HtdY) and expression changes in the virulence factors have been observed at the RNA and protein levels.

A retrospective genomic analysis of drug-resistant strains of M. tuberculosis in a high-burden setting, with an emphasis on comparative diagnostics and reactivation and reinfection status

Background

Recurrence of drug-resistant tuberculosis (DR-TB) after treatment occurs through relapse of the initial infection or reinfection by a new drug-resistant strain. Outbreaks of DR-TB in high burden regions present unique challenges in determining recurrence status for effective disease management and treatment. In the Republic of Moldova the burden of DR-TB is exceptionally high, with many cases presenting as recurrent.

Methods

We performed a retrospective analysis of Mycobacterium tuberculosis from Moldova to better understand the genomic basis of drug resistance and its effect on the determination of recurrence status in a high DR-burden environment. To do this we analyzed genomes from 278 isolates collected from 189 patients, including 87 patients with longitudinal samples. These pathogen genomes were sequenced using Illumina technology, and SNP panels were generated for each sample for use in phylogenetic and network analysis. Discordance between genomic resistance profiles and clinical drug-resistance test results was examined in detail to assess the possibility of mixed infection.

Results

There were clusters of multiple patients with 10 or fewer differences among DR-TB samples, which is evidence of person-to-person transmission of DR-TB. Analysis of longitudinally collected isolates revealed that many infections exhibited little change over time, though 35 patients demonstrated reinfection by divergent (number of differences > 10) lineages. Additionally, several same-lineage sample pairs were found to be more divergent than expected for a relapsed infection. Network analysis of the H3/4.2.1 clade found very close relationships among 61 of these samples, making differentiation of reactivation and reinfection difficult. There was discordance between genomic profile and clinical drug sensitivity test results in twelve samples, and four of these had low level (but not statistically significant) variation at DR SNPs suggesting low-level mixed infections.

Conclusions

Whole-genome sequencing provided a detailed view of the genealogical structure of the DR-TB epidemic in Moldova, showing that reinfection may be more prevalent than currently recognized. We also found increased evidence of mixed infection, which could be more robustly characterized with deeper levels of genomic sequencing.

Characterization of Polymorphisms Associated with Multidrug-Resistant Tuberculosis by Whole Genomic Sequencing: A Preliminary Report from Mexico

Whole genome sequencing (WGS) has been proposed as a tool for the diagnosis of drug resistance in tuberculosis (TB); however, there have been few studies on its effectiveness in countries with significantly high drug resistance rates. This study therefore aimed to evaluate the effectiveness of WGS to identify mutations related to drug resistance in TB isolates from an endemic region of Mexico. The results showed that, of 35 multidrug-resistant isolates analyzed, the values of congruence found between the phenotypic drug susceptibility testing and polymorphisms were 94% for isoniazid, 97% for rifampicin, 90% for ethambutol, and 82% for pyrazinamide. It was also possible to identify eight isolates as potential pre-extensive drug resistant (XDR) and one as XDR. Twenty nine isolates were classified within L4 and two transmission clusters were identified. The results show the potential utility of WGS for predicting resistance against first- and second-line drugs, as well as providing a phylogenetic characterization of TB drug-resistant isolates circulating in Mexico.

Recent updates on drug resistance in Mycobacterium tuberculosis

Tuberculosis (TB) along with acquired immune deficiency syndrome and malaria rank among the top three fatal infectious diseases which pose threat to global public health, especially in middle and low income countries. TB caused by Mycobacterium tuberculosis (Mtb) is an airborne infectious disease and one-third of the world's population gets infected with TB leading to nearly 1·6 million deaths annually. TB drugs are administered in different combinations of four first-line drugs (rifampicin, isoniazid, pyrazinamide and ethambutol) which form the core of treatment regimens in the initial treatment phase of 6-9 months. Several reasons account for the failure of TB therapy such as (i) late diagnosis, (ii) lack of timely and proper administration of effective drugs, (iii) lower availability of less toxic, inexpensive and effective drugs, (iv) long treatment duration, (v) nonadherence to drug regimen and (vi) evolution of drug-resistant TB strains. Drug-resistant TB poses a significant challenge to TB therapy and control programs. In the background of worldwide emergence of 558 000 new TB cases with resistance to rifampicin in the year 2017 and of them, 82% becoming multidrug-resistant TB (MDR-TB), it is essential to continuously update the knowledge on the mechanisms and molecular basis for evolution of Mtb drug resistance. This narrative and traditional review summarizes the progress on the anti-tubercular agents, their mode of action and drug resistance mechanisms in Mtb. The aim of this review is to provide recent updates on drug resistance mechanisms, newly developed/repurposed anti-TB agents in pipeline and international recommendations to manage MDR-TB. It is based on recent literature and WHO guidelines and aims to facilitate better understanding of drug resistance for effective TB therapy and clinical management.

Deciphering drug resistance in Mycobacterium tuberculosis using whole-genome sequencing: progress, promise, and challenges

Tuberculosis (TB) is a global infectious threat that is intensified by an increasing incidence of highly drug-resistant disease. Whole-genome sequencing (WGS) studies of Mycobacterium tuberculosis, the causative agent of TB, have greatly increased our understanding of this pathogen. Since the first M. tuberculosis genome was published in 1998, WGS has provided a more complete account of the genomic features that cause resistance in populations of M. tuberculosis, has helped to fill gaps in our knowledge of how both classical and new antitubercular drugs work, and has identified specific mutations that allow M. tuberculosis to escape the effects of these drugs. WGS studies have also revealed how resistance evolves both within an individual patient and within patient populations, including the important roles of de novo acquisition of resistance and clonal spread. These findings have informed decisions about which drug-resistance mutations should be included on extended diagnostic panels. From its origins as a basic science technique, WGS of M. tuberculosis is becoming part of the modern clinical microbiology laboratory, promising rapid and improved detection of drug resistance, and detailed and real-time epidemiology of TB outbreaks. We review the successes and highlight the challenges that remain in applying WGS to improve the control of drug-resistant TB through monitoring its evolution and spread, and to inform more rapid and effective diagnostic and therapeutic strategies.

Beyond the SNP Threshold: Identifying Outbreak Clusters Using Inferred Transmissions

Whole-genome sequencing (WGS) is increasingly used to aid the understanding of pathogen transmission. A first step in analyzing WGS data is usually to define "transmission clusters," sets of cases that are potentially linked by direct transmission. This is often done by including two cases in the same cluster if they are separated by fewer single-nucleotide polymorphisms (SNPs) than a specified threshold. However, there is little agreement as to what an appropriate threshold should be. We propose a probabilistic alternative, suggesting that the key inferential target for transmission clusters is the number of transmissions separating cases. We characterize this by combining the number of SNP differences and the length of time over which those differences have accumulated, using information about case timing, molecular clock, and transmission processes. Our framework has the advantage of allowing for variable mutation rates across the genome and can incorporate other epidemiological data. We use two tuberculosis studies to illustrate the impact of our approach: with British Columbia data by using spatial divisions; with Republic of Moldova data by incorporating antibiotic resistance. Simulation results indicate that our transmission-based method is better in identifying direct transmissions than a SNP threshold, with dissimilarity between clusterings of on average 0.27 bits compared with 0.37 bits for the SNP-threshold method and 0.84 bits for randomly permuted data. These results show that it is likely to outperform the SNP-threshold method where clock rates are variable and sample collection times are spread out. We implement the method in the R package transcluster.

Mycobacterium tuberculosis Metabolism

Mycobacterium tuberculosis is the cause of tuberculosis (TB), a disease which continues to overwhelm health systems in endemic regions despite the existence of effective combination chemotherapy and the widespread use of a neonatal anti-TB vaccine. For a professional pathogen, M. tuberculosis retains a surprisingly large proportion of the metabolic repertoire found in nonpathogenic mycobacteria with very different lifestyles. Moreover, evidence that additional functions were acquired during the early evolution of the M. tuberculosis complex suggests the organism has adapted (and augmented) the metabolic pathways of its environmental ancestor to persistence and propagation within its obligate human host. A better understanding of M. tuberculosis pathogenicity, however, requires the elucidation of metabolic functions under disease-relevant conditions, a challenge complicated by limited knowledge of the microenvironments occupied and nutrients accessed by bacilli during host infection, as well as the reliance in experimental mycobacteriology on a restricted number of experimental models with variable relevance to clinical disease. Here, we consider M. tuberculosis metabolism within the framework of an intimate host-pathogen coevolution. Focusing on recent advances in our understanding of mycobacterial metabolic function, we highlight unusual adaptations or departures from the better-characterized model intracellular pathogens. We also discuss the impact of these mycobacterial "innovations" on the susceptibility of M. tuberculosis to existing and experimental anti-TB drugs, as well as strategies for targeting metabolic pathways. Finally, we offer some perspectives on the key gaps in the current knowledge of fundamental mycobacterial metabolism and the lessons which might be learned from other systems.

Whole genomic sequencing as a tool for diagnosis of drug and multidrug-resistance tuberculosis in an endemic region in Mexico

Background

Whole genome sequencing (WGS) has been proposed as a tool for diagnosing drug resistance in tuberculosis. However, reports of its effectiveness in endemic countries with important numbers of drug resistance are scarce. The goal of this study was to evaluate the effectiveness of this procedure in isolates from a tuberculosis endemic region in Mexico.

Methods

WGS analysis was performed in 81 tuberculosis positive clinical isolates with a known phenotypic profile of resistance against first-line drugs (isoniazid, rifampin, ethambutol, pyrazinamide and streptomycin). Mutations related to drug resistance were identified for each isolate; drug resistant genotypes were predicted and compared with the phenotypic profile. Genotypes and transmission clusters based on genetic distances were also characterized.

Findings

Prediction by WGS analysis of resistance against isoniazid, rifampicin, ethambutol, pyrazinamide and streptomycin showed sensitivity values of 84%, 96%, 71%, 75% and 29%, while specificity values were 100%, 94%, 90%, 90% and 98%, respectively. Prediction of multidrug resistance showed a sensitivity of 89% and specificity of 97%. Moreover, WGS analysis revealed polymorphisms related to second-line drug resistance, enabling classification of eight and two clinical isolates as pre- and extreme drug-resistant cases, respectively. Lastly, four lineages were identified in the population (L1, L2, L3 and L4). The most frequent of these was L4, which included 90% (77) of the isolates. Six transmission clusters were identified; the most frequent was TC6, which included 13 isolates with a L4.1.1 and a predominantly multidrug-resistant condition.

Conclusions

The results illustrate the utility of WGS for establishing the potential for prediction of resistance against first and second line drugs in isolates of tuberculosis from the region. They also demonstrate the feasibility of this procedure for use as a tool to support the epidemiological surveillance of drug- and multidrug-resistant tuberculosis.

Whole genome sequencing Mycobacterium tuberculosis directly from sputum identifies more genetic diversity than sequencing from culture

Background

Repeated culture reduces within-sample Mycobacterium tuberculosis genetic diversity due to selection of clones suited to growth in culture and/or random loss of lineages, but it is not known to what extent omitting the culture step altogether alters genetic diversity. We compared M. tuberculosis whole genome sequences generated from 33 paired clinical samples using two methods. In one method DNA was extracted directly from sputum then enriched with custom-designed SureSelect (Agilent) oligonucleotide baits and in the other it was extracted from mycobacterial growth indicator tube (MGIT) culture.

Results

DNA directly sequenced from sputum showed significantly more within-sample diversity than that from MGIT culture (median 5.0 vs 4.5 heterozygous alleles per sample, p = 0.04). Resistance associated variants present as HAs occurred in four patients, and in two cases may provide a genotypic explanation for phenotypic resistance.

Conclusions

Culture-free M. tuberculosis whole genome sequencing detects more within-sample diversity than a leading culture-based method and may allow detection of mycobacteria that are not actively replicating.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5782-2) contains supplementary material, which is available to authorized users.